Quantitative Susceptibility Mapping (QSM) Algorithms: Mathematical Rationale and Computational Implementations

5D image reconstruction exploiting space-motion-echo sparsity for accelerated free-breathing quantitative liver MRI

Recent advances in 3D non-Cartesian multi-echo gradient-echo (mGRE) imaging and compressed sensing (CS)-based 4D (3D image space + 1D respiratory motion) motion-resolved image reconstruction, which applies temporal total variation to the respiratory motion dimension, have enabled free-breathing liver tissue MR parameter mapping. This technology now allows for robust reconstruction of high-resolution proton density fat fraction (PDFF), R2∗, and quantitative susceptibility mapping (QSM), previously unattainable with conventional Cartesian mGRE imaging. However, long scan times remain a persistent challenge in free-breathing 3D non-Cartesian mGRE imaging. Recognizing that the underlying dimension of the imaging data is essentially 5D (4D + 1D echo signal evolution), we propose a CS-based 5D motion-resolved mGRE image reconstruction method to further accelerate the acquisition. Our approach integrates discrete wavelet transforms along the echo and spatial dimensions into a CS-based reconstruction model and devises a solution algorithm capable of handling such a 5D complex-valued array. Through phantom and in vivo human subject studies, we evaluated the effectiveness of leveraging unexplored correlations by comparing the proposed 5D reconstruction with the 4D reconstruction (i.e., motion-resolved reconstruction with temporal total variation) across a wide range of acceleration factors. The 5D reconstruction produced more reliable and consistent measurements of PDFF, R2∗, and QSM compared to the 4D reconstruction. In conclusion, the proposed 5D motion-resolved image reconstruction demonstrates the feasibility of achieving accelerated, reliable, and free-breathing liver mGRE imaging for the measurement of PDFF, R2∗, and QSM.

Spiral cardiac quantitative susceptibility mapping for differential cardiac chamber oxygenation-Initial validation in relation to invasive blood sampling

PURPOSE

To develop a breath-hold cardiac quantitative susceptibility mapping (QSM) sequence for noninvasive measurement of differential cardiac chamber blood oxygen saturation (ΔSO2).

METHODS

A non-gated three-dimensional stack-of-spirals QSM sequence was implemented to continuously sample the data throughout the cardiac cycle. Measurements of ΔSO2 between the right and left heart chamber obtained by the proposed sequence and a previously validated navigator Cartesian QSM sequence were compared in three cohorts consisting of healthy volunteers, coronavirus disease 2019 survivors, and patients with pulmonary hypertension. In the pulmonary-hypertension cohort, Bland-Altman plots were used to assess the agreement of ΔSO2 values obtained by QSM and those obtained by invasive right heart catheterization (RHC).

RESULTS

Compared with navigator QSM (average acquisition time 419 ± 158 s), spiral QSM reduced the scan time on average by over 20-fold to a 20-s breath-hold. In all three cohorts, spiral QSM and navigator QSM yielded similar ΔSO2. Among healthy volunteers and coronavirus disease 2019 survivors, ΔSO2 was 17.41 ± 4.35% versus 17.67 ± 4.09% for spiral and navigator QSM, respectively. In pulmonary-hypertension patients, spiral QSM showed a slightly smaller ΔSO2 bias and narrower 95% limits of agreement than that obtained by navigator QSM (1.09% ± 6.47% vs. 2.79% ± 6.99%) when compared with right heart catheterization.

CONCLUSION

Breath-hold three-dimensional spiral cardiac QSM for measuring differential cardiac chamber blood oxygenation is feasible and provides values in good agreement with navigator cardiac QSM and with reference right heart catheterization.

Microstructure-Informed Myelin Mapping (MIMM) from routine multi-echo gradient echo data using multiscale physics modeling of iron and myelin effects and QSM

PURPOSE

Myelin quantification is used in the study of demyelination in neurodegenerative diseases. A novel noninvasive MRI method, Microstructure-Informed Myelin Mapping (MIMM), is proposed to quantify the myelin volume fraction (MVF) from a routine multi-gradient echo sequence (mGRE) using a multiscale biophysical signal model of the effects of microstructural myelin and iron.

THEORY AND METHODS

In MIMM, the effects of myelin are modeled based on the Hollow Cylinder Fiber Model accounting for anisotropy, while iron is considered as an isotropic paramagnetic point source. This model is used to create a dictionary of mGRE magnitude signal evolution and total voxel susceptibility using finite elements of size 0.2 μm. Next, voxel-by-voxel stochastic matching pursuit between acquired mGRE data (magnitude+QSM) and the pre-computed dictionary generates quantitative MVF and iron susceptibility maps. Dictionary matching was evaluated under three conditions: (1) without fiber orientation (basic), (2) with fiber orientation obtained using DTI, and (3) with fiber orientation obtained using an atlas (atlas). MIMM was compared with the three-pool complex fitting (3PCF) using T2-relaxometry myelin water fraction (MWF) map as reference.

RESULTS

The DTI MIMM and atlas MIMM approaches were equally effective in reducing the overestimation of MVF in certain white matter tracts observed in the basic MIMM approach, and they both showed good agreement with T2-relaxometry MWF. MIMM MVF reduced myelin overestimation of globus pallidus observed in 3PCF MWF.

CONCLUSION

MIMM processing of mGRE data can provide MVF maps from routine clinical scans without requiring special sequences.

Association of Increased Brain Iron Levels With Anxiety and Motor Dysfunction in Cerebral Small Vessel Disease

AIMS

This study explored the relationships between brain iron levels, emotion, and cognitive and motor function in cerebral small vessel disease (CSVD) patients using quantitative susceptibility mapping (QSM).

METHODS

A total of 208 subjects were enrolled in this study. A brain QSM map was calculated from multiecho GRE data via morphology-enabled dipole inversion with an automatic uniform cerebrospinal fluid zero reference algorithm (MEDI+0). Multiple linear regression analysis was applied to explore the clinical factors influencing cerebral susceptibility in CSVD patients. Correlation analysis and pathway-specific mediation effects between brain iron levels and motor function were investigated.

RESULTS

There were significant differences in the MoCA scores, depression scores, five-repetition sit-to-stand test (5R-STS) time, and susceptibility values of the caudate nucleus and putamen among the three groups (p < 0.05, FDR correction). Age and history of diabetes played crucial roles in brain iron levels in the caudate nucleus and putamen, which may increase iron levels in the basal ganglia, associated with cognitive decline. Notably, the susceptibility values of the left caudate nucleus and putamen were positively correlated with the 5R-STS time in CSVD subjects, and there were significant mediating effects of anxiety on the prediction of motor dysfunction with respect to iron levels in the left putamen in CSVD patients.

CONCLUSION

Age, diabetes status, and anxiety may serve as effective intervention targets for individuals with CSVD, especially individuals with cognitive and motor dysfunction. A greater brain iron burden may be a quantitative imaging marker of cognitive and motor dysfunction in CSVD patients.

TRIAL REGISTRATION

ISRCTN20008650.

Brain iron deposition and cognitive decline in patients with cerebral small vessel disease : a quantitative susceptibility mapping study

BACKGROUND

Quantitative susceptibility mapping (QSM) can study the susceptibility values of brain tissue which allows for noninvasive examination of local brain iron levels in both normal and pathological conditions.

PURPOSE

Our study compares brain iron deposition in gray matter (GM) nuclei between cerebral small vessel disease (CSVD) patients and healthy controls (HCs), exploring factors that affect iron deposition and cognitive function.

MATERIALS AND METHODS

A total of 321 subjects were enrolled in this study. All subjects had cognitive examination including the Stroop color word test (SCWT) and MRI including multiecho gradient echo (mGRE) sequence. The patients with CSVD were divided into mild to moderate group (CSVD-M, total CSVD score ≤ 1) and severe group (CSVD-S, total CSVD score > 1). Morphology-enabled dipole inversion with an automated uniform cerebrospinal fluid zero reference algorithm (MEDI + 0) was used to generate brain QSM maps from mGRE data. Deep gray regional susceptibility values and cognitive function were compared among three groups (CSVD-S, CSVD-M, and HC) using multiple linear regression analysis and mediation effect analysis.

RESULTS

There were significant differences in the SCWT scores and mean susceptibility values of the globus pallidus (GP), putamen (Put), and caudate nucleus (CN) among the three groups (P < 0.05, FDR correction). Age had a significant positive impact on the susceptibility values of GP (p = 0.018), Put (p < 0.001), and CN (p < 0.001). A history of diabetes had a significant positive influence on the susceptibility values of Put (p = 0.011) and CN (p < 0.001). A smoking history had a significant positive association with the susceptibility values of CN (p = 0.019). Mediation effect analysis demonstrated that iron deposition in the neostriatum partially mediated the relationship between hypertension and cognitive function. Age, diabetes, and smoking may increase iron deposition in the basal ganglia, associated with cognitive decline. The mean susceptibility values of the neostriatum played a mediating role in the association between hypertension and cognitive scores.

CONCLUSIONS

Age, diabetes, and smoking are associated with increased iron deposition in the basal ganglia and also linked to cognitive decline. This can help with understanding CSVD and its prevention and treatment.

Exploring the links among brain iron accumulation, cognitive performance, and dietary intake in older adults: A longitudinal MRI study

This study evaluated longitudinal brain iron accumulation in older adults, its association with cognition, and the role of specific nutrients in mitigating iron accumulation. MRI-based, quantitative susceptibility mapping estimates of brain iron concentration were acquired from seventy-two healthy older adults (47 women, ages 60-86) at a baseline timepoint (TP1) and a follow-up timepoint (TP2) 2.5-3.0 years later. Dietary intake was evaluated at baseline using a validated questionnaire. Cognitive performance was assessed at TP2 using the uniform data set (Version 3) neuropsychological tests of episodic memory (MEM) and executive function (EF). Voxel-wise, linear mixed-effects models, adjusted for longitudinal gray matter volume alterations, age, and several non-dietary lifestyle factors revealed brain iron accumulation in multiple subcortical and cortical brain regions, which was negatively associated with both MEM and EF performance at T2. However, consumption of specific dietary nutrients at TP1 was associated with reduced brain iron accumulation. Our study provides a map of brain regions showing iron accumulation in older adults over a short 2.5-year follow-up and indicates that certain dietary nutrients may slow brain iron accumulation.

XSIM: A structural similarity index measure optimized for MRI QSM

PURPOSE

The structural similarity index measure (SSIM) has become a popular quality metric to evaluate QSM in a way that is closer to human perception than RMS error (RMSE). However, SSIM may overpenalize errors in diamagnetic tissues and underpenalize them in paramagnetic tissues, resulting in biasing. In addition, extreme artifacts may compress the dynamic range, resulting in unrealistically high SSIM scores (hacking). To overcome biasing and hacking, we propose XSIM: SSIM implemented in the native QSM range, and with internal parameters optimized for QSM.

METHODS

We used forward simulations from a COSMOS ground-truth brain susceptibility map included in the 2016 QSM Reconstruction Challenge to investigate the effect of QSM reconstruction errors on the SSIM, XSIM, and RMSE metrics. We also used these metrics to optimize QSM reconstructions of the in vivo challenge data set. We repeated this experiment with the QSM abdominal phantom. To validate the use of XSIM instead of SSIM for QSM quality assessment across a range of different reconstruction techniques/algorithms, we analyzed the reconstructions submitted to the 2019 QSM Reconstruction Challenge 2.0.

RESULTS

Our experiments confirmed the biasing and hacking effects on the SSIM metric applied to QSM. The XSIM metric was robust to those effects, penalizing the presence of streaking artifacts and reconstruction errors. Using XSIM to optimize QSM reconstruction regularization weights returned less overregularization than SSIM and RMSE.

CONCLUSION

XSIM is recommended over traditional SSIM to evaluate QSM reconstructions against a known ground truth, as it avoids biasing and hacking effects and provides a larger dynamic range of scores.

Trajectories and sex differences of brain structure, oxygenation and perfusion functions in normal aging

BACKGROUND

Brain structure, oxygenation and perfusion are important factors in aging. Coupling between regional cerebral oxygen consumption and perfusion also reflects functions of neurovascular unit (NVU). Their trajectories and sex differences during normal aging important for clinical interpretation are still not well defined. In this study, we aim to investigate the relationship between brain structure, functions and age, and exam the sex disparities.

METHOD

A total of 137 healthy subjects between 20∼69 years old were enrolled with conventional MRI, structural three-dimensional T1-weighted imaging (3D-T1WI), 3D multi-echo gradient echo sequence (3D-mGRE), and 3D pseudo-continuous arterial spin labeling (3D-pCASL). Oxygen extraction fraction (OEF) and cerebral blood flow (CBF) were respectively reconstructed from 3D-mGRE and 3D-pCASL images. Cerebral metabolic rate of oxygen (CMRO2) were calculated as follows: CMRO2=CBF·OEF·[H]a, [H]a=7.377 μmol/mL. Brains were segmented into global gray matter (GM), global white matter (WM), and 148 cortical subregions. OEF, CBF, CMRO2, and volumes of GM/WM relative to intracranial volumes (rel_GM/rel_WM) were compared between males and females. Generalized additive models were used to evaluate the aging trajectories of brain structure and functions. The coupling between OEF and CBF was analyzed by correlation analysis. P or PFDR < 0.05 was considered statistically significant.

RESULTS

Females had larger rel_GM, higher CMRO2 and CBF of GM/WM than males (P < 0.05). With control of sex, CBF of GM significantly declined between 20 and 32 years, CMRO2 of GM declined subsequently from 33 to 41 years and rel_GM decreased significantly at all ages (R2 = 0.27, P < 0.001; R2 = 0.17, P < 0.001; R2 = 0.52, P < 0.001). In subregion analysis, CBF declined dispersedly while CMRO2 declined widely across most subregions of the cortex during aging. Robust negative coupling between OEF and CBF was found in most of the subregions (r range = -0.12∼-0.48, PFDR < 0.05).

CONCLUSION

The sex disparities, age trajectories of brain structure and functions as well as the coupling of NVU in healthy individuals provide insights into normal aging which are potential targets for study of pathological conditions.

Reliability of quantitative magnetic susceptibility imaging metrics for cerebral cortex and major subcortical structures

BACKGROUND AND PURPOSE

Susceptibility estimates derived from quantitative susceptibility mapping (QSM) images for the cerebral cortex and major subcortical structures are variably reported in brain magnetic resonance imaging (MRI) studies, as average of all (μ all ${{{{\mu}}}_{{\mathrm{all}}}}$ ), absolute (μ abs ${{{{\mu}}}_{{\mathrm{abs}}}}$ ), or positive- (μ p ${{{{\mu}}}_{\mathrm{p}}}$ ) and negative-only (μ n ${{{{\mu}}}_{\mathrm{n}}}$ ) susceptibility values using a region of interest (ROI) approach. This pilot study presents a reliability analysis of currently used ROI-QSM metrics and an alternative ROI-based approach to obtain voxel-weighted ROI-QSM metrics (μ wp ${{{{\mu}}}_{{\mathrm{wp}}}}$ andμ wn ${{{{\mu}}}_{{\mathrm{wn}}}}$ ).

METHODS

Ten healthy subjects underwent repeated (test-retest) 3-dimensional multi-echo gradient-echo (3DMEGE) 3 Tesla MRI measurements. Complex-valued 3DMEGE images were acquired and reconstructed with slice thicknesses of 1 and 2 mm (3DMEGE1, 3DMEGE2) along with 3DT1-weighted isometric (voxel 1 mm3) images for independent registration and ROI segmentation. Agreement, consistency, and reproducibility of ROI-QSM metrics were assessed through Bland-Altman analysis, intraclass correlation coefficient, and interscan and intersubject coefficient of variation (CoV).

RESULTS

All ROI-QSM metrics exhibited good to excellent consistency and test-retest agreement with no proportional bias. Interscan CoV was higher forμ all ${{{{\mu}}}_{{\mathrm{all}}}}$ in comparison to the other metrics where it was below 15%, in both 3DMEGE1 and 3DMEGE2 datasets. Intersubject CoV forμ all ${{{{\mu}}}_{{\mathrm{all}}}}$ andμ abs ${{{{\mu}}}_{{\mathrm{abs}}}}$ exceeded 50% in all ROIs.

CONCLUSIONS

Among the evaluated ROI-QSM metrics,μ all ${{{{\mu}}}_{{\mathrm{all}}}}$ andμ abs ${{{{\mu}}}_{{\mathrm{abs}}}}$ estimates were less reliable, whereas separating positive and negative values (usingμ p , μ n , μ wp , μ wn ${{{{\mu}}}_{\mathrm{p}}},\ {{{{\mu}}}_{\mathrm{n}}},\ {{{{\mu}}}_{{\mathrm{wp}}}},\ {{{{\mu}}}_{{\mathrm{wn}}}}$ ) improved the reproducibility within, and the comparability between, subjects, even when reducing the slice thickness. These preliminary findings may offer valuable insights toward standardizing ROI-QSM metrics across different patient cohorts and imaging settings in future clinical MRI studies.

SIPAS: A comprehensive susceptibility imaging process and analysis studio

Quantitative susceptibility mapping (QSM) is a rising MRI-based technology and quite a few QSM-related algorithms have been proposed to reconstruct maps of tissue susceptibility distribution from phase images. In this paper, we develop a comprehensive susceptibility imaging process and analysis studio (SIPAS) that can accomplish reliable QSM processing and offer a standardized evaluation system. Specifically, SIPAS integrates multiple methods for each step, enabling users to select algorithm combinations according to data conditions, and QSM maps could be evaluated by two aspects, including image quality indicators within all voxels and region-of-interest (ROI) analysis. Through a sophisticated design of user-friendly interfaces, the results of each procedure are able to be exhibited in axial, coronal, and sagittal views in real-time, meanwhile ROIs can be displayed in 3D rendering visualization. The accuracy and compatibility of SIPAS are demonstrated by experiments on multiple in vivo human brain datasets acquired from 3T, 5T, and 7T MRI scanners of different manufacturers. We also validate the QSM maps obtained by various algorithm combinations in SIPAS, among which the combination of iRSHARP and SFCR achieves the best results on its evaluation system. SIPAS is a comprehensive, sophisticated, and reliable toolkit that may prompt the QSM application in scientific research and clinical practice.

The changes of oxygen extraction fraction in different types of lesions in relapsing-remitting multiple sclerosis: A cross-sectional and longitudinal study

OBJECTIVES

To explore the oxygen metabolism level of different types of lesions in relapsing-remitting multiple sclerosis (RRMS) patients by oxygen extraction fraction (OEF) both cross-sectionally and longitudinally.

METHODS

Forty-six RRMS patients and forty-one healthy controls (HC) went MRI examination. The quantitative susceptibility mapping (QSM) and OEF map were reconstructed from a 3D multi-echo gradient echo sequence. MS lesions in white matter were classified as contrast-enhancing lesions (CELs) on post-gadolinium T1-weighted sequence, paramagnetic rim lesions (PRLs), hyperintense lesions and non-hyperintense lesions on QSM, respectively. The susceptibility and OEF of different types of lesions were compared. The susceptibility and OEF values were measured and compared among different types of lesions. Among these RRMS patients, seventeen had follow-up MRI and 232 lesions, and baseline to follow-up longitudinal changes in susceptibility and OEF were measured.

RESULTS

PRLs had higher susceptibility and lower OEF than CELs, hyperintense lesions, and non-hyperintense lesions. The hyperintense lesions had higher susceptibility and lower OEF than non-hyperintense lesions. In longitudinal changes, PRLs had susceptibility increased (P < 0.001) and OEF decreased (P < 0.001). The hyperintense lesions showed significant decreases in susceptibility (P = 0.020), and non-hyperintense lesions showed significant increases in OEF during follow-up (P = 0.005). Notably, hyperintense lesions may convert to PRLs or non-hyperintense lesions as time progresses, accompanied by changes of OEF and susceptibility in the lesions.

CONCLUSION

This study revealed tissue damage and oxygen metabolism level in different types of MS lesions. The OEF may contribute to further understanding the evolution of MS lesions.

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.

Maximum spherical mean value filtering for whole-brain QSM

PURPOSE

To develop a tissue field-filtering algorithm, called maximum spherical mean value (mSMV), for reducing shadow artifacts in QSM of the brain without requiring brain-tissue erosion.

THEORY AND METHODS

Residual background field is a major source of shadow artifacts in QSM. The mSMV algorithm filters large field-magnitude values near the border, where the maximum value of the harmonic background field is located. The effectiveness of mSMV for artifact removal was evaluated by comparing existing QSM algorithms in numerical brain simulation as well as using in vivo human data acquired from 11 healthy volunteers and 93 patients.

RESULTS

Numerical simulation showed that mSMV reduces shadow artifacts and improves QSM accuracy. Better shadow reduction, as demonstrated by lower QSM variation in the gray matter and higher QSM image quality score, was also observed in healthy subjects and in patients with hemorrhages, stroke, and multiple sclerosis.

CONCLUSION

The mSMV algorithm allows QSM maps that are substantially equivalent to those obtained using SMV-filtered dipole inversion without eroding the volume of interest.

Free-breathing high isotropic resolution quantitative susceptibility mapping (QSM) of liver using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction

PURPOSE

To enable free-breathing and high isotropic resolution liver quantitative susceptibility mapping (QSM) using 3D multi-echo UTE cones acquisition and respiratory motion-resolved image reconstruction.

METHODS

Using 3D multi-echo UTE cones MRI, a respiratory motion was estimated from the k-space center of the imaging data. After sorting the k-space data with estimated motion, respiratory motion state-resolved reconstruction was performed for multi-echo data followed by nonlinear least-squares fitting for proton density fat fraction (PDFF),R 2 * $$ {\mathrm{R}}_2^{\ast } $$ , and fat-corrected B0 field maps. PDFF and B0 field maps were subsequently used for QSM reconstruction. The proposed method was compared with motion-averaged (gridding) reconstruction and conventional 3D multi-echo Cartesian MRI in moving gadolinium phantom and in vivo studies. Region of interest (ROI)-based linear regression analysis was performed on these methods to investigate correlations between gadolinium concentration and QSM in the phantom study and betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM in in vivo study.

RESULTS

Cones with motion-resolved reconstruction showed sharper image quality compared to motion-averaged reconstruction with a substantial reduction of motion artifacts in both moving phantom and in vivo studies. For ROI-based linear regression analysis of the phantom study, susceptibility values from cones with motion-resolved reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.31 × gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.05,R 2 $$ {R}^2 $$  = 0.999) and Cartesian without motion (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.32× gadolinium mM + $$ \times {\mathrm{gadolinium}}_{\mathrm{mM}}+ $$ 0.04,R 2 $$ {R}^2 $$  = 1.000) showed linear relationships with gadolinium concentrations and showed good agreement with each other. For in vivo, motion-resolved reconstruction showed higher goodness of fit (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$  = 0.00261 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.524,R 2 $$ {R}^2 $$  = 0.977) compared to motion-averaged reconstruction (QSM ppm $$ {\mathrm{QSM}}_{\mathrm{ppm}} $$ = 0.0021 × R 2 s - 1 * - $$ \times {\mathrm{R}}_{2_{{\mathrm{s}}^{-1}}}^{\ast }- $$ 0.572,R 2 $$ {R}^2 $$ = 0.723) in ROI-based linear regression analysis betweenR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and QSM.

CONCLUSION

Feasibility of free-breathing liver QSM was demonstrated with motion-resolved 3D multi-echo UTE cones MRI, achieving high isotropic resolution currently unachievable in conventional Cartesian MRI.

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.

Data adaptive regularization with reference tissue constraints for liver quantitative susceptibility mapping

PURPOSE

To improve repeatability and reproducibility across acquisition parameters and reduce bias in quantitative susceptibility mapping (QSM) of the liver, through development of an optimized regularized reconstruction algorithm for abdominal QSM.

METHODS

An optimized approach to estimation of magnetic susceptibility distribution is formulated as a constrained reconstruction problem that incorporates estimates of the input data reliability and anatomical priors available from chemical shift-encoded imaging. The proposed data-adaptive method was evaluated with respect to bias, repeatability, and reproducibility in a patient population with a wide range of liver iron concentration (LIC). The proposed method was compared to the previously proposed and validated approach in liver QSM for two multi-echo spoiled gradient-recalled echo protocols with different acquisition parameters at 3T. Linear regression was used for evaluation of QSM methods against a reference FDA-approvedR 2 $$ {R}_2 $$ -based LIC measure andR 2 ∗ $$ {R}_2^{\ast } $$ measurements; repeatability/reproducibility were assessed by Bland-Altman analysis.

RESULTS

The data-adaptive method produced susceptibility maps with higher subjective quality due to reduced shading artifacts. For both acquisition protocols, higher linear correlation with bothR 2 $$ {R}_2 $$ - andR 2 ∗ $$ {R}_2^{\ast } $$ -based measurements were observed for the data-adaptive method (r 2 = 0 . 74 / 0 . 69 $$ {r}^2=0.74/0.69 $$ forR 2 $$ {R}_2 $$ ,0 . 97 / 0 . 95 $$ 0.97/0.95 $$ forR 2 ∗ $$ {R}_2^{\ast } $$ ) than the standard method (r 2 = 0 . 60 / 0 . 66 $$ {r}^2=0.60/0.66 $$ and0 . 79 / 0 . 88 $$ 0.79/0.88 $$ ). For both protocols, the data-adaptive method enabled better test-retest repeatability (repeatability coefficients 0.19/0.30 ppm for the data-adaptive method, 0.38/0.47 ppm for the standard method) and reproducibility across protocols (reproducibility coefficient 0.28 vs. 0.53ppm) than the standard method.

CONCLUSIONS

The proposed data-adaptive QSM algorithm may enable quantification of LIC with improved repeatability/reproducibility across different acquisition parameters as 3T.

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

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.

Deep Learning-Assisted Quantitative Susceptibility Mapping as a Tool for Grading and Molecular Subtyping of Gliomas

This study aimed to explore the value of deep learning (DL)-assisted quantitative susceptibility mapping (QSM) in glioma grading and molecular subtyping. Forty-two patients with gliomas, who underwent preoperative T2 fluid-attenuated inversion recovery (T2 FLAIR), contrast-enhanced T1-weighted imaging (T1WI + C), and QSM scanning at 3.0T magnetic resonance imaging (MRI) were included in this study. Histopathology and immunohistochemistry staining were used to determine glioma grades, and isocitrate dehydrogenase (IDH) 1 and alpha thalassemia/mental retardation syndrome X-linked gene (ATRX) subtypes. Tumor segmentation was performed manually using Insight Toolkit-SNAP program (www.itksnap.org). An inception convolutional neural network (CNN) with a subsequent linear layer was employed as the training encoder to capture multi-scale features from MRI slices. Fivefold cross-validation was utilized as the training strategy (seven samples for each fold), and the ratio of sample size of the training, validation, and test dataset was 4:1:1. The performance was evaluated by the accuracy and area under the curve (AUC). With the inception CNN, single modal of QSM showed better performance in differentiating glioblastomas (GBM) and other grade gliomas (OGG, grade II-III), and predicting IDH1 mutation and ATRX loss (accuracy: 0.80, 0.77, 0.60) than either T2 FLAIR (0.69, 0.57, 0.54) or T1WI + C (0.74, 0.57, 0.46). When combining three modalities, compared with any single modality, the best AUC/accuracy/F1-scores were reached in grading gliomas (OGG and GBM: 0.91/0.89/0.87, low-grade and high-grade gliomas: 0.83/0.86/0.81), predicting IDH1 mutation (0.88/0.89/0.85), and predicting ATRX loss (0.78/0.71/0.67). As a supplement to conventional MRI, DL-assisted QSM is a promising molecular imaging method to evaluate glioma grades, IDH1 mutation, and ATRX loss.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-022-00087-6.

Iron quantitative analysis of motor combined with bulbar region in M1 cortex may improve diagnosis performance in ALS

OBJECTIVES

To explore whether the combined analysis of motor and bulbar region of M1 on susceptibility-weighted imaging (SWI) can be a valid biomarker for amyotrophic lateral sclerosis (ALS).

METHODS

Thirty-two non-demented ALS patients and 35 age- and gender-matched healthy controls (HC) were retrospectively recruited. SWI and 3D-T1-MPRAGE images were obtained from all individuals using a 3.0-T MRI scan. The bilateral posterior band of M1 was manually delineated by three neuroradiologists on phase images and subdivided into the motor and bulbar regions. We compared the phase values in two groups and performed a stratification analysis (ALSFRS-R score, duration, disease progression rate, and onset). Receiver operating characteristic (ROC) curves were also constructed.

RESULTS

ALS group showed significantly increased phase values in M1 and the two subregions than the HC group, on the all and elderly level (p < 0.001, respectively). On all-age level comparison, negative correlations were found between phase values of M1 and clinical score and duration (p < 0.05, respectively). Similar associations were found in the motor region (p < 0.05, respectively). On both the total (p < 0.01) and elderly (p < 0.05) levels, there were positive relationships between disease progression rate and M1 phase values. In comparing ROC curves, the entire M1 showed the best diagnostic performance.

CONCLUSIONS

Combining motor and bulbar analyses as an integral M1 region on SWI can improve ALS diagnosis performance, especially in the elderly. The phase value could be a valuable biomarker for ALS evaluation.

KEY POINTS

• Integrated analysis of the motor and bulbar as an entire M1 region on SWI can improve the diagnosis performance in ALS. • Quantitative analysis of iron deposition by SWI measurement helps the clinical evaluation, especially for the elderly patients. • Phase value, when combined with the disease progression rate, could be a valuable biomarker for ALS.

Incomplete spectrum QSM using support information

Introduction

Reconstructing a bounded object from incomplete k-space data is a well posed problem, and it was recently shown that this incomplete spectrum approach can be used to reconstruct undersampled MRI images with similar quality to compressed sensing approaches. Here, we apply this incomplete spectrum approach to the field-to-source inverse problem encountered in quantitative magnetic susceptibility mapping (QSM). The field-to-source problem is an ill-posed problem because of conical regions in frequency space where the dipole kernel is zero or very small, which leads to the kernel's inverse being ill-defined. These "ill-posed" regions typically lead to streaking artifacts in QSM reconstructions. In contrast to compressed sensing, our approach relies on knowledge of the image-space support, more commonly referred to as the mask, of our object as well as the region in k-space with ill-defined values. In the QSM case, this mask is usually available, as it is required for most QSM background field removal and reconstruction methods.

Methods

We tuned the incomplete spectrum method (mask and band-limit) for QSM on a simulated dataset from the most recent QSM challenge and validated the QSM reconstruction results on brain images acquired in five healthy volunteers, comparing incomplete spectrum QSM to current state-of-the art-methods: FANSI, nonlinear dipole inversion, and conventional thresholded k-space division.

Results

Without additional regularization, incomplete spectrum QSM performs slightly better than direct QSM reconstruction methods such as thresholded k-space division (PSNR of 39.9 vs. 39.4 of TKD on a simulated dataset) and provides susceptibility values in key iron-rich regions similar or slightly lower than state-of-the-art algorithms, but did not improve the PSNR in comparison to FANSI or nonlinear dipole inversion. With added (ℓ1-wavelet based) regularization the new approach produces results similar to compressed sensing based reconstructions (at sufficiently high levels of regularization).

Discussion

Incomplete spectrum QSM provides a new approach to handle the "ill-posed" regions in the frequency-space data input to QSM.

DeepSPIO: Super Paramagnetic Iron Oxide Particle Quantification Using Deep Learning in Magnetic Resonance Imaging

The susceptibility of super paramagnetic iron oxide (SPIO) particles makes them a useful contrast agent for different purposes in MRI. These particles are typically quantified with relaxometry or by measuring the inhomogeneities they produced. These methods rely on the phase, which is unreliable for high concentrations. We present in this study a novel Deep Learning method to quantify the SPIO concentration distribution. We acquired the data with a new sequence called View Line in which the field map information is encoded in the geometry of the image. The novelty of our network is that it uses residual blocks as the bottleneck and multiple decoders to improve the gradient flow in the network. Each decoder predicts a different part of the wavelet decomposition of the concentration map. This decomposition improves the estimation of the concentration, and also it accelerates the convergence of the model. We tested our SPIO concentration reconstruction technique with simulated images and data from actual scans from phantoms. The simulations were done using images from the IXI dataset, and the phantoms consisted of plastic cylinders containing agar with SPIO particles at different concentrations. In both experiments, the model was able to quantify the distribution accurately.

Global cerebrospinal fluid as a zero-reference regularization for brain quantitative susceptibility mapping

BACKGROUND AND PURPOSE

The objective ofthis study was to demonstrate a global cerebrospinal fluid (CSF) method for a consistent and automated zero referencing of brain quantitative susceptibility mapping (QSM).

METHODS

Whole brain CSF mask was automatically segmented by thresholding the gradient echo transverse relaxation ( R2∗) map, and regularization was employed to enforce uniform susceptibility distribution within the CSF volume in the field-to-susceptibility inversion. This global CSF regularization method was compared with a prior ventricular CSF regularization. Both reconstruction methods were compared in a repeatability study of 12 healthy subjects using t-test on susceptibility measurements, and in patient studies of 17 multiple sclerosis (MS) and 10 Parkinson's disease (PD) patients using Wilcoxon rank-sum test on radiological scores.

RESULTS

In scan-rescan experiments, global CSF regularization provided more consistent CSF volume as well as higher repeatability of QSM measurements than ventricular CSF regularization with a smaller bias: -2.7 parts per billion (ppb) versus -0.13 ppb (t-test p<0.05) and a narrower 95% limits of agreement: [-7.25, 6.99] ppb versus [-16.60, 11.19 ppb] (f-test p<0.05). In PD and MS patients, global CSF regularization reduced smoothly varying shadow artifacts and significantly improved the QSM quality score (p<0.001).

CONCLUSIONS

The proposed whole brain CSF method for QSM zero referencing improves repeatability and image quality of brain QSM compared to the ventricular CSF method.

Streaking artifact suppression of quantitative susceptibility mapping reconstructions via L1-norm data fidelity optimization (L1-QSM)

PURPOSE

The presence of dipole-inconsistent data due to substantial noise or artifacts causes streaking artifacts in quantitative susceptibility mapping (QSM) reconstructions. Often used Bayesian approaches rely on regularizers, which in turn yield reduced sharpness. To overcome this problem, we present a novel L1-norm data fidelity approach that is robust with respect to outliers, and therefore prevents streaking artifacts.

METHODS

QSM functionals are solved with linear and nonlinear L1-norm data fidelity terms using functional augmentation, and are compared with equivalent L2-norm methods. Algorithms were tested on synthetic data, with phase inconsistencies added to mimic lesions, QSM Challenge 2.0 data, and in vivo brain images with hemorrhages.

RESULTS

The nonlinear L1-norm-based approach achieved the best overall error metric scores and better streaking artifact suppression. Notably, L1-norm methods could reconstruct QSM images without using a brain mask, with similar regularization weights for different data fidelity weighting or masking setups.

CONCLUSION

The proposed L1-approach provides a robust method to prevent streaking artifacts generated by dipole-inconsistent data, renders brain mask calculation unessential, and opens novel challenging clinical applications such asassessing brain hemorrhages and cortical layers.

Ironsmith: An Automated Pipeline for QSM-based Data Analyses

Quantitative susceptibility mapping (QSM) is an MRI-based, computational method for anatomically localizing and measuring concentrations of specific biomarkers in tissue such as iron. Growing research suggests QSM is a viable method for evaluating the impact of iron overload in neurological disorders and on cognitive performance in aging. Several software toolboxes are currently available to reconstruct QSM maps from 3D GRE MR Images. However, few if any software packages currently exist that offer fully automated pipelines for QSM-based data analyses: from DICOM images to region-of-interest (ROI) based QSM values. Even less QSM-based software exist that offer quality control measures for evaluating the QSM output. Here, we address these gaps in the field by introducing and demonstrating the reliability and external validity of Ironsmith; an open-source, fully automated pipeline for creating and processing QSM maps, extracting QSM values from subcortical and cortical brain regions (89 ROIs) and evaluating the quality of QSM data using SNR measures and assessment of outlier regions on phase images. Ironsmith also features automatic filtering of QSM outlier values and precise CSF-only QSM reference masks that minimize partial volume effects. Testing of Ironsmith revealed excellent intra- and inter-rater reliability. Finally, external validity of Ironsmith was demonstrated via an anatomically selective relationship between motor performance and Ironsmith-derived QSM values in motor cortex. In sum, Ironsmith provides a freely-available, reliable, turn-key pipeline for QSM-based data analyses to support research on the impact of brain iron in aging and neurodegenerative disease.

QQ-NET - using deep learning to solve quantitative susceptibility mapping and quantitative blood oxygen level dependent magnitude (QSM+qBOLD or QQ) based oxygen extraction fraction (OEF) mapping

PURPOSE

To improve accuracy and speed of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) -based oxygen extraction fraction (OEF) mapping using a deep neural network (QQ-NET).

METHODS

The 3D multi-echo gradient echo images were acquired in 34 ischemic stroke patients and 4 healthy subjects. Arterial spin labeling and diffusion weighted imaging (DWI) were also performed in the patients. NET was developed to solve the QQ model inversion problem based on Unet. QQ-based OEF maps were reconstructed with previously introduced temporal clustering, tissue composition, and total variation (CCTV) and NET. The results were compared in simulation, ischemic stroke patients, and healthy subjects using a two-sample Kolmogorov-Smirnov test.

RESULTS

In the simulation, QQ-NET provided more accurate and precise OEF maps than QQ-CCTV with 150 times faster reconstruction speed. In the subacute stroke patients, OEF from QQ-NET had greater contrast-to-noise ratio (CNR) between DWI-defined lesions and their unaffected contralateral normal tissue than with QQ-CCTV: 1.9 ± 1.3 vs 6.6 ± 10.7 (p = 0.03). In healthy subjects, both QQ-CCTV and QQ-NET provided uniform OEF maps.

CONCLUSION

QQ-NET improves the accuracy of QQ-based OEF with faster reconstruction.

The Effect of Deep Learning-Based QSM Magnetic Resonance Imaging on the Subthalamic Nucleus

In order to study the influence of quantitative magnetic susceptibility mapping (QSM) on them. A 2.5D Attention U-Net Network based on multiple input and multiple output, a method for segmenting RN, SN, and STN regions in high-resolution QSM images is proposed, and deep learning realizes accurate segmentation of deep nuclei in brain QSM images. Experimental results show data first cuts each layer of 0 100 case data, based on the image center, from 384 × 288 to the size of 128 × 128. Image combination: each layer of the image in the layer direction combines with two adjacent images into a 2.5D image, i.e., (It - m It; It + i), where It represents the layer i image. At this time, the size of the image changes from 128 × 128 to 128 × 128 × 3, in which 3 represents three consecutive layers of images. The SNR of SWP I to STN is twice that of SWI. The small deep gray matter nuclei (RN, SN, and STN) in QSM images of the brain and the pancreas with irregular shape and large individual differences in abdominal CT images can be automatically segmented.

QSM reconstruction challenge 2.0: Design and report of results

PURPOSE

The aim of the second quantitative susceptibility mapping (QSM) reconstruction challenge (Oct 2019, Seoul, Korea) was to test the accuracy of QSM dipole inversion algorithms in simulated brain data.

METHODS

A two-stage design was chosen for this challenge. The participants were provided with datasets of multi-echo gradient echo images synthesized from two realistic in silico head phantoms using an MR simulator. At the first stage, participants optimized QSM reconstructions without ground truth data available to mimic the clinical setting. At the second stage, ground truth data were provided for parameter optimization. Submissions were evaluated using eight numerical metrics and visual ratings.

RESULTS

A total of 98 reconstructions were submitted for stage 1 and 47 submissions for stage 2. Iterative methods had the best quantitative metric scores, followed by deep learning and direct inversion methods. Priors derived from magnitude data improved the metric scores. Algorithms based on iterative approaches and total variation (and its derivatives) produced the best overall results. The reported results and analysis pipelines have been made public to allow researchers to compare new methods to the current state of the art.

CONCLUSION

The synthetic data provide a consistent framework to test the accuracy and robustness of QSM algorithms in the presence of noise, calcifications and minor voxel dephasing effects. Total Variation-based algorithms produced the best results among all metrics. Future QSM challenges should assess whether this good performance with synthetic datasets translates to more realistic scenarios, where background fields and dipole-incompatible phase contributions are included.

In vivo assessment of anisotropy of apparent magnetic susceptibility in white matter from a single orientation acquisition

Multiple studies have reported a significant dependence of the effective transverse relaxation rate constant (R₂^*) and the phase of gradient-echo based (GRE) signal on the orientation of white matter fibres in the human brain. It has also been hypothesized that magnetic susceptibility, as obtained by single-orientation quantitative susceptibility mapping (QSM), exhibits such a dependence. In this study, we investigated this hypothesized relationship in a cohort of healthy volunteers. We show that R₂^* follows the predicted orientation dependence consistently across white matter regions, whereas the apparent magnetic susceptibility is related differently to fibre orientation across the brain and often in a complex non-monotonic manner. In addition, we explored the effect of fractional anisotropy measured by diffusion-weighted MRI on the strength of the orientation dependence and observed only a limited influence in many regions. However, with careful consideration of such an impact and the limitations imposed by the ill-posed nature of the dipole inversion process, it is possible to study magnetic susceptibility anisotropy in specific brain regions with a single orientation acquisition.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

Single-step calculation of susceptibility through multiple orientation sampling

Quantitative susceptibility mapping (QSM) was developed to estimate the spatial distribution of magnetic susceptibility from MR signal phase acquired using a gradient echo (GRE) sequence. The field-to-susceptibility inversion in QSM is known to be ill-posed and needs numerical stabilization through either regularization or data oversampling. The calculation of susceptibility through the multiple orientation sampling (COSMOS) method uses phase data acquired at three or more head orientations to achieve a well-conditioned field-to-susceptibility inversion and is often considered the gold standard for in vivo QSM. However, the conventional COSMOS approach, here named multistep COSMOS (MSCOSMOS), solves the dipole inversion from the local field derived from raw GRE phase through multiple steps of phase preprocessing. Error propagations between these consecutive phase processing steps can thus affect the final susceptibility quantification. On the other hand, recently proposed single-step QSM (SSQSM) methods aim to solve an integrated inversion from unprocessed or total phase to mitigate such error propagations but have been limited to single orientation QSM. This study therefore aimed to test the feasibility of using single-step COSMOS (SSCOSMOS) to jointly perform background field removal and dipole inversion with multiple orientation sampling, which could serve as a better standard for gauging SSQSM methods. We incorporated multiple spherical mean value (SMV) kernels of various radii with the dipole inversion in SSCOSMOS. QSM reconstructions with SSCOSMOS and MSCOSMOS were compared using both simulations with a numerical head phantom and in vivo human brain data. SSCOSMOS permitted integrated background removal and dipole inversion without the need to adjust any regularization parameters. In addition, with sufficiently large SMV kernels, SSCOSMOS performed consistently better than MSCOSMOS in all the tested error metrics in our simulations, giving better susceptibility quantification and smaller reconstruction error. Consistent tissue susceptibility values were obtained between SSCOSMOS and MSCOSMOS.

Multiecho complex total field inversion method (mcTFI) for improved signal modeling in quantitative susceptibility mapping

PURPOSE

Typical quantitative susceptibility mapping (QSM) reconstruction steps consist of first estimating the magnetization field from the gradient-echo images, and then reconstructing the susceptibility map from the estimated field. The errors from the field-estimation steps may propagate into the final QSM map, and the noise in the estimated field map may no longer be zero-mean Gaussian noise, thus, causing streaking artifacts in the resulting QSM. A multiecho complex total field inversion (mcTFI) method was developed to compute the susceptibility map directly from the multiecho gradient echo images using an improved signal model that retains the Gaussian noise property in the complex domain. It showed improvements in QSM reconstruction over the conventional field-to-source inversion.

METHODS

The proposed mcTFI method was compared with the nonlinear total field inversion (nTFI) method in a numerical brain with hemorrhage and calcification, the numerical brains provided by the QSM Challenge 2.0, 18 brains with intracerebral hemorrhage scanned at 3T, and 6 healthy brains scanned at 7T.

RESULTS

Compared with nTFI, the proposed mcTFI showed more accurate QSM reconstruction around the lesions in the numerical simulations. The mcTFI reconstructed QSM also showed the best image quality with the least artifacts in the brains with intracerebral hemorrhage scanned at 3T and healthy brains scanned at 7T.

CONCLUSION

The proposed multiecho complex total field inversion improved QSM reconstruction over traditional field-to-source inversion through better signal modeling.

Computed inverse MRI (CIMRI) for intrinsic brain magnetic susceptibility mapping

In magnetic resonance imaging (MRI), tissue magnetization in the main field B₀ is a necessary preparation for magnetic resonance signal formation that imposes an inherent dipole effect on MRI signals, which predisposes an artifact on tissue MRI. In the MRI principle, T2*-weighted MRI can be described by a cascade of data transformations: from the source of tissue magnetic susceptibility (denoted by χ) to the output of complex-valued T2* image (in a magnitude and phase pair). Under the linear approximation of the T2* phase MRI, we can computationally reconstruct the source χ by quantitative susceptibility mapping (QSM), which is an inverse solution that is modeled by computed inverse MRI (CIMRI). For a brain function study using MRI (fMRI), we can reconstruct a timeseries of brain χ images to represent the intrinsic brain function activity called functional QSM (fQSM). This intrinsic depiction is defined as the removal of the artifactual dipole effect and other MRI-introduced distortions from phase data through inverse mapping. With one high-resolution QSM experiment and one group (20 subjects) low-resolution fQSM experiment, we show that the dipole effect manifests as ripples around vessels and a spatial split at a local activation blob and that the dipole effect could be removed by CIMRI. In the context of inverse imaging or undoing MRI transformations (including dipole convolution), we computationally achieve brain intrinsic structural depiction by QSM and intrinsic functional depiction by fQSM.

QSMART: Quantitative susceptibility mapping artifact reduction technique

PURPOSE

Quantitative susceptibility mapping (QSM) is a novel MR technique that allows mapping of tissue susceptibility values from MR phase images. QSM is an ill-conditioned inverse problem, and although several methods have been proposed in the field, in the presence of a wide range of susceptibility sources, streaking artifacts appear around high susceptibility regions and contaminate the whole QSM map. QSMART is a post-processing pipeline that uses two-stage parallel inversion to reduce the streaking artifacts and remove banding artifact at the cortical surface and around the vasculature.

METHOD

Tissue and vein susceptibility values were separately estimated by generating a mask of vasculature driven from the magnitude data using a Frangi filter. Spatially dependent filtering was used for the background field removal step and the two susceptibility estimates were combined in the final QSM map. QSMART was compared to RESHARP/iLSQR and V-SHARP/iLSQR inversion in a numerical phantom, 7T in vivo single and multiple-orientation scans, 9.4T ex vivo mouse data, and 4.7T in vivo rat brain with induced focal ischemia.

RESULTS

Spatially dependent filtering showed better suppression of phase artifacts near cortex compared to RESHARP and V-SHARP, while preserving voxels located within regions of interest without brain edge erosion. QSMART showed successful reduction of streaking artifacts as well as improved contrast between different brain tissues compared to the QSM maps obtained by RESHARP/iLSQR and V-SHARP/iLSQR.

CONCLUSION

QSMART can reduce QSM artifacts to enable more robust estimation of susceptibility values in vivo and ex vivo.

Comparison of parameter optimization methods for quantitative susceptibility mapping

PURPOSE

Quantitative Susceptibility Mapping (QSM) is usually performed by minimizing a functional with data fidelity and regularization terms. A weighting parameter controls the balance between these terms. There is a need for techniques to find the proper balance that avoids artifact propagation and loss of details. Finding the point of maximum curvature in the L-curve is a popular choice, although it is slow, often unreliable when using variational penalties, and has a tendency to yield overregularized results.

METHODS

We propose 2 alternative approaches to control the balance between the data fidelity and regularization terms: 1) searching for an inflection point in the log-log domain of the L-curve, and 2) comparing frequency components of QSM reconstructions. We compare these methods against the conventional L-curve and U-curve approaches.

RESULTS

Our methods achieve predicted parameters that are better correlated with RMS error, high-frequency error norm, and structural similarity metric-based parameter optimizations than those obtained with traditional methods. The inflection point yields less overregularization and lower errors than traditional alternatives. The frequency analysis yields more visually appealing results, although with larger RMS error.

CONCLUSION

Our methods provide a robust parameter optimization framework for variational penalties in QSM reconstruction. The L-curve-based zero-curvature search produced almost optimal results for typical QSM acquisition settings. The frequency analysis method may use a 1.5 to 2.0 correction factor to apply it as a stand-alone method for a wider range of signal-to-noise-ratio settings. This approach may also benefit from fast search algorithms such as the binary search to speed up the process.

Spatially Adaptive Regularization in Total Field Inversion for Quantitative Susceptibility Mapping

Adaptive Total Field Inversion is described for quantitative susceptibility mapping (QSM) reconstruction from total field data through a spatially adaptive suppression of shadow artifacts through spatially adaptive regularization. The regularization for shadow suppression consists of penalizing low-frequency components of susceptibility in regions of small susceptibility contrasts as estimated by R2∗ derived signal intensity. Compared with a conventional local field method and two previously proposed regularized total field inversion methods, improvements were demonstrated in phantoms and subjects without and with hemorrhages. This algorithm, named TFIR, demonstrates the lowest error in numerical and gadolinium phantom datasets. In COSMOS data, TFIR performs well in matching ground truth in high-susceptibility regions. For patient data, TFIR comes close to meeting the quality of the reference local field method and outperforms other total field techniques in both clinical scores and shadow reduction.

•

TFIR's adaptive regularization obtains magnetic susceptibility from magnetic field

•

TFIR has low artifact incidence on both quantitative and clinical scores

•

The error for TFIR is low on various numerical and ground truth tests

•

Clinical applications for TFIR include hemorrhages and whole head mapping

The 2016 QSM Challenge: Lessons learned and considerations for a future challenge design

PURPOSE

The 4th International Workshop on MRI Phase Contrast and QSM (2016, Graz, Austria) hosted the first QSM Challenge. A single-orientation gradient recalled echo acquisition was provided, along with COSMOS and the χ33 STI component as ground truths. The submitted solutions differed more than expected depending on the error metric used for optimization and were generally over-regularized. This raised (unanswered) questions about the ground truths and the metrics utilized.

METHODS

We investigated the influence of background field remnants by applying additional filters. We also estimated the anisotropic contributions from the STI tensor to the apparent susceptibility to amend the χ33 ground truth and to investigate the impact on the reconstructions. Lastly, we used forward simulations from the COSMOS reconstruction to investigate the impact noise had on the metric scores.

RESULTS

Reconstructions compared against the amended STI ground truth returned lower errors. We show that the background field remnants had a minor impact in the errors. In the absence of inconsistencies, all metrics converged to the same regularization weights, whereas structural similarity index metric was more insensitive to such inconsistencies.

CONCLUSION

There was a mismatch between the provided data and the ground truths due to the presence of unaccounted anisotropic susceptibility contributions and noise. Given the lack of reliable ground truths when using in vivo acquisitions, simulations are suggested for future QSM Challenges.

Quantitative Measurement of Metal Accumulation in Brain of Patients With Wilson's Disease

BACKGROUND

Currently, no study has evaluated metal accumulation in the brains of patients with Wilson's disease by using quantitative susceptibility mapping at 3T MRI. The objectives of this study were to qualitatively and quantitatively evaluate changes in magnetic susceptibility and R2* maps in deep gray matter nuclei to discriminate Wilson's disease patients from healthy controls and to evaluate their sensitivities in diagnosing Wilson's disease.

METHODS

Magnetic susceptibility and R2* maps and conventional T1-weighted, T2-weighted, and T2-weighted fluid-attenuated inversion recovery images were obtained from 17 Wilson's disease patients and 14 age-matched healthy controls on a 3T MRI scanner. Differences between Wilson's disease and healthy control groups in susceptibility and R2* values in multiple deep nuclei were evaluated using a Mann-Whitney U test and receiver operating characteristic curves. The correlations of susceptibility and R2* values with Unified Wilson's Disease Rating Scale score were also performed.

RESULTS

Magnetic susceptibility and R2* can effectively distinguish different types of signal abnormalities. Magnetic susceptibility and R2* values in multiple deep nuclei of Wilson's disease patients were significantly higher than those in healthy controls. Magnetic susceptibility value in the substantia nigra had the highest area under the curve (0.888). There were positive correlations of the Unified Wilson's Disease Rating Scale score with susceptibility values in the caudate nucleus (r = 0.757, P = 0.011), putamen (r = 0.679, P = 0.031), and red nucleus (r = 0.638, P = 0.047), as well as R2* values in the caudate nucleus (r = 0.754, P = 0.012).

CONCLUSIONS

Quantitative susceptibility mapping at 3T could be a useful tool to evaluate metal accumulation in deep gray matter nuclei of Wilson's disease patients. © 2020 International Parkinson and Movement Disorder Society.

Fidelity imposed network edit (FINE) for solving ill-posed image reconstruction

Deep learning (DL) is increasingly used to solve ill-posed inverse problems in medical imaging, such as reconstruction from noisy and/or incomplete data, as DL offers advantages over conventional methods that rely on explicit image features and hand engineered priors. However, supervised DL-based methods may achieve poor performance when the test data deviates from the training data, for example, when it has pathologies not encountered in the training data. Furthermore, DL-based image reconstructions do not always incorporate the underlying forward physical model, which may improve performance. Therefore, in this work we introduce a novel approach, called fidelity imposed network edit (FINE), which modifies the weights of a pre-trained reconstruction network for each case in the testing dataset. This is achieved by minimizing an unsupervised fidelity loss function that is based on the forward physical model. FINE is applied to two important inverse problems in neuroimaging: quantitative susceptibility mapping (QSM) and under-sampled image reconstruction in MRI. Our experiments demonstrate that FINE can improve reconstruction accuracy.

Automated adaptive preconditioner for quantitative susceptibility mapping

PURPOSE

To develop an automated adaptive preconditioner for QSM reconstruction with improved susceptibility quantification accuracy and increased image quality.

THEORY AND METHODS

The total field was used to rapidly produce an approximate susceptibility map, which was then averaged and trended over R 2 ∗ binning to generate a spatially varying distribution of preconditioning values. This automated adaptive preconditioner was used to reconstruct QSM via total field inversion and was compared with its empirical counterparts in a numerical simulation, a brain experiment with 5 healthy subjects and 5 patients with intracerebral hemorrhage, and a cardiac experiment with 3 healthy subjects.

RESULTS

Among evaluated preconditioners, the automated adaptive preconditioner achieved the fastest convergence in reducing the RMSE of the QSM in the simulation, suppressed hemorrhage-associated artifacts while preserving surrounding brain tissue contrasts, and provided cardiac chamber oxygenation values consistent with those reported in the literature.

CONCLUSION

An automated adaptive preconditioner allows high-quality QSM from the total field in imaging various anatomies with dynamic susceptibility ranges.

Quantitative susceptibility mapping of the spine using in-phase echoes to initialize inhomogeneous field and R2* for the nonconvex optimization problem of fat-water separation

Quantitative susceptibility mapping (QSM) of human spinal vertebrae from a multi-echo gradient-echo (GRE) sequence is challenging, because comparable amounts of fat and water in the vertebrae make it difficult to solve the nonconvex optimization problem of fat-water separation (R2*-IDEAL) for estimating the magnetic field induced by tissue susceptibility. We present an in-phase (IP) echo initialization of R2*-IDEAL for QSM in the spinal vertebrae. Ten healthy human subjects were recruited for spine MRI. A 3D multi-echo GRE sequence was implemented to acquire out-phase and IP echoes. For the IP method, the R2* and field maps estimated by separately fitting the magnitude and phase of IP echoes were used to initialize gradient search R2*-IDEAL to obtain final R2*, field, water, and fat maps, and the final field map was used to generate QSM. The IP method was compared with the existing Zero method (initializing the field to zero), VARPRO-GC (variable projection using graphcuts but still initializing the field to zero), and SPURS (simultaneous phase unwrapping and removal of chemical shift using graphcuts for initialization) on both simulation and in vivo data. The single peak fat model was also compared with the multi-peak fat model. There was no substantial difference on QSM between the single peak and multi-peak fat models, but there were marked differences among different initialization methods. The simulations demonstrated that IP provided the lowest error in the field map. Compared to Zero, VARPRO-GC and SPURS, the proposed IP method provided substantially improved spine QSM in all 10 subjects.

Clinical Integration of Automated Processing for Brain Quantitative Susceptibility Mapping: Multi-Site Reproducibility and Single-Site Robustness

BACKGROUND AND PURPOSE

Quantitative susceptibility mapping (QSM) of the brain has become highly reproducible and has applications in an expanding array of diseases. To translate QSM from bench to bedside, it is important to automate its reconstruction immediately after data acquisition. In this work, a server system that automatically reconstructs QSM and exchange images with the scanner using the DICOM standard is demonstrated using a multi-site, multi-vendor reproducibility study and a large, single-site, multi-scanner image quality review study in a clinical environment.

METHODS

A single healthy subject was scanned with a 3D multi-echo gradient echo sequence at nine sites around the world using scanners from three manufacturers. A high-resolution (HiRes, .5 × .5 × 1 mm³ reconstructed) and standard-resolution (StdRes, .5 × .5 × 3 mm³ ) protocol was performed. ROI analysis of various white matter and gray matter regions was performed to investigate reproducibility across sites. At one institution, a retrospective multi-scanner image quality review was carried out of all clinical QSM images acquired consecutively in 1 month.

RESULTS

Reconstruction times using a GPU were 29 ± 22 seconds (StdRes) and 55 ± 39 seconds (HiRes). ROI standard deviation across sites was below 24 ppb (StdRes) and 17 ppb (HiRes). Correlations between ROI averages across sites were on average .92 (StdRes) and .96 (HiRes). Image quality review of 873 consecutive patients revealed diagnostic or excellent image quality in 96% of patients.

CONCLUSION

Online QSM reconstruction for a variety of sites and scanner platforms with low cross-site ROI standard deviation is demonstrated. Image quality review revealed diagnostic or excellent image quality in 96% of 873 patients.

Clinical feasibility of brain quantitative susceptibility mapping

PURPOSE

To evaluate the quality of brain quantitative susceptibility mapping (QSM) that is fully automatically reconstructed in clinical MRI of various neurological diseases.

METHODS

393 consecutive patients in one month were recruited for this evaluation study. QSM was reconstructed using Morphology Enabled Dipole Inversion without zero reference regularization (MEDI) and using MEDI with cerebrospinal fluid automatic zero-reference regularization to generate susceptibility values (MEDI+0). Two neuroradiologists independently assessed the image quality of MEDI+0 and MEDI and image concordance between them. Lesion susceptibility values were measured in 20 cases of glioma, 21 cases of ischemic stroke and 43 multiple sclerosis (MS) cases on both MEDI+0 and MEDI images.

RESULTS

The two neuroradiologists rated the MEDI+0 image qualities of the 393 cases as 351 (89.3%) and 362 (92.1%) excellent, 29 (7.4%) and 24 (6.1%) diagnostic, and 13 (3.3%) and 7 (1.8%) poor, and scored the concordances between MEDI+0 and MEDI as 364 (92.6%) and 351 (89.3%) excellent, 13 (3.3%) and 31 (7.9%) good, 14 (3.6%) and 9 (2.3%) intermediate, 2 (0.5%) and 2 (0.5%) poor, and 0 (0%) and 0 (0%) none. There was good correlation between MEDI+0 and MEDI in lesion susceptibility contrast of glioma, ischemic stroke, and MS cases (all p < 0.05). The MS lesion susceptibility time course from this patient cohort was found to be similar to the reported pattern: isointense initially for acute enhancing lesions, and hyperintense over the following years for active chronic lesions.

CONCLUSION

Brain QSM images of various neurological diseases have reliable diagnostic quality in clinical MRI, with MEDI+0 providing susceptibility values automatically referenced to CSF in longitudinal and cross-center studies.

Phase fMRI informs whole-brain function connectivity balance across lifespan with connection-specific aging effects during the resting state

A functional magnetic resonance imaging (fMRI) experiment produces complex-valued images consisting of pairwise magnitude and phase images. As different perspective on the same magnetic source, fMRI magnitude and phase data are complementary for brain function analysis. We collected 600-subject fMRI data during rest, decomposed via group-level independent component analysis (ICA) (mICA and pICA for magnitude and phase respectively), and calculated brain functional network connectivity matrices (mFC and pFC). The pFC matrix shows a fewer of significant connections balanced across positive and negative relationships. In comparison, the mFC matrix contains a positively-biased pattern with more significant connections. Our experiment data analyses also show that human brain maintains a whole-brain connection balance in resting state across an age span from 10 to 76 years, however, phase and magnitude data analyses reveal different connection-specific age effects on significant positive and negative subnetwork couplings.

Weak-harmonic regularization for quantitative susceptibility mapping

PURPOSE

Background-field removal is a crucial preprocessing step for quantitative susceptibility mapping (QSM). Remnants from this step often contaminate the estimated local field, which in turn leads to erroneous tissue-susceptibility reconstructions. The present work aimed to mitigate this undesirable behavior with the development of a new approach that simultaneously decouples background contributions and local susceptibility sources on QSM inversion.

METHODS

Input phase data for QSM can be seen as a composite scalar field of local effects and residual background components. We developed a new weak-harmonic regularizer to constrain the latter and to separate the 2 components. The resulting optimization problem was solved with the alternating directions of multipliers method framework to achieve fast convergence. In addition, for convenience, a new alternating directions of multipliers method-based preconditioned nonlinear projection onto dipole fields solver was developed to enable initializations with wrapped-phase distributions. Weak-harmonic QSM, with and without nonlinear projection onto dipole fields preconditioning, was compared with the original (alternating directions of multipliers method-based) total variation QSM algorithm in phantom and in vivo experiments.

RESULTS

Weak-harmonic QSM returned improved reconstructions regardless of the method used for background-field removal, although the proposed nonlinear projection onto dipole fields method often obtained better results. Streaking and shadowing artifacts were substantially suppressed, and residual background components were effectively removed.

CONCLUSION

Weak-harmonic QSM with field preconditioning is a robust dipole inversion technique and has the potential to be extended as a single-step formulation for initialization with uncombined multi-echo data.

Patents on Quantitative Susceptibility Mapping (QSM) of Tissue Magnetism

BACKGROUND

Quantitative susceptibility mapping (QSM) depicts biodistributions of tissue magnetic susceptibility sources, including endogenous iron and calcifications, as well as exogenous paramagnetic contrast agents and probes. When comparing QSM with simple susceptibility weighted MRI, QSM eliminates blooming artifacts and shows reproducible tissue susceptibility maps independent of field strength and scanner manufacturer over a broad range of image acquisition parameters. For patient care, QSM promises to inform diagnosis, guide surgery, gauge medication, and monitor drug delivery. The Bayesian framework using MRI phase data and structural prior knowledge has made QSM sufficiently robust and accurate for routine clinical practice.

OBJECTIVE

To address the lack of a summary of US patents that is valuable for QSM product development and dissemination into the MRI community.

METHOD

We searched the USPTO Full-Text and Image Database for patents relevant to QSM technology innovation. We analyzed the claims of each patent to characterize the main invented method and we investigated data on clinical utility.

RESULTS

We identified 17 QSM patents; 13 were implemented clinically, covering various aspects of QSM technology, including the Bayesian framework, background field removal, numerical optimization solver, zero filling, and zero-TE phase.

CONCLUSION

Our patent search identified patents that enable QSM technology for imaging the brain and other tissues. QSM can be applied to study a wide range of diseases including neurological diseases, liver iron disorders, tissue ischemia, and osteoporosis. MRI manufacturers can develop QSM products for more seamless integration into existing MRI scanners to improve medical care.

A robust multi-scale approach to quantitative susceptibility mapping

Quantitative Susceptibility Mapping (QSM), best known as a surrogate for tissue iron content, is becoming a highly relevant MRI contrast for monitoring cellular and vascular status in aging, addiction, traumatic brain injury and, in general, a wide range of neurological disorders. In this study we present a new Bayesian QSM algorithm, named Multi-Scale Dipole Inversion (MSDI), which builds on the nonlinear Morphology-Enabled Dipole Inversion (nMEDI) framework, incorporating three additional features: (i) improved implementation of Laplace's equation to reduce the influence of background fields through variable harmonic filtering and subsequent deconvolution, (ii) improved error control through dynamic phase-reliability compensation across spatial scales, and (iii) scalewise use of the morphological prior. More generally, this new pre-conditioned QSM formalism aims to reduce the impact of dipole-incompatible fields and measurement errors such as flow effects, poor signal-to-noise ratio or other data inconsistencies that can lead to streaking and shadowing artefacts. In terms of performance, MSDI is the first algorithm to rank in the top-10 for all metrics evaluated in the 2016 QSM Reconstruction Challenge. It also demonstrated lower variance than nMEDI and more stable behaviour in scan-rescan reproducibility experiments for different MRI acquisitions at 3 and 7 Tesla. In the present work, we also explored new forms of susceptibility MRI contrast making explicit use of the differential information across spatial scales. Specifically, we show MSDI-derived examples of: (i) enhanced anatomical detail with susceptibility inversions from short-range dipole fields (hereby referred to as High-Pass Susceptibility Mapping or HPSM), (ii) high specificity to venous-blood susceptibilities for highly regularised HPSM (making a case for MSDI-based Venography or VenoMSDI), (iii) improved tissue specificity (and possibly statistical conditioning) for Macroscopic-Vessel Suppressed Susceptibility Mapping (MVSSM), and (iv) high spatial specificity and definition for HPSM-based Susceptibility-Weighted Imaging (HPSM-SWI) and related intensity projections.

Serial quantitative neuroimaging of iron in the intracerebral hemorrhage pig model

Iron released after intracerebral hemorrhage (ICH) is damaging to the brain. Measurement of the content and distribution of iron in the hematoma could predict brain damage. In this study, 16 Yorkshire piglets were subjected to autologous blood injection ICH model and studied longitudinally using quantitative susceptibility mapping and R2* relaxivity MRI on day 1 and 7 post-ICH. Phantom calibration of susceptibility demonstrated (1) iron distribution heterogeneity within the hematoma and (2) natural absorption of iron from 154 ± 78 µg/mL (day 1) to 127 ± 33 µg/mL (day 7). R2* in the hematoma decreased at day 7. This method could be adopted for ICH in humans.