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

Depthwise cortical iron relates to functional connectivity and fluid cognition in healthy aging

Age-related differences in fluid cognition have been associated with both the merging of functional brain networks, defined from resting-state functional magnetic resonance imaging (rsfMRI), and with elevated cortical iron, assessed by quantitative susceptibility mapping (QSM). Limited information is available, however, regarding the depthwise profile of cortical iron and its potential relation to functional connectivity. Here, using an adult lifespan sample (n = 138; 18-80 years), we assessed relations among graph theoretical measures of functional connectivity, column-based depthwise measures of cortical iron, and fluid cognition (i.e., tests of memory, perceptual-motor speed, executive function). Increased age was related both to less segregated functional networks and to increased cortical iron, especially for superficial depths. Functional network segregation mediated age-related differences in memory, whereas depthwise iron mediated age-related differences in general fluid cognition. Lastly, higher mean parietal iron predicted lower network segregation for adults younger than 45 years of age. These findings suggest that functional connectivity and depthwise cortical iron have distinct, complementary roles in the relation between age and fluid cognition in healthy adults.

In-vivo high-resolution χ-separation at 7T

A recently introduced quantitative susceptibility mapping (QSM) technique, χ-separation, offers the capability to separate paramagnetic (χpara) and diamagnetic (χdia) susceptibility distribution within the brain. In-vivo high-resolution mapping of iron and myelin distribution, estimated by χ-separation, could provide a deeper understanding of brain substructures, assisting the investigation of their functions and alterations. This can be achieved using 7T MRI, which benefits from a high signal-to-noise ratio and susceptibility effects. However, applying χ-separation at 7T presents difficulties due to the requirement of an R2 map, coupled with issues such as high specific absorption rate (SAR), large B1 transmit field inhomogeneities, and prolonged scan time. To address these challenges, we developed a novel deep neural network, R2PRIMEnet7T, designed to convert a 7T R2* map into a 3T R2' map. Building on this development, we present a new pipeline for χ-separation at 7T, enabling us to generate high-resolution χ-separation maps from multi-echo gradient-echo data. The proposed method is compared with alternative pipelines, such as an end-to-end network and linearly-scaled R2', and is validated against χ-separation maps at 3T, demonstrating its accuracy. The 7T χ-separation maps generated by the proposed method exhibit similar contrasts to those from 3T, while 7T high-resolution maps offer enhanced clarity and detail. Quantitative analysis confirms that the proposed method surpasses the alternative pipelines. The proposed method results well delineate the detailed brain structures associated with iron and myelin. This new pipeline holds promise for analyzing iron and myelin concentration changes in various neurodegenerative diseases through precise structural examination.

χ-sepnet: Deep Neural Network for Magnetic Susceptibility Source Separation

Magnetic susceptibility source separation (χ-separation), an advanced quantitative susceptibility mapping (QSM) method, enables the separate estimation of paramagnetic and diamagnetic susceptibility source distributions in the brain. Similar to QSM, it requires solving the ill-conditioned problem of dipole inversion, suffering from so-called streaking artifacts. Additionally, the method utilizes reversible transverse relaxation (R 2 ' = R 2 * - R 2 $$ {R}_2^{\prime }={R}_2^{\ast }-{R}_2 $$ ) to complement frequency shift information for estimating susceptibility source concentrations, requiring time-consuming data acquisition forR 2 $$ {R}_2 $$ (e.g., multi-echo spin-echo) in addition to multi-echo GRE data forR 2 * $$ {R}_2^{\ast } $$ . To address these challenges, we develop a new deep learning network, χ-sepnet, and propose two deep learning-based susceptibility source separation pipelines, χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ for inputs with multi-echo GRE and multi-echo spin-echo (or turbo spin-echo) and χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ for input with multi-echo GRE only. The neural network is trained using multiple head orientation data that provide streaking artifact-free labels, generating high-quality χ-separation maps. The evaluation of the pipelines encompasses both qualitative and quantitative assessments in healthy subjects, and visual inspection of lesion characteristics in multiple sclerosis patients. The susceptibility source-separated maps of the proposed pipelines delineate detailed brain structures with substantially reduced artifacts compared to those from the conventional regularization-based reconstruction methods. In quantitative analysis, χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ achieves the best outcomes followed by χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ , outperforming the conventional methods. When the lesions of multiple sclerosis patients are classified into subtypes, most lesions are identified as the same subtype in the maps from χ-sepnet-R 2 ' $$ {R}_2^{\prime } $$ and χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ (paramagnetic susceptibility: 99.6% and diamagnetic susceptibility: 98.4%; both out of 250 lesions). The χ-sepnet-R 2 * $$ {R}_2^{\ast } $$ pipeline, which only requires multi-echo GRE data, has demonstrated its potential to offer broad clinical and scientific applications, although further evaluations for various diseases and pathological conditions are necessary.

Atrophy-Independent and Dependent Iron and Myelin Changes in Deep Gray Matter of Multiple Sclerosis: A Longitudinal Study Using χ-Separation Imaging

RATIONALE AND OBJECTIVES

To investigate iron and myelin changes in deep gray matter (DGM) of relapsing-remitting multiple sclerosis (RRMS) patients and their relationship to atrophy by χ-separation imaging.

MATERIALS AND METHODS

33 RRMS patients and 34 healthy controls (HC) were included in this study. The χ-separation map reconstructed from a 3D multi-echo gradient echo scan was used to measure the positive susceptibility (χpos) and negative susceptibility (χneg) of DGM. To take into account the effect of atrophy, susceptibility mass of DGM was calculated by multiplying volume by the mean bulk susceptibility. Differences in MRI metrics between baseline patients, follow-up patients, and HC were compared respectively.

RESULTS

Compared to HC, χpos of basal ganglia were significantly increased in follow-up patients (P < 0.05). The χpos of pallidum was significantly higher in follow-up patients than that in baseline patients (P = 0.006). The χneg of caudate, pallidum and hippocampus in baseline and follow-up patients was significantly higher than that in HC (P < 0.05). When taking into account the effect of atrophy, there was a significant decrease in χpos mass and a significant increase in χneg mass of thalamus, accumbens and amygdala in follow-up patients compared to HC (P < 0.05). The χpos mass of the thalamus was further decreased in follow-up patients compared to baseline patients (P = 0.006).

CONCLUSION

χ-separation imaging could generate independent information on iron and myelin changes in RRMS patients, showing atrophy-dependent iron increase in basal ganglia and atrophy-independent iron and myelin decrease in thalamus.

A human brain atlas of χ-separation for normative iron and myelin distributions

Iron and myelin are primary susceptibility sources in the human brain. These substances are essential for a healthy brain, and their abnormalities are often related to various neurological disorders. Recently, an advanced susceptibility mapping technique, which is referred to as χ-separation (pronounced as "chi"-separation), has been proposed, successfully disentangling paramagnetic iron from diamagnetic myelin. This method provided a new opportunity for generating high-resolution iron and myelin maps of the brain. Utilizing this technique, this study constructs a normative χ-separation atlas from 106 healthy human brains. The resulting atlas provides detailed anatomical structures associated with the distributions of iron and myelin, clearly delineating subcortical nuclei, thalamic nuclei, and white matter fiber bundles. Additionally, susceptibility values in a number of regions of interest are reported along with age-dependent changes. This atlas may have direct applications such as localization of subcortical structures for deep brain stimulation or high-intensity focused ultrasound and also serve as a valuable resource for future research.

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.

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.

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.

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.

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.

Microstructure-Informed Myelin Mapping (MIMM) from Gradient Echo MRI using Stochastic Matching Pursuit

Quantification of the myelin content of the white matter is important for studying demyelination in neurodegenerative diseases such as Multiple Sclerosis (MS), particularly for longitudinal monitoring. A novel noninvasive MRI method, called Microstructure-Informed Myelin Mapping (MIMM), is developed to quantify the myelin volume fraction (MVF) by utilizing a multi gradient echo sequence (mGRE) and a detailed biophysical model of tissue microstructure. Myelin is modeled as anisotropic negative susceptibility source based on the Hollow Cylindrical Fiber Model (HCFM), and iron as isotropic positive susceptibility source in the extracellular region. Voxels with a range of biophysical parameters are simulated to create a dictionary of MR echo time magnitude signals and total susceptibility values. MRI signals measured using a mGRE sequence are then matched voxel-by-voxel to the created dictionary to obtain the spatial distributions of myelin and iron. Three different MIMM versions are presented to deal with the fiber orientation dependent susceptibility effects of the myelin sheaths: a basic variation, which assumes fiber orientation is an unknown to fit, two orientation informed variations, which assume the fiber orientation distribution is available either from a separate diffusion tensor imaging (DTI) acquisition or from a DTI atlas based fiber orientation map. While all showed a significant linear correlation with the reference method based on T2-relaxometry (p < 0.0001), DTI orientation informed and atlas orientation informed variations reduced overestimation at white matter tracts compared to the basic variation. Finally, the implications and usefulness of attaining an additional iron susceptibility distribution map are discussed.

Highlights

novel stochastic matching pursuit algorithm called microstructure-informed myelin mapping (MIMM) is developed to quantify Myelin Volume Fraction (MVF) using Magnetic Resonance Imaging (MRI) and microstructural modeling.utilizes a detailed biophysical model to capture the susceptibility effects on both magnitude and phase to quantify myelin and iron.matter fiber orientation effects are considered for the improved MVF quantification in the major fiber tracts.acquired myelin and iron maps may be utilized to monitor longitudinal disease progress.