Diffusion models reveal white matter microstructural changes with ageing, pathology and cognition.

White matter microstructure undergoes progressive changes during the lifespan, but the neurobiological underpinnings related to ageing and disease remains unclear. We used an advanced diffusion MRI, Neurite Orientation Dispersion and Density Imaging, to investigate the microstructural alterations due to demographics, common age-related pathological processes (amyloid, tau and white matter hyperintensities) and cognition. We also compared Neurite Orientation Dispersion and Density Imaging findings to the older Diffusion Tensor Imaging model-based findings. Three hundred and twenty-eight participants (264 cognitively unimpaired, 57 mild cognitive impairment and 7 dementia with a mean age of 68.3 ± 13.1 years) from the Mayo Clinic Study of Aging with multi-shell diffusion imaging, fluid attenuated inversion recovery MRI as well as amyloid and tau PET scans were included in this study. White matter tract level diffusion measures were calculated from Diffusion Tensor Imaging and Neurite Orientation Dispersion and Density Imaging. Pearson correlation and multiple linear regression analyses were performed with diffusion measures as the outcome and age, sex, education/occupation, white matter hyperintensities, amyloid and tau as predictors. Analyses were also performed with each diffusion MRI measure as a predictor of cognitive outcomes. Age and white matter hyperintensities were the strongest predictors of all white matter diffusion measures with low associations with amyloid and tau. However, neurite density decrease from Neurite Orientation Dispersion and Density Imaging was observed with amyloidosis specifically in the temporal lobes. White matter integrity (mean diffusivity and free water) in the corpus callosum showed the greatest associations with cognitive measures. All diffusion measures provided information about white matter ageing and white matter changes due to age-related pathological processes and were associated with cognition. Neurite orientation dispersion and density imaging and diffusion tensor imaging are two different diffusion models that provide distinct information about variation in white matter microstructural integrity. Neurite Orientation Dispersion and Density Imaging provides additional information about synaptic density, organization and free water content which may aid in providing mechanistic insights into disease progression.

Abbreviated summary

Raghavan et al. report that advanced diffusion MRI can provide specific information about synaptic density, organization, and free water content of white matter microstructure. Age and white matter hyperintensities but not Alzheimer’s disease are key drivers of these changes. Corpus callosum white matter health is fundamental for better cognitive performance.

Introduction

The white matter (WM) architecture of the human brain undergoes substantial changes across the life span. There is clear evidence for the association between WM changes and age as well as neuropathological processes that will lead to cognitive decline. Diffusion MRI is a versatile method that allows us to study these WM microstructural details. Previous findings based on diffusion tensor imaging (DTI) revealed reduced fractional anisotropy (FA) and increased mean diffusivity (MD) in association with amyloid deposition, a hallmark of Alzheimer’s disease, and cerebrovascular disease.

Despite its sensitivity, the clinical utility of DTI is constrained by its inherent limitation in specificity of identifying the different diffusion environments that exist within most individual voxels. Characterizing the different water pools within a voxel with a single diffusion tensor is well known to be problematic in crossing fibre regions of WM, and also confounds the macroscopically isotropic diffusion of grey matter (GM), with that of CSF. The growing availability of multiband excitation allows the acquisition of roughly three times as much data in the same time as a standard DTI scan, making multiple b value (diffusion weighting) shells clinically practical. Distributing the diffusion samples over >2 b values allows the use of more sophisticated and biologically plausible models to characterize the general properties of the microstructural environments inside the axons, between them, and in the extracellular water. In addition, these models can ideally handle the ‘crossing fibre problem’ better than DTI. Neurite orientation dispersion and density imaging (NODDI) is an advanced dMRI technique that uses the additional degrees of freedom from multi-shell data to probe the microstructural complexity of neurites (dendrites and axons), separately from CSF and to a large degree also separately from each other. This biophysical modelling method divides water diffusion in the brain into three microstructural compartments: intracellular space through the Neurite Density Index (NDI), which measures the signal fraction that is due to axons and dendrites; Orientation Dispersion Index (ODI), which measures angular variation or dispersion of the neurites; and the Isotropic Volume Fraction (ISOVF), which measures free water (FW) fraction. More recently, a number of studies have demonstrated the efficiency of NODDI to provide finer granularity, in comparison to DTI metrics, to decipher the intra and extracellular microscopic features of age- and sex-specific diffusion trajectories., In addition, NODDI has been found to be useful for the early detection of neurodegenerative changes and its association with cognitive deficits.,

Recent findings have suggested that amyloid affects DTI measures, and also the greater effect of cerebrovascular disease on diffusion alteration than Alzheimer’s disease in memory clinic patients. It is also well known that WM plays an important role in normal cognition and age-related cognitive decline., However, the efficiency of NODDI models over DTI models to detect Alzheimer’s disease and cerebrovascular disease pathologies, and their contribution to cognitive performance, in population-based studies remains unclear. Given the detailed quantification of biological processes by NODDI models,,, we hypothesized that NODDI measures would provide more sensitive features of microstructural brain changes than conventional FA and MD, and would be more sensitive in capturing disease related processes. The overall goal of the study was to identify the relationships between demographics (age, sex and education/occupation), neuroimaging measures of Alzheimer’s disease and cerebrovascular disease, and cognition with diffusion MRI (NODDI and conventional DTI) in participants from Mayo Clinic Study of Aging (MCSA).

Materials and methods

Selection of participants

We identified 328 participants consisting of 264 cognitively unimpaired (CU), 57 mild cognitive impairment (MCI) and 7 dementia from the MCSA, an epidemiological cohort designed to investigate the prevalence, incidence and risk factors for MCI and dementia among the residents of Olmsted County, Minnesota. The Rochester Epidemiology Project (REP) medical records-linkage system, was used to enumerate the MCSA sample population. The MCI and dementia participants were diagnosed based on the previously published consensus criteria. Our inclusion criteria were participants who had multi-shell diffusion data, fluid attenuated inversion recovery (FLAIR)-MRI, amyloid and tau PET scans. The clinical diagnosis was made at the time of MRI assessment and almost all clinical and imaging visits were within 60–70 days (median was 63 days with a range of 0–124 days). The amyloid negative/positive (A−/A+) and tau negative/positive (T−/T+) proportions in the sample were CU (A − T − n = 171, A − T + n = 23, A + T − n = 41 and A + T + n = 29), MCI (A − T − n = 16, A − T + n = 3, A + T − n = 12 and A + T + n = 26) and dementia (A + T − n = 1 and A + T + n = 6).

Standard protocol approvals, registrations and patient consents: The study was approved by the Mayo Clinic institutional review board and written informed consent was obtained from all participants or their surrogates.

Imaging

MRI acquisition and processing

All participants underwent a 3 T head MRI protocol on one of two 3 T Siemens Prisma scanners running VE11 software with 64 channel receiver head coils. The protocol included a magnetization prepared rapid gradient echo (MPRAGE) sequence (TR/TE/TI = 2300/3.14/945 ms, flip angle 9°, 0.8 mm isotropic resolution), and a diffusion scan using the product VE11 Simultaneous Multi-Slice acceleration with adaptive coil combination. For the diffusion scan the field of view was 232 mm in X and Y and 162 mm in the Z direction, with 2.0 mm isotropic voxels. The echo and repetition times were 71 and 3400 ms respectively. Data consisted of 127 volumes with 13 non-diffusion-weighted images (b = 0 s/mm2), and 114 diffusion-encoding gradient directions (6 b = 500, 48 b = 1000 and 60 b = 2000 s/mm2), evenly spread over the entire spherical shells using an electrostatic repulsion model, and interspersed in time to minimize gradient heating.

The diffusion data were preprocessed using the in-house developed pipeline. After visual inspection, an intracranial mask was made for the diffusion MRI scan and the noise in the raw diffusion images was estimated and removed using random matrix theory. Then FSL’s eddy_cuda was used to correct for head motion and eddy current distortion, followed by the correction of Gibbs ringing and Rician bias. Diffusion tensors were fitted for both the multi-shell and extracted b = 1000 data using a non-linear least-squares fitting algorithm implemented in dipy, from which FA and MD images were generated. The NODDI model was fit by the Accelerated Microstructure Imaging via Convex Optimization (AMICO) implementation in Python. FA, MD, NDI, ODI and ISOVF maps generated from a cognitively unimpaired subject are shown in Fig. 1A.

Figure 1

DTI and NODDI maps from a participant and white matter tracts of interest from the JHU atlas. (A) FA, MD, NDI, ODI and ISOVF maps generated from a 67-year-old cognitively unimpaired female participant. (B) White matter tracts of interest from JHU atlas. FA, fractional anisotropy; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index.

Amyloid and tau assessment from PET scans

The acquisition and processing were described previously. From amyloid PET scans, a global amyloid load measure for each subject [standardized uptake value ratio (SUVR)] was computed by calculating the median uptake in the prefrontal, orbitofrontal, parietal, temporal, anterior cingulate and posterior cingulate/precuneus regions of interest (ROIs) normalized by the median amyloid PET uptake in the cerebellar crus grey matter. From tau PET scans, a composite ratio for each subject was computed by calculating median tau PET uptake in the entorhinal, amygdala, parahippocampal, fusiform, inferior temporal and middle temporal ROIs normalized by the median tau PET uptake in the cerebellar crus grey matter.

WMH assessment from FLAIR scans

The 3D MPRAGE and 3D T2-weighted FLAIR images were used to calculate WMH volume via a fully automated algorithm, updated from a previously described in-house semi-automated method. Briefly, 3D FLAIR images were preprocessed for intensity inhomogeneity correction and de-noising using a non-local means filter. Then, WMH were segmented based on location (spatial priors), intensity relative to the global distribution of GM intensity values, and intensity relative to the local neighbourhood of WM voxels. False-positive WMH segmentations were reduced by applying a white matter mask derived from the 3D MPRAGE segmentation, and by removing isolated single-voxel detections.

Cognitive performance measures

All MCSA participants underwent a detailed neuropsychological test battery that consisted of 9 tests covering 4 cognitive subdomains., The present study utilized a global cognitive z-score that was derived as the z-transformation of the average of all nine tests across the 4 cognitive domains (memory, language, attention/executive and visuo-spatial function). Individual and compound scores from Trail Making Test (Trails) A and B (time to complete the test) were used as a sensitive test for processing speed. The raw scores were transformed into z-scores and averaged to create a composite score.

Image analysis

ROI-based analysis

We performed an ROI analysis in ten WM tracts which were selected based on literature suggesting their association with cognition.,,,: commissural fibres: genu (GCC), body (BCC), and splenium (SCC) of corpus callosum, and fornix (FX); association fibres: cingulum (CGC), parahippocampal cingulum (CGH), superior longitudinal fasciculus (SLF), inferior fronto-occipital fasciculus; and other relevant tracts: inferior temporal WM (ITWM) and anterior limb of internal capsule (Fig. 1B). The median values of FA, MD, NDI, ODI and ISOVF were computed in these tracts by non-linearly registering an in-house modified version of the JHU ‘Eve’ WM atlas to each subject’s image using Advanced Normalization Tools–Symmetric Normalization (ANTS-SyN). In this analysis, we excluded the cuneus, precuneus, fusiform and lingual WM regions since they are too small for reliable registration. The median values of bilateral regions were then averaged, weighting by region size, to produce a single measure for each bilateral structure.

Voxel-based analysis of diffusion metrics

Diffusion images were analysed using an in house developed voxel-based analysis (VBA) pipeline for SPM12 in MATLAB to identify the global brain changes in association with demographics and disease pathologies. Briefly, each subjects FA, MD, NDI, ODI and ISOVF maps were non-linearly registered to a custom-made study-specific template using ANTs-SyN. To reduce partial volume effects and understand the regional results based on tissue classes, additional mask images were made using GM, WM and GM+WM masks from each subject’s segmented T1-weighted image. The masks were registered to the study template using an ANTS-calculated warp from the subject’s T1-weighted image to a T1-like target synthesized from the FA and MD templates. The GM, WM and GM+WM masks were thresholded to include voxels with respective tissue type fractions >0.5. Each of the normalized diffusion images was then smoothed with an 8-mm FWHM isotropic Gaussian kernel and analysed per-voxel within each tissue-class mask, using SPM12.

Statistical analyses

Characteristics of the participants were summarized as mean (standard deviation) for the continuous variables and count (%) for the categorical variables. WMH was presented and analysed as a percentage of total intracranial volume (TIV). The distributions of WMH and amyloid were skewed, and hence log transformed to obtain a more normal distribution. To describe the relationships between NODDI and DTI parameters, we performed a series of unadjusted Pearson correlation analyses associating FA with NDI, FA with ODI, MD with NDI, MD with ODI, and MD with ISOVF across subjects within each WM tract. We also used unadjusted Pearson correlation analyses to describe associations between demographics (age, sex and education/occupation), cerebrovascular disease (WMH) biomarkers, Alzheimer’s disease (amyloid and tau) biomarkers, and ROI-based diffusion (FA, MD, NDI, ODI and ISOVF) measures (corrplot package 0.84).

To assess the contributions of cerebrovascular disease (WMH) and Alzheimer’s disease (amyloid and tau) biomarkers on the WM integrity changes, we fit multiple linear regression models with each ROI diffusion measure as the outcome variable, and with age, sex, education/occupation scores, WMH, amyloid and tau as predictor variables. All the imaging variables were standardized.

We also repeated the above analyses using voxel-wise multiple regression analyses on the smoothed DTI and NODDI images with age, sex, education/occupation, WMH, amyloid and tau as predictor variables. The generated SPM-T maps were corrected for multiple comparisons using family-wise error (FWE) with PFWE <0.05. The voxel level analyses also helped confirm ROI level analyses and provide insights into subtle associations missed by ROI analyses.

Finally, we estimated the association of global cognition with each diffusion variable after adjusting for age, sex, education/occupation, cycle number (the number of times the cognitive battery was administered to each specific subject to adjust for practice effects), and amyloid and tau PET. We repeated the analyses for subdomain scores (memory, attention, language and visuospatial) and processing speed (Trail A, Trail B, composite score) with regional WM microstructural integrity measures. We computed partial Pearson correlations with 95% confidence intervals and report the beta coefficients from the multiple regression analyses.

Data availability

Data used in this study are available upon reasonable request via MCSA/ADRC data sharing website.

Results

The characteristics of the participants are summarized in Table 1. The mean (standard deviation) age was 68.3 (13.1) years, 52% were men, 30% were APOE4 carriers, 35% were amyloid positive and 27% were tau positive. Cognitively unimpaired individuals comprised 80% of this sample.

Table 1

Characteristics table with the mean (SD) listed for the continuous variables and count (%) for the categorical variables

	CU n = 264	MCI n = 57	Dementia n = 7	P-value
Male, n (%)	135 (51%)	30 (53%)	4 (57%)	0.94
Age, years	65.9 (12.8)	77.2 (9.8)	84.7 (4.2)	<0.001
APOE4 carrier, n (%)	63 (27%)	21 (40%)	5 (71%)	0.010
Education/occupation	12.9 (2.2)	11.9 (2.8)	12.1 (2.3)	0.013
CMC	1.5 (1.4)	2.6 (1.8)	3.0 (1.0)	<0.001
MMSE	29.0 (1.1)	26.1 (1.9)	20.3 (3.9)	<0.001
zGlobal	0.63 (0.99)	−1.80 (1.17)	−3.05 (0.49)	<0.001
zMemory	0.64 (0.96)	−1.85 (0.92)	−3.04 (0.15)	<0.001
zAttention	0.41 (1.02)	−1.50 (1.57)	−3.24 (0.78)	<0.001
zLanguage	0.35 (0.99)	−1.48 (1.38)	−2.96 (1.32)	<0.001
zVisual-spatial	0.56 (0.99)	−1.05 (1.23)	−1.87 (1.36)	<0.001
Trails A	32.0 (11.7)	51.4 (24.9)	91.5 (51.3)	<0.001
Trails B	75.6 (40.0)	163.5 (96.2)	262.3 (93.0)	<0.001
Amyloid, SUVR	1.49 (0.34)	1.95 (0.62)	3.04 (0.41)	<0.001
Amyloid positive, n (%)	70 (27%)	38 (67%)	7 (100%)	<0.001
Tau, SUVR	1.18 (0.11)	1.29 (0.21)	1.54 (0.30)	<0.001
Tau positive, n (%)	52 (20%)	29 (51%)	6 (86%)	<0.001
WMH percentage	0.66 (0.88)	1.22 (0.93)	1.98 (1.57)	<0.001

The P-values reported are from an ANOVA for continuous measurements and a chi-squared for categorical variables.

CU, cognitively unimpaired; CMC, cardiovascular and metabolic conditions; MCI, mild cognitive impairment; MMSE, mini mental state examination; SUVR, standardized uptake value ratio; WMH—white matter hyperintensity.

Association between DTI and NODDI metrics in different WM tracts

Pearson correlations between FA, MD, NDI, ODI and ISOVF are shown in Fig. 2. Correlations within the same regions between measures: across the WM tracts, MD and NDI showed the strongest association with each other (r ≤ −0.676) except for in the fornix. In contrast, a modest association was observed between MD and ODI in half of the regions (fornix and association tracts). FA and NDI were associated modestly in most of the WM tracts, while FA and ODI (indicators of dispersion) had strong associations in the association tracts, anterior limb of internal capsule and inferior temporal WM (r ≤ −0.51). Similarly, MD was associated strongly with ISOVF in the corpus callosum.

Figure 2

Correlation matrix. Association between diffusion measures (FA and MD, FA and NDI, FA and ODI, MD and NDI, MD and ODI, MD and ISOVF). Colour legend indicates the range of correlations, the size of the circle indicates the strength of the correlation, and the symbol ‘X’ indicates the non-significant P-value. FA, fractional anisotropy; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index.

Associations with demographics and biomarkers of cerebrovascular disease and Alzheimer’s disease

Univariate associations: the univariate associations using unadjusted Pearson correlations between tract measures and age, sex, education/occupation, WMH, amyloid and tau are shown in Fig. 3. This figure highlights three broad aspects of the data: older age was significantly associated with lower FA, lower NDI, higher MD and higher ISOVF; the sex and education/occupation scores had either modest or no association with diffusion measures; and WMH showed the strongest association with all diffusion metrics in major WM tracts. We also observed associations across pairs of tracts with each DTI and NODDI measure (bottom of each triangle in Fig. 3). One can observe greater variability in the FA, ODI and ISOVF correlations across the tracts but MD and NDI appear to be correlated across all the tracts.

Figure 3

Correlation matrix. Association between demographics (age, sex, education/occupation) or white matter hyperintensity (WMH) or amyloid or tau and diffusion measures. Edu.occ represents education/occupation. Colour legend indicates the range of correlations, the size of the circle indicates the strength of the correlation, and the symbol “X” indicates the non-significant P-value. FA, fractional anisotropy; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index.

Multiple regression models with focus on disease pathologies

The regression models with standardized disease pathologies (WMH, amyloid and tau) as predictors and standardized WM integrity measures of FA, MD, NDI, ODI and ISOVF as outcomes are shown in Fig. 4 and Supplementary Table 1. Associations with WMH: Across all models, WMH had the strongest associations with all dMRI metrics from all tracts. Higher WMH (a surrogate of cerebrovascular disease) was significantly associated with lower FA, higher MD, lower NDI and higher ISOVF. Splenium was the only region where WMH showed a statistically significant association with ODI.

Figure 4

Association of diffusion metrics with white matter hyperintensity (WMH), amyloid and tau after controlling for age, sex and education/occupation. Different symbols below are used for each of the primary predictors. FA, fractional anisotropy; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index.

Associations with amyloid and tau

Higher amyloid was significantly associated with higher MD in parahippocampal cingulum (P = 0.026). Higher amyloid was associated with lower NDI in the same region, but the P-value was 0.053. In addition, higher MD was significantly associated with greater tau in inferior temporal WM (P = 0.014).

Voxel level associations for confirmation of ROI analyses

Similar to the ROI analysis, the voxel-wise analyses found the strongest associations for age and WMH with all diffusion metrics as displayed in Fig. 5. Modest associations were found with amyloid for both DTI and NODDI in the medial temporal lobe regions, specifically at the grey and WM junctions (Supplementary Fig. 1A). The extent and strength of tau associations with dMRI measures was minimal (Supplementary Fig. 1B).

Figure 5

Association between age or white matter hyper intensity (WMH) and diffusion measures. Significance level set at P < 0.05, FWE corrected with an extend threshold K = 100. (+ve and –ve represents the kind of association between variables). FA, fractional anisotropy; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index.

Association of diffusion measures with cognition

Association results from multiple linear regression models of the global cognition and cognitive subdomain z-scores with DTI and NODDI metrics after controlling for age, sex, education/occupation, cycle visit, amyloid and tau are shown in Fig. 6 for the corpus callosum tracts (where the correlations were highest). The regression coefficients for all tracts with global cognition are shown in Supplementary Table 2. Corpus callosum generally had the most significant findings except for analyses of NDI and ODI where superior longitudinal fasciculus and cingulum respectively had the greatest impact. The associations between subdomain scores and diffusion metrics had a similar pattern to that of global cognition with stronger associations with attention. Further analyses revealed significant associations between diffusion metrics and speed scores (Supplementary Fig. 2). As expected, the strongest associations were observed for corpus callosum fibres with MD and ISOVF and these changes were greatest with Trail B as an outcome.

Figure 6

Association of diffusion metrics with cognition after controlling for age, sex and education/occupation, cycle visit, amyloid and tau. Different symbols below are used for each of the diffusion measures. BCC, body of corpus callosum; FA, fractional anisotropy; GCC, genu of corpus callosum; ISOVF, isotropic volume fraction; MD, mean diffusivity; NDI, neurite density index; ODI, orientation dispersion index; SCC, splenium of corpus callosum.

Association of diffusion measures with disease pathologies and cognition in non-demented participants

We also performed sensitivity analyses after excluding dementia participants. As described above, we evaluated (i) the contribution of WMH, amyloid and tau on WM changes (after adjusting for age, sex, education/occupation) and (ii) association of global cognition, subdomain and processing speed scores with corpus callosum WM measures (after adjusting for age, sex, education/occupation, and cycle number and amyloid and tau PET). There were no significant differences observed in these sensitivity analyses as shown in the supplemental material (Supplementary Figs 3–5).

Discussion

We investigated the performance of DTI and NODDI models in capturing the microstructural brain changes associated with demographics and pathological processes and their association with cognition in 328 MCSA subjects aged 33–98 years. The major findings of the study were as follows: (i) NODDI and DTI are two different biophysical models that provide distinct information about variation in WM health. There was complementary information such that only MD and NDI had the strongest correlations with each other across the tracts; (ii) Age and WMH had the strongest associations with DTI and NODDI measures among the measured WM tracts; (iii) After adjusting for demographics, WMH was the strongest predictor of diffusion measures; (iv) Both dMRI measures were able to detect subtle Alzheimer’s disease related WM changes mainly at the medial temporal grey-white matter junctions and also WM tracts in the temporal lobes; and (v) MD and ISOVF in the corpus callosum were strongest predictors of cognitive function. Taken together, NODDI and traditional DTI measures are comparable in their predictive ability of WMH and cognition, but the non-overlapping information provided by each may aid in providing mechanistic insights into disease progression.

Advanced biophysical models versus traditional models

An advanced biophysical model such as NODDI leverages richer multi-shell diffusion gradients to examine the physiological alterations in neurites. In Fig. 2, we directly compared variation in DTI with variation in NODDI signal. The idea is that DTI signals are sensitive to gross anatomical and neuropathological changes associated with WM, but they are inherently non-specific to disentangle the complex tissue properties of a given voxel with crossing, kissing and fanning fibres., On the other hand, NODDI measures demonstrated more putative cell microstructure associations across studies, and have been found to be strongly correlated with neurobiological underpinnings.,, While decreases in NDI and increases in ISOVF are straightforward to understand, ODI changes have been hard to interpret because there is no simple physical mechanism that directly relates them to disease processes, such as demyelination, inflammation, or atrophy. Also, the direction of the correlation between axonal loss and ODI changes depends on region. Consider a hypothetical axon that runs parallel to a bundle for a few centimetres and then perpendicularly to a different bundle. The loss of the axon would increase the ODI in the parallel region and decrease it in the perpendicular region (ODI = 0 for a perfectly aligned bundle and goes to 1 where fibres spread out equally in all directions).

Unfortunately, at the macroscopic scale of a voxel the NODDI measures by themselves do not specify which axons changed. FA also suffers from this ambiguity, but unlike ODI is also directly coupled to demyelination and atrophy. Therefore, contrasting the relationship between DTI (FA and MD) and NODDI (NDI, ODI and ISOVF) can help understand the regional variations in these associations, which are largely unknown in the population.

The most consistent relationship between NODDI and DTI was seen with MD and NDI (but not with ODI and ISOVF), implicating that rather than orientation and geometry of tracts, a higher density may drive more diffusion restriction. We found a positive association of FA with NDI across various WM tracts with primarily a stronger relationship in corpus callosum fibres. This fibre pathway connects the two hemispheres, and the observed positive association between NDI and FA suggests the existence of the same underlying physiological processes (reduced axonal packing and demyelination)., Interestingly, ODI did not show a close association with FA in the corpus callosum, which may be due to their different responses to degeneration when most of the fibres are strongly aligned., Specifically, if the callosal boundary retreats due to atrophy, the edge voxels will be filled in by more CSF, affecting FA but not ODI.

Outside the corpus callosum, the associations between NODDI and DTI measures were inconsistent in the association and temporal WM tracts, which may be explained by the differing sensitivity of neurites to growth/maturational trajectories., Across the regions, the correlation between DTI and NODDI measures was smallest in the fornix. The fornix is part of the limbic system that connects the hippocampus to the subcortical structures and is also well known for partial volume contamination by CSF. This selective weakened association indicates the correction for CSF-contamination effects in the NODDI method.

Diffusion measures and age, sex and education/occupation

Age-associated WM changes in imaging have been widely reported. Past studies demonstrated NODDI as a key marker for studying ageing,, and the association of age and sex. Consistent with the ageing literature, age was associated with a decrease in FA and NDI, increase in MD ISOVF, and increase or decrease in ODI (which depends on the tract tortuosity and the presence of crossing fibres). Among the diffusion parameters, MD had the most sensitive age effects across the majority of tracts,, while the unconstrained diffusivity metric ISOVF demonstrated the greatest age effect in corpus callosum (genu and splenium) and cingulum. The overall widespread increase in ISOVF with age suggests the increase in FW concentration in specific brain regions. However, the key drivers of this increased FW are largely unknown. In addition to cerebrovascular disease and neurodegenerative pathologies, other possible underlying neuropathological factors include an influx of CSF or other factors like cell shrinkage, edoema, and neuroinflammation. Interestingly, past MRI and histology evidence clarified this as age- related increase in interstitial water.,

While most tracts had strong correlations with age, there were some subtle differences. With NDI, the association tracts (especially cingulum, superior longitudinal fasciculus and inferior fronto-occipital fasciculus) had the greatest age associations suggesting the presence of higher neurite density fibres in more metabolically active brain regions, that may be vulnerable to detrimental systemic age effects. ODI exhibited heterogeneous regional variations with age. Although there is reduced tract complexity in the corpus callosum fibres, the higher dispersion in fornix, cingulum and parahippocampal cingulum suggests the greater loosening, fanning and possibly bending of axonal bundles with the advancement of age., As stated above, this could be explained as evidence of continuous remodelling of WM during the life span which is more evident after the sixth decade. Notably, the age-related changes in the hippocampal connections might explain the amnestic changes in the elderly population.,

The sex-specific WM integrity association of DTI and NODDI is sparsely covered in the literature. A few DTI studies evaluated sex differences and reported inconsistent findings. The inconsistency across these findings may be due to the heterogeneity of populations and differences in the analysis methods. In this study, we observed small sex differences in eight tracts, with higher FA in males and greater ODI and ISOVF in females which is consistent with previous studies in healthy adults.,,

Reserve and resilience factors are important modulating parameters between brain injury and cognitive outcomes. Past studies showed education/occupation as a ‘proxy’ for cognitive reserve and related differences in fibre tract integrity., Similarly, a previous study in the MCSA population demonstrated a significant association of intellectual enrichment on FA of the genu. Although this association was modest, the present study showed a positive association of DTI and NODDI measures with corpus callosum, cingulum, parahippocampal cingulum, superior longitudinal fasciculus, inferior fronto-occipital fasciculus and inferior temporal WM, which is consistent with prior DTI studies., The anatomical localization provides insights into the associations between resilience mechanisms and brain maturation and plasticity.

Diffusion measures as markers of cerebrovascular disease

Although WMH is the most commonly used biomarker for cerebrovascular disease, it only represents extensive (and structurally visible) WM damage and fails to measure the disruption or subtle changes of the underlying WM tracts. There is growing evidence supporting the utility of DTI to characterize the WM changes in cerebrovascular disease,, even before the appearance of WMH and cognitive decline. Notably, the observed decrease in diffusion directionality and an increase in the extent of water diffusion in conjunction with WMH are consistent with prior findings.,, Consistent with this idea, a recent study using FW imaging demonstrated a greater contribution of cerebrovascular disease markers than Alzheimer’s disease biomarkers (CSF and PET) in memory clinic patients.

Similar to the age effect, we found a strong association of WMH with corpus callosum and association fibres. Although there is evidence for more vascular damage in the thinly myelinated anterior corpus callosum,,, the present study showed slight variations across the measures. Importantly, conventional DTI performed as well as NDI in detecting cerebrovascular disease changes. The decreased density and dispersion of the neurites and increased FW might contrast the lack of specificity in FA and MD to explain the underlying histological changes associated with WMH. The NDI finding in the genu of the corpus callosum was in accordance with a prior NODDI study that explored the diabetic encephalopathy in subjects with cognitive impairment. The only measure that did not show consistent associations with WMH was ODI. We believe this is due to ODI nominally being a property of only healthy neurites, and thus being more orthogonal to neuronal decay than the other NODDI measures or DTI.

Though the correlations with vascular risk were not a focus of the manuscript, as previously reported we found that WM measures from traditional DTI (FA and MD) and NODDI (NDI specifically) were significantly associated with worsening vascular risk (Results not shown). Recent researchers focussed on using global diffusion MRI as a cerebrovascular disease marker.,, However, this work sheds light on the variability in regional associations suggesting a greater sensitivity and specificity of regional markers. Future work should be undertaken to widely validate and compare diffusion outcomes as cerebrovascular disease measures.

Diffusion measures and neuroimaging Alzheimer’s disease measures

The association between amyloid deposition and WM microstructure is still a matter of debate. A non-monotonic behaviour was found between both measures in GM and WM, in human studies. Consistent with our region level findings, prior DTI studies reported reduced FA in corpus callosum and fornix, in cognitively unimpaired individuals and increased axial diffusivity and accelerated FA decrease in the parahippocampal cingulum of amyloid positive individuals. As expected, the current study identified a significant global association between Aβ deposition and increased MD and ISOVF along with decreased NDI in medial temporal lobe grey-white matter junctions, which are consistent with a more recent study that reported lower neurite density in limbic and association fibres and higher medial temporal FW. The medial temporal lobe is an early region of neuronal changes in Alzheimer’s disease, so the parahippocampal cingulum findings were as expected. We also found associations between tau and non-specific MD and ISOVF association in the inferior temporal WM (Fig. 4). These results are supported by a study of tau and NODDI in a transgenic Alzheimer’s disease model. Our findings in the temporal lobe (hippocampal and parahippocampal regions) and the temporo-occipital fusiform gyrus suggest that NODDI may be able to provide more detailed information about neurite health in the presence of Alzheimer’s disease pathology.

Diffusion measures with cognition

Association between WM DTI alterations and cognitive decline in the CU, MCI and Alzheimer’s disease populations have been reported previously.,,, Although there was decreased FA and increased MD in association with cognitive decline, the exact sources of DTI signal were not studied. The present study is one of the earliest studies to compare DTI and NODDI based on their association with cognitive performance after accounting for amyloid and tau, which allows us to evaluate its utility as a cerebrovascular disease marker. As expected, both DTI and NODDI were significantly associated with global cognition and cognitive subdomain scores after adjusting for age, sex, cycle visit and Alzheimer’s disease biomarkers. The overall pattern suggests that higher coherence and density, and lower FW concentration and tract complexity, both correlate with better cognitive performance., Across the tracts, the strongest association of reduced WM integrity and worse global cognitive performance was observed in the corpus callosum. This is consistent with a previous DTI study in cerebrovascular disease that reported highly significant correlations of genu and splenium with global cognitive performance. Impaired interhemispheric connection pathways contribute to multiple impaired cognitive functions, such as impaired memory, psychomotor speed, frontal lobe mediated attention and executive function., Additionally, these observations replicated our recent study in MCI that showed greater predictability of high FA of genu on better cognitive performance, even after controlling for amyloid and tau PET.

Although there are contributions from other domains, we found that the associations with WM integrity and cognitive performance were mainly driven by attention. Importantly, our detailed investigation indicated that FW fraction in the corpus callosum predicted cognitive decline. In general, reduced neurite density correlated with worse cognitive performance with most of the tracts in all domains. Among these, the stronger association of NDI than FA in superior longitudinal fasciculus may be due to its proximity to the crossing fibres in the centrum semiovale, which corresponds to higher FA and lower tract complexity. Another speculation may be that the superior longitudinal fasciculus is connecting lateral prefrontal to parietal brain areas, which are responsible for the multifaceted processes we studied here. Notably, cingulum performed uniformly well across all domains and diffusion metrics to predict cognition. This bundle is the prominent WM tract that interconnects frontal, parietal, and medial temporal lobe and the posterior cingulate cortex. Surprisingly, the parahippocampal cingulum bundle, which connects the hippocampus to the rest of the brain areas, emerged as an important tract in visuospatial function. In contrast, deteriorations in parahippocampal cingulum have previously been implicated in association with episodic memory in older subjects and Alzheimer’s disease.

Diffusion metrics are suggested to be most strongly associated with processing speed., Therefore, we also tested these hypotheses in the supplemental material and found that both DTI and NODDI strongly predicted processing speed (Trail B and combined). As expected, commissural fibres had the greatest effect size.

The present study has several strengths and limitations. The main strength was the extensive analyses of single and multi-shell diffusion data on WM health and cognition. Also, this is the first study to assess the relationship between NODDI metrics and cerebrovascular disease and Alzheimer’s disease biomarkers together along with associations with cognitive performance. The inclusion of a representative sample population strengthens the generalizability of the findings. Our voxel-wise and regional findings mostly corroborated each other, and the slight differences may be due to partial volume effects, smoothing, and more stringent FWE corrections. The major limitations are the cross-sectional nature of the study and the lack of histological confirmation of the observed associations. Another limitation is the smaller number of subjects in the dementia group, but the results remain the same after excluding them. Furthermore, the regularization scheme used by the AMICO implementation of NODDI acts like a prior that gives a mild preference to some values of NDI, ISOVF, and especially ODI, which could be obscuring some differences between subjects. Future longitudinal research with multiple biophysical models may provide more sensitive and conclusive findings.

In summary, the present study provides evidence of microstructural WM alterations due to ageing and age-related pathological processes, and their impact on cognition. Although NODDI-derived indices perform similar to traditional FA and MD in predicting cognitive performance, NODDI provides additional insights into the underlying synaptic density, organization and FW content which are biological processes that cannot be separated with DTI. Among DTI and NODDI indices, MD and FW fraction provided by ISOVF were the key parameters in predicting cognition. This study also highlights the spatial heterogeneity of tracts across the metrics and which highlights the importance of looking at each diffusion metrics to investigate changes in each WM region of the brain as a function of disease progression.

Supplementary material

Supplementary material is available at Brain Communication online.

Click here for additional data file.