Utility of quantitative MRI metrics in brain ageing research

The advent of new, advanced quantitative MRI metrics allows for in vivo evaluation of multiple biological processes highly relevant for ageing. The presented study combines several MRI parameters hypothesised to detect distinct biological characteristics as myelin density, cellularity, cellular membrane integrity and iron concentration. 116 healthy volunteers, continuously distributed over the whole adult age span, underwent a multi-modal MRI protocol acquisition. Scatterplots of individual MRI metrics revealed that certain MRI protocols offer much higher sensitivity to early adulthood changes while plateauing in higher age (e.g., global functional connectivity in cerebral cortex or orientation dispersion index in white matter), while other MRI metrics provided reverse ability—stable levels in young adulthood with sharp changes with rising age (e.g., T1ρ and T2ρ). Nonetheless, despite the previously published validations of specificity towards microstructural biology based on cytoarchitectonic maps in healthy population or alterations in certain pathologies, several metrics previously hypothesised to be selective to common measures failed to show similar scatterplot distributions, pointing to further confounding factors directly related to age. Furthermore, other metrics, previously shown to detect different biological characteristics, exhibited substantial intercorrelations, be it due to the nature of the MRI protocol itself or co-dependence of relevant biological microstructural processes. All in all, the presented study provides a unique basis for the design and choice of relevant MRI parameters depending on the age group of interest. Furthermore, it calls for caution in simplistic biological inferences in ageing based on one simple MRI metric, even though previously validated under other conditions. Complex multi-modal approaches combining several metrics to extract the shared subcomponent will be necessary to achieve the desired goal of histological MRI.

Introduction

Over the course of ageing, human brains undergo extensive remodelling, a combination of accumulating damage and reparatory processes affecting structure, microstructure and functional patterns. Be it cellular loss, degenerative changes of myelin sheaths, alterations of iron metabolism, accumulation of neurofibrillary tangles or shifts in biochemical processes (Yankner et al., 2008), it’s exceedingly difficult to get a complex, but still non-invasive insight into these processes.

The versatility of MRI has emerged as a viable candidate—the combination of MRI protocols with different sensitivities allows us to exploit the relative contributions of various histological features of tissue. Nonetheless, the true linkage of MRI signal to the underlying biological system poses an exceptional challenge due to its elaborate, complex nature. The multitude of various cell types, their densely packed axons and dendrites of immense importance for the cerebral function despite their small diameters, astrocytes and terminal vasculature are well beyond the current resolution limits of in vivo MRI in humans. The MRI signal from an imaging voxel is therefore an inherently aggregated measure of the constituent parts, where different MRI protocols provide a different sensitivity to primary contrast drivers as macromolecules—mainly myelin (Stüber et al., 2014), cell membranes (Palombo et al., 2020) and paramagnetic iron (Bilgic et al., 2012; Möller et al., 2019). All of these are expected to be closely linked with ageing.

Furthermore, true assessment of normal biological processes, pathologies or even responses to therapeutical interventions intrinsically requires the extraction of quantifiable measures—quantitative imaging biomarkers (Kessler et al., 2015), going beyond the conventional MRI aimed at local anatomical contrast. Ideally, these quantitative MRI (qMRI) parameters should provide calibrated measures directly comparable across time points and imaging sites (Zaidi, 2006). A large number of MRI modalities hypothesised to meet these criteria are available. The measurement of water proton relaxation time/rate constants is one of the most straightforward approaches, providing quantitative information on molecular constituents in the micro-environment and also water molecule movement between different micro-environments (Does, 2018). Diffusion weighted imaging (DWI) offers a wide range of applications from reconstructing white matter pathways to estimating cell sizes and evaluating the axonal integrity (Dyrby et al., 2018). Myelination has been the focus of multiple targeted protocols. The rather controversial, but easy-to-acquire ratio of T1-weighted (T1w) and T2-weighted (T2w) signal has been labelled as “myelin map” based on the strong spatial correlation with histologically-derived myelin maps (Glasser and Van Essen, 2011). More advanced measures as magnetization transfer (Schmierer et al., 2004) and non-adiabatic RAFF4 (Relaxation Along a Fictitious Field of rank 4; Liimatainen et al., 2010) have been validated directly against histological samples. And last, but not least, parameters as volumes of specific structures of interest and certain “semiquantitative” measures of resting-state functional MRI (rs-fMRI) with an extensive body of research backing their utility may be also considered of value in this context.

Ergo, the benefits of multi-modal MRI acquisitions are more than apparent despite the time requirements and costs since they provide multifaceted and complementing views of the under-pinning biological structures (Castiglioni et al., 2019; Zhang et al., 2020; Weiskopf et al., 2021). The presented cross-sectional study was designed to capitalise exactly on this quantitative multiparameter mapping approach, consisting of:

– three relaxometry metrics:

o adiabatic longitudinal relaxation constant T1ρ as a measure of cellular density (Michaeli et al., 2009)

o adiabatic transverse relaxation constant T2ρ as a measure of iron level (Mitsumori et al., 2009)

o non-adiabatic RAFF4 as a measure of myelin (Liimatainen et al., 2010)

– two DWI metrics based on NODDI (Tariq et al., 2016)

o orientation dispersion index (ODI) associated with the spatial configuration of neurite structures

o intracellular volume fraction (fICVF) associated with neurite density

– three rs-fMRI measures

o weighted Degree Centrality (wDeCe) describing whole-brain, global connectivity

o regional homogeneity (ReHo) describing local connectivity, i.e., synchrony of local neuronal populations (Zang et al., 2004)

o fractional amplitude of low-frequency fluctuations (fALFF) reflecting the extent of spontaneous neural activity of the brain (Zang et al., 2007)

– T1w/T2w ratio as a measure of myelin (Glasser and Van Essen, 2011; Filip et al., 2020a, 2020b)

– volumes of structures of interest.

A cohort continuous over the whole adult age span was examined to describe age-related brain tissue alterations, with the following specific aims:

(1) Plot the age-related courses of individual qMRI parameters over main regions of interest (ROIs), enabling us to evaluate the plausibility of their use in future ageing research considering the specific neurobiological hypotheses.

(2) Perform cross-correlations between individual metrics in the main ROIs to assess eventual similarities in Aim 1 outputs quantitatively.

(3) Evaluate the unique contributions of individual protocols in main large ROIs (subcortical grey matter [GM], cerebellar GM, cortical GM and whole white matter [WM]), i.e., the ability of the MRI metric to detect a process not visible with the other MRI metrics.

Our main hypothesis presumed that individual qMRI parameters will exhibit different age curves, with different plateaus and increasing/decreasing trendlines, Furthermore, the ability of rotating frame relaxations to detect further information from much lower frequencies in the kHz range of the effective fields will give them edge over other metrics.

Methods

Subjects

116 volunteers were screened for this study. Basic demographic information was collected, including the Brief Questionnaire Regarding Severity of Memory & Emotional Problems (BQRS-M&E; Reisberg, 2013)—a simple, screening, patient-facing 5-item scale providing subjective assessment of eventual cognitive problems. The following exclusion criteria were utilised: MRI contraindications, claustrophobia, significant space occupying or vascular lesion in the MRI scan, psychiatric or neurologic disorder, subjectively perceived cognitive problems (BQRS-M&E ≥ 3). Every participant provided a written informed consent in accordance with the Declaration of Helsinki. The study protocol was approved by the ethics committee of the University Hospital of St. Anne, Czech Republic.

Imaging protocol and data analysis.

MRI acquisition was performed in a 3-Tesla Siemens MAGNETOM Prisma scanner using a 64-channel head–neck coil at the Central European Institute of Technology in Brno, Czech Republic. Magnetization-prepared rapid gradient echo (MPRAGE) sequence in sagittal orientation with 1.0 mm isotropic resolution, echo time (TE) 2.47 ms, repetition time (TR) 2,150 ms, inversion time (TI) of 1,100 ms, flip angle 8°, generalized autocalibrating partial parallel acquisition (GRAPPA) acceleration factor of 2 was utilised to generate a T1-weighted (T1w) image. T2-weighted (T2w) image was acquired using the SPACE sequence in sagittal orientation, with 1.0 mm isotropic resolution, TE 72.6 ms, TR 2,820 ms, GRAPPA 2. Manufacturer-implemented pre-scan normalisation algorithm was employed for both T1w and T2w scans. Diffusion-weighted imaging (DWI) was based on the following parameters: 1.8 mm isotropic resolution, TR 2,820 ms, TE 726 ms, MB acceleration factor 4, 93 directions with 7 additional non-diffusion weighted (b0) images, two b-shells of 750 and 1,500 s/mm². Two DWI acquisitions were performed with opposite phase encoding polarity [anterior–posterior (AP) followed by posterior–anterior (PA)]. T2*-weighted scans sensitive to the blood oxygenation level dependent contrast were used for the rs-fMRI acquisition, with gradient-recalled echo (GRE) echo-planar imaging (EPI) sequence, 3.0 mm isotropic resolution, TR 900 ms, TE 30.0 ms, flip angle 45°, MB acceleration factor 4, interleaved acquisition, consisting of 502 volumes. Rotating frame relaxation measurements (adiabatic T1ρ, T2ρ and RAFF4) were obtained with voxel size of 1.6 × 1.6 × 3.6 mm³, GRAPPA 3, TE 3.18 ms and TR = 2,000 ms, 36 slices in total. Hyperbolic secant pulses were utilised for adiabatic T1ρ and T2ρ with adiabaticity factor R = 10, pulse duration of 6 ms, bandwidth 1.3 kHz, peak power ω1^max/(2π) of 800 Hz, and number of pulses 0, 4, 8, 12, and 16. In RAFF4 acquisitions, pulse duration was 4.56 ms, number of pulses 0, 4, 8, 12, 16 and ω1^max/(2π) 327 Hz. No medication was administered for the MRI scanning session.

The image processing pipeline for structural T1w and T2w images was based on the human connectome project (HCP) minimal pre-processing pipeline (Glasser et al., 2013) with minor modifications. Specifically, the PreFreeSurfer step did not include the optional gradient non-linearity and B₀ field inhomogeneity correction. The FreeSurfer step was based on the CUDA (Compute Unified Device Architecture)-enabled version of FreeSurfer 6.0¹. PostFreeSurfer step was implemented without any changes.

The analysis of rs-fMRI acquisitions was also performed using the HCP minimal pre-processing pipeline—concatenated transformation consisting of gradient non-linearity correction, motion correction, correction of the B₀ field inhomogeneities using the spin-echo fieldmaps, rigid-body co-registration to the T1w native space with spline interpolation, followed by the mapping to the Connectivity Informatics Technology Initiative (CIFTI) grey ordinates space using ribbon-constrained approach and regularisation with 2-mm full width at half maximum (FWHM) surface and subcortical volume smoothing. The subsequent processing utilised the HCP rsfMRI pipeline (Smith et al., 2013), which includes MELODIC independent component analysis (ICA), automatic artefact components identification via the FIX algorithm (Salimi-Khorshidi et al., 2014) and regressing out motion-related time-courses. HCP training data was used for the automatic classification algorithm with subsequent manual correction of the output and “non-aggressive” regression of the contribution of artefactual components out of the rsfMRI data. The AFNI package (Cox, 1996) was utilised to calculate voxel-wise weighted degree centrality (DeCe) with sparsity threshold of 0.1 (Craddock and Clark, 2016). Furthermore, Regional homogeneity was calculated as Kendall’s W over 28 neighbouring voxels (Zang et al., 2004) and fALFF as the total power within the frequency range between 0.01 and 0.1 Hz (Zang et al., 2007).

Diffusion weighted imaging processing also followed the HCP minimal pre-processing pipeline, including the optional gradient non-linearity correction. Afterwards, the output coregistered to the T1w native space of the subject was utilised for NODDI parameter calculation (Tariq et al., 2016), specifically orientation dispersion index and intracellular volume fraction.

The processing of T1ρ, T2ρ and RAFF4 maps utilised one pipeline: after 3D rigid-body motion correction of all the acquired scans to the first scan of each of these sequences (trilinear interpolation, mutual information as cost function followed by an optimization pass with sinc interpolation as implemented in the FSL 6.0 MCFLIRT), relaxation time constants were calculated utilising 2-parameter non-linear fitting (custom routines in MATLAB R2016a; MathWorks, Inc., Natick, MA). Afterwards, each map was co-registered to the T1w native space of the subject—mri_robust_register-initialized BB-register algorithm.

Visual inspection and evaluation of root-mean-squared voxel displacement for motion-correction in T1ρ, T2ρ and RAFF4 maps reconstruction, rs-fMRI and DWI data was performed to exclude subjects with framewise motion exceeding the extent of two voxels and/or substantial anatomical deviations.

And lastly, predetermined ROI masks derived from automatic FreeSurfer segmentation (cortical GM encompassing all the cerebral lobes; subcortical GM [thalamus, striatum, pallidum, hippocampus, amygdala, nucleus accumbens], cerebellar GM and whole WM) were co-registered to the images with lower resolution to avoid oversampling and inherently related problems. Brainstem was not considered in the presented analyses due to the complex structure combining GM and WM and lack of individual subsegmentation output of FreeSurfer for this area. See Figure 1 for the visualisation of GM ROIs.

FIGURE 1

Figure 1. Visualisation of regions of interest from one representative subject in axial plane. Slice coordinates in the z axis provided, with zero reference level in the middle of the anterior commissure and posterior commissure line. Red = cortical grey matter; blue = subcortical grey matter; green = cerebellar cortical grey matter. White matter region of interest not marked in a specific colour.

All data processing outputs were visually evaluated by trained operators (VK and PF).

Medians of individual MRI metrics of interest (T1ρ, T2ρ, RAFF4, ODI, fICVF, T1w/T2w ratio, wDeCe, ReHo, fALFF) were calculated over individual ROIs in each subject. ROI volumes were normalised using the estimated intracranial volume. For the Aim 2, rank full correlation (Spearman’s correlation coefficient) between individual MRI metrics was calculated, and for Aim 3, both rank full and partial correlation of age and individual MRI metrics was calculated, controlling each time for all the other MRI metrics in the latter case and the effect of sex. False discovery rate (FDR) correction was implemented: in Aim 2 over individual ROIs, and in Aim 3 over full and partial correlation analysis separately (Benjamini and Hochberg, 1995). Alpha of 0.05 was utilised as the type I error threshold.

Results

Seven subjects were excluded in total due to incomplete acquisition and/or excessive motion, leaving 109 subjects for further analyses.

Demographic and BQRS-M&E data was lost in one subject. Furthermore, based on the evaluation of BQRS-M&E data, four more subjects were excluded, leaving 104 subjects (51 female; median age 44.5 years, age range 19–89) in total for the final analysis. For further demographic information, including body weight, height, body mass index, education level, alcohol, tobacco and medication use, see the Supplementary Table S1. For the overview of age distribution, see the Supplementary Figure S1.

The results for Aim 1, scatterplots of individual MRI parameters and ages over predetermined ROIs, are presented in the Figure 2 showing individual scatterplots (axis x representing the relevant metric, axis y representing age in years) with overlaid 2nd order polynomial trendlines. Figure 3 provides a simplified overview of visually estimated ages when individual qMRI metrics detect some underlying changes, ordering the qMRI metrics based on this parameter.

FIGURE 2

Figure 2. Scatterplots of age in years (axis x) and individual MRI metrics (axis y): normalised volume, T1w/T2w ratio, T1ρ, T2ρ, rank 4 Relaxation Along a Fictitious Field (RAFF4), intracellular volume fraction (fICVF), orientation dispersion index (ODI), weighted Degree Centrality (wDeCe), regional homogeneity (ReHo), and fractional amplitude of low-frequency fluctuations (fALFF), for each of the four regions of interest (cerebral cortical grey matter (GM), subcortical GM, cerebellar cortical GM and white matter [not provided for resting-state fMRI metrics]). Second order polynomial trendline is overlaid in each scatterplot.

FIGURE 3

Figure 3. Simplified overview of ages when individual MRI metrics are able to detect underlying changes (in ascending order), for each of the four regions of interest [cerebral cortical grey matter (GM), subcortical GM, cerebellar cortical GM and white matter (not provided for resting-state fMRI metrics)]. Normalised volume, T1w/T2w ratio, T1ρ, T2ρ, rank 4 Relaxation Along a Fictitious Field (RAFF4), intracellular volume fraction (fICVF), orientation dispersion index (ODI), weighted Degree Centrality (wDeCe), regional homogeneity (ReHo), and fractional amplitude of low-frequency fluctuations (fALFF) are presented. Red colour or the bar denotes a descending trend of the MRI metric, green colour an ascending trend and white bar marks plateau/no change in the respective age group.

Figure 4 provides the output for the Aim 2—analysis of cross-corelation between individual MRI metrics—for all the predetermined ROIs. The lower triangle of presented matrices states relevant Spearman’s correlation coefficients, marking positive correlations in blue and negative correlations in red. The upper triangle of each matrix contains FDR-corrected (separately for each of the four ROIs) p-values, marking p-values < 0.001 in green and p-values in the range of 0.001–0.05 yellow for better visibility.

FIGURE 4

Figure 4. Rank cross-correlation matrices between individual MRI metrics across the four predetermined regions of interest [cerebral cortical grey matter (GM), subcortical GM, cerebellar cortical GM and white matter (not provided for resting-state fMRI metrics)]. Numbers in the lower triangle represent Spearman’s correlation coefficients (positive correlations in blue and negative correlations in red). The upper triangle of each matrix contains False discovery rate-corrected p-values (p-values < 0.001 in green, p-values between 0.001 and 0.05 yellow).

And lastly, the results for Aim 3, Spearman’s correlation analysis of age and MRI metrics, are presented in the Table 1. Full rank correlation coefficients, with respective FDR-corrected p-values, and partial rank correlation coefficients after controlling for the remaining MRI metrics of interest, with respective FDR-corrected p-values, are shown.

TABLE 1

Table 1. Full and partial rank correlations between the age and individual MRI metrics in the four predetermined regions of interest [cerebral cortical grey matter (GM), subcortical GM, cerebellar cortical GM and white matter (not provided for resting-state fMRI metrics)].

Discussion

This is the first study providing a complex overview of research-focused qMRI modalities, which are still relatively simple to implement and potentially useful in ageing research. Unsurprisingly, subject age is associated with a substantial effect on MRI metrics and all MRI protocols utilised in this study provided statistically significant results in one or more ROIs. However, the interpretation of their findings is far from straightforward, pointing to the presence of multiple elusive, latent variables.

Previously published histological validations or presumed models of tissue structure which individual MRI metrics (Michaeli et al., 2009; Mitsumori et al., 2009; Hakkarainen et al., 2016; Tariq et al., 2016) are supposed to be sensitive to do not fully explain the shapes of trendlines overlaid over the scatterplots in Figure 2. The use of 2nd order polynomial trendlines is definitely an oversimplification, but a logical one in the least. They fit rather well the presumed plateau corresponding to the relatively stable health and performance in early adulthood, followed by sharper decline in higher age, sometimes hypothesised to reach even exponential rate in some areas (Cox, 1996; Salimi-Khorshidi et al., 2014). For volume trendlines, this hypothesis seems to hold well for subcortical GM and WM, with sharper decline of volumes after the age of 50. However, cerebral and cerebellar cortical GM volume seem to decline more linearly over the whole age span acquired in the study.

While the interpretation of volume alteration is rather straightforward, the situation is more complex in further metrics. T1w/T2w ratio, previously labelled “myelin map” for cerebral cortical GM, since areas with higher T1w/T2w values corresponded to regions with higher myelin content (Glasser and Van Essen, 2011), fails to follow the age-related development of the underlying biology as seen in electron microscopy. While myelination is thought to continue well into early adulthood, normal ageing has been repeatedly associated with the loss of myelin integrity such as formation of splits and myelin balloons directly correlated with cognitive decline (Craddock and Clark, 2016). To add to the complexity of the process, oligodendroglial cells in the cerebral cortex seem to be unable to sustain myelin sheath quality, leading to decrease of internodal lengths of myelin sheaths, structural alterations with the presence of degenerating oligodendroglial cytoplasm in some sheaths and concomitant production of dysfunctional redundant myelin lamellae (Peters, 2002; Peters and Sethares, 2003). Since these very myelin lamellae are thought to be a substantial contributor to MRI signal attributed to myelin, we are hesitant to explain the rise in T1w/T2w ratio in the cerebral cortical GM as clear proof of the increase of myelination in its true meaning. Strangely though, although T1w/T2w ratio as “myelin map” was validated only in cerebral cortex, the shape and position of its trendline in the white matter in the Figure 2 may be considered to correspond much better to the presumed biology, despite well substantiated doubts about the association of this metric with white matter myelin concentration (Uddin et al., 2018). And lastly, while our T1w/T2w ratio age scatterplots are not in complete discord with previous studies (Grydeland et al., 2019), they fail to follow the shapes of other metrics previously validated as myelin markers both in histological (Satzer et al., 2015) and clinical studies in patients with demyelinating disorders (Filip et al., 2020a, 2020b, 2021)—RAFF4 and T1ρ. RAFF4 exhibits generally U-shaped curves, with low inter-individual variability in individual age groups mainly in cerebral cortical GM. However, as seen in the Table 1, this rather conspicuous and believable development of a plateau around the age of 50 is associated with only modest correlation and statistical significance.

The very same issue seems to haunt also MRI metrics hypothesised to be associated with cellular membrane integrity or cellularity in general—fICVF, T1ρ and partly also ODI. In all four ROIs, T1ρ seems to correspond rather well to the overlaid second order polynomial trendline, with steep increase after the age of 50–60, preceded by a mild decrease or a plateau stage. fICVF, i.e., fraction of the intracellular volume over the sum of intra-and extracellular volumes, rises rather steadily in all GM types throughout the age span, without sharp direction changes. Even its cross-correlation with T1ρ is borderline at best. On the other hand, ODI corresponds rather well to the ageing milestones in the cerebral cortical GM detected by T1ρ. However, cerebellar cortical GM exhibits a dissimilar shape in ODI than cerebral cortical GM, pointing to the selectivity of this method to non-negligible differences in cerebral and cerebellar ageing (Filip et al., 2019; Filip and Bareš, 2021). Moreover, ODI, as the only metric of the studied parameters, detects an age-related process leading to steep change in WM characteristics in youth followed by a plateau or mild decrease after the age of 50. This unique sensitivity is also confirmed by partial correlation analysis in WM in the Table 1 and in the Figure 3, where ODI shows no or only minimal cross-correlations with other MRI metrics in WM, making it an important candidate to consider in any ageing-focused research. Nonetheless, the most striking discrepancy in these findings is the failure to replicate the microscopic findings: a mild loss of cortical neurons has been suggested in normal ageing, but nowhere as extensive as suggested in early reports (Peters et al., 1998) or to be inferred from the presented data.

Adiabatic T2ρ was the only MRI metric in this study previously suggested as an iron marker (Michaeli et al., 2007; Mitsumori et al., 2009). Iron content in the cerebral cortical GM is an especially important factor since iron load is a risk factor of several neurodegenerative diseases (Bartzokis et al., 2004). But although about 80% of neocortical iron is stored in ferritin, myelinated fibers also have an increased iron concentration in several cortical areas (Stüber et al., 2014). While this may be one of the reasons for high cross-correlations (see Figure 3) of this metric with other MRI parameters presumed to be sensitive to other biological underpinnings, the extent of similarity casts doubts on the utility of the combined use of these highly intercorrelated MRI metrics in the ageing research. Nonetheless, the very biological processes detected by these MRI protocols may well be interlinked in normal ageing, even though they are distinct in their nature and physiological processes. Ergo, by no means does the Figure 3 prove the combined use of presented MRI metrics to be futile, since other physiological or pathological circumstances may lead to uncoupling of underlying biology and therefore different sensitivity of utilised MRI protocols (Filip et al., 2020a, 2020b). Furthermore, as seen in partial correlations in the Table 1, T2ρ is definitely able to detect a unique process in subcortical GM not visible by any other MRI measures implemented in this study.

And lastly, rsfMRI metrics, by their very nature, are generally difficult to interpret and relate specific biological processes to rsfMRI alterations, since the origin of BOLD (blood oxygen level dependent) signal and its association with the underlying neural activity is still a matter of dispute (Drew, 2019). Nonetheless, a widely accepted model of BOLD and local field potential coupling assumes an indirect entrainment of BOLD by neural activity, mostly post-synaptic potentials, but not the spiking rate (Ekstrom, 2010). Although highly cross-corelated in the cerebral cortical GM, they still provide an interesting addition to the whole picture of brain ageing. The age-related decrease of global connectivity (wDeCe) in cerebral cortical GM is then strangely countered by its increase in subcortical and cerebellar cortical GM. Local connectivity (ReHo) seems to decline generally with age in all the ROIs.

Several limitations to this study need to be acknowledged. Firstly, this study lacks complex formal neuropsychological examination to exclude subjects with non-diagnosed cognitive problems. While the evaluation of subjective cognitive impairment utilised here provides some help in this matter, further studies focused mainly on subjects of higher age should be designed with this possibility in mind. And secondly, this paper generally avoids biological interpretations of the presented findings and cautions against simplistic likening of individual qMRI parameters to one underlying pathophysiology under all circumstances. Hence, even the claims on qMRI parameters utilised in this study, as stated in the introduction, must be taken with a grain of salt. While T1ρ as a compound metric has been previously associated with neuronal density, there are no studies directly distinguishing between the effect of other cells in the central nervous system on this parameter. By the same token, other metals may affect T2ρ as an additional relaxation channel from previously published iron. And lastly, NODDI parameters distinguish in general between intracellular and extracellular space, without any putative selectivity for neurons or glia. While exceedingly difficult, neuroscientific research would greatly benefit from MRI protocols able to disentangle this issue and differentiate between neurons and glia.

Conclusion

All in all, this unique overview provides a quantitative baseline for the design of further studies and choice of relevant MRI parameters depending on the age group of interest. It clearly shows that MRI metrics are generally highly sensitive, but very poorly specific tools. One must even consider the possibility that the main biological underpinnings of some MRI metrics may be age-dependent as well -one biological process may be the primary driver of the detected signal in early age, but later end up overshadowed by another process developing only in higher age–a process which could explain the U-shape of curves in RAFF4. Even partial correlations do not show whether the MRI metrics detect the same biology or interlinked, but distinct microstructural features, since ageing-related cascades often consist of concomitant and co-dependent biological processes of substantially different nature and effect on MRI signal. Therefore, the newly emerging field of in vivo histological MRI aiming to uncover the true underlying biology inherently requires multimodal approach with advanced post-processing models combining several metrics, which would extract the dominant, shared subcomponent with much higher microstructural specificity than each metric individually. Furthermore, very careful multimodal validation studies against biological tissue samples, ranging from simplified phantoms to true biological specimen in combination with extensive histological and microscopic evaluation, must be considered to achieve the full potential of novel MRI approaches.

Data availability statement

Raw or processed data of this study are not publicly available due to the sensitive nature of human data acquired in patients. The data is available upon reasonable request to the corresponding author.

Ethics statement

The studies involving human participants were reviewed and approved by Ethics Committee of the University Hospital of St. Anne, Brno, Czechia. The patients/participants provided their written informed consent to participate in this study.

Author contributions

PF participated in the design of the study, performed the data analysis and wrote the manuscript. VK recruited the study subjects, performed the blinded quality control and reviewed the manuscript. ZV was consulted about the statistical methods and reviewed the manuscript. MB recruited the study subjects and reviewed the manuscript. SiM and ShM participated in the design of the study and reviewed the manuscript. LV designed the acquisition protocol, monitored data acquisition, and reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

Financial support for this project was provided by the General University Hospital in Prague (MH CZ-DRO-VFN64165). We acknowledge the core facility MAFIL of CEITEC supported by the MEYS CR (LM2018129 Czech-BioImaging). We also acknowledge the National Institutes of Health for support to SiM and ShM (Funding P41 EB027061 and R01 AG055591).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the funding agencies, the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2023.1099499/full#supplementary-material

Abbreviations

MRI, magnetic resonance imaging; rsfMRI, resting-state functional MRI; qMRI, quantitative MRI; T1w, T1-weighted; T2w, T2-weighted; DWI, diffusion-weighted imaging; RAFF4, Relaxation Along a Fictitious Field of rank 4; ODI, orientation dispersion index; fICVF, intracellular volume fraction; wDeCe, weighted Degree Centrality; ReHo, regional homogeneity; fALFF, fractional amplitude of low-frequency fluctuations; ROI, region of interest; GM, grey matter; WM, white matter.

Footnotes

1. ^http://surfer.nmr.mgh.harvard.edu

References

Bartzokis, G., Tishler, T. A., Shin, I. S., Lu, P. H., and Cummings, J. L. (2004). Brain ferritin iron as a risk factor for age at onset in neurodegenerative diseases. Ann. N. Y. Acad. Sci. 1012, 224–236. doi: 10.1196/annals.1306.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B (Methodological). 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Bilgic, B., Pfefferbaum, A., Rohlfing, T., Sullivan, E. V., and Adalsteinsson, E. (2012). MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. NeuroImage 59, 2625–2635. doi: 10.1016/j.neuroimage.2011.08.077

PubMed Abstract | CrossRef Full Text | Google Scholar

Castiglioni, I., Gallivanone, F., Soda, P., Avanzo, M., Stancanello, J., Aiello, M., et al. (2019). AI-based applications in hybrid imaging: how to build smart and truly multi-parametric decision models for radiomics. Eur. J. Nucl. Med. Mol. Imaging 46, 2673–2699. doi: 10.1007/s00259-019-04414-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Cox, R. W. (1996). AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. 29, 162–173. doi: 10.1006/cbmr.1996.0014

PubMed Abstract | CrossRef Full Text | Google Scholar

Craddock, R. C., and Clark, D. J. (2016). Optimized implementations of voxel-wise degree centrality and local functional connectivity density mapping in AFNI. Gigascience 5:s13742-016. doi: 10.1186/s13742-016-0147-0-d

CrossRef Full Text | Google Scholar

Does, M. D. (2018). Inferring brain tissue composition and microstructure via MR relaxometry. NeuroImage 182, 136–148. doi: 10.1016/j.neuroimage.2017.12.087

PubMed Abstract | CrossRef Full Text | Google Scholar

Drew, P. J. (2019). Vascular and neural basis of the BOLD signal. Curr. Opin. Neurobiol. 58, 61–69. doi: 10.1016/j.conb.2019.06.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dyrby, T. B., Innocenti, G. M., Bech, M., and Lundell, H. (2018). Validation strategies for the interpretation of microstructure imaging using diffusion MRI. NeuroImage 182, 62–79. doi: 10.1016/j.neuroimage.2018.06.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Ekstrom, A. (2010). How and when the fMRI BOLD signal relates to underlying neural activity: the danger in dissociation. Brain Res. Rev. 62, 233–244. doi: 10.1016/j.brainresrev.2009.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Filip, P., and Bareš, M. (2021). “Cerebellum—aging of the neuronal machine” in Factors Affecting Neurological Aging. eds. C. R. Martin, V. R. Preedy, and R. Rajendram (Netherlands: Elsevier), 281–288.

Google Scholar

Filip, P., Dufek, M., Mangia, S., Michaeli, S., Bareš, M., Schwarz, D., et al. (2021). Alterations in sensorimotor and mesiotemporal cortices and diffuse white matter changes in primary progressive multiple sclerosis detected by adiabatic relaxometry. Front. Neurosci. 15:711067. doi: 10.3389/fnins.2021.711067

CrossRef Full Text | Google Scholar

Filip, P., Gallea, C., Lehéricy, S., Lungu, O., and Bareš, M. (2019). Neural scaffolding as the Foundation for Stable Performance of aging cerebellum. Cerebellum [internet]. Cerebellum 18, 500–510. doi: 10.1007/s12311-019-01015-7

CrossRef Full Text | Google Scholar

Filip, P., Svatkova, A., Carpenter, A. F., Eberly, L. E., Nestrasil, I., Nissi, M. J., et al. (2020a). Rotating frame MRI relaxations as markers of diffuse white matter abnormalities in multiple sclerosis. NeuroImage Clin. 1:102234. doi: 10.1016/j.nicl.2020.102234

CrossRef Full Text | Google Scholar

Filip, P., Vojtíšek, L., Baláž, M., Mangia, S., Michaeli, S., Šumec, R., et al. (2020b). Differential diagnosis of tremor syndromes using MRI relaxometry. Parkinsonism Relat. Disord. 81, 190–193. doi: 10.1016/j.parkreldis.2020.10.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Glasser, M. F., Sotiropoulos, S. N., Wilson, J. A., Coalson, T. S., Fischl, B., Andersson, J. L., et al. (2013). The minimal preprocessing pipelines for the human Connectome project. NeuroImage 80, 105–124. doi: 10.1016/j.neuroimage.2013.04.127

PubMed Abstract | CrossRef Full Text | Google Scholar

Glasser, M. F., and Van Essen, D. C. (2011). Mapping human cortical areas in vivo based on myelin content as revealed by T1-and T2-weighted MRI. J. Neurosci. 31, 11597–11616. doi: 10.1523/JNEUROSCI.2180-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Grydeland, H., Vértes, P. E., Váša, F., Romero-Garcia, R., Whitaker, K., Alexander-Bloch, A. F., et al. (2019). Waves of maturation and senescence in micro-structural MRI markers of human cortical myelination over the lifespan. Cereb. Cortex 29, 1369–1381. doi: 10.1093/cercor/bhy330

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakkarainen, H., Sierra, A., Mangia, S., Garwood, M., Michaeli, S., Gröhn, O., et al. (2016). MRI relaxation in the presence of fictitious fields correlates with myelin content in normal rat brain. Magn. Reson. Med. 75, 161–168. doi: 10.1002/mrm.25590

PubMed Abstract | CrossRef Full Text | Google Scholar

Kessler, L. G., Barnhart, H. X., Buckler, A. J., Choudhury, K. R., Kondratovich, M. V., Toledano, A., et al. (2015). The emerging science of quantitative imaging biomarkers terminology and definitions for scientific studies and regulatory submissions. Stat. Methods Med. Res. 24, 9–26. doi: 10.1177/0962280214537333

PubMed Abstract | CrossRef Full Text | Google Scholar

Liimatainen, T., Sorce, D. J., O’Connell, R., Garwood, M., and Michaeli, S. (2010). MRI contrast from relaxation along a fictitious field (RAFF). Magn. Reson. Med. 64, 983–994. doi: 10.1002/mrm.22372

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaeli, S., Burns, T. C., Kudishevich, E., Harel, N., Hanson, T., Sorce, D. J., et al. (2009). Detection of neuronal loss using T1ρ MRI assessment of 1H2O spin dynamics in the aphakia mouse. J. Neurosci. Methods 177, 160–167. doi: 10.1016/j.jneumeth.2008.10.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Michaeli, S., Oez, G., Sorce, D. J., Garwood, M., Ugurbil, K., Majestic, S., et al. (2007). Assessment of brain iron and neuronal integrity in patients with Parkinson’s disease using novel MRI contrasts. Mov. Disord. 22, 334–340. doi: 10.1002/mds.21227

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitsumori, F., Watanabe, H., and Takaya, N. (2009). Estimation of brain iron concentration in vivo using a linear relationship between regional iron and apparent transverse relaxation rate of the tissue water at 4.7 T. Magn. Reson. Med. 62, 1326–1330. doi: 10.1002/mrm.22097

PubMed Abstract | CrossRef Full Text | Google Scholar

Möller, H. E., Bossoni, L., Connor, J. R., Crichton, R. R., Does, M. D., Ward, R. J., et al. (2019). Iron, myelin, and the brain: neuroimaging meets neurobiology. Trends Neurosci. 42, 384–401. doi: 10.1016/j.tins.2019.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Palombo, M., Ianus, A., Guerreri, M., Nunes, D., Alexander, D. C., Shemesh, N., et al. (2020). SANDI: a compartment-based model for non-invasive apparent soma and neurite imaging by diffusion MRI. NeuroImage 215:116835. doi: 10.1016/j.neuroimage.2020.116835

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, A. (2002). The effects of normal aging on myelin and nerve fibers: a review. J. Neurocytol. 31, 581–593. doi: 10.1023/A:1025731309829

CrossRef Full Text | Google Scholar

Peters, A., Morrison, J. H., Rosene, D. L., and Hyman, B. T. (1998). Are neurons lost from the primate cerebral cortex during normal aging? Cerebral Cortex (New York, NY: 1991) 8, 295–300. doi: 10.1093/cercor/8.4.295

CrossRef Full Text | Google Scholar

Peters, A., and Sethares, C. (2003). Is there remyelination during aging of the primate central nervous system? J. Comp. Neurol. 460, 238–254. doi: 10.1002/cne.10639

CrossRef Full Text | Google Scholar

Reisberg, D. The Oxford Handbook of Cognitive Psychology. Oxford University Press; (2013), England

Google Scholar

Salimi-Khorshidi, G., Douaud, G., Beckmann, C. F., Glasser, M. F., Griffanti, L., and Smith, S. M. (2014). Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. NeuroImage 90, 449–468. doi: 10.1016/j.neuroimage.2013.11.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Satzer, D., DiBartolomeo, C., Ritchie, M. M., Storino, C., Liimatainen, T., Hakkarainen, H., et al. (2015). Assessment of dysmyelination with RAFFn MRI: application to murine MPS I. PLoS One 10:e0116788. doi: 10.1371/journal.pone.0116788

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmierer, K., Scaravilli, F., Altmann, D. R., Barker, G. J., and Miller, D. H. (2004). Magnetization transfer ratio and myelin in postmortem multiple sclerosis brain. Ann. Neurol. 56, 407–415. doi: 10.1002/ana.20202

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M., Beckmann, C. F., Andersson, J., Auerbach, E. J., Bijsterbosch, J., Douaud, G., et al. (2013). Resting-state fMRI in the human connectome project. NeuroImage 80, 144–168. doi: 10.1016/j.neuroimage.2013.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Stüber, C., Morawski, M., Schäfer, A., Labadie, C., Wähnert, M., Leuze, C., et al. (2014). Myelin and iron concentration in the human brain: a quantitative study of MRI contrast. NeuroImage 93, 95–106. doi: 10.1016/j.neuroimage.2014.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Tariq, M., Schneider, T., Alexander, D. C., Gandini Wheeler-Kingshott, C. A., and Zhang, H. (2016). Bingham–NODDI: mapping anisotropic orientation dispersion of neurites using diffusion MRI. NeuroImage 133, 207–223. doi: 10.1016/j.neuroimage.2016.01.046

CrossRef Full Text | Google Scholar

Uddin, M. N., Figley, T. D., Marrie, R. A., and Figley, C. R. (2018). Can T1w/T2w ratio be used as a myelin-specific measure in subcortical structures? Comparisons between FSE-based T1w/T2w ratios, GRASE-based T1w/T2w ratios and multi-echo GRASE-based myelin water fractions. NMR Biomed. 31:e3868. doi: 10.1002/nbm.3868

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiskopf, N., Edwards, L. J., Helms, G., Mohammadi, S., and Kirilina, E. (2021). Quantitative magnetic resonance imaging of brain anatomy and in vivo histology. Nat. Rev. Phys. 3, 570–588. doi: 10.1038/s42254-021-00326-1

CrossRef Full Text | Google Scholar

Yankner, B. A., Lu, T., and Loerch, P. (2008). The aging brain. Annu. Rev. Pathol. 3, 41–66. doi: 10.1146/annurev.pathmechdis.2.010506.092044

CrossRef Full Text | Google Scholar

Zaidi, H. Quantitative Analysis in Nuclear Medicine Imaging. Springer; (2006), Berlin

Google Scholar

Zang, Y. F., He, Y., Zhu, C. Z., Cao, Q. J., Sui, M. Q., Liang, M., et al. (2007). Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain and Development 29, 83–91. doi: 10.1016/j.braindev.2006.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. NeuroImage 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030

CrossRef Full Text | Google Scholar

Zhang, Y. D., Dong, Z., Wang, S. H., Yu, X., Yao, X., Zhou, Q., et al. (2020). Advances in multimodal data fusion in neuroimaging: overview, challenges, and novel orientation. Inf. Fusion. 64, 149–187. doi: 10.1016/j.inffus.2020.07.006

CrossRef Full Text | Google Scholar