Network connectivity determines cortical thinning in early parkinson’s disease progression

ABSTRACT Here we test the hypothesis that the neurodegenerative process in Parkinson’s disease (PD) moves stereotypically along neural networks, possibly reflecting the spread of toxic

alpha-synuclein molecules. PD patients (_n_ = 105) and matched controls (_n_ = 57) underwent T1-MRI at entry and 1 year later as part of the Parkinson’s Progression Markers Initiative. Over

this period, PD patients demonstrate significantly greater cortical thinning than controls in parts of the left occipital and bilateral frontal lobes and right somatomotor-sensory cortex.

Cortical thinning is correlated to connectivity (measured functionally or structurally) to a “disease reservoir” evaluated by MRI at baseline. The atrophy pattern in the ventral frontal

lobes resembles one described in certain cases of Alzheimer’s disease. Our findings suggest that disease propagation to the cortex in PD follows neuronal connectivity and that disease spread

to the cortex may herald the onset of cognitive impairment. SIMILAR CONTENT BEING VIEWED BY OTHERS MICROSTRUCTURAL BUT NOT MACROSTRUCTURAL CORTICAL DEGENERATION OCCURS IN PARKINSON’S

DISEASE WITH MILD COGNITIVE IMPAIRMENT Article Open access 09 November 2022 AKINETIC RIGID SYMPTOMS ARE ASSOCIATED WITH DECLINE IN A CORTICAL MOTOR NETWORK IN PARKINSON’S DISEASE Article

Open access 24 August 2020 ELEVATED CAUDATE CONNECTIVITY IN COGNITIVELY NORMAL PARKINSON’S DISEASE PATIENTS Article Open access 21 October 2020 INTRODUCTION The network-spread hypothesis

proposes that neurodegenerative diseases target intrinsic brain networks. Specifically, it is suggested that the progression of disease follows network connectivity and that neurotoxicity is

caused by the spread and accumulation of toxic agents along the brain’s neuronal connectome1,2. In the case of Parkinson’s disease (PD), the likely agent is misfolded α-synuclein (α-syn).

In its normal conformation, α-syn is an endogenous protein that is soluble and contributes to synaptic vesicular function3; however, in synucleinopathies such as PD, α-syn is present in a

hyperphosphorylated and insoluble state and found to adopt an abnormal configuration into amyloid fibril aggregates4. A model first proposed by Braak suggested that patterns of abnormal

α-syn accumulation observed at postmortem were consistent with a stereotyped propagation from brainstem to subcortical areas to cortex, which paralleled the clinical progression of the

disease5. Numerous lines of evidence support this model: misfolded α-syn aggregates are identified at postmortem in the brains of PD patients, and the protein has been shown to spread

through neuroanatomical pathways causing cell death in animal models due to its neurotoxicity4; however, the model remains controversial and direct human evidence is sparse6,7. There is

currently no imaging technique that allows in vivo measurement of α-syn in the brain. However, we previously observed a pattern of brain atrophy in cross-sectional magnetic resonance imaging

(MRI) scans of de novo PD patients from the Parkinson’s Progression Markers Initiative (PPMI) database8 that was consistent with the network-spread hypothesis9. The regional distribution of

atrophy identified using deformation-based morphometry (DBM) overlapped with an intrinsic connectivity network, and the pattern suggested a disease epicenter in the substantia nigra. The

atrophy distribution was mostly subcortical, involving brainstem, striatum, thalamus, and medial temporal lobe. Recent studies using voxel-based morphometry (VBM) to assess gray matter

density have begun to highlight a pattern of cortical changes in PD10,11, but there have been discrepant findings12,13. Inconsistencies are not only due to heterogeneity in disease

progression and small sample sizes, but also to methodological issues14. Cortical thickness analysis may be better able to delineate cortical morphometrical patterns and may be more

sensitive to age or disease-related regional gray matter changes than VBM15,16. While cortical thinning is normal in the aging process17, deviations from typical patterns can be used as

diagnostic indicators for neurodegenerative disorders (e.g., Alzheimer’s disease (AD)) even at early stages18. Several cross-sectional studies of MRI-derived cortical thickness in

early-stage PD have reported cortical thinning in parietal and premotor regions compared to controls16,19,20,21. Moreover, cortical thinning appears to correlate with poorer cognitive

performance22. However, cortical thickness in PD has yet to be studied in a controlled longitudinal setting, and little is known as to how it may correlate with other disease-related

measures. Disease spread to cortical areas may be accompanied by the start of cognitive impairment. Dementia is now recognized as a common, perhaps inevitable, feature of PD23,24. Cognitive

impairment takes several forms, starting with executive (frontal lobe) dysfunction and apathy, and progressing to visual hallucinations, impulse control, and behavioral disorders, and,

eventually, a global cognitive impairment that resembles AD24,25,26,27. As with the classical dementias, it is likely that the specific pattern of disease spread in PD may shape the clinical

features of the dementing process1. Here we take advantage of follow-up imaging in a large longitudinal, multicenter study that aimed to assess the progression of clinical features,

imaging, and other biomarkers in a cohort of well-characterized de novo PD patients. Using the PPMI database (http://www.ppmi-info.org/data), we test the theory that disease progression at a

1-year follow-up, as assessed by cortical thickness, is best explained by spread from a “disease reservoir” along the brain’s connectome. We further identify certain factors at baseline

that predict severity of disease progression at follow-up and provide evidence that cortical involvement might herald the onset of cognitive decline. RESULTS PARTICIPANTS Participants in the

present study consisted of de novo PD patients (_n_ = 105) and healthy controls (HC, _n_ = 57) enrolled in the PPMI study8 [accessed May 01, 2016), clinicaltrials.gov identifier

NCT01141023] who had follow-up assessment 1-year postenrollment (1.05 years for both groups). Demographic and clinical data are summarized in Table 1. PD patients were unmedicated at the

first evaluation. At 1-year follow-up, 81 of the 105 were on PD medication (e.g., levodopa or dopamine agonist). Among PD patients, performance on the Montreal Cognitive Assessment (MoCA)

was slightly worse at follow-up compared to baseline (from 27.5 to 26.7, _t_(103) = 3.038, _p_ = 0.003). PD patients also scored higher on the Movement Disorders Society Unified Parkinson’s

Disease Rating Scale Part III (UPDRS-III) at follow-up compared to baseline, but this did not reach statistical significance (from 21.7 to 23.1, _t_(95) = −1.544, _p_ = 0.126). CHANGES IN

CORTICAL THICKNESS AT 1-YEAR FOLLOW-UP Over the 1-year period, mean whole-brain cortical thickness significantly decreased among both PD patients (_t_ 1 = 3.055 mm ± 0.013; _t_ 2 = 3.027 mm 

± 0.013) and controls (_t_ 1 = 3.055 mm ± 0.016; _t_ 2 = 3.036 mm ± 0.017) (both _p_ < 0.01). Regional cortical thinning was greater in PD than HC (Fig. 1). Specifically, four significant

clusters were identified: (1) left occipital and inferior and middle temporal gyri, (2) left frontal, (3) right frontal, and (4) right pre- and post-central gyrus and supramarginal gyrus

(Fig. 2). The frontal lobe atrophy mostly involved the inferior parts of lateral and medial frontal cortex, including the orbitofrontal cortex, frontal pole, ventromedial, and rostral

dorsomedial prefrontal cortex. On the lateral surface, inferior and medial frontal gyri were implicated, with extension to the superior frontal gyrus on the left. In order to assess the

likely functional consequences of cortical thinning in PD, we assigned each cortical parcel to one of the seven intrinsic brain networks derived from resting-state functionl MRI (fMRI) by

Yeo et al.28. Cortical regions belonging to the limbic, frontoparietal, and ventral attention networks on average demonstrated the greatest cortical thinning, as well as anterior parts of

the default mode network (DMN; Fig. 3). There were no areas showing greater loss of cortical thickness in HC than in PD. CONNECTIVITY AND DISEASE PROPAGATION The disease propagation

hypothesis predicts that cortical brain atrophy at follow-up should depend on connectivity of cortical areas to affected areas at baseline. We had previously shown predominant subcortical

atrophy in de novo PD in this population. Connectivity between any cortical area and this putative subcortical “disease reservoir” should therefore predict the degree of cortical thinning

relative to age-matched controls when comparing follow-up to baseline. We define the disease reservoir as the subcortical areas showing atrophy in PD patients relative to controls at the

initial evaluation, as determined previously using DBM9. For the current analysis, the brain was segmented into 463 cortical and subcortical regions of interest29. Each surface vertex from

the cortical thickness images was mapped to the nearest volumetric parcel (see Methods section). To examine the propagation hypothesis, the connectivity (Conn _ij_ ) between any cortical

parcel _i_ and any parcel belonging to the disease reservoir _j_ was measured. Here we calculated Conn _ij_ based on healthy brain connectomes (see Methods section). The analysis was

performed for the whole brain and for each hemisphere separately using Spearman’s correlation to examine the relation between disease progression and exposure. Disease exposure was defined

as the product of connectivity and atrophy measure in the corresponding disease reservoir region _j_ at the onset of the disease, summed over all possible connections. $${\mathrm

Disease}\;{\mathrm exposure} \left( i \right) = \mathop {\sum} \limits_j {\mathrm Conn}_{ij} \cdot {\mathrm Atrophy}(j)$$ Both functional and structural connectivity were used to test the

propagation hypothesis (Fig. 4). Connectomes were generated using data from young healthy individuals to calculate the Conn _ij_ . Atrophy(_j_) was set equal to the _z_-value at each region

_j_ of the independent component demonstrating atrophy in PD compared to HC from the previously computed DBM map from the baseline MRI9. Note that for areas not showing significant atrophy

at baseline the value of Atrophy(_j_) is equal to 0. First, a functional connectome was generated using resting-state fMRI to define Conn _ij_ . Significant correlation was observed between

regional cortical thinning and disease exposure (_r_ = 0.19, _p_ < 0.0001). Significance was confirmed by permutation testing where we kept the connectivity structure (i.e.,

subcortical–cortical connections) intact and permuted the cortical thinning values (_n_ = 10,000) and the 95% confidence interval was measured using bootstrapping (_n_ = 10,000) (_r_ ∈ 

[0.11,0.28], _p__perm < 0.0001). Moreover, cortical thickness of both the left (_r_ = 0.22, _p_ = 0.0005, _r_ ∈ [0.09,0.34], _p__perm = 0.0004) and right hemisphere (_r_ = 0.27, _p_ < 

0.0001, _r_ ∈ [0.16,0.38], _p__perm < 0.0001) significantly correlated with the disease exposure from bilateral affected regions. Second, the same hypothesis was tested using structural

connectivity as measured by diffusion-weighted MRI (DW-MRI). Though the correlation was not significant at the whole-brain level (_r_ = −0.03, _p_ = 0.239, _r_ ∈ [−0.12,0.06], _p__perm = 

0.762), change in cortical thickness and disease exposure within the individual hemispheres were significantly correlated. Similar results were observed within the left hemisphere (_r_ = 

0.14, _p_ = 0.016, _r_ ∈ [0.01,0.27], _p__perm = 0.016) and the right hemisphere (_r_ = 0.11, _p_ = 0.044, _r_ ∈ [−0.02,0.25], _p__perm = 0.041). Failure to capture the relationship at the

whole-brain level may be due to inherent limitations in DW-MRI to detect interhemispheric connectivity30,31,32,33 (See Discussion section). For comparison, when only considering pure

Euclidean distance, we found that cortical thickness of neither the whole brain (_r_ = −0.06, _p_ = 0.104, _r_ ∈ [−0.15,0.03], _p__perm = 0.105) nor the left hemisphere (_r_ = −0.09, _p_ = 

0.088, _r_ ∈ [−0.22,0.04], _p__perm = 0.084) correlated with the disease exposure from bilateral-affected regions. However, when only considering the right hemisphere, we observed

significant results (_r_ = −0.21, _p_ = 0.0007, _r_ ∈ [−0.33,−0.09], _p__perm = 0.0006). Furthermore, we tested the relative contribution of functional and structural networks to disease

spread using a multivariate model for each hemisphere. Within the left hemisphere, the analysis shows that both functional connectivity (_β_ = 0.18, _t_(220) = −3.88, _p_ < 0.001) and

structural connectivity (_β_ = 0.04, _t_(220) = −2.27, _p_ = 0.02) significantly relate to cortical atrophy. Within the right hemisphere, again functional connectivity (_β_ = 0.12, _t_(220) 

= −2.56, _p_ = 0.01) and structural connectivity (β = 0.08, _t_(220) = −2.77, _p_ < 0.01) significantly relate to cortical atrophy. These results suggest that both types of connectivity

play a role in the disease-spread model; however, functional connectivity appears to be more significant. As a follow-up exploratory analysis, we tested whether striatal and non-striatal

regions within our disease reservoir may contribute differently to disease exposure. Based on functional connectivity, analysis restricted to a basal ganglia reservoir showed an increased

effect size within the left hemisphere (_r_ = 0.29, _p_ < 0.0001, _r_ ∈ [0.16,0.40], _p__perm < 0.0001), right hemisphere (_r_ = 0.37, _p_ < 0.0001, _r_ ∈ [0.26,0.47], _p__perm <

 0.0001), and across the whole rain (_r_ = 0.29, _p_ < 0.0001, _r_ ∈ [0.21,0.37], _p__perm < 0.0001) compared to the analysis using all affected subcortical regions. When repeating

this analysis for structural connectivity, the effect size observed for a basal ganglia-restricted analysis within the left hemisphere (_r_ = 0.14, _p_ = 0.016, _r_ ∈ [0.01,0.27], _p__perm =

 0.017), right hemisphere (_r_ = 0.15, _p_ = 0.012, _r_ ∈ [0.01,0.28], _p__perm = 0.011), and across the whole brain (_r_ = 0.02, _p_ = 0.35, _r_ ∈ [−0.08,0.11], _p__perm = 0.35) were

similar to the analyses considering all subcortical structures. RELATION TO COGNITION AND CEREBROSPINAL FLUID (CSF) MEASURES Additional analyses restricted to the PD group were performed to

investigate the relationship between cortical thinning of the four identified clusters and other disease-related measures (Fig. 5). To test whether CSF measures may contribute to the

cortical thinning within each of the clusters, we conducted partial correlations (controlling for age and gender). We found that changes in the left frontal (_r_ = −0.16, _p_ = 0.05) and

left occipital cluster (_r_ = −0.18, _p_ = 0.037) both correlated negatively with CSF Aβ42. Changes of the left occipital cluster (_r_ = 0.17, _p = _0.043) also correlated with CSF α-syn. We

then performed a linear regression model to examine whether cortical thinning was related to cognitive decline (as assessed by MoCA) and found the overall model, controlling for covariates,

was significant (_r_ _2_ _ = _0.12, _p_ = 0.03). Further post-hoc analysis revealed that only the cortical thinning of the left frontal cluster (_r_ = −0.19, _p_ = 0.03) significantly

related to MoCA changes. Note that these analyses fail to survive correction for multiple comparisons of the four significant clusters. CSF measures did not correlate with individual time

point measures nor changes in MoCA or motor severity as assessed by UPDRS-III (all _p_ > 0.05). No significant correlation was observed between cortical thickness of any of the four

clusters and measures of UPDRS-III or p-tau181. In the control group, no significant correlations between cortical thickness and measures of UPDRS-III, MoCA, or CSF biomarkers were observed.

DISCUSSION The present study investigated 105 PD and 57 control participants over a 1-year period following initial diagnosis. Early PD progression was associated with significant cortical

changes. While aging was also associated with cortical thinning, PD patients demonstrated significantly greater reduction in cortical thickness than controls (0.028 mm vs 0.019 mm reduction

on average). This was observed mostly in parts of the left occipital and bilateral frontal lobes and right somatomotor-sensory cortex—regions largely belonging to the limbic, frontoparietal,

ventral attention, and default mode networks. We tested the theory that neurodegeneration in PD results from a disease-spreading process that enters the cerebral cortex from initially

affected subcortical areas. Per the network propagation hypothesis, neurodegenerative diseases result from the aggregation and propagation of misfolded proteins, which in turn results in

neuronal death and brain atrophy2,4. The distribution of focal cortical thinning should therefore be contingent on connectivity to already affected areas, which we call a “disease reservoir”

by analogy to the spread of epidemics in populations. In the present study, we tested this hypothesis in PD with longitudinal data. We found that cortical regions with greater connectivity

(measured functionally or structurally) to the subcortical disease reservoir identified at the time of diagnosis demonstrated greater cortical atrophy over the 1-year period. This is in line

with our previous work using DBM in the initial scan from this dataset9, which indicated that regions demonstrating the greatest atrophy were those belonging to a connectivity network with

an epicenter in substantia nigra—the hypothesized origin of misfolded protein spread to the supratentorial cerebrum in PD. Some aspects of our model deserve mention. Predictions based on

physical distance were similar to topographical distance based on DW-MRI in the right hemisphere. This may be a result of the known distance bias of probabilistic tractography31 or that the

brain is organized such that neuronal connectivity correlates with Euclidean distance32. Resting-state connectivity, however, better predicted atrophy than tractography or topographical

distance, thus supporting the notion that connectivity improves our measure of propagation above and beyond consideration of geometrical distance alone. Moreover, there is no described

mechanism by which a protein spreads physically through the neuropil. While functional connectivity models fit the data at a whole-brain level, we only observed a relationship between

connectivity and thinning in individual hemispheres when using a DW-MRI connectome. This is perhaps due to inherent limitations in DW-MRI to capture interhemispheric connections33. Finally,

the disease exposure analysis limited to the basal ganglia only yielded greater effects on cortical thinning than the whole subcortical reservoir model. This might suggest a greater

contribution of basal ganglia to cortical propagation, although our results do not unequivocally demonstrate this. Mathematical models of the dynamic spread through brain networks,

accounting for factors such as protein accumulation and clearance, have been proposed in recent years for AD34,35,36. Such models might be usefully applied to PD. In our post-hoc analysis,

greater reduction in cortical thickness of the left frontal cluster in PD over the 1-year period was associated with worsening cognition (as assessed by the MoCA), although the _p_-value of

0.03 would not survive correction for multiple comparisons. This mirrors results from other studies that have similarly found trending results22,37. While dementia is rare in early PD,

especially as our cohort was selected to exclude dementia at baseline, we still see greater cortical thinning of frontal clusters in the PD group 1-year postdiagnosis. Cortical neuronal

degeneration could be occurring earlier in the clinical course of PD than once thought. Although PD is primarily thought of as a movement disorder, it is recognized to be a brain-wide

neurodegenerative process that spreads up from brainstem into cortex, as originally suggested by Braak38. In this dataset, only cortical thinning of the left frontal cluster related to

cognitive decline. Furthermore, thinning in this frontal cluster was weakly related to reduced Aβ42 at baseline (Fig. 5). Lower CSF Aβ42 and α-syn was also associated with occipital

thinning, which did not correlate with cognitive decline. The frontal clusters demonstrating greater loss of cortical thickness in PD involved lateral and medial frontal cortex,

predominantly in ventral and orbitofrontal areas. These regions of frontal lobe belong to the limbic system, the frontoparietal systems implicated in focused attention, and the DMN

implicated in memory retrieval. This pattern of frontal atrophy is similar to that described in a subgroup of AD patients identified by hierarchical clustering39 and in patients with

behavioral variant frontotemporal dementia1. A notable difference between the pattern of cortical atrophy described here and that commonly seen in AD is the relative absence of more

posterior involvement (precuneus, posterior cingulate, and parietal lobe) in our PD cohort. The effect seen here is much more anterior, possibly reflecting the propagation of disease via

dopaminergic or frontostriatal projections. Nonetheless, the DMN is also disrupted in AD40,41 and a target of amyloid plaques42. The anterior DMN may be an avenue through which dementia

develops in PD. We note also that, in the Braak model of AD, the initial cortical site of amyloid-β involvement is the ventral frontal cortex43 and that the observed ventral frontal atrophy

pattern in PD shows a high degree of similarity to early Aβ42 deposition measured by positron emission tomography in AD44. While cognitive impairment is a recognized consequence of PD23, its

underlying mechanism is not well understood. Our findings suggest that, in the PD group, higher Aβ42 deposition in the brain at baseline (indicated by lower CSF Aβ42 levels45) predicted

greater cortical thinning of the left frontal cluster over the 1-year period. Though our results are not conclusive, they suggest that Aβ42 may make the PD brain more susceptible to cortical

thinning and, in turn, more prone to cognitive decline. Indeed, recent positron emission tomographic studies have revealed that cortical amyloid deposits are associated, in a dose-dependent

manner, with increased risk for cognitive decline in PD46. Additionally, low CSF Aβ42 has previously been associated with the risk of subsequent dementia in PD47. These findings provide

evidence in favor of potentially overlapping mechanisms for dementia in PD and AD. Previous research found that 28% of PD patients have AD pathological features at postmortem48. While it can

be argued that the findings observed in the present study may be driven by coincidental occurrence of AD, we found that only 5 of the 105 PD subjects met a proposed CSF amyloid-positive

criterion49. A possible alternate explanation is that synergy between Aβ42 and α-syn deposition promotes cognitive impairment in PD50. Aβ42 has been found to promote gray matter atrophy in

distributed regions of the parietal and frontal lobe in the normal aging brain51, akin to the atrophy pattern observed in the present study; α-syn may accelerate the neurodegenerative

process associated with normal age-related Aβ42 deposition, especially in regions connected to the PD disease reservoir. The pattern of cortical atrophy in these de novo PD patients may

account for the constellation of mood and cognitive deficits that eventually characterize PD dementia. Involvement of the lateral prefrontal cortex and striatum could lead to executive

dysfunction. Inferior frontal areas also affected in behavioral variant frontotemporal dementia may be implicated in apathy, compulsive behaviors, changes in personality, and social

deficits. Finally, DMN involvement could account for Alzheimer-like memory deficits while degeneration of ventral visual stream areas could explain visual hallucinations. Using the

Neurosynth program and database52, we searched for the terms most often associated in fMRI publications with each anatomical cluster demonstrating atrophy in PD (Fig. 6). The left

occipital–temporal cluster is associated with (visual) motion, observation, perception, and vision. The right somatomotor-sensory cluster is associated with movement. The frontal clusters

are associated with expressions related to language and to the terms theory of mind, self, social, mentalizing, and memory. The atrophic regions observed here subserve the cognitive

functions that fail in PD dementia. Findings from this study should be considered in light of some limitations. While the cellular events leading to changes in cortical thickness are not

clear, cortical thickness as measured with T1-weighted MRI is nonetheless believed to reflect neuronal density along with glial support and dendritic arborization53,54. In PD, reductions in

cortical thickness most likely reflect neurodegeneration. In the current paper, we address the “protein-spread” model, which we note remains controversial and lacking in supportive human

evidence6,7. With our current dataset and imaging modalities, we cannot rule out the possibility of a stochastic attack on the brain as a whole that would give the appearance of network

spread. Areas that are connected could share features (e.g., gene expression55) that make them similarly vulnerable to a specific pathophysiological process. In this case, the map of atrophy

would look like a network but not be driven by connectivity. Additionally, neurodegeneration could propagate due to loss of trophic factors having nothing to do with misfolded protein

spread. Moreover, we cannot measure α-syn accumulation but only its consequence (cortical thinning), and the two may have a complex and nonlinear relationship. However, the present finding

that a connectivity model explains disease progression better than chance lends support to the network-propagation hypothesis. The inclusion criterion of MoCA > 26 may have biased both

the PD and HC pool. Also, while our sample is relatively large, the PPMI sample is a research-based cohort and may not be truly representative of a population-based cohort. Finally, most

patients were prescribed dopaminergic treatments after diagnosis, which may mask the magnitude of cognitive decline and motor symptoms and their relationship with brain atrophy. Future

longitudinal studies and follow-up with patients at later stages of the disease as well as postmortem confirmation are needed to understand the relationship between cortical atrophy and the

underlying disease process. The current study is the largest to date to assess longitudinal measures of cortical thickness using 3 T MRI in a de novo PD cohort. Our findings suggest that

disease propagation to the cortex occurs earlier than previously thought. Focal cortical thinning may be contingent on connectivity to a disease reservoir: the greater the connections, the

higher the disease exposure, and the more the cortical atrophy. PD progression involves subcortical α-syn spreading through intrinsic brain networks to the cortex, possibly in synergy with

Aβ42, which may be the harbinger of cognitive decline. METHODS The PD patients included in this study met the PPMI inclusion criteria8 consisting of asymmetric resting tremor or asymmetric

bradykinesia or any two of bradykinesia, resting tremor, and limb rigidity. Patients had to have been diagnosed within 2 years of enrollment and not have received anti-parkinsonian

medications. Diagnosis was confirmed with dopamine transporter imaging using single-photon emission computed tomography. Control subjects were free of neurological disease. Each

participating PPMI site received approval from their local institutional review board and obtained written informed consent from all subjects. Analysis of this dataset was approved by the

Montreal Neurological Institute Research Ethics Board. Only participants with 3 T MRI data at both initial visit and 1-year follow-up were included in the present study. Of the 237 PD

patients with 3 T acquisition at baseline, 122 had MRI data at the 1-year follow-up. For controls, 118 subjects had 3 T acquisition at baseline, and 62 had MRI data at the 1-year follow-up.

All non-categorical data were normally distributed as per Shapiro–Wilk tests. Group difference in sociodemographic and neuropsychological variables were analyzed using the χ2 test for

categorical data and _t_-test for normally distributed data. MRI ACQUISITION T1-weighted MRI scans were acquired in the sagittal plane on 3 T scanners at each study site using a

magnetization-prepared rapid-acquisition gradient echo sequence (full study protocol: http://www.ppmi-info.org/study-design/research-documents-and-sops/). Acquisition parameters were as

follows: repetition time = 2300/1900 ms; echo time = 2.98/2.96/2.27/2.48/2.52 ms; inversion time = 900 ms; flip angle: 9°; 256 × 256 matrix; and 1 × 1 × 1 mm3 isotropic voxel. Acquisition

site was used as a covariate in all analyses. CORTICAL THICKNESS Cortical thickness models were generated from the T1-weighted 3 T MRI scans using the CIVET2.1 preprocessing pipeline56. A

summary of the steps involved follows; full documentation can be found online (http://www.bic.mni.mcgill.ca/ServicesSoftware/CIVET). The T1-weighted MRIs were linearly registered to the

MNI-ICBM152 volumetric template57,58. Images were then corrected for signal intensity non-uniformity59, and a brain mask was calculated from each input image. Images were then segmented into

gray matter, white matter, CSF, and background56. After the CIVET pipeline, quality control was carried out by two independent reviewers using previously described criteria60 to ensure

adequate quality of the T1-wieghted volume images, such as linear and nonlinear registration, and to exclude distortions in gray and white matter surfaces and motion artifacts. Seventeen of

the 122 PD and 5 of the 62 HC failed quality control. The distance between the corresponding vertices of inner and outer cortical surfaces were evaluated by a fully automated method

(Constrained Laplacian Anatomic Segmentation using Proximity algorithm61) to provide a measure of cortical thickness at each vertex. The white and gray matter interface was fitted across

40,692 vertices in each hemisphere to create the inner cortical surface, then expanded to fit the gray matter and CSF interface to create the outer cortical surface. The surfaces of each

subject were nonlinearly registered to an average surface created from the volumetric ICBM152 template. In order to test the propagation hypothesis, the brain was segmented into 448 cortical

and 15 subcortical regions29, and each of the 40,692 cortical vertices from the cortical thickness analysis was interpolated and assigned to one of the 448 cortical parcels. Changes in

cortical thickness were calculated by subtracting the values (Δ_t = t_ 1 _−t_ 2) at the 1-year follow-up (_t_ 2) from the baseline (_t_ 1) for both groups using the SurfStat software package

(http://www.math.mcgill.ca/keith/surfstat/). Difference between PD and controls [(PD1-PD2)–(HC1-HC2)] were analyzed statistically based on random field theory with a threshold of _p_ < 

0.05 corrected for multiple comparisons over the entire cortical surface62. In order to address the relationship between cognitive decline and cortical thinning over the 1-year period, we

performed a linear regression model. Change in MoCA was the dependent variable, while mean cortical thinning values from each of the four clusters distinguishing PD from controls were the

independent variables. Gender and age at baseline were used as covariates. Observed differences in cortical thickness were then correlated to other potential disease-related measures (i.e.,

UPDRS-III, MoCA, CSF α-syn, and tau and Aβ42 markers) using Pearson partial correlations as post-hoc analysis. All the analyses were performed using MatLab 2015b statistics toolbox. Within

PD patients, the difference in cognitive performance over the 1-year interval was examined using a linear regression model with MoCA as the dependent variable and the thinning of the four

cortical clusters as independent variables. To test whether cortical thickness of a specific cluster was linked to cognitive performance, post-hoc analyses using partial correlations were

conducted. Age at baseline and gender were used as confounding covariates for all analyses including for partial correlations, linear regression model, and for obtaining cortical thickness.

DEFORMATION-BASED MORPHOMETRY Data were analyzed from the baseline visits of newly diagnosed PD patients (_n_ = 232) and an age-matched control group (_n_ = 117) obtained from the PPMI

database. All subjects’ initial visit's 3 T high-resolution T1-weighted MRI scans underwent preprocessing steps including denoising63, intensity non-uniformity correction59, and

intensity range normalization. Next, each subject’s T1-weighted image was first linearly, and then nonlinearly, registered to the MNI-ICBM152 template. Using the resulting nonlinear

transformation fields, we calculated deformation morphometry maps (i.e., the determinant of the Jacobian matrix of transformation). DBM maps were then concatenated and independent component

analysis (ICA) was employed to define independent sources of deformation in the brain. For the resulting ICA components (_n_ = 30), the PD group was compared to healthy controls (unpaired

_t_-test). The atrophy level significantly differed between PD and control groups in only one of these components after correcting for multiple comparisons (Bonferroni correction),

suggesting that this component captures PD-specific alterations in the brain. We used this component (referred to as PD-ICA here) as our atrophy map. The component score (atrophy measure) in

PD-ICA for each subject was significantly correlated with disease severity as measured by UPDRS-III and single-photon emission computed tomography striatum binding ratio. Please see

Zeighami et al.9 for further details. CONNECTIVITY AND DISEASE PROPAGATION The difference in cortical thinning in PD vs controls was used to define disease progression. DBM on the baseline

MRI was employed to assess disease-related brain atrophy and to define a source, or reservoir, for disease propagation. Based on a parcellation that is a subdivision of the common FreeSurfer

implementation of the Desikan–Killiany Atlas, gray matter was segmented into 463 regions of interest of approximately similar size29. Of these parcels, 448 were cortical and 15 were

subcortical (thalamus, caudate, putamen, pallidum, accumbens, hippocampus, amygdala in each hemisphere and brainstem). Each surface vertex from the cortical thickness analysis was then

mapped to the corresponding nearest whole-brain parcel from the Atlas. The “disease reservoir” was defined as the 15 subcortical parcels from the Atlas, and the Atrophy value at each parcel

was obtained from the DBM map. To generate brain networks and connectomes, data-driven brain parcellations from young healthy individuals were used29. The functional connectivity map was

derived from resting-state fMRI data of 40 young healthy subjects (25.3 ± 4.9 years, 24 males)64. The scans were done using Siemens Medical Trio 3 T MRI scanner, consisting of a

high-resolution T1-weighted image, as well as T2*-weighted images with blood-oxygen-level-dependent contrast (3.3 mm3 isotropic voxels, TR 1.9 s). The dataset was preprocessed using

conventional steps consisting of rigid body motion correction, correcting for nuisance variables (white matter, CSF, motion parameters), and low pass filtering. The structural connectivity

map was generated using the Illinois Institute of Technology Human Brain Atlas v.3, (freely available at http://www.iit.edu/~mri/Home.html). The atlas was constructed using data from 72

neurologically healthy subjects (20–40 years, 28.8 ± 5.5 years, 30 males). All subjects were scanned using a 3 T General Electric MRI scanner (GE, Waukesha, WI) with the following parameters

(TR = 5800 ms, TE = 94ms, voxel size = 3 mm3, 256 × 256 final image matrix, FoV = 24 × 24 cm2, _b_ = 900 s mm−2, 45 slices for 12 diffusion directions uniformly distributed in 3D space,

scan time: 21 min 57 s)65. Imposing a maximum of 500 mm trace length and a curvature threshold of ±90, the anatomical connection probability was calculated between all pairs of regions using

a probabilistic graph-based tractography algorithm66. Spearman’s correlation was used to investigate the relationship between disease progression in cortical areas (cortical thinning) and

the ‘disease exposure’ for each region in the PD-ICA (defined by the average values of each region weighted by connectivity). The same procedure was performed to test propagation hypotheses

along the functional and structural connectomes. The propagation analysis was performed for the whole brain and for each hemisphere separately. CSF MEASURES CSF Aβ42, α-syn, t-tau, and

p-tau181 were collected by lumbar puncture at PPMI clinical sites67. Frozen aliquots (at −80 °C) were then shipped and assessed using the MAP Luminex platform (Luminex Corp, Austin, TX) with

Innogenetics (INNO-BIA AlzBio3; Ghent, Belgium; for research use-only reagents) immunoassays used for Aβ42, t-tau and p-tau181. CSF α-syn assay was performed at Covance using a commercially

available ELISA Assay Kit (Covance, Dedham, MA)68. To evaluate the possible contamination of blood in CSF—a factor thought to influence the level of some proteins including α-syn and Aβ42

69,70—the relationship between CSF hemoglobin and other CSF measures was assessed. No significant association was observed (all _p_ < 0.05), suggesting that there is no effect of added

blood on concentration of CSF measures. Owing to the ongoing nature of the study and incomplete data, only CSF data from baseline were used in this analysis68. CODE AVAILABILITY All code

written by the authors of this manuscript is available upon request. DATA AVAILABILITY The DBM maps are available at https://neurovault.org/collections/860/. Cortical thickness maps will be

made available on request. The CIVET pipeline is available at http://www.bic.mni.mcgill.ca/ServicesSoftware/CIVET. The functional and structural connectome data are from cited sources and

available upon request from the authors of the related publications. REFERENCES * Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L. & Greicius, M. D. Neurodegenerative diseases

target large-scale human brain networks. _Neuron_ 62, 42–52 (2009). Article  CAS  PubMed  PubMed Central  Google Scholar  * Warren, J. D. et al. Molecular nexopathies: a new paradigm of

neurodegenerative disease. _Trends Neurosci._ 36, 561–569 (2013). Article  CAS  PubMed  PubMed Central  Google Scholar  * Burre, J. et al. Alpha-synuclein promotes SNARE-complex assembly in

vivo and in vitro. _Science_ 329, 1663–1667 (2010). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Luk, K. C. et al. Pathological alpha-synuclein transmission initiates

Parkinson-like neurodegeneration in nontransgenic mice. _Science_ 338, 949–953 (2012). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Braak, H. et al. Staging of brain

pathology related to sporadic Parkinson’s disease. _Neurobiol. Aging_ 24, 197–211 (2003). Article  PubMed  Google Scholar  * Walsh, D. M. & Selkoe, D. J. A critical appraisal of the

pathogenic protein spread hypothesis of neurodegeneration. _Nat. Rev. Neurosci._ 17, 251–260 (2016). Article  CAS  PubMed  Google Scholar  * Surmeier, D. J., Obeso, J. A. & Halliday, G.

M. Selective neuronal vulnerability in Parkinson disease. _Nat. Rev. Neurosci._ 18, 101–113 (2017). Article  CAS  PubMed  PubMed Central  Google Scholar  * Marek, K. et al. The Parkinson

Progression Marker Initiative (PPMI). _Prog. Neurobiol._ 95, 629–635 (2011). Article  Google Scholar  * Zeighami, Y. et al. Network structure of brain atrophy in de novo Parkinson’s disease.

_eLife_ 4, e08440 (2015). Article  PubMed Central  Google Scholar  * Reetz, K. et al. Structural imaging in the presymptomatic stage of genetically determined parkinsonism. _Neurobiol.

Dis._ 39, 402–408 (2010). Article  CAS  PubMed  Google Scholar  * Fioravanti, V. et al. MRI correlates of Parkinson’s disease progression: a voxel based morphometry study. _Parkinsons Dis._

2015, 8 (2015). Google Scholar  * Camicioli, R. et al. Voxel-based morphometry reveals extra-nigral atrophy patterns associated with dopamine refractory cognitive and motor impairment in

parkinsonism. _Parkinsonism Relat. Disord._ 15, 187–195 (2009). Article  PubMed  Google Scholar  * Menke, R. A. et al. Comprehensive morphometry of subcortical grey matter structures in

early‐stage Parkinson’s disease. _Hum. Brain Mapp._ 35, 1681–1690 (2014). Article  PubMed  Google Scholar  * Davatzikos, C. Why voxel-based morphometric analysis should be used with great

caution when characterizing group differences. _Neuroimage_ 23, 17–20 (2004). Article  PubMed  Google Scholar  * Hutton, C., Draganski, B., Ashburner, J. & Weiskopf, N. A comparison

between voxel-based cortical thickness and voxel-based morphometry in normal aging. _Neuroimage_ 48, 371–380 (2009). Article  PubMed  PubMed Central  Google Scholar  * Pereira, J. B. et al.

Assessment of cortical degeneration in patients with Parkinson’s disease by voxel-based morphometry, cortical folding, and cortical thickness. _Hum. Brain Mapp._ 33, 2521–2534 (2012).

Article  PubMed  Google Scholar  * Khundrakpam, B. S., Lewis, J. D., Kostopoulos, P., Carbonell, F. & Evans, A. C. Cortical thickness abnormalities in autism spectrum disorders through

late childhood, adolescence, and adulthood: a large-scale MRI study. _Cereb. Cortex_ 27, 1721–1731 (2017). Article  PubMed  Google Scholar  * Dickerson, B. C. et al. The cortical signature

of Alzheimer’s disease: regionally specific cortical thinning relates to symptom severity in very mild to mild AD dementia and is detectable in asymptomatic amyloid-positive individuals.

_Cereb. Cortex_ 19, 497–510 (2009). Article  PubMed  Google Scholar  * Ibarretxe-Bilbao, N. et al. Progression of cortical thinning in early Parkinson’s disease. _Mov. Disord._ 27, 1746–1753

(2012). Article  PubMed  Google Scholar  * Jubault, T. et al. Patterns of cortical thickness and surface area in early Parkinson’s disease. _Neuroimage_ 55, 462–467 (2011). Article  PubMed

  Google Scholar  * Zarei, M. et al. Cortical thinning is associated with disease stages and dementia in Parkinson’s disease. _J. Neurol. Neurosurg. Psychiatry_ 84, 875–882 (2013). Article 

PubMed  PubMed Central  Google Scholar  * Pereira, J. B. et al. Initial cognitive decline is associated with cortical thinning in early Parkinson disease. _Neurology_ 82, 2017–2025 (2014).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Hely, M. A., Reid, W. G., Adena, M. A., Halliday, G. M. & Morris, J. G. The Sydney multicenter study of Parkinson’s disease: the

inevitability of dementia at 20 years. _Mov. Disord._ 23, 837–844 (2008). Article  PubMed  Google Scholar  * Aarsland, D. et al. Cognitive decline in Parkinson disease. _Nat. Rev. Neurol._

13, 217–231 (2017). Article  PubMed  Google Scholar  * Dirnberger, G. & Jahanshahi, M. Executive dysfunction in Parkinson’s disease: a review. _J. Neuropsychol._ 7, 193–224 (2013).

Article  PubMed  Google Scholar  * Callesen, M. B., Weintraub, D., Damholdt, M. F. & Moller, A. Impulsive and compulsive behaviors among Danish patients with Parkinson’s disease:

prevalence, depression, and personality. _Parkinsonism Relat. Disord._ 20, 22–26 (2014). Article  CAS  PubMed  Google Scholar  * Ffytche, D. H. et al. The psychosis spectrum in Parkinson

disease. _Nat. Rev. Neurol._ 13, 81–95 (2017). Article  PubMed  Google Scholar  * Yeo, B. T. et al. The organization of the human cerebral cortex estimated by intrinsic functional

connectivity. _J. Neurophysiol._ 106, 1125–1165 (2011). Article  PubMed  Google Scholar  * Cammoun, L. et al. Mapping the human connectome at multiple scales with diffusion spectrum MRI. _J.

Neurosci. Methods_ 203, 386–397 (2012). Article  PubMed  Google Scholar  * Neher, P. F., Descoteaux, M., Houde, J.-C., Stieltjes, B. & Maier-Hein, K. H. Strengths and weaknesses of

state of the art fiber tractography pipelines—a comprehensive in-vivo and phantom evaluation study using tractometer. _Med. Image Anal._ 26, 287–305 (2015). Article  PubMed  Google Scholar 

* Jbabdi, S., Sotiropoulos, S. N., Haber, S. N., Van Essen, D. C. & Behrens, T. E. Measuring macroscopic brain connections in vivo. _Nat. Neurosci._ 18, 1546–1555 (2015). Article  CAS 

PubMed  Google Scholar  * Ercsey-Ravasz, M. et al. A predictive network model of cerebral cortical connectivity based on a distance rule. _Neuron_ 80, 184–197 (2013). Article  CAS  PubMed 

PubMed Central  Google Scholar  * Jones, D. K., Knosche, T. R. & Turner, R. White matter integrity, fiber count, and other fallacies: the do’s and don’ts of diffusion MRI. _Neuroimage_

73, 239–254 (2013). Article  PubMed  Google Scholar  * Raj, A., Kuceyeski, A. & Weiner, M. A network diffusion model of disease progression in dementia. _Neuron_ 73, 1204–1215 (2012).

Article  CAS  PubMed  PubMed Central  Google Scholar  * Iturria-Medina, Y., Sotero, R. C., Toussaint, P. J. & Evans, A. C. & Alzheimer’s Disease Neuroimaging Initiative. Epidemic

spreading model to characterize misfolded proteins propagation in aging and associated neurodegenerative disorders. _PLoS Comput. Biol._ 10, e1003956 (2014). Article  PubMed  PubMed Central

  Google Scholar  * Zhou, J., Gennatas, E. D., Kramer, J. H., Miller, B. L. & Seeley, W. W. Predicting regional neurodegeneration from the healthy brain functional connectome. _Neuron_

73, 1216–1227 (2012). Article  CAS  PubMed  PubMed Central  Google Scholar  * Caspell-Garcia, C. et al. Multiple modality biomarker prediction of cognitive impairment in prospectively

followed de novo Parkinson disease. _PLoS ONE_ 12, e0175674 (2017). Article  PubMed  PubMed Central  Google Scholar  * Braak, H. & Braak, E. Pathoanatomy of Parkinson’s disease. _J.

Neurol._ 247, II3–II10 (2000). Article  PubMed  Google Scholar  * Hwang, J. et al. Prediction of Alzheimer’s disease pathophysiology based on cortical thickness patterns. _Alzheimers Dement.

(Amst.)_ 2, 58–67 (2016). Google Scholar  * Greicius, M. D., Srivastava, G., Reiss, A. L. & Menon, V. Default-mode network activity distinguishes Alzheimer’s disease from healthy aging:

evidence from functional MRI. _Proc. Natl. Acad. Sci. USA_ 101, 4637–4642 (2004). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Koch, W. et al. Diagnostic power of default

mode network resting state fMRI in the detection of Alzheimer’s disease. _Neurobiol. Aging_ 33, 466–478 (2012). Article  PubMed  Google Scholar  * Sheline, Y. I. et al. Amyloid plaques

disrupt resting state default mode network connectivity in cognitively normal elderly. _Biol. Psychiatry_ 67, 584–587 (2010). Article  CAS  PubMed  Google Scholar  * Braak, H. & Braak,

E. Frequency of stages of Alzheimer-related lesions in different age categories. _Neurobiol. Aging_ 18, 351–357 (1997). Article  CAS  PubMed  Google Scholar  * Villain, N. et al. Regional

dynamics of amyloid-β deposition in healthy elderly, mild cognitive impairment and Alzheimer’s disease: a voxelwise PiB–PET longitudinal study. _Brain_ 135, 2126–2139 (2012). Article  PubMed

  Google Scholar  * Palmqvist, S. et al. Accuracy of brain amyloid detection in clinical practice using cerebrospinal fluid beta-amyloid 42: a cross-validation study against amyloid positron

emission tomography. _JAMA Neurol._ 71, 1282–1289 (2014). Article  PubMed  Google Scholar  * Gomperts, S. N. et al. PET radioligands reveal the basis of dementia in Parkinson’s disease and

dementia with Lewy bodies. _Neurodegener. Dis._ 16, 118–124 (2016). Article  CAS  PubMed  Google Scholar  * Alves, G. et al. CSF Aβ42 predicts early-onset dementia in Parkinson disease.

_Neurology_ 82, 1784–1790 (2014). Article  CAS  PubMed  Google Scholar  * Irwin, D. J. et al. Neuropathologic substrates of Parkinson disease dementia. _Ann. Neurol._ 72, 587–598 (2012).

Article  CAS  PubMed  PubMed Central  Google Scholar  * McMillan, C. T. & Wolk, D. A. Presence of cerebral amyloid modulates phenotype and pattern of neurodegeneration in early

Parkinson’s disease. _J. Neurol._ 87, 1112–1122 (2016). Google Scholar  * Frey, K. A. & Petrou, M. Imaging amyloidopathy in Parkinson disease and parkinsonian dementia syndromes. _Clin.

Transl. Imaging_ 3, 57–64 (2015). Article  PubMed  PubMed Central  Google Scholar  * Sepulcre, J. et al. In vivo tau, amyloid, and gray matter profiles in the aging brain. _J. Neurosci._ 36,

7364–7374 (2016). Article  CAS  PubMed  PubMed Central  Google Scholar  * Yarkoni, T., Poldrack, R. A., Nichols, T. E., Van Essen, D. C. & Wager, T. D. Large-scale automated synthesis

of human functional neuroimaging data. _Nat. Methods_ 8, 665–670 (2011). Article  CAS  PubMed  PubMed Central  Google Scholar  * Chklovskii, D. B., Mel, B. & Svoboda, K. Cortical

rewiring and information storage. _Nature_ 431, 782–788 (2004). Article  ADS  CAS  PubMed  Google Scholar  * la Fougère, C. et al. Where in-vivo imaging meets cytoarchitectonics: the

relationship between cortical thickness and neuronal density measured with high-resolution [18F]flumazenil-PET. _Neuroimage_ 56, 951–960 (2011). Article  PubMed  Google Scholar  * Richiardi,

J. et al. Correlated gene expression supports synchronous activity in brain networks. _Science_ 348, 1241–1244 (2015). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  *

Ad-Dab’bagh, Y. et al. The CIVET image-processing environment: a fully automated comprehensive pipeline for anatomical neuroimaging research. In _Proc. Organization for Human Brain Mapping_

(2006). * Fonov, V. S., Evans, A. C., McKinstry, R. C., Almli, C. & Collins, D. Unbiased nonlinear average age-appropriate brain templates from birth to adulthood. _Neuroimage_ 47, S102

(2009). Article  Google Scholar  * Fonov, V. et al. Unbiased average age-appropriate atlases for pediatric studies. _Neuroimage_ 54, 313–327 (2011). Article  PubMed  Google Scholar  * Sled,

J. G., Zijdenbos, A. P. & Evans, A. C. A nonparametric method for automatic correction of intensity nonuniformity in MRI data. _IEEE Trans. Med. Imaging_ 17, 87–97 (1998). Article  CAS 

PubMed  Google Scholar  * Simmons, A. et al. The AddNeuroMed framework for multi‐centre MRI assessment of Alzheimer’s disease: experience from the first 24 months. _Int. J. Geriatr.

Psychiatry_ 26, 75–82 (2011). Article  PubMed  Google Scholar  * Kim, J. S. et al. Automated 3-D extraction and evaluation of the inner and outer cortical surfaces using a Laplacian map and

partial volume effect classification. _Neuroimage_ 27, 210–221 (2005). Article  PubMed  Google Scholar  * Worsley, K. J. et al. A unified statistical approach for determining significant

signals in images of cerebral activation. _Hum. Brain Mapp._ 4, 58–73 (1996). Article  CAS  PubMed  Google Scholar  * Coupé, P. et al. An optimized blockwise nonlocal means denoising filter

for 3-D magnetic resonance images. _IEEE Trans. Med. Imaging_ 27, 425–441 (2008). Article  PubMed  PubMed Central  Google Scholar  * Mišić, B. et al. Cooperative and competitive spreading

dynamics on the human connectome. _Neuron_ 86, 1518–1529 (2015). Article  PubMed  Google Scholar  * Varentsova, A., Zhang, S. & Arfanakis, K. Development of a high angular resolution

diffusion imaging human brain template. _Neuroimage_ 91, 177–186 (2014). Article  PubMed  PubMed Central  Google Scholar  * Iturria-Medina, Y. et al. Characterizing brain anatomical

connections using diffusion weighted MRI and graph theory. _Neuroimage_ 36, 645–660 (2007). Article  CAS  PubMed  Google Scholar  * Kang, J., Irwin, D. J. & Chen-Plotkin, A. S. et al.

Association of cerebrospinal fluid β-amyloid 1-42, t-tau, p-tau181, and α-synuclein levels with clinical features of drug-naive patients with early parkinson disease. _JAMA Neurol._ 70,

1277–1287 (2013). PubMed  PubMed Central  Google Scholar  * Mollenhauer, B. et al. Total CSF alpha-synuclein is lower in de novo Parkinson patients than in healthy subjects. _Neurosci.

Lett._ 532, 44–48 (2013). Article  CAS  PubMed  Google Scholar  * Bjerke, M. et al. Confounding factors influencing amyloid beta concentration in cerebrospinal fluid. _Int. J. Alzheimers

Dis._ 2010, 986310 (2010). PubMed  PubMed Central  Google Scholar  * Hong, Z. et al. DJ-1 and α-synuclein in human cerebrospinal fluid as biomarkers of Parkinson’s disease. _Brain_ 133,

713–726 (2010). Article  PubMed  PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS PPMI—a public–private partnership—is funded by the Michael J. Fox Foundation for

Parkinson’s Research and funding partners, including AbbVie, Avid, Biogen, Bristol-Myers Squibb, Covance, GE Healthcare, Genentech, GlaxoSmithKline, Lilly, Lundbeck, Merck, Meso Scale

Discovery, Pfizer, Piramal, Roche, Sanofi Genzyme, Servier, Teva, and UCB. This research was supported by grants from the Michael J. Fox Foundation for Parkinson’s Research, the W. Garfield

Weston Foundation, and the Alzheimer’s Association; the Canadian Institutes of Health Research; and the Natural Sciences and Engineering Research Council of Canada to A.D. Y.Y. is a Vanier

Scholar and receives funding from the Canadian Institute of Health Research. Y.Z. holds the Jeanne Timmins Costello Fellowship. AUTHOR INFORMATION Author notes * Y. Yau and Y. Zeighami

contributed equally to this work. AUTHORS AND AFFILIATIONS * Montreal Neurological Institute, McGill University, 3801 University Street, Montreal, QC, H3A 2B4, Canada Y. Yau, Y. Zeighami, T.

E. Baker, K. Larcher, U. Vainik, M. Dadar, V. S. Fonov, B. Mišić, D. L. Collins & A. Dagher * Center for Molecular and Behavioral Neuroscience, Rutgers University, 197 University

Avenue, Newark, NJ, 07102, USA T. E. Baker * Institute of Psychology, Faculty of Social Sciences, University of Tartu, Näituse 2, 50409, Tartu, Estonia U. Vainik * Department of Radiology,

Lausanne University Hospital and University of Lausanne, Rue du Bugnon 21, 1011, Lausanne, Switzerland P. Hagmann * Brain Center Rudolf Magnus, UMC Utrecht, Heidelberglaan 100, A01.126, 3508

GA, Utrecht, The Netherlands A. Griffa Authors * Y. Yau View author publications You can also search for this author inPubMed Google Scholar * Y. Zeighami View author publications You can

also search for this author inPubMed Google Scholar * T. E. Baker View author publications You can also search for this author inPubMed Google Scholar * K. Larcher View author publications

You can also search for this author inPubMed Google Scholar * U. Vainik View author publications You can also search for this author inPubMed Google Scholar * M. Dadar View author

publications You can also search for this author inPubMed Google Scholar * V. S. Fonov View author publications You can also search for this author inPubMed Google Scholar * P. Hagmann View

author publications You can also search for this author inPubMed Google Scholar * A. Griffa View author publications You can also search for this author inPubMed Google Scholar * B. Mišić

View author publications You can also search for this author inPubMed Google Scholar * D. L. Collins View author publications You can also search for this author inPubMed Google Scholar * A.

Dagher View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS Y.Y., Y.Z., and A.D. contributed to the study design and interpretation. Y.Y.,

Y.Z., T.B., K.L., U.V., and M.D. analyzed the data. V.F., B.M., D.L.C., P.H., and A.G. contributed to methods. Y.Y. drafted the initial manuscript; all authors contributed to writing of this

manuscript. CORRESPONDING AUTHOR Correspondence to A. Dagher. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing financial interests. ADDITIONAL INFORMATION

PUBLISHER'S NOTE: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ELECTRONIC SUPPLEMENTARY MATERIAL PEER REVIEW

FILE RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and

reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes

were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If

material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain

permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS

ARTICLE Yau, Y., Zeighami, Y., Baker, T.E. _et al._ Network connectivity determines cortical thinning in early Parkinson’s disease progression. _Nat Commun_ 9, 12 (2018).

https://doi.org/10.1038/s41467-017-02416-0 Download citation * Received: 06 June 2017 * Accepted: 23 November 2017 * Published: 02 January 2018 * DOI:

https://doi.org/10.1038/s41467-017-02416-0 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not

currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative