Cohort Profile: the Oxford Parkinson’s Disease Centre Discovery Cohort Magnetic Resonance Imaging sub-study (OPDC-MRI)

Clinical data

Participants receive extensive assessment in designated research clinics as part of their participation in the OPDC Discovery cohort. The assessment, performed by a nurse and neurologist, includes a structured general medical interview, detailed characterisation of motor and non-motor features, and cognitive assessment (see ² for details). Patients are followed up clinically every 18 months, while controls only have a baseline clinic visit. In a research clinic, we cannot formally diagnose patients who convert from RBD to Parkinson’s (or another neurodegenerative disorder). However, where history and examination suggest conversion to neurodegenerative disease, we alert the treating clinicians who then establish the diagnosis and ensure clinical management is in place. On the day of scanning, an additional UPDRS III assessment was performed (“off” in Parkinson’s). Table 1 summarises the main demographic and clinical characteristics of the OPDC-MRI sub-study cohort.

Imaging data

Scanning was performed at the Oxford Centre for Clinical Magnetic Resonance Research (OCMR) using a 3T Siemens Trio MRI scanner (Siemens, Erlangen, Germany) equipped with a 12-channel receive-only head coil. The neuroimaging protocol includes both structural and functional sequences and lasts approximately 45-50 minutes.

Within the allocated time, five modalities were always acquired (core sequences), while the remaining time was used to experiment with novel sequences, which changed during the study. We report here the details of two of them, introduced at a later stage based on new information becoming available in the field, and which showed promising preliminary results. We describe the other additional sequences we experimented in the supplementary material.

For each modality, pipelines for preprocessing and analysis of MRI data were developed and applied to the data to extract single-subject imaging variables (both summary and voxel-wise measurements). Here below we describe for each modality the rationale for acquisition, as well as the main preprocessing steps performed on the images to derive the measures of interest.

The main acquisition parameters and the number of available scans per group for each of the core sequences are listed in Table 2.

T1-weighted MRI

The T1w MPRAGE (Fig. 1a-1d) offers very good contrast across tissue classes: grey matter (GM), white matter (WM), and cerebrospinal fluid (CSF). It is primarily used to study GM structural macroscopic tissue in both cortical and subcortical regions. In Parkinson’s, cortical morphology in cognitively intact patients is generally reported to be normal or mildly altered, while impaired cognition and dementia in Parkinson’s have been found to be associated with more severe patterns of cortical atrophy^6,13-15. Non-motor symptoms have been also associated with structural changes in specific related brain networks (for a review see ⁶). With this sequence we aimed to investigate potential differences in GM density across groups or relationships between GM and clinical variables, as well as taking into account the effect of structural changes when analysing functional MRI data.

• Download figure
• Open in new tab

Figure 1.

Structural sequences and related processing. T1weighted MRI (a) is brain extracted and bias field corrected (b) to perform tissue-type segmentation (GM in red, WM in yellow, CSF in blue) (c) and subcortical structures segmentation (d). T1-weighted white matter nulled (e) allows better contrast in the subcortical structures, which allows, for example, the segmentation of the Globus Pallidus (zoom) into its internal (GPi, yellow) and external (GPe, light blue) portions (f). FLAIR images (g) are used to detect and quantify white matter hyperintensities (red) (h).

The fsl_anat pipeline was used to perform brain-extraction, bias-field correction and 3-class (GM, WM, CSF) tissue-type segmentation with FAST¹⁶. The resulting images were then used to perform voxel-based morphometry analyses¹⁷ and to generate voxel-wise confound regressors for resting state functional MRI statistical analyses. Subcortical structure segmentation was performed with FIRST¹⁸, obtaining 3D meshes for each structure for each subject (used to perform vertex analysis), as well as volumetric images used to calculate volumes or as masks to extract average values from other modalities.

T1-weighted images with suppressed white matter signal (T1 WM nulled) show enhanced contrast across the basal ganglia structures (Fig.1e). This allows the detection of, for example, the boundary between the internal and external globus pallidus and a finer segmentation of this structure using MIST¹⁹ (Fig.1f).

T2-weighted Fluid-attenuated inversion recovery (FLAIR)

This sequence is commonly used to detect white matter abnormalities like leukoencephalopathies, demyelinating diseases and abnormalities of vascular origin. Regarding vascular pathology, it is used to detect white matter hyperintensities (WMH) (Fig. 1g,1h). With this sequence we wanted to characterise the amount and distribution of WMH in the cohort, look for differences across groups and assess the possible relationship with cardiovascular risk.

FLAIR images were brain extracted and bias field corrected using FAST¹⁶. WMHs were then automatically segmented on FLAIR images with BIANCA²⁰, a supervised segmentation tool which assigns to each voxel a probability of being a lesion, based on their intensity in FLAIR and T1 and their location (details and training dataset are online²¹). The probabilistic output was thresholded and restricted to the voxels located within a mask created from the T1-weighted scans (using the command make_bianca_mask), which excluded the cortex and subcortical structures. The total WMH volume was then calculated for each scan.

Diffusion-weighted MRI (dMRI)

Diffusion MRI is used to study the microstructure of brain WM in vivo. Also, it serves to virtually reconstruct putative WM tracts. The basis of dMRI is measurement of the random motion of water molecules (diffusion), which has a preferential orientation in the WM (i.e. is less restricted along than across the axons). It follows that the preferential direction of water diffusion is related to fibre orientation. The amount of diffusion directionality and restriction can inform about the microstructural environment of a voxel under study.

Diffusion-weighted images were acquired along 60 isotropically distributed diffusion directions (b-value of 1000 s/mm²). Five additional images were acquired without diffusion weighting (b=0 s/mm²). B0 inhomogeneity for diffusion imaging was measured using a dual-echo GRE sequence and the resulting phase and magnitude images were processed to produce field maps for correction of inhomogeneity-induced distortions.

Correction for b0-associated and eddy current-related distortion, as well as participant’s movement, were performed using EDDY²². EDDY uses a generative probabilistic model to estimate inter- and intra-volume movements, displacements caused by field inhomogeneity, and distortions caused by eddy currents induced by the diffusion gradients. Additionally, automatic artefact rejection replaces slice drop-outs with model estimates^23,24. The resultant 4D diffusion data were then fed into dtifit, which fits a diffusion tensor model at each voxel^25,26 and generates maps of tensor-derived measures to assess WM microstructural integrity: fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD). Finally, we ran bedpostX to produce fibre orientation estimates and their respective uncertainties. This model can be used to perform probabilistic tractography for reconstructing WM pathways and assess their properties such as structural connectivity.

Resting-state functional MRI (rfMRI)

rfMRI is used to investigate brain function without requiring the subject to undertake a specific task. Although potentially harder to interpret than task fMRI, rfMRI is not affected by subject’s performance or compliance and allows to study different resting state networks, i.e. sets of brain regions sharing a common time-course of spontaneous fluctuations that have been associated with specific brain functions²⁷. Among those, the basal ganglia network²⁸ (supplementary Figure S1) is of particular interest for the OPDC-MRI cohort. In 28 iPD patients we also repeated the rfMRI sequence in the “on-state” 60–90 min after taking their own dopaminergic medication¹¹.

Firstly, images were motion corrected with MCFLIRT and the six rigid-body parameter time series extracted for each subject were used for subsequent cleaning. Mean relative displacement was also calculated to potentially exclude subjects with excessive motion and as possible confound metric in further analyses. Images were brain extracted, corrected for B0 inhomogeneities using field maps, spatially smoothed using a Gaussian kernel of FWHM of 6 mm, and temporally filtered using a high-pass filtering of 150 s. Single-subject probabilistic independent component analysis (ICA) was then performed with MELODIC²⁹ with automated dimensionality estimation, followed by automatic component classification with FSL-FIX^30,31 to identify and regress out the contribution of the artefactual components reflecting non-neuronal fluctuations (FIX training dataset available online²¹). The pre-processed functional data were registered to the individual’s structural scan and standard space images using FLIRT and FNIRT^32,33, using boundary-based-registration. Single-subject resting state networks were derived with dual regression³⁴ and compared across subjects. The spatial maps used as spatial regressors in dual regression can be derived either from group-level ICA on the study-specific data¹¹ or using an external template (for more details regarding this choice see ³⁵). Supplementary Figure S1 shows the basal ganglia network map, part of the template used for the analyses. The full templates used for the studies published so far with data from the OPDC cohort are available online²¹.

Regarding non-core sequences acquired in the OPDC-MRI sub-study, Table 3 shows the main acquisition parameters and the number of available subjects per group for two sequences introduced in 2016, which showed particularly promising results. In the Supplementary material we provide more information on the other sequences experimented in this study, which include quantitative T1 and T2 mapping, diffusion weighted imaging of the substantia nigra and multi-echo T2*-weighted images of the substantia nigra.

Neuromelanin-sensitive MRI (NM-MRI)

The monaminergic neurons in the substantia nigra (SN) and locus coeruleus (LC) are rich in neuromelanin, a dark pigment that gives these structures their distinct colour. Neuromelanin is detectable with MRI as a hyperintense signal using modified T1-weighted sequences^36,37, which exploit the paramagnetic properties of the pigment, due to its iron content. According to pathological studies in Parkinson’s³⁸, the SN and LC are affected early in the neurodegeneration process³⁹, making them an interesting target for the development of neuroimaging biomarkers. Studies using NM-MRI found a reduction of the hyperintense signal on NM-MRI in SN and/or LC in patients in patients with established Parkinson’s^40,41 and RBD^42-44. These promising findings drove the inclusion of this sequence in our protocol.

NM-MRI images acquired in OPDC-MRI were bias field corrected using FAST¹⁶ and the transformations to the individual’s structural scan (T1w) and standard space were calculated using FLIRT and FNIRT^32,33. We also defined two reference regions of interest (ROIs) in MNI space (one for SN and one for LC) and used the average intensities within the ROI (registered in individual subject space) as normalization factors in subsequent analyses.

To also extract quantitative information on these images, we developed a segmentation method to automatically quantify the hyperintense signal in SN and LC (Figure 2a and 2b, preliminary results in ⁴⁵). The analysis on the whole sample is currently ongoing and support data related to the analyses (e.g. reference ROIs) will be available online²¹.

• Download figure
• Open in new tab

Figure 2.

Novel sequences: Neuromelanin-sensitive (NM-MRI, left) and susceptibility weighted imaging (SWI, right). Examples of segmentation of (a) the substantia nigra and (b) the locus coeruleus as hyperintense areas on NM-MRI. (c) Examples of the presence (top) absence (bottom) of the ‘swallow tail’ sign on SWI.

T2*-weighted images / susceptibility weighted imaging (SWI)

Susceptibility-weighted imaging (SWI) uses tissue magnetic susceptibility differences to enhance contrast in MRI. This is achieved by using the phase image in addition to the magnitude of T2*-weighted images. The phase image contains information about local susceptibility changes between tissues, which can be useful in measuring iron content. There are numerous neurologic disorders that can benefit from a sensitive method that monitors the amount of iron in the brain, whether in the form of deoxyhemoglobin, ferritin, or hemosiderin⁴⁶. In Parkinson’s, SWI has recently emerged as a promising sequence for evaluating the integrity of the substantia nigra⁴⁷. In healthy subjects, the dorsolateral SN shows an area of signal hyperintensity, corresponding to nigrosome-I. This feature, described as a ‘swallow tail’ appearance⁴⁸, is lost in Parkinson’s, since nigrosome-I is affected early by synuclein degeneration. Given the promising evidence for the ‘swallow tail’ to be a candidate biomarker for Parkinson’s, we included SWI in our protocol.

The following pre-processing was applied to the T2* images to obtain the final SWI image: macroscopic phase artefacts removal was performed by high-pass filtering the phase images using a 50×50 FWHM window in Fourier space (window size selected empirically to suppress artefacts in the midbrain caused by nearby aerated structures). Then, paramagnetic phase components only were taken to the 4th power and multiplied with the magnitude images.

The presence/absence of the dorsal nigral hyperintensity, the ‘swallow tail’ sign, was visually rated for each subject (Figure 2c). We also calculated the transformations from the individual’s SWI to structural scan (T1w) and standard space using FLIRT and FNIRT^32,33 to perform group comparisons. To extract quantitative information on these images, we also developed a method to automatically quantify the nigrosome-I hyperintense signal⁴⁹. The analysis on the whole sample is currently ongoing and support data related to the analyses will be available online²¹.

Long-term follow-up

These baseline data already represent a rich source of data from a deeply phenotyped cohort. A key aspect of the OPDC-MRI cohort, however, is the longitudinal follow-up, which is ongoing. Clinical longitudinal data are acquired in the Discovery cohort every 18 months, and therefore, we will be able to use them to relate baseline imaging with clinical progression. This will potentially allow us to stratify patients and at-risk individuals and predict their progression. Information about conversion to Parkinson’s of at-risk individuals will also be available, providing the ultimate validation of potential biomarkers. In 2015 we also commenced the acquisition of longitudinal follow-up MRI after 5 years from baseline for Parkinson’s and HC, and after 2.5-3 years for RBD. This difference in time after baseline was chosen because we predict that a number of the RBD patients will convert to Parkinson’s, so we aim to collect multiple data points prior to conversion. In this way we hope to assess their trajectory during the prodromal phase and capture the associated brain changes.