Multi-echo dipole inversion for magnetic susceptibility mapping

PURPOSE

Reconstructing tissue magnetic susceptibility (QSM) from MRI phase data involves solving multiple consecutive ill-posed inverse problems such as phase unwrapping, background field removal, and field-to-source inversion. Multi-echo acquisitions present an additional challenge, as the magnetization field is typically computed from the multiple phase data prior to reconstructing the susceptibility map. Processing the multiple phase data introduces errors during the field estimation, violating assumptions of the subsequent inverse problems, manifesting as streaking artifacts in the susceptibility map. To address this challenge, we propose a multi-echo field-to-source forward model that forgoes the field estimation step. Moreover, we propose a fully general underestimation correction step to recover susceptibility sources that were regularized away during the field-to-source inversion.

METHODS

The multi-echo forward model and correction step were validated on the QSM Challenge 2.0 datasets and compared to the standard single field-to-source model in in vivo human brains using different types of deconvolution algorithms.

RESULTS

On the QSM Challenge 2.0 datasets the multi-echo forward model and correction step attain state-of-the-art results on all metrics by a wide margin. Experiments in in vivo brains show that the multi-echo model is in agreement with the single field-to-source model and that the proposed forward model and correction step can be used with any available dipole inversion method.

CONCLUSION

A multi-echo field-to-source forward model forgoes the need to fit multi-echo phase data and achieves state-of-the-art results on the QSM Challenge 2.0 data. Underestimated low-frequency susceptibility distributions can be partially recovered using a correction step.

An optimized and highly repeatable MRI acquisition and processing pipeline for quantitative susceptibility mapping in the head-and-neck region

PURPOSE

Quantitative Susceptibility Mapping (QSM) is an emerging technique sensitive to disease-related changes including oxygenation. It is extensively used in brain studies and has increasing clinical applications outside the brain. Here we present the first MRI acquisition protocol and QSM pipeline optimized for the head-and-neck region together with a repeatability analysis performed in healthy volunteers.

METHODS

We investigated both the intrasession and the intersession repeatability of the optimized method in 10 subjects. We also implemented two, Tikhonov-regularisation-based susceptibility calculation techniques that were found to have higher contrast-to-noise than existing methods in the head-and-neck region. Repeatability was evaluated by calculating the distributions of susceptibility differences between repeated scans and the corresponding minimum detectable effect sizes (MDEs).

RESULTS

Deep brain regions had higher QSM repeatability than neck regions. As expected, intrasession repeatability was generally better than intersession repeatability. Susceptibility maps calculated using projection onto dipole fields for background field removal were more repeatable than using the Laplacian boundary value method in the head-and-neck region. Small (short-axis diameter <5 mm) lymph nodes had the lowest repeatability (MDE = 0.27 ppm) as imperfect segmentation included some of the surrounding paramagnetic fatty fascia, highlighting the importance of accurate region delineation. MDEs in the larger lymph nodes (0.16 ppm), submandibular glands (0.10 ppm), and especially the parotid glands (0.06 ppm) were much lower, comparable to those of the brain regions.

CONCLUSIONS

The high repeatability of the acquisition and pipeline optimized for QSM will facilitate clinical studies in the head-and-neck region.

Longitudinal Magnetic Resonance Imaging of Cerebral Microbleeds in Multiple Sclerosis Patients

We hypothesized that cerebral microbleeds (CMBs) in multiple sclerosis (MS) patients will be detected with higher prevalence compared to healthy controls (HC) and that quantitative susceptibility mapping (QSM) will help remove false positives seen in susceptibility weighted imaging (SWI). A cohort of 100 relapsing remitting MS subjects scanned at 3T were used to validate a set of CMB detection guidelines specifically using QSM. A second longitudinal cohort of 112 MS and 25 HCs, also acquired at 3T, was reviewed across two time points. Both cohorts were imaged with SWI and fluid attenuated inversion recovery. Fourteen subjects in the first cohort (14%, 95% CI 8-21%) and twenty-one subjects in the second cohort (18.7%, 95% CI 11-27%) had at least one CMB. The combined information from SWI and QSM allowed us to discern stable CMBs and new CMBs from potential mimics and evaluate changes over time. The longitudinal results demonstrated that longer disease duration increased the chance to develop new CMBs. Higher age was also associated with increased CMB prevalence for MS and HC. We observed that MS subjects developed new CMBs between time points, indicating the need for longitudinal quantitative imaging of CMBs.

Quantitative susceptibility mapping (QSM) with an extended physical model for MRI frequency contrast in the brain: a proof-of-concept of quantitative susceptibility and residual (QUASAR) mapping

Quantitative susceptibility mapping (QSM) aims to calculate the tissue's magnetic susceptibility distribution from its perturbing effect on the MRI static main magnetic field. The method is increasingly being applied to study iron and myelin in clinical and preclinical settings. However, recent experimental and theoretical findings have challenged the fundamental theoretical assumptions that form the basis of current numerical implementations of QSM algorithms. The present work introduces a new class of susceptibility mapping algorithms, termed quantitative susceptibility and residual mapping (QUASAR), which takes into account frequency contributions not related to the spatial variation of bulk magnetic susceptibility in the Lorentz sphere model. We present a simple proof-of-concept QUASAR algorithm that, unlike most of the QSM algorithms currently used widely, results in an improved anatomical accuracy of the susceptibility distribution without any a priori assumptions about the susceptibility distribution during the field-to-source inversion. The algorithm was evaluated both in silico and in vivo in the preclinical setting. Our preliminary application of QUASAR in rodents provides the first in vivo evidence that the susceptibility-field model traditionally used in the QSM field cannot fully explain the frequency contrast in brain tissues. Only when an additional local frequency contribution is added to the physical model can the frequency contrast in the brain be related properly to the underlying anatomy.

On the sensitivity of quantitative susceptibility mapping for measuring trabecular bone density

Purpose

To develop a methodological framework to simultaneously measure R2* and magnetic susceptibility in trabecularized yellow bone marrow and to investigate the sensitivity of Quantitative Susceptibility Mapping (QSM) for measuring trabecular bone density using a non‐UTE multi‐gradient echo sequence.

Methods

The ankle of 16 healthy volunteers and two patients was scanned using a time‐interleaved multi‐gradient‐echo (TIMGRE) sequence. After field mapping based on water–fat separation methods and background field removal based on the Laplacian boundary value method, three different QSM dipole inversion schemes were implemented. Mean susceptibility values in regions of different trabecular bone density in the calcaneus were compared to the corresponding values in the R2* maps, bone volume to total volume ratios (BV/TV) estimated from high resolution imaging (in 14 subjects), and CT attenuation (in two subjects). In addition, numerical simulations were performed in a simplified trabecular bone model of randomly positioned spherical bone inclusions to verify and compare the scaling of R2* and susceptibility with BV/TV.

Results

Differences in calcaneus trabecularization were well depicted in susceptibility maps, in good agreement with high‐resolution MR and CT images. Simulations and in vivo scans showed a linear relationship of measured susceptibility with BV/TV and R2*. The ankle in vivo results showed a strong linear correlation between susceptibility and R2* (R² = 0.88, p < 0.001) with a slope and intercept of −0.004 and 0.2 ppm, respectively.

Conclusions

A method for multi‐paramteric mapping, including R2*‐mapping and QSM was developed for measuring trabecularized yellow bone marrow, showing good sensitivity of QSM for measuring trabecular bone density.

An illustrated comparison of processing methods for phase MRI and QSM: removal of background field contributions from sources outside the region of interest

The elimination of so-called background fields is an essential step in phase MRI and quantitative susceptibility mapping (QSM). Background fields, which are caused by sources outside the region of interest (ROI), are often one to two orders of magnitude stronger than tissue-related field variations from within the ROI, hampering quantitative interpretation of field maps. This paper reviews the current literature on background elimination algorithms for QSM and provides insights into similarities and differences between the many algorithms proposed. We discuss the basic theoretical foundations and derive fundamental limitations of background field elimination. Copyright © 2016 John Wiley & Sons, Ltd.

Metallic implant geometry and susceptibility estimation using multispectral B0 field maps

PURPOSE

To estimate the susceptibility and the geometry of metallic implants from multispectral imaging (MSI) information, to separate the metal implant region from the surrounding signal loss region.

THEORY AND METHODS

The susceptibility map of signal-void regions is estimated from MSI B₀ field maps using total variation (TV) regularized inversion. Voxels with susceptibility estimates above a predetermined threshold are identified as metal. The accuracy of the estimated susceptibility and implant geometry was evaluated in simulations, phantom, and in vivo experiments.

RESULTS

The proposed method provided more accurate susceptibility estimation compared with a previous method without TV regularization, in both simulations and phantom experiments. In the phantom experiment where the actual implant was 40% of the signal-void region, the mean estimated susceptibility was close to the susceptibility in literature, and the precision and recall of the estimated geometry was 85% and 93%. In vivo studies in subjects with hip implants also demonstrated that the proposed method can distinguish implants from surrounding low-signal tissues, such as cortical bone.

CONCLUSION

The proposed method can improve the delineation of metallic implant geometry by distinguishing metal voxels from artificial signal voids and low-signal tissues by estimating the susceptibility maps. Magn Reson Med 77:2402-2413, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Hemisphere, gender and age-related effects on iron deposition in deep gray matter revealed by quantitative susceptibility mapping

The purpose of this work was to investigate the effects of hemispheric location, gender and age on susceptibility value, as well as the association between susceptibility value and diffusional metrics, in deep gray matter. Iron content was estimated in vivo using quantitative susceptibility mapping. Microstructure was probed using diffusional kurtosis imaging. Regional susceptibility and diffusional metrics were measured for the putamen, caudate nucleus, globus pallidus, thalamus, substantia nigra and red nucleus in 42 healthy adults (age range 25-78 years). Susceptibility value was significantly higher in the left than the right side of the caudate nucleus (P = 0.043) and substantia nigra (P < 0.001). Women exhibited lower susceptibility values than men in the thalamus (P < 0.001) and red nucleus (P = 0.032). Significant age-related increases of susceptibility were observed in the putamen (P < 0.001), red nucleus (P < 0.001), substantia nigra (P = 0.004), caudate nucleus (P < 0.001) and globus pallidus (P = 0.017). The putamen exhibited the highest rate of iron accumulation with aging (slope of linear regression = 0.73 × 10(-3) ppm/year), which was nearly twice those in substantia nigra (slope = 0.40 × 10(-3) ppm/year) and caudate nucleus (slope = 0.39 × 10(-3) ppm/year). Significant positive correlations between the susceptibility value and diffusion measurements were observed for fractional anisotropy (P = 0.045) and mean kurtosis (P = 0.048) in the putamen without controlling for age. Neither correlation was significant after controlling for age. Hemisphere, gender and age-related differences in iron measurements were observed in deep gray matter. Notably, the putamen exhibited the highest rate of increase in susceptibility with aging. Correlations between susceptibility value and microstructural measurements were inconclusive. These findings could provide new clues for unveiling mechanisms underlying iron-related neurodegenerative diseases.

Fast image reconstruction with L2-regularization

PURPOSE

We introduce L2-regularized reconstruction algorithms with closed-form solutions that achieve dramatic computational speed-up relative to state of the art L1- and L2-based iterative algorithms while maintaining similar image quality for various applications in MRI reconstruction.

MATERIALS AND METHODS

We compare fast L2-based methods to state of the art algorithms employing iterative L1- and L2-regularization in numerical phantom and in vivo data in three applications; (i) Fast Quantitative Susceptibility Mapping (QSM), (ii) Lipid artifact suppression in Magnetic Resonance Spectroscopic Imaging (MRSI), and (iii) Diffusion Spectrum Imaging (DSI). In all cases, proposed L2-based methods are compared with the state of the art algorithms, and two to three orders of magnitude speed up is demonstrated with similar reconstruction quality.

RESULTS

The closed-form solution developed for regularized QSM allows processing of a three-dimensional volume under 5 s, the proposed lipid suppression algorithm takes under 1 s to reconstruct single-slice MRSI data, while the PCA based DSI algorithm estimates diffusion propagators from undersampled q-space for a single slice under 30 s, all running in Matlab using a standard workstation.

CONCLUSION

For the applications considered herein, closed-form L2-regularization can be a faster alternative to its iterative counterpart or L1-based iterative algorithms, without compromising image quality.

Increased iron accumulation occurs in the earliest stages of demyelinating disease: an ultra-high field susceptibility mapping study in Clinically Isolated Syndrome

OBJECTIVE

To determine, using ultra-high field magnetic resonance imaging (MRI), whether changes in iron content occur in the earliest phases of demyelinating disease, by quantifying the magnetic susceptibility of deep grey matter structures in patients with Clinically Isolated Syndrome (CIS) that is suggestive of multiple sclerosis (MS), as compared with age-matched healthy subjects.

METHODS

We compared 19 CIS patients to 20 age-matched, healthy controls. Scanning of the study subjects was performed on a 7T Philips Achieva system, using a 3-dimensional, T2*-weighted gradient echo acquisition. Phase data were first high-pass filtered, using a dipole fitting method, and then inverted to produce magnetic susceptibility maps. Region of interest (ROI) analysis was used to estimate magnetic susceptibility values for deep grey matter structures (caudate nucleus, putamen, globus pallidus, the thalamus and its pulvinar).

RESULTS

Significantly increased relative susceptibilities were found in the CIS group, compared with controls, for the caudate nucleus (p = < 0.01), putamen (p < 0.01), globus pallidus (p < 0.01) and pulvinar (p < 0.05). We found no significant nor consistent trends in the relationship between susceptibility and age for either the study controls or CIS patients, in any ROI (r(2) < 0.5; p > 0.05). In CIS patients, the time elapsed since the clinical event and the Expanded Disability Status Scale (EDSS) scores were not correlated with iron levels in any ROI (r(2) < 0.5; p > 0.05); however, a moderate correlation (r(2) = 0.3; p < 0.01) was found between the T1 lesion load and the mean susceptibility of the caudate nucleus.

CONCLUSION

CIS patients showed an increased iron accumulation, as measured using susceptibility mapping of the deep grey matter, suggesting that iron changes did occur at the earlier stages of CIS disease.

Susceptibility mapping as a means to visualize veins and quantify oxygen saturation

PURPOSE

To create an orientation-independent, 3D reconstruction of the veins in the brain using susceptibility mapping.

MATERIALS AND METHODS

High-resolution, high-pass filtered phase images usually used for susceptibility weighted imaging (SWI) were used as a source for local magnetic field behavior. These images were subsequently postprocessed using an inverse procedure to generate susceptibility maps of the veins. Regularization and interpolation of the data in k-space of the phase images were used to reduce reconstruction artifacts. To understand the effects of artifacts, and to fine-tune the methodology, simulations of blood vessels were performed with and without noise.

RESULTS

With sufficient resolution, major veins in the brain could be visualized with this approach. The usual geometry-dependent phase dipole effects are removed by this processing, leaving basically images of the veins. Different sized vessels show a different level of contrast depending on their partial volume effects. Vessels that are 8 mm or 16 mm in size show quantitative values expected for normal oxygen saturation levels. Smaller vessels show smaller values due to errors in the methodology and due to partial volume effects. Larger vessels show a bias toward a reduced susceptibility approaching 90% of the expected value. Limitations of the method and artifacts related to different sources of errors are demonstrated.

CONCLUSION

Susceptibility maps can successfully create venograms of the brain with varying levels of contrast-to-noise depending on the size of the vessel. Partial volume effects render this approach more useful as an imaging tool or a visualization tool, although certain larger vessels have measured susceptibilities close to expected values associated with normal blood oxygen saturation levels.