Introduction
Parkinson’s disease (PD) is the second most common neurodegenerative disease [
1]. Clinically, PD is characterized by progressive motor symptoms, such as bradykinesia, rigidity, resting tremor, and postural instability, and various non-motor symptoms [
2]. The pathophysiology of PD has yet to be fully elucidated, and this disease has no preventative or curative treatments. Recently, psychoradiology is emerging as the study of brain diseases through a variety of non-invasive imaging techniques. Diffusion tensor imaging (DTI), a widely used magnetic resonance imaging (MRI) technique that visualizes and measures water diffusion, is particularly useful because of its high sensitivity to changes in white matter (WM) integrity [
3]. Fractional anisotropy (FA), the most commonly used DTI parameter, is positively associated with WM anisotropy and can detect microstructural abnormalities of WM at very early stages [
4,
5].
Many DTI studies of PD have identified decreased FA in various brain regions, including corpus callosum (CC) [
6–
9], cerebellum [
8,
10–
12], and superior and inferior longitudinal fasciculus (ILF) [
9,
10], but others have found no significant FA differences between patients with PD and healthy controls (HC) [
13–
15]. These discrepancies may be attributable to factors such as small or heterogeneous study samples and methodological differences. To gain insights into the pathophysiology, a meta-analysis has a good chance of detecting consistent regions of FA alterations in PD.
Earlier reviews and meta-analyses confirmed that symptoms of PD are a result of the degeneration of the substantia nigra [
16–
18]. Although the pathologic process is mainly attributed to disruptions in the nigrostriatal dopamine system, previous neuroimaging studies in PD have focused on the substantia nigra using a seed-point method. By examining pre-defined regions of interest (ROI), some important cerebral alterations may have been missed. To overcome region placement preference bias, we chose a widely used whole-brain method to analyze the WM changes in PD. Previous DTI meta-analyses [
19] did not directly test the hypothesis of WM changes in patients with PD compared with HC, and confounding factors such as medication status, disease stage, and methodological differences were not fully considered. A comprehensive meta-analysis is needed to explore the WM microstructure abnormalities and to investigate the effects of symptom severity and other clinical characteristics on regional WM alterations.
The aim of this meta-analysis was to define the most consistent WM microstructural changes, reflected by FA, in patients with PD using published whole brain DTI studies. To reduce the effects of heterogeneity, several subgroup analyses were conducted. Meta-regression analyses were used to evaluate the influence of clinical and demographic characteristics.
Materials and methods
Data sources, study selection, and quality assessment
This retrospective study was approved by the local institutional review committee. In accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement [
20], a comprehensive search of studies published before May 14, 2018 was conducted in the PubMed, Web of Science, and Embase databases using the keywords “Parkinson” OR “Parkinson’s disease” OR “PD”; AND “diffusion tensor” OR “diffusion tensor imaging” OR “DTI.” The references of these studies and relevant review articles were checked for additional relevant studies. Studies that satisfied the following conditions were included in the meta-analysis: (1) patients were diagnosed with idiopathic PD; (2) FA was compared between patients and HC; (3) three-dimensional coordinates (Montreal Neurological Institute (MNI) or Talairach) were reported for the whole-brain FA analysis; (4) significant results were reported using thresholds for significance corrected for multiple comparisons or uncorrected with spatial extent thresholds; and (5) the study was published as an original article (not as a letter or abstract) in a peer-reviewed English-language journal. Data sets were excluded if they explicitly included patients diagnosed with comorbid neurological or psychiatric diseases (e.g., cognitive impairment or depression) [
21]. To avoid sample overlaps: (1) for longitudinal studies, only baseline data were included [
21]; (2) for studies reporting both on- and off-state results, only off-state data sets were included [
21]; (3) for multiple studies using the same patient group, only the largest sample was included. Finally, the meta-analysis included 27 data sets from 24 studies with 808 patients with PD (429 men and 379 women) and 760 HC (393 men and 367 women) (Fig. 1).
Two of the authors independently rated each included study for quality and completeness using a 12-point checklist adapted from previous published meta-analyses [
22] (Table S1). Any discrepancies were resolved by a third investigator.
Voxel-wise meta-analysis
We analyzed FA differences in WM between patients with PD and HC using anisotropic effect size–signed differential mapping (AES-SDM), a voxel-based meta-analytic approach (see Supplementary Material). As detailed elsewhere [
23–
25], we extracted peak coordinates and effect sizes (e.g., t-values) of FA differences between patients with PD and HC from each data set. For each data set, we recreated a standard MNI map of FA differences using an anisotropic Gaussian kernel, from which the mean map was generated by voxel-wise calculation, weighted by sample size, intra-data set variability, and between-data set heterogeneity. To optimally balance false positive and negative findings, we used the default SDM kernel size and thresholds (full width at half maximum (FWHM) = 20 mm, uncorrected
P = 0.005 was used as the main threshold, peak height
Z = 1, cluster extent= 10 voxels) [
23,
24]. This FWHM kernel is intended to assign indicators of proximity to reported coordinates but not to smooth any image that is different in nature [
25]. We adopted the tract-based spatial statistics (TBSS) template included in AES-SDM to allow the combination of voxel-based analysis (VBA) and TBSS studies [
26]. We used MRIcron software to visualize AES-SDM maps overlaid onto a high-resolution brain image template created by the International Consortium for Brain Mapping.
Jackknife sensitivity analysis
Following preprocessing, we performed a whole-brain voxel-based jackknife sensitivity analysis to test the robustness of the findings by iteratively repeating the analysis, excluding one data set each time. If a brain region remains significant in all or most of the combinations of studies, the finding is considered highly replicable [
23].
Analyses of heterogeneity and publication bias
We conducted a heterogeneity analysis using a random effects model with Q statistics to explore unexplained between-study variability in the results; heterogeneous brain regions were obtained using the default SDM kernel size and thresholds [
25]. We also performed Egger’s test for publication bias by extracting the values from statistically significant relevant peaks between patients with PD and HC [
27].
Subgroup analyses
Several subgroup meta-analyses were conducted to analyze clinical and methodological differences between these studies. This subgroup analysis included only studies with defined clinical and methodological homogeneity and was conducted for studies that investigated medicated and medication-free (including both medication-naïve and off-state) patients, studies that investigated early stage (Hoehn and Yahr (H&Y) stage 1–2.5) patients [
28], studies where the number of acquisition diffusion directions was≥30 and <30 [
29], and studies that used VBA or TBSS approaches.
Meta-regression analyses
We used meta-regression analyses to examine the effects of age, percentage of female patients, illness duration, H&Y stage, Unified Parkinson’s Disease Rating Scale (UPDRS) Part III score, and levodopa equivalent daily dose (LEDD) of patients with PD. We decreased the probability threshold to 0.0005 and cluster extent= 10 voxels, as described previously [
26]. Regression plots were visually inspected to discard the fits driven by too few studies [
23,
24].
Results
Included studies and sample characteristics
The meta-analysis included 24 studies [
6–
15,
30–
43], and 3 of them included 2 independent databases [
8,
13,
15], resulting in 27 data sets in all. Table 1 summarizes the demographic, clinical, and imaging information of these studies, which provided a total of 41 coordinates of altered FA in patients with PD versus HC. Ten data sets from 8 studies and 13 data sets from 12 studies explicitly indicated patients in the medication off-state and on-state, and 1 study included medication-naïve patients. Three studies did not provide any information about the medication status (on- or off-state) at DTI acquisition. The patients recruited in 6 data sets from 5 studies were at their early stage. The number of acquisition diffusion directions was≥30 in 11 data sets from 10 studies and <30 in 16 data sets from 14 studies. Seventeen studies used the TBSS approach, 5 used VBA, 1 used whole-brain-based spatial statistics, and 1 used gray-matter-based spatial statistics. Five of the 24 studies were performed on 1.5T scanners and 19 on 3T scanners. The quality scores, ranging from 9.5 to 12 (mean 10.7), demonstrated that the included studies were of high quality (Table 1).
FA differences in the main meta-analysis
Meta-analysis revealed decreased FA in patients with PD relative to HC, in the left middle cerebellar peduncles (MCP), CC, left inferior fronto-occipital fasciculus (IFOF), and right ILF (Fig. 2, Table 2).
Jackknife sensitivity analysis
Whole-brain jackknife sensitivity analysis showed that the most robust results were obtained for FA decreases in the WM regions of the left MCP and left CC, which were preserved throughout all combinations; right ILF and right CC were relatively replicable, being preserved throughout all but one combination; FA decreases in the left IFOF remained significant in all but two combinations (Table 3).
Analyses of heterogeneity and publication bias
The left CC with altered FA had significant statistical heterogeneity among studies (P<0.005; Table 4). Egger’s tests for publication bias were not significant with respect to FA differences in the left MCP (P = 0.107), left CC (P = 0.122), left IFOF (P = 0.126), right CC (P = 0.059), or right inferior longitudinal fasciculus (P = 0.057) (Table 2).
Subgroup analyses
The medication-free subgroup (comprising off-state (n = 10 data sets) and medication-naïve patients (n = 1 data set)) included 11 data sets comprising 374 patients with PD and 328 HC. The analysis revealed FA decreases in the right olfactory cortex. The medicated group (n = 13 data sets) comprised 340 patients with PD and 347 HC. The main results remained unchanged. We then used a threshold of LEDD= 400 mg/day to define two medicated subgroups. In the LEDD < 400 mg/day subgroup analysis (n = 5 data sets), comprising 205 patients with PD and 144 HC, the main results remained unchanged; by contrast, the LEDD≥400 mg/day subgroup (n = 7 data sets), comprising 122 patients with PD and 157 HC, revealed FA increases in the bilateral CC.
In the early stage disease (n = 6 data sets) and VBA (n = 5 data sets) subgroups, the main results remained unchanged, but the TBSS (n = 20 data sets) subgroup analysis failed to detect between-group differences. In the < 30 diffusion directions subgroup analysis (n = 16 data sets), four of the main results remain unchanged, whereas the≥30 diffusion directions subgroup (n = 11 data sets) only revealed one main result in the left IFOF. Detailed results of subgroup analyses are shown in Table 5 and Figs. S1–S6.
Meta-regression analyses
The mean age of patients, percentage of women with PD (in the whole PD group in each study), mean illness duration, mean H&Y stage, mean UPDRS-III score, and mean LEDD of patients with PD were not significantly associated with PD-related WM FA changes, at least linearly.
Discussion
The main abnormalities and their potential significance
Our pooled meta-analysis of DTI studies revealed significant WM microarchitecture alterations in patients with PD compared with HC as indicated by the decreased FA values. A consistent decrease in FA was observed in the left MCP, CC, left IFOF, and right ILF. Jackknife sensitivity analysis, heterogeneity analysis, and Egger’s tests confirmed the robustness of these findings. Although the subgroup analysis of medication-free patients did not reveal decreased FA as in the main analysis, it detected an additional cluster in the right olfactory cortex. Meta-regression analyses revealed no significant association between WM alterations and relevant demographic and clinical variables.
Given the role of the MCP as the main afferent pathway from the cortex to the cerebellum and additional findings of diffusion changes in brainstem and pons [
44], these data emphasize the role of the cortico-ponto-cerebellar loop in PD pathology. This loop contains afferent axons from the pontine nuclei [
45] and mediates complex brain functions concerning movement and cognition. PD features motor dysfunction such as tremor, rigidity, bradykinesia, gait disturbance, and postural instability [
46]. Studies of patients with freezing of gait found WM abnormalities in the pontine-cerebellar tracts [
8,
47,
48], which might explain their reduced ability to recruit the cerebellum. Functional studies in PD found hypoactivation in the cerebellum during a motor task [
49]. This neurodegeneration in MCP is consistent with the critical role of the cerebellum for motor performance [
46,
50]. Moreover, the subgroup analysis of PD patients with early stage disease detected decreased FA in MCP, an area with a significant decline in diffusivity in early PD, correlating with motor progression [
51].
Being the largest WM bundle in the human brain, the CC connects the bilateral cerebral hemispheres [
52]; therefore, this finding, suggesting hemispheric involvement [
53], is consistent with resting-state studies reporting decreased interhemispheric functional connectivity in PD [
54–
56]. In previous DTI studies in PD, altered FA was reported in different subdivisions of the CC, such as the genu, which is usually associated with cognitive impairment [
37,
38,
57]. Microstructural impairment of CC has been reported in cognitively normal patients with PD, which points toward an affected interhemispheric connectivity that possibly shapes the onset of cognitive deficits [
58]. In our meta-analysis, the abnormal microstructure of the posterior CC was detected in the pooled analysis and the subgroups of patients with early stage disease, medicated patients, VBA method, and <30 diffusion directions; this finding is consistent with reports of a smaller posterior CC volume in patients with PD [
59]. Transcallosal motor fiber bundles cross the CC in the posterior part [
60]. The anatomic projections found in this study are compatible with previous observations that CC abnormalities are related to postural instability and freezing of gait in PD [
60], which becomes more severe with motor worsening [
61]. Thus, the importance of neurodegeneration of CC in PD is also reflected through its role in motor features of the disease.
The ILF is an associative bundle connecting the occipital and temporal lobes. Its long fibers connect visual areas to the amygdala and hippocampus, the main components of the limbic system related to emotional behavior [
62]. It is involved in face recognition [
63], visual perception [
64], and visual memory [
65]. The IFOF is a ventral associative bundle connecting the ventral occipital lobe and the orbitofrontal cortex [
66]. Its occipital component runs parallel to the ILF. On approaching the anterior temporal lobe, its fibers gather together and enter the external capsule dorsally to the fasciculus. The IFOF is involved in reading [
67] and visual processing [
68]. The functions of these two bundles are not yet completely understood in PD. Guan and colleagues found that the ILF is associated with motor impairment and global disease severity in PD [
30], which has also been reported in relation to the tremor-dominant (TD) subtype of PD, indicating probably different underlying pathologies between TD and other motor subtypes [
69]. FA values in the IFOF were negatively correlated with bradykinesia of PD [
70], and disruption of this structure might contribute to freezing of gait due to impaired visuospatial processing [
71]. However, a recent DTI study showed significant differences in these two bundles between patients with postural instability gait disorder and HC [
72]. The role of these two bundles in differentiating WM changes in PD motor subtypes requires further investigation. Notably, previous studies also found relationships between these two bundles and non-motor symptoms of PD, such as depression [
73] and cognitive dysfunction [
74]. Although studies with comorbid depression and cognitive dysfunction were excluded in the current meta-analysis, the potentially confounding effects could not be completely ruled out. Moreover, although detailed depression scales were not used in these included studies, depression is a frequent non-motor symptom in patients with PD, and its prevalence can be as high as 90% [
75].
The quantitative statistical meta-analyses have been carried out on DTI studies in PD [
17]. Although these earlier meta-analyses provide valuable insights into the underlying WM pathology, differences in the ROI placement likely result in significant heterogeneity between studies. Relative to ROI-based analyses, coordinate-based analyses are less susceptible to errors relating to the lack of spatial distinction [
76]. There is relatively poor consistency of results between previous ROI meta-analyses and our current coordinate-based meta-analysis, which has implications for neuroimaging meta-analyses of other disorders, such as amnestic mild cognitive impairment [
76] and post-traumatic stress disorder [
77]. Several factors might contribute to these differences. First, ROI studies typically report on a subset of brain regions and are likely to suffer from publication bias. Second, where seed-based d mapping uses coordinate data, the effect size is biased toward zero in brain regions where there are no significant clusters.
Subgroup analyses
The subgroup analysis of medication-free patients detected a cluster in the right olfactory cortex, which is considered an early affected site in the Braak model of temporal degeneration in PD [
78]. A recent DTI study suggested that olfactory regions were particularly efficient at distinguishing drug-naïve patients with PD from HC [
79]. There is evidence of a strong link between olfactory discrimination and orbitofrontal cortex [
80]. Decreased FA in the olfactory cortex may be related to the olfactory dysfunction, which is present in ~90% of early stage PD [
81]. However, we did not observe this alteration in the olfactory cortex in the pooled whole-group results. The reason for this discrepancy is unclear, although changes associated with medication may normalize, or at least obscure, the intrinsic changes in these regions. Medication exposure is an important potential confounder, and understanding the effect of medications on white-matter abnormalities in patients with PD is critical for the interpretation of results. In functional MRI studies, levodopa can normalize altered functional connectivity in PD [
82]. Medications have also been shown to affect cerebral WM; for example, benztropine used in the treatment of tremor in PD promotes re-myelination [
83]. In the current meta-analysis, the main results of the LEDD < 400 mg/day subgroup analysis remain unchanged, whereas the results of the LEDD≥400 mg/day subgroup show increased FA in the bilateral CC, sharing these two clusters with the results of the pooled meta-analysis. This observation seems to fit a “normalization” effect. Although other diffusion MRI studies have shown that medication
per se has little influence on WM integrity or does not alter diffusion [
38,
84], the effect of dopaminergic medication on WM is still unclear. The differential influence of dopamine medication is likely determined by additional factors. Inter-individual differences and disease severity should also be considered. Further studies that are specifically designed to detect the effect of medication exposure on WM changes in patients with PD are needed.
Scan acquisition parameters also affect FA measurements individually and in group comparisons. The minimum requirement for a qualitative DTI calculation is 6 non-collinear diffusion directions, whereas robust estimation of FA and tensor orientation requires at least 30 unique and evenly distributed sampling orientations [
29]. DTI studies with less than 30 diffusion directions tend to overestimate FA values, possibly affecting group-wise differences. However, our subgroup analysis of≥30 diffusion directions showed the same robust decrease in FA in the left IFOF as in the main analysis. Nevertheless, future DTI meta-analyses should consider the number of diffusion directions.
Subgroup analysis of the VBA studies revealed clusters that were identified in the pooled analysis; however, the subgroup analysis of the TBSS studies did not. This difference may be a methodological issue. TBSS was developed to analyze diffusion data, and it isolates the central core of WM tracts with the highest FA and reports significant clusters within that WM skeleton [
85]; however, it may miss abnormalities in the peripheral WM [
86]. VBA investigates the integrity of all WM voxels, including peripheral regions of WM and fiber crossings [
87]. Differences in findings between the two methods are thus to be expected. Moreover, given the small number of VBA studies (
n = 5), we cannot completely rule out the possibility of chance findings.
Limitations of this study
First, methodological heterogeneity cannot be completely excluded: it was impossible to perform analyses or subgroup analyses of all relevant variables. Second, the medication-free patients of subgroup meta-analysis were mostly medicated patients in the off-state at scanning. Clearly, the best way to minimize the effects of treatment is to study medication-naïve patients; however, this subgroup was small. The results of the subgroup meta-analysis for medication-free data sets should therefore be interpreted with caution, considering the possible methodological differences between the studies, such as diffusion directions and analysis methods. Third, as is inherent in all coordinate-based methods, our meta-analysis used the coordinates from the included studies rather than the original t-statistic maps, which somewhat limits the accuracy of the results. Fourth, although we performed several subgroup meta-analyses, including VBA and different LEDD medicated data sets, for the early stage, the results should be interpreted with caution. The small sample size of these subgroup meta-analyses limits the generalizability of the results.
In conclusion, this meta-analysis demonstrates a consistent pattern of decreased FA in the left MCP, CC, IFOF, and ILF. In particular, FA decreases in the MCP and CC in PD were robustly present in the sensitivity analyses. These results suggest that the altered WM microstructures in the MCP and CC may serve as a potential neuroimaging biomarker of PD. Furthermore, the subgroup analyses show that medication status, analysis approaches, and the number of diffusion directions have an important impact on the findings, which should be treated with caution in future meta-analysis studies.