Alterations of Interhemispheric Functional Connectivity in Parkinson’s Disease With Depression: A Resting-State Functional MRI Study

Background: Depression is the most common non-motor symptom in patients with Parkinson’s disease (PD) with unknown mechanisms, but the diagnostic criteria of PD with depression (PDD) are not uniform.

Purpose: The aim of the study was to investigate interhemispheric interactions between PDD patients and patients with PD without depression (PDND).

Methods: The voxel-mirrored homotopic connectivity (VMHC) combined with the seed-based method was used to investigate intrinsic resting-state functional connectivity (RSFC) in 33 PDD patients, 60 PDND, and 47 healthy controls (HCs).

Results: PDD patients exhibited a decreased VMHC in the bilateral medial frontal gyrus and paracentral lobule (MFG/PCL) than did PDND patients. Parkinson’s disease with depression had a decreased VMHC in the bilateral precentral gyrus than had PDND and HC (p < 0.05). Parkinson’s disease with depression had a decreased homotopic RSFC from the medial frontal gyrus (MFG)/PCL to the contralateral supplementary motor area (SMA) than had PDND (p < 0.05). The decreased homotopic RSFC from the right MFG/PCL to the left SMA was negatively correlated with Hamilton Depression Rating Scale scores (p < 0.05), but not with illness duration, Beck’s Depression Inventory, and Unified Parkinson’s Disease Rating Scale in PD patients.

Conclusions: Our findings indicated that the occurrence of depression in Parkinson’s disease is associated with the dysfunctional connectivity from the MFG/PCL to the contralateral SMA, which could be used as potential neuroimaging markers for the diagnosis of depression in PD patients.

Introduction

Depression is the most common and severe non-motor symptom in patients with Parkinson’s disease (PD), with a prevalence of around 30–50% (Ossowska and Lorenc-Koci, 2013). The incidence of PD with depression (PDD) increases with the progression of the disease (Rickards, 2005). It is commonly thought that depression is associated with reduced functioning and cognitive impairment; this is due to how such impairment may worsen motor and other non-motor symptoms, as well as the general quality of life in PD patients (Aarsland et al., 2011). Cumulative evidence suggests that PDD may be a consequence of the neurodegenerative process rather than simply a reaction to psychological distress (McDonald et al., 2003; Langston, 2006). Depression could be the initial symptom in PD (Postuma et al., 2012), associating with a fast decline of cognition and motor abilities (Ravina et al., 2007). However, the pathophysiological mechanisms of depression in PD patients are not fully understood. In addition, depression in PD patients is frequently unrecognized and untreated (Pachana et al., 2013) owing to the overlap of the signs and symptoms of depression with motor and other non-motor symptoms. Thus, a comparative study of brain structures and functions in PDD patients with and patients with PD without depression (PDND) could uncover the mechanisms of depression occurring in PD patients.

With the advanced neuroimaging approaches, many previous studies found that PDD patients have specific functional and structural abnormalities in the prefrontal, basal ganglia, and limbic regions. For example, studies on neural metabolic activity in the brain of PD patients showed abnormal metabolic activity in the frontal lobe, striatum, and subcortical or limbic regions including the thalamus, amygdala, hippocampus, anterior cingulate cortex, and insula in PDD patients (Huang et al., 2013; Vriend et al., 2014; Frosini et al., 2015; Qamhawi et al., 2015). Structural neuroimaging studies showed a gray matter loss and white matter reduction in the prefrontal, temporal, and some limbic regions (Matsui et al., 2007; Feldmann et al., 2008; Kostić et al., 2010). Studies on neural connectivity showed decreased connectivity in the bilateral anterior cingulate cortex, thalamus, and multiple tracts connecting to the left frontal and deep temporal lobes (Li et al., 2010; Huang et al., 2014).

Recently, resting-state functional magnetic resonance imaging (rs-fMRI) has been suggested as the most promising method to study the neural network mechanisms. However, most previous studies in PDD patients focus on regional deficits (Wen et al., 2013) and/or functional connectivity abnormalities of a specific network (Hu et al., 2015). In contrast, rare studies have focused on the integrity of the neural network. Voxel-mirrored homotopic connectivity (VMHC) is a data-driven voxel-based intrinsic functional connectivity (iFC) approach. Voxel-mirrored homotopic connectivity combines both the amplitude of low frequency fluctuation and regional homogeneity that can reflect the specific patterns of interhemispheric disconnection. Using this method, Zhu et al. (2016) found reduced VMHC values in the bilateral dorsolateral prefrontal cortex and calcarine cortex in PDD patients. However, Zhu et al. (2016) study did not use the seed-based iFC analysis.

We therefore hypothesized that PDD and PDND patients may have altered functional connectivity in some networks. In this study, we used a seed-based iFC analysis to reveal the heterotopic connectivity patterns in seed regions with altered VMHC in PDD patients. In addition, we also expected that VMHC measurements and resting-state functional connectivity (RSFC) would be clinically relevant, so the relationship between the altered VMHC values/RSFC and clinical variables was analyzed. This study can provide a more complete understanding on the abnormal interhemispheric synchronization in PDD and PDND patients.

Materials and Methods

Participants

Participants were recruited from The Second Xiangya Hospital, Central South University. Sixty-four PDND patients and 38 PDD patients were recruited from December 2015 to May 2019. Motor symptom was evaluated by the motor part of the Unified Parkinson’s Disease Rating Scale (UPDRS-III) and Hoehn and Yahr (H&Y) Scale. Parkinson’s disease with depression patients were diagnosed with Statistical Manual of Mental Disorders, Fifth Edition (DSM-V) by an experienced psychiatrist. The severity of depression was quantified with the 17-item Hamilton Depression Rating Scale (HDRS-17) and the 21-item Beck’s Depression Inventory (BDI). All participants underwent a semi-structured interview and were assessed with the Mini-Mental State Exam (MMSE). An MMSE score >17 for illiterate subjects, >20 for primary school literates, and >23 for junior high school and higher education literates were defined as a normal MMSE score in our patients. The inclusion criteria were as follows: (1) patient met the UK Parkinson’s Disease Society Brain Bank Criteria (Daniel and Lees, 1993); (2) patient was right-handed Han Chinese; (3) patient did not use anti-PD or antidepressant medications for 12 h; (4) patient had a motor symptoms of H&Y staging 1–3; (5) patient had a HDRS-17 score higher than 7 (PDD patients were diagnosed by DSM-V (Wang et al., 2018); and (6) patient had an MMSE score not lower than the corresponding education level (Li et al., 2016). Patients were excluded if they had (1) a history of head injury, stroke, or other neurologic or psychiatric disease; (2) abnormal MMSE scores; and (3) any disorders that may interfere with the assessment of the manifestation of PD. Fifty right-handed healthy Han Chinese with matched age, education, and sex were recruited for controls.

Detailed neuropsychological examinations including MMSE, HDRS-17, and BDI were used to exclude dementia and depression in healthy controls (HCs). Twelve participants (five PDD, four PDND, and three HCs) were excluded from analyses. This study was approved by Medical Research Ethical Committee of The Second Xiangya Hospital. Written informed consent was obtained from all participants prior to enrollment.

Imaging Data Acquisition

Imaging data were acquired on a Siemens Skyra 3T MRI scanner at The Second Xiangya Hospital of Central South University. All participants were instructed to lie supine with eyes closed and no thinking but to avoid falling asleep. Their heads were fixed snugly with foam pads and straps to minimize head movement. Resting-state functional magnetic resonance imaging images were obtained using an echo planar imaging with 39 axial slices in 3.5-mm thickness, no gap, 2,500-ms repetition time, 25-ms echo time, 3.8 × 3.8 × 3.5-mm voxel size, 90° flip angle, 240-mm field of view, 64 × 64 data matrix, and 200 volumes. In addition, three dimensional T1-weighted magnetization-prepared rapid gradient echo sagittal images were acquired using the follow parameters: 176 slices, 1,900-ms repetition time, 2.01-ms echo time, 1.00-mm slice thickness, 1.0 × 1.0 × 1.0-mm voxel size, 9° flip angle, 900-ms inversion time, 256-mm field of view, and 256 × 256 matrix.

Magnetic Resonance Imaging Data Preprocessing

The rs-fMRI data were analyzed using the Data Processing and Analysis of Brain Imaging (DPABI) v2.1 toolbox (Yan et al., 2016), based on the Statistical Parametric Mapping software (SPM12)¹ for the MATLAB platform (The MathWorks Inc., Natick, MA, United States) as previously described (Zhou et al., 2018). Preprocessing comprised the following: (1) removal of the first 10 time points; (2) slice-time and motion corrections; (3) registration of the high-resolution individual T1-weighted images to the mean functional data; (4) segmentation of the high-resolution T1-weighted images with the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) toolkit and generation of a group template; (5) transformation and normalization of the resulting aligned data to the Montreal Neurological Institute (MNI) space with the segmented information from DARTEL; (6) resampling of the normalized rs-fMRI data to 3 × 3 × 3-mm voxels and spatial smoothing with an isotropic 6-mm full width at half maximum (FWHM) Gaussian kernel; (7) nuisance linear regression with the white matter, cerebrospinal fluid, and Friston-24 head motion parameters (including 6 head motion parameters and their historical effects as well as the 12 corresponding squared items) (Yan et al., 2013); and (8) removal of the linear trend and temporal bandpass (0.01–0.1 Hz) filtering. Finally, each subject was limited to a maximum displacement in any cardinal direction (x, y, z) of less than 1.5 mm and a maximum rotation (x, y, z) of less than 1.5° during each rs-fMRI scan. The group differences in head motion between the PDD patients, PDND patients, and the HC group were evaluated according to the published criteria (Van Dijk et al., 2012; Table 1). Six participants who displayed excessive head motion (translation >1.5 mm or rotation >1.5°) were excluded from further analysis.

TABLE 1

Table 1. Demographic and clinical characteristics of the participants (mean ± SD).

Voxel-Mirrored Homotopic Connectivity

To obtain the interhemispheric connectivity at the single-voxel level, VMHC was performed using the DPABI toolkit. A group-specific symmetrical left–right hemisphere template was generated from the normalized T1 images of all subjects to minimize the geometric differences between hemispheres (Zuo et al., 2010). The same transformation of the T1-weighted images was then applied to the individual preprocessed resting-state functional images. For each subject, the homotopic connectivity coefficient was computed as the Pearson correlation between every pair of mirrored interhemispheric voxels’ time series. The resultant values constituted the VMHC and were used for group analyses.

Seed-Based Intrinsic Functional Connectivity

We used a seed-region approach to evaluate iFC. To examine the specific networks that were influenced by homotopic connectivity, regions with altered VMHC among PDD, PDND, and HC groups were selected as regions of interest (ROIs) for seeds. Correlation functional analyses were performed by computing temporal correlation between each seed and the rest of brain in a voxel-wise manner. An entire brain z value map was created for each subject using the Fisher r-to-z transformation. A preceding united mask was used for whole-brain correction of multiple comparisons.

Statistical Analyses

Data were analyzed using SPSS software (version 19.0; IBM, Armonk, NY, United States). One-way analyses of variance (ANOVAs) and post hoc two-sample t test in a pair-wise manner were used to identify the differential brain regions among the three groups with gender, age, and education as covariates. The clusters showing significant differences (voxel-level p < 0.01, cluster size >15 voxels, corresponding to a corrected p < 0.05 as determined by AlphaSim correction) in VMHC among three groups were extracted. The mean VMHC z values of these clusters were calculated. The Pearson correlation between the mean VMHC z values and the neuropsychological measurements was then calculated. Finally, the clusters with significant differences were used as seeds to perform a secondary seed-based FC analysis. The clusters showing significant differences in FC between PDD and PDND groups were extracted. The mean FC z values of these clusters were calculated. The Pearson correlation between the mean z-FC values and the neuropsychological measurements was analyzed. The mean FC values within these ROIs were further used for receiver operating characteristic (ROC) curves analyses. A cutoff was selected using Youden index and minimum distance from the coordinate (0, 1) on the ROC curve methods. All tests were two-tailed, and p < 0.05 was considered statistically significant.

Results

Sample Characteristics

Demographic and clinical characteristics of the participants are presented in Table 1. No significant differences in age, gender, education level, MMSE scores, relative gray matter volume, relative white matter volume, and mean framewise displacement were observed between any two of the three groups. Gender, age, and education level were further regressed, but no significant differences were observed. The relative gray matter volume and relative white matter volume used in the data analysis were measured from the segmented T1 images with CAT12 (A Computational Anatomy Toolbox for SPM 12). However, the HDRS-17 and the BDI scores were significantly different among the three groups (p < 0.001). There were no significant differences in PD duration, UPDRS-III, and H&Y staging between the PDD and PDND patients.

Voxel-Mirrored Homotopic Connectivity Differences

Significant differences in the VMHC values of the bilateral precentral gyrus (PCG), medial frontal gyrus, and paracentral lobule (PCL) were observed among PDD, PDND, and HC groups (Figure 1A and Table 2). Voxel-mirrored homotopic connectivity in the PCG was significantly decreased in PDND compared with HCs, whereas VMHC in the PCG was notably decreased in PDD patients compared with HCs (Figure 1B and Table 2). Voxel-mirrored homotopic connectivity values in the bilateral medial frontal gyrus and PCL were significantly decreased in PDD compared with PDND and HCs, but there was no significant difference between PDND and HC groups (Figure 1C and Table 2).

FIGURE 1

Figure 1. VMHC differences among PDD, PDND, and HC groups. The ANCOVA thresholds were set at a voxel-level p < 0.01, cluster size >15 voxels, corresponding to a corrected p < 0.05 as determined by AlphaSim correction. (A) Representative images. The yellow/red color bar indicates the F value from the ANCOVAs. (B) VMHC values of the PCG. *p < 0.05, **p < 0.01 vs. HC. (C) VMHC values of the MFG/PCL. Bar plots representing the mean VMHC values (and standard deviation) between three groups. **p < 0.01 vs. HC, ^#p < 0.05 vs. PDND. VMHC, voxel-mirrored homotopic connectivity; PDD, Parkinson’s disease with depression; PDND, Parkinson’s disease without depression; HC, healthy control; ANCOVA, analysis of covariance; PCG, precentral gyrus; PCL, paracentral lobule.

TABLE 2

Table 2. VMHC differences among PDD, PDND, and HCs.

Secondary Seed-Based Homotopic Resting-State Functional Connectivity

To evaluate the effects of disrupted interhemispheric connectivity on the relevant resting-state networks, the clusters involved in the right and left bilateral PCG and the right and left medial frontal gyrus (MFG)/PCL were used as the four ROI (including all the voxels in these clusters, Figure 2A) to investigate patterns of iFC throughout the entire brain among three groups. With the use of the right PCG as seed 1, a decrease in iFC was observed between the right PCG and the right fusiform gyrus, the left inferior occipital gyrus, the left lingual gyrus, the right calcarine fissure and surrounding cortex, and the bilateral postcentral gyrus (PoCG) among three groups (Figure 2B and Table 3) (AlphaSim corrected, combined with a different cluster size of >56 voxels).

FIGURE 2

Figure 2. RSFC from the right/left PCG, and the right/left MFG/PCL seeds among PDD, PDND, and HC groups. (A) Seed regions with decreased voxel-mirrored homotopic connectivity. (B) RSFC patterns connected with the right PCG as the seed 1 from the entire brain among the three groups. (C) RSFC patterns connected with the left PCG as the seed 2 from the entire brain among the three groups. (D) RSFC patterns connected with the right MFG/PCL as the seed 3 from the entire brain among the three groups. (E) RSFC patterns connected with the left MFG/PCL as the seed 4 from the entire brain among the three groups. The color bars indicate the F values. RSFC, resting-state functional connectivity; PCG, precentral gyrus; PCL, paracentral lobule; PDD, Parkinson’s disease with depression; PDND, Parkinson’s disease without depression; HC, healthy control.

TABLE 3

Table 3. RSFC among PDD, PDND, and HC groups.

With the use of the left PCG as seed 2, a significant difference in iFC was observed between the left PCG and the left PoCG among three groups (AlphaSim corrected, combined with a different cluster size of >74 voxels). Both PDD and PDND patients showed a decreased connectivity between seed 2 and the left PoCG compared with HCs (Figure 2C and Table 3).

With the use of the right MFG/PCL as seed 3, a decreased iFC was observed between the right MFG/PCL and the right lingual gyrus, and the bilateral PoCG, whereas an increased iFC between the right MFG/PCL and the left supplementary motor area (SMA) was observed among PDD, PDND, and HC groups (AlphaSim corrected, combined with a different cluster size of >49 voxels) (Figure 2D and Table 3). A decreased functional connectivity in the right MFG/PCL and the left SMA was observed in PDD patients compared with PDND patients. An increased functional connectivity from the right MFG/PCL to the left SMA was observed in PDND patients compared with HC.

With the use of the left MFG/PCL as seed 4, a decreased iFC was observed between the seed 4 and the bilateral lingual gyrus, the bilateral PoCG, the bilateral precuneus, the left middle occipital gyrus, the left caudate, the left PCG, the right cuneus, and the right SMA (AlphaSim corrected, combined with a different cluster size of >10 voxels) (Figure 2E and Table 3). A significant decrease in functional connectivity between the left MFG/PCL and the right SMA was observed in PDD patients compared with PDND patients.

Relations Between the Voxel-Mirrored Homotopic Connectivity/Resting-State Functional Connectivity Values and Clinical Variables

PCG and MFG/PCL were selected as ROIs, because significant VMHC differences were observed among PDD, PDND, and HC groups. Linear correlation analysis was performed between the mean VMHC values within the ROIs and clinical variables in the PDD and PDND patients. No correlation was found between VMHC and illness duration, HDRS, BDI, and UPDRS-III in either PDD or PDND group. The clusters showing significant differences in seed-based FC between PDD and PDND groups were extracted for correlation analysis. The RSFC values of right MFG/PCL to the left SMA and RSFC values of the left MFG/PCL to the right SMA were then extracted. The decreased homotopic RSFC from the right MFG/PCL to the left SMA was negatively correlated with HDRS scores (r = –0.408, p = 0.012) in only PDD patients (Table 4).

TABLE 4

Table 4. Correlations between the VMHC/RSFC and clinical variables.

Discriminatory Analysis of the Resting-State Functional Connectivity From the Right Medial Frontal Gyrus/Paracentral Lobule to the Left Supplementary Motor Area and the Left Medial Frontal Gyrus/Paracentral Lobule to the Right Supplementary Motor Area

Because seed 3 (the right MFG/PCL to the left SMA) and seed 4 (the left MFG/PCL to the right SMA) exhibited significant RSFC differences between PDD and PDND groups, they were selected as ROIs. Mean RSFC values were extracted for further ROC analyses. The results showed that the area under the curve (AUC) of RSFC from seed 3-cluster 4 was 0.691 [95% confidence interval (CI): 0.569–0.814, p = 0.002] and the AUC of RSFC from seed 4-cluster 11 was 0.680 (95% CI: 0.563–0.798, p = 0.004). When separating PDD patients from PDND patients, the sensitivity and specificity for separating PDD from PDND were relatively high (Table 5).

TABLE 5

Table 5. ROC analyses for PDD from PDND.

Discussion

To our knowledge, this is the first study that used the VMHC method combined with the seed-based homotopic RSFC to explore interhemispheric functional coordination in PDD and PDND patients. Using homotopic RSFC to measure interhemispheric synchrony or homotopy during the resting state is more important than other functional connectivity analyses (Gee et al., 2011). This study found that PDD patients had a decreased VMHC in the bilateral MFG/PCL; a decreased homotopic RSFC from the MFG/PCL to the contralateral SMA when compared with PD patients without depression PDND; and a decreased VMHC in the bilateral PCG when compared with PD patients without depression PDND and to HCs. The decreased homotopic RSFC from the right MFG/PCL to the left SMA was negatively correlated with the HDRS scores. The decreased RSFC in the right MFG/PCL to the left SMA could distinguish PDD patients from PDND.

The medial frontal gyrus (MFG) is located in the medial prefrontal cortex (mPFC). MFG is thought to be involved in cognitive regulation and affective processing (Hänsel and von Känel, 2008), whereas mPFC is widely known to play an important role in cognitive and limbic functions (Liang et al., 2015) and can modulate the expression of reward-seeking behavior through regulating the reward-related neural circuit (Ferenczi et al., 2016; Belujon and Grace, 2017). The abnormal spontaneous neural activity in the MFG has been demonstrated to be an important factor for the development of depressive symptoms in PD patients (Cardoso et al., 2009), whereas the right MFG activation is proposed as a biomarker for the occurrence and the severity of depression in PD patients (Zhi et al., 2019). Previous studies also observed a decreased activation in the mPFC of PDD patients compared with PDND patients (Cardoso et al., 2009). In contrast, the PCL controls motor and sensory innervation and plays an important role in somatosensation (Jellinger, 2008). Voxel-mirrored homotopic connectivity reflects iFC. This study found that PDD had a decreased VMHC in the bilateral MFG/PCL and a decreased homotopic RSFC from the MFG/PCL to the SMA compared with PDND. Taken together, the reduced VMHC and RSFC within the bilateral MFG/PCL may impede the down-regulation of the negative emotion and cause a deficiency in the reward-related neural circuit leading to anhedonia in PDD patients.

Our study showed that the decreased homotopic RSFC from the right MFG/PCL to the left SMA was negatively correlated with HDRS scores (r = -0.41, p = 0.02), but not with illness duration, BDI, and UPDRS-III in PD patients; this suggests that aberrant interactions between the right MFG/PCL and the left SMA are related to the severity of depression in PD patients. The results of ROC curve analyses indicated that the RSFC from the right MFG/PCL to the left SMA may have the ability to identify the PDD patients from PDND, which may reduce misdiagnosis. Thus, the disrupted right MFG/PCL–left SMA homotopic RSFC is not only related to the severity of depression in PD, but it can also distinguish PDD from PDND patients.

The midline core of the default mode network (DMN) comprises the posterior cingulate cortex, adjacent precuneus, and the mPFC (Sampaio et al., 2014). Functionally, DMN is involved in major human brain functions including spontaneous cognition (Xu et al., 2014). The motor network (MN) regions consist of the bilateral sensorimotor area, SMA, and bilateral cerebellum. Parkinson’s disease is traditionally considered as the most common movement disorder characterized by cardinal motor symptoms (Luo et al., 2014). However, in this study, PDD and PDND patients showed significant differences in homotopic RSFC from the MFG/PCL to the contralateral SMA. It could be said that disrupted RSFC in the midline structure of DMN to the MN can be used to distinguish PDD from PDND. However, using the bilateral PCG and the bilateral MFG/PCL for four seeds, we found a decreased homotopic RSFC from the right PCG to the bilateral PoCG, the left PCG to the left PoCG, the right MFG/PCL to the bilateral PoCG and the left SMA, the left MFG/PCL to the bilateral PoCG, and the left PCG and the right SMA among PDD, PDND, and HC groups. Overall, our findings suggest that PDD might depend on the damage to specific neural networks rather than on the dysfunction of single, discrete brain region.

We acknowledge several limitations of this study. First, our sample size is relatively small with the clinically heterogeneous group. Studies with a larger and clinically homogenous sample are needed to further validate our findings. Second, the human brain is not structurally symmetrical. Therefore, we tried to resolve this issue by smoothing the functional data and normalizing them to a symmetric template (Zuo et al., 2010). Hence, morphometric asymmetry could not account for the reduced VMHC (Kelly et al., 2011). Third, the current study mainly focused on the differences of interhemispheric synchrony between PD with and without depression PDND, only three neuropsychological tests (MMSE, HDRS-17, and BDI-21) were used, and we did not observe anxiety or apathy among the PD patients. A broader spectrum of tests may be helpful to evaluate the emotional and cognitive functions of PD patients. Fourth, the patients in the current study were not drug naïve. Although the anti-PD or antidepressant medications were stopped for 12 h before imaging and neuropsychological testing, the potentially confounding effects of chronic medications could not be avoided. Fifth, the lack of a primary depression group left the question of whether PDD shares a common neurobiological substrate with that of primary depression unanswered. Sixth, the multiple-comparison correction is not strict. Seventh, the seed areas were defined in the same sample that further analysis was conducted in. Finally, this is a cross-sectional study. It remains unsolved whether the found alterations represent state or trait characteristics of our patient cohort.

In conclusion, the occurrence of depression in PD might arise from the disturbance in the MFG and PCL, and disrupted RSFC in the midline structure of DMN to the MN can be used to diagnose depression in PD patients.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Ethics Statement

The studies involving human participants were reviewed and approved by Medical Research Ethical Committee of The Second Xiangya Hospital, Central South University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

CT and XZ contributed to the conception and design of the study. HL, SC, CW, HZ, and MW contributed to data collection. HL, JF, and QS contributed to data analysis. HL contributed to writing the manuscript. JL contributed to the English-language revision.

Funding

Funding from the National Natural Science Foundation of China (No. 81471646) and the Science Foundation of Hunan Province (No. 14JJ2027) is gratefully acknowledged.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Footnotes

1. ^ http://www.fil.ion.ucl.ac.uk/spm/software/spm12/

References

Aarsland, D., Påhlhagen, S., Ballard, C. G., Ehrt, U., and Svenningsson, P. (2011). Depression in Parkinson disease–epidemiology, mechanisms and management. Nat. Rev. Neurol. 8, 35–47. doi: 10.1038/nrneurol.2011.189

PubMed Abstract | CrossRef Full Text | Google Scholar

Belujon, P., and Grace, A. A. (2017). Dopamine system dysregulation in major depressive disorders. Int. J. Neuropsychopharmacol. 20, 1036–1046. doi: 10.1093/ijnp/pyx056

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, E. F., Maia, F. M., Fregni, F., Myczkowski, M. L., Melo, L. M., Sato, J. R., et al. (2009). Depression in Parkinson’s disease: convergence from voxel-based morphometry and functional magnetic resonance imaging in the limbic thalamus. Neuroimage 47, 467–472. doi: 10.1016/j.neuroimage.2009.04.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Daniel, S. E., and Lees, A. J. (1993). Parkinson’s disease society brain bank, London: overview and research. J. Neural Transm. Suppl. 39, 165–172.

Google Scholar

Feldmann, A., Illes, Z., Kosztolanyi, P., Illes, E., Mike, A., Kover, F., et al. (2008). Morphometric changes of gray matter in Parkinson’s disease with depression: a voxel-based morphometry study. Mov. Disord. 23, 42–46. doi: 10.1002/mds.21765

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferenczi, E. A., Zalocusky, K. A., Liston, C., Grosenick, L., Warden, M. R., Amatya, D., et al. (2016). Prefrontal cortical regulation of brainwide circuit dynamics and reward-related behavior. Science 351:aac9698. doi: 10.1126/science.aac9698

PubMed Abstract | CrossRef Full Text | Google Scholar

Frosini, D., Unti, E., Guidoccio, F., Del Gamba, C., Puccini, G., Volterrani, D., et al. (2015). Mesolimbic dopaminergic dysfunction in Parkinson’s disease depression: evidence from a 123I-FP-CIT SPECT investigation. J. Neural Transm. 122, 1143–1147. doi: 10.1007/s00702-015-1370-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Gee, D. G., Biswal, B. B., Kelly, C., Stark, D. E., Margulies, D. S., Shehzad, Z., et al. (2011). Low frequency fluctuations reveal integrated and segregated processing among the cerebral hemispheres. Neuroimage 54, 517–527. doi: 10.1016/j.neuroimage.2010.05.073

PubMed Abstract | CrossRef Full Text | Google Scholar

Hänsel, A., and von Känel, R. (2008). The ventro-medial prefrontal cortex: a major link between the autonomic nervous system, regulation of emotion, and stress reactivity? Biol. Psychol. Soc. Med. 2:21.

Google Scholar

Hu, X., Song, X., Yuan, Y., Li, E., Liu, J., Liu, W., et al. (2015). Abnormal functional connectivity of the amygdala is associated with depression in Parkinson’s disease. Mov. Disord. 30, 238–244. doi: 10.1002/mds.26087

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Ravdin, L. D., Nirenberg, M. J., Piboolnurak, P., Severt, L., Maniscalco, J. S., et al. (2013). Neuroimaging markers of motor and nonmotor features of Parkinson’s disease: an 18f fluorodeoxyglucose positron emission computed tomography study. Dement. Geriatr. Cogn. Disord. 35, 183–196. doi: 10.1159/000345987

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, P., Xu, X., Gu, Q., Xuan, M., Yu, X., Luo, W., et al. (2014). Disrupted white matter integrity in depressed versus non-depressed Parkinson’s disease patients: a tract-based spatial statistics study. J. Neurol. Sci. 346, 145–148. doi: 10.1016/j.jns.2014.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Jellinger, K. A. (2008). Neuroscience in Medicine. New York, NY: Humana Press.

Google Scholar

Kelly, C., Zuo, X.-N., Gotimer, K., Cox, C. L., Lynch, L., Brock, D., et al. (2011). Reduced interhemispheric resting state functional connectivity in cocaine addiction. Biol. Psychiatry 69, 684–692. doi: 10.1016/j.biopsych.2010.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostić, V. S., Agosta, F., Petrović, I., Galantucci, S., Spica, V., Jecmenica-Lukic, M., et al. (2010). Regional patterns of brain tissue loss associated with depression in Parkinson disease. Neurology 75, 857–863. doi: 10.1212/wnl.0b013e3181f11c1d

PubMed Abstract | CrossRef Full Text | Google Scholar

Langston, J. W. (2006). The Parkinson’s complex: parkinsonism is just the tip of the iceberg. Ann. Neurol. 59, 591–596. doi: 10.1002/ana.20834

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Jia, J., and Yang, Z. (2016). Mini-mental state examination in elderly chinese: a population-based normative study. J. Alzheimers. Dis. 53, 487–496. doi: 10.3233/jad-160119

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Liu, J., Skidmore, F., Liu, Y., Tian, J., and Li, K. (2010). White matter microstructure changes in the thalamus in Parkinson disease with depression: a diffusion tensor MR imaging study. A.J.N.R. Am. J. Neuroradiol. 31, 1861–1866. doi: 10.3174/ajnr.a2195

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, Z., Watson, G. D., Alloway, K. D., Lee, G., Neuberger, T., and Zhang, N. (2015). Mapping the functional network of medial prefrontal cortex by combining optogenetics and fMRI in awake rats. Neuroimage 117, 114–123. doi: 10.1016/j.neuroimage.2015.05.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, C., Chen, Q., Song, W., Chen, K., Guo, X., Yang, J., et al. (2014). Resting-state fMRI study on drug-naive patients with Parkinson’s disease and with depression. J. Neurol. Neurosurg. Psychiatry 85, 675–683. doi: 10.1136/jnnp-2013-306237

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsui, H., Nishinaka, K., Oda, M., Niikawa, H., Komatsu, K., Kubori, T., et al. (2007). Depression in Parkinson’s disease. Diffusion tensor imaging study. J. Neurol. 254, 1170–1173. doi: 10.1007/s00415-006-0236-6

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald, W. M., Richard, I. H., and DeLong, M. R. (2003). Prevalence, etiology, and treatment of depression in Parkinson’s disease. Biol. Psychiatry 54, 363–375. doi: 10.1016/s0006-3223(03)00530-4

CrossRef Full Text | Google Scholar

Ossowska, K., and Lorenc-Koci, E. (2013). Depression in Parkinson’s disease. Parmacol. Rep. 65, 1545–1557.

Google Scholar

Pachana, N. A., Egan, S. J., Laidlaw, K., Dissanayaka, N., Byrne, G. J., Brockman, S., et al. (2013). Clinical issues in the treatment of anxiety and depression in older adults with Parkinson’s disease. Mov. Disord. 28, 1930–1934. doi: 10.1002/mds.25689

PubMed Abstract | CrossRef Full Text | Google Scholar

Postuma, R. B., Aarsland, D., Barone, P., Burn, D. J., Hawkes, C. H., Oertel, W., et al. (2012). Identifying prodromal Parkinson’s disease: pre-motor disorders in Parkinson’s disease. Mov. Disord. 27, 617–626. doi: 10.1002/mds.24996

PubMed Abstract | CrossRef Full Text | Google Scholar

Qamhawi, Z., Towey, D., Shah, B., Pagano, G., Seibyl, J., Marek, K., et al. (2015). Clinical correlates of raphe serotonergic dysfunction in early Parkinson’s disease. Brain 138(Pt 10), 2964–2973. doi: 10.1093/brain/awv215

PubMed Abstract | CrossRef Full Text | Google Scholar

Ravina, B., Camicioli, R., Como, P. G., Marsh, L., Jankovic, J., Weintraub, D., et al. (2007). The impact of depressive symptoms in early Parkinson disease. Neurology 69, 342–347. doi: 10.1212/01.wnl.0000268695.63392.10

CrossRef Full Text | Google Scholar

Rickards, H. (2005). Depression in neurological disorders: Parkinson’s disease, multiple sclerosis, and stroke. J. Neurol. Neurosur. Psychiatr. 76, I48–I52.

Google Scholar

Sampaio, A., Soares, J. M., Coutinho, J., Sousa, N., and Gonçalves, ÓF. (2014). The Big Five default brain: functional evidence. Brain Struct. Funct. 219, 1913–1922. doi: 10.1007/s00429-013-0610-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Dijk, K. R., Sabuncu, M. R., and Buckner, R. L. (2012). The influence of head motion on intrinsic functional connectivity MRI. Neuroimage 59, 431–438. doi: 10.1016/j.neuroimage.2011.07.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Vriend, C., Raijmakers, P., Veltman, D. J., van Dijk, K. D., van der Werf, Y. D., Foncke, E. M., et al. (2014). Depressive symptoms in Parkinson’s disease are related to reduced [123I]FP-CIT binding in the caudate nucleus. J. Neurol. Neurosurg. Psychiatry 85, 159–164. doi: 10.1136/jnnp-2012-304811

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Chen, H., Wu, J., Tao, L., Pang, Y., Gu, M., et al. (2018). Altered resting-state voxel-level whole-brain functional connectivity in depressed Parkinson’s disease. Parkinsonism Relat. Disord. 50, 74–80. doi: 10.1016/j.parkreldis.2018.02.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, X., Wu, X., Liu, J., Li, K., and Yao, L. (2013). Abnormal baseline brain activity in non-depressed Parkinson’s disease and depressed Parkinson’s disease: a resting-state functional magnetic resonance imaging study. PLoS One 8:e63691. doi: 10.1371/journal.pone.0063691

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Vik, A., Groote, I. R., Lagpopoulos, J., Holen, A., and Ellingsen, Ø, et al. (2014). Nondirective meditation activates default mode network and areas associated with memory retrieval and emotional processing. Front. Hum. Neurosci. 8:86. doi: 10.3389/fnhum.2014.00086

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C. G., Craddock, R. C., Zuo, X. N., Zang, Y. F., and Milham, M. P. (2013). Standardizing the intrinsic brain: towards robust measurement of inter-individual variation in 1000 functional connectomes. Neuroimage 80, 246–262. doi: 10.1016/j.neuroimage.2013.04.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C. G., Wang, X. D., Zuo, X. N., and Zang, Y. F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhi, Y., Yuan, Y., Si, Q., Wang, M., Shen, Y., Wang, L., et al. (2019). The association between DRD3 Ser9Gly polymorphism and depression severity in Parkinson’s disease. Parkinsons Dis. 2019:1642087.

Google Scholar

Zhou, F., Zhao, Y., Huang, M., Zeng, X., Wang, B., and Gong, H. (2018). Disrupted interhemispheric functional connectivity in chronic insomnia disorder: a resting-state fMRI study. Neuropsychiatr. Dis. Treat. 14, 1229–1240. doi: 10.2147/ndt.s162325

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Y., Song, X., Xu, M., Hu, X., Li, E., Liu, J., et al. (2016). Impaired interhemispheric synchrony in Parkinson’s disease with depression. Sci. Rep. 6:27477.

Google Scholar

Zuo, X. N., Kelly, C., Di Martino, A., Mennes, M., Margulies, D. S., Bangaru, S., et al. (2010). Growing together and growing apart: regional and sex differences in the lifespan developmental trajectories of functional homotopy. J. Neurosci. 30, 15034–15043. doi: 10.1523/jneurosci.2612-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar