MRI estimates of brain iron concentration in normal aging using quantitative susceptibility mapping

Quantifying tissue iron concentration in vivo is instrumental for understanding the role of iron in physiology and in neurological diseases associated with abnormal iron distribution. Herein, we use recently-developed Quantitative Susceptibility Mapping (QSM) methodology to estimate the tissue magnetic susceptibility based on MRI signal phase. To investigate the effect of different regularization choices, we implement and compare ℓ[subscript 1] and ℓ[subscript 2] norm regularized QSM algorithms. These regularized approaches solve for the underlying magnetic susceptibility distribution, a sensitive measure of the tissue iron concentration, that gives rise to the observed signal phase. Regularized QSM methodology also involves a pre-processing step that removes, by dipole fitting, unwanted background phase effects due to bulk susceptibility variations between air and tissue and requires data acquisition only at a single field strength. For validation, performances of the two QSM methods were measured against published estimates of regional brain iron from postmortem and in vivo data. The in vivo comparison was based on data previously acquired using Field-Dependent Relaxation Rate Increase (FDRI), an estimate of MRI relaxivity enhancement due to increased main magnetic field strength, requiring data acquired at two different field strengths. The QSM analysis was based on susceptibility-weighted images acquired at 1.5 T, whereas FDRI analysis used Multi-Shot Echo-Planar Spin Echo images collected at 1.5 T and 3.0 T. Both datasets were collected in the same healthy young and elderly adults. The in vivo estimates of regional iron concentration comported well with published postmortem measurements; both QSM approaches yielded the same rank ordering of iron concentration by brain structure, with the lowest in white matter and the highest in globus pallidus. Further validation was provided by comparison of the in vivo measurements, ℓ[subscript 1]-regularized QSM versus FDRI and ℓ[subscript 2]-regularized QSM versus FDRI, which again yielded perfect rank ordering of iron by brain structure. The final means of validation was to assess how well each in vivo method detected known age-related differences in regional iron concentrations measured in the same young and elderly healthy adults. Both QSM methods and FDRI were consistent in identifying higher iron concentrations in striatal and brain stem ROIs (i.e., caudate nucleus, putamen, globus pallidus, red nucleus, and substantia nigra) in the older than in the young group. The two QSM methods appeared more sensitive in detecting age differences in brain stem structures as they revealed differences of much higher statistical significance between the young and elderly groups than did FDRI. However, QSM values are influenced by factors such as the myelin content, whereas FDRI is a more specific indicator of iron content. Hence, FDRI demonstrated higher specificity to iron yet yielded noisier data despite longer scan times and lower spatial resolution than QSM. The robustness, practicality, and demonstrated ability of predicting the change in iron deposition in adult aging suggest that regularized QSM algorithms using single-field-strength data are possible alternatives to tissue iron estimation requiring two field strengths.National Institutes of Health (U.S.) (Grant NIH R01 EB007942)National Institutes of Health (U.S.) (Grant AG019717)National Institutes of Health (U.S.) (Grant AA005965)National Institutes of Health (U.S.) (Grant AA017168)National Institutes of Health (U.S.) (Grant EB008381)National Science Foundation (U.S.) (Grant 0643836)Siemens CorporationSiemens-MIT AllianceMIT-Center for Integration of Medicine and Innovative Technology (Medical Engineering Fellowship

A new discrete dipole kernel for quantitative susceptibility mapping

PURPOSE: Most approaches for quantitative susceptibility mapping (QSM) are based on a forward model approximation that employs a continuous Fourier transform operator to solve a differential equation system. Such formulation, however, is prone to high-frequency aliasing. The aim of this study was to reduce such errors using an alternative dipole kernel formulation based on the discrete Fourier transform and discrete operators. METHODS: The impact of such an approach on forward model calculation and susceptibility inversion was evaluated in contrast to the continuous formulation both with synthetic phantoms and in vivo MRI data. RESULTS: The discrete kernel demonstrated systematically better fits to analytic field solutions, and showed less over-oscillations and aliasing artifacts while preserving low- and medium-frequency responses relative to those obtained with the continuous kernel. In the context of QSM estimation, the use of the proposed discrete kernel resulted in error reduction and increased sharpness. CONCLUSION: This proof-of-concept study demonstrated that discretizing the dipole kernel is advantageous for QSM. The impact on small or narrow structures such as the venous vasculature might by particularly relevant to high-resolution QSM applications with ultra-high field MRI - a topic for future investigations. The proposed dipole kernel has a straightforward implementation to existing QSM routines

To mask or not to mask? Improving QSM quality by accounting for spatial frequency distributions and susceptibility sources

Estimating magnetic susceptibility using MRI depends on inverting a forward relationship between the susceptibility and measured Larmor frequency. However, an often-overlooked constraint in susceptibility fitting is that the Larmor frequency is only measured inside the sample, and after background field removal, susceptibility sources should only reside inside the same sample. Here we test the impact of accounting for such effects in susceptibility fitting and demonstrate that such effects should not be ignored.Comment: 22 pages, 5 figure

Learning-based Single-step Quantitative Susceptibility Mapping Reconstruction Without Brain Extraction

Quantitative susceptibility mapping (QSM) estimates the underlying tissue magnetic susceptibility from MRI gradient-echo phase signal and typically requires several processing steps. These steps involve phase unwrapping, brain volume extraction, background phase removal and solving an ill-posed inverse problem. The resulting susceptibility map is known to suffer from inaccuracy near the edges of the brain tissues, in part due to imperfect brain extraction, edge erosion of the brain tissue and the lack of phase measurement outside the brain. This inaccuracy has thus hindered the application of QSM for measuring the susceptibility of tissues near the brain edges, e.g., quantifying cortical layers and generating superficial venography. To address these challenges, we propose a learning-based QSM reconstruction method that directly estimates the magnetic susceptibility from total phase images without the need for brain extraction and background phase removal, referred to as autoQSM. The neural network has a modified U-net structure and is trained using QSM maps computed by a two-step QSM method. 209 healthy subjects with ages ranging from 11 to 82 years were employed for patch-wise network training. The network was validated on data dissimilar to the training data, e.g. in vivo mouse brain data and brains with lesions, which suggests that the network has generalized and learned the underlying mathematical relationship between magnetic field perturbation and magnetic susceptibility. AutoQSM was able to recover magnetic susceptibility of anatomical structures near the edges of the brain including the veins covering the cortical surface, spinal cord and nerve tracts near the mouse brain boundaries. The advantages of high-quality maps, no need for brain volume extraction and high reconstruction speed demonstrate its potential for future applications.Comment: 26 page

Characterizing aging in the human brainstem using quantitative multimodal MRI analysis.

Aging is ubiquitous to the human condition. The MRI correlates of healthy aging have been extensively investigated using a range of modalities, including volumetric MRI, quantitative MRI (qMRI), and diffusion tensor imaging. Despite this, the reported brainstem related changes remain sparse. This is, in part, due to the technical and methodological limitations in quantitatively assessing and statistically analyzing this region. By utilizing a new method of brainstem segmentation, a large cohort of 100 healthy adults were assessed in this study for the effects of aging within the human brainstem in vivo. Using qMRI, tensor-based morphometry (TBM), and voxel-based quantification (VBQ), the volumetric and quantitative changes across healthy adults between 19 and 75 years were characterized. In addition to the increased R2* in substantia nigra corresponding to increasing iron deposition with age, several novel findings were reported in the current study. These include selective volumetric loss of the brachium conjunctivum, with a corresponding decrease in magnetization transfer and increase in proton density (PD), accounting for the previously described “midbrain shrinkage.” Additionally, we found increases in R1 and PD in several pontine and medullary structures. We consider these changes in the context of well-characterized, functional age-related changes, and propose potential biophysical mechanisms. This study provides detailed quantitative analysis of the internal architecture of the brainstem and provides a baseline for further studies of neurodegenerative diseases that are characterized by early, pre-clinical involvement of the brainstem, such as Parkinson’s and Alzheimer’s diseases

A post-mortem study of the anatomical region differences and age-related changes on Ca and Mg levels in the human brain

Calcium and magnesium levels in 14 different areas of the human brain [frontal cortex, superior and middle temporal gyri, caudate nucleus, putamen, globus pallidus, cingulated gyrus, hippocampus, inferior parietal lobule, visual cortex of the occipital lobe, midbrain, pons-locus coeruleus, medulla and cerebellum-dentate nucleus] of adult individuals (n = 42; 71 +/- 12, range: 50-103 years old) without a known history of neurodegenerative, neurological or psychiatric disorder were studied. Considering the mean values for the 14 regions, Mg was present at ca. 2.3-fold higher levels than Ca (mean +/-: 527 +/- 34 mu g/g versus 226 +/- 53 mu g/g). Calcium distribution within the brain showed to be quite heterogeneous: highest levels were found in the occipital (306 +/- 156 mu g/g) and frontal cortex (287 +/- 78 ug/g), while lowest levels were found in the medulla (186 +/- 70 mu g/g) and cerebellum (145 +/- 42 mu g/g). Higher Ca levels were found in women than in men (248 +/- 59 mu g/g versus 213 +/- 46 mu g/g; p < 0.05). A tendency for Ca levels to increase with age in all studied brain regions and in both genders was also observed. On the contrary, Mg presented a highly homogeneous distribution and seems to remain quite unchanged irrespective of aging.The authors thank to Universidade do Porto and Santander Totta for financial support through the project "TRAIN: Trace elements in human brain: age-related changes and anatomic region specific differences" (PP_IJUP 2011 342)