Substantia nigra ferric overload and neuromelanin loss in Parkinson’s disease measured with 7T MRI

: 1 2 Background: Vulnerability of the substantia nigra dopaminergic neurons in Parkinson’s 3 disease is associated with ferric overload, leading to neurodegeneration with cognitive and 4 motor decline. Here, we quantify iron and neuromelanin-related markers in vivo using ultra- 5 high field 7-Tesla MRI, and examine the clinical correlates of these imaging assessments. 6 7 Methods: Twenty-five people with mild-to-moderate Parkinson’s disease and twenty-six 8 healthy controls underwent high-resolution imaging at 7-Tesla with a T 2 *-weighted sequence 9 (measuring susceptibility-χ and R 2 *, sensitive to iron) and a magnetization transfer-weighted 10 sequence (MT-w, sensitive to neuromelanin). From an independent control group (N=29), we 11 created study-specific regions-of-interest for five neuromelanin- and/or iron-rich subregions 12 within the substantia nigra. Mean R 2 *, susceptibility-χ and their ratio, as well as the MT-w 13 contrast-to-noise ratio (MT-CNR) were extracted from these regions and compared between 14 groups. We then tested the relationships between these imaging metrics and clinical severity. 16 Results: People with Parkinson’s disease showed a significant ~50% reduction in MT-CNR 17 compared to healthy controls. They also showed a 1.2-fold increase in ferric iron loading 18 (elevation of the ∆𝑅 !∗ /∆𝜒 ratio from 0.19±0.058ms/ppm to 0.22±0.059ms/ppm) in an area of 19 the substantia nigra identified as having both high neuromelanin and susceptibility MRI signal 20 in healthy controls. In this region, the ferric-to-ferrous iron loading was associated with 21 disease duration ( b =0.0072, p FDR =0.048) and cognitive impairment ( b =-0.0115, p FDR =0.048). Conclusions:


Introduction 31 32
Iron accumulation and loss of pigmented neuromelanin cells in the substantia nigra are key 33 1 (N1-sign) in healthy controls, and its loss in patients with Parkinson's disease (24). 62 Susceptibility imaging is sensitive to any molecule that generates an MR phase shift. In 63 contrast to SWI, both QSM and R2* provide quantitative measures of susceptibility effects. 64 However, they are unable on their own to differentiate the sources of susceptibility, and 65 therefore an adjunctive method is required. One approach is to capitalise on the fact that 66 ferrous and ferric irons show differences in relaxivity per unit concentration (25), which can 67 be estimated from the R2*-to-χ ratio. We therefore take the R2*-to-χ ratio as a proxy of ferric 68 vs. ferrous loading, and use this ratio as a biomarker for iron accumulation in Parkinson's 69

disease. 70
In this study, we used high resolution dedicated MRI sequences for in vivo imaging of 71 neuromelanin with high resolution (0.08 mm 3 ) and iron (0.34 mm 3 ) using ultra-high field 7T, 72 based on magnetisation transfer (26) and multi-echo T2*, respectively, with nigral regions of 73 interest derived from an independent cohort of age-matched healthy controls (27). The 74 contrast-to-noise ratio (CNR) of the MT signal was interpreted as a measure of neuromelanin 75 content, whereas χ and R2* were used to quantify iron in the tissue. We investigated (1) the 76 molecular changes of neuromelanin and iron in sub-regions of the substantia nigra in 77 Parkinson's disease compared to healthy controls, and (2) the association between imaging 78 markers and cognitive and motor functions in the patient group. We hypothesised that people 79 with Parkinson's disease would show decreased MT signal contrast and increased χ and R2* 80 in the SNpc. Since ferric iron demonstrates higher susceptibility and R2* values than ferrous 81 iron at the same concentration in tissue (25,28) we also expected an increase of the R2*-to-χ 82 ratio in the SNpc.  (33)), Montreal cognitive assessment (MoCA; (34)). Self-rated questionnaires were used to 117 assess anxiety and depression (Hospital Anxiety and Depression Scale; HADS (35)), impulsivity 118 (Barratt Impulsiveness Scale; BIS-11 (36)) and apathy (Apathy Scale; (37)). 119 To create a template and select regions-of-interest, an independent sample of 30 additional 120 age-, sex-and education-matched healthy controls (HCi, 17 male, age 67±8 years) was 121 collected as part of the same study protocol. with ANTs v2.2.0 and FSL v5.0.10 were used for image pre-processing (processing scripts 147 available upon request). QSM maps were estimated from the T2*-w scans using the multi-148 scale dipole inversion algorithm in QSMbox (40) with the pipeline described in (23) but with 149 the regularization parameter set to 10 3.3 . R2* was estimated with the ARLO algorithm (41). To 150 detect the relaxivity effects of ferrous (Fe 2+ ) and ferric (Fe 3+ ) components (25), the ∆ ! * /∆ 151 ratio was measured in each voxel , & , by the ratio of the normalised R2* maps 152 (referenced to the whole brain R2* value ( % * '( )) and the ∆ measured by QSM (MSDI QSM-153 processing references c to the whole brain c value). 154 To calculate SWI maps, for each echo the raw phase was filtered with a homodyne high-pass 155 filter (filter strength: 6.2%). All echoes were then temporally fitted, and a normalized phase 156 mask, % , was generated according to equation 2 of (42) under the assumption that the 157 phase of interest is negative. SWI maps were computed by multiplying the T2*-map with the 158 fourth power of the normalized phase mask ( All MT scans (MT-on and MT-off) were pre-processed as previously described (27). Briefly, 160 data was first bias field corrected. Then the MT-on images (2 repeats for HC and 3 repeats for 161 HCi groups) were combined with the customised antsMultivariateTemplateContruction2 162 function setting the Laplacian sharpening off in the averaging step. 163 The T1-w MP2RAGE image was generated offline from the raw magnitude and phase of the 164 two inversion times (43). Then, were bias-field corrected with ANTs and segmented with 165 SPM12 (v7219) and skull-stripped. 166 167

Analysis of independent control dataset 168
A T1-w driven, cross-modality co-registration pipeline was built to warp the MT and the QSM 169 data from the HCi group into the isotropic 0.5mm ICBM152 T1-w asymmetric brain template 170 (44) (Figure 1). On a within-subject level, the T1-w images were first independently registered 171 to the combined MT-on image (rigid transform, via the MT-off image (27)), and to the first 172 echo of the T2*-w scan (rigid transform using the cross-correlation metric with weight 1 and 173 window radius 4 mm). A T1-w structural template was created from the individual T1-w data 174 using a three step registration process: 1) rigid, 2) affine, and 3) hierarchical nonlinear 175 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Generation of substantia nigra regions-of-interest 198
We applied a semi-automated threshold-based method to segment the substantia nigra 199 independently in the MT-on and QSM images, as it has been validated for the substantia nigra 200 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  The goodness-of-fit for the parametric distributions were calculated and the distribution 258 which showed the lowest Akaike Information Criterion was then used on a general linear 259 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

281
Demographic and clinical characteristics are shown in  image (red outline - Figure 3D; approx. 25% and 31% shared voxels, respectively; Dice=0.3). 292 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

MT images show an overall reduction of signal in the substantia nigra in participants with 293
Parkinson's disease compared to healthy controls. For iron-related measures, in participants 294 with Parkinson's disease, the N1 sign is not evident in the SWI images ( Figure 3D), which is 295 associated with a local increase in iron concentration as observed by the increased signal in 296 QSM and R2* in the Overlap ROI ( Figure 3B,C).  Table 2). The 313 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

332
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  4. Discussion 369 This study reveals the impact of Parkinson's disease on neuromelanin (by MT contrast) and 370 iron (by susceptibility quantities -QSM and R2*), using ultra-high field high-resolution 371 magnetic resonance imaging of the substantia nigra. We defined two main regions-of-372 interest. In healthy individuals, MT-CNR was highest in the inferoventral region of the 373 substantia nigra, whereas increased susceptibility was greatest in the superolateral region. 374 As reported in previous imaging studies (48,52,53), these two regions represent to some 375 degree the difference between the classically defined pars compacta and pars reticulata. 376 The pars compacta has closely packed pigmented neurons, while the pars reticulata shows 377 . We found that the Overlap ROI shares a similar location to the N1-sign 391 obtained with SWI (Dice=0.3), but not completely (red and green outlines in Figure 3). 392 Segmentation was performed on an independent healthy control group using the same 393 imaging strategy as in the group study in MNI space with the purpose of reducing ROI 394 selection bias. Distinguishing types of iron in the brain with MRI in vivo is challenging. Here we referred to 410 R2*-to-χ ratio as an indirect estimate of the ferric (Fe 3+ ) to ferrous (Fe 2+ ) load in the brain. Our 411 results show that, in patients, there is a significant increase of this ratio within the ventral 412 portion of the SNpc. We differentiated iron markers in MRI in Parkinson's patients by 413 combining both magnitude and phase data from a single MRI acquisition. The R2*-to-χ ratio 414 as a proxy for ferric to ferrous loading has been previously validated in phantom work using 415 MRI in (25). However, in future, post-mortem MRI imaging together with histological iron 416 staining is still required to validate this measure in healthy and diseased tissue. 417 Reports of the correlation of imaging and clinical outcomes in Parkinson's disease are diverse. 418 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 20, 2021. ; https://doi.org/10.1101/2021.04.13.21255416 doi: medRxiv preprint We demonstrate the usefulness of a neuromelanin-iron multimodal imaging approach. 451 However, the addition of diffusion scalars similar to (73) at high-resolution would provide 452 additional quantitative values of the microstructural integrity (e.g. myelin) and intra-to extra-453 cellular water diffusion variations within the nigra compartments that could aid 454 characterization of the neurodegeneration in Parkinsonian disorders. 455 Hemispheric effects were found on both control and patient data, particularly evident on 456 smaller ROIs (e.g. results for the SWIROI in Figure 5). Yet, we found no associations between 457 imaging and clinical motor symptom laterality. The hemispheric asymmetry on the imaging 458 results might be related to the known asymmetry in the B1 + transmit profile which occurs 459 when using transmit-receive coils such as the coil used in our study in circular-polarised mode 460 (74,75). In our analysis, we controlled for hemisphere effects by either adding a fixed variable 461 of no interest on the general linear modelling statistics or by averaging the data across both 462

hemispheres. 463
Magnetization transfer imaging shows specificity to the dopaminergic neurons of the SNpc 464 but its biological relationship with neuromelanin is still elusive (76-78) as a contrast-to-noise 465 measurement is only semi-quantitative. In order to better quantify nonexchangeable protons 466 attached to macromolecules like neuromelanin, kinetic model parameters describing 467 magnetization dynamics via a two-pool model could be used (79,80). Susceptibility and R2* 468 mapping have the advantage of quantifying iron-related microscopic field changes, and so 469 provide a more biologically meaningful clinical indicator of the pathogenesis. Therefore, the 470 combination of these two imaging techniques provide a more holistic understanding of the 471 neurodegeneration in the dopaminergic pathway within the substantia nigra and a more 472 sensitive approach in differentiating patients with Parkinson's disease from healthy controls. 473

474
In conclusion, the complementary information from iron and neuromelanin high-resolution 475 ultra-high field 7T imaging allowed the detection of ferric (Fe 3+ ) overload and neuromelanin 476 loss on Parkinson's patients in an area of the substantia nigra with both high neuromelanin 477 and susceptibility MRI signal in normal controls. In addition, patients at later stages of the 478 disease show higher ferric loading and lower cognitive function which is in accordance with 479 histopathological findings. Our results suggest that ultra-high field 7T MT-weighted and T2*-480 weighted imaging are sensitive to the detection of molecular change in Parkinson's disease, 481 and may therefore support future experimental medicines strategies. 482 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 20, 2021.

646
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 20, 2021. ; https://doi.org/10.1101/2021.04.13.21255416 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 20, 2021. ; https://doi.org/10.1101/2021.04.13.21255416 doi: medRxiv preprint