3.1 Posture/balance disorder group data
3.1.1 Clinical features
The sample comprises 20 participants with PD, of whom 10 subjects were men (50%),average age (65.3±10.75 years old ); 10(50%) women with average age of (62.9±9.31 years old). The score of posture/balance was 3.12+3.13 according to MDS-UPDRS (3.12 posture balance, 3.13 posture). Nine patients with the total score of 2 and lower were labled as “Group1”. Eleven patients with the total score of 3 and more were labled as “Group2”. Group2 patients showed higher scores of UPDRS-II and UPDRS-III, being accepted higher levodopa dose, comparing with Group 1(See Table1). Moreover, Group2 possessed higher average age, a higher score of HAMA, HAMD and PSQI than those of Group1, but they showed no statistical difference.
Table1 Clinical characteristics between the groups according to posture/balance.
|
Group 1(n=9)
|
Group 2(n=11)
|
T/Z
|
P
|
Levodopa Dose
|
177.78±164.1
|
404.55±155.65
|
-3.164
|
0.005
|
UPDRSII
|
6.22±4
|
13.55±7.11
|
-2.776
|
0.012
|
UPDRSIII
|
19.27±7.2
|
29.22±10.56
|
-2.401
|
0.027
|
MMSE
|
29(28-29.5)
|
29(27-29)
|
-0.926
|
0.539
|
Disease duration
|
6.22±3.03
|
5.91±2.98
|
0.232
|
0.819
|
Age
|
60.56±7.54
|
67±10.9
|
-1.5
|
0.151
|
HAMA
|
3(1.5-17)
|
10(3-25)
|
-0.761
|
0.446
|
HAMD
|
9.67±9.54
|
13.55±9.77
|
-0.893
|
0.384
|
PSQI
|
7.33±5.59
|
8.1±4.28
|
0.258
|
0.74
|
The value of data that conforms to the normal distribution is displayed as the mean ± standard deviation, and the value of the data that does not conform to the normal distribution is displayed as the median(quartile). MMSE: Mini-Mental State Examination; HAMA: Hamilton Anxiety Scale; HAMD: Hamilton Depression Scale; PSQI: Pittsburgh sleep quality index.
3.1.2 DTI information
Fiber tracking was performed by DTI, Group2 showed lower value of FA and FA skeleton in ACC, PFC, and basal ganglia. For the parameters of RD and RD skeleton, MD and MD skeleton, AD and AD skeleton, there were higher value of the first five terms in ACC; and higher value of all parameters in PFC; higher value of RD skeleton, MD skeleton and AD skeleton in basal ganglia(P<0.05). They overlap with the difference in FA. There were no significant difference between the two groups in the comparison of the AD skeleton of ACC, the AD, MD, and RD of the basal ganglia (P>0.05).( Table 2)
Table 2: The difference between the groups in DTI index
|
Group 1(n=9)
|
Group 2(n=11)
|
T/Z
|
P
|
ACC FA
|
0.178(0.168-0.181)
|
0.16099(0.15-0.177)
|
-2.697
|
0.007
|
ACC FA skeleton
|
0.502±0.044
|
0.454±0.049
|
2.295
|
0.034
|
ACC RD
|
0.00101±0.0001
|
0.0012±0.0002
|
-2.903
|
0.009
|
ACC RD skeleton
|
0.00055±0.00005
|
0.00064±0.000097
|
-2.483
|
0.023
|
ACC MD
|
0.0011±0.0001
|
0.00129±0.00017
|
-2.887
|
0.01
|
ACC MD skeleton
|
0.0008(0.00079-0.00083)
|
0.00085(0.00082-0.0009)
|
-2.089
|
0.037
|
ACC AD
|
0.00128±0.00012
|
0.00148±0.00017
|
-2.843
|
0.011
|
ACC AD skeleton
|
0.00129(0.00123-0.00134)
|
0.00131(0.00126-0.00142)
|
-1.178
|
0.239
|
PFC FA
|
0.121(0.116-0.127)
|
0.112(0.107-0.121)
|
-2.545
|
0.011
|
PFC FA skeleton
|
0.4(0.378-0.411)
|
0.3695(0.338-0.39)
|
-2.621
|
0.009
|
PFC MD
|
0.00073±0.00006
|
0.000815±0.00008
|
-2.58
|
0.019
|
PFC MD skeleton
|
0.000576(0.0007-0.00079)
|
0.00081(0.00079-0.00086)
|
-2.773
|
0.006
|
PFC AD
|
0.00085±0.00007
|
0.00093±0.00008
|
-2.589
|
0.019
|
PFC AD skeleton
|
0.0011(0.00107-0.00113)
|
0.00115(0.00112-0.00116)
|
-2.393
|
0.017
|
PFC RD
|
0.00067±0.00006
|
0.00076±0.00008
|
-2.564
|
0.019
|
PFC RD skeleton
|
0.00057(0.00056-0.0006)
|
0.00063(0.0006-0.0007)
|
-2.545
|
0.011
|
Basal ganglia FA
|
0.247(0.236-0.256)
|
0.2235(0.212-0.242)
|
-2.089
|
0.037
|
Basal ganglia FA skeleton
|
0.419(0.407-0.448)
|
0.397(0.379-0.413)
|
-2.317
|
0.02
|
Basal ganglia AD
|
0.00146±0.0001
|
0.00162±0.0002
|
-1.743
|
0.098
|
Basal ganglia AD skeleton
|
0.00117(0.001-0.0012)
|
0.00122(0.0012-0.00127)
|
-2.241
|
0.025
|
Basal ganglia MD
|
0.0012±0.00012
|
0.00135±0.0002
|
-1.778
|
0.092
|
Basal ganglia MD skeleton
|
0.00078(0.00076-0.0008)
|
0.00085(0.0008-0.00088)
|
-2.925
|
0.003
|
Basal ganglia RD
|
0.00106±0.0001
|
0.0012±0.0002
|
-1.794
|
0.09
|
Basal ganglia RD skeleton
|
0.00059(0.00057-0.00062)
|
0.00066(0.00062-0.00069)
|
-3.001
|
0.005
|
The value of data that conforms to the normal distribution is displayed as the mean±standard deviation, and the value of the data that does not conform to the normal distribution is displayed as the median (quartile).The skeleton is a white matter skeleton, and differences between groups indicate significant results. If the difference of the corresponding data parameters is not on the skeleton, the meaning is not as good as the skeleton result. ACC: anterior cingulate cortex; PFC: prefrontal cortex; FA: fractional anisotropy; MD: mean diffusivity; RD: radial diffusivity; AD: axial diffusivity.
3.1.3 TCD ultrasound parameters between the two groups
Both groups of patients according to posture/balance disorder were detected with TCD. Group 2 showed significantly higher value of the PI and RI of bilateral MCA, ACA and ICA than those of the Group 1(P<0.01). But there were no statistic difference in the peak and mean velocities in MCA, ACA and ICA.(Table 3)
Table3: The difference between the groups in the blood supply
|
Group 1(n=9)
|
Group 2(n=11)
|
T/Z
|
P
|
Left MCA PI
|
0.925±0.058
|
1.069±0.18
|
-2.543
|
0.025
|
Left MCA RI
|
0.568±0.021
|
0.618±0.065
|
-2.427
|
0.031
|
Right MCA PI
|
0.87(0.815-1)
|
1.28(0.99-1.35)
|
-3.003
|
0.005
|
Right MCA RI
|
0.545(0.525-0.595)
|
0.69(0.62-0.71)
|
-3.202
|
0.002
|
Right ACA PI
|
0.834±0.14
|
1.19±0.2
|
-4.286
|
0.001
|
Right ACA RI
|
0.53±0.056
|
0.66±0.067
|
-4.27
|
0.001
|
Left ACA PI
|
0.85(0.81-0.92)
|
1.29(0.95-1.46)
|
-3.309
|
0.002
|
Left ACA RI
|
0.55±0.03
|
0.67±0.078
|
-5.043
|
<0.001
|
Right ICA terminal end PI
|
0.89±0.098
|
1.177±0.22
|
-3.426
|
0.003
|
Right ICA terminal end RI
|
0.55±0.04
|
0.65±0.07
|
-3.34
|
0.004
|
Left ICA terminal end PI
|
0.92±0.074
|
1.13±0.17
|
-3.152
|
0.006
|
Left ICA terminal end RI
|
0.57±0.028
|
0.64±0.056
|
-3.382
|
0.004
|
PI of left ICA siphon segment
|
0.925(0.89-0.95)
|
1.2(1.01-1.35)
|
-3.423
|
0.001
|
RI of left ICA siphon segment
|
0.57(0.555-0.58)
|
0.66(0.6-0.71)
|
-3.397
|
0.001
|
PI of right ICA siphon segment
|
0.92±0.1
|
1.23±0.22
|
-3.674
|
0.002
|
RI of right ICA siphon segment
|
0.57±0.04
|
0.67±0.07
|
-3.627
|
0.002
|
The value of data that conforms to the normal distribution is displayed as the mean ± standard deviation, and the value of the data that does not conform to the normal distribution is displayed as the median (quartile). MCA: middle cerebral arteries; ACA: anterior cerebral artery; ICA: internal carotid arteries; PI: pulsation index;RI: resistance index
3.1.4 Correlation analysis of TCD data and DTI parameters
Analysing the DTI data of the three ROI (ACC, PFC and basal ganglia), the DTI data in basal ganglia showed being associated with the PI and RI of MCA and ICA(Table4). FA skeleton negatively correlated with PI and RI of MCA and ICA, MD skeleton, both with AD and RD skeleton, showed positively correlating with PI and RI of MCA and ICA(P<0.01). As to the PFC and ACC, FA skeleton still showed negatively correlating with PI and RI of ACA and ICA, otherwise, MD skeleton and RD skeleton positively correlated with the PI and RI of ACA and ICA(Table5,6).
Table 4:Correlation analysis of basal ganglia DTI parameters and TCD data
|
Basal ganglia FA skeleton
|
Basal ganglia MD skeleton
|
Basal ganglia AD skeleton
|
Basal ganglia RD skeleton
|
|
R
|
P
|
R
|
P
|
R
|
P
|
R
|
P
|
Left ICA terminal end PI
|
-0.678
|
0.001
|
0.683
|
0.001
|
0.549
|
0.012
|
0.728
|
<0.001
|
Left ICA terminal end RI
|
-0.658
|
0.002
|
0.698
|
0.001
|
0.588
|
0.006
|
0.737
|
<0.001
|
Left MCA RI
|
-0.631
|
0.003
|
0.639
|
0.002
|
0.459
|
0.042
|
0.692
|
0.001
|
Left MCA PI
|
-0.592
|
0.006
|
0.604
|
0.005
|
0.425
|
0.062
|
0.654
|
0.002
|
Right MCA PI
|
-0.58
|
0.007
|
0.743
|
<0.001
|
0.658
|
0.002
|
0.723
|
<0.001
|
Right MCA RI
|
-0.601
|
0.005
|
0.757
|
<0.001
|
0.668
|
0.001
|
0.742
|
<0.001
|
Right ICA terminal end PI
|
-0.591
|
0.006
|
0.751
|
<0.001
|
0.614
|
0.004
|
0.744
|
<0.001
|
Right ICA terminal end RI
|
-0.59
|
0.006
|
0.76
|
<0.001
|
0.635
|
0.003
|
0.747
|
<0.001
|
FA: Fractional anisotropy; MD: average diffusivity; RD: radial diffusivity; AD: axial diffusivity; MCA: middle cerebral arteries; ICA: internal carotid arteries; PI: pulsation index; RI: resistance index.
Table 5:Correlation analysis of PFC DTI parameters and TCD data
|
PFC FA skeleton
|
PFC MD skeleton
|
PFC RD skeleton
|
|
R
|
P
|
R
|
P
|
R
|
P
|
Left ICA terminal end PI
|
-0.653
|
0.002
|
0.49
|
0.028
|
0.515
|
0.02
|
Left ICA terminal end RI
|
-0.655
|
0.002
|
0.483
|
0.031
|
0.512
|
0.021
|
Left ACA PI
|
-0.509
|
0.022
|
0.536
|
0.015
|
0.501
|
0.025
|
Left ACA RI
|
-0.513
|
0.021
|
0.536
|
0.015
|
0.501
|
0.025
|
Right ACA PI
|
-0.563
|
0.01
|
0.531
|
0.016
|
0.49
|
0.028
|
Right ACA RI
|
-0.556
|
0.011
|
0.523
|
0.018
|
0.481
|
0.032
|
Right ICA terminal end PI
|
-0.659
|
0.002
|
0.573
|
0.008
|
0.577
|
0.008
|
Right ICA terminal end RI
|
-0.658
|
0.002
|
0.569
|
0.009
|
0.58
|
0.007
|
PFC: prefrontal cortex; FA: Fractional anisotropy; MD: average diffusivity; RD: radial diffusivity; ACA: anterior cerebral artery; ICA: internal carotid arteries; PI: pulsation index; RI: resistance index
Table 6:Correlation analysis of ACC DTI parameters and TCD data
|
ACC FA skeleton
|
ACC MD skeleton
|
ACC RD skeleton
|
|
R
|
P
|
R
|
P
|
R
|
p
|
Left ICA terminal end PI
|
-0.477
|
0.034
|
0.655
|
0.002
|
0.6
|
0.005
|
Left ICA terminal end RI
|
-0.525
|
0.017
|
0.696
|
0.001
|
0.637
|
0.003
|
Left ACA PI
|
-0.493
|
0.027
|
0.543
|
0.013
|
0.514
|
0.02
|
Left ACA RI
|
-0.495
|
0.026
|
0.55
|
0.012
|
0.519
|
0.019
|
Right ACA PI
|
-0.473
|
0.035
|
0.521
|
0.018
|
0.49
|
0.028
|
Right ACA RI
|
-0.468
|
0.037
|
0.526
|
0.017
|
0.491
|
0.028
|
Right ICA terminal end PI
|
-0.557
|
0.011
|
0.69
|
0.001
|
0.62
|
0.004
|
Right ICA terminal end RI
|
-0.56
|
0.01
|
0.695
|
0.001
|
0.625
|
0.003
|
ACC: anterior cerebral artery; FA: Fractional anisotropy; MD: average diffusivity; RD: radial diffusivity; ACA: anterior cerebral artery; ICA: internal carotid arteries; PI: pulsation index; RI: resistance index
3.2 Analysis of PD-related pain in the patients
3.2.1 Clinical features of patients
The patients were classified as the group with PD-related pain and the group without PD-related pain by use of criteria of Ford classification method[6]. The patients with PD-related pain complained burning pain, acupuncture sensation, tightening sensation, knife cutting sensation. The patients were unable to identify the pain area or the pain area vacillates. Patients with PD-related pain (N=6, 30%) and patients without PD-related pain(N=14, 70%) are closely matched in age (65±8.94 vs.63.71±10.54, p=0.798). No significant differences in score of MDS-UPDRS Ⅱ and III either in “off” or “on” medication states between patients with or without PD-related pain. Again, there were no statistical difference in levodopa dose, HAMA, HAMD, age and disease duration. However, there was a statistically significant difference in PSQI scores(P<0.05) (Table7). Moreover, it showed no correlation between posture/balance disorder and PD-related pain using Chi-square test(P=1).
Table 7: Comparison of clinical information between patients with and without PD-related pain
|
PD-related pain
|
Without PD-related pain
|
T/Z
|
P
|
Levodopa Dose
|
291.67±182.8
|
307.14±204.62
|
0.16
|
0.875
|
PSQI
|
11.67±5.28
|
5.92±3.45
|
-2.856
|
0.011
|
HAMA
|
15.5(2.25-38.25)
|
6(1.75-16.25)
|
-1.116
|
0.264
|
HAMD
|
15.50±10.82
|
10.21±9.01
|
-1.135
|
0.271
|
UPDRSⅡ
|
4.5(3.75-15)
|
10(5.75-14.5)
|
-1.283
|
0.2
|
UPDRSⅢ
|
22.33±9.56
|
25.78±10.79
|
0.676
|
0.507
|
Age
|
65.00±8.94
|
63.71±10.54
|
-0.26
|
0.798
|
Disease duration
|
6.67±2.88
|
5.79±3.02
|
-0.606
|
0.552
|
The value of data that conforms to the normal distribution is displayed as the mean±standard deviation, and the value of the data that does not conform to the normal distribution is displayed as the median (quartile). HAMA: Hamilton Anxiety Scale; HAMD: Hamilton Depression Scale; PSQI: Pittsburgh sleep quality index.
3.2.2 Analysis on fiber bundle tracking
The fiber bundles were tracked from midbrain to the basal ganglia and PFC.The average FA value of the tracked fiber bundle from the midbrain to basal ganglia was significantly higher in the patients with PD-related pain than that in patients without PD-related pain. Moreover, there were more fibers from midbrain to basal ganglia in patients without PD-related pain (p<0.05)(Table8 and Fig.2-A). No statistic difference were found in fibers amount from midbrain to PFC between the two groups, they all showed less fibers left there(Fig.2-B).
Analysis of the white matter fiber bundles in the locus coeruleus area showed that there were no significant differences in the number of fiber bundles, the maximum FA value and the maximum fiber length between the two groups(p>0.05)(Table 8 and Fig.3).
Table 8:The difference between the groups in fiber bundle tracking
|
With PD-related pain
|
Without PD-related pain
|
T/Z
|
P
|
Number of midbrain to basal ganglia fibers
|
176.5(71.75-339)
|
30(6.5-93.5)
|
-2.187
|
0.029
|
Average fiber length from midbrain to basal ganglia
|
81.2886(67.91-87.9)
|
53.2199(26.52-84.18)
|
-1.32
|
0.187
|
Average FA value of midbrain to basal ganglia fibers
|
0.4947(0.47-0.51)
|
0.42265(0.31-0.47)
|
-2.434
|
0.015
|
Number of fiber bundles passing through locus coeruleus
|
111.5±23.72
|
105.6±40.4
|
-0.333
|
0.743
|
The maximum FA value of the fiber bundle walking in the locus coeruleus nucleus
|
0.614±0.074
|
0.6±0.039
|
-0.566
|
0.578
|
Maximum length of fiber bundles running in the locus coeruleus nucleus
|
83.59±27.75
|
70.98±22.57
|
-1.071
|
0.298
|