Altered Functional Connectivity Density in Subtypes of Parkinson’s Disease

Introduction

Parkinson’s disease (PD) is a neurodegenerative disorder that is typically characterized by motor signs that include tremor, akinesia, and rigidity (Lees et al., 2009). It can be classified as the akinetic-rigid (AR) or tremor-dominant (TD) subtype based on the predominant motor signs (Zaidel et al., 2009). This symptomatic heterogeneity of the different PD subtypes may have unique pathophysiological substrates. Uncovering the neural substrates that account for this heterogeneity is critically important for advancing our knowledge of PD and to identify effective therapies. However, the brain mechanisms that underlie the specific PD subtypes are not well-established.

Previous functional neuroimaging studies have demonstrated that the neurobiological substrates of parkinsonian tremor differ from those of akinesia and rigidity. Previous studies have documented structural and/or functional changes in the two PD subtypes. A voxel-based morphometry study showed that TD patients had reduced gray matter volume in the cerebellum (Herman et al., 2013). Using R-fMRI, our research team compared TD and akinetic/rigid-predominant subtypes of PD and found distinct regional homogeneity patterns of spontaneous brain activity in the striato-thalamo-cortical loop and the cerebello-thalamo-cortical loop between these two subtypes of PD (Zhang et al., 2015c). However, these approaches have examined connectivity within a single network [e.g., the default mode network (DMN)] (Tan et al., 2015) or in local regions (e.g., the ReHo) (Zhang et al., 2015c) and have not considered alterations in whole-brain networks. Therefore, in the present study, we sought to fully delineate the abnormal functional connectivity (FC) networks of the two subtypes by searching the entire brain connectome at a refined voxel level. Functional connectivity density (FCD) is a recently developed, data-driven (i.e., with no need for a priori hypotheses) method for identifying the number of functional connections of each brain voxel (Tomasi and Volkow, 2010, 2012; Luo et al., 2014). High FCD values for particular voxels indicate that those voxels are functionally connected to a greater number of other brain voxels and play more important roles (Hub) than others in information processing (Tomasi and Volkow, 2010, 2011; Chen et al., 2015). Based on the neighboring relationships between brain voxels, FCD can be further divided into local and global FCD (Tomasi and Volkow, 2010). Local FCD of a voxel reflects local functional integration, and global FCD reflects functional integration across the whole brain. Data-driven global FCD is useful for exploratory studies and that it can guide seed-voxel correlation analyses when we do not have clear hypotheses for the seed regions (Tomasi and Volkow, 2010, 2012). Although the FCD analysis allows us to identify brain sites that exhibit abnormal FC in two subtypes, it does not provide insight into the locations to which these abnormal connections are linked. To trace these locations, seed-based FC analyses were performed.

In the current study, combined FCD and seed-based FC analyses were performed to fully characterize the abnormal brain networks in the two subtypes. Finally, we correlated the network changes with clinical characteristics of the patients and used receiver operating characteristic (ROC) analysis to evaluate the ability of the FCD measures to discriminate the PD subtypes. In a previously published study, we observed differences in interhemispheric coordination between PD subtypes (Zhang et al., 2015b). Based on prior studies (Lewis et al., 2011; Zhang et al., 2015b) we hypothesized that the functional integration in the cerebellum and motor-related cortical regions would be altered in the two subtypes.

Materials and Methods

Subjects

The present study was approved by the Medical Research Ethics Committee of the Third Military Medical University (Chongqing, China), and written informed consent was obtained from all subjects. Fifty right-handed PD patients [25 AR patients (M:9/F:16, Age(y): 65.12 ± 10.11)] and 25 TD patients [M:9/F:17, Age(y): 60.28 ± 11.14] were recruited from the movement disorders outpatient clinic of Southwest Hospital. Each PD subject was diagnosed by an experienced movement disorder specialist according to the United Kingdom Parkinson’s Disease Society Brain Bank criteria (Hughes et al., 1992). Importantly, all PD patients were screened to exclude patients with visual hallucination symptoms.

In the current study, only those patients who could tolerate long periods of immobility in the MRI scanner were included. Subsequently, motor performance was quantitatively and accurately assessed in detail by using the motor section (part III) of the Unified Parkinson’s Disease Rating Scale (UPDRS) (Fahn et al., 1987) and used to classify PD subjects by subtype. Specifically, the PD subjects were classified as either TD or AR based on the ratio of the mean tremor score (TD score) to the mean akinetic-rigidity score (AR score) (Eggers et al., 2011; Lewis et al., 2011). Just like the previously published study (Zhang et al., 2015b), for each patient, we computed a mean tremor score for 9 tremor items (right and left arm tremor from patient history, and lips or chin tremor, tremor in all 4 limbs, and activity in both arms or postural tremor, on examination) and a mean score for 5 AR items (falling, freezing, and walking difficulty, from patient history, and gait and postural instability, on exam-ination). Briefly, patients were classified as having TD-PD when the ratio of the mean tremor score to the mean AR score was ≥1.5 and as having AR-PD when this ratio was ≤1 (Eggers et al., 2011; Lewis et al., 2011). In addition, a Mini Mental State Examination (MMSE) score > 28 and a Montreal Cognitive Assessment (MoCA) score > 26 were required for all subjects to ensure that no patients met the criteria for dementia (Dubois et al., 2007). The patient exclusion criteria included secondary parkinsonism, atypical parkinsonian disease, and pre-existing neurological or psychiatric disorders (including seizures, aphasia, neglect, substantial sensory disturbances, severe depression, or claustrophobia). Meanwhile, patients with severe tremor which might disturb the MRI procedure were excluded.

All patients were substituted with dopaminergic treatment but none were treated with antidepressants. Three of 25 TD patients were also receiving low-dose anticholinergics (benzhexol, 0.5 mg/3 times daily).

In addition, 26 healthy controls (HCs) [M:9/F:17, Age(y): 59.84 ± 7.01] were recruited from among local individuals who volunteered to participate in scientific studies. All of the control subjects had normal neurological examination results and underwent neuropsychological testing prior to the MRI scan. The exclusion criteria were the same as those applied to the patients.

Data Acquisition

Functional images were acquired using a 3.0 T Siemens Tim Trio whole-body MRI system (Siemens Medical Solutions, Erlangen, Germany). Subjects were instructed to stay awake and close their eyes, and to try not to think of anything. All subjects confirmed they did not fall asleep during the MRI procedures. Foam padding and earplugs were used to reduce head motion and scanner noise. Imaging data were collected transversely by using an echo-planar imaging (EPI) sequence with the following settings: TR = 2000 ms, TE = 30 ms, flip angle = 90°, FOV = 192 mm × 192 mm, in-plane matrix = 64 × 64, thickness = 3 mm, voxel size = 3.0 mm × 3.0 mm × 3.0 mm. For each subject, a total of 240 volumes were acquired with scan time of 480 s. To minimize the impact of dopaminergic medication, the MRI scans were performed during a relatively hypodopaminergic state (12 h after a last dose of dopaminergic treatment). Other medications (anticholinergics) were withdrawn for 12 h prior to scanning (ALI–MELKKILÄ et al., 1993; Khor and Hsu, 2007).

Data Preprocessing

Data preprocessing was carried out using Neuroscience Information Toolbox (NIT¹). A series of conventional preprocessing steps were performed, which include (1) discarding first five volumes; (2) slice-timing correction; (3) head-motion correction; (4) normalization with an EPI template in the Montreal Neurological Institute (MNI) space; (5) linear detrending; (6) temporal filtering (band-pass 0.01–0.08 Hz); and (7) Regressing out six motion parameters, white matter and cerebrospinal fluid. Participants were excluded from further analyses if head motion exceeded 2.0 mm or 2.0°. We further evaluated the mean absolute displacement of each volume as compared with the previous volume (Zhang et al., 2015b). The average/standard deviation of frame-wise displacement (FD) for each group is present at Supplementary Table 1. All the three groups showed very slight head motion (mean FD < 0.1). In addition, there is no differences among the three groups (F = 1.302, p = 0.278, ANOVA). Considering that the censoring may lead to fewer degrees of freedom and introduce noisy (Power et al., 2015), we did not apply the volume censoring in the preprocessing (Supplementary Data).

Functional Connectivity Density Analysis

According to the approach introduced by Tomasi and Volkow (2010) we used custom-written software in Neuroscience Information Toolbox (NIT¹) to evaluate the individual voxel-wise FCD maps. The global FCD of a given voxel was defined as the number of functional connections between that voxel and all other voxels. The local FCD was obtained using a “growing” algorithm. In this algorithm, given a voxel χ0, an additional voxel χj was added to the list of neighbors of χ0 if it was adjacent to a voxel that was linked to χ0 by a continuous path of functionally connected voxels and the correlation coefficient between χ0 and χj was >0.6. This calculation was repeated for all voxels that were adjacent to the neighbors of χ0 in an iterative manner until no new neighbors could be added to the list. The local FCD (intraregional) of χ0 was defined as the number of elements in the list of neighbors. Therefore, the global FCD reflected the whole brain FC, whereas the local FCD reflected local FC. The threshold of the correlation coefficient is a key parameter in FCD analysis. Based on prior knowledge, the threshold was set at 0.6 (Supplementary Data and Figure 1; Tomasi and Volkow, 2010, 2011). For each subject, a global FCD map and a local FCD map were computed. For standardization purposes, individual FCD maps were divided by the mean FCD value. Prior to the group-level analyses, we applied 6-mm FWHM spatial smoothing to the normalized FCD maps.

Seed-Based Functional Connectivity Analysis

The brain sites identified by FCD group comparisons between the two subtypes were used as seeds in subsequent FC analyses. Specifically, a sphere ROI was generated for each of the brain sites, the center of which corresponded to the peak voxel (radius = 3 mm); the mean BOLD time series was then extracted and correlated to the time series of all of the voxels in the brain (Li et al., 2015). In this way, FC maps of each seed were produced; these maps were transformed by Fisher’s r-to-z transformation to improve normality.

Statistical Analysis

One-way analyses of variance and Chi-squared tests were used to analyze demographic characteristics among the three groups, and two-sample t-tests were used to perform post hoc analysis. All tests were two-tailed, and p < 0.05 was considered statistically significant. These analyses were performed with SPSS software (version 20.0; SPSS, Chicago, IL, United States).

To determine the effect of group, a voxel-based one-way analysis of covariance (ANCOVA) was used with age, gender, and education level as covariates in both the FCD and FC analyses using REST software. The mean FCD and FC values were extracted for post hoc analyses using the least significant difference (LSD) t-test. Multiple comparisons correction was performed using a height threshold (p < 0.005) of individual voxel and a cluster size based on Monte Carlo simulations (Ledberg et al., 1998), which corresponds to cluster-level p < 0.05 by AlphaSim correction. The AFNI AlphaSim program² was used to correct for multiple comparisons.

We evaluated the clinical data for normality and then performed Spearman rank correlation analyses to investigate the relationships between each of the clinical variables (duration, age of onset, UPDRS-III, H-Y scale score, TD score, and AR score) and the altered functional brain properties separately (the normalized FCD values of the ROIs and the average z-score of the altered FC). ROC curves were computed using the mean global FCD differences between the TD patients and the AR patients. The discriminatory performance of each index was evaluated separately. Additionally, the indices were combined using a binary logistic regression model, and the performance of this model was also assessed. In each case, the area under the ROC curve (AUC) and the 95% confidence interval (CI) were calculated using non-parametric methods. Spearman rank correlation analyses and ROC analysis were performed with SPSS software.

Results

Demographic and Clinical Data of Participants

There were no significant differences among the three groups in terms of gender, age, or educational level and no significant differences between the TD and AR patients in terms of the duration of disease, Hoehn and Yahr scale score or UPDRS-III score (Table 1).

TABLE 1

TABLE 1. Demographic information and clinical characteristics of study subjects.

One-way analysis of variance did not reveal any significant differences in head motion when measured by either the mean or maximum relative displacement [F(2,69), translational mean: F = 0.039, p = 0.961; translational max: F = 0.493, p = 0.613; rotational mean: F = 2.93, p = 0.061; rotational max: F = 1.012, p = 0.369]. In addition, there is no differences among three groups in term of FD [F = 1.302, p = 0.278, ANOVA, F(2,69)], indicating that the PD patients and controls exhibited similar head motion characteristics.

Functional Connectivity Density Analysis

Using the averaged local and global FCD maps, we found that those regions with high local and global FCD (Hubs) were mainly distributed bilaterally, with maximal magnitude in the posterior cingulate, occipital and prefrontal cortices (Figure 1A and Table 2).

FIGURE 1

FIGURE 1. Results of the FCD analyses. (A) Shows the spatial distribution of the average global functional connectivity density (FCD) for TD, AR, and HC. (B) Show result of the significant main effect of group on FCD (one-way ANCOVA, p < 0.005 corrected for multiple comparison at the cluster level, AlphSim corrected). CP, cerebellum posterior lobe; MTG, middle temporal gyrus; MFG, middle frontal gyrus; SFG, superior frontal gyrus.

TABLE 2

TABLE 2. The difference of global FCD among TDs, ARs, and HCs.

For comparison among three groups of global FCD, the brain regions that showed a significant main effect (one-way ANCOVA, p < 0.005, AlphaSim corrected) predominantly included prefrontal regions (superior/middle frontal gyrus), the temporal cortex (middle temporal gyrus), and the cerebellum (Inferior and posterior lobe) (Figure 1B and Table 2).

Post hoc analyses were applied to identify significant pairwise differences between three pairs of groups (post hoc two-sample t-test, p < 0.005, AlphaSim corrected). (1) TD vs. AR. The TD patients had significantly decreased global FCD in the left inferior frontal gyrus, right middle frontal gyrus and right superior frontal gyrus, and increased global FCD in the cerebellum anterior lobe relative to the AR patients (Figure 2 and Table 2). (2) TD vs. HC. Compared with HCs, the TD patients had significantly increased global FCD in the right cerebellum posterior lobe, left superior frontal gyrus and right fusiform gyrus (Figure 2 and Table 2). (3) AR vs. HC. Compared with the HCs, the AR patients exhibited significantly increased global FCD in the left posterior lobe of the cerebellum, the right middle temporal gyrus, the bilateral middle frontal gyrus, and the right inferior frontal gyrus (Figure 2 and Table 2).

FIGURE 2

FIGURE 2. Statistical group differences (post hoc analyses) of FCD. The statistical differences (t score) of post hoc analysis (post hoc two-sample t-test, p < 0.005 corrected for multiple comparison at the cluster level, AlphSim corrected). Color bar represents the t-value.

No significant differences in local FCD were found among groups. Thus, the four regions that exhibited significant differences in global FCD between TD and AR (the left cerebellum anterior lobe, right superior frontal gyrus, right middle frontal gyrus, and left inferior frontal gyrus) were defined as seeds in the following FC analysis.

We correlated the clinical variables with the normalized FCD values of the four regions just mentioned (the left cerebellum anterior lobe, right superior frontal gyrus, right middle frontal gyrus, and left inferior frontal gyrus) for the TD and AR patients separately. In the PD group (TD and AR patients), a significant negative correlation was identified between H-Y scale score and the FCD value of the cerebellum anterior lobe (Figure 3A).

FIGURE 3

FIGURE 3. (A) Show the correlations between regions with altered global FCD and clinical data in PD group (TDs and ARs) (spearman rank correlation analyses, p < 0.05). (B) Shows receiver operating characteristic (ROC) curve of the combination of the four FCD indexes for distinguishing PD subtypes patients.

Discriminatory Performance of the Functional Connectivity Density

The results of the ROC analysis indicated that the discriminatory performance of each FCD index was far lower than the inclusion of the four FCD indices (Table 3) and the inclusion of the four FCD indices in a binary logistic regression model had the highest power to discriminate the AR patients from the TD patients, with an AUC of 0.96 (95% CI = 0.91–1.00; Figure 3B and Table 3). At a cutoff of 0.23, the sensitivity and specificity were 95.2 and 80.9%, respectively (Table 3).

TABLE 3

TABLE 3. Discriminatory performance of the functional connectivity density.

Functional Connectivity Analysis of Regions of Different FCD

For the seed of the left cerebellum anterior lobe, TD patients showed significantly increased FCs in the left ventral and dorsal post-central gyrus, left superior frontal gyrus, right insula and orbital inferior frontal gyrus compared with AR patients (Figure 4 and Table 4).

FIGURE 4

FIGURE 4. Results of the functional connectivity (FC) of the left cerebellum. (A) Shows the main effect of group on FC [one-way ANCOVA, p < 0.005 corrected for multiple comparison at the cluster level, AlphSim corrected, F(2,69)]. Color bar represents the F-value. (B) Show the differences on FC of the left cerebellum between TD, AR, and HC (post hoc two-sample t-test, p < 0.005 corrected for multiple comparison at the cluster level, AlphSim corrected). Color bar represents the t-value.

TABLE 4

TABLE 4. The altered functional connections of seeds in the group level.

For the remaining seeds, no significant difference in FC was found between TD and AR patients (Table 4).

Discussion

Our findings obtained using a combination of FCD and seed-based FC analyses provide consistent evidence for that the network disorganization of the brains in the two PD subtypes were different. These findings have important implications for understanding the neural substrates that underlie these disparate manifestations of PD.

These results indicate that AR and TD subgroups both represent altered network disorganization. The global FCD abnormalities between the three groups involved the premotor cortex, supplementary motor area (SMA), prefrontal region (e.g., middle frontal gyrus and superior frontal gyrus) and the cerebellum. The premotor cortex plays an important role in the temporal organization of sequential movements, the selection of movements and the generation of motor sequences from memory that fit into a precise plan (Deiber et al., 1991; Halsband et al., 1993). The SMA is critical for the planning and initiation of movements, particularly internally generated, self-paced movement sequences (Jenkins et al., 2000; Tanji and Hoshi, 2001). The prefrontal cortex plays an important role in the learning of new motor sequences (Jueptner et al., 1997). The cerebellum receives information from the motor cortex, and cerebellar output influences various neuronal populations in the motor cortex (Manto et al., 2012). Therefore, we postulated that the FCD of these sites differed between the PD group (AR and TD patients) and the HC group due to the importance of these sites for motor function.

Intriguingly, we found that there was increased global FCD in the frontal lobe in the AR patients compared with the TD patients. Previous studies have demonstrated that these related regions are associated with the AR subtype (Hanakawa et al., 1999). The SMA and the pre-SMA also play a crucial role in the pathogenesis of bradykinesia in PD (Cunnington et al., 2002; Esposito et al., 2013). The dysfunction of the SMA has been recognized as a crucial reason for akinesia (Zhang et al., 2016). Thus, we speculate that increased global FCD in the frontal lobe in AR patients may be the manifestation of an attempt to compensate for the weaker information processing function due to the pathogenesis of akinetic-rigidity even before the onset of gait symptoms. In addition, we found that global FCD in left middle temporal gyrus was increased in AR patients compared with the HC group and a significant negative correlation between age of disease onset and the FCD value in middle frontal gyrus of the AR patients. Several studies have confirmed that patients with late-onset PD have a greater tendency for early development of bradykinesia and rigidity (Wickremaratchi et al., 2011; Thenganatt and Jankovic, 2014), so this functional integration change may suggest that the AR subtype shows more aggressive symptom, manifested by earlier and more rapid physical decline (Thenganatt and Jankovic, 2014). What’s more, it is generally believed that the bilateral prefrontal cortices as well as middle temporal gyri are crucial nodes of the DMN (Raichle et al., 2001; Tessitore et al., 2012), and the AR subtype shows a more rapid development of cognitive decline (Karunanayaka et al., 2016). This result suggests that there is an early functional disruption of the DMN in AR prior to clinical evidence of cognitive impairment, and the AR patients were likely to invest more neural resources to compensate for the cognitive decline (Marras et al., 2002; Rajput et al., 2009).

We also found increased global FCD in the cerebellum of TD patients compared with the AR patients. Previous studies have shown that structural and functional changes in the cerebellum are highly related to tremor-related diseases, such as cortical myoclonic tremor (Buijink et al., 2013) essential tremor (Louis et al., 2006; Popa et al., 2011) and resting tremor in PD (Mure et al., 2011) these changes may also be related to the symptoms. Mounting evidence suggests that cerebellar glucose metabolism and blood flow are associated with the severity of tremor in PD (Fukuda et al., 2004; Benninger et al., 2009), especially the correlations of changes in tremor frequency within condition with rCBF in the contralateral cerebellum has been observed (Fukuda et al., 2004). Wang’s research showed that the enhanced Vim–cerebellum/dentate nucleus connectivity is associated with parkinsonian tremor; this connectivity is more strengthened as tremor becomes more severe (Wang et al., 2016). We also found that alterations in interhemispheric correlations in the cerebellum may be involved in the neuronal basis of the resting tremor symptoms of PD (Zhang et al., 2015b). Furthermore, we found a significant negative correlation between H-Y stage and the FCD value of the cerebellum anterior lobe in the PD group. We suspect that in TD patients, abnormal elevations of blood perfusion and greater mobilization of the cerebellum may constitute a potential compensatory mechanism for the cerebellar dysfunction of PD. One likely explanation is that this phenomenon presents a functional compensation for the defective basal ganglia (Wu and Hallett, 2005). Previous study showed that increased coupling between the subthalamic nucleus (STN) and cerebellum might underlie the neural substrate of PD tremors. Lower FC between the STN and putamen might underpin PD gait and posture disturbances, while higher FC between the STN and visual cortex might play a compensatory role (Zhang et al., 2016). But we find no significant difference in the basal ganglia and substantia nigra between two subtypes, just like our previously published study (Supplementary Data and Figure 2).

Previously, Zhang et al. (2015b) demonstrated abnormal FCD in PD and speculated that the observed FCD reductions in the ventral visual pathway might be responsible for the non-psychotic visual hallucinations in PD patients. They also speculated that higher FCD in the posterior cingulate cortex (PCC) and precuneus might be associated with the inability of PD patients to redirect attentional resources. To identify FC alteration related to PD motor symptoms and subtypes, we strictly controlled the inclusion criteria, and patients with visual hallucinations or cognitive dysfunction were carefully excluded. We found no altered FCD in the ventral visual pathway, PCC or precuneus. However, compared with controls, we did find increased global FCD in the fusiform gyrus in the AR patients. Because the fusiform gyrus is involved in visual information processing and recognition, such as face recognition and face hallucinations (Benninger et al., 2009; Blom, 2009) we speculate that non-psychotic visual hallucinations could affect FCD in PD earlier than we thought, even before the appearance of these symptoms.

We further investigated FC in the four altered FCD regions. FC analysis also showed altered FC between the frontal lobe and cerebellum in the TD patients compared with the AR patients, consistent with the notion that PD tremor is most likely related to combined impairment of the cerebello-cortical circuit. Considering these results together with the abnormal FCD that was observed in the two subtypes and the relationships with clinical data, we speculate that as part of the compensatory effects required for the generation of movements, more cerebellar resources are recruited in the TD subtype, whereas more frontal lobe resources are recruited in the AR subtype. Further research is needed to further explore the mechanisms underlying these changes.

We also found that the FCD provided good discrimination between the AR and TD patients. The diagnosis of different subtypes depended on UPDRS score (Jankovic et al., 1990; Stebbins et al., 2013) and we achieved relatively objective discrimination of the different subtypes, with good sensitivity (95%) and specificity (80%) for classifying PD patients. Previous studies have only attempted to use FC metrics to distinguish between PD and healthy individuals (Fernandez-Seara et al., 2015; Zhang et al., 2015a). The results demonstrate the importance of altered global FCD in discriminating the two subtypes. As we relied on an experienced neurologist for the diagnosis of the two subtypes in this study and as these findings have not been reported previously, the classification results presented here are preliminary, and further studies are required to corroborate our interpretation. The diagnostic power of FCD for PD subtypes also requires further investigation.

Limitations

Some methodological limitations of this study should be mentioned. First, the small sample size and the relatively weak correction strategy employed (AlphaSim program) limited the conclusions that could be drawn; thus, a large number of participants should be included in future studies. Second, most of the PD patients were administered dopaminergic medications and a few with anticholinergics which might affect brain condition (Ehrt et al., 2010). Studies of drug-naïve individuals are warranted to exclude the effects of dopaminergic and anticholinergics medications on the findings. Longitudinal data are also required to further study the neuroimaging changes observed in the two subtypes. Third, the threshold (Tc) of the correlation coefficient used to define the functional connection was set at 0.6 in the FCD analysis based on prior knowledge, and this fixed value may lead to some false positive or negative findings. Flexible thresholds, such as a set of continuous thresholds, might be an alternative option for overcoming this limitation. We also evaluated the difference in FCD among groups using five thresholds (from 0.3 to 0.7 at intervals of 0.1), and similar findings were observed across several calculations using different thresholds (Supplementary Data and Figure 1). Lastly, considering other literature on PD pathophysiology, there are certainly other ways to subtype PD patients, i.e., according to different genetic backgrounds and/or clinical presentations. Our study, however, focused on TD versus AR subtype, one of the first clinical sub-types described in Jankovic et al. (1990). Future studies, therefore, are warranted to understand the relevance of resting state FC differences between PD subtypes other than TD versus AR.

Author Contributions

XH: project execution; MRI data acquisition; healthy control recruitment and demographic information acquisition; writing of manuscript. YJ: MRI data analysis, statistical analysis; writing of manuscript; revision of this manuscript. XJ: MRI data analysis. JZ: manuscript review and revision. ML: MRI data acquisition. JL: participant recruitment and clinical data acquisition. YZ: participant recruitment and clinical data acquisition. DY: manuscript review and revision. CL: project conception, manuscript review and critique. JW: project conception, manuscript review and critique.

Funding

This study was funded by grants from the National Nature Science Foundation of China (81701664, 81330032, 81471638), the Technology Innovation Program in Southwest Hospital (SWH2016JCYB-30), the PCSIRT Project (IRT0910), and the Special-Funded Program on National Key Scientific Instruments and Equipment Development of China (2013YQ49085908).

Conflict of Interest Statement

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnhum.2017.00458/full#supplementary-material

Footnotes

1. ^http://www.neuro.uestc.edu.cn/NIT.html
2. ^https://afni.nimh.nih.gov/afni/doc/manual/AlphaSim.pdf

References

ALI–MELKKILÄ, T., Kanto, J., and Iisalo, E. (1993). Pharmacokinetics and related pharmacodynamics of anticholinergic drugs. Acta Anaesthesiol. Scand. 37, 633–642. doi: 10.1111/j.1399-6576.1993.tb03780.x

CrossRef Full Text | Google Scholar

Benninger, D. H., Thees, S., Kollias, S. S., Bassetti, C. L., and Waldvogel, D. (2009). Morphological differences in Parkinson’s disease with and without rest tremor. J. Neurol. 256, 256–263. doi: 10.1007/s00415-009-0092-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Blom, J. D. (2009). A Dictionary of Hallucinations. New York, NY: Springer.

Google Scholar

Buijink, A. G., Caan, M. A., Tijssen, M. J., Hoogduin, J., Maurits, N., and van Rootselaar, A.-F. (2013). Decreased cerebellar fiber density in cortical myoclonic tremor but not in essential tremor. Cerebellum 12, 199–204. doi: 10.1007/s12311-012-0414-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Duan, M., Xie, Q., Lai, Y., Dong, L., Cao, W., et al. (2015). Functional disconnection between the visual cortex and the sensorimotor cortex suggests a potential mechanism for self-disorder in schizophrenia. Schizophr. Res. 166, 151–157. doi: 10.1016/j.schres.2015.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Cunnington, R., Windischberger, C., Deecke, L., and Moser, E. (2002). The preparation and execution of self-initiated and externally-triggered movement: a study of event-related fMRI. Neuroimage 15, 373–385. doi: 10.1006/nimg.2001.0976

PubMed Abstract | CrossRef Full Text | Google Scholar

Deiber, M. P., Passingham, R. E., Colebatch, J. G., Friston, K. J., Nixon, P. D., and Frackowiak, R. S. (1991). Cortical areas and the selection of movement: a study with positron emission tomography. Exp. Brain Res. 84, 393–402. doi: 10.1007/BF00231461

CrossRef Full Text | Google Scholar

Dubois, B., Burn, D., Goetz, C., Aarsland, D., Brown, R. G., Broe, G. A., et al. (2007). Diagnostic procedures for Parkinson’s disease dementia: recommendations from the movement disorder society task force. Mov. Disord. 22, 2314–2324. doi: 10.1002/mds.21844

PubMed Abstract | CrossRef Full Text | Google Scholar

Eggers, C., Kahraman, D., Fink, G. R., Schmidt, M., and Timmermann, L. (2011). Akinetic-rigid and tremor-dominant Parkinson’s disease patients show different patterns of FP-CIT single photon emission computed tomography. Mov. Disord. 26, 416–423. doi: 10.1002/mds.23468

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrt, U., Broich, K., Larsen, J. P., Ballard, C., and Aarsland, D. (2010). Use of drugs with anticholinergic effect and impact on cognition in Parkinson’s disease: a cohort study. J. Neurol. Neurosurg. Psychiatry 81, 160–165. doi: 10.1136/jnnp.2009.186239

PubMed Abstract | CrossRef Full Text | Google Scholar

Esposito, F., Tessitore, A., Giordano, A., De Micco, R., Paccone, A., Conforti, R., et al. (2013). Rhythm-specific modulation of the sensorimotor network in drug-naïve patients with Parkinson’s disease by levodopa. Brain 136, 710–725. doi: 10.1093/brain/awt007

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahn, S., Elton, R., and motUDC (1987). The Unified Parkinson’s Disease Rating Scale. Florham Park, NJ: McMillan Healthcare Information.

Google Scholar

Fernandez-Seara, M. A., Mengual, E., Vidorreta, M., Castellanos, G., Irigoyen, J., Erro, E., et al. (2015). Resting state functional connectivity of the subthalamic nucleus in Parkinson’s disease assessed using arterial spin-labeled perfusion fMRI. Hum. Brain Mapp. 36, 1937–1950. doi: 10.1002/hbm.22747

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukuda, M., Barnes, A., Simon, E. S., Holmes, A., Dhawan, V., Giladi, N., et al. (2004). Thalamic stimulation for parkinsonian tremor: correlation between regional cerebral blood flow and physiological tremor characteristics. Neuroimage 21, 608–615. doi: 10.1016/j.neuroimage.2003.09.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Halsband, U., Ito, N., Tanji, J., and Freund, H. J. (1993). The role of premotor cortex and the supplementary motor area in the temporal control of movement in man. Brain 116(Pt 1), 243–266. doi: 10.1093/brain/116.1.243

CrossRef Full Text | Google Scholar

Hanakawa, T., Fukuyama, H., Katsumi, Y., Honda, M., and Shibasaki, H. (1999). Enhanced lateral premotor activity during paradoxical gait in Parkinson’s disease. Ann. Neurol. 45, 329–336. doi: 10.1002/1531-8249(199903)45:3<329::AID-ANA8<3.0.CO;2-S

CrossRef Full Text | Google Scholar

Herman, T., Giladi, N., and Hausdorff, J. (2013). Neuroimaging as a window into gait disturbances and freezing of gait in patients with Parkinson’s disease. Curr. Neurol. Neurosci. Rep. 13, 1–13. doi: 10.1007/s11910-013-0411-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, A. J., Daniel, S. E., Kilford, L., and Lees, A. J. (1992). Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatry 55, 181–184. doi: 10.1136/jnnp.55.3.181

CrossRef Full Text | Google Scholar

Jankovic, J., McDermott, M., Carter, J., Gauthier, S., Goetz, C., Golbe, L., et al. (1990). Variable expression of Parkinson’s disease: a base-line analysis of the DATATOP cohort. The Parkinson study group. Neurology 40, 1529–1529. doi: 10.1212/WNL.40.10.1529

CrossRef Full Text | Google Scholar

Jenkins, I. H., Jahanshahi, M., Jueptner, M., Passingham, R. E., and Brooks, D. J. (2000). Self-initiated versus externally triggered movements. II. The effect of movement predictability on regional cerebral blood flow. Brain 123(Pt 6), 1216–1228. doi: 10.1093/brain/123.6.1216

PubMed Abstract | CrossRef Full Text | Google Scholar

Jueptner, M., Stephan, K. M., Frith, C. D., Brooks, D. J., Frackowiak, R. S., and Passingham, R. E. (1997). Anatomy of motor learning. I. Frontal cortex and attention to action. J. Neurophysiol. 77, 1313–1324.

Google Scholar

Karunanayaka, P. R., Lee, E.-Y., Lewis, M. M., Sen, S., Eslinger, P. J., Yang, Q. X., et al. (2016). Default mode network differences between rigidity-and tremor-predominant Parkinson’s disease. Cortex 81, 239–250. doi: 10.1016/j.cortex.2016.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Khor, S.-P., and Hsu, A. (2007). The pharmacokinetics and pharmacodynamics of levodopa in the treatment of Parkinson’s disease. Curr. Clin. Pharmacol. 2, 234–243. doi: 10.2174/157488407781668802

CrossRef Full Text | Google Scholar

Ledberg, A., Akerman, S., and Roland, P. E. (1998). Estimation of the probabilities of 3D clusters in functional brain images. Neuroimage 8, 113–128. doi: 10.1006/nimg.1998.0336

PubMed Abstract | CrossRef Full Text | Google Scholar

Lees, A. J., Hardy, J., and Revesz, T. (2009). Parkinson’s disease. Lancet 373, 2055–2066. doi: 10.1016/S0140-6736(09)60492-X

CrossRef Full Text | Google Scholar

Lewis, M. M., Du, G., Sen, S., Kawaguchi, A., Truong, Y., Lee, S., et al. (2011). Differential involvement of striato- and cerebello-thalamo-cortical pathways in tremor- and akinetic/rigid-predominant Parkinson’s disease. Neuroscience 177, 230–239. doi: 10.1016/j.neuroscience.2010.12.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., He, H., Huang, M., Zhang, X., Lu, J., Lai, Y., et al. (2015). Identifying enhanced cortico-basal ganglia loops associated with prolonged dance training. Sci. Rep. 5:10271. doi: 10.1038/srep10271

PubMed Abstract | CrossRef Full Text | Google Scholar

Louis, E. D., Vonsattel, J. P., Honig, L. S., Lawton, A., Moskowitz, C., Ford, B., et al. (2006). Essential tremor associated with pathologic changes in the cerebellum. Arch. Neurol. 63, 1189–1193. doi: 10.1001/archneur.63.8.1189

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, C., Tu, S., Peng, Y., Gao, S., Li, J., Dong, L., et al. (2014). Long-term effects of musical training and functional plasticity in salience system. Neural Plast. 2014:180138. doi: 10.1155/2014/180138

PubMed Abstract | CrossRef Full Text | Google Scholar

Manto, M., Bower, J. M., Conforto, A. B., Delgado-Garcia, J. M., da Guarda, S. N., Gerwig, M., et al. (2012). Consensus paper: roles of the cerebellum in motor control–the diversity of ideas on cerebellar involvement in movement. Cerebellum 11, 457–487. doi: 10.1007/s12311-011-0331-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Marras, C., Rochon, P., and Lang, A. E. (2002). Predicting motor decline and disability in Parkinson disease: a systematic review. Arch. Neurol. 59, 1724–1728. doi: 10.1001/archneur.59.11.1724

CrossRef Full Text | Google Scholar

Mure, H., Hirano, S., Tang, C. C., Isaias, I. U., Antonini, A., Ma, Y., et al. (2011). Parkinson’s disease tremor-related metabolic network: characterization, progression, and treatment effects. Neuroimage 54, 1244–1253. doi: 10.1016/j.neuroimage.2010.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Popa, T., Gallea, C., Marais, L., Roze, E., Vidailhet, M., Lehericy, S., et al. (2011). Abnormal cerebello-cortical connectivity in essential tremor is corrected by non invasive repeated stimulation of the cerebellum. Mov. Disord. 26, S250–S250. doi: 10.1016/S1388-2457(11)60296-4

CrossRef Full Text | Google Scholar

Power, J. D., Schlaggar, B. L., and Petersen, S. E. (2015). Recent progress and outstanding issues in motion correction in resting state fMRI. Neuroimage 105, 536–551. doi: 10.1016/j.neuroimage.2014.10.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U.S.A. 98, 676–682. doi: 10.1073/pnas.98.2.676

PubMed Abstract | CrossRef Full Text | Google Scholar

Rajput, A. H., Voll, A., Rajput, M. L., Robinson, C. A., and Rajput, A. (2009). Course in Parkinson disease subtypes: a 39-year clinicopathologic study. Neurology 73, 206–212. doi: 10.1212/WNL.0b013e3181ae7af1

PubMed Abstract | CrossRef Full Text | Google Scholar

Stebbins, G. T., Goetz, C. G., Burn, D. J., Jankovic, J., Khoo, T. K., and Tilley, B. C. (2013). How to identify tremor dominant and postural instability/gait difficulty groups with the movement disorder society unified Parkinson’s disease rating scale: comparison with the unified Parkinson’s disease rating scale. Mov. Disord. 28, 668–670. doi: 10.1002/mds.25383

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, Y., Tan, J., Deng, J., Cui, W., He, H., Yang, F., et al. (2015). Alteration of basal ganglia and right frontoparietal network in early drug-naïve Parkinson’s disease during heat pain stimuli and resting state. Front. Hum. Neurosci. 9:467. doi: 10.3389/fnhum.2015.00467

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanji, J., and Hoshi, E. (2001). Behavioral planning in the prefrontal cortex. Curr. Opin. Neurobiol. 11, 164–170. doi: 10.1016/S0959-4388(00)00192-6

CrossRef Full Text | Google Scholar

Tessitore, A., Esposito, F., Vitale, C., Santangelo, G., Amboni, M., Russo, A., et al. (2012). Default-mode network connectivity in cognitively unimpaired patients with Parkinson disease. Neurology 79, 2226–2232. doi: 10.1212/WNL.0b013e31827689d6

PubMed Abstract | CrossRef Full Text | Google Scholar

Thenganatt, M. A., and Jankovic, J. (2014). Parkinson disease subtypes. JAMA Neurol. 71, 499–504. doi: 10.1001/jamaneurol.2013.6233

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomasi, D., and Volkow, N. D. (2010). Functional connectivity density mapping. Proc. Natl. Acad. Sci. U.S.A. 107, 9885–9890. doi: 10.1073/pnas.1001414107

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomasi, D., and Volkow, N. D. (2011). Functional connectivity hubs in the human brain. Neuroimage 57, 908–917. doi: 10.1016/j.neuroimage.2011.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomasi, D., and Volkow, N. D. (2012). Aging and functional brain networks. Mol. Psychiatry 471, 549–558. doi: 10.1038/mp.2011.81

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Z., Chen, H., Ma, H., Ma, L., Wu, T., and Feng, T. (2016). Resting-state functional connectivity of subthalamic nucleus in different Parkinson’s disease phenotypes. J. Neurol. Sci. 371, 137–147. doi: 10.1016/j.jns.2016.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickremaratchi, M. M., Knipe, M. D. W., Sastry, B. S. D., Morgan, E., Jones, A., Salmon, R., et al. (2011). The motor phenotype of Parkinson’s disease in relation to age at onset. Mov. Disord. 26, 457–463. doi: 10.1002/mds.23469

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., and Hallett, M. (2005). A functional MRI study of automatic movements in patients with Parkinson’s disease. Brain 128, 2250–2259. doi: 10.1093/brain/awh569

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaidel, A., Arkadir, D., Israel, Z., and Bergman, H. (2009). Akineto-rigid vs. tremor syndromes in Parkinsonism. Curr. Opin. Neurol. 22, 387–393. doi: 10.1097/WCO.0b013e32832d9d67

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Liu, X., Chen, J., Liu, B., and Wang, J. (2015a). Widespread increase of functional connectivity in Parkinson’s disease with tremor: a resting-state FMRI study. Front. Aging Neurosci. 7:6. doi: 10.3389/fnagi.2015.00006

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Bi, W., Zhang, Y., Zhu, M., Zhang, Y., Feng, H., et al. (2015b). Abnormal functional connectivity density in Parkinson’s disease. Behav. Brain Res. 280, 113–118. doi: 10.1016/j.bbr.2014.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Wei, L., Hu, X., Xie, B., Zhang, Y., Wu, G.-R., et al. (2015c). Akinetic-rigid and tremor-dominant Parkinson’s disease patients show different patterns of intrinsic brain activity. Parkinsonism Relat. Disord. 21, 23–30. doi: 10.1016/j.parkreldis.2014.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. R., Feng, T., Hou, Y. N., Chan, P., and Wu, T. (2016). Functional connectivity of vim nucleus in tremor-and akinetic-/rigid-dominant Parkinson’s disease. CNS Neurosci. Ther. 22, 378–386. doi: 10.1111/cns.12512

PubMed Abstract | CrossRef Full Text | Google Scholar