Quantitative susceptibility mapping of small objects using volume constraints

Brain iron imaging markers in the presence of white matter hyperintensities

PURPOSE

To investigate the relationship between pathological brain iron deposition and white matter hyperintensities (WMHs) in cerebral small vessel disease (CSVD), via Monte Carlo simulations of magnetic susceptibility imaging and the development of a novel imaging marker called the Expected Iron Coefficient (EIC).

METHODS

A synthetic pathological model of a different number of impenetrable spheres at random locations was employed to represent pathological iron deposition. The diffusion process was simulated with a Monte Carlo method with adjustable parameters to manipulate sphere size, distribution, and extracellular properties. Quantitative susceptibility mapping (QSM) was performed in a clinical dataset to study CSVD to derive and evaluate QSM, R2*, the iron microenvironment coefficient (IMC), and the EIC in the presence of WMHs.

RESULTS

The simulations show that QSM signals increase in the presence of increased tissue iron, confirming that the EIC increases with pathology. Clinical results demonstrate that while QSM, R2*, and the IMC do not show significant differences in brain iron, the EIC does in the context of CSVD.

CONCLUSION

The EIC is more sensitive to subtle changes in brain iron deposition caused by pathology, even when QSM, R2*, and the IMC fail.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Longitudinal ultra-high field MRI of brain lesions in neuromyelitis optica spectrum disorders

BACKGROUND

In neuromyelitis optica spectrum disorder (NMOSD), clinical disability in NMOSD patients is relapse-related and progressive phase is rare. This observation raises the question whether there is any radiographic disease activity. The aim of present study was to determine the longitudinal changes in cerebral lesion number, lesion size, lesion-to-venule relationship, and morphological patterns of lesions in NMOSD using multiparametric 7T MR imaging. We also aimed to assess brain volume changes in NMOSD.

METHODS

A cohort of 22 patients with NMOSD underwent high-resolution 3D-susceptibility weighted imaging (SWI) and 2D-gradient-echo (GRE-T2*) weighted imaging on 7T MRI of brain at baseline and after ~2.8 years of follow-up. Morphologic imaging characteristics, and signal intensity patterns of lesions were recorded at both time points. Lesions were classified as "iron-laden" if they demonstrated hypointense signal on GRE-T2* images and/or SWI as well as hyperintense signal on quantitative susceptibility mapping (QSM). Lesions were considered "non-iron-laden" if they were hyperintense on GRE-T2*/SWI and isointense or hyperintense on QSM. Additionally, fractional brain parenchymal volume (fBPV) was computed at both time points.

RESULTS

A total of 169 lesions were observed at baseline. At follow-up, 6 new lesions were found in 5 patients. In one patient, a single lesion could not be detected on the follow-up scan. No appreciable change in lesion size and vessel-lesion relationship was observed at follow up. All lesions demonstrated hyperintense signal intensity on GRE-T2* weighted images and isointense signal on QSM at both time points. Therefore, these lesions were considered as non-associated with iron pathology. Additionally, no significant change in brain volume was observed: fBPV 0.78 ± 0.06 at baseline vs. 0.77 ± 0.05 at follow up, p>0.05.

CONCLUSION

Cerebral lesions in NMOSD patients remain 'inert' and do not show any substantial variations in morphological characteristics during a 2-3-year follow-up period.

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.

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.

Susceptibility weighted imaging and quantitative susceptibility mapping of the cerebral vasculature using ferumoxytol

PURPOSE

To demonstrate the potential of imaging cerebral arteries and veins with ferumoxytol using susceptibility weighted imaging (SWI) and quantitative susceptibility mapping (QSM).

MATERIALS AND METHODS

The relationships between ferumoxytol concentration and the apparent susceptibility at 1.5T, 3T, and 7T were determined using phantom data; the ability of visualizing subvoxel vessels was evaluated using simulations; and the feasibility of using ferumoxytol to enhance the visibility of small vessels was confirmed in three healthy volunteers at 7T(with doses 1 mg/kg to 4 mg/kg). The visualization of the lenticulostriate arteries and the medullary veins was assessed by two raters and the contrast-to-noise ratios (CNRs) of these vessels were measured.

RESULTS

The relationship between ferumoxytol concentration and susceptibility was linear with a slope 13.3 ± 0.2 ppm·mg^-1 ·mL at 7T. Simulations showed that SWI data with an increased dose of ferumoxytol, higher echo time (TE), and higher imaging resolution improved the detection of smaller vessels. With 4 mg/kg ferumoxytol, voxel aspect ratio = 1:8, TE = 10 ms, the diameter of the smallest detectable artery was approximately 50μm. The rating score for arteries was improved from 1.5 ± 0.5 (precontrast) to 3.0 ± 0.0 (post-4 mg/kg) in the in vivo data and the apparent susceptibilities of the arteries (0.65 ± 0.02 ppm at 4 mg/kg) agreed well with the expected susceptibility (0.71 ± 0.05 ppm).

CONCLUSION

The CNR for cerebral vessels with ferumoxytol can be enhanced using SWI, and the apparent susceptibilities of the arteries can be reliably quantified using QSM. This approach improves the imaging of the entire vascular system outside the capillaries and may be valuable for a variety of neurodegenerative diseases which involve the microvasculature.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;47:621-633.

Determination of detection sensitivity for cerebral microbleeds using susceptibility-weighted imaging

Cerebral microbleeds (CMBs) are small brain hemorrhages caused by the break down or structural abnormalities of small vessels of the brain. Owing to the paramagnetic properties of blood degradation products, CMBs can be detected in vivo using susceptibility-weighted imaging (SWI). SWI can be used not only to detect iron changes and CMBs, but also to differentiate them from calcifications, both of which may be important MR-based biomarkers for neurodegenerative diseases. Moreover, SWI can be used to quantify the iron in CMBs. SWI and gradient echo (GE) imaging are the two most common methods for the detection of iron deposition and CMBs. This study provides a comprehensive analysis of the number of voxels detected in the presence of a CMB on GE magnitude, phase and SWI composite images as a function of resolution, signal-to-noise ratio (SNR), TE, field strength and susceptibility using in silico experiments. Susceptibility maps were used to quantify the bias in the effective susceptibility value and to determine the optimal TE for CMB quantification. We observed a non-linear trend with susceptibility for CMB detection from the magnitude images, but a linear trend with susceptibility for CMB detection from the phase and SWI composite images. The optimal TE values for CMB quantification were found to be 3 ms at 7 T, 7 ms at 3 T and 14 ms at 1.5 T for a CMB of one voxel in diameter with an SNR of 20: 1. The simulations of signal loss and detectability were used to generate theoretical formulae for predictions. Copyright © 2016 John Wiley & Sons, Ltd.

Evaluating the Role of Reduced Oxygen Saturation and Vascular Damage in Traumatic Brain Injury Using Magnetic Resonance Perfusion-Weighted Imaging and Susceptibility-Weighted Imaging and Mapping

The cerebral vasculature, along with neurons and axons, is vulnerable to biomechanical insult during traumatic brain injury (TBI). Trauma-induced vascular injury is still an underinvestigated area in TBI research. Cerebral blood flow and metabolism could be important future treatment targets in neural critical care. Magnetic resonance imaging offers a number of key methods to probe vascular injury and its relationship with traumatic hemorrhage, perfusion deficits, venous blood oxygen saturation changes, and resultant tissue damage. They make it possible to image the hemodynamics of the brain, monitor regional damage, and potentially show changes induced in the brain's function not only acutely but also longitudinally following treatment. These methods have recently been used to show that even mild TBI (mTBI) subjects can have vascular abnormalities, and thus they provide a major step forward in better diagnosing mTBI patients.

Iron and Non-Iron-Related Characteristics of Multiple Sclerosis and Neuromyelitis Optica Lesions at 7T MRI

BACKGROUND AND PURPOSE

Characterization of iron deposition associated with demyelinating lesions of multiple sclerosis and neuromyelitis optica has not been well studied. Our aim was to investigate the potential of ultra-high-field MR imaging to distinguish MS from neuromyelitis optica and to characterize tissue injury associated with iron pathology within lesions.

MATERIALS AND METHODS

Twenty-one patients with MS and 21 patients with neuromyelitis optica underwent 7T high-resolution 2D-gradient-echo-T2* and 3D-susceptibility-weighted imaging. An in-house-developed algorithm was used to reconstruct quantitative susceptibility mapping from SWI. Lesions were classified as "iron-laden" if they demonstrated hypointensity on gradient-echo-T2*-weighted images and/or SWI and hyperintensity on quantitative susceptibility mapping. Lesions were considered "non-iron-laden" if they were hyperintense on gradient-echo-T2* and isointense or hyperintense on quantitative susceptibility mapping.

RESULTS

Of 21 patients with MS, 19 (90.5%) demonstrated at least 1 quantitative susceptibility mapping-hyperintense lesion, and 11/21 (52.4%) had iron-laden lesions. No quantitative susceptibility mapping-hyperintense or iron-laden lesions were observed in any patients with neuromyelitis optica. Iron-laden and non-iron-laden lesions could each be further characterized into 2 distinct patterns based on lesion signal and morphology on gradient-echo-T2*/SWI and quantitative susceptibility mapping. In MS, most lesions (n = 262, 75.9% of all lesions) were hyperintense on gradient-echo T2* and isointense on quantitative susceptibility mapping (pattern A), while a small minority (n = 26, 7.5% of all lesions) were hyperintense on both gradient-echo-T2* and quantitative susceptibility mapping (pattern B). Iron-laden lesions (n = 57, 16.5% of all lesions) were further classified as nodular (n = 22, 6.4%, pattern C) or ringlike (n = 35, 10.1%, pattern D).

CONCLUSIONS

Ultra-high-field MR imaging may be useful in distinguishing MS from neuromyelitis optica. Different patterns related to iron and noniron pathology may provide in vivo insight into the pathophysiology of lesions in MS.

Susceptibility Weighted Imaging and Mapping of Micro-Hemorrhages and Major Deep Veins after Traumatic Brain Injury

Micro-hemorrhages are a common result of traumatic brain injury (TBI), which can be quantified with susceptibility weighted imaging and mapping (SWIM), a quantitative susceptibility mapping approach. A total of 23 TBI patients (five women, 18 men; median age, 41.25 years old; range, 21.69-67.75 years) with an average Glasgow Coma Scale score of 7 (range, 3-15) at admission were recruited at mean 149 d (range, 57-366) after injury. Susceptibility-weighted imaging data were collected and post-processed to create SWIM images. The susceptibility value of small hemorrhages (diameter ≤10 mm) and major deep veins (right septal, left septal, central septal, right thalamostriate, left thalamostriate, internal cerebral, right basal vein of Rosenthal, left basal vein of Rosenthal, and pial veins) were evaluated. Different susceptibility thresholds were tested to determine SWIM's sensitivity and specificity for differentiating hemorrhages from the veins. A total of 253 deep veins and 173 small hemorrhages were identified and evaluated. The mean susceptibility of hemorrhages was 435±206 parts per billion (ppb) and the mean susceptibility of deep veins was 108±56 ppb. Hemorrhages showed a significantly higher susceptibility than all deep veins (p<0.001). With different thresholds (250, 227 and 200 ppb), the specificity was 97%, 95%, and 92%, and the sensitivity was 84%, 90%, and 92%, respectively. These results show that SWIM could be used to differentiate hemorrhages from veins in TBI patients in a semi-automated manner with reasonable sensitivity and specificity. A larger cohort will be needed to validate these findings.

An improved method for susceptibility and radius quantification of cylindrical objects from MRI

A new method is developed to measure the magnetic susceptibilities and radii of small cylinder-like objects at arbitrary orientations accurately. This method for most biological substances only requires a standard gradient echo sequence with one or two echo times, depending on the orientation of an object relative to the main magnetic field. For objects oriented at the magic angle, however, this method is not applicable. As a byproduct of this method, the cross-sectional area as well as signals inside and outside the object can be determined. The uncertainty of each measurement is estimated from the error propagation method. Partial volume, dephasing, and phase aliasing effects are naturally included in the equations of this method. A number of simulations, phantom, and pilot in-vivo human studies are carried out to validate the theory. When the maximal phase value at the boundary of a given cylindrical object is larger than 3 radians, and the phase inside the object is more than 1 radian, the susceptibility can be accurately quantified within 15%. The radius of the object can be determined to subpixel accuracy. This is the case when the signal-to-noise ratio inside the object is about 6:1 or higher and the radius of the object is about one pixel or larger. These conditions are realistic when considering medullary and pial veins for example.

Magnetic moment quantifications of small spherical objects in MRI

PURPOSE

The purpose of this work is to develop a method for accurately quantifying effective magnetic moments of spherical-like small objects from magnetic resonance imaging (MRI). A standard 3D gradient echo sequence with only one echo time is intended for our approach to measure the effective magnetic moment of a given object of interest.

METHODS

Our method sums over complex MR signals around the object and equates those sums to equations derived from the magnetostatic theory. With those equations, our method is able to determine the center of the object with subpixel precision. By rewriting those equations, the effective magnetic moment of the object becomes the only unknown to be solved. Each quantified effective magnetic moment has an uncertainty that is derived from the error propagation method. If the volume of the object can be measured from spin echo images, the susceptibility difference between the object and its surrounding can be further quantified from the effective magnetic moment. Numerical simulations, a variety of glass beads in phantom studies with different MR imaging parameters from a 1.5T machine, and measurements from a SQUID (superconducting quantum interference device) based magnetometer have been conducted to test the robustness of our method.

RESULTS

Quantified effective magnetic moments and susceptibility differences from different imaging parameters and methods all agree with each other within two standard deviations of estimated uncertainties.

CONCLUSION

An MRI method is developed to accurately quantify the effective magnetic moment of a given small object of interest. Most results are accurate within 10% of true values, and roughly half of the total results are accurate within 5% of true values using very reasonable imaging parameters. Our method is minimally affected by the partial volume, dephasing, and phase aliasing effects. Our next goal is to apply this method to in vivo studies.

Quantitative susceptibility mapping: current status and future directions

Quantitative susceptibility mapping (QSM) is a new technique for quantifying magnetic susceptibility. It has already found various applications in quantifying in vivo iron content, calcifications and changes in venous oxygen saturation. The accuracy of susceptibility mapping is dependent on several factors. In this review, we evaluate the entire process of QSM from data acquisition to individual data processing steps. We also show preliminary results of several new concepts introduced in this review in an attempt to improve the quality and accuracy for certain steps. The uncertainties in estimating susceptibility differences using susceptibility maps, phase images, and T2* maps are analyzed and compared. Finally, example clinical applications are presented. We conclude that QSM holds great promise in quantifying iron and becoming a standard clinical tool.

Measuring venous blood oxygenation in fetal brain using susceptibility-weighted imaging

PURPOSE

To evaluate fetal cerebral venous blood oxygenation, Yv, using principles of MR susceptometry.

MATERIALS AND METHODS

A cohort of 19 pregnant subjects, with a mean gestational age of 31.6 ± 4.7 weeks were imaged using a modified susceptibility-weighted imaging (SWI) sequence. Data quality was first assessed for feasibility of oxygen saturation measurement, and data from five subjects (mean ± std gestational age of 33.7 ± 3.6 weeks) were then chosen for further quantitative analysis. SWI phase in the superior sagittal sinus was used to evaluate oxygen saturation using the principles of MR susceptometry. Systematic error in the measured Y(v) values was studied through simulations.

RESULTS

Simulations showed that the systematic error in Yv depended upon the assumed angle of the vessel, θ, relative to the main magnetic field and the error in that vessel angle δθ. For the typical vessel angle of θ = 30° encountered in the fetal data analyzed, a δθ as large as ±20° led to an absolute error, δYv, of less than 11%. The measured mean oxygen saturation across the five fetuses was 66% ± 9.4%. This average cerebral venous blood oxygenation value is in close agreement with values in the published literature.

CONCLUSION

We have reported the first in vivo measurement of human fetal cerebral venous oxygen saturation using MRI.