Quantitative susceptibility mapping for susceptibility source separation with adaptive relaxometric constant estimation (QSM-ARCS) from solely gradient-echo data.

To separate the contributions of paramagnetic and diamagnetic sources within a voxel, a magnetic susceptibility source separation method based solely on gradient-echo data has been developed. To measure the opposing susceptibility sources more accurately, we propose a novel single-orientation quantitative susceptibility mapping method with adaptive relaxometric constant estimation (QSM-ARCS) for susceptibility source separation. Moreover, opposing susceptibilities and their anisotropic effects were determined in healthy volunteers in the white matter. Multiple spoiled gradient echo and diffusion tensor imaging of ten healthy volunteers was obtained using a 3 T magnetic resonance scanner. After the opposing susceptibility and fractional anisotropy (FA) maps had been reconstructed, the parametric maps were spatially normalized. To evaluate the agreements of QSM-ARCS against the susceptibility source separation method using R2 and R2* maps (χ-separation) by Bland-Altman plots, the opposing susceptibility values were measured using white and deep gray matter atlases. We then evaluated the relationships between the opposing susceptibilities and FAs in the white matter and used a field-to-fiber angle to assess the fiber orientation dependencies of the opposing susceptibilities. The susceptibility maps in QSM-ARCS were successfully reconstructed without large artifacts. In the Bland-Altman analyses, the opposing QSM-ARCS susceptibility values excellently agreed with the χ-separation maps. Significant inverse and proportional correlations were observed between FA and the negative and positive susceptibilities estimated by QSM-ARCS. The fiber orientation dependencies of the negative susceptibility represented a nonmonotonic feature. Conversely, the positive susceptibility increased linearly with the fiber angle with respect to the B0 field. The QSM-ARCS could accurately estimate the opposing susceptibilities, which were identical values of χ-separation, even using gradient echo alone. The opposing susceptibilities might offer direct biomarkers for assessment of the myelin and iron content in glial cells and, through the underlying magnetic sources, provide biologic insights toward clinical transition.

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.

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.

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 R2* with R2' in QSM. The magnitude decay kernel, describing the proportionality coefficient between R2' 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 R2* 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 R2' 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.

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 R2*, 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.

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 approximated R 2 ' map from acquired R 2 * map. The linear relationship between R 2 * and R 2 ' maps was established from the regression coefficients on a small cohort data acquired with both 3D gradient recalled echo data and R 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.

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

Obtaining a histological fingerprint from the in-vivo brain has been a long-standing target of magnetic resonance imaging (MRI). In particular, non-invasive imaging of iron and myelin, which are involved in normal brain functions and are histopathological hallmarks in neurodegenerative diseases, has practical utilities in neuroscience and medicine. Here, we propose a biophysical model that describes the individual contribution of paramagnetic (e.g., iron) and diamagnetic (e.g., myelin) susceptibility sources to the frequency shift and transverse relaxation of MRI signals. Using this model, we develop a method, χ-separation, that generates the voxel-wise distributions of the two sources. The method is validated using computer simulation and phantom experiments, and applied to ex-vivo and in-vivo brains. The results delineate the well-known histological features of iron and myelin in the specimen, healthy volunteers, and multiple sclerosis patients. This new technology may serve as a practical tool for exploring the microstructural information of the brain.

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.

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.

Distinguishing microgliosis and tau deposition in the mouse brain using paramagnetic and diamagnetic susceptibility source separation.

Tauopathies, including Alzheimer's disease (AD), are neurodegenerative disorders characterized by hyperphosphorylated tau protein aggregates in the brain. In addition to protein aggregates, microglia-mediated inflammation and iron dyshomeostasis are other pathological features observed in AD and other tauopathies. It is known that these alterations at the subcellular level occur much before the onset of macroscopic tissue atrophy or cognitive deficits. The ability to detect these microstructural changes with MRI therefore has substantive importance for improved characterization of disease pathogenesis. In this study, we demonstrate that quantitative susceptibility mapping (QSM) with paramagnetic and diamagnetic susceptibility source separation has the potential to distinguish neuropathological alterations in a transgenic mouse model of tauopathy. 3D multi-echo gradient echo data were acquired from fixed brains of PS19 (Tau) transgenic mice and age-matched wild-type (WT) mice (n = 5 each) at 11.7 T. The multi-echo data were fit to a 3-pool complex signal model to derive maps of paramagnetic component susceptibility (PCS) and diamagnetic component susceptibility (DCS). Group-averaged signal fraction and composite susceptibility maps showed significant region-specific differences between the WT and Tau mouse brains. Significant bilateral increases in PCS and |DCS| were observed in specific hippocampal and cortical sub-regions of the Tau mice relative to WT controls. Comparison with immunohistological staining for microglia (Iba1) and phosphorylated-tau (AT8) further indicated that the PCS and DCS differences corresponded to regional microgliosis and tau deposition in the PS19 mouse brains, respectively. The results demonstrate that quantitative susceptibility source separation may provide sensitive imaging markers to detect distinct pathological alterations in tauopathies.

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 R2 . 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.

A novel phantom with dia- and paramagnetic substructure for quantitative susceptibility mapping and relaxometry.

PURPOSE: A phantom is presented in this study that allows for an experimental evaluation of QSM reconstruction algorithms. The phantom contains susceptibility producing particles with dia- and paramagnetic properties embedded in an MRI visible medium and is suitable to assess the performance of algorithms that attempt to separate isotropic dia- and paramagnetic susceptibility at the sub-voxel level. METHODS: The phantom was built from calcium carbonate (diamagnetic) and tungsten carbide particles (paramagnetic) embedded in gelatin and surrounded by agarose gel. Different mass fractions and mixing ratios of both susceptibility sources were used. Gradient echo data were acquired at 1.5 T, 3 T and 7 T. Susceptibility maps were calculated using the MEDI toolbox and relaxation rates ΔR2∗ were determined using exponential fitting. RESULTS: Relaxation rates as well as susceptibility values generally coincide with the theoretical values for particles fulfilling the assumptions of the the static dephasing regime with stronger deviations for relaxation rates at higher field strength and for high susceptibility values. MRI raw data are available for free academic use as supplementary material. CONCLUSIONS: In this study, a susceptibility phantom is presented that might find its application in the development and quantitative validation of current and future QSM reconstruction algorithms which aim to separate the influence of isotropic dia- and paramagnetic substructure in quantitative susceptibility mapping.