Distinct Role of Striatal Functional Connectivity and Dopaminergic Loss in Parkinson’s Symptoms

Degeneration of dopaminergic neurons is a hallmark of Parkinson’s disease. However, its link to Parkinson’s disease symptoms remains unclear. Striatal resting state functional connectivity differentiates between Parkinson’s disease patients and healthy controls and might be a potential mediator of the effects of striatal dopaminergic degeneration onto Parkinson’s disease symptoms. Here, we evaluated the relationship between dopaminergic deficits, striatal functional connectivity (SFC) at rest and different Parkinson’s disease clinical symptoms in the largest currently established cohort of de novo Parkinson’s disease patients. We show that SFC is an independent predictor of symptom severity in Parkinson’s disease in addition to striatal dopaminergic deficits. Furthermore, we find that distinct SFC networks are associated with symptoms reflecting the ability to perform daily routine automatized motor tasks and clinician-rated Parkinson’s disease motor symptoms. We find that reduced SFC is a major and independent predictor of Parkinson’s disease symptoms going beyond the mere reflection of striatal dopaminergic input loss. These findings indicate the high value of SFC as a clinically relevant biomarker in Parkinson’s disease.

Introduction

Parkinson’s disease (PD) is a progressive neurodegenerative disorder that is associated with loss of dopaminergic projections from substantia nigra to the striatum. This decline of dopamine signaling initiates a pathological cascade encompassing numerous brain regions (Lotharius and Brundin, 2002; Dauer and Przedborski, 2003). Though loss of dopaminergic neurons as measured for example by dopamine transporter single photon emission tomography (DAT-SPECT) is consistently found to be correlated with PD motor symptoms at advanced disease stages, data on its association to symptom severity in de novo PD is rather limited (Pirker, 2003; McGhee et al., 2013). Also the contribution of different pathology levels and brain anatomical networks to the observed motor and non-motor symptoms covered by the Unified Parkinson’s Disease Rating Scale (UPDRS) and its revised version (MDS-UPDRS) – the most established PD rating scale – remains rather unclear (Fahn and Elton, 1987; Ramaker et al., 2002; Goetz et al., 2008). Subscales of the original and revised UPDRS scale have been repeatedly shown to reflect differential clinical aspects observed in PD patients (Martinez-Martin et al., 1994; Stebbins and Goetz, 1998; Goetz et al., 2008). UPDRS part I comprises symptoms related to mood and physiology – functions that are rather associated with prefrontal, insular and brain stem circuitry (Murai et al., 2001; Wolters, 2009; Benoit and Robert, 2011). UPDRS part II covers self-evaluation of highly automatized motor aspects of daily living. Retention of such automatized motor functions is closely associated with a striato-cerebellar circuitry (Doyon et al., 1998; Lang and Bastian, 2002). In contrast, UPDRS part III provides a direct evaluation of current motor deficits by the clinician. For this subscale, a stronger involvement of striatal, motor, premotor and prefrontal regions involved in execution and planning is therefore expected (Goldberg, 1985; Doyon et al., 1998; Petrides, 2005). A common mechanism underlying all symptoms captured by UPDRS subscales appears therefore unlikely.

The current knowledge regarding these potentially differential mechanisms underlying PD symptoms is rather limited. Evidence from single photon and positron emission tomography studies indicates that striatal dopaminergic signal accounts for up to 25% of variance in UPDRS total and its motor related part III subscale in advanced PD [for a detailed review see (McGhee et al., 2013)]. In contrast, the evidence for a link to UPDRS part I and II remains sparse. More recently, resting state functional magnetic resonance imaging (rsMRI) studies in PD reported significant associations between UPDRS part III and different striatal functional connectivity (SFC) metrics, explaining a similar proportion of variance as found for dopaminergic tracers (Wu et al., 2011; Hacker et al., 2012). At the same time several studies found a correlation between striatal different striatal connectivity metrics and striatal DAT binding (Lebedev et al., 2014; Rieckmann et al., 2015). These findings suggest that SFC desynchronization could be a further more downstream key pathological mechanism contributing to the observed PD motor symptoms (Ham et al., 2015). However, these studies did not evaluate if SFC provides any additional contribution to PD symptoms when accounting for loss of dopaminergic neurons. Correspondingly, the observed correlations could be due to the dopaminergic signal loss inducing both reductions in SFC and increases in symptom severity without a direct contribution of SFC to PD symptoms. The generalizability of these findings is also rather limited due to small sample sizes, inclusion of mostly advanced PD patients and insufficient control for potential treatment and atrophy effects onto the extracted resting state measures and the unclear link to other PD symptom domains.

Here we aim to address the question of differential contribution of dopaminergic deficits and SFC to PD clinical symptoms in a de novo PD population. We evaluate associations between dopaminergic loss, SFC and clinical symptoms as measured by the revised UPDRS (further referred to as UPDRS). Based on the above arguments, we hypothesize that reductions in SFC contribute to PD symptoms beyond the mere loss of dopaminergic neurons in the striatum. We further hypothesize that differential SFC networks indicated above are associated with specific clusters of symptoms observed in PD.

Materials and Methods

Subjects and Imaging Data

We extracted rsMRI data of de novo PD patients (n = 87) from the Parkinson’s Disease Progression Marker Initiative (PPMI)¹ database as available on 14th January, 2015. For 75 PD patients DAT-SPECT striatal binding ratios were available for the same visits as rsMRI. These data were used for cross-modality correlations described below. A detailed description of the PD cohort used for functional connectivity analyses is provided in Table 1. Additionally, we downloaded the first available (screening or baseline) structural MRI (sMRI) and DAT-SPECT scans for all available PD patients [sMRI: n = 148 (55 female, mean age ± standard deviation: 61.4 ± 9.4); DAT-SPECT: n = 419 (132 female, mean age ± standard deviation: 62.3 ± 9.4)] and age and sex matched healthy controls [sMRI: n = 69 (24 female, mean age ± standard deviation: 60.2 ± 11.2); DAT-SPECT: n = 198 (63 female, mean age ± standard deviation: 61.4 ± 10.8)]. As resting state was added after the study start some of the PD patients were scanned on PD medication. For some PD patients receiving drug treatment double evaluations of UPDRS part III both in the on and off medication state were available. As all other UPDRS subscales for these patients were acquired in the on state we restricted our analyses to the on medication scores. Written informed consent was obtained from all subjects. The study was approved by Institutional Review Boards/Independent Ethics Committees. For more details on the study please see http://www.ppmi-info.org/wp-content/uploads/2013/02/PPMI-Protocol-AM5-Final-27Nov2012v6-2.pdf.

TABLE 1

TABLE 1. Subject group characteristics for resting state SFC analyses.

Image Processing

All image pre-processing steps were performed using the Statistical Parametric Mapping 12 software package (SPM12)² and Matlab R2013.b (MathWorks). In brief, for rsMRI it comprised motion correction, for sMRI data segmentation using NewSegment, spatial normalization to the Montreal Neurological Institute (MNI) space using sMRI based parameters (without modulation for rsMRI to preserve the original signal and with modulation for gray matter probability maps to preserve volumetric information – GMV). Preprocessing of DAT-SPECT data comprised normalization to an average size DAT-SPECT template with subsequent normalization into MNI space. All imaging data were smoothed with Gaussian kernel of 8 mm full-width at half maximum. All analyses were restricted to a gray matter mask derived from the gray matter template provided by SPM (probability of gray matter >0.2).

Group Comparisons of DAT-SPECT and sMRI Data

The pre-processed DAT-SPECT and GMV images PD patients and healthy controls were entered into voxel-wise general linear models (GLMs) with diagnosis as a factor and controlled for age and sex (Table 1). In GMV analyses, total intracranial volume was additionally controlled for. The contrast showing reduced signal in PD was evaluated for both. A voxel-wise family-wise error (FWE) corrected threshold of p < 0.05 was applied for these analyses.

Functional Connectivity of the Striatum

As motion has been reported to be a critical factor affecting resting state signal (Van Dijk et al., 2012), we regressed out the effects of motion for each subject using a 24 parameter model (6 translational and rotational motion parameters, 6 translational and rotational acceleration parameters, and the squared terms of both) that allows for a rigid control of such data (Yan et al., 2013). As PPMI ensured that only data passing pre-specified quality checks are uploaded to the database no patient had to be excluded due to extensive motion (above 1 voxel-size). Functional connectivity maps were computed using the REST toolbox (Song et al., 2011) implemented in Matlab (MathWorks). A sphere with a radius of 5 mm around the striatal peak voxel providing for DAT-SPECT the maximum difference to healthy controls was used as a seed for computation of Fisher’s z-transformed Pearson correlational maps of the average signal within this region and the rest of the brain (SFC) (Figure 1A). We used the default settings in the REST toolbox including removing of a linear trend and application of a low- and high-pass filter restricting the computation to frequencies between 0.01 and 0.08 Hz to reduce contribution of physiological noise.

FIGURE 1

FIGURE 1. Results of dopamine transporter single photon emission tomography (DAT-SPECT) and functional connectivity analyses. (A) DAT-SPECT results showing reduced striatal activity in Parkinson’s disease patients relative to healthy controls and the seed region used for subsequent functional connectivity analyses. (B) Regions showing significant negative correlations with the Unified Parkinson’s Disease Rating Scale (UPDRS) total score in the functional connectivity analysis. (C) Regions showing significant negative correlations with different parts of the UPDRS in the functional connectivity analysis with all subscales. (D) Regions showing significant negative correlations with different parts of the UPDRS in the functional connectivity analysis with separate subscores. All images are displayed in neurological convention.

Statistical Analysis of rsMRI Data

To facilitate the interpretation of potential correlations between UPDRS subscales and SFC maps, we first aimed to quantify the shared variance across these subscales by computing pair-wise Pearson correlation coefficients. We then aimed to understand how the SFC maps are linked to the UPDRS subscales with and without accounting for the variance shared by the subscales. To address this question, the computed SFC maps were first entered into a voxel-wise GLM including UPDRS parts I, II and III as regressors in the same model therewith accounting for each other’s contribution. To understand the link between UPDRS subscales and SFC maps without controlling for the shared variance, the analysis was then recomputed including each of the subscales and the total score in separate GLMs. All GLMs were tested for negative correlations between SFC and the respective UPDRS regressors. All analyses were further controlled for age, sex, dominant side of symptoms and 3 binary medication status covariates (levodopa, dopamine agonists or other PD medications: all yes/no, no information on dose available in the data base). As we were more interested in the identification of anatomical networks underlying the UPDRS subscales rather than the peak regions associated with each of them, an uncorrected voxel-wise threshold of p < 0.01 was applied combined with a whole-brain FWE corrected cluster threshold of p < 0.05 adjusted for non-stationarity of smoothness. Recently, Eklund et al. (2016), reported that this type of threshold may lead to an increased false positive rate for fMRI type of analyses. To address this issue we ran 1000 permutations of our parametric design randomly assigning the PD clinical data to individual scans and recomputing the contrasts testing for significant correlations with randomized UPDRS scores at the above voxel- and cluster-wise threshold. This evaluation resulted in 4.9% [99% confidence interval: 4.2–5.8%] of permutations showing significant clusters at the chosen statistical threshold strongly suggesting the validity of the applied multiple comparison correction for our data. We also report the exact p-values for the obtained clusters as derived from this null distribution (Table 2).

TABLE 2

TABLE 2. Regions showing a negative correlation with UPDRS subscales in the multiple regression analysis including all subscales.

To quantify the contribution of SFC within the identified networks to the observed symptom severity, we extracted its eigenvariate from all significant clusters adjusted for variance explained by covariates of no interest. To visualize the identified relationships between clinical scales and SFC measures, the eigenvariates from each cluster were plotted against each UPDRS subscale (Figure 2).

FIGURE 2

FIGURE 2. Results of multiple regression analyses based on functional connectivity regions identified in correlations with the UPDRS subscales. (A) Observed UPDRS subscales and total scores vs. those predicted in a leave-one-out cross-validation using regions identified as being correlated with UPDRS part II. (B) Observed UPDRS subscales and total scores vs. those predicted in a leave-one-out cross-validation using regions identified as being correlated with UPDRS part III. (C) Observed UPDRS subscales and total scores vs. those predicted in a leave-one-out cross-validation using regions identified as being correlated with UPDRS total scores.

Combined Analysis of rsMRI and DAT-SPECT Data

To test if DAT-SPECT and rsMRI indeed provide a differential contribution to clinical symptoms, we performed a mediator analysis. A mediator analysis tests if an initially significant correlation (p < 0.05) between two variables is no longer significant after controlling for a third variable that is correlated with both. The third variable is then considered as a mediator between both. In contrast, if the correlation remains significant after controlling for the third variable the correlation between both is considered as independent of the respective variable. For this, partial correlations were computed to test if SFC networks identified in regression analyses with UPDRS subscales and DAT striatal binding explain independent variance in PD symptom severity going beyond the effect of striatal dopaminergic signal loss. The following partial correlations were computed for this purpose (all controlling for age, sex and PD medication status: levodopa, dopamine agonists, other):

(1) Between DAT striatal binding (ipsi- and contralateral putamen and caudate nucleus), UPDRS subscales (UPDRS I, II and III) and SFC cluster eigenvariates identified in voxel-wise regression analyses (SFC clusters showing significant correlations with UPDRS II, III or the total score),

(2) Between UPDRS subscales and SFC eigenvariates controlling for DAT striatal binding,

(3) Between UPDRS subscales and DAT striatal binding controlling for SFC.

Additionally, we formally compared the partial correlation coefficients obtained in (1) between DAT striatal binding and UPDRS subscales and SFC and UPDRS subscales using Steiger’s z-test for correlation coefficients (p < .05 two-sided) computed on overlapping variables (Wuensch, 2007).

Further, we aimed to better understand the directionality of relationships identified in partial correlation analyses between striatal dopaminergic uptake, SFC and UPDRS. For this, we performed a Bayesian Network analyses using the three phase dependency algorithm implemented in the Matlab-based Causal Explorer toolbox (Cheng et al., 1998; Statnikov et al., 2010). In brief, this algorithm creates a Bayesian network testing for significant causal relationships between any of the variables (p < 0.05) in presence of all other included variables. If for example a variable A can be regarded in a causal relationship as a parent or a child of variable B a directed graph is assigned to this relationship. In case no directionality can be derived despite a significant relationship a bidirectional graph is assigned to such connection. The magnitude of the identified associations was quantified using determination coefficients.

Results

Demographic and clinical description of the rsMRI cohort including the status of PD medication and clinical severity are provided in Table 1.

Results of Group Comparisons for DAT-SPECT and sMRI Data

No atrophic regions were identified in GMV comparisons of PD patients and healthy control subjects. DAT binding in PD was significantly reduced in all striatal regions beside posterior caudate with most pronounced reductions in the right posterior putamen (Figure 1A). This region was used as seed for subsequent computations of SFC maps.

Results of rsMRI Data Analyses

As expected, significant correlations were observed between the UPDRS subscales (UPDRS I and II: r = 0.56; p < 0.001, UPDRS I and III: r = 0.33; p < 0.001, UPDRS II and III: r = 0.46; p < 0.001) indicating that the subscales though sharing up to 31% of variance also provide sufficiently distinct information on disease severity. In the multiple regression analysis of imaging measures with the UPRDS subscales accounting for their shared variance, we found significant negative correlations between SFC and UPDRS part II and III but not with UPDRS part I. Extensive correlations with the UPDRS part II were observed in the vermis and bilateral cerebellar, precuneal, occipital and right temporal and posterior insular regions (this cluster is further referred to as SFC II) (Table 2; Figure 1C). We also identified significant negative correlations with the UPDRS part III in bilateral premotor, supplementary motor, anterior and middle cingulate, and dorsolateral prefrontal regions (SFC III) (Table 2; Figure 1C). In correlational analysis with single UPDRS subscales not controlling for the shared variance we found similar but more wide-spread anatomical network for both UPDRS II and III (Figure 1D). In regression analysis with UPDRS total scores, significant negative correlations with SFC were detected in regions covering an extensive cortico-subcortical-cerebellar network and including all regions identified for both UPDRS part II and III (SFC total) (Table 2; Figure 1B). In subsequent correlational analyses, eigenvariates from SFC II, SFC III and SFC total explained up to 19% of variance in these subscales (Figure 2).

Results of Combined rsMRI and DAT-SPECT Data Analyses

In partial correlation analyses only controlling for covariates of no interest but not for shared variance between DAT striatal binding and SFC, the UPRDS I subscale showed a marginally significant correlation with SFC total regions (r = -0.22, p = 0.071) (Figure 3A). As expected, the UPDRS II subscale was significantly correlated with SFC II regions (r = -0.46, p < 0.001) but also with SFC total (r = -0.47, p < 0.001), SFC III (r = -0.27, p = 0.021), with DAT striatal binding in ipsilateral putamen (r = -0.31, p = 0.009) and in bilateral caudate nucleus (ipsilateral: r = -0.36, p = 0.002; contralateral: r = -0.29, p = 0.015). The UPDRS III subscale was significantly correlated with SFC III regions (r = -0.42, p < 0.001) but also with SFC total (r = -0.42, p < 0.001), SFC II (r = -0.26, p = 0.029) and with ipsilateral caudate DAT binding (r = -0.25, p = 0.034). No significant correlations were observed between SFC and DAT striatal uptake. There was no significant difference between DAT-SPECT and SFC with respect to correlation strength with UPDRS I (Z = 0.15, p = 0.878), II (Z = 0.82, p = 0.414) and III (Z = 1.17, p = 0.244).

FIGURE 3

FIGURE 3. Results of cross-modal partial correlation and Bayesian network analyses. (A) Results of partial correlation analyses between imaging and clinical measures when controlling for demographic factors are displayed. Significant partial correlations between the corresponding imaging measures and UPDRS subscales are indicated by solid lines. (B) Results of partial correlation analyses between striatal functional connectivity (SFC) and clinical measures when controlling for demographic factors and DAT striatal binding are displayed. Significant partial correlations between the corresponding imaging measures and UPDRS subscales are indicated by solid lines. (C) Results of partial correlation analyses between DAT striatal binding and clinical measures when controlling for demographic factors and SFC are displayed. Significant partial correlations between the corresponding imaging measures and UPDRS subscales are indicated by solid lines. (D) Results of Bayesian network analyses are displayed. DAT, dopamine transporter; caud, caudate nucleus; put, putamen; SFC I, II and III, striatal functional connectivity regions associated with UPDRS I, II and III, respectively; R², determination coefficient.

In partial correlations between UPDRS subscales and SFC controlling for DAT striatal uptake all significant correlations reported above, except the non-significant correlation between SFC total and UPDRS I (r = -0.18, p = 0.146), remained significant or marginally significant (UPDRS II and SFC II: r = -0.45, p < 0.001; UPDRS II and SFC III: r = -0.22, p = 0.071; UPDRS II and SFC total: r = -0.44, p < 0.001; UPDRS III and SFC III: r = -0.40, p = 0.001; UPDRS III and SFC total: r = -0.39, p = 0.001; UPDRS III and SFC II: r = -0.24, p = 0.050) (Figure 3B).

Similarly, in partial correlations between UPDRS subscales and DAT striatal binding controlling for SFC, all significant correlations between UPDRS II and DAT striatal binding remained significant (UPDRS II and ipsilateral putamen DAT: r = -0.27, p = 0.028; UPDRS II and ipsilateral caudate DAT: r = -0.32, p = 0.007; UPDRS II and contralateral caudate DAT: r = -0.24, p = 0.045) (Figure 3C). The initially significant correlation between UPDRS III and ipsilateral caudate DAT binding was no longer significant after controlling for SFC indicating a mediator effect of SFC onto the initially observed association between both (r = -0.17, p = 0.164).

The Bayesian network analyses revealed an independent contribution of DAT striatal binding and II and III to UDPRS symptoms (Figure 3D). More specifically, as expected SFC II was linked to UPDRS II, and SFC III to UPDRS III. DAT striatal binding in ipsilateral caudate nucleus was linked to UPDRS II. All other DAT striatal binding regions were only indirectly linked to UPDRS through this relationship. Interestingly, UPDRS I and III were directionally contributing to UPDRS II scale. No significant link was identified between any DAT striatal binding and SFC regions nor between UPDRS I and any of the imaging measures.

Discussion

Here we analyzed the association between PD symptom domains, striatal dopaminergic deficits and SFC at rest in the largest currently available cohort of de novo PD patients. For the first time, we show that SFC provides an independent and comparably strong link to PD clinical symptoms as the striatal dopaminergic deficits. We also find distinct anatomical SFC networks to contribute to different PD symptom domains as captured by UPDRS part II and III; the first reflecting the ability to perform daily routine automatized motor functions and the second measuring directly clinician-rated motor symptoms associated with PD.

The mediator and Bayesian network analyses confirm the independent contribution of the SFC and dopaminergic deficits to PD clinical symptoms as measured by UPDRS II. More specifically, the Bayesian network analyses suggest that dopaminergic signal and SFC independently contribute to the UPDRS II activities of daily living subscale. In contrast, both the mediator and the Bayesian network analysis suggest a mediator effect of SFC onto the observed relationship between UPDRS III and DAT striatal binding reported in previous studies and reviews (Benamer et al., 2000; Eshuis et al., 2006; Nobili et al., 2010; McGhee et al., 2013). In line with this, the lack of a significant correlation between DAT striatal binding and SFC supports the idea of independent pathological mechanisms underlying both types of measures. Conceptually, DAT binding reflects the loss of dopaminergic projections from substantia nigra into the striatum. However, it does not provide direct information on the resulting downstream perturbations of striato-cortical loops which are presumably more directly reflected by SFC. In example, DAT striatal binding is already substantially reduced in recently diagnosed PD with mild symptom severity (Tissingh et al., 1998). This finding supports the idea of striatal dopaminergic projections not being directly associated with the manifestation of clinical symptoms, i.e., due to compensational mechanisms (Lloyd, 1977).

We find the UPDRS part II being closely related to SFC to vermis, cerebellum and posterior cortical regions, including the temporo-parietal junction. These relationships match well-known brain networks for retention of highly automatized motor functions and bodily self-awareness (Doyon et al., 1998; Lang and Bastian, 2002; Ionta et al., 2011; Serino et al., 2013). In contrast, UPDRS part III is more closely linked to SFC to prefrontal regions associated with motor and executive functions. The lack of significant correlations with UPDRS part I may be related to the higher heterogeneity of functions covered by this subscale and/or involvement of other disease mechanisms (Murai et al., 2001; Goetz et al., 2008). Both identified anatomical networks are strikingly similar to regions reported in previous studies in PD and healthy controls as being functionally connected to the posterior putamen (Helmich et al., 2010; Tziortzi et al., 2014). The prefrontal regions identified in our study are also consistent with regions showing diffusion alterations in PD patients (Kendi et al., 2008). The anatomical networks are also consistent with recent findings reporting these regions to be associated with altered activity and connectivity in PD patients during execution and attention to novel versus automatic movements (Wu et al., 2014). The cerebellar and vermis SFC to the striatum as well as the within-striatum connectivity has been shown to be in general decreased in PD and increased after levodopa administration, supporting the hypothesized involvement of this striato-cerebellar network in generation of some of the PD symptoms (Hacker et al., 2012; Jech et al., 2013; Bell et al., 2015). Also importantly, several studies suggested a potential compensatory role of both cerebellar and premotor synchronization in PD (Wu et al., 2009; Hacker et al., 2012). Our findings are fully in line with these hypothesized mechanisms as one would expect that a reduction or deficit in the compensatory synchronization would lead to increased PD symptoms.

As the same striatal seed region was used for all analyses the identified differences in underlying networks can be therefore regarded as symptom specific desynchronization of posterior putamen activity with corresponding SFC networks. This notion of differential processes involved in generation of PD symptoms is also supported by the relatively low to moderate correlations observed between UPDRS subscales. By controlling our statistical design for different medication and for dominant side of symptoms and testing for potential effects of atrophy we further account for known potential factors affecting functional connectivity in this patient group (Hacker et al., 2012; Jech et al., 2013).

Overall, these findings suggest the potential usability of resting state SFC as a surrogate biomarker endpoint for PD symptoms. However, several further studies will be needed to first validate this identified relationship and to establish longitudinal and potentially mechanistic links between SFC networks and the respective symptoms. Nonetheless, the distinct anatomical networks identified in our study indicate the existence of differential neuronal mechanisms underlying the corresponding PD symptoms. These findings are also in line with results of previous factorial analyses, suggesting contribution of different factors to UPDRS subscales (Stebbins and Goetz, 1998; Goetz et al., 2008). Notably, we find that both anatomical networks associated with UPDRS part II and III also correlate with the UPDRS total score. This finding supports the use of the total score as a composite measure reflecting several PD related pathological processes. However, our results also suggest that caution might be required when using the UPDRS total score for efficacy evaluation of focal treatment procedures such as transcranial magnetic stimulation (Shimamoto et al., 1999; Pal et al., 2010). Based on our results, one would for example expect that a stimulation focus on the premotor region would predominantly result in an improvement of the UPDRS part III but not of the other subscales. This prediction is consistent with the results of a recent large double-blind sham-controlled transcranial magnetic stimulation study reporting an improvement in UPDRS part III but not in UPDRS part I and II after stimulation of the supplementary motor area (Shirota et al., 2013).

We do not identify any statistically significant atrophy in our study in the comparison of de novo PD patients and control subjects suggesting a neglectable contribution of atrophy to the observed correlations between functional connectivity and UPDRS subscales (Dukart and Bertolino, 2014). Considering the large sample size included for structural analyses, the lack of findings in our study is unlikely to be a result of insufficient power to detect such effects. The lack of significant atrophy is in contrast to a previous study reporting striatal atrophy in subsample of the cohort included here using a different volumetric estimation algorithm combined with an independent component based approach leading to a more lenient statistical threshold (Zeighami et al., 2015). In a recent meta-analysis of voxel-based morphometry studies in PD the orbitofrontal gray matter volume was identified as consistently reduced in comparisons with healthy controls (Pan et al., 2012). However, all studies included in this meta-analysis targeted more advanced PD populations. Considering the neurodegenerative nature of PD, occurrence of more pronounced atrophy in later stages is therefore likely and might explain the lack of findings in our study. Nonetheless, in line with several previous studies our results strongly question the usability of gray matter atrophy as detected by voxel-based morphometry as an early biomarker of PD (Tessa et al., 2008; Dalaker et al., 2009; Menke et al., 2014).

Author Contributions

JD designed the study, analyzed the data and wrote the manuscript. FS and AB designed the study and reviewed the manuscript.

Funding

JD, FS, and AB are current or former full-time employees of F. Hoffmann-La Roche, Basel Switzerland. The authors received no specific funding for this work. F. Hoffmann-La Roche provided financial contribution in the form of salary for all authors but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

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.

Acknowledgments

Data used in the preparation of this article were obtained from the Parkinson’s Progression Markers Initiative (PPMI) database (www.ppmi-info.org/data). For up-to-date information on the study, visit www.ppmi-info.org. PPMI – a public-private partnership – is funded by the Michael J. Fox Foundation for Parkinson’s Research and funding partners, including Abbvie, Avid Radiopharmaceuticals, Biogen Idec, Bristol-Myers Squibb, Covance, GE Healthcare, Genentech, GlaxoSmithKline, Lilly, Lundbeck, Merck, Meso Scale Discovery, Pfizer, Piramal, Roche and UCB.

Footnotes

1. ^ http://www.ppmi-info.org/
2. ^ http://www.fil.ion.ucl.ac.uk/spm

References

Bell, P. T., Gilat, M., O’Callaghan, C., Copland, D. A., Frank, M. J., Lewis, S. J. G., et al. (2015). Dopaminergic basis for impairments in functional connectivity across subdivisions of the striatum in Parkinson’s disease. Hum. Brain Mapp. 36, 1278–1291. doi: 10.1002/hbm.22701

PubMed Abstract | CrossRef Full Text | Google Scholar

Benamer, H. T., Patterson, J., Wyper, D. J., Hadley, D. M., Macphee, G. J. A., and Grosset, D. G. (2000). Correlation of Parkinson’s disease severity and duration with 123I-FP-CIT SPECT striatal uptake. Mov. Disord. 15, 692–698. doi: 10.1002/1531-8257(200007)15:4<692::AID-MDS1014>3.0.CO;2-V

CrossRef Full Text | Google Scholar

Benoit, M., and Robert, P. H. (2011). Imaging correlates of apathy and depression in Parkinson’s disease. J. Neurol. Sci. 310, 58–60. doi: 10.1016/j.jns.2011.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, J., Bell, D., Greiner, R., and Liu, W. (1998). Learning Bayesian networks from data: an information-theory based approach. Artif. Intell. 137, 43–90. doi: 10.1016/S0004-3702(02)00191-1

CrossRef Full Text | Google Scholar

Dalaker, T. O., Larsen, J. P., Bergsland, N., Beyer, M. K., Alves, G., Dwyer, M. G., et al. (2009). Brain atrophy and white matter hyperintensities in early Parkinson’s disease. Mov. Disord. 24, 2233–2241. doi: 10.1002/mds.22754

PubMed Abstract | CrossRef Full Text | Google Scholar

Dauer, W., and Przedborski, S. (2003). Parkinson’s disease: mechanisms and models. Neuron 39, 889–909. doi: 10.1016/S0896-6273(03)00568-3

CrossRef Full Text | Google Scholar

Doyon, J., Laforce, R., Bouchard, G., Gaudreau, D., Roy, J., Poirier, M., et al. (1998). Role of the striatum, cerebellum and frontal lobes in the automatization of a repeated visuomotor sequence of movements. Neuropsychologia 36, 625–641. doi: 10.1016/S0028-3932(97)00168-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dukart, J., and Bertolino, A. (2014). When structure affects function–the need for partial volume effect correction in functional and resting state magnetic resonance imaging studies. PLoS ONE 9:e114227. doi: 10.1371/journal.pone.0114227

PubMed Abstract | CrossRef Full Text | Google Scholar

Eklund, A., Nichols, T. E., and Knutsson, H. (2016). Cluster failure: why fMRI inferences for spatial extent have inflated false-positive rates. Proc. Natl. Acad. Sci. U.S.A. 113, 7900–7905. doi: 10.1073/pnas.1602413113

PubMed Abstract | CrossRef Full Text | Google Scholar

Eshuis, S. A., Maguire, R. P., Leenders, K. L., Jonkman, S., and Jager, P. L. (2006). Comparison of FP-CIT SPECT with F-DOPA PET in patients with de novo and advanced Parkinson’s disease. Eur. J. Nucl. Med. Mol. Imaging 33, 200–209. doi: 10.1007/s00259-005-1904-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Fahn, S., and Elton, R. L. (1987). “Members of the updrs development committee,” in Recent Developments in Parkinson’s Disease, Vol. 2, eds S. Fahn, C. D. Marsden, D. B. Calne, and M. Goldstein (Florham Park, NJ: Macmillan Health Care Information), 153–163.

Google Scholar

Goetz, C. G., Tilley, B. C., Shaftman, S. R., Stebbins, G. T., Fahn, S., Martinez-Martin, P., et al. (2008). Movement Disorder Society-sponsored revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): scale presentation and clinimetric testing results. Mov. Disord. 23, 2129–2170. doi: 10.1002/mds.22340

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, G. (1985). Supplementary motor area structure and function: review and hypotheses. Behav. Brain Sci. 8, 567–588. doi: 10.1017/S0140525X00045167

CrossRef Full Text | Google Scholar

Hacker, C., Perlmutter, J., Criswell, S., Ances, B., and Snyder, A. (2012). Resting state functional connectivity of the striatum in Parkinson’s disease. Brain 135(Pt 12), 3699–3711. doi: 10.1093/brain/aws281

PubMed Abstract | CrossRef Full Text | Google Scholar

Ham, J. H., Cha, J., Lee, J. J., Baek, G.-M., Sunwoo, M. K., Hong, J. Y., et al. (2015). Nigrostriatal dopamine-independent resting-state functional networks in Parkinson’s disease. Neuroimage 119, 296–304. doi: 10.1016/j.neuroimage.2015.06.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Helmich, R. C., Derikx, L. C., Bakker, M., Scheeringa, R., Bloem, B. R., and Toni, I. (2010). Spatial remapping of cortico-striatal connectivity in Parkinson’s disease. Cereb. Cortex 20, 1175–1186. doi: 10.1093/cercor/bhp178

PubMed Abstract | CrossRef Full Text | Google Scholar

Ionta, S., Heydrich, L., Lenggenhager, B., Mouthon, M., Fornari, E., Chapuis, D., et al. (2011). Multisensory mechanisms in temporo-parietal cortex support self-location and first-person perspective. Neuron 70, 363–374. doi: 10.1016/j.neuron.2011.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Jech, R., Mueller, K., Schroeter, M. L., and Růžička, E. (2013). Levodopa increases functional connectivity in the cerebellum and brainstem in Parkinson’s disease. Brain 136, e234. doi: 10.1093/brain/awt015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kendi, A. K., Lehericy, S., Luciana, M., Ugurbil, K., and Tuite, P. (2008). Altered diffusion in the frontal lobe in Parkinson disease. Am. J. Neuroradiol. 29, 501–505. doi: 10.3174/ajnr.A0850

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, C. E., and Bastian, A. J. (2002). Cerebellar damage impairs automaticity of a recently practiced movement. J. Neurophysiol. 87, 1336–1347.

PubMed Abstract | Google Scholar

Lebedev, A. V., Westman, E., Simmons, A., Lebedeva, A., Siepel, F. J., Pereira, J. B., et al. (2014). Large-scale resting state network correlates of cognitive impairment in Parkinson’s disease and related dopaminergic deficits. Front. Syst. Neurosci. 8:45. doi: 10.3389/fnsys.2014.00045

PubMed Abstract | CrossRef Full Text | Google Scholar

Lloyd, K. G. (1977). “CNS compensation to dopamine neuron loss in Parkinson’s Disease,” in Parkinson’s Disease: Neurophysiological, Clinical, and Related Aspects, eds F. S. Messiha and A. D. Kenny (New York, NY: Springer), 255–266.

Google Scholar

Lotharius, J., and Brundin, P. (2002). Pathogenesis of Parkinson’s disease: dopamine, vesicles and α-synuclein. Nat. Rev. Neurosci. 3, 932–942. doi: 10.1038/nrn983

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Martin, P., Gil-Nagel, A., Gracia, L. M., Gómez, J. B., Martinez-Sarriés, J., and Bermejo, F. (1994). Unified Parkinson’s disease rating scale characteristics and structure. Mov. Disord. 9, 76–83. doi: 10.1002/mds.870090112

PubMed Abstract | CrossRef Full Text | Google Scholar

McGhee, D. J., Royle, P. L., Thompson, P. A., Wright, D. E., Zajicek, J. P., and Counsell, C. E. (2013). A systematic review of biomarkers for disease progression in Parkinson’s disease. BMC Neurol. 13:35. doi: 10.1186/1471-2377-13-35

PubMed Abstract | CrossRef Full Text | Google Scholar

Menke, R. A., Szewczyk-Krolikowski, K., Jbabdi, S., Jenkinson, M., Talbot, K., Mackay, C. E., et al. (2014). Comprehensive morphometry of subcortical grey matter structures in early-stage Parkinson’s disease. Hum. Brain Mapp. 35, 1681–1690. doi: 10.1002/hbm.22282

PubMed Abstract | CrossRef Full Text | Google Scholar

Murai, T., Müller, U., Werheid, K., Sorger, D., Reuter, M., Becker, T., et al. (2001). In vivo evidence for differential association of striatal dopamine and midbrain serotonin systems with neuropsychiatric symptoms in Parkinson’s disease. J. Neuropsychiatry Clin. Neurosci. 13, 222–228. doi: 10.1176/jnp.13.2.222

PubMed Abstract | CrossRef Full Text | Google Scholar

Nobili, F., Campus, C., Arnaldi, D., De Carli, F., Cabassi, G., Brugnolo, A., et al. (2010). Cognitive-nigrostriatal relationships in de novo, drug-naïve Parkinson’s disease patients: a [I-123] FP-CIT SPECT study. Mov. Disord. 25, 35–43. doi: 10.1002/mds.22899

PubMed Abstract | CrossRef Full Text | Google Scholar

Pal, E., Nagy, F., Aschermann, Z., Balazs, E., and Kovacs, N. (2010). The impact of left prefrontal repetitive transcranial magnetic stimulation on depression in Parkinson’s disease: a randomized, double-blind, placebo-controlled study. Mov. Disord. 25, 2311–2317. doi: 10.1002/mds.23270

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, P. L., Song, W., and Shang, H. F. (2012). Voxel-wise meta-analysis of gray matter abnormalities in idiopathic Parkinson’s disease. Eur. J. Neurol. 19, 199–206. doi: 10.1111/j.1468-1331.2011.03474.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrides, M. (2005). Lateral prefrontal cortex: architectonic and functional organization. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360, 781–795. doi: 10.1098/rstb.2005.1631

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirker, W. (2003). Correlation of dopamine transporter imaging with parkinsonian motor handicap: how close is it? Mov. Disord. 18, S43–S51. doi: 10.1002/mds.10579

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramaker, C., Marinus, J., Stiggelbout, A. M., and van Hilten, B. J. (2002). Systematic evaluation of rating scales for impairment and disability in Parkinson’s disease. Mov. Disord. 17, 867–876. doi: 10.1002/mds.10248

PubMed Abstract | CrossRef Full Text | Google Scholar

Rieckmann, A., Gomperts, S. N., Johnson, K. A., Growdon, J. H., and Van Dijk, K. R. A. (2015). Putamen–midbrain functional connectivity is related to striatal dopamine transporter availability in patients with Lewy body diseases. Neuroimage Clin. 8, 554–559. doi: 10.1016/j.nicl.2015.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Serino, A., Alsmith, A., Costantini, M., Mandrigin, A., Tajadura-Jimenez, A., and Lopez, C. (2013). Bodily ownership and self-location: components of bodily self-consciousness. Conscious. Cogn. 22, 1239–1252. doi: 10.1016/j.concog.2013.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimamoto, H., Morimitsu, H., Sugita, S., Nakahara, K., and Shigemori, M. (1999). Therapeutic effect of repetitive transcranial magnetic stimulation in Parkinson’s disease. Rinsho Shinkeigaku 39, 1264–1267. doi: 10.1007/s00405-013-2851-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirota, Y., Ohtsu, H., Hamada, M., Enomoto, H., and Ugawa, Y. (2013). Supplementary motor area stimulation for Parkinson disease A randomized controlled study. Neurology 80, 1400–1405. doi: 10.1212/WNL.0b013e31828c2f66

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X.-W., Dong, Z.-Y., Long, X.-Y., Li, S.-F., Zuo, X.-N., Zhu, C.-Z., et al. (2011). REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS ONE 6:e25031. doi: 10.1371/journal.pone.0025031

PubMed Abstract | CrossRef Full Text | Google Scholar

Statnikov, A., Tsamardinos, I., Brown, L. E., and Aliferis, C. F. (2010). Causal explorer: a matlab library of algorithms for causal discovery and variable selection for classification. Causation Predict. Chall. Chall. Mach. Learn. 2, 267.

Google Scholar

Stebbins, G. T., and Goetz, C. G. (1998). Factor structure of the Unified Parkinson’s Disease Rating Scale: motor examination section. Mov. Disord. 13, 633–636. doi: 10.1002/mds.870130404

PubMed Abstract | CrossRef Full Text | Google Scholar

Tessa, C., Giannelli, M., Della Nave, R., Lucetti, C., Berti, C., Ginestroni, A., et al. (2008). A whole-brain analysis in de novo Parkinson disease. Am. J. Neuroradiol. 29, 674–680. doi: 10.3174/ajnr.A0900

PubMed Abstract | CrossRef Full Text | Google Scholar

Tissingh, G., Booij, J., Bergmans, P., Winogrodzka, A., Janssen, A. G., van Royen, E. A., et al. (1998). Iodine-123-N-(omega)-fluoropropyl-2 (beta)-carbomethoxy-3 (beta)-(4-iodophenyl) tropane SPECT in healthy controls and early-stage, drug-naive Parkinson’s disease. J. Nucl. Med. 39, 1143–1148.

Google Scholar

Tziortzi, A. C., Haber, S. N., Searle, G. E., Tsoumpas, C., Long, C. J., Shotbolt, P., et al. (2014). Connectivity-based functional analysis of dopamine release in the striatum using diffusion-weighted MRI and positron emission tomography. Cereb. Cortex 24, 1165–1177. doi: 10.1093/cercor/bhs397

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

Wolters, E. C. (2009). Non-motor extranigral signs and symptoms in Parkinson’s disease. Parkinsonism Relat. Disord. 15, S6–S12. doi: 10.1016/S1353-8020(09)70770-9

CrossRef Full Text | Google Scholar

Wu, T., Liu, J., Zhang, H., Hallett, M., Zheng, Z., and Chan, P. (2014). Attention to automatic movements in Parkinson’s Disease: modified automatic mode in the Striatum. Cereb. Cortex 25, 3330–3342. doi: 10.1093/cercor/bhu135

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Long, X., Wang, L., Hallett, M., Zang, Y., Li, K., et al. (2011). Functional connectivity of cortical motor areas in the resting state in Parkinson’s disease. Hum. Brain Mapp. 32, 1443–1457. doi: 10.1002/hbm.21118

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, T., Long, X., Zang, Y., Wang, L., Hallett, M., Li, K., et al. (2009). Regional homogeneity changes in patients with Parkinson’s disease. Hum. Brain Mapp. 30, 1502–1510. doi: 10.1002/hbm.20622

PubMed Abstract | CrossRef Full Text | Google Scholar

Wuensch, K. L. (2007). Comparing Correlation Coefficients, Slopes, and Intercepts. Stat. Lessons. Available at: http://core.ecu.edu/psyc/wuenschk/docs30/CompareCorrCoeff.pdf

Yan, C.-G., Craddock, R. C., He, Y., and Milham, M. P. (2013). Addressing head motion dependencies for small-world topologies in functional connectomics. Front. Hum. Neurosci. 7:910. doi: 10.3389/fnhum.2013.00910

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeighami, Y., Ulla, M., Iturria-Medina, Y., Dadar, M., Zhang, Y., Larcher, K. M.-H., et al. (2015). Network structure of brain atrophy in de novo Parkinson’s disease. Elife 4:e08440. doi: 10.7554/eLife.08440

PubMed Abstract | CrossRef Full Text | Google Scholar