Neuroimaging biomarkers differentiate Parkinson disease with and without cognitive impairment and dementia

Conor Owens-Walton, David Jakabek, Brian D. Power, Mark Walterfang, Sara Hall, Danielle van Westen, Jeffrey C.L. Looi, Marnie Shaw, Oskar Hansson
doi: https://doi.org/10.1101/19002147

Abstract

Mild cognitive impairment in Parkinson disease places a high burden on patients and is likely a precursor to Parkinson disease-related dementia. Studying the functional connectivity and morphology of subcortical structures within basal ganglia-thalamocortical circuits may uncover neuroimaging biomarkers of cognitive dysfunction in PD. We used an atlas-based seed region-of-interest approach to investigate resting-state functional connectivity of important subdivisions of the caudate nucleus, putamen and thalamus, between controls (n = 33), cognitively unimpaired Parkinson disease subjects (n = 33), Parkinson disease subjects with mild cognitive impairment (n = 22) and Parkinson disease subjects with dementia (n = 17). We then investigated how the morphology of the caudate, putamen and thalamus structures and differed between groups. Results indicate that cognitively unimpaired Parkinson disease subjects, compared to controls, display increased functional connectivity of the dorsal caudate, anterior putamen and mediodorsal thalamic subdivisions with areas across the frontal lobe, as well as reduced functional connectivity of the dorsal caudate with posterior cortical and cerebellar regions. Compared to cognitively unimpaired subjects, Parkinson disease subjects with mild cognitive impairment demonstrated reduced functional connectivity of the mediodorsal thalamus with midline nodes within the executive-control network. Compared to subjects with mild cognitive impairment, subjects with dementia demonstrated reduced functional connectivity of the mediodorsal thalamus with the posterior cingulate cortex, a key node within the default-mode network. Extensive volumetric and surface-based deflation was found in Parkinson disease subjects with dementia. Our research demonstrates how functional connectivity of the caudate, putamen and thalamus are implicated in the pathophysiology of cognitive impairment and dementia in Parkinson disease, with mild cognitive impairment and dementia in Parkinson disease associated with a breakdown in functional connectivity of the mediodorsal thalamus with para- and posterior cingulate regions of the brain.

1. Introduction

Parkinson disease (PD) is the second most common neurodegenerative disorder in the world, affecting 2-3% of the population over the age of 65 (Poewe et al., 2017). Traditionally categorized as a disorder of movement, it is now widely recognised that patients experience significant neuropsychiatric symptoms relating to executive, memory, attentional and visual disturbances (Litvan et al., 2012). Recent recognition of the high prevalence and substantial impact of cognitive impairment in PD has drawn research attention to PD-related mild cognitive impairment (Goldman et al., 2018) which is associated with a significantly higher likelihood of developing dementia (Broeders et al., 2013).

The neural processes that distinguish PD patients with cognitive impairment from those without it are poorly understood (Williams-Gray et al., 2007). The pathological hallmark of PD is the presence of α-synuclein-immunopositive Lewy bodies and neurites (Obeso et al., 2000) which results in a loss of dopaminergic neurons in the substantia nigra pars compacta, and subsequent depletion of dopamine at the striatum (caudate nucleus and putamen) (Kish et al., 1988). While nigrostriatal dopamine loss is a core feature of PD, evidence suggests that pathology within the thalamus also contributes to the abnormal neural activity associated with the disorder (Halliday, 2009). The combined neural activity in these nodes within basal ganglia-thalamocortical circuits in PD results in brain network abnormalities being as key components to PD pathophysiology (Strafella, 2013).

To understand the functioning of brain networks it is necessary to study both the constituent neuronal elements of networks, as well as their interconnections (Sporns et al., 2005). Investigating neuronal elements can be done by quantifying disease-related effects on the morphology of ‘hubs’ (Looi et al., 2014) while interconnections can be investigated via resting-state functional connectivity (FC) analyses (Damoiseaux et al., 2006). Due to the importance of the caudate nucleus, putamen and the thalamus to the abnormal neural activity associated with PD, coupled with the fact that these structures are considered important hubs in brain networks (Hwang et al., 2017, Looi and Walterfang, 2013), an analysis of the morphology and FC of these structures may yield neuroimaging biomarkers of brain network abnormalities in PD that relate to cognitive dysfunction.

There is a significant variability in the field of neuroimaging biomarkers in PD, with reviews of structural and functional studies revealing few clear patterns (Khan et al., 2018). Research in this field has yielded important information about how FC of widescale intrinsic connectivity networks is impacted in PD relative to cognitive status, however the specific role played by important hubs within basal ganglia-thalamocortical circuitry is yet to be fully elucidated. Studies have found cognitively unimpaired PD is associated with both increases (Gorges et al., 2015) and decreases (Bell et al., 2015) in FC of basal ganglia-thalamocortical circuitry, while PD subjects with MCI have been shown to generally display decreases in FC (Amboni et al., 2015, Gorges et al., 2015). Gorges et al. (2015) argue that the evidence suggests that PD may involve a transient process from the hyper- to the hypo-connected state which is linked to cognitive decline. Data on the morphological changes to the caudate, putamen and thalamus in PD subjects without dementia is also varied. Research groups have demonstrated both the presence (Garg et al., 2015, Mckeown et al., 2008) and absence (Menke et al., 2014, Messina et al., 2011) of morphological changes to the thalamus, while others have shown either the presence (Pitcher et al., 2012, Sterling et al., 2013) or absence (Menke et al., 2014, Owens-Walton et al., 2018) of morphological changes to the caudate or putamen. Studies suggest that PD with dementia is associated with atrophic changes to the caudate nucleus, putamen and thalamus (Summerfield et al., 2005, Owens-Walton et al., 2018), however the presence and extent of such changes in PD-MCI is uncertain (Chen et al., 2016, Melzer et al., 2012). Using a well-defined clinical cohort, and a targeted rs-fMRI FC approach, the current study aims to better understand how FC and morphology of the caudate nucleus, putamen and thalamus are impacted in PD, and how these factors may vary between cognitively unimpaired PD patients, PD-MCI and PD patients with dementia.

A common response to neurological disruption is hyper-connectivity of brain circuits, and it is suggested that such a response reflects compensatory mechanisms in the brain to maintain normal levels of neuronal functioning (Hillary et al., 2015). This mechanism has been used to explain increases in FC of structures within basal ganglia-thalamocortical circuits in PD (Gorges et al., 2015, Helmich et al., 2010) and also in disorders such as mild cognitive impairment and Alzheimer’s disease (Mevel et al., 2011, Sheline and Raichle, 2013). Hillary et al. (2015) argue that hyper-connectivity relies on the availability of sufficient neuronal resources, however these resources can be diminished via the loss of grey matter as diseases processes advance, thus depleting and ultimately exhausting any redundant FC capacity which might otherwise underpin a compensatory response. When this neuronal resource loss reaches a hypothetical ‘critical threshold,’ increased FC is no longer a viable response to situational demand, and subsequent reductions in FC ensue. While the presence and extent of grey matter loss in cognitively unimpaired PD is questionable, PD-MCI and PDD more consistently display grey matter atrophy (Chen et al., 2016, Melzer et al., 2012). Accordingly, we hypothesized that cognitively unimpaired patients with PD would have sufficient neural resources to support functional compensation, and display increased FC, concomitant with no structural (morphological) alterations to the caudate, putamen and thalamus. However, we anticipated the presence of morphological alterations in PD-MCI, representing a ‘critical threshold,’ where FC begins to decrease along with reductions in morphology, evidenced by reduced volumes and surface deflations. Figure 1 is adapted from Hillary et al. (2015) and outlines these hypothesised structural and functional changes in a theoretical framework which describes the onset of FC changes in PD, and a critical threshold where both FC and morphology changes in PD-MCI and PDD can be found.

Figure 1:
Figure 1: Hypothetical ‘critical threshold’ framework.

This figure demonstrates onset of a ‘critical threshold,’ indicated by red crosses, after which point functional connectivity and morphology (volume and surface shape) changes can be found (Hillary et al., 2015).

2. Materials and Methods

2.1 Participants

Participants in this study were from the Swedish BioFINDER study (www.biofinder.se) and gave informed written consent. All elements of the research were performed in accordance with the World Medical Association’s Declaration of Helsinki and ethical approval was obtained through the Ethical Review Board of Lund, Sweden, and the Human Research Ethics Committee at the Australian National University, Canberra, Australia.

2.2 PD disease groups: cognitive impairment

Diagnosis of PD was made by a neurologist using the National Institute of Neurological and Stroke Diagnostic Criteria (Gelb et al., 1999). The total PD cohort was subdivided into three cognitive impairment disease groups. The first group comprises cognitively unimpaired PD patients (PD-CU; n = 33). The second group comprises PD patients with mild cognitive impairment (PD-MCI; n = 22), based on the Movement Disorder Society Task force guidelines (Litvan et al., 2012). This criterion is met when patients score at least 1 standard deviation below the normative mean in at least two cognitive tests, including executive function, attention, visuospatial, memory and language domains. The third subgroup consists of PD participants with dementia (PDD; n = 17) diagnosed using the Clinical Diagnostic Criteria for Dementia Associated with PD (Emre et al., 2007). A healthy control group (Controls) was included for comparative purposes (n = 26). Exclusion criteria for this study included refusing MRI, the presence of significant alcohol or substance misuse and/or a significant unstably systematic illness or organ failure. All participants underwent a cognitive and neurological examination by a medical doctor with extensive experience with movement disorders, while PD patients remained on their usual medication regimes. Clinical functioning of participants was quantified using the Hoehn and Yahr staging, assessing the progression of the disorder, the ‘Unified Parkinson’s Disease Rating Scale Part-III’ test (UPDRS-III), assessing the motor signs of PD (Fahn and Elton, 1987), and the ‘Mini Mental State Examination (MMSE), assessing cognitive mental state (Folstein et al., 1975).

2.3 MRI acquisition

Structural MRI data was obtained on a 3T scanner (Trio, Siemens Magnetom, Erlangen, Germany) with a 20-channel head-coil. High resolution T1-weighted three-dimensional anatomical brain images were acquired using a magnetization-prepared rapid acquisition technique with gradient-echo sequence (repetition time = 7 ms; echo time = 3 ms; flip angle = 90 degrees; voxel size = isotropic 1mm3). Image matrix size was 356 voxels in the coronal and sagittal planes and 176 voxels in the axial plane. Functional MRI data consists of 256 T2*-weighted echo planar imaging volumes (repetition time = 1850 ms; echo time 30 ms; flip angle = 90 degrees; voxel size 3 × 3 × 3.75 mm3) for each participant. Matrix size was 64 voxels in the coronal and sagittal planes and 36 voxels in the axial plane. Subjects were instructed to lie still with their eyes closed, not to fall asleep and not to think of anything in particular during the scan, which lasted for approximately 8 minutes.

2.4 Preprocessing of T1 anatomical MRI data

Preprocessing of T1-weighted structural data was performed using the FSL software package (FSL, Oxford, UK; version 5.0.10) (Jenkinson et al., 2012). Images underwent brain extraction using the FSL-BET tool (Smith, 2002) which were then processed using FSL- FAST, segmenting the T1-image into brain tissue types whilst correcting for image intensity fluctuations caused by inhomogeneities in the radio frequency field (Zhang et al., 2001).

2.5 Preprocessing of resting-state fMRI data

Preprocessing of resting-state functional images was performed using FSL-FEAT (v6.0). Steps included deletion of first 6 volumes, motion correction using FSL-MCFLIRT (Jenkinson et al., 2002) and removal of non-brain tissue using FSL-BET (Smith, 2002). Registration of functional images to subject-specific high resolution T1 weighted structural images using a boundary-based registration (Greve and Fischl, 2009) within FSL-FLIRT linear registration tool (Jenkinson et al., 2002, Jenkinson and Smith, 2001). Registration of functional images to MNI 152 T1 2mm3 standard space was performed with FSL-FLIRT with 12 degrees of freedom, further refined using FSL-FNIRT nonlinear registration (Andersson et al., 2007a, Andersson et al., 2007b) with a warp resolution of 10mm and a resampling resolution of 4mm. Data denoising was then performed with FSL-MELODIC (v3.15) independent component analysis (Beckmann and Smith, 2004) to identify any potential influence of head motion, scanner and cerebrospinal fluid artefacts. FSL-FIX (v1.06) (Salimi- Khorshidi et al., 2014) was then used to remove components classified as ‘noise’ and also correct for motion confounds. Spatial smoothing was then performed (FWHM Gaussian kernel = 5mm) followed by a multiplicative mean intensity normalization and a high-pass temporal filtering Gaussian-weighted least-squares straight line fitting with sigma = 50.0s (Woolrich et al., 2001). For a full explanation of the denoising procedure see ‘1. Supplementary Materials and Methods: 1.1 Independent component analysis denoising.’

2.6 Resting-state fMRI: Seed-based region-of-interest approach

To investigate how the caudate, putamen and thalamus are impacted in PD groups with varying levels of cognitive impairment, we used an atlas-based method to subdivide the structures into functionally relevant regions-of-interest (ROI) (Figure 2A). We performed this step using the Oxford Thalamic Connectivity Atlas (Behrens et al., 2003) and the Oxford- GSK-Imanova Connectivity Atlas (Tziortzi et al., 2014) within FSLEyes. Regions of the dorsal and central caudate nucleus (Figure 1) and the anterior putamen have the highest probability of connectivity with the anterior and dorsolateral prefrontal cortices in the brain (BA9, 9/46 and 10) and are said to be involved in cognitive processes, including perception, memory, reasoning and judgement (Haber, 2003). The mediodorsal and anterior thalamic nuclei have the highest probability of connectivity to the dorsolateral and prefrontal cortices (Behrens et al., 2003) and these thalamic nuclear regions are said to be involved in limbic, episodic memory and cognitive functions (Power and Looi, 2015). For ease of reference, we will refer to these three seed-ROIs as the ‘dorsal caudate’, ‘anterior putamen’ and ‘mediodorsal thalamus’ respectively. These three bilateral standard space MNI 152 2mm3 seed-ROI masks were thresholded to only include voxels that had a greater than 50% likelihood of inclusion within the mask (Figure 2A). Masks were then registered to subject- specific structural T1 space where we used them to ‘cookie-cut’ ROIs from the FSL-FIRST segmentation masks. This was done by overlaying the seed mask over the FSL-FIRST segmentation mask, with the seed-ROI being the neuroanatomical area covered by both masks. The created seed-ROIs were then converted to subject-specific functional MRI space where we performed a final erosion to avoid partial volume effects (Caballero-Gaudes and Reynolds, 2017). Mean activation data from within the three seed-ROIs at each timepoint was used as explanatory variables for subsequent statistical analyses. Before investigating the differences in FC between groups we conducted three one-sample t-tests on all of the Control subjects to look at the FC of our seed-ROIs with the rest of the brain. This step indicates the areas of highest FC with our seed-ROIs and provides a backdrop on which to interpret our results (Figure 2B).

Figure 2:
Figure 2: Seed regions-of-interest and one-sample t-tests demonstrating putative functional connectivity pathways in Control subjects.

A) Subject specific seed regions-of-interest masks were converted to MNI 152 T1 standard space and combined to produce this image. Lighter areas of each colour are indicative of a greater proportion of voxels in that area being included in the relevant mask. B) Z-statistic (Gaussianised) images from Control participant’s individual-level GLM statistical maps, thresholded non-parametrically using clusters determined by Z > 3.1 and a (corrected) cluster significance of p < 0.05 (Worsley, 2001).

2.7 Resting-state fMRI: Statistical analyses

Three individual-level FC analyses were performed for each participant using a general linear model (GLM) in a mass univariate voxelwise whole-brain analysis (Woolrich et al., 2004). This GLM approach measures the correlation between activity within each of the three seed- ROIs and the rest of the brain. We used a FILM pre-whitening correction (Woolrich et al., 2001), shifted the model with a temporal derivative to account for slice timing effects and variations in the hemodynamic response function, removed volumes potentially corrupted by large motion and regressed out average signal from within white matter, ventricle and whole- brain masks. For a full explanation these steps, please see ‘1. Supplementary Materials and Methods: 1.2 Individual-level GLM functional connectivity analyses.’

The following higher-level analyses of FC differences were conducted: PD-CU compared to Controls (section 3.2.1); PD-MCI compared to PD-CU (section 3.2.2); PDD compared to PD- CU (section 3.2.3) and PDD compared to PD-MCI (section 3.3.4). These analyses were performed using FSL-FLAME in a mixed-effects design (Woolrich et al., 2004). Age and gender data were included as covariates. When investigating FC between the PD-CU, PD- MCI and PDD groups we also included levodopa-equivalent daily dosage (LEDD) data as a covariate. Z-statistic (Gaussianized) images from participants individual-level GLM statistical maps were thresholded non-parametrically using clusters determined by Z > 3.1 and a (corrected) cluster significance of p < 0.05 (Worsley, 2001).

2.8 Volumetrics: Statistical analyses

Volumes of the caudate, putamen and thalamus were estimated using the segmentation and registration tool FSL-FIRST (Patenaude et al., 2011). This uses a model-based approach to segment subcortical structures producing output meshes and volumetric data. Investigating differences between experimental groups was performed with SPSS 22.0 (IBM Corporation, Somers, New York, USA) utilising two-tailed multivariate analysis of covariance models controlling for intracranial volume (derived from FSL-FAST) and age. Multivariate partial eta squared values (η2) are used to report effect sizes (Cohen, 1992).

2.9 SPHARM-PDM surface shape: Statistical analyses

Volumetric analyses can investigate overall atrophic differences between groups, however they cannot highlight where this atrophy takes place. Accordingly, we used spherical harmonic parametrization three-dimensional point distribution model (SPHARM-PDM) to investigate localized surface changes to the segmentations of the caudate, putamen and thalamus between experimental groups (Styner et al., 2006). Segmentations were preprocessed to ensure a spherical topology, then described by spherical harmonic functions and sampled onto surfaces of 1002 points. Surfaces were then aligned using a rigid-body Procrustes alignment to a mean template created from a group sample. Comparisons between groups were performed using multivariate analysis of covariance models with a Hotelling statistic in R (version 3.2.1, R Development Core Team, 2014). All analyses included covariates of total intracranial volume, age and sex. SPHARM-PDM produces mean displacement maps which show the magnitude of surface displacement (deflation or inflation in mm) between corresponding points on the mean surface of one experimental group compared to another. Significant surface changes are indicated by regions with a p value < 0.05, corrected for multiple comparisons with a false-discovery rate (FDR) bound q of 5% (Genovese et al., 2002). Results are presented as a composite map of each structure with displacement only displayed in regions that are FDR significant.

3. Results

3.1 Participant characteristics

The characteristics of participants are displayed in Table 1. A chi-square test for independence found no significant difference in sex between the groups [?2 (3, n = 98) = 0.89, p = 0.86]. One-way analyses of variance found no significant difference in age [F (3, 94) = 2.169, p = 0.10], ICV [F (3, 94) = 0.323, p = 0.81], years of education [F (3, 79) = 1.255, p = 0.30] or average DVARS [F (3, 94) = 2.460, p = 0.07] between groups (Table 1).

View this table:
Table 1: Participant characteristics

3.2 Resting-state seed-based functional connectivity analyses

3.2.1 PD-CU compared to Controls

Dorsal caudate

Analysis of the dorsal caudate in PD-CU compared to Controls found clusters of increased FC with the right superior frontal gyrus (BA8), left frontal pole (BA10) and right middle frontal gyrus (BA8) (Figure 3A, Table 2). These clusters extended to include local maxima at the left postcentral gyrus (BA3), left precentral gyrus (BA6), left superior frontal gyrus (BA6), left paracingulate gyrus (BA32), right frontal pole (BA10), left anterior cingulate (BA24) and the left middle frontal gyrus (BA8) (Figure 2A; Supplementary Table 1). We found clusters of decreased FC at the right lateral occipital cortex (BA19), left cerebellar crus I, left angular gyrus (BA39) and left cerebellar lobule VI (Figure 3A; Table 2). These clusters extended to include local maxima at the left supramarginal gyrus (BA39) and left cerebellar vermis VI (Figure 3A; Supplementary Table 1).

View this table:
Table 2: Functional connectivity differences in PD-CU compared to Controls
Figure 3:
Figure 3: Significant between-group differences in functional connectivity of the dorsal caudate (A), anterior putamen (B) and mediodorsal thalamus (C) in PD-CU compared to Controls.

Z- score statistic (Gaussianized) images thresholded non-parametrically using clusters determined by Z > 3.1 and a corrected cluster significance of p < 0.05, overlaid on MNI 152 T1 1mm standard axial brain image (Z = -26mm ascending to 64mm; 6mm interslice distance). Warm colours represent increased in functional connectivity in PD-CU, while shades of blue represent decreased functional connectivity in PD-CU. Abbreviations: Crus I, Crus I of cerebellum; Lobe VI, posterior cerebellar lobule VI; AG, angular gyrus; LOC, lateral occipital gyrus; A, anterior; R, right hemisphere; MFG, middle frontal gyrus; FP, frontal pole; SFG, superior frontal gyrus; PCG, paracingulate gyrus; FPM, frontal pole (medial); PostCG, postcentral gyrus.

Anterior putamen

Analysis of the anterior putamen in PD-CU compared to Controls found clusters of increased FC at the bilateral middle frontal gyri (BA8, 9), left frontal pole (BA10), bilateral postcentral gyri (BA2, 3), left paracingulate gyrus (BA32) and right inferior frontal gyrus (BA44) (Figure 3B; Table 2). These clusters extended to include local maxima at the bilateral anterior cingulate cortex (BA24), right frontal pole (BA10), left superior frontal gyrus (BA9), bilateral precentral gyri (BA4, 6) and bilateral central operculum (BA42) (Figure 3B; Supplementary Table 2). We also found clusters of decreased FC at the left lateral occipital cortex (BA7) and right precuneus (BA7) (Figure 3B; Table 2). These clusters extended to include a local maxima at the left precuneus (Figure 3B; Supplementary Table 2).

Mediodorsal thalamus

Analysis of the mediodorsal thalamus in PD-CU compared to Controls found clusters of increased FC at the left postcentral gyrus (BA1) and the right middle frontal gyrus (BA8) (Figure 3C; Table 2). These clusters extended to include local maxima at the left postcentral gyrus (BA3) and the left precentral gyrus (BA4, 6) (Figure 3C; Supplementary Table 3). No significant clusters of decreased FC were found for the mediodorsal thalamus.

View this table:
Table 3: Functional connectivity differences in PD-MCI compared to PD-CU

3.2.2 PD-MCI compared to PD-CU

Mediodorsal thalamus

Analysis of the mediodorsal thalamus in PD-MCI compared to PD- CU found one significant cluster of decreased FC at the right paracingulate gyrus (BA32) (Figure 4; Table 3). This cluster extended to include local maxima at the medial frontal gyrus (BA9) and right frontal pole (BA8) (Figure 4; Supplementary Table 4). No significant clusters of increased FC were found. There were no significant FC differences to the dorsal caudate or anterior putamen when comparing PD-MCI and PD-CU.

View this table:
Table 4: Functional connectivity differences in PDD compared to PD-MCI
Figure 4:
Figure 4: Significant between-group differences in functional connectivity of the mediodorsal thalamus in PD-MCI compared to PD-CU.

Z-score statistic (Gaussianized) images thresholded non-parametrically using clusters determined by Z > 3.1 and a corrected cluster significance of p < 0.05, overlaid on MNI 152 T1 1mm axial standard brain image (Z = -30.5mm ascending to 45.5mm; 5mm interslice distance). Shades of blue show areas of decreased in functional connectivity in PD- MCI. Abbreviations: A, anterior; R, right hemisphere; PCG, paracingulate gyrus.

3.2.3. PDD compared to PD-CU

There were no significant FC differences to the dorsal caudate, anterior putamen or mediodorsal thalamus when comparing PD-CU subjects to PDD subjects.

3.2.4 PDD compared to PD-MCI

Mediodorsal thalamus

Analysis of the mediodorsal thalamus in PDD compared to PD-MCI found one cluster of decreased FC at the right posterior cingulate cortex (BA23) (Figure 5; Table 4; Supplementary Table 5). There were no significant FC differences of the dorsal caudate or anterior putamen when comparing PDD with PD-MCI.

Figure 5:
Figure 5: Significant between-group differences in functional connectivity of the mediodorsal thalamus in PDD compared to PD-MCI.

Z-score statistic (Gaussianized) images thresholded non- parametrically using clusters determined by Z > 3.1 and a corrected cluster significance of p < 0.05, overlaid on MNI 152 T1 1mm axial standard brain image (Z = -3mm ascending to 57mm; 4mm interslice distance). Shades of blue show areas of decreased in functional connectivity in PDD. Abbreviations: PCC, Posterior cingulate cortex; R, right hemisphere; A, anterior.

3.3 Volumetric analyses

Pairwise comparisons of volumes (Supplementary Table 6 and 7) revealed the following significant differences between experimental groups: the left caudate was significantly reduced in PDD compared to PD-CU (p = 0.009) and Controls (p < 0.001) with a large effect size (η2 = 0.177) (Figure 6A). The right caudate was significantly reduced in PDD compared to PD-CU (p = 0.021) and Controls (p = 0.001) also with a large effect size (η2 = 0.143) (Figure 6A). The left putamen was significantly reduced in PDD compared to PD-CU (p = 0.005) and Controls (p = 0.05) with a large effect size (η2 = 0.193) (Figure 6B). The right putamen was significantly reduced in PDD compared to Controls (p = 0.016) with a medium effect size (η2 = 0.107) (Figure 6B). The left thalamus was significantly reduced in PDD compared to PD-CU (p = 0.011) and Controls (p = 0.037) with a medium effect size (η2 = 0.108) (Figure 6C). The right thalamus was significantly reduced in PDD compared to PD- CU (p = 0.002) with a medium effect size (η2 = 0.136) (Figure 6C).

Figure 6:
Figure 6: Pairwise comparison of volumes between PD subgroups and Controls.

(A) Caudate (B) putamen and (C) thalamus. Volumes presented are estimated marginal mean values for each experimental group derived from MANCOVA models controlling for age and ICV. *, indicates a pairwise comparison with a significant p-value <0.05 confidence interval correction: Bonferroni.

3.4 SPHARM-PDM shape analyses

SPHARM-PDM shape analysis of the caudate, putamen and thalamus found widespread deflation to the surface of these structures in PDD compared to PD-CU subjects. Surface deflation can be observed at the bilateral rostral as well as dorsomedial surfaces areas of the caudate (Figure 7). Widespread deflation was also found on the surface of the left putamen, at the anterior, dorsal and posterior regions. Significant and widespread deflation was found across the medial surface of the bilateral thalami, potentially corresponding to surface regions of the anterior and mediodorsal thalamic nuclei. Surface deflation was also found across the posterior-ventral surfaces of the thalami, potentially corresponding to the pulvinar as well as medial and lateral geniculate bodies. No FDR-corrected significant shape changes were found when comparing PD-CU to PD-MCI or PD-MCI to PDD.

Figure 7:
Figure 7: SPHARM-PDM shape analysis of the caudate nucleus, putamen and thalamus in PDD compared to PD-CU subjects.

Neuroanatomical surface areas of FDR-corrected mean difference displacement maps. Displacement colour scale corresponds to the millimetres of deflation/inflation of the surface in that region in mm; warmer colours correspond to surface expansion in PDD while cooler colours correspond to surface deflation in PDD. The bilateral thalami sit medially, the caudate nuclei rostrally and the putamen are the most lateral of the three structures. Abbreviations: R, right hemisphere; FDR, false-discovery rate; Left, looking at the structures from the left side of the brain; R, looking at the structures from the right side of the brain.

4. Discussion

4.1 Study overview

This study investigates the pattern of caudate nucleus, putamen and thalamus functional connectivity (FC) and morphology alterations associated with cognitive impairment and dementia in PD. Our study contributes new knowledge by showing how important functional subdivisions of these structures are impacted across cognitive disease stages. Our results indicate that PD-CU subjects display increased FC of the dorsal caudate, anterior putamen and mediodorsal thalamus with neuroanatomical areas across the frontal lobe, as well as reduced FC of the dorsal caudate with posterior cortical and cerebellar regions, compared to Controls. PD-MCI subjects demonstrated reduced FC of the mediodorsal thalamus with midline structures within the executive control network while PDD subjects demonstrated reduced FC of the mediodorsal thalamus with the posterior cingulate cortex (PCC), a key node within the default mode network (DMN). Significant volumetric and surface-based changes to the caudate nucleus, putamen and thalamus were found exclusively in our PDD cohort, suggesting that FC changes in PD-CU and PD-MCI take place independently of observable morphological alterations.

4.2 Functional connectivity analyses: PD-CU compared to Controls

Our findings of increased FC between the dorsal caudate, anterior putamen and mediodorsal thalamus with the pre- and postcentral gyri support a recent meta-analysis indicating that PD is associated with increased FC of these brain regions (Ji et al., 2018). After weighing the relative contributions of seed-ROIs, Ji et al. (2018) were able to demonstrate that seeds placed within the caudate, putamen and thalamus contributed most to this result. More specifically, our finding supports work by Gorges et al. (2015) who found increases in FC between a basal ganglia-thalamic network and areas of the prefrontal cortex and anterior/paracingulate gyri. These researchers argue that such hyper-connectivity in cognitively unimpaired PD subjects may be representative of the recruitment of additional resources in the brain to maintain cognitive performance (Gorges et al., 2015). Our results support our hypotheses and may signal how resources can be recruited to maintain ‘normal’ cognitive functioning. These increases in FC between the dorsal caudate and regions across the prefrontal cortex and anterior/paracingulate gyri may be indicative of compensatory changes taking place within the executive control network (ECN). This network is comprised of dorsal and ventrolateral prefrontal, frontoinsular, lateral parietal, middle prefrontal cortices, anterior and paracingulate gyri, as well as subcortical sites at the mediodorsal and anterior thalamic nuclei and the dorsal caudate (Beckmann et al., 2005, Seeley et al., 2007). The ECN has been validated on a large scale and has shown that activity within this network corresponds strongly with several paradigms of cognition in task-based FC analyses (Smith et al., 2009).

Our data also indicate that PD-CU is associated with a reduction in FC between the dorsal caudate and the cerebellum. This result supports the findings of Hacker et al. (2012) who demonstrated decreases in FC between the caudate nucleus and an ‘extended brainstem’ region encompassing the cerebellum. These findings are significant because loss of FC between the caudate and the cerebellum may serve as a viable candidate for PD related gait difficulties, postural instability and freezing, which can be linked to dysfunction in cerebellar circuitry (Hacker et al., 2012). Our data also indicate that PD-CU subjects have decreased FC between the dorsal caudate/anterior putamen and the precuneus and angular gyrus. This finding supports the work of Amboni et al. (2015) who found decreases in FC between nodes in the DMN, including lateral and medial parietal regions.

While our results support a number of previous research findings, they contrast with similar work, necessitating explication. Bell et al. (2015) demonstrated no significant differences in FC of the anterior striatum in a PD cohort (when placed on medication) compared to Control subjects. They also found no differences in FC of seed-ROIs in the thalamus with intrinsic connectivity networks. Similarly, our data contrast with the work of Agosta et al. (2014) who found decreased FC between the caudate and putamen and brain regions across the frontal lobe, as well as hypo-connectivity of the thalamus with insular cortices. There are a number of methodological considerations that may explain these discrepancies. The work of Bell and colleagues (2015) involved a younger PD cohort to the current study, which brings into question the impact of age related changes in FC (Fair et al., 2008), while the work of Agosta et al. (2014) focussed on a PD cohort with unilateral disease presentation, studying FC of the ‘more affected’ and ‘less affected’ hemispheres. This approach increases the specificity of their analysis, however it makes comparisons with the present data difficult as we averaged FC signals from seeds within both hemispheres at each timepoint.

4.3 Functional connectivity analyses: PD-MCI compared to PD-CU

PD-MCI subjects showed decreased FC of the mediodorsal thalamus with the paracingulate gyrus, medial frontal gyrus and frontal pole compared to cognitively unimpaired PD subjects. Our findings of decreased FC of the mediodorsal thalamus support the work of Gorges et al. (2015) who showed that PD-MCI is associated with lower FC between a basal ganglia-thalamic intrinsic connectivity network and the anterior/paracingulate gyri compared to a cognitively unimpaired PD cohort. Decreased FC between the mediodorsal thalamus and medial prefrontal brain regions provides further evidence for the important link between midline nodes within the ECN and cognitive function. As outlined, structures within this intrinsic connectivity network are related to cognitive function (Beckmann et al., 2005, Seeley et al., 2007), as they work together to provide bias signals to other areas of the brain to implement and maintain cognitive control processes (Miller and Cohen, 2001). While we have demonstrated increased FC of nodes within this network in PD-CU subjects, here we demonstrate that PD-MCI subjects have a reduced FC between the mediodorsal thalamus and paracingulate cortices. While traditionally thought to be subserve affective functions, the cingulate has emerged as an important structure associated with higher cognitive processes, and has been shown to be particularly vulnerable to disease processes where executive control is impaired (Carter et al., 2000). Seeley et al. (2007) demonstrated that stronger FC between nodes of the ECN correlated positively with superior executive task performance of subjects. If the corollary holds true, poorer cognitive performance, as our PD-MCI subjects demonstrate, would be associated with reduced FC within this network, which we have demonstrated.

4.4 Functional connectivity analyses: PDD compared to PD-MCI

PDD subjects showed decreased FC of the mediodorsal thalamus with the posterior PCC compared to PD-MCI subjects, directly supporting recent graph theoretical research (Zhan et al., 2018). Decreases in FC between these regions is significant due to the crucial role played by the PCC in cognitive functioning as well as the role of the mediodorsal thalamus and PCC within the DMN. The functional role played by the PCC involves elements of memory consolidation, attention and the control required to balance internally and externally focused thoughts (Leech and Sharp, 2014). The PCC is affected in neurodegenerative diseases like Alzheimer’s disease (Buckner et al., 2009, Greicius et al., 2004, Mevel et al., 2011) and atrophy of the PCC has been demonstrated in PD subjects with dementia (Melzer et al., 2012). The PCC has important structural and also functional connectivity with the thalamus (Cunningham et al., 2017) and Greicius et al. (2003) suggested that the structure may play an important role balancing the interaction between brain regions in the prefrontal cortex that are involved cognitive processes and posterior nodes within the DMN. These researchers suggested that the thalamus is uniquely positioned to act as an intermediary structure in the brain that balances the interaction between these regions due to its extensive cortico-cortical connectivity (Greicius et al., 2003). Our research supports this idea, highlighting how the FC between the mediodorsal thalamus and the PCC is impacted in PD subjects with dementia.

4.5 Morphological analyses

Our research found no evidence of morphological (structural) alterations to the caudate, putamen or thalamus in PD-CU subjects compared to Controls, supporting the findings from a number of studies (Lee et al., 2014, Mak et al., 2014, Mckeown et al., 2008, Menke et al., 2014, Messina et al., 2011, Nemmi et al., 2015, Tinaz et al., 2011). We also found no evidence of atrophy to the caudate, putamen or thalamus in PD-MCI compared to Controls, supporting work by Amboni et al. (2015). Our data do indicate that PDD subjects have reduced caudate volumes (bilaterally) in PDD compared to PD-CU and Controls. Left putamen and left thalamus volumes were reduced in PDD compared to PD-CU and Controls, while the right putamen and right thalamus was reduced in PDD compared to Controls and PD-CU, respectively. These results support studies showing reduced caudate, putamen and thalamic grey matter changes in PD subjects with dementia compared to controls (Burton et al., 2004, Melzer et al., 2012, Summerfield et al., 2005). Shape analysis of the caudate, putamen and thalamus found widespread surface deflation in PDD compared to PD-CU subjects, fine-tuning these volumetric results. Regarding the caudate, surface deflation was found primarily at the rostral and dorsomedial surfaces, areas with putative connectivity to the executive regions of the frontal lobe as well as the rostral motor regions (Tziortzi et al., 2014). Surface deflation was also found at the anterior, dorsal and posterior putamen surface regions areas with putative connectivity to the executive regions of the frontal lobe as well as caudal motor regions (Tziortzi et al., 2014). Widespread surface deflation was found primarily across the medial surface of the bilateral thalami, corresponding to surface regions of the anterior and mediodorsal thalamic nuclei with connections to the dorsolateral and prefrontal cortices (Behrens et al., 2003) as well as posterior cingulate, precuneus and hippocampal brain regions (Cunningham et al., 2017). No significant surface changes were found when comparing the other experimental groups. Our findings support the work of Mak et al. (2014) who found no significant difference between PD-MCI and cognitively unimpaired PD subjects using a vertex-wise shape analysis method. Our morphological data on PD-CU subjects are also consistent with previous studies which indicate that PD is not associated with localized surface alterations (Lee et al., 2014, Menke et al., 2014, Messina et al., 2011) however they contrast with the results of two groups who found localised surface deflation in PD (Garg et al., 2015, Mckeown et al., 2008). As raised previously, methodological inconsistencies between studies may be a significant factor behind these differences.

4.6 A hypothetical model of functional connectivity and morphological alterations associated with cognitive impairment in PD

Our findings provide only partial support for our model regarding the relationship between FC and morphology across PD cognitive disease states (Figure 1). PD-CU patients showed no observable morphological changes compared to Controls, however they displayed significant increases in FC with prefrontal regions of the brain. This may be evidence of a compensatory mechanism whereby the brain recruits additional neuronal resources to maintain normal cognitive function, due to impact of disease-related pathology. We hypothesised that these recruited neuronal resources would be diminished in PD-MCI and PDD, evidenced by morphological alterations to the caudate nucleus, putamen and thalamus, and would accompany reductions in FC. Whilst we found that PD-MCI subjects showed reductions in FC compared to PD-CU subjects, this was despite any observable morphological differences between the groups. PDD subjects displayed significant reductions in surface morphology as well as reductions in seed-based FC, supporting the proposed model. However, our research was unable to find any significant differences in FC between PD-CU and PDD subjects, which we would expect considering a difference in FC was found between the PDD and PD-MCI groups.

4.7 Limitations and future research

This research was performed in a cross-sectional manner, logically assuming a sequential temporal relationship between PD-CU, PD-MCI and PDD groupings, though we acknowledge this is not a longitudinal analysis. Categorisation of PD subjects in this way is a potential limitation, as we cannot be sure that all participants with PD-CU will progress to PD-MCI, then to PDD. While there is strong evidence that PD-MCI does represent an intermediate cognitive disease stage between cognitively unimpaired PD and PDD (Litvan et al., 2012), future research should follow a longitudinal design to better study how functional connectivity and morphology of basal ganglia-thalamocortical circuits are related to cognitive impairment in PD. This would have the added benefit of reducing the influence of the significant disease heterogeneity in PD which has been shown to have a significant impact on the reproducibility of results in this field of research (Badea et al., 2017).

4.8 Conclusion

This study demonstrated how the functional connectivity of subcortical hubs within basal ganglia-thalamocortical circuits are implicated in PD between cognitive disease stages. Functional connectivity changes are found in cognitively unimpaired PD subjects, who display no morphological changes, while PD-MCI subjects and PDD subjects both display decreases in functional connectivity of the mediodorsal thalamus with the anterior and posterior cingulate cortices respectively. Our research highlights the importance of the caudate, putamen and thalamus as core structures within key intrinsic connectivity networks in the brain, and indicates that a breakdown in functional connectivity of the mediodorsal thalamus is tied to mild cognitive impairment and dementia in PD.

Data Availability

Anonymized data will be shared by request from any qualified investigator for the sole purpose of replicating procedures and results presented in the article and as long as data transfer is in agreement with EU legislation on the general data protection regulation. Data underlying the results described in our manuscript are available from at Attention: Conor Owens-Walton Academic Unit of Psychiatry & Addiction Medicine Australian National University Medical School Building 4, Level 2, Canberra Hospital Woden, A.C.T. 2605 AUSTRALIA or (conor.owens-walton@anu.edu.au)

Funding

CO-W would like to acknowledge the Australian National University for their funding support via the University Research Scholarship. Work at the authors’ research center was supported by the European Research Council, the Swedish Research Council, the Knut and Alice Wallenberg foundation, the Marianne and Marcus Wallenberg foundation, the Strategic Research Area MultiPark (Multidisciplinary Research in Parkinson’s disease) at Lund University, the Swedish Alzheimer Foundation, the Swedish Brain Foundation, The Parkinson foundation of Sweden, The Parkinson Research Foundation, the Skåne University Hospital Foundation, and the Swedish federal government under the ALF agreement. DJ reports no funding sources. MW has received consulting fees and honoraria from Actelion Pharmaceuticals and receives royalties from ACER (Australian Council of Educational research) for the NUCOG, a pencil and paper cognitive assessment tool (not used in this study). Funding sources had no role in the design and conduct of the study, in the collection, analysis, interpretation of the data or in the preparation, review or approval of the manuscript.

Competing interests

OH has acquired research support (for the institution) from Roche, GE Healthcare, Biogen, AVID Radiopharmaceuticals and Euroimmun. In the past 2 years, he has received consultancy/speaker fees (paid to the institution) from Biogen and Roche.

Acknowledgements

The authors are indebted to all patients and control subjects who participated in this study. This project is an initiative of the Swedish BioFINDER Study, of whom DvW and OH are steering committee members, and also the AUSSIE network coordinated by JCLL at the Australian National University Medical School, who self-funds related expenses.

Footnotes

  • ↵ Denotes co-senior authors

  • Declaration of interest: Declaration of interest: None

References

  1. ↵
    Agosta, F., Caso, F., Stankovic, I., Inuggi, A., Petrovic, I., Svetel, M., Kostic, V. S. & Filippi,M. 2014. Cortico-striatal-thalamic network functional connectivity in hemiparkinsonism. Neurobiology of Aging, 35, 2592–2602.
    OpenUrlPubMed
  2. ↵
    Amboni, M., Tessitore, A., Esposito, F., Santangelo, G., Picillo, M., Vitale, C., Giordano, A., Erro, R., De Micco, R. & Corbo, D. 2015. Resting-state functional connectivity associated with mild cognitive impairment in Parkinson’s disease. Journal of neurology, 262, 425–434.
    OpenUrlCrossRefPubMed
  3. ↵
    Andersson, J. L. R., Jenkinson, M. & Smith, S. M. 2007a. Non-linear optimisation. FMRIB technical report TR07JA1.
  4. ↵
    Andersson, J. L. R., Jenkinson, M. & Smith, S. M. 2007b. Non-linear optimisation, aka spatial normalisation. FMRIB technical report TR07JA2.
  5. ↵
    Badea, L., Onu, M., Wu, T., Roceanu, A. & Bajenaru, O. 2017. Exploring the reproducibility of functional connectivity alterations in Parkinson’s disease. PLoS ONE, 12, e0188196.
    OpenUrl
  6. ↵
    Beckmann, C. F., Deluca, M., Devlin, J. T. & Smith, S. M. 2005. Investigations into restingstate connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci, 360, 1001–13.
    OpenUrlCrossRefPubMed
  7. ↵
    Beckmann, C. F. & Smith, S. M. 2004. Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging, 23, 137–52.
    OpenUrlCrossRefPubMedWeb of Science
  8. ↵
    Behrens, T. E., Johansen-Berg, H., Woolrich, M. W., Smith, S. M., Wheeler-Kingshott, C. A., Boulby, P. A., Barker, G. J., Sillery, E. L., Sheehan, K., Ciccarelli, O., Thompson, A. J., Brady, J. M. & Matthews, P. M. 2003. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nature Neuroscience, 6, 750–7.
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    Bell, P. T., Gilat, M., O’callaghan, C., Copland, D. A., Frank, M. J., Lewis, S. J. G. & Shine, J. M. 2015. Dopaminergic basis for impairments in functional connectivity across subdivisions of the striatum in Parkinson’s disease. Human Brain Mapping, 36, 1278–1291.
    OpenUrlCrossRefPubMed
  10. ↵
    Broeders, M., De Bie, R., Velseboer, D. C., Speelman, J. D., Muslimovic, D. & Schmand, B. 2013. Evolution of mild cognitive impairment in Parkinson disease. Neurology, 81, 346–352.
    OpenUrlCrossRef
  11. ↵
    Buckner, R. L., Sepulcre, J., Talukdar, T., Krienen, F. M., Liu, H., Hedden, T., Andrews-Hanna, J. R., Sperling, R. A. & Johnson, K. A. 2009. Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer’s disease. Journal of Neuroscience, 29, 1860–1873.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    Burton, E. J., Mckeith, I. G., Burn, D. J., Williams, E. D. & O’brien, J. T. 2004. Cerebral atrophy in Parkinson’s disease with and without dementia: a comparison with Alzheimer’s disease, dementia with Lewy bodies and controls.
  13. ↵
    Caballero-Gaudes, C. & Reynolds, R. C. 2017. Methods for cleaning the BOLD fMRI signal. Neuroimage, 154, 128–149.
    OpenUrl
  14. ↵
    Carter, C. S., Macdonald, A. M., Botvinick, M., Ross, L. L., Stenger, V. A., Noll, D. & Cohen, J. D. 2000. Parsing executive processes: strategic vs. evaluative functions of the anterior cingulate cortex. Proceedings of the National Academy of Sciences of the United States of America, 97, 1944–1948.
    OpenUrlAbstract/FREE Full Text
  15. ↵
    Chen, F. X., Kang, D. Z., Chen, F. Y., Liu, Y., Wu, G., Li, X., Yu, L. H., Lin, Y. X. & Lin, Z. Y. 2016. Gray matter atrophy associated with mild cognitive impairment in Parkinson’s disease. Neuroscience Letters, 617, 160–5.
    OpenUrl
  16. ↵
    Cohen, J. 1992. Statistical power analysis. Current directions in psychological science, 1, 98–101.
    OpenUrlCrossRef
  17. ↵
    Cunningham, S. I., Tomasi, D. & Volkow, N. D. 2017. Structural and functional connectivity of the precuneus and thalamus to the default mode network. Hum Brain Mapp, 38, 938–956.
    OpenUrlCrossRef
  18. ↵
    Damoiseaux, J. S., Rombouts, S. a. R. B., Barkhof, F., Scheltens, P., Stam, C. J., Smith, S. M. & Beckmann, C. F. 2006. Consistent resting-state networks across healthy subjects. Proceedings of the National Academy of Sciences of the United States of America, 103, 13848–13853.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    Emre, M., Aarsland, D., Brown, R., Burn, D. J., Duyckaerts, C., Mizuno, Y., Broe, G. A., Cummings, J., Dickson, D. W., Gauthier, S., Goldman, J., Goetz, C., Korczyn, A., Lees, A., Levy, R., Litvan, I., Mckeith, I., Olanow, W., Poewe, W., Quinn, N., Sampaio, C., Tolosa, E. & Dubois, B. 2007. Clinical diagnostic criteria for dementia associated with Parkinson’s disease. Movment Disorders, 22, 1689-707; quiz 1837.
    OpenUrl
  20. ↵
    Fahn, S. & Elton, R. (eds.) 1987. Members of the UPDRS development committee: The unified Parkison disease rating scale, Florham Park, NJ.: Macmillan Health Care Information.
  21. ↵
    Fair, D. A., Cohen, A. L., Dosenbach, N. U., Church, J. A., Miezin, F. M., Barch, D. M., Raichle, M. E., Petersen, S. E. & Schlaggar, B. L. 2008. The maturing architecture of the brain’s default network. Proceedings of the National Academy of Sciences, 105, 4028–4032.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    Folstein, M. F., Folstein, S. E. & Mchugh, P. R. 1975. “Mini-mental state”: a practical method for grading the cognitive state of patients for the clinician. Journal of psychiatric research, 12, 189–198.
    OpenUrlCrossRefPubMedWeb of Science
  23. ↵
    Garg, A., Appel-Cresswell, S., Popuri, K., Mckeown, M. J. & Beg, M. F. 2015. Morphological alterations in the caudate, putamen, pallidum, and thalamus in Parkinson’s disease. Frontiers in Neuroscience, 9.
  24. ↵
    Gelb, D. J., Oliver, E. & Gilman, S. 1999. Diagnostic criteria for Parkinson disease. Archives of neurology, 56, 33–39.
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    Genovese, C. R., Lazar, N. A. & Nichols, T. 2002. Thresholding of statistical maps in functional neuroimaging using the false discovery rate. Neuroimage, 15, 870–878.
    OpenUrlCrossRefPubMedWeb of Science
  26. ↵
    Goldman, J. G., Holden, S. K., Litvan, I., Mckeith, I., Stebbins, G. T. & Taylor, J. P. 2018. Evolution of diagnostic criteria and assessments for Parkinson’s disease mild cognitive impairment. Movement Disorders, 33, 503–510.
    OpenUrl
  27. ↵
    Gorges, M., Muller, H. P., Lule, D., Pinkhardt, E. H., Ludolph, A. C. & Kassubek, J. 2015. To rise and to fall: functional connectivity in cognitively normal and cognitively impaired patients with Parkinson’s disease. Neurobiology of Aging, 36, 1727–1735.
    OpenUrl
  28. ↵
    Greicius, M. D., Krasnow, B., Reiss, A. L. & Menon, V. 2003. Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc Natl Acad Sci U S A, 100, 253–8.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    Greicius, M. D., Srivastava, G., Reiss, A. L. & Menon, V. 2004. Default-mode network activity distinguishes Alzheimer’s disease from healthy aging: Evidence from functional MRI. Proceedings of the National Academy of Sciences of the United States of America, 101, 4637–4642.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    Greve, D. N. & Fischl, B. 2009. Accurate and robust brain image alignment using boundary-based registration. Neuroimage, 48, 63–72.
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    Haber, S. N. 2003. The primate basal ganglia: parallel and integrative networks. Journal of chemical neuroanatomy, 26, 317–330.
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    Hacker, C. D., Perlmutter, J. S., Criswell, S. R., Ances, B. M. & Snyder, A. Z. 2012. Resting state functional connectivity of the striatum in Parkinson’s disease. Brain, 135, 3699–3711.
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    Halliday, G. M. 2009. Thalamic changes in Parkinson’s disease. Parkinsonism & Related Disorders, 15, Supplement 3, S152–S155.
    OpenUrl
  34. ↵
    Helmich, R. C., Derikx, L. C., Bakker, M., Scheeringa, R., Bloem, B. R. & Toni, I. 2010. Spatial remapping of cortico-striatal connectivity in Parkinson’s disease. Cerebral cortex, 20, 1175–1186.
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    Hillary, F. G., Roman, C. A., Venkatesan, U., Rajtmajer, S. M., Bajo, R. & Castellanos, N. D. 2015. Hyperconnectivity is a fundamental response to neurological disruption. Neuropsychology, 29, 59–75.
    OpenUrlCrossRefPubMed
  36. ↵
    Hwang, K., Bertolero, M. A., Liu, W. B. & D’esposito, M. 2017. The human thalamus is an integrative hub for functional brain networks. Journal of Neuroscience, 37, 5594–5607.
    OpenUrlAbstract/FREE Full Text
  37. ↵
    Jenkinson, M., Bannister, P., Brady, M. & Smith, S. 2002. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage, 17, 825–41.
    OpenUrlCrossRefPubMedWeb of Science
  38. ↵
    Jenkinson, M., Beckmann, C. F., Behrens, T. E. J., Woolrich, M. W. & Smith, S. M. 2012. FSL. Neuroimage, 62, 782–790.
    OpenUrlCrossRefPubMedWeb of Science
  39. ↵
    Jenkinson, M. & Smith, S. 2001. A global optimisation method for robust affine registration of brain images. Medical image analysis, 5, 143–156.
    OpenUrlCrossRefPubMedWeb of Science
  40. ↵
    Ji, G.-J., Hu, P., Liu, T.-T., Li, Y., Chen, X., Zhu, C., Tian, Y., Chen, X. & Wang, K. 2018. Functional Connectivity of the Corticobasal Ganglia–Thalamocortical Network in Parkinson Disease: A Systematic Review and Meta-Analysis with Cross-Validation. Radiology, 287, 172183.
    OpenUrl
  41. ↵
    Khan, A. R., Hiebert, N. M., Vo, A., Wang, B. T., Owen, A. M., Seergobin, K. N. & Macdonald, P. A. 2018. Biomarkers of Parkinson’s disease: Striatal sub-regional structural morphometry and diffusion MRI. Neuroimage: Clinical.
  42. ↵
    Kish, S. J., Shannak, K. & Hornykiewicz, O. 1988. Uneven pattern of dopamine loss in the striatum of patients with idiopathic Parkinson’s disease. Pathophysiologic and clinical implications. New England Journal of Medicine, 318, 876–80.
    OpenUrlCrossRefPubMedWeb of Science
  43. ↵
    Lee, H. M., Kwon, K.-Y., Kim, M.-J., Jang, J.-W., Suh, S.-I., Koh, S.-B. & Kim, J. H. 2014. Subcortical grey matter changes in untreated, early stage Parkinson’s disease without dementia. Parkinsonism & related disorders, 20, 622–626.
    OpenUrl
  44. ↵
    Leech, R. & Sharp, D. J. 2014. The role of the posterior cingulate cortex in cognition and disease. Brain, 137, 12–32.
    OpenUrlCrossRefPubMedWeb of Science
  45. ↵
    Litvan, I., Goldman, J. G., Troster, A. I., Schmand, B. A., Weintraub, D., Petersen, R. C., Mollenhauer, B., Adler, C. H., Marder, K., Williams-Gray, C. H., Aarsland, D., Kulisevsky, J., Rodriguez-Oroz, M. C., Burn, D. J., Barker, R. A. & Emre, M. 2012. Diagnostic criteria for mild cognitive impairment in Parkinson’s disease: Movement Disorder Society Task Force guidelines. Movment Disorders, 27, 349–56.
    OpenUrl
  46. ↵
    Looi, J. C. & Walterfang, M. 2013. Striatal morphology as a biomarker in neurodegenerative disease. Molecular Psychiatry, 18, 417–24.
    OpenUrlCrossRefPubMed
  47. ↵
    Looi, J. C., Walterfang, M., Nilsson, C., Power, B. D., Van Westen, D., Velakoulis, D., Wahlund, L. O. & Thompson, P. M. 2014. The subcortical connectome: Hubs, spokes and the space between - a vision for further research in neurodegenerative disease. Australian and New Zealand Journal of Psychiatry.
  48. ↵
    Mak, E., Bergsland, N., Dwyer, M. G., Zivadinov, R. & Kandiah, N. 2014. Subcortical atrophy is associated with cognitive impairment in mild Parkinson disease: a combined investigation of volumetric changes, cortical thickness, and vertex-based shape analysis. American Journal of Neuroradiology, 35, 2257–64.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    Mckeown, M. J., Uthama, A., Abugharbieh, R., Palmer, S., Lewis, M. & Huang, X. 2008. Shape (but not volume) changes in the thalami in Parkinson disease. BMC neurology, 8, 8.
    OpenUrl
  50. ↵
    Melzer, T. R., Watts, R., Macaskill, M. R., Pitcher, T. L., Livingston, L., Keenan, R. J., Dalrymple-Alford, J. C. & Anderson, T. J. 2012. Grey matter atrophy in cognitively impaired Parkinson’s disease. Journal of Neurology, Neurosurgery and Psychiatry, 83, 188.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    Menke, R. A., Szewczyk□Krolikowski, K., Jbabdi, S., Jenkinson, M., Talbot, K., Mackay, C. E. & Hu, M. 2014. Comprehensive morphometry of subcortical grey matter structures in early□stage Parkinson’s disease. Human brain mapping, 35, 1681–1690.
    OpenUrlCrossRefPubMed
  52. ↵
    Messina, D., Cerasa, A., Condino, F., Arabia, G., Novellino, F., Nicoletti, G., Salsone, M., Morelli, M., Lanza, P. L. & Quattrone, A. 2011. Patterns of brain atrophy in Parkinson’s disease, progressive supranuclear palsy and multiple system atrophy. Parkinsonism & related disorders, 17, 172–176.
    OpenUrl
  53. ↵
    Mevel, K., Chetelat, G., Eustache, F. & Desgranges, B. 2011. The default mode network in healthy aging and Alzheimer’s disease. International journal of Alzheimer’s disease, 2011, 535816.
    OpenUrl
  54. ↵
    Miller, E. K. & Cohen, J. D. 2001. An integrative theory of prefrontal cortex function. Annu Rev Neurosci, 24, 167–202.
    OpenUrlCrossRefPubMedWeb of Science
  55. ↵
    Nemmi, F., Sabatini, U., Rascol, O. & Peran, P. 2015. Parkinson’s disease and local atrophy in subcortical nuclei: insight from shape analysis. Neurobiology of Aging, 36, 424–33.
    OpenUrlCrossRefPubMed
  56. ↵
    Obeso, J. A., Rodriguez-Oroz, M. C., Rodriguez, M., Lanciego, J. L., Artieda, J., Gonzalo, N. & Olanow, C. W. 2000. Pathophysiology of the basal ganglia in Parkinson’s disease. Trends in neurosciences, 23, S8–S19.
    OpenUrlCrossRefPubMedWeb of Science
  57. ↵
    Owens-Walton, C., Jakabek, D., Li, X., Wilkes, F. A., Walterfang, M., Velakoulis, D., Van Westen, D., Looi, J. C. L. & Hansson, O. 2018. Striatal changes in Parkinson disease: An investigation of morphology, functional connectivity and their relationship to clinical symptoms. Psychiatry Research: Neuroimaging, 275, 5–13.
    OpenUrl
  58. ↵
    Patenaude, B., Smith, S. M., Kennedy, D. N. & Jenkinson, M. 2011. A Bayesian model of shape and appearance for subcortical brain segmentation. Neuroimage, 56, 907–922.
    OpenUrlCrossRefPubMedWeb of Science
  59. ↵
    Pitcher, T. L., Melzer, T. R., Macaskill, M. R., Graham, C. F., Livingston, L., Keenan, R. J., Watts, R., Dalrymple-Alford, J. C. & Anderson, T. J. 2012. Reduced striatal volumes in Parkinson’s disease: a magnetic resonance imaging study. Translational Neurodegeneration, 1, 68–76.
    OpenUrl
  60. ↵
    Poewe, W., Seppi, K., Tanner, C. M., Halliday, G. M., Brundin, P., Volkmann, J., Schrag, A.-E. & Lang, A. E. 2017. Parkinson disease. Nature Reviews Disease Primers, 3, 17013.
    OpenUrl
  61. ↵
    Power, B. D. & Looi, J. C. 2015. The thalamus as a putative biomarker in neurodegenerative disorders. Australian and New Zealand Journal of Psychiatry, 0004867415585857.
  62. Salimi-Khorshidi, G., Douaud, G., Beckmann, C. F., Glasser, M. F., Griffanti, L. & Smith, S. M. 2014. Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. Neuroimage, 90, 449–468.
    OpenUrlCrossRefPubMedWeb of Science
  63. ↵
    Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., Reiss, A. L. & Greicius, M. D. 2007. Dissociable intrinsic connectivity networks for salience processing and executive control. J Neurosci, 27, 2349–56.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    Sheline, Y. I. & Raichle, M. E. 2013. Resting state functional connectivity in preclinical Alzheimer’s disease. Biological Psychiatry, 74, 340–7.
    OpenUrlCrossRefPubMedWeb of Science
  65. ↵
    Smith, S. M. 2002. Fast robust automated brain extraction. Human brain mapping, 17, 143–155.
    OpenUrlCrossRefPubMedWeb of Science
  66. ↵
    Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., Filippini, N., Watkins, K. E., Toro, R., Laird, A. R. & Beckmann, C. F. 2009. Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci U S A, 106, 13040–5.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    Sporns, O., Tononi, G. & Kötter, R. 2005. The Human Connectome: A Structural Description of the Human Brain. PLoS Comput Biol, 1.
  68. ↵
    Sterling, N. W., Du, G., Lewis, M. M., Dimaio, C., Kong, L., Eslinger, P. J., Styner, M. & Huang, X. 2013. Striatal shape in Parkinson’s disease. Neurobiology of Aging, 34, 2510–2516.
    OpenUrl
  69. ↵
    Strafella, A. P. 2013. Anatomical and functional connectivity as a tool to study brain networks in Parkinson’s disease. Movement Disorders, 28, 411–412.
    OpenUrlCrossRefPubMed
  70. ↵
    Styner, M., Oguz, I., Xu, S., Brechbuhler, C., Pantazis, D., Levitt, J. J., Shenton, M. E. & Gerig, G. 2006. Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM. Insight J, 242–250.
  71. ↵
    Summerfield, C., Junque, C., Tolosa, E., Salgado-Pineda, P., Gomez-Anson, B., Marti, M. J., Pastor, P., Ramirez-Ruiz, B. & Mercader, J. 2005. Structural brain changes in Parkinson disease with dementia: a voxel-based morphometry study. Arch Neurol, 62, 281–5.
    OpenUrlCrossRefPubMedWeb of Science
  72. ↵
    Tinaz, S., Courtney, M. G. & Stern, C. E. 2011. Focal cortical and subcortical atrophy in early Parkinson’s disease. Movement Disorders, 26, 436–441.
    OpenUrlCrossRefPubMed
  73. ↵
    Tziortzi, A. C., Haber, S. N., Searle, G. E., Tsoumpas, C., Long, C. J., Shotbolt, P., Douaud, G., Jbabdi, S., Behrens, T. E. J., Rabiner, E. A., Jenkinson, M. & Gunn, R. N. 2014. Connectivity-Based Functional Analysis of Dopamine Release in the Striatum Using Diffusion-Weighted MRI and Positron Emission Tomography. Cerebral Cortex, 24, 1165–1177.
    OpenUrlCrossRefPubMed
  74. ↵
    Williams-Gray, C. H., Foltynie, T., Brayne, C. E. G., Robbins, T. W. & Barker, R. A. 2007. Evolution of cognitive dysfunction in an incident Parkinson’s disease cohort. Brain, 130, 1787–1798.
    OpenUrlCrossRefPubMedWeb of Science
  75. ↵
    Woolrich, M. W., Behrens, T. E. J., Beckmann, C. F., Jenkinson, M. & Smith, S. M. 2004. Multilevel linear modelling for FMRI group analysis using Bayesian inference. Neuroimage, 21, 1732–1747.
    OpenUrlCrossRefPubMedWeb of Science
  76. ↵
    Woolrich, M. W., Ripley, B. D., Brady, M. & Smith, S. M. 2001. Temporal autocorrelation in univariate linear modeling of FMRI data. Neuroimage, 14, 1370–1386.
    OpenUrlCrossRefPubMedWeb of Science
  77. ↵
    Worsley, K. 2001. Statistical analysis of activation images. Functional MRI: an introduction to methods, 14, 251–270.
    OpenUrl
  78. ↵
    Zhan, Z.-W., Lin, L.-Z., Yu, E.-H., Xin, J.-W., Lin, L., Lin, H.-L., Ye, Q.-Y., Chen, X.-C. & Pan, X.-D. 2018. Abnormal resting-state functional connectivity in posterior cingulate cortex of Parkinson’s disease with mild cognitive impairment and dementia. CNS Neuroscience & Therapeutics, 24, 897–905.
    OpenUrl
  79. ↵
    Zhang, Y., Brady, M. & Smith, S. 2001. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Transactions on Medical Imaging, 20, 45–57.
    OpenUrlCrossRefPubMedWeb of Science
Back to top
Posted August 26, 2019.
Download PDF

Supplementary Material

Data/Code
Email
Neuroimaging biomarkers differentiate Parkinson disease with and without cognitive impairment and dementia
Conor Owens-Walton, David Jakabek, Brian D. Power, Mark Walterfang, Sara Hall, Danielle van Westen, Jeffrey C.L. Looi, Marnie Shaw, Oskar Hansson
medRxiv 19002147; doi: https://doi.org/10.1101/19002147
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area