On the separation of susceptibility sources in quantitative susceptibility mapping: Theory and phantom validation with an in vivo application to multiple sclerosis lesions of different age

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.

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.

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.

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.

Optimal echo times for quantitative susceptibility mapping: A test-retest study on basal ganglia and subcortical brain nuclei

Quantitative Susceptibility Mapping (QSM) is a recent MRI-technique able to quantify the bulk magnetic susceptibility of myelin, iron, and calcium in the brain. Its variability across different acquisition parameters has prompted the need for standardisation across multiple centres and MRI vendors. However, a high level of agreement between repeated imaging acquisitions is equally important. With this study we aimed to assess the inter-scan repeatability of an optimised multi-echo GRE sequence in 28 healthy volunteers. We extracted and compared the susceptibility measures from the scan and rescan acquisitions across 7 bilateral brain regions (i.e., 14 regions of interest (ROIs)) relevant for neurodegeneration. Repeatability was first assessed while reconstructing QSM with a fixed number of echo times (i.e., 8). Excellent inter-scan repeatability was found for putamen, globus pallidus and caudate nucleus, while good performance characterised the remaining structures. An increased variability was instead noted for small ROIs like red nucleus and substantia nigra. Secondly, we assessed the impact exerted on repeatability by the number of echoes used to derive QSM maps. Results were impacted by this parameter, especially in smaller regions. Larger brain structures, on the other hand, showed more consistent performance. Nevertheless, with either 8 or 7 echoes we managed to obtain good inter-scan repeatability on almost all ROIs. These findings indicate that the designed acquisition/reconstruction protocol has wide applicability, particularly in clinical or research settings involving longitudinal acquisitions (e.g. rehabilitation studies).

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.

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.

From pathology to MRI and back: Clinically relevant biomarkers of multiple sclerosis lesions

Focal lesions in both white and gray matter are characteristic of multiple sclerosis (MS). Histopathological studies have helped define the main underlying pathological processes involved in lesion formation and evolution, serving as a gold standard for many years. However, histopathology suffers from an intrinsic bias resulting from over-reliance on tissue samples from late stages of the disease or atypical cases and is inadequate for routine patient assessment. Pathological-radiological correlative studies have established advanced MRI's sensitivity to several relevant MS-pathological substrates and its practicality for assessing dynamic changes and following lesions over time. This review focuses on novel imaging techniques that serve as biomarkers of critical pathological substrates of MS lesions: the central vein, chronic inflammation, remyelination and repair, and cortical lesions. For each pathological process, we address the correlative value of MRI to MS pathology, its contribution in elucidating MS pathology in vivo, and the clinical utility of the imaging biomarker.

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.