Age-related magnetic susceptibility changes in deep grey matter and cerebral cortex of normal young and middle-aged adults depicted by whole brain analysis
Original Article

Age-related magnetic susceptibility changes in deep grey matter and cerebral cortex of normal young and middle-aged adults depicted by whole brain analysis

Romana Burgetova1,2, Petr Dusek1,3, Andrea Burgetova1^, Adam Pudlac1, Manuela Vaneckova1, Dana Horakova3, Jan Krasensky1, Zsoka Varga3, Lukas Lambert1

1Department of Radiology, First Faculty of Medicine, Charles University and General University Hospital in Prague, Prague, Czech Republic; 2Department of Radiology, Third Faculty of Medicine, Charles University and University Hospital Královské Vinohrady, Prague, Czech Republic; 3Department of Neurology, First Faculty of Medicine, Charles University and General University Hospital in Prague, Prague, Czech Republic

Contributions: (I) Conception and design: P Dusek, A Burgetova, L Lambert; (II) Administrative support: A Burgetova, L Lambert, A Pudlac; (III) Provision of study materials or patients: D Horakova, Z Varga, P Dusek, M Vaneckova, J Krasensky; (IV) Collection and assembly of data: R Burgetova, M Vaneckova, P Dusek, J Krasensky, D Horakova; (V) Data analysis and interpretation: P Dusek, L Lambert, R Burgetova, A Burgetova; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

^ORCID: 0000-0002-9975-2338.

Correspondence to: Assoc. prof. Andrea Burgetova, MD, PhD.First Faculty of Medicine Charles University and General Faculty Hospital in Prague, U nemocnice 2, 128 08 Prague 2, Czech Republic. Email: andrea.burgetova@vfn.cz.

Background: Iron accumulates in brain tissue in healthy subjects during aging. Our goal was to conduct a detailed analysis of iron deposition patterns in the cerebral deep grey matter and cortex using region-based and whole-brain analyses of brain magnetic susceptibility.

Methods: Brain MRI was performed in 95 healthy individuals aged between 21 and 58 years on a 3T scanner. MRI protocol included T1-weighted (T1W) magnetization-prepared rapid acquisition with gradient echo images and 3D flow-compensated multi-echo gradient-echo images for quantitative susceptibility mapping (QSM). In the region-based analysis, QSM and T1W images entered an automated multi-atlas segmentation pipeline and regional mean bulk susceptibility values were calculated. The whole-brain analysis included a non-linear transformation of QSM images to the standard MNI template. For the whole-brain analysis voxel-wise maps of linear regression slopes β and P values were calculated. Regional masks of cortical voxels with a significant association between susceptibility and age were created and further analyzed.

Results: In cortical regions, the highest increase of susceptibility values with age was found in areas involved in motor functions (precentral and postcentral areas, premotor cortex), in cognitive processing (prefrontal cortex, superior temporal gyrus, insula, precuneus), and visual processing (occipital gyri, cuneus, posterior cingulum, fusiform, calcarine and lingual gyrus). Thalamic susceptibility increased until the fourth decade and decreased thereafter with the exception of the pulvinar where susceptibility increase was observed throughout the adult lifespan. Deep grey matter structures with the highest increase of susceptibility values with age included the red nucleus, putamen, substantia nigra, dentate nucleus, external globus pallidus, caudate nucleus, and the subthalamic nucleus in decreasing order.

Conclusions: Accumulation of iron in basal ganglia follows a linear pattern whereas in the thalamus, pulvinar, precentral cortex, and precuneus, it follows a quadratic or exponential pattern. Age-related changes of iron content are different in the pulvinar and the rest of the thalamus as well as in internal and external globus pallidus. In the cortex, areas involved in motor and cognitive functions and visual processing show the highest iron increase with aging. We suggest that the departure from normal patterns of regional brain iron trajectories during aging may be helpful in the detection of subtle neurodegenerative and neuroinflammatory processes.

Keywords: Magnetic susceptibility; iron; aging; brain; deep grey matter; cerebral cortex


Submitted Jan 23, 2021. Accepted for publication Apr 19, 2021.

doi: 10.21037/qims-21-87


Introduction

During the lifespan, the human brain undergoes microstructural changes including grey matter atrophy, myelin degeneration, or iron accumulation that can be picked in vivo by various quantitative MRI techniques (1-4). Knowledge about the spectrum of age-related changes in brain structure and associated MR parameters is important for understanding the neurobiology of aging and for differentiating physiological aging from disease.

Magnetic susceptibility, a measure that indicates the extent to which a material may be magnetized in external magnetic field, can be quantified in tissues by a recent MR technique called quantitative susceptibility mapping (QSM) (5). While magnetic susceptibility of cerebral tissue may be influenced by myelin, calcium, or deoxyhemoglobin, it is mainly determined by ferritin-iron content in the basal ganglia and other grey matter structures (6). QSM is thus considered as a surrogate measure of iron content in the grey matter (7,8).

Cerebral accumulation of non-heme iron occurs during normal aging. This increase occurs heterogeneously in certain cell types and brain regions such as basal ganglia, hippocampus, motor cortex, cerebellar nuclei, and other subcortical brain regions (9-13). Globus pallidus, substantia nigra, red nuclei, and putamen are regions with the highest iron deposition (14-16).

Iron is involved in numerous pathways in the brain such as oxygen transportation, oxidative phosphorylation, DNA synthesis, mitochondrial respiration, myelin synthesis, antioxidant enzyme activation, and metabolism of neurotransmitters (13). On the other hand, disturbance of iron homeostasis can lead to brain neurotoxicity by different mechanisms such as membrane damage, protoxin bioactivation, aberrant cell signaling, bioenergetic failure, proteosomal and mitochondrial dysfunction, protein aggregation and inclusion formation, electrophysiological derangements, and synaptolysis, conveying all of these factors on apoptosis, necrosis or a specific cell death triggered by iron called ferroptosis (13,17,18). Increased subcortical iron is associated with decreased cognitive and motor functions in the elderly population (14,19,20) suggesting that basal ganglia may be more vulnerable to iron accumulation throughout aging, exposing them to an increased risk of neurodegeneration or/and neuroinflammation. Understanding normal patterns of iron accumulation with aging might thus assist in identifying subjects who are at greater risk of accelerated neurodegeneration.

The majority of previous studies examining iron accumulation during normal aging were based on region-of-interest analysis of deep grey matter on R2* relaxometry and magnetic susceptibility maps (21-27). Whole-brain analysis, which is not limited to a priori predefined regions, has been used only sparsely for studying age-related changes of magnetic susceptibility and, by extension, of iron levels (28,29).

The objective of this study was to conduct a detailed analysis of iron deposition patterns not only in the cerebral deep grey matter but also in the cortex by a voxel-wise analysis of QSM at the whole-brain level and using region-based analysis.


Methods

The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the Ethics Committee of the General University Hospital in Prague (ID1018/17), and all subjects signed informed consent.

Study participants

Brain MRI was performed in 95 healthy individuals (57 females and 38 males) aged between 21 and 58 (mean ± SD = 37±10) years, who responded to an announcement presented at the First Faculty of Medicine, Charles University and General University Hospital in Prague and met the following inclusion criteria: (I) willing to sign an informed consent in accordance with local regulations, (II) age 18–60 years (chosen to be comparable to that of adult patients with multiple sclerosis). The exclusion criteria were: (I) unable to undergo MRI examination, (II) unable to be examined 4 times during the next three years at 0, 12, 24, 36 months, (III) pregnancy at the time of enrollment into the study, (IV) any diseases, medical conditions, or brain MRI abnormalities that could influence the structure of the brain or spinal cord.

Imaging protocol

The examination was performed on a 3T MRI scanner (Siemens Skyra 3T, Siemens Healthcare, Erlangen, Germany) with a 32-channel birdcage head coil. MRI protocol included T1-weighted 3D Magnetization-Prepared Rapid Acquisition with Gradient Echo (MPRAGE) images (T1WI) in sagittal plane with the following acquisition parameters: repetition time (TR), 2,300 ms; echo time (TE), 2.96 ms; inversion time (TI) 900 ms; flip angle (FA) 9°; field of view (FOV) 176×256×256 mm; spatial resolution 1.0×1.0×1.0 mm3 for anatomical imaging, and segmentation and 3D flow-compensated multi-echo gradient-echo (GRE) images in axial plane (TR =33 ms; first TE =4.5 ms, uniform echo spacing =5 ms, last TE =29.5 ms, number of echoes =6; FA =18°; FOV =195×240×164 mm; spatial resolution =0.94×0.94×0.94 mm3) for QSM.

QSM processing

QSM was processed by a multi-scale dipole inversion algorithm implemented in QSMbox software package (https://gitlab.com/acostaj/QSMbox) (7,30). T1WI (i.e., MPRAGE) images were transferred into the QSM space by rigid-body registering to the first echo magnitude image from the GRE pulse sequence, using the SPM12 (Statistical Parametric Mapping) tool, (www.fil.ion.ucl.ac.uk/spm/software/spm12) (31) and skull stripped by multiplying with an SPM-based brain binary mask calculated by QSMbox. Co-registered skull-stripped QSM and T1WI entered an automated multi-atlas segmentation pipeline using dual, i.e., QSM/T1, contrast for delineation of DGM nuclei (32) implemented at a cloud-based platform (www.mricloud.org) (33). Subsequently, segmented volumes were eroded by one voxel, and mean bulk susceptibility values of the following regions of interest (ROI) were extracted: internal globus pallidus, external globus pallidus, putamen, caudate nucleus, thalamus, pulvinar, substantia nigra, red nucleus, dentate nucleus, subthalamic nucleus, and total white matter (Figure 1). Magnetic susceptibility values from both hemispheres were averaged and further analyzed as (I) raw uncorrected values and (II) values referenced to the mean bulk magnetic susceptibility of the total white matter. Notably, the total white matter was segmented by the automatic algorithm and its mean bulk susceptibility thus represents a highly reproducible operator independent measurement.

Figure 1 Illustration of MRI data processing. Example (A) magnitude and (B) raw phase (TE =24.5 ms) GRE images; (C) reconstructed QSM image and (D) corregistered skull-stripped T1-weighted anatomical image that were used for automatic multiatlas segmentation; (E) automatically segmented regions of interest overlaid onto anatomical image. ppb, parts per billion; GRE, gradient recalled echo; TE, echo time; QSM, quantitative susceptibility mapping.

T2* maps processing

T2* relaxation maps were calculated from multi-echo GRE magnitude images using non-linear least squares fitting according to the Levenberg-Marquardt algorithm (MRI Processor v.1.1.6, ImageJ 1.51k; T2* values were capped at 100 ms, the maximum number of iterations was 100, and forced no bias). T2* relaxation maps were used for whole-brain analysis.

Whole-brain analysis

The analysis was based on a procedure described previously (28). Briefly, a linear pre-registration (antsRegistration) and non-linear transformation (antsApplyTransformation) of susceptibility and T2* relaxation maps to the standard MNI template (MNI152 “ICBM 2009c Nonlinear Symmetric”) was performed using ANTS package ver. 2.1.0 (34). To compensate for a small mismatch in registration, we applied a median filter followed by a maximum intensity filter, both calculated across a sphere with a diameter of 3 mm. This approach was chosen empirically—we have visually compared the output images after eleven different combinations of filters (median spherical, median cubic, Gaussian spherical) with different kernel sizes applied. The selected combination of filters yielded images with a reasonable level of noise and detail and it appeared to be the most suitable approach for co-registration of susceptibility and T2* relaxation maps; the median filter removed both positive and negative signal values caused by vessels while the maximum intensity spherical filter (3 voxels in diameter) compensated for imperfections in co-registration of narrow structures (especially the gyri) to the MNI152 template.

In a post-hoc analysis, a mask of cortical voxels with a significant association of magnetic susceptibility and age was created. This mask was overlaid onto an anatomical T1 template and nine ROIs consisting of voxels from this mask were manually segmented using ITK-SNAP (www.itksnap.org) (35): precentral gyrus, postcentral gyrus, dorsolateral prefrontal cortex, medial prefrontal cortex, insular cortex, precuneus, occipital cortex, occipito-temporal cortex, and cerebellar vermis (Figure S1). Individual mean bulk magnetic susceptibility values from these ROIs were retrieved and further analyzed as described for deep grey matter ROIs (25).

Statistics

Univariate regression analysis for regional age-related magnetic susceptibility changes was performed with linear and non-linear (quadratic and exponential) models and the model with the best fit based on Akaike information criterion with correction for small sample sizes (AICc), R2, and root mean squared error (RMSE) was selected. For linear relationships, regression slope (β) was expressed in parts per billion (ppb)/year and its difference from zero was tested by the F test.

The exponential model was based on one-phase association and the susceptibility values were fitted with the following equation: susceptibility = A*(1−exp[−B*age]) + C, where A, B, and C are tissue-specific parameters (25). Holm-Bonferroni method was applied to correct for multiple hypothesis testing. Statistical analyses were carried out using Prism 8 (GraphPad Software, San Diego, CA, USA).

Whole-brain analyses of age-related magnetic susceptibility and T2* relaxation time changes were performed using SPM12. Statistical maps were generated using a cluster definition threshold of P<0.005 (uncorrected) and clusterwise family-wise error (FWE)-corrected threshold of P<0.05. Voxel-wise maps of linear regression slopes β were visualized using MRIcroGL (http://www.nitrc.org/projects/mricrogl); anatomical locations of significant clusters were retrieved from the Atlas of Intrinsic Connectivity of Homotopic Areas (AICHA) (36).


Results

ROI-based analysis

The age-related increase of magnetic susceptibility in the basal ganglia followed a linear pattern, but in the thalamus and pulvinar, it followed a quadratic and exponential pattern, respectively (Figure 2, Table S1). Deep grey matter structures with a significant linear increase of raw bulk magnetic susceptibility included the red nucleus (β=1.30, P<0.001), putamen (β=0.89, P<0.001), substantia nigra (β=0.76, P<0.001), dentate nucleus (β=0.73, P=0.002), external globus pallidus (β=0.62, P<0.001), caudate nucleus (β=0.46, P<0.001), and subthalamic nucleus (β=0.43, P=0.024). The effect of age was not significant for the internal globus pallidus (β=0.20, P=0.15). In the thalamus, magnetic susceptibility increased until the age of 40 years and decreased thereafter. In the pulvinar, susceptibility increase was followed by a relative flattening of the curve from the age of 40 years onwards. Analysis of magnetic susceptibility referenced to the total white matter value returned similar results (Figures S2,S3), whereby mean magnetic susceptibility in the total white matter showed a mild linear positive effect of age (β=0.03, P=0.002).

Figure 2 Bulk magnetic susceptibility values in ppb in deep gray matter structures as a function of age with their respective regression fits and coefficients. Blue dotted lines represent 95% confidence bands of the regression fit and red dotted lines represent 95% population prediction bands. Please note different Y-axis range for thalamus and pulvinar. ppb, parts per billion; GP, globus pallidus, STN, subthalamic nucleus.

Whole-brain analyses

Results of the whole-brain QSM analysis documenting magnetic susceptibility increase in ppb per year are shown in Figure 3 and Table 1. In addition to deep grey matter structures, multiple cortical regions with a positive association of magnetic susceptibility and age were identified in the precentral and postcentral areas, medial and dorsolateral prefrontal cortex, superior temporal gyrus and insula, posteriorly in the superior, lateral, and medial occipital gyri, precuneus, cuneus, posterior cingulate, fusiform, calcarine and lingual gyrus (P<0.05, FWE-corrected). In the cerebellum, the vermis and tonsils showed a significant positive correlation with age (P<0.05, FWE-corrected). Spatial coordinates and anatomical labeling of clusters with significant age-related magnetic susceptibility increase are listed in Table 1.

Figure 3 Statistical map of age-related susceptibility changes in the brain of normal adults represented as mean bulk susceptibility increase in ppb per year (thresholded at P<0.05 at cluster level and corrected for family-wise error; initial cluster definition threshold was set to P<0.005 uncorrected). Z-coordinate in the MNI space in [mm] is indicated for each slice. ppb, parts per billion; MNI, Montreal Neurological Institute.
Table 1
Table 1 Clusters with significant age-related susceptibility increase identified by whole brain analysis with corresponding anatomical locations
Full table

Several regions showed magnetic susceptibility decrease with age, although the effect was much smaller compared to areas with susceptibility increase reaching β not greater than −0.2 ppb per year (P<0.05, FWE-corrected). These regions include distinct white matter areas, particularly the rostrum of the corpus callosum, external capsulae, mesencephalic tectum and tegmentum, and frontal and temporal subcortical regions. Several extracerebral regions, the choroid plexus, falx cerebri, and venous sinuses showed age-related susceptibility decrease with β up to −1.0 ppb per year (P<0.05, FWE-corrected).

Post-hoc analysis of magnetic susceptibility in cortical regions with a significant effect of age showed its linear increase in the postcentral, dorsolateral prefrontal, medial prefrontal, insular, occipital cortices, precuneus, and cerebellar vermis (β≈0.2, P<0.001). In the precentral gyrus and temporo-occipital gyrus, magnetic susceptibility followed a quadratic pattern with a relatively steep increase until the age of 40 years and flattening of the curve afterwards (Figure 4).

Figure 4 Uncorrected bulk magnetic susceptibility values in ppb in cortical regions of interest as a function of age with their respective regression fits. Blue dotted lines represent 95% confidence bands of the regression fit and red dotted lines represent 95% population prediction bands. ppb, parts per billion; DLPFC, dorsolateral prefrontal cortex; MPFC, medial prefrontal cortex.

Results of T2* analysis are shown in Figure S4 and Table S2. The results are similar to QSM however, a larger number of areas correlate with age. Compared to QSM, associations are evident also in the white matter.


Discussion

In this study, we quantified age-related changes of magnetic susceptibility in deep grey matter structures and the cerebral cortex in a cohort of adult healthy individuals using whole-brain and ROI-based analyses. We found a linear association between age and magnetic susceptibility in deep grey matter structures including the red nucleus, putamen, substantia nigra, dentate nucleus, external globus pallidus, and caudate nucleus. In the thalamus, quadratic association with age was observed with magnetic susceptibility increasing before and decreasing after the age of 40 years. Whole-brain QSM analysis documented that beyond these regions, there is a significant positive association between age and magnetic susceptibility in extensive areas of the cerebral cortex—in the frontal lobes, Rolandic area, insular area, precuneus, cuneus, posterior cingulum, and lingual gyrus. Whole-brain T2* analysis showed a similar pattern of age-related changes. In opposite to susceptibility, which showed a positive effect of age only in grey matter, areas of age-related T2* relaxation times shortening overlapped grey and white matter structures. This is likely caused by the fact that in QSM, the paramagnetic effect of accumulating iron and the diamagnetic effect of maturating myelin counteract each other in the white matter, while in T2* maps, both, iron accumulation and myelin maturation with aging, cause local magnetic field inhomogeneities and thus contribute to the shortening of T2* relaxation times. This effect was particularly noted in subcortical white matter regions which are known to accumulate iron. Comparison of QSM and T2* map analyses indicates that the former may be more suitable for the assessment of age-related iron changes in the cortical regions.

It is well known that the trajectory of brain iron accumulation during the lifespan is not linear. A sharp increase of tissue iron concentration (15) and subsequently of magnetic susceptibility (25,27,37) during the first two decades is followed by flattening of the curve or even decrease at old age. Results of this study indicate the evolution of susceptibility values is linear between the age of 20 and 60 years for most deep grey matter and cortical structures. It is likely that regions that showed an exponential relationship between age and magnetic susceptibility in the current study, i.e., pulvinar, motor and occipito-temporal cortices, exhibit less abrupt initial iron accumulation with flattening of the curve shifted to an older age.

We found the largest effect of age on magnetic susceptibility (the steepest regression slope) in the red nucleus and putamen, which is consistent with previous studies (11,27,38). The red nucleus has been previously less investigated and, in some studies, the absolute magnitude of susceptibility increase was higher for the red nucleus than for the putamen (22,27,29), similarly as in our study. Our observation of milder gradual increase of magnetic susceptibility in the caudate nucleus, substantia nigra, and dentate nucleus is also concordant with post-mortem evidence (16) and several previous in vivo MRI studies (9,13,14,27,28).

In the pallidum, previous works reported contradictory findings including no age-related susceptibility changes (23,27,28), or mild iron accumulation during adulthood (22). A meta-analysis of 20 MRI studies that estimated iron content in the caudate nucleus, globus pallidus, putamen, red nucleus, and substantia nigra, documented the lowest age-related differences in globus pallidus although it was generally rich in iron in absolute value (11). Studies that included infants showed an exponential pattern with a steeply increasing iron concentration as well as magnetic susceptibility from birth and a plateau from the 3rd decade onwards (13,14,16,38-40). However, these studies did not aim at analyzing specifically the pallidal iron content trajectory during the adult lifespan. Our results indicate that the reason for contradictory findings may be the differential contribution of external and internal globus pallidus depending on segmentation strategy. In contrast to the internal part, the external part, which is adjacent to the putamen, showed a linear increase of susceptibility during aging in the current study. A similar effect was observed for R2* transverse relaxivity in another study that segmented internal and external globus pallidus (29). These results indicate distinct regulation of iron concentration in pallidal subregions and suggest that internal and external globus pallidus should be analyzed separately for proper assessment of disease-related changes in this structure.

Thalamus presents an overall low iron content and it was suggested that myelin contributes by a large degree to its susceptibility masking the effect of small age-related changes in iron content (41). On the other hand, the bidirectional pattern with an initial increase until the age of 40 years followed later by a decrease of magnetic susceptibility observed in the current and a previous study (42) conspicuously matches age-related changes in the thalamic iron concentration measured in post mortem tissue (15). It is thus likely that thalamic susceptibility is mostly driven by the effect of iron, similarly as in the basal ganglia. Contrary to some previous studies (27), we and others have found a positive age-related magnetic susceptibility increase in the pulvinar (42,43). It is interesting to point up the differences between pulvinar and other thalamic nuclei. The pulvinar, although being a posterior part of the thalamus, is visually discernible on QSM and behaves unlike other thalamic nuclei. Zhang et al. (42) suggests that in the medial and lateral nuclei of the thalamus susceptibility follows a bidirectional pattern while anterior nuclei and pulvinar display gradual susceptibility increase with age. Overall, there is converging evidence that MRI metrics sensitive to iron content should be analyzed separately for pulvinar and other thalamic nuclei.

Age-related changes in iron distribution in the cerebral and cerebellar cortex in vivo have been much less investigated. It may be related to the fact that the cerebral cortex has more complex and variable anatomy than deep grey matter nuclei and its automatic analysis at a group level requires non-linear transformation to match corresponding regions among subjects. Further challenges include the thinning of the cortex and widening of subarachnoid cerebrospinal fluid spaces in cerebral sulci with cerebral atrophy that accompanies aging. Results of the current whole-brain analysis largely corroborate findings from a single previous voxel-based QSM study in that cortical magnetic susceptibility evolves unevenly during aging (28,29). The strongest effect of age was found in areas involved in the motor (Rolandic area, premotor cortex, cerebellar vermis), cognitive (dorsolateral prefrontal cortex, insular region, precuneus), and visual (cuneus, posterior cingulum, lingual and fusiform gyri) cortical areas. This is also in line with the direct measurement of iron concentration in post mortem tissue (15) which indicates that iron has an important contribution to magnetic susceptibility in the cortex. Overall, structures that accumulate iron with aging contain dopaminergic iron-dependent pathways and it was hypothesized that both, changes in iron content and dopaminergic transmission may be responsible for behavioral changes during aging (27,28,44).

In the white matter, magnetic susceptibility is predominantly affected by diamagnetic myelin and it was previously shown to follow a trajectory of initial susceptibility decrease reaching its minimum between 25 and 45 years of age, followed by an increase, consistent with myelin maturation and decay in the course of normal brain development (25). We found a small but significant positive effect of age on the magnetic susceptibility of the total white matter, similar to a previous study (30) suggesting that susceptibility in most white matter regions was beyond its minimum in the examined age group. However, the whole-brain QSM analysis revealed white matter regions with a small negative effect of age corresponding to areas with an age-related decrease of T2* relaxation time, predominantly in the fronto-temporal subcortical areas and rostral corpus callosum which are regions where myelin maturates in late adulthood (45).

In line with previous studies (25,28) which demonstrated that the effect of age on magnetic susceptibility outweighs the impact of reference adjustment, we have used raw unreferenced values for the main analysis in order to avoid bias introduced by manual reference region placement. Secondary analysis with magnetic susceptibility referenced to the mean value of total white matter confirmed a negligible effect of susceptibility adjustment.

This study has several limitations. Firstly, its cross-sectional nature does not allow direct analysis of individual regional magnetic susceptibility values evolution throughout aging. A longitudinal study will be necessary to confirm the relationship between aging and magnetic susceptibility. Secondly, the age spectrum of participants spanning 20–60 years does not include adolescents and seniors; the results are thus not generalizable to these populations. The primary aim of this study was to provide reference values for multiple sclerosis patients, but the results may be as well used for other diseases with an adult onset, e.g., Huntington disease, Wilson disease, or neuromyelitis optica (46,47). Thirdly, magnetic susceptibility is only a surrogate marker of iron concentration and is also influenced by other substances such as diamagnetic myelin. However, an excellent match of susceptibility values with previously published directly measured iron concentrations indicates that for grey matter structures, magnetic susceptibility is dominantly influenced by iron and the confounding effect of myelin is rather small. Fourthly, voxel-based analysis of the brain involves non-rigid warping of anatomy to match the template. This is inherently associated with registration errors that may lead to false-negative age-related associations for some brain areas. Lastly, susceptibility changes are presumably influenced not only by absolute iron content but also by changes in the relative iron concentration secondary to age-related reductions in volumes of evaluated structures (48). However, such complex relationships cannot be inferred from this cross-sectional study.

In conclusion, we have demonstrated that iron accumulates in specific brain regions with different patterns. Regions most prone to accumulate iron are the ones involved in the motor, visual, and cognitive functions. Iron accumulation follows a linear pattern during adult lifespan in the basal ganglia and most cortical areas but in the thalamus, pulvinar, precentral and occipito-temporal cortices, it follows a quadratic or exponential pattern. Changes in iron content during aging are different in the pulvinar and the rest of the thalamus. We suggest that whole brain analysis of magnetic susceptibility and segmentation of thalamic and pallidal subregions is beneficial for proper assessment of aging and disease-related changes in cerebral iron content. Departure from normal patterns of regional brain iron trajectories during aging might be helpful in the detection of subtle neurodegenerative and neuroinflammatory processes.


Acknowledgments

This work was presented as a poster at ECTRIMS virtual congress 2020.

Funding: This study was supported by the Czech Ministry of Health grant (NV18-08-00062, and RVO VFN64165), by the Charles University in Prague grant (PROGRES Q27) and by Roche company.


Footnote

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/qims-21-87). Dr. PD reports funding from Czech Ministry of Health, grant No. NU21-04-00535 and European Union’s Horizon 2020 research and innovation programme, grant No. 633190, and advisory board payment from Alexion Pharmaceuticals. Dr. MV reports funding from Czech Ministry of Health, grant No. NV 18-04-00168, honoraria for lectures and presentation from Biogen Idec, Novartis, Sanofi Genzyme, Merck Serono and Teva and support for attending meetings and/or travel from Biogen Idec, Novartis, Sanofi Genzyme. Dr. DH reports honoraria for lectures and presentation from Biogen Idec, Novartis, Sanofi, Roche and Merck and support for attending meetings and/or travel from Biogen Idec, Novartis, Sanofi, Roche and Merck. Dr. JK reports funding from Czech Ministry of Health, grant No. NV 18-04-00168 and support for attending meetings and/or travel Biogen Idec, Novartis, Sanofi Genzyme. Dr. LL serves as an unpaid editorial board member of Quantitative Imaging in Medicine and Surgery. The other authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). This study was approved by the Ethics Committee of the General University Hospital in Prague (ID1018/17), and all subjects signed informed consent.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Callaghan MF, Freund P, Draganski B, Anderson E, Cappelletti M, Chowdhury R, Diedrichsen J, FitzGerald THB, Smittenaar P, Helms G, Lutti A, Weiskopf N. Widespread age-related differences in the human brain microstructure revealed by quantitative magnetic resonance imaging. Neurobiol Aging 2014;35:1862-72. [Crossref] [PubMed]
  2. Lorio S, Lutti A, Kherif F, Ruef A, Dukart J, Chowdhury R, Frackowiak RS, Ashburner J, Helms G, Weiskopf N, Draganski B. Disentangling in vivo the effects of iron content and atrophy on the ageing human brain. NeuroImage 2014;103:280-9. [Crossref] [PubMed]
  3. Pfefferbaum A, Adalsteinsson E, Rohlfing T, Sullivan EV. MRI estimates of brain iron concentration in normal aging: Comparison of field-dependent (FDRI) and phase (SWI) methods. NeuroImage 2009;47:493-500. [Crossref] [PubMed]
  4. Taubert M, Roggenhofer E, Melie-Garcia L, Muller S, Lehmann N, Preisig M, Vollenweider P, Marques-Vidal P, Lutti A, Kherif F, Draganski B. Converging patterns of aging-associated brain volume loss and tissue microstructure differences. Neurobiol Aging 2020;88:108-18. [Crossref] [PubMed]
  5. Wang Y, Liu T. Quantitative susceptibility mapping (QSM): Decoding MRI data for a tissue magnetic biomarker. Magn Reson Med 2015;73:82-101. [Crossref] [PubMed]
  6. Langkammer C, Bredies K, Poser BA, Barth M, Reishofer G, Fan AP, Bilgic B, Fazekas F, Mainero C, Ropele S. Fast quantitative susceptibility mapping using 3D EPI and total generalized variation. NeuroImage 2015;111:622-30. [Crossref] [PubMed]
  7. Acosta-Cabronero J, Milovic C, Mattern H, Tejos C, Speck O, Callaghan MF. A robust multi-scale approach to quantitative susceptibility mapping. NeuroImage 2018;183:7-24. [Crossref] [PubMed]
  8. Dusek P, Dezortova M, Wuerfel J. Imaging of iron. Int Rev Neurobiol 2013;110:195-239. [Crossref] [PubMed]
  9. Burgetova A, Seidl Z, Krasensky J, Horakova D, Vaneckova M. Multiple sclerosis and the accumulation of iron in the Basal Ganglia: quantitative assessment of brain iron using MRI t(2) relaxometry. Eur Neurol 2010;63:136-43. [Crossref] [PubMed]
  10. Connor JR, Menzies SL, Martin SMS, Mufson EJ. Cellular distribution of transferrin, ferritin, and iron in normal and aged human brains. J Neurosci Res 1990;27:595-611. [Crossref] [PubMed]
  11. Daugherty A, Raz N. Age-related differences in iron content of subcortical nuclei observed in vivo: a meta-analysis. NeuroImage 2013;70:113-21. [Crossref] [PubMed]
  12. Schipper HM. Neurodegeneration with brain iron accumulation - clinical syndromes and neuroimaging. Biochim Biophys Acta 2012;1822:350-60. [Crossref] [PubMed]
  13. Ward RJ, Zucca FA, Duyn JH, Crichton RR, Zecca L. The role of iron in brain ageing and neurodegenerative disorders. Lancet Neurol 2014;13:1045-60. [Crossref] [PubMed]
  14. Ghadery C, Pirpamer L, Hofer E, Langkammer C, Petrovic K, Loitfelder M, Schwingenschuh P, Seiler S, Duering M, Jouvent E, Schmidt H, Fazekas F, Mangin J-F, Chabriat H, Dichgans M, Ropele S, Schmidt R. R2* mapping for brain iron: associations with cognition in normal aging. Neurobiol Aging 2015;36:925-32. [Crossref] [PubMed]
  15. Hallgren B, Sourander P. The effect of age on the non-haemin iron in the human brain. J Neurochem 1958;3:41-51. [Crossref] [PubMed]
  16. Ramos P, Santos A, Pinto NR, Mendes R, Magalhães T, Almeida A. Iron levels in the human brain: a post-mortem study of anatomical region differences and age-related changes. J Trace Elem Med Biol 2014;28:13-7. [Crossref] [PubMed]
  17. Li J, Cao F, Yin H, Huang Z, Lin Z, Mao N, Sun B, Wang G. Ferroptosis: past, present and future. Cell Death Dis 2020;11:88. [Crossref] [PubMed]
  18. Schipper HM. Brain iron deposition and the free radical-mitochondrial theory of ageing. Ageing Res Rev 2004;3:265-301. [Crossref] [PubMed]
  19. Daugherty AM, Haacke EM, Raz N. Striatal Iron Content Predicts Its Shrinkage and Changes in Verbal Working Memory after Two Years in Healthy Adults. J Neurosci 2015;35:6731-43. [Crossref] [PubMed]
  20. Kalpouzos G, Garzón B, Sitnikov R, Heiland C, Salami A, Persson J, Bäckman L. Higher Striatal Iron Concentration is Linked to Frontostriatal Underactivation and Poorer Memory in Normal Aging. Cereb Cortex 2017;27:3427-36. [Crossref] [PubMed]
  21. Aquino D, Bizzi A, Grisoli M, Garavaglia B, Bruzzone MG, Nardocci N, Savoiardo M, Chiapparini L. Age-related Iron Deposition in the Basal Ganglia: Quantitative Analysis in Healthy Subjects. Radiology 2009;252:165-72. [Crossref] [PubMed]
  22. Bilgic B, Pfefferbaum A, Rohlfing T, Sullivan EV, Adalsteinsson E. MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping. NeuroImage 2012;59:2625-35. [Crossref] [PubMed]
  23. Gong NJ, Wong CS, Hui ES, Chan CC, Leung LM. Hemisphere, gender and age-related effects on iron deposition in deep gray matter revealed by quantitative susceptibility mapping. NMR Biomed 2015;28:1267-74. [Crossref] [PubMed]
  24. Haacke EM, Miao Y, Liu M, Habib CA, Katkuri Y, Liu T, Yang Z, Lang Z, Hu J, Wu J. Correlation of putative iron content as represented by changes in R2* and phase with age in deep gray matter of healthy adults. J Magn Reson Imaging 2010;32:561-76. [Crossref] [PubMed]
  25. Li W, Wu B, Batrachenko A, Bancroft‐Wu V, Morey RA, Shashi V, Langkammer C, Bellis MDD, Ropele S, Song AW, Liu C. Differential developmental trajectories of magnetic susceptibility in human brain gray and white matter over the lifespan. Hum Brain Mapp 2014;35:2698-713. [Crossref] [PubMed]
  26. Liu M, Liu S, Ghassaban K, Zheng W, Dicicco D, Miao Y, Habib C, Jazmati T, Haacke EM. Assessing global and regional iron content in deep gray matter as a function of age using susceptibility mapping. J Magn Reson Imaging 2016;44:59-71. [Crossref] [PubMed]
  27. Persson N, Wu J, Zhang Q, Liu T, Shen J, Bao R, Ni M, Liu T, Wang Y, Spincemaille P. Age and Sex Related Differences in Subcortical Brain Iron Concentrations among Healthy Adults. NeuroImage 2015;122:385-98. [Crossref] [PubMed]
  28. Acosta-Cabronero J, Betts MJ, Cardenas-Blanco A, Yang S, Nestor PJ. In Vivo MRI Mapping of Brain Iron Deposition across the Adult Lifespan. J Neurosci 2016;36:364-74. [Crossref] [PubMed]
  29. Betts MJ, Acosta-Cabronero J, Cardenas-Blanco A, Nestor PJ, Düzel E. High-resolution characterisation of the aging brain using simultaneous quantitative susceptibility mapping (QSM) and R2* measurements at 7T. NeuroImage 2016;138:43-63. [Crossref] [PubMed]
  30. Milovic C, Prieto C, Bilgic B, Uribe S, Acosta-Cabronero J, Irarrazaval P, Tejos C. Comparison of parameter optimization methods for quantitative susceptibility mapping. Magn Reson Med 2021;85:480-94. [Crossref] [PubMed]
  31. Huhdanpaa H, Hwang DH, Gasparian GG, Booker MT, Cen Y, Lerner A, Boyko OB, Go JL, Kim PE, Rajamohan A, Law M, Shiroishi MS. Image Coregistration: Quantitative Processing Framework for the Assessment of Brain Lesions. J Digit Imaging 2014;27:369-79. [Crossref] [PubMed]
  32. Li X, Chen L, Kutten K, Ceritoglu C, Li Y, Kang N, Hsu J, Qiao Y, Wei H, Liu C, Miller M, Mori S, Yousem DM, van Zijl P, Faria AV. Multi-Atlas Tool for Automated Segmentation of Brain Gray Matter Nuclei and Quantification of their Magnetic Susceptibility. NeuroImage 2019;191:337-49. [Crossref] [PubMed]
  33. Mori S, Wu D, Ceritoglu C, Li Y, Kolasny A, Vaillant MA, Faria AV, Oishi K, Miller MI. MRICloud: Delivering High-Throughput MRI Neuroinformatics as Cloud-Based Software as a Service. Comput Sci Eng 2016;18:21-35. [Crossref]
  34. Avants BB, Tustison NJ, Song G, Gee JC. Ants: Open-source tools for normalization and neuroanatomy. HeanetIe 2009;10:1-11.
  35. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, Gee JC, Gerig G. User-guided 3D active contour segmentation of anatomical structures: significantly improved efficiency and reliability. NeuroImage 2006;31:1116-28. [Crossref] [PubMed]
  36. Joliot M, Jobard G, Naveau M, Delcroix N, Petit L, Zago L, Crivello F, Mellet E, Mazoyer B, Tzourio-Mazoyer N. AICHA: An atlas of intrinsic connectivity of homotopic areas. J Neurosci Methods 2015;254:46-59. [Crossref] [PubMed]
  37. van der Weijden MCM, van Laar PJ, Lambrechts RA, Verbeek DS, Tijssen MAJ. Cortical pencil lining on SWI MRI in NBIA and healthy aging. BMC Neurol 2019;19:233. [Crossref] [PubMed]
  38. Liu C, Wei H, Gong NJ, Cronin M, Dibb R, Decker K. Quantitative Susceptibility Mapping: Contrast Mechanisms and Clinical Applications. Tomography 2015;1:3-17. [Crossref] [PubMed]
  39. Ning N, Liu C, Wu P, Hu Y, Zhang W, Zhang L, Li M, Gho SM, Kim DH, Guo H, Yang J, Jin C. Spatiotemporal variations of magnetic susceptibility in the deep gray matter nuclei from 1 month to 6 years: A quantitative susceptibility mapping study. J Magn Reson Imaging 2019;49:1600-9. [Crossref] [PubMed]
  40. Peterson ET, Kwon D, Luna B, Larsen B, Prouty D, Bellis MDD, Voyvodic J, Liu C, Li W, Pohl KM, Sullivan EV, Pfefferbaum A. Distribution of brain iron accrual in adolescence: Evidence from cross-sectional and longitudinal analysis. Hum Brain Mapp 2019;40:1480-95. [Crossref] [PubMed]
  41. Hagemeier J, Dwyer MG, Bergsland N, Schweser F, Magnano CR, Heininen-Brown M, Ramasamy DP, Carl E, Kennedy C, Melia R, Polak P, Deistung A, Geurts JJG, Reichenbach JR, Zivadinov R. Effect of age on MRI phase behavior in the subcortical deep gray matter of healthy individuals. Am J Neuroradiol 2013;34:2144-51. [Crossref] [PubMed]
  42. Zhang Y, Wei H, Cronin MJ, He N, Yan F, Liu C. Longitudinal atlas for normative human brain development and aging over the lifespan using quantitative susceptibility mapping. NeuroImage 2018;171:176-89. [Crossref] [PubMed]
  43. Burgetova A, Dusek P, Vaneckova M, Horakova D, Langkammer C, Krasensky J, Sobisek L, Matras P, Masek M, Seidl Z. Thalamic Iron Differentiates Primary-Progressive and Relapsing-Remitting Multiple Sclerosis. Am J Neuroradiol 2017;38:1079-86. [Crossref] [PubMed]
  44. Steiger TK, Weiskopf N, Bunzeck N. Iron Level and Myelin Content in the Ventral Striatum Predict Memory Performance in the Aging Brain. J Neurosci 2016;36:3552-8. [Crossref] [PubMed]
  45. Branson HM. Normal Myelination: A Practical Pictorial Review. Neuroimaging Clin N Am 2013;23:183-95. [Crossref] [PubMed]
  46. Dezortova M, Lescinskij A, Dusek P, Herynek V, Acosta‐Cabronero J, Bruha R, Jiru F, Robinson SD, Hajek M. Multiparametric Quantitative Brain MRI in Neurological and Hepatic Forms of Wilson’s Disease. J Magn Reson Imaging 2020;51:1829-35. [Crossref] [PubMed]
  47. Pudlac A, Burgetova A, Dusek P, Nytrova P, Vaneckova M, Horakova D, Krasensky J, Lambert L. Deep Gray Matter Iron Content in Neuromyelitis Optica and Multiple Sclerosis. Biomed Res Int 2020;2020:6492786 [Crossref] [PubMed]
  48. Schweser F, Hagemeier J, Dwyer MG, Bergsland N, Hametner S, Weinstock-Guttman B, Zivadinov R. Decreasing brain iron in multiple sclerosis: The difference between concentration and content in iron MRI. Hum Brain Mapp 2021;42:1463-74. [Crossref] [PubMed]
Cite this article as: Burgetova R, Dusek P, Burgetova A, Pudlac A, Vaneckova M, Horakova D, Krasensky J, Varga Z, Lambert L. Age-related magnetic susceptibility changes in deep grey matter and cerebral cortex of normal young and middle-aged adults depicted by whole brain analysis. Quant Imaging Med Surg 2021;11(9):3906-3919. doi: 10.21037/qims-21-87

Download Citation