χ-separation: Magnetic susceptibility source separation toward iron and myelin mapping in the brain

Assessing brain integrity in patients with long-term and well-functioning metal-based hip implants

Production of metal debris from implant wear and corrosion processes is now a well understood occurrence following hip arthroplasty. Evidence has shown that metal ions can enter the bloodstream and travel to distant organs including the brain, and in extreme cases, can induce sensorial and neurological diseases. Our objective was tosimultaneously analyze brain anatomy and physiology in patients with long-term and well-functioning implants. Included were subjects who had received total hip or hip resurfacing arthroplastywith an implantation time of a minimum of 7 years (n = 28) and age- and sex-matched controls (n = 32). Blood samples were obtained to measure ion concentrations of cobalt and chromium, and the Montreal Cognitive Assessment was performed. 3T MRI brain scans were completed with an MPRAGE sequence for ROI segmentation and multiecho gradient echo sequences to generate QSM and R2* maps. Mean QSM and R2* values were recorded for five deep brain and four middle and cortical brain structures on both hemispheres: pallidum, putamen, caudate, amygdala, hippocampus, anterior cingulate, inferior temporal, and cerebellum. No differences in QSM or R2* or cognition scores were found between both groups (p > 0.6654). No correlation was found between susceptibility and blood ion levels for cobalt or chromium in any region of the brain. No correlation was found between blood ion levels and cognition scores. Clinical significance: Results suggest that metal ions released by long-term and well-functioning implants do not affect brain integrity.

Value of Quantitative Susceptibility Mapping in Clinical Neuroradiology

Quantitative susceptibility mapping (QSM) is a unique technique for providing quantitative information on tissue magnetic susceptibility using phase image data. QSM can provide valuable information regarding physiological and pathological processes such as iron deposition, hemorrhage, calcification, and myelin. QSM has been considered for use as an imaging biomarker to investigate physiological status and pathological changes. Although various studies have investigated the clinical applications of QSM, particularly regarding the use of QSM in clinical practice, have not been examined well. This review provides on an overview of the basics of QSM and its clinical applications in neuroradiology. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Iron Deposition and Distribution Across the Hippocampus Is Associated with Pattern Separation and Pattern Completion in Older Adults at Risk for Alzheimer's Disease

Elevated iron deposition in the brain has been observed in older adult humans and persons with Alzheimer's disease (AD), and has been associated with lower cognitive performance. We investigated the impact of iron deposition, and its topographical distribution across hippocampal subfields and segments (anterior, posterior) measured along its longitudinal axis, on episodic memory in a sample of cognitively unimpaired older adults at elevated familial risk for AD (N = 172, 120 females, 52 males; mean age = 68.8 ± 5.4 years). MRI-based quantitative susceptibility maps were acquired to derive estimates of hippocampal iron deposition. The Mnemonic Similarity Task was used to measure pattern separation and pattern completion, two hippocampally mediated episodic memory processes. Greater hippocampal iron load was associated with lower pattern separation and higher pattern completion scores, both indicators of poorer episodic memory. Examination of iron levels within hippocampal subfields across its long axis revealed topographic specificity. Among the subfields and segments investigated here, iron deposition in the posterior hippocampal CA1 was the most robustly and negatively associated with the fidelity memory representations. This association remained after controlling for hippocampal volume and was observed in the context of normal performance on standard neuropsychological memory measures. These findings reveal that the impact of iron load on episodic memory performance is not uniform across the hippocampus. Both iron deposition levels as well as its spatial distribution, must be taken into account when examining the relationship between hippocampal iron and episodic memory in older adults at elevated risk for AD.

Comparison between R2'-based and R2*-based χ-separation methods: A clinical evaluation in individuals with multiple sclerosis

Susceptibility source separation, or χ-separation, estimates diamagnetic (χdia) and paramagnetic susceptibility (χpara) signals in the brain using local field and R2' (= R2* - R2) maps. Recently proposed R2*-based χ-separation methods allow for χ-separation using only multi-echo gradient echo (ME-GRE) data, eliminating the need for additional data acquisition for R2 mapping. Although this approach reduces scan time and enhances clinical utility, the impact of missing R2 information remains a subject of exploration. In this study, we evaluate the viability of two previously proposed R2*-based χ-separation methods as alternatives to their R2'-based counterparts: model-based R2*-χ-separation versus χ-separation and deep learning-based χ-sepnet-R2* versus χ-sepnet-R2'. Their performances are assessed in individuals with multiple sclerosis (MS), comparing them with their corresponding R2'-based counterparts (i.e., R2*-χ-separation vs. χ-separation and χ-sepnet-R2* vs. χ-sepnet-R2'). The evaluations encompass qualitative visual assessments by experienced neuroradiologists and quantitative analyses, including region of interest analyses and linear regression analyses. Qualitatively, R2*-χ-separation tends to report higher χpara and χdia values compared with χ-separation, leading to less distinct lesion contrasts, while χ-sepnet-R2* closely aligns with χ-sepnet-R2'. Quantitative analysis reveals a robust correlation between both R2*-based methods and their R2'-based counterparts (r ≥ 0.88). Specifically, in the whole-brain voxels, χ-sepnet-R2* exhibits higher correlation and better linearity than R2*-χ-separation (χdia/χpara from R2*-χ-separation: r = 0.88/0.90, slope = 0.79/0.86; χdia/χpara from χ-sepnet-R2*: r = 0.90/0.92, slope = 0.99/0.97). In MS lesions, both R2*-based methods display comparable correlation and linearity (χdia/χpara from R2*-χ-separation: r = 0.90/0.91, slope = 0.98/0.91; χdia/χpara from χ-sepnet-R2*: r = 0.88/0.88, slope = 0.91/0.95). Notably, χ-sepnet-R2* demonstrates negligible offsets, whereas R2*-χ-separation exhibits relatively large offsets (0.02 ppm in the whole brain and 0.01 ppm in the MS lesions), potentially indicating the false presence of myelin or iron in MS lesions. Overall, both R2*-based χ-separation methods demonstrated their viability as alternatives to their R2'-based counterparts. χ-sepnet-R2* showed better alignment with its R2'-based counterpart with minimal susceptibility offsets, compared with R2*-χ-separation that reported higher χpara and χdia values compared with R2'-based χ-separation.

Recommended implementation of quantitative susceptibility mapping for clinical research in the brain: A consensus of the ISMRM electro-magnetic tissue properties study group

This article provides recommendations for implementing QSM for clinical brain research. It is a consensus of the International Society of Magnetic Resonance in Medicine, Electro-Magnetic Tissue Properties Study Group. While QSM technical development continues to advance rapidly, the current QSM methods have been demonstrated to be repeatable and reproducible for generating quantitative tissue magnetic susceptibility maps in the brain. However, the many QSM approaches available have generated a need in the neuroimaging community for guidelines on implementation. This article outlines considerations and implementation recommendations for QSM data acquisition, processing, analysis, and publication. We recommend that data be acquired using a monopolar 3D multi-echo gradient echo (GRE) sequence and that phase images be saved and exported in Digital Imaging and Communications in Medicine (DICOM) format and unwrapped using an exact unwrapping approach. Multi-echo images should be combined before background field removal, and a brain mask created using a brain extraction tool with the incorporation of phase-quality-based masking. Background fields within the brain mask should be removed using a technique based on SHARP or PDF, and the optimization approach to dipole inversion should be employed with a sparsity-based regularization. Susceptibility values should be measured relative to a specified reference, including the common reference region of the whole brain as a region of interest in the analysis. The minimum acquisition and processing details required when reporting QSM results are also provided. These recommendations should facilitate clinical QSM research and promote harmonized data acquisition, analysis, and reporting.

Revealing the MRI-Contrast in Optically Cleared Brains

The current consensus holds that optically-cleared specimens are unsuitable for Magnetic Resonance Imaging (MRI); exhibiting absence of contrast. Prior studies combined MRI with tissue-clearing techniques relying on the latter's ability to eliminate lipids, thereby fostering the assumption that lipids constitute the primary source of ex vivo MRI-contrast. Nevertheless, these findings contradict an extensive body of literature that underscores the contribution of other features to contrast. Furthermore, it remains unknown whether non-delipidating clearing methods can produce MRI-compatible specimens or whether MRI-contrast can be re-established. These limitations hinder the development of multimodal MRI-light-microscopy (LM) imaging approaches. This study assesses the relation between MRI-contrast, and delipidation in optically-cleared whole brains following different tissue-clearing approaches. It is demonstrated that uDISCO and ECi-brains are MRI-compatible upon tissue rehydration, despite both methods' substantial delipidating-nature. It is also demonstrated that, whereas Scale-clearing preserves most lipids, Scale-cleared brain lack MRI-contrast. Furthermore, MRI-contrast is restored to lipid-free CLARITY-brains without introducing lipids. Our results thereby dissociate between the essentiality of lipids to MRI-contrast. A tight association is found between tissue expansion, hyperhydration and loss of MRI-contrast. These findings then enabled us to develop a multimodal MRI-LM-imaging approach, opening new avenues to bridge between the micro- and mesoscale for biomedical research and clinical applications.

So You Want to Image Myelin Using MRI: Magnetic Susceptibility Source Separation for Myelin Imaging

In MRI, researchers have long endeavored to effectively visualize myelin distribution in the brain, a pursuit with significant implications for both scientific research and clinical applications. Over time, various methods such as myelin water imaging, magnetization transfer imaging, and relaxometric imaging have been developed, each carrying distinct advantages and limitations. Recently, an innovative technique named as magnetic susceptibility source separation has emerged, introducing a novel surrogate biomarker for myelin in the form of a diamagnetic susceptibility map. This paper comprehensively reviews this cutting-edge method, providing the fundamental concepts of magnetic susceptibility, susceptibility imaging, and the validation of the diamagnetic susceptibility map as a myelin biomarker that indirectly measures myelin content. Additionally, the paper explores essential aspects of data acquisition and processing, offering practical insights for readers. A comparison with established myelin imaging methods is also presented, and both current and prospective clinical and scientific applications are discussed to provide a holistic understanding of the technique. This work aims to serve as a foundational resource for newcomers entering this dynamic and rapidly expanding field.

Vascular Aging in the Choroid Plexus: A 7T Ultrasmall Superparamagnetic Iron Oxide (USPIO)-MRI Study

BACKGROUND

The choroid plexus (ChP), a densely vascularized structure, has drawn increasing attention for its involvement in brain homeostasis and waste clearance. While the volumetric changes have been explored in many imaging studies, few studies have investigated the vascular degeneration associated with aging in the ChP.

PURPOSE

To investigate the sub-structural characteristics of the ChP, particularly the vascular compartment using high-resolution 7T imaging enhanced with Ferumoxytol, an ultrasmall super-paramagnetic iron oxide, which greatly increase the susceptibility contrast for vessels.

STUDY TYPE

Prospective.

SUBJECTS

Forty-nine subjects without neurological disorders (age: 21-80 years; 42 ± 17 years; 20 females).

FIELD STRENGTH/SEQUENCE

7-T with 2D and 3D T2* GRE, 3D MPRAGE T1, 2D TSE T2, and 2D FLAIR.

ASSESSMENT

The vascular and stromal compartments of the ChP were segmented using K-means clustering on post-contrast 2D GRE images. Visual and qualitative assessment of ChP vascular characteristics were conducted independently by three observers. Vascular density (Volvessel/VolChP ratio) and susceptibility change (Δχ) induced by Ferumoxytol were analyzed on 3D GRE-derived susceptibility-weighted imaging and quantitative susceptibility mapping, respectively.

STATISTICAL TESTS

Independent t-test, Mann-Whitney U test, and Chi-square test were utilized for group comparisons. The relationship between age and ChP's vascular alterations was examined using Pearson's correlation. Intra-class coefficient was calculated for inter-observer agreement. A P value <0.05 was considered statistically significant.

RESULTS

2D GRE images demonstrated superior contrast and accurate delineation of ChP substructures (ICC = 0.86). Older subjects exhibited a significantly smaller vascular density (16.5 ± 4.34%) and lower Δχ (22.10 ± 12.82 ppb) compared to younger subjects (24.85 ± 6.84% and 34.64 ± 12.69 ppb). Vascular density and mean Δχ within the ChP negatively correlated with age (r = -0.48, and r = -0.45).

DATA CONCLUSION

Ferumoxytol-enhanced 7T images can demonstrate ChP alterations in elderly with decreased vascular density and expansion of nonvascular compartment.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

Evaluation of brain iron deposition in different cerebral arteries of acute ischaemic stroke patients using quantitative susceptibility mapping

AIM

To investigate differences in iron deposition between infarct and normal cerebral arterial regions in acute ischaemic stroke (AIS) patients using quantitative susceptibility mapping (QSM).

MATERIALS AND METHODS

Forty healthy controls and 40 AIS patients were recruited, and their QSM images were obtained. There were seven regions of interest (ROIs) in AIS patients, including the infarct regions of responsible arteries (R1), the non-infarct regions of responsible arteries (R2), the contralateral symmetrical sites of lesions (R3), and the non-responsible cerebral arterial regions (R4, R5, R6, R7). For the healthy controls, the cerebral arterial regions corresponding to the AIS patient group were selected as ROIs. The differences in corresponding ROI susceptibilities between AIS patients and healthy controls and the differences in susceptibilities between infarcted and non-infarct regions in AIS patients were compared.

RESULTS

The susceptibilities of infarct regions in AIS patients were significantly higher than those in healthy controls (p<0.0001). There was no significant difference in non-infarct regions between the two groups (p>0.05). The susceptibility of the infarct regions in AIS patients was significantly higher than those of the non-infarct region of responsible artery and non-responsible cerebral arterial regions (p<0.01).

CONCLUSIONS

Abnormal iron deposition detected by QSM in the infarct regions of AIS patients may not affect iron levels in the non-infarct regions of responsible arteries and normal cerebral arteries, which may open the door for potential new diagnostic and treatment strategies.

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.

MR Susceptibility Separation for Quantifying Lesion Paramagnetic and Diamagnetic Evolution in Relapsing-Remitting Multiple Sclerosis

BACKGROUND

Multiple sclerosis (MS) lesion evolution may involve changes in diamagnetic myelin and paramagnetic iron. Conventional quantitative susceptibility mapping (QSM) can provide net susceptibility distribution, but not the discrete paramagnetic and diamagnetic components.

PURPOSE

To apply susceptibility separation (χ separation) to follow lesion evolution in MS with comparison to R2 */R2 ' /QSM.

STUDY TYPE

Longitudinal, prospective.

SUBJECTS

Twenty relapsing-remitting MS subjects (mean age: 42.5 ± 9.4 years, 13 females; mean years of symptoms: 4.3 ± 1.4 years).

FIELD STRENGTH/SEQUENCE

Three-dimensional multiple echo gradient echo (QSM and R2 * mapping), two-dimensional dual echo fast spin echo (R2 mapping), T2 -weighted fluid attenuated inversion recovery, and T1-weighted magnetization prepared gradient echo sequences at 3 T.

ASSESSMENT

Data were analyzed from two scans separated by a mean interval of 14.4 ± 2.0 months. White matter lesions on fluid-attenuated inversion recovery were defined by an automatic pipeline, then manually refined (by ZZ/AHW, 3/25 years' experience in MRI), and verified by a radiologist (MN, 25 years' experience in MS). Susceptibility separation yielded the paramagnetic and diamagnetic susceptibility content of each voxel. Lesions were classified into four groups based on the variation of QSM/R2 * or separated into positive/negative components from χ separation.

STATISTICAL TESTS

Two-sample paired t tests for assessment of longitudinal differences. Spearman correlation coefficients to assess associations between χ separation and R2 */R2 ' /QSM. Significant level: P < 0.005.

RESULTS

A total of 183 lesions were quantified. Categorizing lesions into groups based on χ separation demonstrated significant annual changes in QSM//R2 */R2 ' . When lesions were grouped based on changes in QSM and R2 *, both changing in unison yielded a significant dominant paramagnetic variation and both opposing yielded a dominant diamagnetic variation. Significant Spearman correlation coefficients were found between susceptibility-sensitive MRI indices and χ separation.

DATA CONCLUSION

Susceptibility separation changes in MS lesions may distinguish and quantify paramagnetic and diamagnetic evolution, potentially providing additional insight compared to R2 * and QSM alone.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Incorporating the effect of white matter microstructure in the estimation of magnetic susceptibility in ex vivo mouse brain

PURPOSE

To extend quantitative susceptibility mapping to account for microstructure of white matter (WM) and demonstrate its effect on ex vivo mouse brain at 16.4T.

THEORY AND METHODS

Previous studies have shown that the MRI measured Larmor frequency also depends on local magnetic microstructure at the mesoscopic scale. Here, we include effects from WM microstructure using our previous results for the mesoscopic Larmor frequencyΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ of cylinders with arbitrary orientations. We scrutinize the validity of our model and QSM in a digital brain phantom includingΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from a WM susceptibility tensor and biologically stored iron with scalar susceptibility. We also apply susceptibility tensor imaging to the phantom and investigate how the fitted tensors are biased fromΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ . Last, we demonstrate how to combine multi-gradient echo and diffusion MRI images of ex vivo mouse brains acquired at 16.4T to estimate an apparent scalar susceptibility without sample rotations.

RESULTS

Our new model improves susceptibility estimation compared to QSM for the brain phantom. Applying susceptibility tensor imaging to the phantom withΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ from WM axons with scalar susceptibility produces a highly anisotropic susceptibility tensor that mimics results from previous susceptibility tensor imaging studies. For the ex vivo mouse brain we find theΩ ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ due to WM microstructure to be substantial, changing susceptibility in WM up to 25% root-mean-squared-difference.

CONCLUSION

Ω ‾ Meso $$ {\overline{\Omega}}^{\mathrm{Meso}} $$ impacts susceptibility estimates and biases susceptibility tensor imaging fitting substantially. Hence, it should not be neglected when imaging structurally anisotropic tissue such as brain WM.

Depth- and curvature-based quantitative susceptibility mapping analyses of cortical iron in Alzheimer's disease

In addition to amyloid beta plaques and neurofibrillary tangles, Alzheimer's disease (AD) has been associated with elevated iron in deep gray matter nuclei using quantitative susceptibility mapping (QSM). However, only a few studies have examined cortical iron, using more macroscopic approaches that cannot assess layer-specific differences. Here, we conducted column-based QSM analyses to assess whether AD-related increases in cortical iron vary in relation to layer-specific differences in the type and density of neurons. We obtained global and regional measures of positive (iron) and negative (myelin, protein aggregation) susceptibility from 22 adults with AD and 22 demographically matched healthy controls. Depth-wise analyses indicated that global susceptibility increased from the pial surface to the gray/white matter boundary, with a larger slope for positive susceptibility in the left hemisphere for adults with AD than controls. Curvature-based analyses indicated larger global susceptibility for adults with AD versus controls; the right hemisphere versus left; and gyri versus sulci. Region-of-interest analyses identified similar depth- and curvature-specific group differences, especially for temporo-parietal regions. Finding that iron accumulates in a topographically heterogenous manner across the cortical mantle may help explain the profound cognitive deterioration that differentiates AD from the slowing of general motor processes in healthy aging.

Sub-voxel quantitative susceptibility mapping for assessing whole-brain magnetic susceptibility from ages 4 to 80

The evolution of magnetic susceptibility of the brain is mainly determined by myelin in white matter (WM) and iron deposition in deep gray matter (DGM). However, existing imaging techniques have limited abilities to simultaneously quantify the myelination and iron deposition within a voxel throughout brain development and aging. For instance, the temporal trajectories of iron in the brain WM and myelination in DGM have not been investigated during the aging process. This study aimed to map the age-related iron and myelin changes in the whole brain, encompassing myelin in DGM and iron deposition in WM, using a novel sub-voxel quantitative susceptibility mapping (QSM) method. To achieve this, a cohort of 494 healthy adults (18-80 years old) was studied. The sub-voxel QSM method was employed to obtain the paramagnetic and diamagnetic susceptibility based on the approximatedR 2 ' map from acquiredR 2 * map. The linear relationship betweenR 2 * andR 2 ' maps was established from the regression coefficients on a small cohort data acquired with both 3D gradient recalled echo data andR 2 mapping. Large cohort sub-voxel susceptibility maps were used to create longitudinal and age-specific atlases via group-wise registration. To explore the differential developmental trajectories in the DGM and WM, we employed nonlinear models including exponential and Poisson functions, along with generalized additive models. The constructed atlases reveal the iron accumulation in the posterior part of the putamen and the gradual myelination process in the globus pallidus with aging. Interestingly, the developmental trajectories show that the rate of myelination differs among various DGM regions. Furthermore, the process of myelin synthesis is paralleled by an associated pattern of iron accumulation in the primary WM fiber bundles. In summary, our study offers significant insights into the distinctive developmental trajectories of iron in the brain's WM and myelination/demyelination in the DGM in vivo. These findings highlight the potential of using sub-voxel QSM to uncover new perspectives in neuroscience and improve our understanding of whole-brain myelination and iron deposit processes across the lifespan.

Multi-parametric MRI can detect enhanced myelination in the Gli1 -/- mouse brain

Purpose

The ability of MRI-based markers to detect myelin in the brain is limited. This study investigated the potential of combining multiple MRI markers, each targeting distinct myelin properties, to improve myelin characterization.

Methods

We acquired ex vivo multiparametric MRI data at 7 Tesla from control and Gli1 -/- mouse brains at postnatal day 10 (P10), which exhibits enhanced myelination in the corpus callosum, followed by myelin basic protein (MBP) stained immunohistochemistry.

Results

Although most MRI markers included in this study showed significant differences in the corpus callosum between control and Gli1 -/- , only fractional anisotropy (FA), mean diffusivity (MD), and T 2 had strong correlations with MBP signals. Partial least square regression (PSLR) based on MRI and MBP values from white matter regions suggested that T 2 had the highest contributions to myelin estimation. When both white and gray matter regions were included, inhomogeneous MT ratio and FA showed strong contributions.

Conclusion

This study demonstrates the efficacy of multi-parametric MRI in detecting enhanced myelination in the Gli1 -/- mouse brain. T 2 and diffusion MRI parameters showed strong correlation with MBP signals in the genu of the corpus callosum at P10. The contribution of individual MRI parameter for detecting myelin can be evaluated using PLSR.

Quantitative susceptibility mapping of brain iron in healthy aging and cognition

Quantitative susceptibility mapping (QSM) is a magnetic resonance imaging (MRI) technique that can assess the magnetic properties of cerebral iron in vivo. Although brain iron is necessary for basic neurobiological functions, excess iron content disrupts homeostasis, leads to oxidative stress, and ultimately contributes to neurodegenerative disease. However, some degree of elevated brain iron is present even among healthy older adults. To better understand the topographical pattern of iron accumulation and its relation to cognitive aging, we conducted an integrative review of 47 QSM studies of healthy aging, with a focus on five distinct themes. The first two themes focused on age-related increases in iron accumulation in deep gray matter nuclei versus the cortex. The overall level of iron is higher in deep gray matter nuclei than in cortical regions. Deep gray matter nuclei vary with regard to age-related effects, which are most prominent in the putamen, and age-related deposition of iron is also observed in frontal, temporal, and parietal cortical regions during healthy aging. The third theme focused on the behavioral relevance of iron content and indicated that higher iron in both deep gray matter and cortical regions was related to decline in fluid (speed-dependent) cognition. A handful of multimodal studies, reviewed in the fourth theme, suggest that iron interacts with imaging measures of brain function, white matter degradation, and the accumulation of neuropathologies. The final theme concerning modifiers of brain iron pointed to potential roles of cardiovascular, dietary, and genetic factors. Although QSM is a relatively recent tool for assessing cerebral iron accumulation, it has significant promise for contributing new insights into healthy neurocognitive aging.

High angular resolution susceptibility imaging and estimation of fiber orientation distribution functions in primate brain

Uncovering brain-tissue microstructure including axonal characteristics is a major neuroimaging research focus. Within this scope, anisotropic properties of magnetic susceptibility in white matter have been successfully employed to estimate primary axonal trajectories using mono-tensorial models. However, anisotropic susceptibility has not yet been considered for modeling more complex fiber structures within a voxel, such as intersecting bundles, or an estimation of orientation distribution functions (ODFs). This information is routinely obtained by high angular resolution diffusion imaging (HARDI) techniques. In applications to fixed tissue, however, diffusion-weighted imaging suffers from an inherently low signal-to-noise ratio and limited spatial resolution, leading to high demands on the performance of the gradient system in order to mitigate these limitations. In the current work, high angular resolution susceptibility imaging (HARSI) is proposed as a novel, phase-based methodology to estimate ODFs. A multiple gradient-echo dataset was acquired in an entire fixed chimpanzee brain at 61 orientations by reorienting the specimen in the magnetic field. The constant solid angle method was adapted for estimating phase-based ODFs. HARDI data were also acquired for comparison. HARSI yielded information on whole-brain fiber architecture, including identification of peaks of multiple bundles that resembled features of the HARDI results. Distinct differences between both methods suggest that susceptibility properties may offer complementary microstructural information. These proof-of-concept results indicate a potential to study the axonal organization in post-mortem primate and human brain at high resolution.

DeepSTI: Towards tensor reconstruction using fewer orientations in susceptibility tensor imaging

Susceptibility tensor imaging (STI) is an emerging magnetic resonance imaging technique that characterizes the anisotropic tissue magnetic susceptibility with a second-order tensor model. STI has the potential to provide information for both the reconstruction of white matter fiber pathways and detection of myelin changes in the brain at mm resolution or less, which would be of great value for understanding brain structure and function in healthy and diseased brain. However, the application of STI in vivo has been hindered by its cumbersome and time-consuming acquisition requirement of measuring susceptibility induced MR phase changes at multiple head orientations. Usually, sampling at more than six orientations is required to obtain sufficient information for the ill-posed STI dipole inversion. This complexity is enhanced by the limitation in head rotation angles due to physical constraints of the head coil. As a result, STI has not yet been widely applied in human studies in vivo. In this work, we tackle these issues by proposing an image reconstruction algorithm for STI that leverages data-driven priors. Our method, called DeepSTI, learns the data prior implicitly via a deep neural network that approximates the proximal operator of a regularizer function for STI. The dipole inversion problem is then solved iteratively using the learned proximal network. Experimental results using both simulation and in vivo human data demonstrate great improvement over state-of-the-art algorithms in terms of the reconstructed tensor image, principal eigenvector maps and tractography results, while allowing for tensor reconstruction with MR phase measured at much less than six different orientations. Notably, promising reconstruction results are achieved by our method from only one orientation in human in vivo, and we demonstrate a potential application of this technique for estimating lesion susceptibility anisotropy in patients with multiple sclerosis.

QSM Throughout the Body

Magnetic materials in tissue, such as iron, calcium, or collagen, can be studied using quantitative susceptibility mapping (QSM). To date, QSM has been overwhelmingly applied in the brain, but is increasingly utilized outside the brain. QSM relies on the effect of tissue magnetic susceptibility sources on the MR signal phase obtained with gradient echo sequence. However, in the body, the chemical shift of fat present within the region of interest contributes to the MR signal phase as well. Therefore, correcting for the chemical shift effect by means of water-fat separation is essential for body QSM. By employing techniques to compensate for cardiac and respiratory motion artifacts, body QSM has been applied to study liver iron and fibrosis, heart chamber blood and placenta oxygenation, myocardial hemorrhage, atherosclerotic plaque, cartilage, bone, prostate, breast calcification, and kidney stone.

Multimodal comparisons of QSM and PET in neurodegeneration and aging

Quantitative susceptibility mapping (QSM) has been used to study susceptibility changes that may occur based on tissue composition and mineral deposition. Iron is a primary contributor to changes in magnetic susceptibility and of particular interest in applications of QSM to neurodegeneration and aging. Iron can contribute to neurodegeneration through inflammatory processes and via interaction with aggregation of disease-related proteins. To better understand the local susceptibility changes observed on QSM, its signal has been studied in association with other imaging metrics such as positron emission tomography (PET). The associations of QSM and PET may provide insight into the pathophysiology of disease processes, such as the role of iron in aging and neurodegeneration, and help to determine the diagnostic utility of QSM as an indirect indicator of disease processes typically evaluated with PET. In this review we discuss the proposed mechanisms and summarize prior studies of the associations of QSM and amyloid PET, tau PET, TSPO PET, FDG-PET, 15O-PET, and F-DOPA PET in evaluation of neurologic diseases with a focus on aging and neurodegeneration.

APART-QSM: an improved sub-voxel quantitative susceptibility mapping for susceptibility source separation using an iterative data fitting method

The brain tissue phase contrast in MRI sequences reflects the spatial distributions of multiple substances, such as iron, myelin, calcium, and proteins. These substances with paramagnetic and diamagnetic susceptibilities often colocalize in one voxel in brain regions. Both opposing susceptibilities play vital roles in brain development and neurodegenerative diseases. Conventional QSM methods only provide voxel-averaged susceptibility value and cannot disentangle intravoxel susceptibilities with opposite signs. Advanced susceptibility imaging methods have been recently developed to distinguish the contributions of opposing susceptibility sources for QSM. The basic concept of separating paramagnetic and diamagnetic susceptibility proportions is to include the relaxation rate R₂^* with R₂^' in QSM. The magnitude decay kernel, describing the proportionality coefficient between R₂^' and susceptibility, is an essential reconstruction coefficient for QSM separation methods. In this study, we proposed a more comprehensive complex signal model that describes the relationship between 3D GRE signal and the contributions of paramagnetic and diamagnetic susceptibility to the frequency shift and R₂^* relaxation. The algorithm is implemented as a constrained minimization problem in which the voxel-wise magnitude decay kernel and sub-voxel susceptibilities are determined alternately in each iteration until convergence. The calculated voxel-wise magnitude decay kernel could realistically model the relationship between the R₂^' relaxation and the volume susceptibility. Thus, the proposed method effectively prevents the errors of the magnitude decay kernel from propagating to the final susceptibility separation reconstruction. Phantom studies, ex vivo macaque brain experiments, and in vivo human brain imaging studies were conducted to evaluate the ability of the proposed method to distinguish paramagnetic and diamagnetic susceptibility sources. The results demonstrate that the proposed method provides state-of-the-art performances for quantifying brain iron and myelin compared to previous QSM separation methods. Our results show that the proposed method has the potential to simultaneously quantify whole brain iron and myelin during brain development and aging. The proposed model was also deployed with multiple-orientation complex GRE data input measurements, resulting in high-quality QSM separation maps with more faithful tissue delineation between brain structures compared to those reconstructed by single-orientation QSM separation methods.

Simulated and experimental phantom data for multi-center quality assurance of quantitative susceptibility maps at 3 T, 7 T and 9.4 T

PURPOSE

To develop methods for quality assurance of quantitative susceptibility mapping (QSM) using MRI at different magnetic field strengths, and scanners, using different MR-sequence protocols, and post-processing pipelines.

METHODS

We built a custom phantom based on iron in two forms: homogeneous susceptibility ('free iron') and with fine-scaled variations in susceptibility ('clustered iron') at different iron concentrations. The phantom was measured at 3.0 T (two scanners), 7.0 T and 9.4 T using multi-echo, gradient echo acquisition sequences. A digital phantom analogue to the iron-phantom, tailored to obtain similar results as in experimentation was developed, with similar geometry and susceptibility values. Morphology enabled dipole inversion was applied to the phase images to obtain QSM for experimental and simulated data using the MEDI + 0 approach for background regularization.

RESULTS

Across all scanners, QSM-values showed a linear increase with iron concentrations. The QSM-relaxivity was 0.231 ± 0.047 ppm/mM for free and 0.054 ± 0.013 ppm/mM for clustered iron, with adjusted determination coefficients (DoC) ≥ 0.87. Similarly, the simulations yielded linear increases (DoC ≥ 0.99). In both the experimental and digital phantoms, the estimated molar susceptibility was lower with clustered iron, because clustering led to highly localized field effects.

CONCLUSION

Our iron phantom can be used to evaluate the capability of QSM to detect local variations in susceptibility across different field strengths, when using different MR-sequence protocols. The devised simulation method captures the effect of iron clustering in QSM as seen experimentally and could be used in the future to optimize QSM processing pipelines and achieve higher accuracy for local field effects, as also seen in Alzheimer's beta-amyloid plaques.

Effect of motion, cortical orientation and spatial resolution on quantitative imaging of cortical R2* and magnetic susceptibility at 0.3 mm in-plane resolution at 7 T

MR images of the effective relaxation rate R2* and magnetic susceptibility χ derived from multi-echo T2*-weighted (T2*w) MRI can provide insight into iron and myelin distributions in the brain, with the potential of providing biomarkers for neurological disorders. Quantification of R2* and χ at submillimeter resolution in the cortex in vivo has been difficult because of challenges such as head motion, limited signal to noise ratio, long scan time, and motion related magnetic field fluctuations. This work aimed to improve the robustness for quantifying intracortical R2* and χ and analyze the effects from motion, spatial resolution, and cortical orientation. T2*w data was acquired with a spatial resolution of 0.3 × 0.3 × 0.4 mm3 at 7 T and downsampled to various lower resolutions. A combined correction for motion and B0 changes was deployed using volumetric navigators. Such correction improved the T2*w image quality rated by experienced image readers and test-retest reliability of R2* and χ quantification with reduced median inter-scan differences up to 10 s-1 and 5 ppb, respectively. R2* and χ near the line of Gennari, a cortical layer high in iron and myelin, were as much as 10 s-1 and 10 ppb higher than the region at adjacent cortical depth. In addition, a significant effect due to the cortical orientation relative to the static field (B0) was observed in χ with a peak-to-peak amplitude of about 17 ppb. In retrospectively downsampled data, the capability to distinguish different cortical depth regions based on R2* or χ contrast remained up to isotropic 0.5 mm resolution. This study highlights the unique characteristics of R2* and χ along the cortical depth at submillimeter resolution and the need for motion and B0 corrections for their robust quantification in vivo.

Patterns of inflammation, microstructural alterations, and sodium accumulation define multiple sclerosis subtypes after 15 years from onset

Introduction

Conventional MRI is routinely used for the characterization of pathological changes in multiple sclerosis (MS), but due to its lack of specificity is unable to provide accurate prognoses, explain disease heterogeneity and reconcile the gap between observed clinical symptoms and radiological evidence. Quantitative MRI provides measures of physiological abnormalities, otherwise invisible to conventional MRI, that correlate with MS severity. Analyzing quantitative MRI measures through machine learning techniques has been shown to improve the understanding of the underlying disease by better delineating its alteration patterns.

Methods

In this retrospective study, a cohort of healthy controls (HC) and MS patients with different subtypes, followed up 15 years from clinically isolated syndrome (CIS), was analyzed to produce a multi-modal set of quantitative MRI features encompassing relaxometry, microstructure, sodium ion concentration, and tissue volumetry. Random forest classifiers were used to train a model able to discriminate between HC, CIS, relapsing remitting (RR) and secondary progressive (SP) MS patients based on these features and, for each classification task, to identify the relative contribution of each MRI-derived tissue property to the classification task itself.

Results and discussion

Average classification accuracy scores of 99 and 95% were obtained when discriminating HC and CIS vs. SP, respectively; 82 and 83% for HC and CIS vs. RR; 76% for RR vs. SP, and 79% for HC vs. CIS. Different patterns of alterations were observed for each classification task, offering key insights in the understanding of MS phenotypes pathophysiology: atrophy and relaxometry emerged particularly in the classification of HC and CIS vs. MS, relaxometry within lesions in RR vs. SP, sodium ion concentration in HC vs. CIS, and microstructural alterations were involved across all tasks.

Implications of quantitative susceptibility mapping at 7 Tesla MRI for microbleeds detection in cerebral small vessel disease

Background

Cerebral microbleeds (MBs) are a hallmark of cerebral small vessel disease (CSVD) and can be found on T2*-weighted sequences on MRI. Quantitative susceptibility mapping (QSM) is a postprocessing method that also enables MBs identification and furthermore allows to differentiate them from calcifications.

Aims

We explored the implications of using QSM at submillimeter resolution for MBs detection in CSVD.

Methods

Both 3 and 7 Tesla (T) MRI were performed in elderly participants without MBs and patients with CSVD. MBs were quantified on T2*-weighted imaging and QSM. Differences in the number of MBs were assessed, and subjects were classified in CSVD subgroups or controls both on 3T T2*-weighted imaging and 7T QSM.

Results

48 participants [mean age (SD) 70.9 (8.8) years, 48% females] were included: 31 were healthy controls, 6 probable cerebral amyloid angiopathy (CAA), 9 mixed CSVD, and 2 were hypertensive arteriopathy [HA] patients. After accounting for the higher number of MBs detected at 7T QSM (Median = Mdn; Mdn7T−QSM = 2.5; Mdn3T−T2 = 0; z = 4.90; p &lt; 0.001) and false positive MBs (6.1% calcifications), most healthy controls (80.6%) demonstrated at least one MB and more MBs were discovered in the CSVD group.

Conclusions

Our observations suggest that QSM at submillimeter resolution improves the detection of MBs in the elderly human brain. A higher prevalence of MBs than so far known in healthy elderly was revealed.

Depth-wise profiles of iron and myelin in the cortex and white matter using χ-separation: A preliminary study

The in-vivo profiling of iron and myelin across cortical depths and underlying white matter has important implications for advancing knowledge about their roles in brain development and degeneration. Here, we utilize χ-separation, a recently-proposed advanced susceptibility mapping that creates positive (χ_pos) and negative (χ_neg) susceptibility maps, to generate the depth-wise profiles of χ_pos and χ_neg as surrogate biomarkers for iron and myelin, respectively. Two regional sulcal fundi of precentral and middle frontal areas are profiled and compared with findings from previous studies. The results show that the χ_pos profiles peak at superificial white matter (SWM), which is an area beneath cortical gray matter known to have the highest accumulation of iron within the cortex and white matter. On the other hand, the χ_neg profiles increase in SWM toward deeper white matter. These characteristics in the two profiles are in agreement with histological findings of iron and myelin. Furthermore, the χ_neg profiles report regional differences that agree with well-known distributions of myelin concentration. When the two profiles are compared with those of QSM and R₂*, different shapes and peak locations are observed. This preliminary study offers an insight into one of the possible applications of χ-separation for exploring microstructural information of the human brain, as well as clinical applications in monitoring changes of iron and myelin in related diseases.

Quantitative magnetic resonance imaging biomarkers for cortical pathology in multiple sclerosis at 7 T

Substantial cortical gray matter tissue damage, which correlates with clinical disease severity, has been revealed in multiple sclerosis (MS) using advanced magnetic resonance imaging (MRI) methods at 3 T and the use of ultra-high field, as well as in histopathology studies. While clinical assessment mainly focuses on lesions using T 1 - and T 2 -weighted MRI, quantitative MRI (qMRI) methods are capable of uncovering subtle microstructural changes. The aim of this ultra-high field study is to extract possible future MR biomarkers for the quantitative evaluation of regional cortical pathology. Because of their sensitivity to iron, myelin, and in part specifically to cortical demyelination, T 1 , T 2 , R 2 * , and susceptibility mapping were performed including two novel susceptibility markers; in addition, cortical thickness as well as the volumes of 34 cortical regions were computed. Data were acquired in 20 patients and 16 age- and sex-matched healthy controls. In 18 cortical regions, large to very large effect sizes (Cohen's d ≥ 1) and statistically significant differences in qMRI values between patients and controls were revealed compared with only four regions when using more standard MR measures, namely, volume and cortical thickness. Moreover, a decrease in all susceptibility contrasts ( χ , χ + , χ - ) and R 2 * values indicates that the role of cortical demyelination might outweigh inflammatory processes in the form of iron accumulation in cortical MS pathology, and might also indicate iron loss. A significant association between susceptibility contrasts as well as R 2 * of the caudal middle frontal gyrus and disease duration was found (adjusted R² : 0.602, p = 0.0011). Quantitative MRI parameters might be more sensitive towards regional cortical pathology compared with the use of conventional markers only and therefore may play a role in early detection of tissue damage in MS in the future.

Biophysical mechanisms underlying the effects of static magnetic fields on biological systems

With the widespread use of static magnetic fields (SMFs) in medicine, it is imperative to explore the biological effects of SMFs and the mechanisms underlying their effects on biological systems. The presence of magnetic materials within cells and organisms could affect various biological metabolism and processes, including stress responses, proliferation, and structural alignment. SMFs were generally found to be safe at the organ and organism levels. However. human subjects exposed to strong SMFs have reported side effects. In this review, we combined the magnetic properties of biological samples to illustrate the mechanism of action of SMFs on biological systems from a biophysical point of view. We suggest that the mechanisms of action of SMFs on biological systems mainly include the induction of electric fields and currents, generation of magnetic effects, and influence of electron spins. An electrolyte flowing in a static magnetic field generates an induced current and an electric field. Magnetomechanical effects include orientation effects upon subjecting biological samples to SMFs and movement of biological samples in strong field gradients. SMFs are thought to affect biochemical reaction rates and yields by influencing electron spin. This paper helps people how can harness the favorable biological effects of SMFs.

χ-Separation Imaging for Diagnosis of Multiple Sclerosis versus Neuromyelitis Optica Spectrum Disorder

Background Use of χ-separation imaging can provide surrogates for iron and myelin that relate closely to abnormal changes in multiple sclerosis (MS) lesions. Purpose To evaluate the appearances of MS and neuromyelitis optica spectrum disorder (NMOSD) brain lesions on χ-separation maps and explore their diagnostic value in differentiating the two diseases in comparison with previously reported diagnostic criteria. Materials and Methods This prospective study included individuals with MS or NMOSD who underwent χ-separation imaging from October 2017 to October 2020. Positive (χ_pos) and negative (χ_neg) susceptibility were estimated separately by using local frequency shifts and calculating R2' (R2' = R2* - R2). R2 mapping was performed with a machine learning approach. For each lesion, presence of the central vein sign (CVS) and paramagnetic rim sign (PRS) and signal characteristics on χ_neg and χ_pos maps were assessed and compared. For each participant, the proportion of lesions with CVS, PRS, and hypodiamagnetism was calculated. Diagnostic performances were assessed using receiver operating characteristic (ROC) curve analysis. Results A total of 32 participants with MS (mean age, 34 years ± 10 [SD]; 25 women, seven men) and 15 with NMOSD (mean age, 52 years ± 17; 14 women, one man) were evaluated, with a total of 611 MS and 225 NMOSD brain lesions. On the χ_neg maps, 80.2% (490 of 611) of MS lesions were categorized as hypodiamagnetic versus 13.8% (31 of 225) of NMOSD lesions (P < .001). Lesion appearances on the χ_pos maps showed no evidence of a difference between the two diseases. In per-participant analysis, participants with MS showed a higher proportion of hypodiamagnetic lesions (83%; IQR, 72-93) than those with NMOSD (6%; IQR, 0-14; P < .001). The proportion of hypodiamagnetic lesions achieved excellent diagnostic performance (area under the ROC curve, 0.96; 95% CI: 0.91, 1.00). Conclusion On χ-separation maps, multiple sclerosis (MS) lesions tend to be hypodiamagnetic, which can serve as an important hallmark to differentiate MS from neuromyelitis optica spectrum disorder. © RSNA, 2022 Supplemental material is available for this article.

Susceptibility networks reveal independent patterns of brain iron abnormalities in multiple sclerosis

Brain iron homeostasis is necessary for healthy brain function. MRI and histological studies have shown altered brain iron levels in the brains of patients with multiple sclerosis (MS), particularly in the deep gray matter (DGM). Previous studies were able to only partially separate iron-modifying effects because of incomplete knowledge of iron-modifying processes and influencing factors. It is therefore unclear to what extent and at which stages of the disease different processes contribute to brain iron changes. We postulate that spatially covarying magnetic susceptibility networks determined with Independent Component Analysis (ICA) reflect, and allow for the study of, independent processes regulating iron levels. We applied ICA to quantitative susceptibility maps for 170 individuals aged 9-81 years without neurological disease ("Healthy Aging" (HA) cohort), and for a cohort of 120 patients with MS and 120 age- and sex-matched healthy controls (HC; together the "MS/HC" cohort). Two DGM-associated "susceptibility networks" identified in the HA cohort (the Dorsal Striatum and Globus Pallidus Interna Networks) were highly internally reproducible (i.e. "robust") across multiple ICA repetitions on cohort subsets. DGM areas overlapping both robust networks had higher susceptibility levels than DGM areas overlapping only a single robust network, suggesting that these networks were caused by independent processes of increasing iron concentration. Because MS is thought to accelerate brain aging, we hypothesized that associations between age and the two robust DGM-associated networks would be enhanced in patients with MS. However, only one of these networks was altered in patients with MS, and it had a null age association in patients with MS rather than a stronger association. Further analysis of the MS/HC cohort revealed three additional disease-related networks (the Pulvinar, Mesencephalon, and Caudate Networks) that were differentially altered between patients with MS and HCs and between MS subtypes. Exploratory regression analyses of the disease-related networks revealed differential associations with disease duration and T2 lesion volume. Finally, analysis of ROI-based disease effects in the MS/HC cohort revealed an effect of disease status only in the putamen ROI and exploratory regression analysis did not show associations between the caudate and pulvinar ROIs and disease duration or T2 lesion volume, showing the ICA-based approach was more sensitive to disease effects. These results suggest that the ICA network framework increases sensitivity for studying patterns of brain iron change, opening a new avenue for understanding brain iron physiology under normal and disease conditions.

Quantitative susceptibility mapping (QSM) and R2* of silent cerebral infarcts in sickle cell anemia

Silent cerebral infarction (SCI) is the most commonly reported radiological abnormality in patients with sickle cell anemia (SCA) and is associated with future clinical stroke risk. To date, there have been few histological and quantitative MRI studies of SCI and multiple radiological definitions exist. As a result, the tissue characteristics and composition of SCI remain elusive. The objective of this work was therefore to investigate the composition of segmented SCI lesions using quantitative MRI for R2 * and quantitative magnetic susceptibility mapping (QSM). 211 SCI lesions were segmented from 32 participants with SCA and 6 controls. SCI were segmented according to two definitions (FLAIR+/-T1w-based threshold) using a semi-automated pipeline. Magnetic susceptibility (χ) and R2 * maps were calculated from a multi-echo gradient echo sequence and mean SCI values were compared to an equivalent region of interest in normal appearing white matter (NAWM). SCI χ and R2 * were investigated as a function of SCI definition, patient demographics, anatomical location, and cognition. Compared to NAWM, SCI were significantly less diamagnetic (χ = -0.0067 ppm vs. -0.0153 ppm, p < 0.001) and had significantly lower R2 * (16.7 s-1 vs. 19.2 s-1, p < 0.001). SCI definition had a significant effect on the mean SCI χ and R2 * , with lesions becoming significantly less diamagnetic and having significantly lower R2 * after the application of a more stringent T1w-based threshold. SCI-NAWM R2 * decrease was significantly greater in patients with SCA compared with controls (-2.84 s-1 vs. -0.64 s-1, p < 0.0001). No significant association was observed between mean SCI-NAWM χ or R2* differences and subject age, lesion anatomical location, or cognition. The increased χ and decreased R2 * in SCI relative to NAWM observed in both patients and controls is indicative of lower myelin or increased water content within the segmented lesions. The significant SCI-NAWM R2 * differences observed between SCI in patients with SCA and controls suggests there may be differences in tissue composition relative to NAWM in SCI in the two populations. Quantitative MRI techniques such as QSM and R2 * mapping can be used to enhance our understanding of the pathophysiology and composition of SCI in patients with SCA as well as controls.

Quantitative susceptibility mapping of basal ganglia iron is associated with cognitive and motor functions that distinguish spinocerebellar ataxia type 6 and type 3

Background

In spinocerebellar ataxia type 3 (SCA3), volume loss has been reported in the basal ganglia, an iron-rich brain region, but iron content has not been examined. Recent studies have reported that patients with SCA6 have markedly decreased iron content in the cerebellar dentate, coupled with severe volume loss. Changing brain iron levels can disrupt cognitive and motor functions, yet this has not been examined in the SCAs, a disease in which iron-rich regions are affected.

Methods

In the present study, we used quantitative susceptibility mapping (QSM) to measure tissue magnetic susceptibility (indicating iron concentration), structural volume, and normalized susceptibility mass (indicating iron content) in the cerebellar dentate and basal ganglia in people with SCA3 (n = 10) and SCA6 (n = 6) and healthy controls (n = 9). Data were acquired using a 7T Philips MRI scanner. Supplemental measures assessed motor, cognitive, and mood domains.

Results

Putamen volume was lower in both SCA groups relative to controls, replicating prior findings. Dentate susceptibility mass and volume in SCA6 was lower than in SCA3 or controls, also replicating prior findings. The novel finding was that higher basal ganglia susceptibility mass in SCA6 correlated with lower cognitive performance and greater motor impairment, an association that was not observed in SCA3. Cerebellar dentate susceptibility mass, however, had the opposite relationship with cognition and motor function in SCA6, suggesting that, as dentate iron is depleted, it relocated to the basal ganglia, which contributed to cognitive and motor decline. By contrast, basal ganglia volume loss, rather than iron content, appeared to drive changes in motor function in SCA3.

Conclusion

The associations of higher basal ganglia iron with lower motor and cognitive function in SCA6 but not in SCA3 suggest the potential for using brain iron deposition profiles beyond the cerebellar dentate to assess disease states within the cerebellar ataxias. Moreover, the role of the basal ganglia deserves greater attention as a contributor to pathologic and phenotypic changes associated with SCA.

Quantitative susceptibility mapping versus phase imaging to identify multiple sclerosis iron rim lesions with demyelination

BACKGROUND AND PURPOSE

To compare quantitative susceptibility mapping (QSM) and high-pass-filtered (HPF) phase imaging for (1) identifying chronic active rim lesions with more myelin damage and (2) distinguishing patients with increased clinical disability in multiple sclerosis.

METHODS

Eighty patients were scanned with QSM for paramagnetic rim detection and Fast Acquisition with Spiral Trajectory and T2prep for myelin water fraction (MWF). Chronic lesions were classified based on the presence/absence of rim on HPF and QSM images. A lesion-level linear mixed-effects model with MWF as the outcome was used to compare myelin damage among the lesion groups. A multiple patient-level linear regression model was fit to establish the association between Expanded Disease Status Scale (EDSS) and the log of the number of rim lesions.

RESULTS

Of 2062 lesions, 188 (9.1%) were HPF rim+/QSM rim+, 203 (9.8%) were HPF rim+/QSM rim-, and the remainder had no rim. In the linear mixed-effects model, HPF rim+/QSM rim+ lesions had significantly lower MWF than both HPF rim+/QSM rim- (p < .001) and HPF rim-/QSM rim- (p < .001) lesions, while the MWF difference between HPF rim+/QSM rim- and HPF rim-/QSM rim- lesions was not statistically significant (p = .130). Holding all other factors constant, the log number of QSM rim+ lesion was associated with EDSS increase (p = .044). The association between the log number of HPF rim+ lesions and EDSS was not statistically significant (p = .206).

CONCLUSIONS

QSM identifies paramagnetic rim lesions that on average have more myelin damage and stronger association with clinical disability than those detected by phase imaging.

Linking Plasma Amyloid Beta and Neurofilament Light Chain to Intracortical Myelin Content in Cognitively Normal Older Adults

Evidence suggests that lightly myelinated cortical regions are vulnerable to aging and Alzheimer’s disease (AD). However, it remains unknown whether plasma markers of amyloid and neurodegeneration are related to deficits in intracortical myelin content, and whether this relationship, in turn, is associated with altered patterns of resting-state functional connectivity (rs-FC). To shed light into these questions, plasma levels of amyloid-β fragment 1–42 (Aβ_1–42) and neurofilament light chain (NfL) were measured using ultra-sensitive single-molecule array (Simoa) assays, and the intracortical myelin content was estimated with the ratio T1-weigthed/T2-weighted (T1w/T2w) in 133 cognitively normal older adults. We assessed: (i) whether plasma Aβ_1–42 and/or NfL levels were associated with intracortical myelin content at different cortical depths and (ii) whether cortical regions showing myelin reductions also exhibited altered rs-FC patterns. Surface-based multiple regression analyses revealed that lower plasma Aβ_1–42 and higher plasma NfL were associated with lower myelin content in temporo-parietal-occipital regions and the insular cortex, respectively. Whereas the association with Aβ_1–42 decreased with depth, the NfL-myelin relationship was most evident in the innermost layer. Older individuals with higher plasma NfL levels also exhibited altered rs-FC between the insula and medial orbitofrontal cortex. Together, these findings establish a link between plasma markers of amyloid/neurodegeneration and intracortical myelin content in cognitively normal older adults, and support the role of plasma NfL in boosting aberrant FC patterns of the insular cortex, a central brain hub highly vulnerable to aging and neurodegeneration.

Magnetic Susceptibility Source Separation Solely from Gradient Echo Data: Histological Validation

Quantitative susceptibility mapping (QSM) facilitates mapping of the bulk magnetic susceptibility of tissue from the phase of complex gradient echo (GRE) MRI data. QSM phase processing combined with an R2* model of magnitude of multiecho gradient echo data (R2*QSM) allows separation of dia- and para-magnetic components (e.g., myelin and iron) that contribute constructively to R2* value but destructively to the QSM value of a voxel. This R2*QSM technique is validated against quantitative histology—optical density of myelin basic protein and Perls’ iron histological stains of rim and core of 10 ex vivo multiple sclerosis lesions, as well as neighboring normal appearing white matter. We found that R2*QSM source maps are in good qualitative agreement with histology, e.g., showing increased iron concentration at the edge of the rim+ lesions and myelin loss in the lesions’ core. Furthermore, our results indicate statistically significant correlation between paramagnetic and diamagnetic tissue components estimated with R2*QSM and optical densities of Perls’ and MPB stains. These findings provide direct support for the use of R2*QSM magnetic source separation based solely on GRE complex data to characterize MS lesion composition.

Susceptibility source separation from gradient echo data using magnitude decay modeling

BACKGROUND AND PURPOSE

The objective is to demonstrate feasibility of separating magnetic sources in quantitative susceptibility mapping (QSM) by incorporating magnitude decay rates R 2 ∗ $R_2^{\rm{*}}$ in gradient echo (GRE) MRI.

METHODS

Magnetic susceptibility source separation was developed using R 2 ∗ $R_2^{\rm{*}}$ and compared with a prior method using R 2 ' = R 2 ∗ - R 2 ${R^{\prime}_2} = R_2^* - {R_2}$ that required an additional sequence to measure the transverse relaxation rate R₂ . Both susceptibility separation methods were compared in multiple sclerosis (MS) patients (n = 17). Susceptibility values of negative sources estimated with R 2 ∗ $R_2^{\rm{*}}$ -based source separation in a set of enhancing MS lesions (n = 44) were correlated against longitudinal myelin water fraction (MWF) changes.

RESULTS

In in vivo data, linear regression of the estimated χ + ${\chi}^{+}$ and χ - ${\chi}^{-}$ susceptibility values between the R 2 ∗ $R_2^*$ - and the R 2 ' ${R^{\prime}_2}$ -based separation methods performed across 182 segmented lesions revealed correlation coefficient r = .96 and slope close .99. Correlation analysis in enhancing lesions revealed a significant positive association between the χ - ${\chi}^{-}$ increase at 1-year post-onset relative to 0 year and the MWF increase at 1 year relative to 0 year (β = -0.144, 95% confidence interval: [-0.199, -0.1], p = .0008) and good agreement between R 2 ' ${R^{\prime}_2}$ and R 2 ∗ $R_2^*$ methods (r = .79, slope = .95).

CONCLUSIONS

Separation of magnetic sources based solely on GRE complex data is feasible by combining magnitude decay rate modeling and phase-based QSM and χ - ${\chi}^{-}$ change may serve as a biomarker for myelin recovery or damage in acute MS lesions.

Phenotypic and genetic associations of quantitative magnetic susceptibility in UK Biobank brain imaging

A key aim in epidemiological neuroscience is identification of markers to assess brain health and monitor therapeutic interventions. Quantitative susceptibility mapping (QSM) is an emerging magnetic resonance imaging technique that measures tissue magnetic susceptibility and has been shown to detect pathological changes in tissue iron, myelin and calcification. We present an open resource of QSM-based imaging measures of multiple brain structures in 35,273 individuals from the UK Biobank prospective epidemiological study. We identify statistically significant associations of 251 phenotypes with magnetic susceptibility that include body iron, disease, diet and alcohol consumption. Genome-wide associations relate magnetic susceptibility to 76 replicating clusters of genetic variants with biological functions involving iron, calcium, myelin and extracellular matrix. These patterns of associations include relationships that are unique to QSM, in particular being complementary to T2* signal decay time measures. These new imaging phenotypes are being integrated into the core UK Biobank measures provided to researchers worldwide, creating the potential to discover new, non-invasive markers of brain health.

Magnetic susceptibility anisotropy in normal appearing white matter in multiple sclerosis from single-orientation acquisition

•

Orientation dependence of QSM is studied in a large cohort of MS patients.

•

Apparent magnetic susceptibility anisotropy (MSA) obtained from single-orientation QSM.

•

Apparent MSA found decreased in optic radiation (OR) of MS patients.

•

Apparent MSA decreases with lesion load in OR and with disease duration in splenium.

•

Negative apparent MSA observed in SLF indicates limitations of the proposed method.

Quantitative susceptibility mapping (QSM) has been successfully applied to study changes in deep grey matter nuclei as well as in lesional tissue, but its application to white matter has been complicated by the observed orientation dependence of gradient echo signal. The anisotropic susceptibility tensor is thought to be at the origin of this orientation dependence, and magnetic susceptibility anisotropy (MSA) derived from this tensor has been proposed as a marker of the state and integrity of the myelin sheath and may therefore be of particular interest for the study of demyelinating pathologies such as multiple sclerosis (MS). Reconstruction of the susceptibility tensor, however, requires repeated measurements with multiple head orientations, rendering the approach impractical for clinical applications.

In this study, we combined single-orientation QSM with fibre orientation information to assess apparent MSA in three white matter tracts, i.e., optic radiation (OR), splenium of the corpus callosum (SCC), and superior longitudinal fascicle (SLF), in two cohorts of 64 healthy controls and 89 MS patients. The apparent MSA showed a significant decrease in optic radiation in the MS cohort compared with healthy controls. It decreased in the MS cohort with increasing lesion load in OR and with disease duration in the splenium. All of this suggests demyelination in normal appearing white matter. However, the apparent MSA observed in the SLF pointed to potential systematic issues that require further exploration to realize the full potential of the presented approach. Despite the limitations of such single-orientation ROI-specific estimation, we believe that our clinically feasible approach to study degenerative changes in WM is worthy of further investigation.

[Brain Iron Imaging in Aging and Cognitive Disorders: MRI Approaches]

Iron has a vital role in the human body, including the central nervous system. Increased deposition of iron in the brain has been reported in aging and important neurodegenerative diseases. Owing to the unique magnetic resonance properties of iron, MRI has great potential for in vivo assessment of iron deposition, distribution, and non-invasive quantification. In this paper, we will review the MRI methods for iron assessment and their changes in aging and neurodegenerative diseases, focusing on Alzheimer's disease. In addition, we will summarize the limitations of current approaches and introduce new areas and MRI methods for iron imaging that are expected in the future.

Blip up-down acquisition for spin- and gradient-echo imaging (BUDA-SAGE) with self-supervised denoising enables efficient T2 , T2 *, para- and dia-magnetic susceptibility mapping

PURPOSE

To rapidly obtain high resolution T₂ , T₂ *, and quantitative susceptibility mapping (QSM) source separation maps with whole-brain coverage and high geometric fidelity.

METHODS

We propose Blip Up-Down Acquisition for Spin And Gradient Echo imaging (BUDA-SAGE), an efficient EPI sequence for quantitative mapping. The acquisition includes multiple T₂ *-, T₂ '-, and T₂ -weighted contrasts. We alternate the phase-encoding polarities across the interleaved shots in this multi-shot navigator-free acquisition. A field map estimated from interim reconstructions was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to eliminate distortion. A self-supervised neural network (NN), MR-Self2Self (MR-S2S), was used to perform denoising to boost SNR. Using Slider encoding allowed us to reach 1 mm isotropic resolution by performing super-resolution reconstruction on volumes acquired with 2 mm slice thickness. Quantitative T₂ (=1/R₂ ) and T₂ * (=1/R₂ *) maps were obtained using Bloch dictionary matching on the reconstructed echoes. QSM was estimated using nonlinear dipole inversion on the gradient echoes. Starting from the estimated R₂ /R₂ * maps, R₂ ' information was derived and used in source separation QSM reconstruction, which provided additional para- and dia-magnetic susceptibility maps.

RESULTS

In vivo results demonstrate the ability of BUDA-SAGE to provide whole-brain, distortion-free, high-resolution, multi-contrast images and quantitative T₂ /T₂ * maps, as well as yielding para- and dia-magnetic susceptibility maps. Estimated quantitative maps showed comparable values to conventional mapping methods in phantom and in vivo measurements.

CONCLUSION

BUDA-SAGE acquisition with self-supervised denoising and Slider encoding enables rapid, distortion-free, whole-brain T₂ /T₂ * mapping at 1 mm isotropic resolution under 90 s.

Subsecond accurate myelin water fraction reconstruction from FAST-T2 data with 3D UNET

PURPOSE

To develop a 3D UNET convolutional neural network for rapid extraction of myelin water fraction (MWF) maps from six-echo fast acquisition with spiral trajectory and T₂ -prep data and to evaluate its accuracy in comparison with multilayer perceptron (MLP) network.

METHODS

The MWF maps were extracted from 138 patients with multiple sclerosis using an iterative three-pool nonlinear least-squares algorithm (NLLS) without and with spatial regularization (srNLLS), which were used as ground-truth labels to train, validate, and test UNET and MLP networks as a means to accelerate data fitting. Network testing was performed in 63 patients with multiple sclerosis and a numerically simulated brain phantom at SNR of 200, 100 and 50.

RESULTS

Simulations showed that UNET reduced the MWF mean absolute error by 30.1% to 56.4% and 16.8% to 53.6% over the whole brain and by 41.2% to 54.4% and 21.4% to 49.4% over the lesions for predicting srNLLS and NLLS MWF, respectively, compared to MLP, with better performance at lower SNRs. UNET also outperformed MLP for predicting srNLLS MWF in the in vivo multiple-sclerosis brain data, reducing mean absolute error over the whole brain by 61.9% and over the lesions by 67.5%. However, MLP yielded 41.1% and 51.7% lower mean absolute error for predicting in vivo NLLS MWF over the whole brain and the lesions, respectively, compared with UNET. The whole-brain MWF processing time using a GPU was 0.64 seconds for UNET and 0.74 seconds for MLP.

CONCLUSION

Subsecond whole-brain MWF extraction from fast acquisition with spiral trajectory and T₂ -prep data using UNET is feasible and provides better accuracy than MLP for predicting MWF output of srNLLS algorithm.