Abstract
Background Patients with Parkinson’s disease undergo a loss of melanized neurons in substantia nigra pars compacta and locus coeruleus. There are very few in vivo studies of locus coeruleus pathology in Parkinson’s disease with magnetic resonance imaging. Existing studies have used varying methodologies, and reproducibility has not been established for any approach.
Objectives We used neuromelanin-sensitive magnetic resonance imaging to examine the reproducibility of locus coeruleus and substantia nigra pars compacta degeneration in discovery and validation cohorts.
Methods Two cohorts, discovery and validation, recruited from the Emory Movement Disorders Clinic were scanned on two different magnetic resonance imaging scanners. In cohort 1, imaging data from 19 controls and 22 Parkinson’s disease patients were acquired with a Siemens Trio 3 Tesla scanner using a 2D gradient echo sequence with magnetization transfer preparation pulse. Cohort 2 consisted of 33 controls and 39 Parkinson’s disease patients were scanned on a Siemens Prisma 3 Tesla scanner with a similar imaging protocol. Nigral and locus coeruleus volumes were segmented in both cohorts.
Results Nigral volume (Cohort 1: p=0.0148; Cohort 2: p=0.0011) and locus coeruleus volume (Cohort 1: p=0.0412; Cohort 2: p=0.0056) were significantly reduced in the Parkinson’s disease group as compared to controls in both cohorts.
Interpretation Substantia nigra pars compacta volume and locus coeruleus volume were significantly reduced in Parkinson’s disease compared to controls in both cohorts. This imaging approach robustly detects Parkinson’s disease effects on these structures suggesting its use as a potential diagnostic biomarker of Parkinson’s disease.
1. Introduction
Parkinson’s disease (PD) is a heterogeneous neurodegenerative disorder with a variety of motor and non-motor symptoms that can be clinically challenging to diagnose and manage. By the time of symptom onset and clinical diagnosis, there is an estimated 50-70% of dopaminergic neuronal loss in the substantia nigra compacta (SNc).1-3 There are currently no effective neuroprotective interventions to stop PD neurodegeneration. Historically, SNc degeneration is known as the pathological hallmark of PD, but it has been proposed that alpha synuclein protein aggregates, in the form of Lewy bodies, spread across the brain in multiple stages affecting other structures, including locus coeruleus (LC), before the SNc.4, 5 Despite insights gained through postmortem study of PD neuropathology, there is a lack of reliable biomarkers that can predict, diagnose, and track progression in PD. Neuroprotection trial designs could be improved with brain imaging markers of neurodegeneration, which might assist participant selection and serve as surrogate outcome measures. Therefore, dependable and reproducible neuroimaging biomarkers are needed to detect and quantify PD neuropathology in vivo.
Melanized neurons in SNc and LC can be visualized in vivo with neuromelanin-sensitive MRI (NM-MRI) sequences using either incidental magnetization transfer effects from an interleaved multi-slice turbo spin echo acquisition6 or explicit magnetization transfer effects generated by magnetization transfer preparation pulses,7-9 Magnetization transfer contrast (MTC) colocalizes with melanized neurons in LC and SNc.10, 11 NM-MRI approaches using incidental12 or explicit13-15 magnetization transfer effects have been found to exhibit high scan-rescan reproducibility in controls and gradient echo-based approaches with explicit magnetization transfer effects exhibit the highest reproducibility.15
NM-MRI can be used to assess PD-related reductions in tissue neuromelanin content in SNc and LC. Application of NM-MRI approaches based on incidental magnetization transfer effects to image depigmentation in SNc and LC has revealed PD-related reductions in NM-MRI contrast ratios in SNc or LC,6, 16, 17 nigral volume,18 or area of SNc in a single slice.19, 20 Similar reductions have been observed in nigral volume7 and contrast.21 Nigral regions of interest, derived from NM-MRI images, have also been used to examine PD-related microstructural changes22, 23 or iron deposition24, 25 in SNc.
Replication of imaging markers in multiple cohorts is a crucial step in biomarker development and reliable detection of neuromelanin loss in SNc and LC may provide diagnostic imaging markers for PD. Reductions in nigral volume have been replicated in separate cohorts using explicit magnetization transfer effect-based NM-MRI approaches.26 However, PD-related reductions in LC volume loss have not been replicated and this deficiency may be due to the size and stature of LC, a small rod-shaped structure approximately 2 mm in diameter and 15 mm long.27 Here, we aim to remedy this deficiency and apply a NM-MRI approach based on explicit magnetization transfer effects to measure volume loss in LC as well as in SNc in PD patients and controls in separate discovery and validation cohorts acquired on different MRI scanner models.
2. Methods
2.1 Research Participants
Two cohorts, discovery and validation, were used in this study and were scanned on two different MRI scanners as described below. Cohort 1 data was collected from 2012-2014 and included 19 controls and 22 PD patients. Cohort 2 consisted of 33 healthy controls (HC) and 39 PD patients with data collected from 2015-2016. All research participants were recruited from the Emory University Movement Disorders Clinic under an institutional review board approved protocol with informed written consent. Controls were recruited from the community and the Emory Alzheimer’s Disease Research Center control population. PD patients fulfilled the Movement Disorders Society clinical diagnostic criteria,28 and diagnosis was established by a fellowship-trained movement disorders neurologist at the Emory Movement Disorders Clinic. PD participants had early to moderate disease with Unified Parkinson’s Disease
Rating Scale Part III (UPDRS-III) motor score :−25 in the practically defined ON state. Disease duration in years and levodopa equivalent daily dose (LEDD) were also determined for PD participants. Cognition was assessed using the Montreal Cognitive Assessment (MoCA).29 The Non-motor Symptoms Questionnaire (NMSQ)30 was used to assess non-motor symptoms. Symptoms of rapid eye movement (REM) behavior disorder were assessed using the REM Sleep Behavior Disorder Screening Questionnaire (RBD-SQ).31
2.2 MRI Acquisition
MRI data for Cohort 1 were acquired with a Siemens Trio 3 Tesla scanner (Siemens Medical Solutions, Malvern, PA, USA) at Emory University with a 12-channel receive-only head coil. NM-MRI data was acquired using a 2D magnetization transfer (MT) prepared gradient echo (GRE) sequence8, 9: echo time (TE) / repetition time (TR) = 2.68 ms / 337 ms, slice thickness 3 mm, in plane resolution 0.39×0.39 mm2, field of view (FOV) = 162 mm × 200 mm, flip angle (FA) = 40°, 470 Hz/pixel bandwidth, 15 contiguous slices, and MTC preparation pulse (300°, 1.2 kHz off resonance, 10 ms duration), 7 measurements, scan time 16 minutes 17 seconds. For registration from subject space to common space, T1 magnetization-prepared rapid gradient echo (MP-RAGE) sequence was acquired with the following parameters: TE/TR= 3.02 ms/2600 ms, inversion time = 800 ms, FA=8°, voxel size = 1.0×1.0×1.0 mm3.
Cohort 2 was scanned with a Siemens Prisma 3 Tesla scanner using a 64-channel receive-only coil. NM-MRI data were acquired using a 2D GRE sequence with a MT preparation pulse: TE/TR = 3.10 ms /354 ms, 15 contiguous slices, FOV = 162×200 mm2, in plane resolution = 0.39 × 0.39 ×mm2), 7 measurements, FA = 40°, 470 Hz/pixel receiver bandwidth, and MTC pulses (300°, 1.2 kHz off resonance, 10 ms duration), scan time 17 minutes 12 seconds. For registration, structural images were acquired with a MP-RAGE sequence: TE/TR = 2.46 ms/1900 ms, inversion time = 900 ms, FA = 9°, voxel size = 0.8×0.8×0.8 mm3.
On the sagittal T1 images for both cohorts, the NM-MRI scan slices were positioned perpendicular to the dorsal edge of the brainstem at midline along the fourth ventricle, starting from the lower pons (below the most caudal extent of LC), with slices covering SNc and LC.
2.3 Image Processing
MRI data was processed using the FMRIB Software Library (FSL). A transformation was derived between each individual’s T1-weighted image and 2 mm Montreal Neurological Institute (MNI) T1-space using FMRIB’s Linear Image Registration Tool (FLIRT) and FMRIB’s Nonlinear Image Registration Tool (FNIRT) in the FSL software package using the following steps.32, 33 The T1-weighted image was brain extracted using the brain extraction tool (BET). Next, an affine transform was used to align the brain extracted T1-weighted images with the MNI brain extracted image. Finally, a nonlinear transformation was used to generate a transformation from individual T1-weighted images to T1-weighted MNI T1-space.
For each participant individual NM-MRI measurements were corrected for motion by registering the seven measurements to the first image using a rigid-body transform in FMRIB Linear Image Registration Tool (FLIRT) and then averaged. Next, a transform was derived between each individual’s T1-weighted image and the averaged NM-MRI image with a boundary-based registration cost function.
SNc and LC volumes were segmented using an automated thresholding method. To ensure consistent placement of reference regions of interest (ROIs), a reference ROI in the cerebral peduncle was created using the MNI template and, for each subject, the cerebral peduncle ROI was transformed to individual NM-MRI images using transforms described in previous paragraphs. This process ensured that the reference ROI was placed in similar locations for each subject. The mean (denoted μref), and standard deviation (σref) of the signal intensities were measured in the reference ROI. Next, SNc and LC were segmented using probabilistic standard space masks.34 These atlases were transformed from standard space to individual NM-MRI images, thresholded at a level of 5%, and dilated. Voxels in SNc and LC ROIs with intensity > μref+3.9 σref and > μref+2.8σref were considered to be part of LC and SNc, respectively.
2.4 Statistical Analysis
All statistical analyses were performed using IBM SPSS Statistics software version 24 (IBM Corporation, Somers, NY, USA) and results are reported as mean ± standard error. A p value of 0.05 was considered significant for all statistical tests performed in this work. Normality of SNc and LC volumes was assessed using the Shapiro-Wilk test for each group and all data was found to be normal.
For demographic data, independent samples t-test was used to assess differences in age and education and Chi square was used to examine differences in gender and race between controls and PD in each cohort. Group means for UPDRS-III score, disease duration, levodopa equivalents, MoCA, NMSQ, and RBD-SQ were compared using a two-tailed Welch’s t-test in each cohort.
In cohort 1, the effect of group (PD, control) was tested with separate analysis of covariance (ANCOVA) for SNc and LC volume controlling for age and education. In cohort 2, differences in SNc volume and LC volume between control and PD groups were compared using a two-tailed Welch’s t-test. Pearson’s correlation was used to assess the relationship between SNc volume with UPDRS-III, LEDD, MoCA and LC volume in PD for cohort 1, cohort 2, and both cohorts combined. Pearson’s correlation was used to assess the relationship between LC volume with UPDRS-III and MoCA for cohort 1, cohort 2, and both cohorts combined. Age was a treated as a covariate in all correlations.
3. Results
In cohort 1, significant differences were in age (p<10−3) and education (p=0.015) with controls being older and having higher levels of education than the PD group. No differences were seen between groups in gender (p=0.453) or in MoCA score (p=0.556). The PD group had significantly higher UPDRS-III (p <10−3), RBD-SQ (p=0.002), and NMSQ (p<10−3) scores compared to the control group. In cohort 2, no difference was seen between PD and control groups in age (p=0.894), race (p=0.402), gender (p=0.500), education (p=0.661), RBD-SQ score (p=0.087), or MoCA score (p=0.621). A significant difference was seen in NMSQ (p=0.0006) and UPDRS-III (p<10−3) scores between the two groups. Demographic information is summarized in Table 1.
Figures 1 and 2 show LC and SNc in mean MTC images of both cohorts. In cohort 1, the effect of group on volume was assessed with separate ANCOVAs for each ROI (SNc, LC) with age and education as covariates. In SNc, results revealed a significantly smaller volume in the PD group relative to the control group (Control: 474 mm3 ± 31 mm3; PD: 340 mm3 ± 28 mm3; F=8.031; p=0.007). Similarly, a reduction in LC volume was seen in the PD group relative to the control group (Control: 6.9 mm3 ± 0.7 mm3; PD: 4.4 mm3 ± 0.7 mm3; F=8.031; p=0.033). In cohort 2, Welch’s t-test was used to examine group differences in SNc volume and LC volume. SNc (Control: 429 mm3 ± 20 mm3; PD: 329 mm3 ± 17 mm3; t=3.370; p=0.0002) and LC (Control: 8.0 mm3 ± 0.6 mm3; PD: 5.2 mm3 ± 0.6 mm3; t=3.306; p=0.0008) volumes were significantly lower in the PD group compared to the control group in cohort 2. These comparisons are shown in Figure 3.
In the PD group, MoCA showed a significant positive association with SNc volume in cohort 1 (p=0.012, r=0.478) and both cohorts (p=0.040, r=0.228) but not in cohort 2 (p=0.319, r=0.079). SNc volume showed no significant associations with UPDRS-III or disease duration in the PD group in cohort 1, cohort 2 or both cohorts (ps>0.436). There were no significant correlations between MoCA, UPDRS-III and disease duration with LC volume in cohort 1, cohort 2, or both cohorts combined for the PD group (ps>0.272). SNc volume showed a significant positive association with LC volume in PD in both cohorts (p=0.012, r=0.296) but no significance in HC (p=0.434, r=-0.025). In the PD group, SNc and LC volumes had a significant positive correlation in cohort 2 (p=0.024, r=0.333) but not in cohort 1 (p=0.138, r=0.243). These correlations are shown in Figure 4.
In cohort 1 SNc volume outperformed LC volume as a diagnostic imaging marker of PD. The area under the receiver operating characteristic (ROC) curve (AUC) for SNc volume was 0.756 [standard error (SE): 0.078, 95% confidence interval (CI): 0.603-0.909, p=0.005] while the AUC for LC volume was 0.644 [SE: 0.088, 95% CI: 0.471-0.816, p=0.117]. SNc volume and LC volume performed similarly in cohort 2 as diagnostic imaging markers of PD. LC volume had an AUC of 0.752 [SE: 0.063, 95% CI: 0.629-0.876, p=0.001] and SNc volume had an AUC of 0.749 [SE: 0.062, 95% CI: 0.627-0.871, p=0.001].
4. Discussion
This study examined PD-related loss of NM-MRI contrast in LC and SNc in separate discovery and validation cohorts, using two different scanner models. We observed significant volume loss in both LC and SNc in the PD groups as compared to controls in both cohorts. In cohort 2, SNc volume and LC volume performed similarly as diagnostic imaging markers of PD. This is the first report, to our knowledge, of reproducible detection of PD-associated LC volume loss using the same NM-MRI approach in multiple cohorts. It is also the first report of simultaneous imaging of LC and SNc detecting PD effects in both structures in discovery and validation cohorts. The NM-MRI pulse sequence and image processing methods that were used in this study have previously established high scan-rescan reproducibility.13, 15 To increase reproducibility and ensure consistent placement of ROIs, reference regions and thresholded regions were defined using standard space ROIs and transformed to each individual’s NM-MRI image. These ROIs were then used to define thresholded and reference regions in the SNc and LC segmentation procedure. The results of this study, therefore, support the utility of these NM-MRI methods as a robust and reproducible approach to measure PD neuropathology in vivo.
The two Parkinsonian cohorts used in this study had similar motor impairment, but different disease durations. The mean and standard error of nigral volume was strikingly consistent across both cohorts. Similarly, no difference was observed in LC volume between cohorts (p=0.081, nonsignificant). The consistency in SNc and LC volume loss may be due to similarity of NM-MRI protocols used in this study. This is in contrast to an earlier study which found LC volume to be dependent on scanner,26 which used different implementations of NM-MRI protocols to image LC and the variability in scanning protocols may explain the differences in LC volume.
The ROC analysis found AUC of nigral volume in both cohorts and LC volume in cohort 2 to have comparable AUCs to previously published nigral diagnostic imaging markers of PD.21, 26, 35-38 However, the performance of cohort’s 1 LC volume as a diagnostic marker of PD was significantly below AUCs reported in earlier studies.21, 26, 35-38 Reduced performance of LC volume in cohort 1 may be due to several factors. First, the Siemens Trio MRI scanner used in cohort 1 has a higher noise profile than the Siemens Prisma MRI scanner used in cohort 2 and elevated noise likely reduced efficacy of LC volume as a diagnostic marker in cohort 1. Second, the control group was significantly older than the PD group in cohort 1. Neuromelanin peaks in LC between age 50-60 and declines after age 60.39-41 Thus, this age difference likely reduced the effect size by reducing the mean LC volume in the control group.
Nigrosome-1 is the subregion of SNc that experiences the greatest reductions in melanized neurons42 and is found in the lateral-ventral portion of SNc.1 This region can be visualized in T2- or T2*-weighted images as a hyperintense region in the substantia nigra and the loss of hyperintensity of nigrosome-1 in T2- or T2*-weighted images of PD patients is used as a diagnostic imaging marker.43 In NM-MRI images, nigrosome-1 has been localized to the posterior portion of SNc.44 We observed reductions in NM-sensitive contrast in the posterior regions of SNc (see Figure 2) in both cohorts. This is in agreement with earlier studies that found volume loss or reductions in neuromelanin-sensitive contrast in the posterior portion of SNc45 or in the lateral-ventral tier of SNc.21 Studies examining PD-related changes using other MRI contrasts have found increases in free water,46, 47 associated with neuronal loss, in posterior nigral ROIs as well as elevated iron levels in the lateral-ventral tier.48, 49 Taken together, these changes may manifest from neurodegeneration in nigrosome-1.
Animal studies using 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) to induce Parkinsonism in mouse and primate models found that lesioning LC increases the amount of dopaminergic neurons lost in SNc50-53 and histology in humans found elevated levels of dopamine loss in PD correlates with reduced noradrenaline.54 Thus, integrity of LC neurons may protect dopaminergic neurons in SNc. In cohort 2, SNc and LC volume were significantly correlated in the PD group and this result may indicate that LC integrity plays a key role in Parkinsonian neurodegeneration in SNc. However, given the small sample size, different disease durations between cohorts, and the correlation not being replicated in cohort 1, these results should be interpreted with caution.
There are several caveats to this study. First, UPDRS-III was measured in the on-medication state. Earlier studies found a significant relationship between nigral volume and UPDRS score26 and lack of correlation between nigral volume and UPDRS-III score may be due to measurement in the on medication state or lack of inclusion of PD patients with more severe motor symptoms. Second, the control group in cohort 1 was significantly older and more educated than the PD group. This may have reduced the effect size for LC volume.
5. Conclusion
NM-MRI approaches based on explicit magnetization transfer effects have already demonstrated high scan-rescan reproducibility.13-15 The current findings provide additional evidence that NM-MRI robustly detects PD-related reductions in volume in both LC and SNc. Therefore, further investigation of these SNc and LC volumes from NM-MRI as candidate PD biomarkers is warranted, both individually and as part of multivariate marker profiles using iron-sensitive and diffusion-sensitive MRI modalities.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Acknowledgements
This work was supported by 1K23NS105944-01A1 (Huddleston) from the National Institutes of Health/National Institute of Neurological Diseases and Stroke, Michael J. Fox Foundation grants (MJFF-010556 and MJFF-010854), and the American Parkinson’s Disease Foundation Center for Advanced Research at Emory University (Huddleston).
Footnotes
E-mail: daniel.huddleston{at}emory.edu