MR parameter quantification with magnetization-prepared double echo steady-state (MP-DESS)

Primary Multiparametric Quantitative Brain MRI: State-of-the-Art Relaxometric and Proton Density Mapping Techniques

This review on brain multiparametric quantitative MRI (MP-qMRI) focuses on the primary subset of quantitative MRI (qMRI) parameters that represent the mobile ("free") and bound ("motion-restricted") proton pools. Such primary parameters are the proton densities, relaxation times, and magnetization transfer parameters. Diffusion qMRI is also included because of its wide implementation in complete clinical MP-qMRI application. MP-qMRI advances were reviewed over the past 2 decades, with substantial progress observed toward accelerating image acquisition and increasing mapping accuracy. Areas that need further investigation and refinement are identified as follows: (a) the biologic underpinnings of qMRI parameter values and their changes with age and/or disease and (b) the theoretical limitations implicitly built into most qMRI mapping algorithms that do not distinguish between the different spatial scales of voxels versus spin packets, the central physical object of the Bloch theory. With rapidly improving image processing techniques and continuous advances in computer hardware, MP-qMRI has the potential for implementation in a wide range of clinical applications. Currently, three emerging MP-qMRI applications are synthetic MRI, macrostructural qMRI, and microstructural tissue modeling.

Designing contrasts for rapid, simultaneous parameter quantification and flow visualization with quantitative transient-state imaging

Magnetic resonance imaging (MRI) has evolved into an outstandingly versatile diagnostic modality, as it has the ability to non-invasively produce detailed information on a tissue's structure and function. Complementary data is normally obtained in separate measurements, either as contrast-weighted images, which are fast and simple to acquire, or as quantitative parametric maps, which offer an absolute quantification of underlying biophysical effects, such as relaxation times or flow. Here, we demonstrate how to acquire and reconstruct data in a transient-state with a dual purpose: 1 - to generate contrast-weighted images that can be adjusted to emphasise clinically relevant image biomarkers; exemplified with signal modulation according to flow to obtain angiography information, and 2 - to simultaneously infer multiple quantitative parameters with a single, highly accelerated acquisition. This is achieved by introducing three novel elements: a model that accounts for flowing blood, a method for sequence design using smooth flip angle excitation patterns that incorporates both parameter encoding and signal contrast, and the reconstruction of temporally resolved contrast-weighted images. From these images we simultaneously obtain angiography projections and multiple quantitative maps. By doing so, we increase the amount of clinically relevant data without adding measurement time, creating new dimensions for biomarker exploration and adding value to MR examinations for patients and clinicians alike.

Accelerated 3D T2 mapping with dictionary-based matching for prostate imaging

PURPOSE

To develop a fast and accurate method for 3D T2 mapping of prostate cancer using undersampled acquisition and dictionary-based fitting.

METHODS

3D high-resolution T2 -weighted images (0.9 × 0.9 × 3 mm3 ) were obtained with a multishot T2 -prepared balanced steady-state free precession (T2 -prep-bSSFP) acquisition sequence using a 3D variable density undersampled Cartesian trajectory. Each T2 -weighted image was reconstructed using total variation regularized sensitivity encoding. A flexible simulation framework based on extended phase graphs generated a dictionary of magnetization signals, which was customized to the proposed sequence. The dictionary was matched to the acquired T2 -weighted images to retrieve quantitative T2 values, which were then compared to gold-standard spin echo acquisition values using monoexponential fitting. The proposed approach was validated in simulations and a T1 /T2 phantom, and feasibility was tested in 8 healthy subjects.

RESULTS

The simulation analysis showed that the proposed T2 mapping approach is robust to noise and typically observed T1 variations. T2 values obtained in the phantom with T2 prep-bSSFP and the acquisition-specific, dictionary-based matching were highly correlated with the gold-standard spin echo method (r = 0.99). Furthermore, no differences were observed with the accelerated acquisition compared to the fully sampled acquisition (r = 0.99). T2 values obtained in prostate peripheral zone, central gland, and muscle in healthy subjects (age, 26 ± 6 years) were 97 ± 14, 76 ± 7, and 36 ± 3 ms, respectively.

CONCLUSION

3D quantitative T2 mapping of the whole prostate can be achieved in 3 minutes.

Simultaneous multi-parametric mapping of total sodium concentration, T1, T2 and ADC at 7 T using a multi-contrast unbalanced SSFP

PURPOSE

Quantifying multiple NMR properties of sodium could be of benefit to assess changes in cellular viability in biological tissues. A proof of concept of Quantitative Imaging using Configuration States (QuICS) based on a SSFP sequence with multiple contrasts was implemented to extract simultaneously 3D maps of applied flip angle (FA), total sodium concentration, T₁, T₂, and Apparent Diffusion Coefficient (ADC).

METHODS

A 3D Cartesian Gradient Recalled Echo (GRE) sequence was used to acquire 11 non-balanced SSFP contrasts at a 6 × 6 × 6 mm³ isotropic resolution with carefully-chosen gradient spoiling area, RF amplitude and phase cycling, with TR/TE = 20/3.2 ms and 25 averages, leading to a total acquisition time of 1 h 18 min. A least-squares fit between the measured and the analytical complex signals was performed to extract quantitative maps from a mono-exponential model. Multiple sodium phantoms with different compositions were studied to validate the ability of the method to measure sodium NMR properties in various conditions.

RESULTS

Flip angle maps were retrieved. Relaxation times, ADC and sodium concentrations were estimated with controlled precision below 15%, and were in accordance with measurements from established methods and literature.

CONCLUSION

The results illustrate the ability to retrieve sodium NMR properties maps, which is a first step toward the estimation of FA, T₁, T₂, concentration and ADC of ²³Na for clinical research. With further optimization of the acquired QuICS contrasts, scan time could be reduced to be suitable with in vivo applications.

Optimizing MR Scan Design for Model-Based ${T}_{1}$ , ${T}_{2}$ Estimation From Steady-State Sequences

Rapid, reliable quantification of MR relaxation parameters T1 and T2 is desirable for many clinical applications. Steady-state sequences such as Spoiled Gradient-Recalled Echo (SPGR) and Dual-Echo Steady-State (DESS) are fast and well-suited for relaxometry because the signals they produce are quite sensitive to T1 and T2 variation. However, T1, T2 estimation with these sequences typically requires multiple scans with varied sets of acquisition parameters. This paper describes a systematic framework for selecting scan types (e.g., combinations of SPGR and DESS scans) and optimizing their respective parameters (e.g., flip angles and repetition times). The method is based on a Cramér-Rao Bound (CRB)-inspired min-max optimization that finds scan parameter combinations that robustly enable precise object parameter estimation. We apply this technique to optimize combinations of SPGR and DESS scans for T1, T2 relaxometry in white matter (WM) and grey matter (GM) regions of the human brain at 3T field strength. Phantom accuracy experiments show that SPGR/DESS scan combinations are in excellent agreement with reference measurements. Phantom precision experiments show that trends in T1,T2 pooled sample standard deviations reflect CRB-based predictions. In vivo experiments show that in WM and GM, T1 and T2 estimates from a pair of optimized DESS scans exhibit precision (but not necessarily accuracy) comparable to that of optimized combinations of SPGR and DESS scans. To our knowledge, T1 maps from DESS acquisitions alone are new. This example application illustrates that scan optimization may help reveal new parameter mapping techniques from combinations of established pulse sequences.

Diffusion-weighted DESS protocol optimization for simultaneous mapping of the mean diffusivity, proton density and relaxation times at 3 Tesla

PURPOSE

To design a general framework for the optimization of an MRI protocol based on the the diffusion-weighted dual-echo steady-state (DW-DESS) sequence, enabling quantitative and simultaneous mapping of proton density (PD), relaxation times T1 and T2 and diffusion coefficient D.

METHODS

A parameterization of the DW-DESS sequence minimizing the Cramér-Rao lower bound of each parameter estimate was proposed and tested in a phantom experiment. An extension of the protocol was implemented for brain imaging to return the rotationally invariant mean diffusivity (MD).

RESULTS

In an NiCl₂ -doped agar gel phantom wherein T1/T2=920/65 ms, the parameter estimation errors were below 3% for PD and T1 and below 7% for T2 and D while the measured signal-to-noise ratio always exceeded 20. In the human brain, the in vivo parametric maps obtained were overall in reasonable agreement with gold standard measurements, despite a broadening of the distributions due to physiological motion.

CONCLUSION

Within the optimization framework presented here, DW-DESS images can be quantitatively interpreted to yield four intrinsic parameters of the tissue. Currently, the method is limited by the sensitivity of the DW-DESS sequence in terms of physiological motion. Magn Reson Med 78:130-141, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Dual-pathway multi-echo sequence for simultaneous frequency and T2 mapping

PURPOSE

To present a dual-pathway multi-echo steady state sequence and reconstruction algorithm to capture T2, T2(∗) and field map information.

METHODS

Typically, pulse sequences based on spin echoes are needed for T2 mapping while gradient echoes are needed for field mapping, making it difficult to jointly acquire both types of information. A dual-pathway multi-echo pulse sequence is employed here to generate T2 and field maps from the same acquired data. The approach might be used, for example, to obtain both thermometry and tissue damage information during thermal therapies, or susceptibility and T2 information from a same head scan, or to generate bonus T2 maps during a knee scan.

RESULTS

Quantitative T2, T2(∗) and field maps were generated in gel phantoms, ex vivo bovine muscle, and twelve volunteers. T2 results were validated against a spin-echo reference standard: A linear regression based on ROI analysis in phantoms provided close agreement (slope/R(2)=0.99/0.998). A pixel-wise in vivo Bland-Altman analysis of R2=1/T2 showed a bias of 0.034 Hz (about 0.3%), as averaged over four volunteers. Ex vivo results, with and without motion, suggested that tissue damage detection based on T2 rather than temperature-dose measurements might prove more robust to motion.

CONCLUSION

T2, T2(∗) and field maps were obtained simultaneously, from the same datasets, in thermometry, susceptibility-weighted imaging and knee-imaging contexts.

MR fingerprinting using the quick echo splitting NMR imaging technique

PURPOSE

The purpose of the study is to develop a quantitative method for the relaxation properties with a reduced radio frequency (RF) power deposition by combining magnetic resonance fingerprinting (MRF) technique with quick echo splitting NMR imaging technique (QUEST).

METHODS

A QUEST-based MRF sequence was implemented to acquire high-order echoes by increasing the gaps between RF pulses. Bloch simulations were used to calculate a dictionary containing the range of physically plausible signal evolutions using a range of T₁ and T₂ values based on the pulse sequence. MRF-QUEST was evaluated by comparing to the results of spin-echo methods. The specific absorption rate (SAR) of MRF-QUEST was compared with the clinically available methods.

RESULTS

MRF-QUEST quantifies the relaxation properties with good accuracy at the estimated head SAR of 0.03 W/kg. T₁ and T₂ values estimated by MRF-QUEST are in good agreement with the traditional methods.

CONCLUSIONS

The combination of the MRF and the QUEST provides an accurate quantification of T₁ and T₂ simultaneously with reduced RF power deposition. The resulting lower SAR may provide a new acquisition strategy for MRF when RF energy deposition is problematic. Magn Reson Med 77:979-988, 2017. © 2016 International Society for Magnetic Resonance in Medicine.