Quantitative susceptibility mapping of multiple sclerosis lesions at various ages

Multi-atlas tool for automated segmentation of brain gray matter nuclei and quantification of their magnetic susceptibility

Quantification of tissue magnetic susceptibility using MRI offers a non-invasive measure of important tissue components in the brain, such as iron and myelin, potentially providing valuable information about normal and pathological conditions during aging. Despite many advances made in recent years on imaging techniques of quantitative susceptibility mapping (QSM), accurate and robust automated segmentation tools for QSM images that can help generate universal and sharable susceptibility measures in a biologically meaningful set of structures are still not widely available. In the present study, we developed an automated process to segment brain nuclei and quantify tissue susceptibility in these regions based on a susceptibility multi-atlas library, consisting of 10 atlases with T1-weighted images, gradient echo (GRE) magnitude images and QSM images of brains with different anatomic patterns. For each atlas in this library, 10 regions of interest in iron-rich deep gray matter structures that are better defined by QSM contrast were manually labeled, including caudate, putamen, globus pallidus internal/external, thalamus, pulvinar, subthalamic nucleus, substantia nigra, red nucleus and dentate nucleus in both left and right hemispheres. We then tested different pipelines using different combinations of contrast channels to bring the set of labels from the multi-atlases to each target brain and compared them with the gold standard manual delineation. The results showed that the segmentation accuracy using dual contrasts QSM/T1 pipeline outperformed other dual-contrast or single-contrast pipelines. The dice values of 0.77 ± 0.09 using the QSM/T1 multi-atlas pipeline rivaled with the segmentation reliability obtained from multiple evaluators with dice values of 0.79 ± 0.07 and gave comparable or superior performance in segmenting subcortical nuclei in comparison with standard FSL FIRST or recent multi-atlas package of volBrain. The segmentation performance of the QSM/T1 multi-atlas was further tested on QSM images acquired using different acquisition protocols and platforms and showed good reliability and reproducibility with average dice of 0.79 ± 0.08 to manual labels and 0.89 ± 0.04 in an inter-protocol manner. The extracted quantitative magnetic susceptibility values in the deep gray matter nuclei also correlated well between different protocols with inter-protocol correlation constants all larger than 0.97. Such reliability and performance was ultimately validated in an external dataset acquired at another study site with consistent susceptibility measures obtained using the QSM/T1 multi-atlas approach in comparison to those using manual delineation. In summary, we designed a susceptibility multi-atlas tool for automated and reliable segmentation of QSM images and for quantification of magnetic susceptibilities. It is publicly available through our cloud-based platform (www.mricloud.org). Further improvement on the performance of this multi-atlas tool is expected by increasing the number of atlases in the future.

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.

Reproducibility of quantitative susceptibility mapping in lumbar vertebra

Background

To evaluate the reliability and reproducibility of quantitative susceptibility mapping (QSM) in the lumbar vertebra.

Methods

From May 2017 to September 2017, 61 subjects who underwent QSM MRI and quantitative computed tomography (QCT) were consecutively enrolled in this prospective study. QSM examination was performed two times with an interval of less than 1 week for each subject. For each data set, the QSM and QCT values on L1-L4 vertebral bodies were measured independently by two radiologists. The correlation coefficient between QSM and QCT values was calculated on L1-L4 vertebral bodies. The intraclass correlation coefficient (ICC) and Bland-Altman plots were used to evaluate the inter-observer reliability and the inter-scan reproducibility on QSM.

Results

A total of 61 subjects (mean age, 55.5±13.7 years) with 244 vertebral bodies were analyzed. Overall, QSM and QCT showed good correlation in the L1-L4 vertebral body, especially in the L3 (R=-0.75). QSM value showed excellent inter-observer reliability (ICC, 0.992, 95% CI: 0.985-0.996) with a mean difference of 0.35 and 95% limits of agreements of within -22.74 to 23.45 ppb, and very good inter-scan reproducibility (ICC, 0.932, 95% CI: 0.886-0.959) with a mean difference of -7.60 ppb and 95% limits of agreements of within of -92.85 to 77.62 ppb.

Conclusions

QSM in the lumbar vertebra is a reliable and reproducible technique for evaluating bone mineral density.

Ultrashort Echo Time Quantitative Susceptibility Mapping (UTE-QSM) of Highly Concentrated Magnetic Nanoparticles: A Comparison Study about Different Sampling Strategies

The ability to accurately and non-invasively quantify highly concentrated magnetic nanoparticles (MNPs) is desirable for many emerging applications. Ultrashort echo time quantitative susceptibility mapping (UTE-QSM) has demonstrated the capability to detect high iron concentrations. In this study, we aimed to investigate the effect of different sampling trajectories on the accuracy of quantification based on MNPs acquired through UTE-QSM. A phantom with six different MNP concentrations was prepared for UTE-QSM study with different UTE sampling trajectories, including radial acquisition, continuous single point imaging (CSPI), and Cones with four different gradient stretching factors of 1.0, 1.2, 1.4, and 1.6. No significant differences were found in QSM values derived from the different UTE sampling strategies, suggesting that the UTE-QSM technique could be accelerated with extended Cones sampling.

Imaging outcome measures of neuroprotection and repair in MS: A consensus statement from NAIMS

OBJECTIVE

To summarize current and emerging imaging techniques that can be used to assess neuroprotection and repair in multiple sclerosis (MS), and to provide a consensus opinion on the potential utility of each technique in clinical trial settings.

METHODS

Clinicians and scientists with expertise in the use of MRI in MS convened in Toronto, Canada, in November 2016 at a North American Imaging in Multiple Sclerosis (NAIMS) Cooperative workshop meeting. The discussion was compiled into a manuscript and circulated to all NAIMS members in attendance. Edits and feedback were incorporated until all authors were in agreement.

RESULTS

A wide spectrum of imaging techniques and analysis methods in the context of specific study designs were discussed, with a focus on the utility and limitations of applying each technique to assess neuroprotection and repair. Techniques were discussed under specific themes, and included conventional imaging, magnetization transfer ratio, diffusion tensor imaging, susceptibility-weighted imaging, imaging cortical lesions, magnetic resonance spectroscopy, PET, advanced diffusion imaging, sodium imaging, multimodal techniques, imaging of special regions, statistical considerations, and study design.

CONCLUSIONS

Imaging biomarkers of neuroprotection and repair are an unmet need in MS. There are a number of promising techniques with different strengths and limitations, and selection of a specific technique will depend on a number of factors, notably the question the trial seeks to answer. Ongoing collaborative efforts will enable further refinement and improved methods to image the effect of novel therapeutic agents that exert benefit in MS predominately through neuroprotective and reparative mechanisms.

Iron related changes in MS lesions and their validity to characterize MS lesion types and dynamics with Ultra-high field magnetic resonance imaging

Iron accumulates with age in the normal human brain. This process is altered at several levels in the brain of multiple sclerosis (MS) patients. Since iron is mainly stored in oligodendrocytes and myelin in the normal brain, its liberation in demyelinating lesions may amplify tissue damage in demyelinating lesions and its uptake in macrophages and microglia may help to more precisely define activity stages of the lesions. In addition, glia cells change their iron import, export and storage properties in MS lesions, which is reflected by alterations in the expression of iron transport molecules. Changes of iron distribution in the brain can be reliably detected by MRI, particularly upon application of Ultra‐high magnetic field (7 Tesla). Iron‐sensitive MRI allows to more accurately distinguish the lesions in MS from those in other inflammatory brain diseases, to visualize a subset of slowly expanding lesions in the progressive stage of MS and to increase the sensitivity for lesion detection in the gray matter, such as the cerebral cortex or deep gray matter nuclei.

Quantitative susceptibility mapping identifies inflammation in a subset of chronic multiple sclerosis lesions

Chronic active multiple sclerosis lesions, characterized by a hyperintense rim of iron-enriched, activated microglia and macrophages, have been linked to greater tissue damage. Post-mortem studies have determined that chronic active lesions are primarily related to the later stages of multiple sclerosis; however, the occurrence of these lesions, and their relationship to earlier disease stages may be greatly underestimated. Detection of chronic active lesions across the patient spectrum of multiple sclerosis requires a validated imaging tool to accurately identify lesions with persistent inflammation. Quantitative susceptibility mapping provides efficient in vivo quantification of susceptibility changes related to iron deposition and the potential to identify lesions harbouring iron-laden inflammatory cells. The PET tracer 11C-PK11195 targets the translocator protein expressed by activated microglia and infiltrating macrophages. Accordingly, this study aimed to validate that lesions with a hyperintense rim on quantitative susceptibility mapping from both relapsing and progressive patients demonstrate a higher level of innate immune activation as measured on 11C-PK11195 PET. Thirty patients were enrolled in this study, 24 patients had relapsing remitting multiple sclerosis, six had progressive multiple sclerosis, and all patients had concomitant MRI with a gradient echo sequence and PET with 11C-PK11195. A total of 406 chronic lesions were detected, and 43 chronic lesions with a hyperintense rim on quantitative susceptibility mapping were identified as rim+ lesions. Susceptibility (relative to CSF) was higher in rim+ (2.42 ± 17.45 ppb) compared to rim- lesions (-14.6 ± 19.3 ppb, P < 0.0001). Among rim+ lesions, susceptibility within the rim (20.04 ± 14.28 ppb) was significantly higher compared to the core (-5.49 ± 14.44 ppb, P < 0.0001), consistent with the presence of iron. In a mixed-effects model, 11C-PK11195 uptake, representing activated microglia/macrophages, was higher in rim+ lesions compared to rim- lesions (P = 0.015). Validating our in vivo imaging results, multiple sclerosis brain slabs were imaged with quantitative susceptibility mapping and processed for immunohistochemistry. These results showed a positive translocator protein signal throughout the expansive hyperintense border of rim+ lesions, which co-localized with iron containing CD68+ microglia and macrophages. In conclusion, this study provides evidence that suggests that a hyperintense rim on quantitative susceptibility measure within a chronic lesion is a correlate for persistent inflammatory activity and that these lesions can be identified in the relapsing patients. Utilizing quantitative susceptibility measure to differentiate chronic multiple sclerosis lesion subtypes, especially chronic active lesions, would provide a method to assess the impact of these lesions on disease progression.

Slowly expanding/evolving lesions as a magnetic resonance imaging marker of chronic active multiple sclerosis lesions

Background:

Chronic lesion activity driven by smoldering inflammation is a pathological hallmark of progressive forms of multiple sclerosis (MS).

Objective:

To develop a method for automatic detection of slowly expanding/evolving lesions (SELs) on conventional brain magnetic resonance imaging (MRI) and characterize such SELs in primary progressive MS (PPMS) and relapsing MS (RMS) populations.

Methods:

We defined SELs as contiguous regions of existing T2 lesions showing local expansion assessed by the Jacobian determinant of the deformation between reference and follow-up scans. SEL candidates were assigned a heuristic score based on concentricity and constancy of change in T2- and T1-weighted MRIs. SELs were examined in 1334 RMS patients and 555 PPMS patients.

Results:

Compared with RMS patients, PPMS patients had higher numbers of SELs (p = 0.002) and higher T2 volumes of SELs (p < 0.001). SELs were devoid of gadolinium enhancement. Compared with areas of T2 lesions not classified as SEL, SELs had significantly lower T1 intensity at baseline and larger decrease in T1 intensity over time.

Conclusion:

We suggest that SELs reflect chronic tissue loss in the absence of ongoing acute inflammation. SELs may represent a conventional brain MRI correlate of chronic active MS lesions and a candidate biomarker for smoldering inflammation in MS.

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.

Characterization of normal-appearing white matter in multiple sclerosis using quantitative susceptibility mapping in conjunction with diffusion tensor imaging

PURPOSE

Quantitative susceptibility mapping (QSM) is influenced by iron as well as myelin, which makes interpretation of pathologic changes challenging. Concurrent acquisition of MR sequences that are sensitive to axonal/myelin integrity, such as diffusion tensor imaging (DTI), may provide context for interpreting quantitative susceptibility (QS) signal. The purpose of our study was to investigate alterations in normal-appearing white matter (NAWM) in multiple sclerosis (MS) using QSM in conjunction with DTI.

METHODS

Twenty relapsing-remitting MS patients and 20 age-matched healthy controls (HC) were recruited for this prospective study. QS, radial diffusivity (RD), fractional anisotropy (FA), and R2* maps within the whole brain as well as individual tracts were generated for comparison between NAWM and HC white matter (HCWM).

RESULTS

MS lesions demonstrated significant differences in QS, FA, RD, and R2* compared to HCWM (p < 0.03). These metrics did not show a significant difference between whole-brain NAWM and HCWM. Among NAWM tracts, the cingulate gyri demonstrated significantly decreased QS compared to HCWM (p = 0.004). The forceps major showed significant differences in FA and RD without corresponding changes in QS (p < 0.01).

CONCLUSION

We found discordant changes in QSM and DTI metrics within the cingulate gyri and forceps major. This may potentially reflect the influence of paramagnetic substrates such as iron, which could be decreased along these NAWM tracts. Our results point to the potential role of QSM as a unique biomarker, although additional validation studies are needed.

A Radiodensity Histogram Study of the Brain in Multiple Sclerosis

Multiple sclerosis (MS) is a progressive neurodegenerative disease, affecting 1 million Americans and 2.5 million people globally. Although the diagnosis is made clinically, imaging plays a major role in diagnosing and monitoring disease progression and treatment response. Magnetic resonance imaging (MRI) has proven sensitive in imaging MS lesions, but the characterization offered by routine clinical MRI remains qualitative and with discrepancies between imaging and clinical findings. We investigated the ability of digital analysis of noncontrast head computed tomography (CT) images to detect global brain changes of MS. All routine diagnostic head CTs obtained on patients with known MS obtained from 1 of 2 scan platforms from 6/1/2011 to 6/1/2015 were reviewed. Head CT images from 54 patients with MS met inclusion criteria. Head CT images were processed and histogram metrics were compared to age- and gender- matched control subjects from the same CT scanners during the same time interval. Histogram metrics were correlated with plaque burden as seen on MRI studies. Compared with control subjects, patients had increased total brain radiodensity (P < .0001), further characterized as an increased histogram modal radiodensity (P < .0001) with decrease in histogram skewness (P < .0001). Radiodensity decreased with increasing plaque burden. Similar findings were seen in the patients with only mild plaque burden sub- group. Radiodensity is a unique tissue metric that is not measured by other imaging techniques. Our study finds that brain radiodensity histogram metrics highly correlate with MS, even in cases with minimal plaque burden.

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.

True phase quantitative susceptibility mapping using continuous single-point imaging: a feasibility study

PURPOSE

In this study, we explore the feasibility of a new imaging scheme for quantitative susceptibility mapping (QSM): continuous single-point imaging (CSPI), which uses a pure phase encoding strategy to achieve true phase imaging and improve QSM accuracy.

METHODS

The proposed CSPI is a modification of conventional SPI to allow acquisition of multiple echoes in a single scan. Immediately following a phase encoding gradient, the free induction decay is continuously sampled with extremely high temporal resolution to obtain k-space data at a fixed spatial frequency (i.e., at a fixed k-space coordinate). By having near-0 readout duration, CSPI results in a true snapshot of the transverse magnetization at each TE. Additionally, parallel imaging with autocalibration is utilized to reduce scan time, and an optional temporal averaging strategy is presented to improve signal-to-noise ratio for objects with low proton density or short T2* decay. The reconstructed CSPI images were input to a QSM framework based on morphology enabled dipole inversion.

RESULT

In an experiment performed using iron phantoms, susceptibility estimated using CSPI showed high linearity (R² = 0.9948) with iron concentration. Additionally, reconstructed CSPI phase images showed much reduced ringing artifact compared with phase images obtained using a frequency encoding strategy. In an ex vivo experiment performed using human tibia samples, estimated susceptibilities ranged from -1.6 to -2.1 ppm, in agreement with values reported in the literature (ranging from -1.2 to -2.2 ppm).

CONCLUSION

We have demonstrated the feasibility of using CSPI to obtain true phase images for QSM.

Region-specific susceptibility change in cognitively impaired patients with diabetes mellitus

Emerging evidence suggests that diabetes mellitus (DM) is associated with iron and calcium metabolism. However, few studies have investigated the presence of DM in cognitively impaired patients and its effect on brain iron and calcium accumulation. Therefore, we assessed the effects of DM on cognitively impaired patients using quantitative susceptibility mapping (QSM). From June 2012 to Feb 2014, 92 eligible cognitively impaired patients underwent 3T magnetic resonance imaging (MRI). There were 46 patients with DM (DM+) and 46 aged matched patients without DM (DM-). QSM was obtained from gradient echo data and analyzed by drawing regions of interest around relevant anatomical structures. Clinical factors and vascular pathology were also evaluated. Measurement differences between DM+ and DM- patients were assessed by t tests. A multiple regression analysis was performed to identify independent predictors of magnetic susceptibility. DM+ patients showed lower susceptibility values, indicative of lower brain iron content, than DM- patients, which was significant in the hippocampus (4.80 ± 8.31 ppb versus 0.22 ± 10.60 ppb, p = 0.024) and pulvinar of the thalamus (36.30 ± 19.88 ppb versus 45.90 ± 20.02 ppb, p = 0.023). On multiple regression analysis, microbleed number was a predictor of susceptibility change in the hippocampus (F = 4.291, beta = 0.236, p = 0.042) and DM was a predictor of susceptibility change in the pulvinar of the thalamus (F = 4.900, beta = - 0.251, p = 0.030). In cognitively impaired patients, presence of DM was associated with lower susceptibility change in the pulvinar of the thalamus and hippocampus. This suggests that there may be region-specific alterations of calcium deposition in cognitively impaired subjects with DM.

Background field removal for susceptibility mapping of human brain with large susceptibility variations

PURPOSE

In quantitative susceptibility mapping (QSM) of human brain, the background field induced by air-tissue interface varies significantly with respect to the rotation angle between the head and the static field, which may result in substantial error in the estimated magnetic susceptibility values. The goal of this study was to develop a strategy to better remove such orientation dependent background field.

METHODS

An improved background field removal method is proposed based on the sophisticated harmonic artifact reduction for phase data using a region adaptive kernel (R-SHARP), named iRSHARP. It uses a spatially weighted spherical Gaussian kernel exploiting the amplitude, gradient, and wrap count of the phase map. The method was validated using both numerical simulations and in vivo human brain data at multiple head orientations. Performance was compared with the variable kernel (V-SHARP) and R-SHARP methods.

RESULTS

The proposed iRSHARP method showed improved background removal over R-SHARP while cutting the computational time in half. As compared to V-SHARP and R-SHARP, the iRSHARP generated local field and susceptibility maps showed fewer artifacts in regions of large susceptibility variations, and for the in vivo human brain, the susceptibilities of the deep gray matter nuclei were consistent with the in vivo gold-standard "Calculation of Susceptibility through Multiple Orientation Sampling" (COSMOS) values.

CONCLUSION

iRSHARP can remove the orientation dependent background field effectively. Using iRSHARP, the paranasal sinus regions can be preserved in the brain mask and the brain integrity was conserved, which may facilitate further data analysis and clinical application.

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.

Longitudinal study of multiple sclerosis lesions using ultra-high field (7T) multiparametric MR imaging

Pathophysiology of multiple sclerosis (MS) lesions is dynamic and changes over time. The purpose of this exploratory study was to determine the longitudinal changes in MS lesions over time on ultra-high field MR imaging. Nine patients with MS underwent high-resolution 3D-susceptibility weighted imaging (SWI) and 2D-gradient-echo-T2*-weighted imaging on 7T MRI at baseline and after ~2.4 years of follow-up. Morphologic imaging characteristics, signal intensity patterns and quantitative susceptibility mapping (QSM) values of lesions were recorded at both time points. Lesions were classified as "iron-laden" if they demonstrated hypointense signal on T2*-weighted images and/or SWI as well as hyperintense signal on QSM. Lesions were considered "non-iron-laden" if they were hyperintense on T2*/SWI and isointense or hyperintense on QSM. Total of 162 non-iron-laden and 29 iron-laden lesions were observed at baseline. No change in baseline lesion size during follow up was recorded in 92.7%; no change in lesion-vessel relationship in 86.5%; and no change in signal intensity pattern in 96.9% of lesions. Three lesions which were non-iron-laden at baseline, exhibited iron at follow-up. In two iron-laden lesions, redistribution of iron content was observed at follow-up. Two-thirds of these iron-laden lesions showed an increase in QSM at follow-up relative to baseline, and the remaining one-third exhibited decrease in QSM. Most of the newly formed lesions (11/13, 84.6%) at follow-up were iron-laden. 7T multiparametric MRI is a useful tool for tracking the evolution of MS lesions, especially with regard to changes in iron content.

Mask-Adapted Background Field Removal for Artifact Reduction in Quantitative Susceptibility Mapping of the Prostate

We propose an alternative processing method for quantitative susceptibility mapping of the prostate that reduces artifacts and enables better visibility and quantification of calcifications and other lesions. Three-dimensional gradient-echo magnetic resonance data were obtained from 26 patients at 3 T who previously received a planning computed tomography of the prostate. Phase images were unwrapped using Laplacian-based phase unwrapping. The background field was removed with the V-SHARP method using tissue masks for the entire abdomen (Method 1) and masks that excluded bone and the rectum (Method 2). Susceptibility maps were calculated with the iLSQR method. The quality of susceptibility maps was assessed by one radiologist and two physicists who rated the data for visibility of lesions and data quality on a scale from 1 (poor) to 4 (good). The readers rated susceptibility maps computed with Method 2 to be, on average, better for visibility of lesions with a score of 2.9 ± 1.1 and image quality with a score of 2.8 ± 0.8 compared with maps computed with Method 1 (2.4 ± 1.2/2.3 ± 1.0). Regarding strong artifacts, these could be removed using adapted masks, and the susceptibility values seemed less biased by the artifacts. Thus, using an adapted mask for background field removal when calculating susceptibility maps of the prostate from phase data reduces artifacts and improves visibility of lesions.

Iron deposition quantification: Applications in the brain and liver

Iron has long been implicated in many neurological and other organ diseases. It is known that over and above the normal increases in iron with age, in certain diseases there is an excessive iron accumulation in the brain and liver. MRI is a noninvasive means by which to image the various structures in the brain in three dimensions and quantify iron over the volume of the object of interest. The quantification of iron can provide information about the severity of iron-related diseases as well as quantify changes in iron for patient follow-up and treatment monitoring. This article provides an overview of current MRI-based methods for iron quantification, specifically for the brain and liver, including: signal intensity ratio, R₂ , R2*, R2', phase, susceptibility weighted imaging and quantitative susceptibility mapping (QSM). Although there are numerous approaches to measuring iron, R₂ and R2* are currently preferred methods in imaging the liver and QSM has become the preferred approach for imaging iron in the brain.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2018. J. MAGN. RESON. IMAGING 2018;48:301-317.

Signal Change of Acute Cortical and Juxtacortical Microinfarction on Follow-Up MRI

BACKGROUND AND PURPOSE

Although the clinical importance of cortical microinfarcts has become well-recognized recently, the evolution of cortical microinfarcts on MR imaging is not fully understood. The aim of this study was to examine the temporal changes in acute cortical microinfarcts using susceptibility-weighted imaging and conventional MR imaging.

MATERIALS AND METHODS

Patients with acute infarcts located in the cortical and/or juxtacortical region measuring ≤10 mm in axial diameter based on diffusion-weighted imaging who had a follow-up 3T MR imaging were retrospectively included in the study. All lesions did not show hypointensity on initial T2*WI. For cortical and/or juxtacortical microinfarcts detected on initial DWI, 2 neuroradiologists evaluated the follow-up MR imaging (T2WI, FLAIR, T2*WI, and SWI) and assessed lesion signal intensities and locations (cortical microinfarcts or microinfarcts with juxtacortical white matter involvement).

RESULTS

On initial DWI, 2 radiologists observed 180 cortical and/or juxtacortical microinfarcts in 35 MR imaging examinations in 25 patients; on follow-up, the neuroradiologists identified 29 cortical microinfarcts (16%) on T2WI, 9 (5%) on FLAIR, 4 (2%) on T2*, and 97 (54%) on SWI. All cortical microinfarcts detected with any follow-up MR imaging showed hyperintensity on T2WI/FLAIR and/or hypointensity on T2*WI and SWI.

CONCLUSIONS

SWI revealed conversion (paramagnetic susceptibility changes) of acute cortical microinfarcts, suggesting that a substantial number of cortical microinfarcts may contain hemorrhagic components.

Quantifying iron content in magnetic resonance imaging

Measuring iron content has practical clinical indications in the study of diseases such as Parkinson's disease, Huntington's disease, ferritinopathies and multiple sclerosis as well as in the quantification of iron content in microbleeds and oxygen saturation in veins. In this work, we review the basic concepts behind imaging iron using T2, T2*, T2', phase and quantitative susceptibility mapping in the human brain, liver and heart, followed by the applications of in vivo iron quantification in neurodegenerative diseases, iron tagged cells and ultra-small superparamagnetic iron oxide (USPIO) nanoparticles.

A new discrete dipole kernel for quantitative susceptibility mapping

PURPOSE

Most approaches for quantitative susceptibility mapping (QSM) are based on a forward model approximation that employs a continuous Fourier transform operator to solve a differential equation system. Such formulation, however, is prone to high-frequency aliasing. The aim of this study was to reduce such errors using an alternative dipole kernel formulation based on the discrete Fourier transform and discrete operators.

METHODS

The impact of such an approach on forward model calculation and susceptibility inversion was evaluated in contrast to the continuous formulation both with synthetic phantoms and in vivo MRI data.

RESULTS

The discrete kernel demonstrated systematically better fits to analytic field solutions, and showed less over-oscillations and aliasing artifacts while preserving low- and medium-frequency responses relative to those obtained with the continuous kernel. In the context of QSM estimation, the use of the proposed discrete kernel resulted in error reduction and increased sharpness.

CONCLUSION

This proof-of-concept study demonstrated that discretizing the dipole kernel is advantageous for QSM. The impact on small or narrow structures such as the venous vasculature might by particularly relevant to high-resolution QSM applications with ultra-high field MRI - a topic for future investigations. The proposed dipole kernel has a straightforward implementation to existing QSM routines.

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.

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

Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.Quantitative susceptibility mapping (QSM) solves the magnetic field-to-magnetization (tissue susceptibility) inverse problem under conditions of noisy and incomplete field data acquired using magnetic resonance imaging. Therefore, sophisticated algorithms are necessary to treat the ill-posed nature of the problem and are reviewed here. The forward problem is typically presented as an integral form, where the field is the convolution of the dipole kernel and tissue susceptibility distribution. This integral form can be equivalently written as a partial differential equation (PDE). Algorithmic challenges are to reduce streaking and shadow artifacts characterized by the fundamental solution of the PDE. Bayesian maximum a posteriori estimation can be employed to solve the inverse problem, where morphological and relevant biomedical knowledge (specific to the imaging situation) are used as priors. As the cost functions in Bayesian QSM framework are typically convex, solutions can be robustly computed using a gradient-based optimization algorithm. Moreover, one can not only accelerate Bayesian QSM, but also increase its effectiveness at reducing shadows using prior knowledge based preconditioners. Improving the efficiency of QSM is under active development, and a rigorous analysis of preconditioning needs to be carried out for further investigation.

Significance and In Vivo Detection of Iron-Laden Microglia in White Matter Multiple Sclerosis Lesions

Microglia are resident immune cells that fulfill protective and homeostatic functions in the central nervous system (CNS) but may also promote neurotoxicity in the aged brain and in chronic disease. In multiple sclerosis (MS), an autoimmune demyelinating disease of the CNS, microglia and macrophages contribute to the development of white matter lesions through myelin phagocytosis, and possibly to disease progression through diffuse activation throughout myelinated white matter. In this review, we discuss an additional compartment of myeloid cell activation in MS, i.e., the rim and normal adjacent white matter of chronic active lesions. In chronic active lesions, microglia and macrophages may contain high amounts of iron, express markers of proinflammatory polarization, are activated for an extended period of time (years), and drive chronic tissue damage. Iron-positive myeloid cells can be visualized and quantified with quantitative susceptibility mapping (QSM), a magnetic resonance imaging technique. Thus, QSM has potential as an in vivo biomarker for chronic inflammatory activity in established white matter MS lesions. Reducing chronic inflammation associated with iron accumulation using existing or novel MS therapies may impact disease severity and progression.

Mapping of thalamic magnetic susceptibility in multiple sclerosis indicates decreasing iron with disease duration: A proposed mechanistic relationship between inflammation and oligodendrocyte vitality

Recent advances in susceptibility MRI have dramatically improved the visualization of deep gray matter brain regions and the quantification of their magnetic properties in vivo, providing a novel tool to study the poorly understood iron homeostasis in the human brain. In this study, we used an advanced combination of the recent quantitative susceptibility mapping technique with dedicated analysis methods to study intra-thalamic tissue alterations in patients with clinically isolated syndrome (CIS) and multiple sclerosis (MS). Thalamic pathology is one of the earliest hallmarks of MS and has been shown to correlate with cognitive dysfunction and fatigue, but the mechanisms underlying the thalamic pathology are poorly understood. We enrolled a total of 120 patients, 40 with CIS, 40 with Relapsing Remitting MS (RRMS), and 40 with Secondary Progressive MS (SPMS). For each of the three patient groups, we recruited 40 controls, group matched for age- and sex (120 total). We acquired quantitative susceptibility maps using a single-echo gradient echo MRI pulse sequence at 3 T. Group differences were studied by voxel-based analysis as well as with a custom thalamus atlas. We used threshold-free cluster enhancement (TFCE) and multiple regression analyses, respectively. We found significantly reduced magnetic susceptibility compared to controls in focal thalamic subregions of patients with RRMS (whole thalamus excluding the pulvinar nucleus) and SPMS (primarily pulvinar nucleus), but not in patients with CIS. Susceptibility reduction was significantly associated with disease duration in the pulvinar, the left lateral nuclear region, and the global thalamus. Susceptibility reduction indicates a decrease in tissue iron concentration suggesting an involvement of chronic microglia activation in the depletion of iron from oligodendrocytes in this central and integrative brain region. Not necessarily specific to MS, inflammation-mediated iron release may lead to a vicious circle that reduces the protection of axons and neuronal repair.

Diagnostic accuracy of semiautomatic lesion detection plus quantitative susceptibility mapping in the identification of new and enhancing multiple sclerosis lesions

Purpose

To evaluate the diagnostic accuracy of a novel non-contrast brain MRI method based on semiautomatic lesion detection using T2w FLAIR subtraction image, the statistical detection of change (SDC) algorithm (T2w + SDC), and quantitative susceptibility mapping (QSM). This method identifies new lesions and discriminates between enhancing and nonenhancing lesions in multiple sclerosis (MS).

Methods

Thirty three MS patients who had MRIs at two different time points with at least one new Gd-enhancing lesion on the 2nd MRI were included in the study. For a reference standard, new lesions were identified by two neuroradiologists on T2w and post-Gd T1w images with the help of T2w + SDC. The diagnostic accuracy of the proposed method based on QSM and T2w + SDC lesion detection (T2w + SDC + QSM) for assessing lesion enhancement status was determined. Receiver operating characteristic (ROC) analysis was performed to compute the optimal lesion susceptibility cutoff value.

Results

A total of 165 new lesions (54 enhancing, 111 nonenhancing) were identified. The sensitivity and specificity of T2w + SDC + QSM in predicting lesion enhancement status were 90.7% and 85.6%, respectively. For lesions ≥50 mm³, ROC analysis showed an optimal QSM cutoff value of 13.5 ppb with a sensitivity of 88.4% and specificity of 88.6% (0.93, 95% CI, 0.87–0.99). For lesions ≥15 mm³, the optimal QSM cutoff was 15.4 ppb with a sensitivity of 77.9% and specificity of 94.0% (0.93, 95% CI, 0.89–0.97).

Conclusion

The proposed T2w + SDC + QSM method is highly accurate for identifying and predicting the enhancement status of new MS lesions without the use of Gd injection.

•

T2w + SDC has high sensitivity and accuracy in detecting new MS lesions.

•

T2w + SDC + QSM is highly accurate in discriminating between new enhancing and new nonenhancing lesions.

•

T2w + SDC + QSM can form the basis of an imaging protocol without Gadolinium injection for routine surveillance of MS patients.

Combining Quantitative Susceptibility Mapping with Automatic Zero Reference (QSM0) and Myelin Water Fraction Imaging to Quantify Iron-Related Myelin Damage in Chronic Active MS Lesions

BACKGROUND AND PURPOSE

A hyperintense rim on susceptibility in chronic MS lesions is consistent with iron deposition, and the purpose of this study was to quantify iron-related myelin damage within these lesions as compared with those without rim.

MATERIALS AND METHODS

Forty-six patients had 2 longitudinal quantitative susceptibility mapping with automatic zero reference scans with a mean interval of 28.9 ± 11.4 months. Myelin water fraction mapping by using fast acquisition with spiral trajectory and T2 prep was obtained at the second time point to measure myelin damage. Mixed-effects models were used to assess lesion quantitative susceptibility mapping and myelin water fraction values.

RESULTS

Quantitative susceptibility mapping scans were on average 6.8 parts per billion higher in 116 rim-positive lesions compared with 441 rim-negative lesions (P < .001). All rim-positive lesions retained a hyperintense rim over time, with increasing quantitative susceptibility mapping values of both the rim and core regions (P < .001). Quantitative susceptibility mapping scans and myelin water fraction in rim-positive lesions decreased from rim to core, which is consistent with rim iron deposition. Whole lesion myelin water fractions for rim-positive and rim-negative lesions were 0.055 ± 0.07 and 0.066 ± 0.04, respectively. In the mixed-effects model, rim-positive lesions had on average 0.01 lower myelin water fraction compared with rim-negative lesions (P < .001). The volume of the rim at the initial quantitative susceptibility mapping scan was negatively associated with follow-up myelin water fraction (P < .01).

CONCLUSIONS

Quantitative susceptibility mapping rim-positive lesions maintained a hyperintense rim, increased in susceptibility, and had more myelin damage compared with rim-negative lesions. Our results are consistent with the identification of chronic active MS lesions and may provide a target for therapeutic interventions to reduce myelin damage.

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.

The Role of Advanced Magnetic Resonance Imaging Techniques in Multiple Sclerosis Clinical Trials

Magnetic resonance imaging has been crucial in the development of anti-inflammatory disease-modifying treatments. The current landscape of multiple sclerosis clinical trials is currently expanding to include testing not only of anti-inflammatory agents, but also neuroprotective, remyelinating, neuromodulating, and restorative therapies. This is especially true of therapies targeting progressive forms of the disease where neurodegeneration is a prominent feature. Imaging techniques of the brain and spinal cord have rapidly evolved in the last decade to permit in vivo characterization of tissue microstructural changes, connectivity, metabolic changes, neuronal loss, glial activity, and demyelination. Advanced magnetic resonance imaging techniques hold significant promise for accelerating the development of different treatment modalities targeting a variety of pathways in MS.

Susceptibility-sensitive MRI of multiple sclerosis lesions and the impact of normal-appearing white matter changes

Susceptibility-sensitive magnetic resonance imaging (MRI) has gained importance in multiple sclerosis (MS) research because of its versatility, high resolution and excellent sensitivity to changes in tissue structure and composition. In particular, mapping of the resonance frequency of the MR signal and quantitative susceptibility mapping (QSM) have been explored for the description of MS lesions. Many current studies utilizing these techniques attribute increases in the MR frequency or QSM to elevated tissue iron content, in addition to myelin loss. However, this common interpretation is inconsistent with recent histopathological studies. Here, we investigate the nature of MR frequency shifts related to MS lesions by comparing post-mortem MRI data with histology, and contrast them with numerical simulations of the MR signal. We demonstrate that iron accumulation is not the driving source of the MR frequency or QSM image contrast in our sample; rather, most chronic MS lesions are characterized by advanced loss of both myelin and iron. Moreover, our results suggest that the appearance of MS lesions on MR frequency maps and QSM depends on changes in the non-lesional white matter surrounding the lesions. Understanding and accounting for these changes is essential for the quantitative interpretation of MR frequency or QSM data in white matter.

Quantifying the Susceptibility Variation of Normal-Appearing White Matter in Multiple Sclerosis by Quantitative Susceptibility Mapping

OBJECTIVE

The purpose of this study is to evaluate the magnetic susceptibility of normal-appearing white matter (NAWM) in patients with multiple sclerosis (MS) using quantitative susceptibility mapping.

MATERIALS AND METHODS

Ninety-four patients with relapse-remitting MS (RRMS) (37 with gadolinium-enhancing lesions and 57 with only gadolinium-nonenhancing lesions) and 55 healthy control subjects were included in this retrospective study. The susceptibility values of NAWM relative to CSF in patients with MS were compared with those of white matter (WM) in healthy control subjects and were correlated with the patient status of gadolinium-enhancing lesions, disease duration, and expanded disability status scale scores.

RESULTS

All 37 patients with RRMS and gadolinium-enhancing lesions also had gadolinium-nonenhancing lesions. Susceptibility values of NAWM in patients with MS and only gadolinium-nonenhancing lesions (-18.29 ± 8.03 parts per billion [ppb]) were higher than those for WM in healthy control subjects (-25.81 ± 6.02 ppb; p < 0.001) and NAWM in patients with gadolinium-enhancing lesions (-25.64 ± 6.55 ppb; p < 0.001). Susceptibility values of NAWM in patients with MS with gadolinium-enhancing lesions were similar to those for WM in healthy control subjects (p = 0.91). This trend was dependent on neither NAWM region nor disease duration when the data were controlled for age. NAWM susceptibility was not correlated with either disease duration or expanded disability status scale (p > 0.05).

CONCLUSION

In patients with RRMS and gadolinium-nonenhancing lesions, the susceptibility values of NAWM decrease when gadolinium-enhancing lesions appear, approaching values similar to those of WM in healthy control subjects, suggesting that NAWM may contribute to the iron accumulation observed in early active MS lesions.

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.

Heterogeneity of Cortical Lesion Susceptibility Mapping in Multiple Sclerosis

BACKGROUND AND PURPOSE

Quantitative susceptibility mapping has been used to characterize iron and myelin content in the deep gray matter of patients with multiple sclerosis. Our aim was to characterize the susceptibility mapping of cortical lesions in patients with MS and compare it with neuropathologic observations.

MATERIALS AND METHODS

The pattern of microglial activation was studied in postmortem brain tissues from 16 patients with secondary-progressive MS and 5 age-matched controls. Thirty-six patients with MS underwent 3T MR imaging, including 3D double inversion recovery and 3D-echo-planar SWI.

RESULTS

Neuropathologic analysis revealed the presence of an intense band of microglia activation close to the pial membrane in subpial cortical lesions or to the WM border of leukocortical cortical lesions. The quantitative susceptibility mapping analysis revealed 131 cortical lesions classified as hyperintense; 33, as isointense; and 84, as hypointense. Quantitative susceptibility mapping hyperintensity edge found in the proximity of the pial surface or at the white matter/gray matter interface in some of the quantitative susceptibility mapping-hyperintense cortical lesions accurately mirrors the microglia activation observed in the neuropathology analysis.

CONCLUSIONS

Cortical lesion susceptibility maps are highly heterogeneous, even at individual levels. Quantitative susceptibility mapping hyperintensity edge found in proximity to the pial surface might be due to the subpial gradient of microglial activation.

Dentate nucleus iron deposition is a potential biomarker for tremor-dominant Parkinson's disease

Parkinson's disease (PD) is a heterogeneous neurodegenerative disorder with variable clinicopathologic phenotypes and underlying neuropathologic mechanisms. Each clinical phenotype has a unique set of motor symptoms. Tremor is the most frequent initial motor symptom of PD and is the most difficult symptom to treat. The dentate nucleus (DN) is a deep iron-rich nucleus in the cerebellum and may be involved in PD tremor. In this study, we test the hypothesis that DN iron may be elevated in tremor-dominant PD patients using quantitative susceptibility mapping. Forty-three patients with PD [19 tremor dominant (TD)/24 akinetic rigidity (AR) dominant] and 48 healthy gender- and age-matched controls were recruited. Multi-echo gradient echo data were collected for each subject on a 3.0-T MR system. Inter-group susceptibility differences in the bilateral DN were investigated and correlations of clinical features with susceptibility were also examined. In contrast with the AR-dominant group, the TD group was found to have increased susceptibility in the bilateral DN when compared with healthy controls. In addition, susceptibility was positively correlated with tremor score in drug-naive PD patients. These findings indicate that iron load within the DN may make an important contribution to motor phenotypes in PD. Moreover, our results suggest that TD and AR-dominant phenotypes of PD can be differentiated on the basis of the susceptibility of the DN, at least at the group level. Copyright © 2016 John Wiley & Sons, Ltd.

Quantitative susceptibility mapping at 3 T: comparison of acquisition methodologies

Quantitative susceptibility mapping (QSM) is fast becoming a routine clinical tool in the evaluation and assessment of neurological diseases. Unfortunately, there is currently no established standard scanning protocol, and it is uncertain whether different acquisition strategies alter the derived estimates of magnetic susceptibility. Here, we compare some key deep grey matter susceptibility values in healthy adults acquired from various QSM sequences using either unipolar or bipolar readout gradients, accelerated imaging or not, and gradient-warp correction or not. Four healthy adult volunteers were scanned three times each within 4 days at 3 T. The eight different QSM combinations were acquired in different randomised order for each session, and then co-registered to an anatomical atlas. The average and standard deviations of magnetic susceptibilities in the caudate, putamen, red nucleus, internal and external globus pallidus were used in a linear mixed effects model to determine the influence of the various acquisition parameters. Gradient-warp correction was the only statistically significant fixed effect (p < 0.01), but its impact was small (~5% change) compared with the overall fixed effects. The random effects coefficients (i.e. the various tissues) were statistically significant. Based on our limited multiple observations in healthy adult volunteers, the susceptibilities in deep grey matter are statistically equivalent when QSM source data are acquired with or without accelerated imaging using either unipolar or bipolar readout gradients. There is, however, a statistically meaningful, but small, difference if gradient-warp correction is used or not. Copyright © 2016 John Wiley & Sons, Ltd.

Effects of biological tissue structural anisotropy and anisotropy of magnetic susceptibility on the gradient echo MRI signal phase: theoretical background

Quantitative susceptibility mapping is a potentially powerful technique for mapping tissue magnetic susceptibility from gradient recalled echo (GRE) MRI signal phase. In this review, we present up-to-date theoretical developments in analyzing the relationships between GRE signal phase and the underlying tissue microstructure and magnetic susceptibility at the cellular level. Two important phenomena contributing to the GRE signal phase are at the focus of this review - tissue structural anisotropy (e.g. cylindrical axonal bundles in white matter) and magnetic susceptibility anisotropy. One of the most intriguing and challenging problems in this field is calculating the so-called Lorentzian contribution to the phase shift induced by the local environment - magnetized tissue structures that have dimensions smaller than the imaging voxel (e.g. cells, cellular components, blood capillaries). In this review, we briefly discuss a "standard" approach to this problem, based on introduction of an imaginary Lorentzian cavity, as well as a more recent method - the generalized Lorentzian tensor approach (GLTA) - that is based on a statistical approach and a direct solution of the magnetostatic Maxwell equations. The latter adequately accounts for both types of anisotropy: the anisotropy of magnetic susceptibility and the structural tissue anisotropy. In the GLTA the frequency shift due to the local environment is characterized by the Lorentzian tensor L^, which has a substantially different structure than the susceptibility tensor χ^. While the components of χ^ are compartmental susceptibilities "weighted" by their volume fractions, the components of L^ are weighted by specific numerical factors depending on tissue geometrical microsymmetry. In multi-compartment structures, the components of the Lorentzian tensor also depend on the compartmental relaxation properties, hence the MR pulse sequence settings. Copyright © 2016 John Wiley & Sons, Ltd.

Susceptibility tensor imaging (STI) of the brain

Susceptibility tensor imaging (STI) is a recently developed MRI technique that allows quantitative determination of orientation-independent magnetic susceptibility parameters from the dependence of gradient echo signal phase on the orientation of biological tissues with respect to the main magnetic field. By modeling the magnetic susceptibility of each voxel as a symmetric rank-2 tensor, individual magnetic susceptibility tensor elements as well as the mean magnetic susceptibility and magnetic susceptibility anisotropy can be determined for brain tissues that would still show orientation dependence after conventional scalar-based quantitative susceptibility mapping to remove such dependence. Similar to diffusion tensor imaging, STI allows mapping of brain white matter fiber orientations and reconstruction of 3D white matter pathways using the principal eigenvectors of the susceptibility tensor. In contrast to diffusion anisotropy, the main determinant factor of the susceptibility anisotropy in brain white matter is myelin. Another unique feature of the susceptibility anisotropy of white matter is its sensitivity to gadolinium-based contrast agents. Mechanistically, MRI-observed susceptibility anisotropy is mainly attributed to the highly ordered lipid molecules in the myelin sheath. STI provides a consistent interpretation of the dependence of phase and susceptibility on orientation at multiple scales. This article reviews the key experimental findings and physical theories that led to the development of STI, its practical implementations, and its applications for brain research. Copyright © 2016 John Wiley & Sons, Ltd.

The clinical utility of QSM: disease diagnosis, medical management, and surgical planning

Quantitative susceptibility mapping (QSM) is an MR technique that depicts and quantifies magnetic susceptibility sources. Mapping iron, the dominant susceptibility source in the brain, has many important clinical applications. Herein, we review QSM applications in the diagnosis, medical management, and surgical treatment of disease. To assist in early disease diagnosis, QSM can identify elevated iron levels in the motor cortex of amyotrophic lateral sclerosis patients, in the substantia nigra of Parkinson's disease (PD) patients, in the globus pallidus, putamen, and caudate of Huntington's disease patients, and in the basal ganglia of Wilson's disease patients. Additionally, QSM can distinguish between hemorrhage and calcification, which could prove useful in tumor subclassification, and can measure microbleeds in traumatic brain injury patients. In guiding medical management, QSM can be used to monitor iron chelation therapy in PD patients, to monitor smoldering inflammation of multiple sclerosis (MS) lesions after the blood-brain barrier (BBB) seals, to monitor active inflammation of MS lesions before the BBB seals without using gadolinium, and to monitor hematoma volume in intracerebral hemorrhage. QSM can also guide neurosurgical treatment. Neurosurgeons require accurate depiction of the subthalamic nucleus, a tiny deep gray matter nucleus, prior to inserting deep brain stimulation electrodes into the brains of PD patients. QSM is arguably the best imaging tool for depiction of the subthalamic nucleus. Finally, we discuss future directions, including bone QSM, cardiac QSM, and using QSM to map cerebral metabolic rate of oxygen. Copyright © 2016 John Wiley & Sons, Ltd.

Overview of quantitative susceptibility mapping

Magnetic susceptibility describes the magnetizability of a material to an applied magnetic field and represents an important parameter in the field of MRI. With the recently introduced method of quantitative susceptibility mapping (QSM) and its conceptual extension to susceptibility tensor imaging (STI), the non-invasive assessment of this important physical quantity has become possible with MRI. Both methods solve the ill-posed inverse problem to determine the magnetic susceptibility from local magnetic fields. Whilst QSM allows the extraction of the spatial distribution of the bulk magnetic susceptibility from a single measurement, STI enables the quantification of magnetic susceptibility anisotropy, but requires multiple measurements with different orientations of the object relative to the main static magnetic field. In this review, we briefly recapitulate the fundamental theoretical foundation of QSM and STI, as well as computational strategies for the characterization of magnetic susceptibility with MRI phase data. In the second part, we provide an overview of current methodological and clinical applications of QSM with a focus on brain imaging. Copyright © 2016 John Wiley & Sons, Ltd.

Stability of R2* and quantitative susceptibility mapping of the brain tissue in a large scale multi-center study

Multi-center studies are advantageous for enrolling participants of varying pathological and demographical conditions, and especially in neurological studies. Hence stability of the obtained quantitative R2* and susceptibility in multicenter studies is a key issue for their widespread applications. In this work, the stabilities of simultaneously obtained R2* and susceptibility are investigated and compared across 10 sites that are equipped with the same scanner and receiver coil, the same post-processing process was used to achieve consistent experiment setup. Two healthy adult volunteers (one male and female) participated in this study. High intraclass correlation coefficient was obtained for both susceptibility (0.94) and R2* (0.96). The coefficients of variance for all measurements obtained were smaller than 0.1, the largest variations of measurements in all the chosen ROIs fall within ±20% from the median value. Higher level of stability was obtained in R2* as compared to susceptibility at 1 mm resolution (P < 0.05) and at 1.5 mm (P < 0.01).

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.

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-Based Neuroimaging: Standard Methods, Clinical Applications, and Future Directions

The evaluation of neuropathologies using MRI methods that leverage tissue susceptibility have become standard practice, especially to detect blood products or mineralization. Additionally, emerging MRI techniques have the ability to provide new information based on tissue susceptibility properties in a robust and quantitative manner. This paper discusses these advanced susceptibility imaging techniques and their clinical applications.

Exploring the origins of echo-time-dependent quantitative susceptibility mapping (QSM) measurements in healthy tissue and cerebral microbleeds

Quantitative susceptibility mapping (QSM) is increasingly used to measure variation in tissue composition both in the brain and in other areas of the body in a range of disease pathologies. Although QSM measurements were originally believed to be independent of the echo time (TE) used in the gradient-recalled echo (GRE) acquisition from which they are derived; recent literature (Sood et al., 2016) has shown that these measurements can be highly TE-dependent in a number of brain regions. In this work we systematically investigate possible causes of this effect through analysis of apparent frequency and QSM measurements derived from data acquired at multiple TEs in vivo in healthy brain regions and in cerebral microbleeds (CMBs); QSM data acquired in a gadolinium-doped phantom; and in QSM data derived from idealized simulated phase data. Apparent frequency measurements in the optic radiations (OR) and central corpus callosum (CC) were compared to those predicted by a 3-pool white matter model, however the model failed to fully explain contrasting frequency profiles measured in the OR and CC. Our results show that TE-dependent QSM measurements can be caused by a failure of phase unwrapping algorithms in and around strong susceptibility sources such as CMBs; however, in healthy brain regions this behavior appears to result from intrinsic non-linear phase evolution in the MR signal. From these results we conclude that care must be taken when deriving frequency and QSM measurements in strong susceptibility sources due to the inherent limitations in phase unwrapping; and that while signal compartmentalization due to tissue microstructure and content is a plausible cause of TE-dependent frequency and QSM measurements in healthy brain regions, better sampling of the MR signal and more complex models of tissue are needed to fully exploit this relationship.

Iron Mapping in Multiple Sclerosis

Increased iron deposition in cerebral deep gray matter has been considered a global marker for neurodegeneration in multiple sclerosis (MS); it scales with disease duration and severity. Iron accumulation in white matter and MS lesions might be more directly related to disease activity and has been discussed as a contributor to the inflammatory and neurodegenerative cascade. New insights into iron and MS are expected from MR imaging. We discuss findings from MR iron mapping proposed. Because of the confounding magnetic properties of myelin, iron mapping in white matter remains an unresolved issue.