Quantitative Susceptibility Mapping: Contrast Mechanisms and Clinical Applications.

Quantitative susceptibility mapping (QSM) is a recently developed MRI technique for quantifying the spatial distribution of magnetic susceptibility within biological tissues. It first uses the frequency shift in the MRI signal to map the magnetic field profile within the tissue. The resulting field map is then used to determine the spatial distribution of the underlying magnetic susceptibility by solving an inverse problem. The solution is achieved by deconvolving the field map with a dipole field, under the assumption that the magnetic field is a result of the superposition of the dipole fields generated by all voxels and that each voxel has its unique magnetic susceptibility. QSM provides improved contrast to noise ratio for certain tissues and structures compared to its magnitude counterpart. More importantly, magnetic susceptibility is a direct reflection of the molecular composition and cellular architecture of the tissue. Consequently, by quantifying magnetic susceptibility, QSM is becoming a quantitative imaging approach for characterizing normal and pathological tissue properties. This article reviews the mechanism generating susceptibility contrast within tissues and some associated applications

Advances in Quantitative MRI: Acquisition, Estimation, and Application

Quantitative magnetic resonance imaging (QMRI) produces images of potential MR biomarkers: measurable tissue properties related to physiological processes that characterize the onset and progression of specific disorders. Though QMRI has potential to be more diagnostic than conventional qualitative MRI, QMRI poses challenges beyond those of conventional MRI that limit its feasibility for routine clinical use. This thesis first seeks to address two of those challenges. It then applies these solutions to develop a new method for myelin water imaging, a challenging application that may be specifically indicative of certain white matter (WM) disorders. One challenge that presently precludes widespread clinical adoption of QMRI involves long scan durations: to disentangle potential biomarkers from nuisance MR contrast mechanisms, QMRI typically requires more data than conventional MRI and thus longer scans. Even allowing for long scans, it has previously been unclear how to systematically tune the "knobs" of MR acquisitions to reliably enable precise biomarker estimation. Chapter 4 formalizes these challenges as a min-max optimal acquisition design problem and solves this problem to design three fast steady-state (SS) acquisitions for precise T1/T2 estimation, a popular QMRI application. The resulting optimized acquisition designs illustrate that acquisition design can enable new biomarker estimation techniques from established MR pulse sequences, a fact that subsequent chapters exploit. Another QMRI challenge involves the typically nonlinear dependence of MR signal models on the underlying biomarkers: these nonlinearities cause conventional likelihood-based estimators to either scale very poorly with the number of unknowns or risk producing suboptimal estimates due to spurious local minima. Chapter 5 instead introduces a fast, general method for dictionary-free QMRI parameter estimation via regression with kernels (PERK). PERK first uses prior distributions and the nonlinear MR signal model to simulate many parameter-measurement pairs. Inspired by machine learning, PERK then takes these pairs as labeled training points and learns from them a nonlinear regression function using kernel functions and convex optimization. Chapter 5 demonstrates PERK for T1/T2 estimation using one of the acquisitions optimized in Chapter 4. Simulations as well as single-slice phantom and in vivo experiments demonstrated that PERK and two well-suited maximum-likelihood (ML) estimators produce comparable T1/T2 estimates, but PERK is consistently at least 140x faster. Similar comparisons to an ML estimator in a more challenging problem (Chapter 6) suggest that this 140x acceleration factor will increase by several orders of magnitude for full-volume QMRI estimation problems involving more latent parameters per voxel. Chapter 6 applies ideas developed in previous chapters to design a new fast method for imaging myelin water content, a potential biomarker for healthy myelin. It first develops a two-compartment dual-echo steady-state (DESS) signal model and then uses a Bayesian variation of acquisition design (Chapter 4) to optimize a new DESS acquisition for precise myelin water imaging. The precision-optimized acquisition is as fast as conventional SS myelin water imaging acquisitions, but enables 2-3x better expected coefficients of variation in fast-relaxing fraction estimates. Simulations demonstrate that PERK (Chapter 5) and ML fast-relaxing fraction estimates from the proposed DESS acquisition exhibit comparable root mean-squared errors, but PERK is more than 500x faster. In vivo experiments are to our knowledge the first to demonstrate lateral WM myelin water content estimates from a fast (3m15s) SS acquisition that are similar to conventional estimates from a slower (32m4s) MESE acquisition.PHDElectrical and Computer EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/147486/1/gnataraj_1.pd

Quantitative Susceptibility Imaging of Tissue Microstructure Using Ultra-High Field MRI

This thesis has used ultra-high field (UHF) magnetic resonance imaging (MRI) to investigate the fundamental relationships between tissue microstructure and such susceptibility-based contrast parameters as the apparent transverse relaxation rate (R2*), the local Larmor frequency shift (LFS) and quantitative volume magnetic susceptibility (QS). The interaction of magnetic fields with biological tissues results in shifts in the LFS which can be used to distinguish underlying cellular architecture. The LFS is also linked to the relaxation properties of tissues in a gradient echo MRI sequence. Equally relevant, histological analysis has identified iron and myelin as two major sources of the LFS. As a result, computation of LFS and the associated volume magnetic susceptibility from MRI phase data may serve as a significant method for in vivo monitoring of changes in iron and myelin associated with normal, healthy aging, as well as neurological disease processes. In this research, the cellular level underpinnings of the R2* and LFS signals were examined in a model rat brain system using 9.4 T MRI. The study was carried out using biophysical modeling and correlation with quantitative histology. For the first time, multiple biophysical modeling schemes were compared in both gray and white matter of excised rat brain tissue. Suprisingly, R2* dependence on tissue orientation has not been fully understood. Accordingly, scaling relations were derived for calculating the reversible, mesoscopic magnetic field component, R2\u27, of the apparent transverse relaxation rate from the orientation dependence in gray and white matter. Our results demonstrate that the orientation dependence of R2* and LFS in both white and cortical gray matter has a sinusoidal dependence on tissue orientation and a linear dependence on the volume fraction of myelin in the tissue. A susceptibility processing pipeline was also developed and applied to the calculation of phase-combined LFS and QS maps. The processing pipeline was subsequently used to monitor myelin and iron changes in multiple sclerosis (MS) patients compared to healthy, age and gender-matched controls. With the use of QS and R2* mapping, evidence of statistically significant increases in iron deposition in sub-cortical gray matter, as well as myelin degeneration along the white matter skeleton, were identified in MS patients. The magnetic susceptibility-based MRI methods were then employed as potential clinical biomarkers for disease severity monitoring of MS. It was demonstrated that the combined use of R2* and QS, obtained from multi-echo gradient echo MRI, could serve as an improved metric for monitoring both gray and white matter changes in early MS

Advances in noninvasive myelin imaging

Myelin is important for the normal development and healthy function of the nervous system. Recent developments in MRI acquisition and tissue modeling aim to provide a better characterization and more specific markers for myelin. This allows for specific monitoring of myelination longitudinally and noninvasively in the healthy brain as well as assessment of treatment and intervention efficacy. Here, we offer a nontechnical review of MRI techniques developed to specifically monitor myelin such as magnetization transfer (MT) and myelin water imaging (MWI). We further summarize recent studies that employ these methods to measure myelin in relation to development and aging, learning and experience, and neuropathology and psychiatric disorders

Fat fraction mapping using bSSFP Signal Profile Asymmetries for Robust multi-Compartment Quantification (SPARCQ)

Purpose: To develop a novel quantitative method for detection of different tissue compartments based on bSSFP signal profile asymmetries (SPARCQ) and to provide a validation and proof-of-concept for voxel-wise water-fat separation and fat fraction mapping. Methods: The SPARCQ framework uses phase-cycled bSSFP acquisitions to obtain bSSFP signal profiles. For each voxel, the profile is decomposed into a weighted sum of simulated profiles with specific off-resonance and relaxation time ratios. From the obtained set of weights, voxel-wise estimations of the fractions of the different components and their equilibrium magnetization are extracted. For the entire image volume, component-specific quantitative maps as well as banding-artifact-free images are generated. A SPARCQ proof-of-concept was provided for water-fat separation and fat fraction mapping. Noise robustness was assessed using simulations. A dedicated water-fat phantom was used to validate fat fractions estimated with SPARCQ against gold-standard 1H MRS. Quantitative maps were obtained in knees of six healthy volunteers, and SPARCQ repeatability was evaluated in scan rescan experiments. Results: Simulations showed that fat fraction estimations are accurate and robust for signal-to-noise ratios above 20. Phantom experiments showed good agreement between SPARCQ and gold-standard (GS) fat fractions (fF(SPARCQ) = 1.02*fF(GS) + 0.00235). In volunteers, quantitative maps and banding-artifact-free water-fat-separated images obtained with SPARCQ demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ fat fraction was 0.0512. Conclusion: The SPARCQ framework was proposed as a novel quantitative mapping technique for detecting different tissue compartments, and its potential was demonstrated for quantitative water-fat separation.Comment: 20 pages, 7 figures, submitted to Magnetic Resonance in Medicin

White matter changes measured by multi-component MR Fingerprinting in multiple sclerosis

T2-hyperintense lesions are the key imaging marker of multiple sclerosis (MS). Previous studies have shown that the white matter surrounding such lesions is often also affected by MS. Our aim was to develop a new method to visualize and quantify the extent of white matter tissue changes in MS based on relaxometry properties. We applied a fast, multi-parametric quantitative MRI approach and used a multi-component MR Fingerprinting (MC-MRF) analysis. We assessed the differences in the MRF component representing prolongedrelaxation time between patients with MS and controls and studied the relation between this component's volume and structural white matter damage identified on FLAIR MRI scans in patients with MS. A total of 48 MS patients at two different sites and 12 healthy controls were scanned with FLAIR and MRF-EPI MRI scans. MRF scans were analyzed with a joint-sparsity multi-component analysis to obtain magnetization fraction maps of different components, representing tissues such as myelin water, white matter, gray matter and cerebrospinal fluid. In the MS patients, an additional component was identified with increased transverse relaxation times compared to the white matter, likely representing changes in free water content. Patients with MS had a higher volume of the long- component in the white matter of the brain compared to healthy controls (B (95%-CI) = 0.004 (0.0006–0.008), p = 0.02). Furthermore, this MRF component had a moderate correlation (correlation coefficient R 0.47) with visible structural white matter changes on the FLAIR scans. Also, the component was found to be more extensive compared to structural white matter changes in 73% of MS patients. In conclusion, our MRF acquisition and analysis captured white matter tissue changes in MS patients compared to controls. In patients these tissue changes were more extensive compared to visually detectable white matter changes on FLAIR scans. Our method provides a novel way to quantify the extent of white matter changes in MS patients, which is underestimated using only conventional clinical MRI scans.</p

In vivo myelin water quantification using diffusion–relaxation correlation MRI: A comparison of 1D and 2D methods

Multidimensional Magnetic Resonance Imaging (MRI) is a versatile tool for microstructure mapping. We use a diffusion weighted inversion recovery spin echo (DW-IR-SE) sequence with spiral readouts at ultra-strong gradients to acquire a rich diffusion–relaxation data set with sensitivity to myelin water. We reconstruct 1D and 2D spectra with a two-step convex optimization approach and investigate a variety of multidimensional MRI methods, including 1D multi-component relaxometry, 1D multi-component diffusometry, 2D relaxation correlation imaging, and 2D diffusion-relaxation correlation spectroscopic imaging (DR-CSI), in terms of their potential to quantify tissue microstructure, including the myelin water fraction (MWF). We observe a distinct spectral peak that we attribute to myelin water in multi-component T1 relaxometry, T1-T2 correlation, T1-D correlation, and T2-D correlation imaging. Due to lower achievable echo times compared to diffusometry, MWF maps from relaxometry have higher quality. Whilst 1D multi-component T1 data allows much faster myelin mapping, 2D approaches could offer unique insights into tissue microstructure and especially myelin diffusion