APOE-ε4-related differences in left thalamic microstructure in cognitively healthy adults

APOE-ε4 is a main genetic risk factor for developing late onset Alzheimer’s disease (LOAD) and is thought to interact adversely with other risk factors on the brain. However, evidence regarding the impact of APOE-ε4 on grey matter structure in asymptomatic individuals remains mixed. Much attention has been devoted to characterising APOE-ε4-related changes in the hippocampus, but LOAD pathology is known to spread through the whole of the Papez circuit including the limbic thalamus. Here, we tested the impact of APOE-ε4 and two other risk factors, a family history of dementia and obesity, on grey matter macro- and microstructure across the whole brain in 165 asymptomatic individuals (38–71 years). Microstructural properties of apparent neurite density and dispersion, free water, myelin and cell metabolism were assessed with Neurite Orientation Density and Dispersion (NODDI) and quantitative magnetization transfer (qMT) imaging. APOE-ε4 carriers relative to non-carriers had a lower macromolecular proton fraction (MPF) in the left thalamus. No risk effects were present for cortical thickness, subcortical volume, or NODDI indices. Reduced thalamic MPF may reflect inflammation-related tissue swelling and/or myelin loss in APOE-ε4. Future prospective studies should investigate the sensitivity and specificity of qMT-based MPF as a non-invasive biomarker for LOAD risk.

www.nature.com/scientificreports/ more promising 32,33 . However, it is also possible that limbic regions other than the hippocampus may play an important role in the development of LOAD. Notably, it has been recognised for a while that LOAD pathology may spread through the whole of the Papez circuit and may critically involve the limbic thalamus 4 . For instance, neurofibrillary accumulations in the anterodorsal thalamic nucleus have been found at the same time as those in the hippocampus in LOAD brains 34 and reduced thalamic MRI volume has been observed in amnestic Mild Cognitive Impairment (MCI) 35 , LOAD 36 and presymptomatic presenilin 1 mutation carriers 37 . Similarly, Positron Emission Tomography (PET) studies have found APOE-ε4 state to accelerate longitudinal reductions in glucose metabolism in the thalamus and frontal, parietal, and posterior cingulate regions in MCI 38 . Reduced glucose metabolism in anterior and posterior cingulate cortices, retrosplenial, precuneus, parietal cortex, hippocampus and thalamus was also observed in cognitively healthy middle-aged APOE-ε4 carriers 39 , suggesting that metabolic tissue changes in regions beyond the hippocampus can already occur at asymptomatic stages 40 . While PET imaging is sensitive to metabolic changes and can identify amyloid-β and tau burden 41 , it is invasive and expensive and, therefore, difficult to scale up. Recent advances in non-invasive multi-parametric quantitative MRI (qMRI) methods can reveal subtle microstructural brain changes and promise to provide alternative imaging markers that may be sensitive to early risk-related changes. Up to now qMRI measurements have primarily been studied in LOAD patients and animal models, thus evidence with regards to the effects of risk factors in asymptomatic individuals is sparse.
To address this gap in the literature, we went beyond morphological analyses by employing multi-parametric qMRI to study the effects of APOE-ε4, Family History (FH) of dementia and obesity on cortical and subcortical grey matter in 165 asymptomatic individuals from the Cardiff Ageing and Risk of Dementia Study (CARDS) [42][43][44] (Table 1). More specifically we applied indices sensitive to neurite dispersion and density, free water, myelin and cell metabolism from Neurite Orientation Density and Dispersion Imaging (NODDI) 45 , quantitative magnetization transfer (qMT) [46][47][48][49] and T 1 -relaxometry 50 (Table 2).
NODDI fits a three-compartment biophysical tissue model to diffusion-weighted data acquired with a twoshell (b-values of 1200 s/mm 2 and 2400 s/mm 2 ) High Angular Resolution Diffusion Imaging (HARDI) 51 protocol to separate isotropic from intra-and extracellular diffusion compartments 45 . This allows the calculation of the isotropic signal fraction (ISOSF), an estimate of free water, and the intracellular signal fraction (ICSF), i.e. the fraction of the tissue comprised of neurites. In addition, NODDI yields the orientation dispersion index (ODI) that reflects the spatial configuration of neurite structures ( Table 2). Recent studies reported ICSF and ODI reductions in grey and white matter of patients with MCI, LOAD and young onset AD [52][53][54] . For instance, Fu et al. (2019) found decreased ICSF and ODI in the corpus callosum in MCI and LOAD patients, while Colgan et al. 55 reported positive correlations between ICSF and histological measurements of hyperphosphorylated tau protein in the hippocampus of rTg4510 mice.
The qMT method models the exchange rate between macromolecular protons and protons in surrounding free water when macromolecular protons are selectively saturated by a radiofrequency pulse with a frequency that is off-resonance for protons in free water [46][47][48][49] . This allows the quantification of a number of parameters including the macromolecular proton fraction (MPF) and the magnetization transfer exchange rate k f 49 . In combined neuroimaging and histology studies of Shiverer mice and puppies [56][57][58] , MPF has been shown to be highly sensitive to the myelin content in white matter such that MPF increases with the amount of myelin. MPF in the anterior hippocampus was also found to distinguish healthy controls from MCI and LOAD patients 59 . Furthermore, MCI Table 1. Summary of demographic, genetic, and lifestyle risk information of CARDS participants. APOE = Apolipoprotein-E based on DNA extraction and APOE genotyping of saliva samples using TaqMan genotyping of single nucleotide polymorphism (SNP) rs7412 and KASP genotyping of SNP rs429358. FH = Family History of a first degree relative affected by Alzheimer's or Lewy body disease or vascular dementia. MMSE = Mini Mental State Exam (maximum score = 30) 42 , NART-IQ = National Adult Reading Test-Intelligence Quotient 66 , PHQ-9 = Patient Health Questionnaire (maximum score = 27) 109 . WHR = Waistto-Hip-Ratio.

Sample size n 165
Age (in years) 55 www.nature.com/scientificreports/ and LOAD patients exhibit a reduced rate of magnetization transfer k f in grey and white matter 59-61 suggesting reduced cell metabolism 60 . Finally, indices from relaxometry imaging such as the longitudinal relaxation rate R 1 have been proposed as non-invasive biomarkers of LOAD 62 . R 1 values are influenced by microstructural characteristics such as tissue density, macromolecular, protein and lipid composition, and paramagnetic atoms. A number of patient and preclinical studies have reported increases in R 1 that may reflect LOAD pathology, although the precise mechanisms underpinning these changes remain unknown (see for review 62 ).
Here, we characterised age and risk-related differences in mean values of ICSF, ISOSF, ODI, MPF, k f and R 1 across cortical and subcortical grey matter regions that were segmented from T 1 -weighted images with the Free-Surfer image analysis suite (version 5.3) 63 . Microstructural changes were compared with differences in standard morphological metrics of cortical thickness and subcortical volumes. We expected to see risk effects in brain regions known to be early affected in LOAD including limbic regions of the hippocampus, parahippocampus, entorhinal cortex, posterior cingulate cortex as well as thalamus 2,4,34,64 . We hypothesised that APOE-ε4, a positive FH, and central obesity [measured with the Waist-Hip-Ratio (WHR)] would be associated with reduced ICSF, R 1 , MPF and k f as well as with increased ISOSF and ODI but with no differences in cortical thickness and/or subcortical volume. In addition, we expected to see the largest differences in those individuals at greatest risk, i.e. in obese APOE-ε4 carriers with a positive FH.

Results
Microstructural and morphological dependent variables were fitted to a general linear model in SPSS version 26 65 . All data were examined for outliers defined as above or below three times of the interquartile range (75th percentile value-25th percentile value). This led to an exclusion of 0.6% of the microstructural but no exclusions of the morphological data.
Separate multivariate analyses of covariance (MANCOVA) were carried out to test for the effects of APOE genotype (ε4 + , ε4-), FH (FH + , FH-) and WHR (WHR + , WHR-) on brain morphology (cortical thickness and subcortical volume measures) and on each of the microstructural indices (MPF, k f , R 1 , ISOSF, ICSF, ODI) across 68 cortical and 14 subcortical regions of interest, whilst controlling for age, sex, and IQ estimates from the revised National Adult Reading Test (NART-R) 66 . Significant omnibus effects were further investigated with post-hoc comparisons across all outcome measures. All first and post-hoc models were corrected for multiple comparisons with a False Discovery Rate (FDR) of 5% using the Benjamini-Hochberg procedure 67 (p BHadj ). As the aim of the study was to explore microstructural indices that could potentially provide novel biomarkers of dementia risk in future studies, a false positive rate of below 5% was regarded as an acceptable threshold to control for false positives while minimising the risk of missing any true risk-related microstructural differences. Information about effects sizes was provided with the partial eta squared index η p 2 for MANCOVA analyses, Cohen's d z for group comparisons and Pearson's r for correlational analyses. MPF   www.nature.com/scientificreports/ MPF post-hoc effects. APOE-ε4 carriers relative to non-carriers had lower MPF in the left thalamus (Table 3) (Fig. 1). Women had higher MPF than men in the left and right rostral middle frontal cortices, in the left superior temporal cortex and the right transverse temporal cortex ( ISOSF post-hoc effects. Ageing was associated with bilateral increases in ISOSF in medial regions including the cingulate, precuneus and cuneus cortices and in lateral regions including superior temporal, supramarginal, postcentral, pars opercularis and insula cortices. Age-related increases in ISOSF were also observed in left middle temporal and pars triangularis regions as well as in subcortical hippocampi, thalami, nuclei accumbens and right putamen (Table 4) (Fig. 3). Men relative to women had higher ISOSF in widespread frontal, temporal, parietal and cingulate cortices and in caudate nuclei, hippocampi, thalami and right nucleus accumbens (Table 4) (Fig. 2). In addition, NART-IQ correlated positively with ISOSF in the superior temporal sulci (left: r = 0.253, p BHadj = 0.008; right: r = 0.241, p BHadj = 0.006), left superior parietal (r = 0.227, p BHadj = 0.006), and right lingual (r = 0.182, p BHadj = 0.026) cortices (Table 4). After partialling out of age only correlations on the left hemisphere remained significant [superior parietal cortex [(r = 0.206, p BHadj = 0.048), superior temporal sulcus (r = 0.197, p BHadj = 0.032)] but those on the right did not [superior temporal sulcus (p BHadj = 0.053), lingual (p BHadj = 0.08)]. ODI post-hoc effects. Age-related increases in ODI were observed in left hippocampus, amygdala, caudate and right transverse temporal cortex (Table 5) (Fig. 3).

MANCOVAs of microstructural qMT metrics.
Post-hoc effects for the interaction between FH and WHR did not survive 5% FDR correction (Supplementary Table 2). Post-hoc effects. Ageing was associated with widespread thinning in bilateral frontal, temporal, and parietal cortical regions as well as with volume loss in subcortical structures, i.e. in the left hippocampus, left nucleus accumbens, bilateral thalami and putamen (Table 6) (Fig. 3). Women relative to men had larger volumes in left hippocampus, left nucleus accumbens, left putamen, right caudate and right pallidum. They also had larger cortical thickness in the right isthmus cingulate but lower cortical thickness in the left insula ( Exploring moderator effects of blood pressure and markers of inflammation. We then explored with two separate analyses of covariances whether controlling for differences in (i) systolic and diastolic blood pressure (BP) and (ii) inflammation-related measures of C-Reactive Protein (CRP), Interleukin-8 (IL-8) and leptin/adiponectin ratio (LAR) would account for the effect of APOE on left thalamus MPF. www.nature.com/scientificreports/   121    www.nature.com/scientificreports/

Discussion
Here, we investigated whether qMRI indices of apparent neurite density and dispersion, free water, myelin, and cell metabolism were sensitive to grey matter differences related to LOAD risk in cognitively healthy individuals. Such microstructural measurements hold the potential for novel imaging biomarkers to identify asymptomatic individuals at heightened risk of developing LOAD. As such they may provide non-invasive and cheaper alternatives to PET and cerebrospinal fluid (CSF)-based biomarkers, that are currently employed in clinical trials, in the future. The only significant difference between asymptomatic APOE-ε4 carriers relative to non-carriers was in the qMT measure MPF in the left thalamus with APOE-ε4 related reductions in MPF (Fig. 1). This effect was observed independently of age, sex, and verbal intelligence. Reduced MPF may arise from processes that lead to an increase in free water and/or a reduction in the macromolecular content of grey matter including changes in myelin, proteins, and and/or iron concentrations 68,69 . Such changes may be consistent with the presence of inflammatory processes leading to tissue swelling associated with glia activation 70 and/or with a deficit in cholesterol transport in APOE-ε4 carriers [70][71][72] . Consistent with this interpretation we observed that the effect of APOE genotype on left thalamus MPF was moderated by plasma markers of inflammation (CRP, IL-8, LAR). Furthermore, evidence suggests that APOE-ε4 carriage may increase susceptibility to inflammation 22,23 and that inflammatory processes contribute significantly to the pathogenesis of LOAD [73][74][75] .
Notably these APOE-ε4-related differences in MPF were only observed in the left thalamus but not in any other cortical or subcortical region. The limbic thalamic nuclei maintain dense reciprocal connections with the hippocampal formation and the retrosplenial cortex 76,77 , which, together with the fornix, mamillary bodies and posterior cingulate cortex, comprise the Papez circuit important for episodic memory function 78 . As outlined above it is increasingly recognised that the Papez circuit, including the anterior thalamus, can be affected early in LOAD 4 . Neurofibrillary accumulations are found in the anterodorsal thalamic nucleus at the same time as those in the hippocampus in LOAD brains 34 and neuroimaging studies have revealed reduced thalamic volume in both amnestic MCI 35 and LOAD 36 . Furthermore, studies into the effects of APOE in middle-aged asymptomatic adults found reduced glucose metabolism in the thalamus, hippocampus and cingulate cortex 39 as well as increased metabolism in bilateral thalami and superior temporal gyrus in amyloid-β positive APOE-ε4 carriers with a maternal history of LOAD 79 . Cacciaglia et al. 80 studied the effects of APOE on grey matter volume in over 500 middle-aged asymptomatic individuals and identified reduced hippocampus, caudate, precentral gyrus, and cerebellum volumes but increased volumes in the thalamus, superior frontal and middle occipital gyri in APOE-ε4 carriers. While it remains unknown why APOE-ε4 may be related to increased thalamic volume it was suggested that this could reflect brain swelling associated with glial activation in response to larger amyloid-β www.nature.com/scientificreports/ burden 81 . As mentioned above, the here observed pattern of APOE-ε4-related reductions in MPF in the left thalamus is consistent with this interpretation 56,82 . One other study investigated the impact of APOE-ε4 on qMT white matter metrics in young adults and did not find any differences 83 . This suggests that such risk-related glial dysfunction may accumulate with age and may only become apparent from midlife onwards.
The question arises why we did not observe any risk-related effects in brain regions that have previously been reported to be affected by LOAD risk factors 10,84,85 . Reports with regards to the impact of APOE-ε4 on grey matter structures in healthy young and middle-aged adults have been mixed 10,84 , with some studies reporting no changes in hippocampal grey matter volume in APOE-ε4 carriers 31,86 . Studies assessing the impact of APOE-ε4 on tissue microstructure have primarily focused on diffusion tensor imaging (DTI) of white matter. While some reported widespread white matter differences in DTI measures 83,87,88 , this has not been replicated in all studies 30,89 . These discrepancies may arise due to DTI indices not being sufficiently sensitive and/or specific to detect early risk-related tissue abnormalities 90 . Direct comparisons between DTI and NODDI indices revealed that although fractional anisotropy (FA) was sensitive to white matter differences between healthy controls and patients with metabolic disease, FA was less anatomically specific and did not identify all brain regions that were captured by ICSF and ODI 91 . Thus we employed NODDI and qMT measurements to study risk effects on grey matter here and on white matter in a previous CARDS analysis 92 . In the previously published white matter analysis 92 we did not observe any main effects of risk but found that individuals with the highest genetic risk (obese FH + and APOE-ε4) exhibited obesity-related reductions in MPF and ICSF in the right parahippocampal cingulum.
Taken together, our previous and here reported findings demonstrate that MPF from qMT can identify risk-related microstructural differences in limbic grey and white matter that were not apparent in conventional volumetric or cortical thickness measurements. We propose that these differences may reflect subtle changes related to neuroglia activation and that limbic structures including the thalamus are particularly susceptible to adverse effects of APOE-ε4 on glia cells. Inconsistencies in previous studies may have arisen from standard morphological and DTI measurements not being sensitive and/or specific enough to detect such glia-related changes.
It is important to note that while we did not find any risk-related effects on brain morphology we did replicate the well-established pattern of widespread age-related thinning in frontal, temporal and parietal regions 93 as well as volume loss in subcortical structures including the hippocampi and thalami (Fig. 3). The subcortical volume loss was accompanied by age-related increases in ISOSF in bilateral hippocampi and thalami but effects on cortical regions were more localised: increased ISOSF was apparent along medial regions of the cingulate and parietal cortices including the precuneus as well as in superior temporal and lateral and orbito prefrontal cortices. Age-related increases in ISOSF have been previously observed 94 and most likely reflect lost tissue being replaced by CSF. Consistent with a previous study 95 we also observed a positive correlation between age and ODI, an estimate of neurite dispersion, in the hippocampus and the left caudate and amygdala. In contrast to Nazari et al. 95 however, we did not find any effects in cortical regions, while they reported reduced ODI with age in fronto-parietal regions. These opposing patterns in cortical and subcortical regions may reflect age-related reductions of neocortical dendritic spine density 96 with accompanying compensatory increases in the dendritic extent of dentate gyrus granular cells 97,98 . Similar age-related increases in the dendritic tree have also been reported in the basolateral nucleus of the amygdala of rats 99 .
Furthermore, we observed positive correlations between ISOSF and NART-IQ in superior temporal, parietal and lingual cortices that were partly driven by age. NART requires the reading of irregularly pronounced words and older relative to younger adults tended to perform better in the NART. However, positive albeit weak correlations between NART-IQ and ISOSF remained for the left superior temporal sulcus and left superior parietal cortex. Developmental imaging studies have revealed cortical thinning during adolescence 100 that may be due to increased myelination 100 or synaptic pruning and dendritic arborization 101,102 . It may therefore be possible that childhood developmental differences in cortical maturation as well as in education may have contributed to this effect. For instance, childhood cognitive abilities have been found to account for relationships between cognitive performance and brain cortical thickness decades later in older adults from the Lothian birth cohort 103 .
Consistent with previous reports 104 we did not observe widespread sex-differences in brain morphology measurements with the exception of larger volumes in the left hippocampus in women than men 105 . However, qMRI indices revealed the following pattern: Women compared to men, had lower ISOSF in widespread cortical and subcortical regions and larger MPF in frontal and temporal regions. Previously we also reported higher MPF and lower ISOSF for white matter in women than men 44 . Overall this pattern of sex differences suggests higher cortical myelination and lower free water signal in women as they tended to be overall in better health i.e. were less obese, had lower systolic BP, and reported drinking less alcohol than men 44 . All of these factors may have contributed to women showing "healthier" grey and white matter in the CARDS cohort.
Finally, some study limitations need to be considered. First of all, CARDS is a cross-sectional study that cannot answer whether the observed APOE effects on left thalamus MPF are predictive of accelerated development of LOAD pathology, cognitive, or neuronal decline. Future prospective longitudinal studies are required to address this question. We also propose that our findings require replication in larger samples that can control for possible interactions between APOE and other LOAD risk genes such as variants of TREM2 and polygenic risk hazards as the number of participants in the CARDS study was too small to do so. It is also worth mentioning that other qMRI measurements, that were not included in the current study, may prove helpful in characterising risk effects on the brain. Notably quantitative T 2 and T 2 * measurements have been proposed to be sensitive to neurodegenerative processes. For instance, prolonged T 2 relaxometry has been reported in the hippocampus of LOAD patients 106 and has been proposed to increase the sensitivity and specificity of MCI and LOAD detection 107 . Finally, it should be noted that we only studied the thalamus as a whole structure while neuropathological evidence suggests a specific vulnerability of the anterodorsal thalamic nucleus to LOAD pathology. Future studies may investigate www.nature.com/scientificreports/ risk-related effects on specific subthalamic nuclei, which was beyond the scope of the current study as we were focusing on risk effects across the whole brain. In summary, we have shown APOE-ε4 related reductions in the qMT measure MPF in the left thalamus that were moderated by peripheral markers of inflammation. This effect occurred independently of age, sex and NART-IQ and was not observed in morphological or microstructural indices from diffusion-weighted imaging. In addition, the effect was specific to the left thalamus and was not present in other cortical and subcortical grey matter regions. We propose that MPF reductions may reflect the effects of glia-mediated inflammatory and demyelination processes in APOE-ε4 carriers. As such qMT measurements hold the potential for non-invasive and cheaper biomarker alternatives to PET, that may aid our understanding of the pathological processes leading to LOAD. In addition, qMT may help with the identification of asymptomatic individuals at heightened risk of LOAD for stratification into clinical trials for future preventative therapeutics.

Materials and methods
The Cardiff Ageing and Risk of Dementia Study (CARDS) has been described previously including a detailed description of the participant sample 43,92 , assessment of genetic and metabolic risk factors 44,92 and the acquisition and processing of the MRI data 43,44,92,108 . Here we provide a brief summary of the most important points. CARDS received ethical approval from the School of Psychology Research Ethics Committee at Cardiff University (EC.14.09.09.3843R2) and all participants provided written informed consent in accordance with the Declaration of Helsinki. All research methods were performed in line with Cardiff University's Research Integrity and Governance Code of Practice and relevant data protection regulations.
Participants. The CARDS cohort comprised 166 community-dwelling individuals between the age of 38 and 71 years who underwent cognitive and health assessment as well as MRI scanning (Table 1). Exclusion criteria were a history of neurological and/or psychiatric disease, head injury, drug/alcohol dependency, high risk cardio-embolic source, large-vessel disease or MRI incompatibility due to pacemaker, stents or other surgical implants. As a group, participants intellectual functioning was above average as assessed with the National Adult Reading Test (NART) 66 . All but one participant scored > 26 on the Mini Mental State Exam (MMSE) 42 thus the remaining 165 participants were classified as cognitively healthy. Eight participants scored ≥ 10 in the Patient Health Questionnaire (PHQ)-9 109 , suggesting moderate levels of depression but no participant was severely depressed.
Assessment of risk factors. Saliva samples were collected with the Genotek Oragene-DNA kit (OG-500) and APOE genotypes ε2, ε3, and ε4 were determined with TaqMan genotyping of single nucleotide polymorphism (SNP) rs7412 and KASP genotyping of SNP rs429358. Participants self-reported their family history of dementia, i.e., whether a first-grade relative was affected by Alzheimer's disease, vascular dementia or any other type of dementia.
Central obesity was assessed from the waist-hip ratio (WHR) 44 with abdominal obesity defined as a WHR ≥ 0.9 for males and ≥ 0.85 for females. Resting systolic and diastolic blood pressure (BP) readings were taken with a digital blood pressure monitor (Model UA-631; A&D Medical, Tokyo, Japan) and the means of three readings were calculated. Participants self-reported other metabolic risk factors, including diabetes mellitus, high levels of blood cholesterol controlled with statin medication, history of smoking, and weekly alcohol intake. There were only few diabetics, smokers, and individuals on statins and, hence, these variables were not included in the analyses.
Blood plasma analysis. As previously reported 44,92 , venous blood samples were drawn into 9 ml heparin coated plasma tubes after 12 h overnight fasting and were centrifuged for 10 min at 2000 × g within 1 h from blood collection. Plasma samples were then transferred into 0.5 ml polypropylene microtubes and stored in a freezer at − 80 °C. Circulating levels of high-sensitivity C-Reactive Protein (CRP) in mg/dL were assayed using a human CRP Quantikine enzyme-linked immunosorbent assay (ELISA) kit (R & D Systems, Minneapolis, USA). Six individuals had a CRP value > 10 mg/ml indicative of acute infection and were, therefore, excluded from the statistical analyses testing for moderating effects of inflammation. Leptin concentrations in pg/ml were determined with the DRP300 Quantikine ELISA kit (R & D Systems) and adiponectin in ng/ml with the human total adiponectin/Acrp30 Quantitkine ELISA kit (R & D Systems). Leptin/adiponectin ratios for each participant were calculated. Interleukin IL-8 levels in pg/mL were determined using a high sensitivity CXCL8/ INTER-LEUKIN-8 Quantikine ELISA kit (R & D Systems). Determination of interleukin-1β, interleukin-6 and Tumor Necrosis Factor α (TNFα) levels were trialled with high-sensitivity Quantikine ELISA kits but did not result in reliable measurements consistently above the level of detection for each assay.
HARDI and qMT data processing. As described in 43,44,92,108 , the dual-shell HARDI data were split and b = 1200 and 2400 s/mm 2 data were corrected separately for distortions induced by the diffusion-weighted gradients and motion artifacts with appropriate reorientation of the encoding vectors 112 in ExploreDTI (Version 4.8.3) 113 . EPI-induced geometrical distortions were corrected by warping the diffusion-weighted image volumes to the T 1 -weighted anatomical images 114 . After pre-processing, the NODDI model 45 was fitted to the HARDI data with the fast, linear model fitting algorithms of the Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework 115 to gain ISOSF, ICSF, and ODI maps.
Using Elastix 116 , MT-weighted GRE volumes were co-registered to the MT-volume with the most contrast using a rigid body (6 degrees of freedom) registration to correct for inter-scan motion. Data from the 11 MT-weighted GRE images and T 1 -maps were fitted by a two-pool model using the Ramani pulsed-MT approximation 117 . This approximation provided MPF and k f maps. To remove voxels with noise-only data, MPF maps were thresholded to an upper intensity limit of 0.3 and k f maps to an upper limit of 3.0 using the fslmaths imaging calculator from the Functional Magnetic Resonance Imaging of the Brain (FMRIB) library (version 6).
All image modality maps were spatially aligned to the T 1 -weighted anatomical volume as reference image with linear affine registration (12 degrees of freedom) in within-subject space using FMRIB's Linear Image Registration Tool (FLIRT) 118,119 . Cortical and subcortical grey matter region segmentation. Grey matter cortical and subcortical regions were automatically segmented from T 1 -weighted images with the Freesurfer image analysis suite (version 5.3), which is documented online (https ://surfe r.nmr.mgh.harva rd.edu/) 64 . The images were processed by running the "recon-all" script using the default analysis settings. In brief, the images were registered to the Montreal Neurological Institute standard space and intensity normalization was performed. This was followed by automatic skull stripping to remove extracerebral structures, the cerebellum and the brain stem, followed by segmentation into grey matter, white matter and CSF and separation of the hemispheres. Pial surfaces were obtained by tessellating the grey and white matter boundary and by surface deformation following intensity gradients for optimal placement of grey and white matter and grey matter and CSF boundaries 120 . Surface inflation and registration to a spherical atlas were then performed and the cerebral cortex was parcellated into 34 regions per hemisphere based on gyral and sulcal structures following the Desikan-Killiany atlas 121 . Cortical thickness measurements were estimated as the average shortest distance between the pial surface and the white matter boundary 122 . For each hemisphere, seven deep grey matter structures (hippocampus, amygdala, thalamus, caudate, putamen, pallidum, and nucleus accumbens) were automatically parcellated using a probabilistic atlas so that average volumetric measurements could be determined 123,124 . Mean intracranial volume fractions (ICV) were extracted for each brain as estimates of individual differences in head sizes and all volumetric measurements were adjusted for ICV by dividing each participant's subcortical volume by their ICV.
Finally, the mean values of all microstructural indices were extracted from each participants' cortical and subcortical region of interests. Mean measurements were taken in each participants' native space. This was done by first converting each participants' cortical and subcortical masks from the FreeSurfer Massachusetts General Hospital volume file format (MGZ) into the Neuroimaging Informations Technology Initiative (NIfTI) analyze-style data format and then uploading the microstructural maps onto each region of interest mask using the fslmaths command from the FMRIB library. Mean values of each index for each mask were then extracted using the FMRIB fslstats command. NODDI and qMT indices of ISOSF, ICSF, ODI, MPF and k f , could not be extracted from bilateral caudal middle frontal, left isthmus cingulate and left pericalcarine regions and R 1 could not be extracted from the right postcentral region.