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.

Quantitative Susceptibility Mapping of Time-Dependent Susceptibility Changes in Multiple Sclerosis Lesions

BACKGROUND AND PURPOSE

MR imaging studies have demonstrated that magnetic susceptibility in multiple sclerosis lesions is dependent on lesion age. The objective of this study was to use quantitative susceptibility mapping to determine whether lesions with a hyperintense rim, indicative of iron-laden inflammatory cells (rim+), follow a unique time-dependent trajectory of susceptibility change compared with those without (rim-).

MATERIALS AND METHODS

We studied patients with MS with at least 1 new gadolinium-enhancing lesion and at least 3 longitudinal quantitative susceptibility mapping scans obtained between 1.1 and 6.1 years. Lesions were classified as rim+ if a hyperintense rim appeared on quantitative susceptibility mapping at any time. A multilevel growth curve model compared longitudinal susceptibility among rim+ and rim- lesions.

RESULTS

Thirty-two new gadolinium-enhancing lesions from 19 patients with MS were included, and 16 lesions (50%) were identified as rim+. Quantitative susceptibility mapping rim+ lesions were larger than rim- lesions with gadolinium enhancement (P < .001). Among all lesions, susceptibility increased sharply after enhancement to a peak between 1 and 2 years followed by a decrease. The overall susceptibility curve height for rim- lesions was 4.27 parts per billion lower than that for rim+ lesions (P = .01). Rim- lesions demonstrated a higher linear slope relative to rim+ lesions (P = .023) but faster cubic decay relative to rim+ lesions (P = .005). Rim- lesions started decaying approximately 2 years earlier compared with rim+ lesions.

CONCLUSIONS

There was a marked difference in the susceptibility temporal trajectory between rim+ and rim- lesions during the first 6 years of lesion formation. Most rim+ lesions retain iron for years after the initial lesion appearance.

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.

Rapid automated liver quantitative susceptibility mapping

BACKGROUND

Accurate measurement of the liver iron concentration (LIC) is needed to guide iron-chelating therapy for patients with transfusional iron overload. In this work, we investigate the feasibility of automated quantitative susceptibility mapping (QSM) to measure the LIC.

PURPOSE

To develop a rapid, robust, and automated liver QSM for clinical practice.

STUDY TYPE

Prospective.

POPULATION

13 healthy subjects and 22 patients.

FIELD STRENGTH/SEQUENCES

1.5 T and 3 T/3D multiecho gradient-recalled echo (GRE) sequence.

ASSESSMENT

Data were acquired using a 3D GRE sequence with an out-of-phase echo spacing with respect to each other. All odd echoes that were in-phase (IP) were used to initialize the fat-water separation and field estimation (T₂ *-IDEAL) before performing QSM. Liver QSM was generated through an automated pipeline without manual intervention. This IP echo-based initialization method was compared with an existing graph cuts initialization method (simultaneous phase unwrapping and removal of chemical shift, SPURS) in healthy subjects (n = 5). Reproducibility was assessed over four scanners at two field strengths from two manufacturers using healthy subjects (n = 8). Clinical feasibility was evaluated in patients (n = 22).

STATISTICAL TESTS

IP and SPURS initialization methods in both healthy subjects and patients were compared using paired t-test and linear regression analysis to assess processing time and region of interest (ROI) measurements. Reproducibility of QSM, R₂ *, and proton density fat fraction (PDFF) among the four different scanners was assessed using linear regression, Bland-Altman analysis, and the intraclass correlation coefficient (ICC).

RESULTS

Liver QSM using the IP method was found to be ~5.5 times faster than SPURS (P < 0.05) in initializing T₂ *-IDEAL with similar outputs. Liver QSM using the IP method were reproducibly generated in all four scanners (average coefficient of determination 0.95, average slope 0.90, average bias 0.002 ppm, 95% limits of agreement between -0.06 to 0.07 ppm, ICC 0.97).

DATA CONCLUSION

Use of IP echo-based initialization enables robust water/fat separation and field estimation for automated, rapid, and reproducible liver QSM for clinical applications.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:725-732.

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.

Discontinuity Preserving Liver MR Registration with 3D Active Contour Motion Segmentation

OBJECTIVE

The sliding motion of the liver during respiration violates the homogeneous motion smoothness assumption in conventional non-rigid image registration and commonly results in compromised registration accuracy. This paper presents a novel approach, registration with 3D active contour motion segmentation (RAMS), to improve registration accuracy with discontinuity-aware motion regularization.

METHODS

A Markov random field-based discrete optimization with dense displacement sampling and self-similarity context metric is used for registration, while a graph cuts-based 3D active contour approach is applied to segment the sliding interface. In the first registration pass, a mask-free L1 regularization on an image-derived minimum spanning tree is performed to allow motion discontinuity. Based on the motion field estimates, a coarse segmentation finds the motion boundaries. Next, based on MR signal intensity, a fine segmentation aligns the motion boundaries with anatomical boundaries. In the second registration pass, smoothness constraints across the segmented sliding interface are removed by masked regularization on a minimum spanning forest and masked interpolation of the motion field.

RESULTS

For in vivo breath-hold abdominal MRI data, the motion masks calculated by RAMS are highly consistent with manual segmentations in terms of Dice similarity and bidirectional local distance measure. These automatically obtained masks are shown to substantially improve registration accuracy for both the proposed discrete registration as well as conventional continuous non-rigid algorithms.

CONCLUSION/SIGNIFICANCE

The presented results demonstrated the feasibility of automated segmentation of the respiratory sliding motion interface in liver MR images and the effectiveness of using the derived motion masks to preserve motion discontinuity.

Multicenter reproducibility of quantitative susceptibility mapping in a gadolinium phantom using MEDI+0 automatic zero referencing

PURPOSE

To determine the reproducibility of quantitative susceptibility mapping at multiple sites on clinical and preclinical scanners (1.5 T, 3 T, 7 T, and 9.4 T) from different vendors (Siemens, GE, Philips, and Bruker) for standardization of multicenter studies.

METHODS

Seven phantoms distributed from the core site, each containing 5 compartments with gadolinium solutions with fixed concentrations between 0.625 mM and 10 mM. Multi-echo gradient echo scans were performed at 1.5 T, 3 T, 7 T, and 9.4 T on 12 clinical and 3 preclinical scanners. DICOM images from the scans were processed into quantitative susceptibility maps using the Laplacian boundary value (LBV) and MEDI+0 automatic uniform reference algorithm. Region of interest (ROI) analyses were performed by a physicist to determine agreement between results from all sites. Measurement reproducibility was assessed using regression, Bland-Altman plots, and the intra-class correlation coefficient (ICC).

RESULTS

Quantitative susceptibility mapping (QSM) from all scanners had similar, artifact-free visual appearance. Regression analysis showed a linear relationship between gadolinium concentrations and average QSM measurements for all phantoms (y = 350x - 0.0346, r² >0.99). The SD of measurements increased almost linearly from 32 ppb to 230 ppb as the measured susceptibility increased from 0.26 ppm to 3.56 ppm. A Bland-Altman plot showed the bias, upper, and lower limits of agreement for all comparisons were -10, -210, and 200 ppb, respectively. The ICC was 0.991 with a 95% CI (0.973, 0.99).

CONCLUSIONS

QSM shows excellent multicenter reproducibility for a large range of susceptibility values encountered in cranial and extra-cranial applications on a diverse set of scanner platforms.

Clinical Integration of Quantitative Susceptibility Mapping Magnetic Resonance Imaging into Neurosurgical Practice

OBJECTIVE

To introduce quantitative susceptibility mapping (QSM), a novel magnetic resonance imaging sequence, to the field of neurosurgery.

METHODS

QSM is introduced both in its historical context and by providing a brief overview of the physics behind the technique tailored to a neurosurgical audience. Its application to clinical neurosurgery is then highlighted using case examples.

RESULTS

QSM offers a quantitative assessment of susceptibility (previously considered as an artifact) via a single, straightforward gradient echo acquisition. QSM differs from standard susceptibility weighted imaging in its ability to both quantify and precisely localize susceptibility effects. Clinical applications of QSM are wide reaching and include precise localization of the deep nuclei for deep brain stimulation electrode placement, differentiation between blood products and calcification within brain lesions, and enhanced sensitivity of cerebral micrometastasis identification.

CONCLUSIONS

We present this diverse range of QSM's clinical applications to neurosurgical care via case examples. QSM can be obtained in all patients able to undergo magnetic resonance imaging and is easily integratable into busy neuroradiology programs because of its short acquisition time and straightforward, automated offline postprocessing workflow. Clinical integration of QSM may help clinicians better identify and characterize neurosurgical lesions, thereby improving patient care.

Brain Iron Distribution after Multiple Doses of Ultra-small Superparamagnetic Iron Oxide Particles in Rats

The purpose of this study is to determine the effects of high cumulative doses of ultra-small paramagnetic iron oxide (USPIO) used in neuroimaging studies. We intravenously administered 8 mg/kg of 2 USPIO compounds daily for 4 wk to male Sprague-Dawley rats (Crl:SD). Multiecho gradient-echo MRI, serum iron levels, and histology were performed at the end of dosing and after a 7-d washout period. R2* maps and quantitative susceptibility maps (QSM) were generated from multiecho gradient-echo data. R2* maps and QSM showed iron accumulation in brain ventricles on MR images acquired at the 4- and 5-wk time points. Estimates from QSM data showed ventricular iron concentration was equal to or higher than serum iron concentration. Histologic analysis revealed choroid plexus hemosiderosis and midbrain vacuolation, without iron deposition in brain parenchyma. Serum iron levels increased with administration of both compounds, and a 7-d washout period effectively reduced serum iron levels of one but not both of the compounds. High cumulative doses from multiple, frequent administrations of USPIO can lead to iron deposition in brain ventricles, resulting in persistent signal loss on T2*-weighted images. Techniques such as QSM are helpful in quantifying iron biodistribution in this situation.

Brain Injury Lesion Imaging Using Preconditioned Quantitative Susceptibility Mapping without Skull Stripping

BACKGROUND AND PURPOSE

Identifying cerebral microhemorrhage burden can aid in the diagnosis and management of traumatic brain injury, stroke, hypertension, and cerebral amyloid angiopathy. MR imaging susceptibility-based methods are more sensitive than CT for detecting cerebral microhemorrhage, but methods other than quantitative susceptibility mapping provide results that vary with field strength and TE, require additional phase maps to distinguish blood from calcification, and depict cerebral microhemorrhages as bloom artifacts. Quantitative susceptibility mapping provides universal quantification of tissue magnetic property without these constraints but traditionally requires a mask generated by skull-stripping, which can pose challenges at tissue interphases. We evaluated the preconditioned quantitative susceptibility mapping MR imaging method, which does not require skull-stripping, for improved depiction of brain parenchyma and pathology.

MATERIALS AND METHODS

Fifty-six subjects underwent brain MR imaging with a 3D multiecho gradient recalled echo acquisition. Mask-based quantitative susceptibility mapping images were created using a commonly used mask-based quantitative susceptibility mapping method, and preconditioned quantitative susceptibility images were made using precondition-based total field inversion. All images were reviewed by a neuroradiologist and a radiology resident.

RESULTS

Ten subjects (18%), all with traumatic brain injury, demonstrated blood products on 3D gradient recalled echo imaging. All lesions were visible on preconditioned quantitative susceptibility mapping, while 6 were not visible on mask-based quantitative susceptibility mapping. Thirty-one subjects (55%) demonstrated brain parenchyma and/or lesions that were visible on preconditioned quantitative susceptibility mapping but not on mask-based quantitative susceptibility mapping. Six subjects (11%) demonstrated pons artifacts on preconditioned quantitative susceptibility mapping and mask-based quantitative susceptibility mapping; they were worse on preconditioned quantitative susceptibility mapping.

CONCLUSIONS

Preconditioned quantitative susceptibility mapping MR imaging can bring the benefits of quantitative susceptibility mapping imaging to clinical practice without the limitations of mask-based quantitative susceptibility mapping, especially for evaluating cerebral microhemorrhage-associated pathologies, such as traumatic brain injury.

Magnetic susceptibility increases as diamagnetic molecules breakdown: Myelin digestion during multiple sclerosis lesion formation contributes to increase on QSM

BACKGROUND

The pathological processes in the first weeks of multiple sclerosis (MS) lesion formation include myelin digestion that breaks chemical bonds in myelin lipid layers. This can increase lesion magnetic susceptibility, which is a potentially useful biomarker in MS patient management, but not yet investigated.

PURPOSE

To understand and quantify the effects of myelin digestion on quantitative susceptibility mapping (QSM) of MS lesions.

STUDY TYPE

Histological and QSM analyses on in vitro models of myelin breakdown and MS lesion formation in vivo.

POPULATION/SPECIMENS

Acutely demyelinating white matter lesions from MS autopsy tissue were stained with the lipid dye oil red O. Myelin basic protein (MBP), a major membrane protein of myelin, was digested with trypsin. Purified human myelin was denatured with sodium dodecyl sulfate (SDS). QSM was performed on phantoms containing digestion products and untreated controls. In vivo QSM was performed on five MS patients with newly enhancing lesions, and then repeated within 2 weeks.

FIELD STRENGTH/SEQUENCE

3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> -weighted spoiled multiecho gradient echo scans performed at 3T.

ASSESSMENT

Region of interest analyses were performed by a biochemist and a neuroradiologist to determine susceptibility changes on in vitro and in vivo QSM images.

STATISTICAL TESTS

Not applicable.

RESULTS

MBP degradation by trypsin increased the QSM measurement by an average of 112 ± 37 ppb, in excellent agreement with a theoretical estimate of 111 ppb. Degradation of human myelin by SDS increased the QSM measurement by 23 ppb. As MS lesions changed from gadolinium enhancing to nonenhancing over an average of 15.8 ± 3.7 days, their susceptibility increased by an average of 7.5 ± 6.3 ppb.

DATA CONCLUSION

Myelin digestion in the early stages of MS lesion formation contributes to an increase in tissue susceptibility, detectable by QSM, as a lesion evolves from gadolinium enhancing to nonenhancing.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2018;47:1281-1287.

Cerebral metabolic rate of oxygen (CMRO2 ) mapping by combining quantitative susceptibility mapping (QSM) and quantitative blood oxygenation level-dependent imaging (qBOLD)

PURPOSE

To map the cerebral metabolic rate of oxygen (CMRO₂ ) by estimating the oxygen extraction fraction (OEF) from gradient echo imaging (GRE) using phase and magnitude of the GRE data.

THEORY AND METHODS

3D multi-echo gradient echo imaging and perfusion imaging with arterial spin labeling were performed in 11 healthy subjects. CMRO₂ and OEF maps were reconstructed by joint quantitative susceptibility mapping (QSM) to process GRE phases and quantitative blood oxygen level-dependent (qBOLD) modeling to process GRE magnitudes. Comparisons with QSM and qBOLD alone were performed using ROI analysis, paired t-tests, and Bland-Altman plot.

RESULTS

The average CMRO₂ value in cortical gray matter across subjects were 140.4 ± 14.9, 134.1 ± 12.5, and 184.6 ± 17.9 μmol/100 g/min, with corresponding OEFs of 30.9 ± 3.4%, 30.0 ± 1.8%, and 40.9 ± 2.4% for methods based on QSM, qBOLD, and QSM+qBOLD, respectively. QSM+qBOLD provided the highest CMRO₂ contrast between gray and white matter, more uniform OEF than QSM, and less noisy OEF than qBOLD.

CONCLUSION

Quantitative CMRO₂ mapping that fits the entire complex GRE data is feasible by combining QSM analysis of phase and qBOLD analysis of magnitude.

MEDI+0: Morphology enabled dipole inversion with automatic uniform cerebrospinal fluid zero reference for quantitative susceptibility mapping

PURPOSE

To develop a quantitative susceptibility mapping (QSM) method with a consistent zero reference using minimal variation in cerebrospinal fluid (CSF) susceptibility.

THEORY AND METHODS

The ventricular CSF was automatically segmented on the R2* map. An L₂ -regularization was used to enforce CSF susceptibility homogeneity within the segmented region, with the averaged CSF susceptibility as the zero reference. This regularization for CSF homogeneity was added to the model used in a prior QSM method (morphology enabled dipole inversion [MEDI]). Therefore, the proposed method was referred to as MEDI+0 and compared with MEDI in a numerical simulation, in multiple sclerosis (MS) lesions, and in a reproducibility study in healthy subjects.

RESULTS

In both the numerical simulations and in vivo experiments, MEDI+0 not only decreased the susceptibility variation within the ventricular CSF, but also suppressed the artifact near the lateral ventricles. In the simulation, MEDI+0 also provided more accurate quantification compared to MEDI in the globus pallidus, substantia nigra, corpus callosum, and internal capsule. MEDI+0 measurements of MS lesion susceptibility were in good agreement with those obtained by MEDI. Finally, both MEDI+0 and MEDI showed good and similar intrasubject reproducibility.

CONCLUSION

QSM with a minimal variation in ventricular CSF is viable to provide a consistent zero reference while improving image quality. Magn Reson Med 79:2795-2803, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Vastly accelerated linear least-squares fitting with numerical optimization for dual-input delay-compensated quantitative liver perfusion mapping

PURPOSE

To propose an efficient algorithm to perform dual input compartment modeling for generating perfusion maps in the liver.

METHODS

We implemented whole field-of-view linear least squares (LLS) to fit a delay-compensated dual-input single-compartment model to very high temporal resolution (four frames per second) contrast-enhanced 3D liver data, to calculate kinetic parameter maps. Using simulated data and experimental data in healthy subjects and patients, whole-field LLS was compared with the conventional voxel-wise nonlinear least-squares (NLLS) approach in terms of accuracy, performance, and computation time.

RESULTS

Simulations showed good agreement between LLS and NLLS for a range of kinetic parameters. The whole-field LLS method allowed generating liver perfusion maps approximately 160-fold faster than voxel-wise NLLS, while obtaining similar perfusion parameters.

CONCLUSIONS

Delay-compensated dual-input liver perfusion analysis using whole-field LLS allows generating perfusion maps with a considerable speedup compared with conventional voxel-wise NLLS fitting. Magn Reson Med 79:2415-2421, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Cardiac quantitative susceptibility mapping (QSM) for heart chamber oxygenation

PURPOSE

To demonstrate the feasibility of in vivo quantitative susceptibility mapping (QSM) in cardiac MRI and to show that mixed-venous oxygen saturation (SvO₂ ) can be measured non-invasively using QSM.

METHODS

Electrocardiographic-gated multi-echo 2D gradient echo data were collected at 1.5 T from 14 healthy volunteers during successive breath-holds. Phase wraps and fat chemical shift were removed using a graph-cut-based phase analysis and IDEAL in an iterative approach. The large susceptibility range from air in the lungs to blood in the heart was addressed by using the preconditioning approach in the dipole field inversion. SvO₂ was calculated based on the difference in blood susceptibility between the right ventricle (RV) and left ventricle (LV). Cardiac QSM quality was assessed by two independent readers.

RESULTS

Nine out of fourteen volunteers (64%) yielded interpretable cardiac QSM. QSM maps showed strong differential contrast between RV and LV blood with RV blood having higher susceptibility values (291.5 ± 32.4 ppb), which correspond to 78.3 ± 2.3% SvO₂ .

CONCLUSION

In vivo cardiac QSM is feasible and can be used to measure SvO₂ , but improvements in data acquisition are needed. Magn Reson Med 79:1545-1552, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Coherence enhancement in quantitative susceptibility mapping by means of anisotropic weighting in morphology enabled dipole inversion

PURPOSE

To investigate an anisotropic structural prior in morphology enabled dipole inversion (MEDI) for improving accuracy in quantitative susceptibility mapping (QSM).

THEORY AND METHODS

Anisotropic weighting (AW) was devised and implemented to incorporate orientation information into the edge agreement in the MEDI method. AW performance was compared with isotropic weighting by testing and validating on in vivo brain multiple orientation MRI data using COSMOS and the (33) component of the susceptibility tensor as reference.

RESULTS

Suppressing streaking artifacts, AW improved not only QSM image quality but also accuracy in terms of RMSE (root mean square error), HFEN (high frequency error norm), SSIM (structural similarity index), and GDA (gradient direction agreement). In addition, it outperformed isotropic weighting in region of interest-based analysis. From a computational perspective, AW was as fast as isotropic weighting, taking approximately the same central processing unit times.

CONCLUSION

Using AW in MEDI improves QSM accuracy compared with isotropic weighting. Magn Reson Med 79:1172-1180, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Clinical quantitative susceptibility mapping (QSM): Biometal imaging and its emerging roles in patient care

Quantitative susceptibility mapping (QSM) has enabled magnetic resonance imaging (MRI) of tissue magnetic susceptibility to advance from simple qualitative detection of hypointense blooming artifacts to precise quantitative measurement of spatial biodistributions. QSM technology may be regarded to be sufficiently developed and validated to warrant wide dissemination for clinical applications of imaging isotropic susceptibility, which is dominated by metals in tissue, including iron and calcium. These biometals are highly regulated as vital participants in normal cellular biochemistry, and their dysregulations are manifested in a variety of pathologic processes. Therefore, QSM can be used to assess important tissue functions and disease. To facilitate QSM clinical translation, this review aims to organize pertinent information for implementing a robust automated QSM technique in routine MRI practice and to summarize available knowledge on diseases for which QSM can be used to improve patient care. In brief, QSM can be generated with postprocessing whenever gradient echo MRI is performed. QSM can be useful for diseases that involve neurodegeneration, inflammation, hemorrhage, abnormal oxygen consumption, substantial alterations in highly paramagnetic cellular iron, bone mineralization, or pathologic calcification; and for all disorders in which MRI diagnosis or surveillance requires contrast agent injection. Clinicians may consider integrating QSM into their routine imaging practices by including gradient echo sequences in all relevant MRI protocols.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2017;46:951-971.

Quantitative susceptibility mapping-based cerebral metabolic rate of oxygen mapping with minimum local variance

PURPOSE

The objective of this study was to demonstrate the feasibility of a cerebral metabolic rate of oxygen (CMRO₂ ) mapping method based on its minimum local variance (MLV) without vascular challenge using quantitative susceptibility mapping (QSM) and cerebral blood flow (CBF).

METHODS

Three-dimensional multi-echo gradient echo imaging and arterial spin labeling were performed in 11 healthy subjects to calculate QSM and CBF. Minimum local variance was used to compute whole-brain CMRO₂ map from QSM and CBF. The MLV method was compared with a reference method using the caffeine challenge. Their agreement within the cortical gray matter (CGM) was assessed on CMRO₂ and oxygen extraction fraction (OEF) maps at both baseline and challenge states.

RESULTS

Mean CMRO₂ (in µmol/100 g/min) obtained in CGM using the caffeine challenge and MLV were 142 ± 16.5 and 139 ± 14.8 µmol/100 g/min, respectively; the corresponding baseline OEF were 33.0 ± 4.0% and 31.8 ± 3.2%, respectively. The MLV and caffeine challenge methods showed no statistically significant differences across subjects with small ( < 4%) biases in CMRO₂ and OEF values.

CONCLUSIONS

Minimum local variance-based CMRO₂ mapping without vascular challenge using QSM and arterial spin labeling is feasible in healthy subjects. Magn Reson Med 79:172-179, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Bone quantitative susceptibility mapping using a chemical species-specific R2* signal model with ultrashort and conventional echo data

PURPOSE

To develop quantitative susceptibility mapping (QSM) of bone using an ultrashort echo time (UTE) gradient echo (GRE) sequence for signal acquisition and a bone-specific effective transverse relaxation rate ( R2*) to model water-fat MR signals for field mapping.

METHODS

Three-dimensional radial UTE data (echo times ≥ 40 μs) was acquired on a 3 Tesla scanner and fitted with a bone-specific signal model to map the chemical species and susceptibility field. Experiments were performed ex vivo on a porcine hoof and in vivo on healthy human subjects (n = 7). For water-fat separation, a bone-specific model assigning R2* decay mostly to water was compared with the standard models that assigned the same decay for both fat and water. In the ex vivo experiment, bone QSM was correlated with CT.

RESULTS

Compared with standard models, the bone-specific R2* method significantly reduced errors in the fat fraction within the cortical bone in all tested data sets, leading to reduced artifacts in QSM. Good correlation was found between bone CT and QSM values in the porcine hoof (R²  = 0.77). Bone QSM was successfully generated in all subjects.

CONCLUSIONS

The QSM of bone is feasible using UTE with a conventional echo time GRE acquisition and a bone-specific R2* signal model. Magn Reson Med 79:121-128, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Primal-dual and forward gradient implementation for quantitative susceptibility mapping

PURPOSE

To investigate the computational aspects of the prior term in quantitative susceptibility mapping (QSM) by (i) comparing the Gauss-Newton conjugate gradient (GNCG) algorithm that uses numerical conditioning (ie, modifies the prior term) with a primal-dual (PD) formulation that avoids this, and (ii) carrying out a comparison between a central and forward difference scheme for the discretization of the prior term.

THEORY AND METHODS

A spatially continuous formulation of the regularized QSM inversion problem and its PD formulation were derived. The Chambolle-Pock algorithm for PD was implemented and its convergence behavior was compared with that of GNCG for the original QSM. Forward and central difference schemes were compared in terms of the presence of checkerboard artifacts. All methods were tested and validated on a gadolinium phantom, ex vivo brain blocks, and in vivo brain MRI data with respect to COSMOS.

RESULTS

The PD approach provided a faster convergence rate than GNCG. The GNCG convergence rate slowed considerably with smaller (more accurate) values of the conditioning parameter. Using a forward difference suppressed the checkerboard artifacts in QSM, as compared with the central difference. The accuracy of PD and GNCG were validated based on excellent correlation with COSMOS.

CONCLUSIONS

The PD approach with forward difference for the gradient showed improved convergence and accuracy over the GNCG method using central difference. Magn Reson Med 78:2416-2427, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Susceptibility underestimation in a high-susceptibility phantom: Dependence on imaging resolution, magnitude contrast, and other parameters

PURPOSE

We assessed the accuracy of quantitative susceptibility mapping in a gadolinium balloon phantom with a large range of susceptibility values and imaging resolutions at 1.5 and 3 Tesla (T).

THEORY AND METHODS

The phantom contained sources with susceptibility values of 0.4, 0.8, 1.6, and 3.2 ppm and was imaged at isotropic resolutions of 0.7, 0.8, 1.2, and 1.8 mm. Numerical simulations were performed to match the experimental findings. Voxel sensitivity effects were used to explain the susceptibility underestimations.

RESULTS

Both phantom data and simulation demonstrated that systematic underestimation of the susceptibility values increased with voxel size, field strength, and object susceptibility.

CONCLUSION

The underestimation originates from the signal formation in a voxel, which can be described by the voxel sensitivity function. The amount of underestimation is thus affected by imaging resolution, magnitude contrast, image filtering, and details of the susceptibility inclusions such as the susceptibility value and geometry. High-resolution imaging is therefore needed for accurate reconstruction of QSM values, especially at higher susceptibilities. Magn Reson Med 78:1080-1086, 2017. © 2016 International Society for Magnetic Resonance in Medicine.