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

Synthetic MR: Clinical applications in neuroradiology

PURPOSE

Synthetic MR is a quantitative MRI method that measures tissue relaxation times and generates multiple contrast-weighted images using suitable algorithms. The present article principally discusses the multiple dynamic multiple echo (MDME) technique of synthetic MR and briefly describes other quantitative MR sequences.

METHODS

Using illustrative cases, various applications of the MDME sequence in neuroradiology are explained. The MDME sequence allows rapid quantification of tissue relaxation times in a scan duration of 5-7 minutes for full brain coverage. It also has the additional advantages of myelin quantification and automatic segmentation of brain volumes.

RESULTS

Applications including reducing scan time, improved detection of demyelinating plaques in Multiple Sclerosis (MS), objective assessment and follow-up for brain atrophy in neurodegenerative MS and dementia cases, and applications in stroke imaging and neuro-oncology are discussed. Uses in the pediatric population, including assessment of brain development and progression of myelination in children, evaluation of white matter disorders, and evaluation of pediatric and adult epilepsy, are elaborated. Quantitative evaluation by synthetic MR is discussed, which allows homogenization and objectification of the radiology data and can serve as a valuable source for artificial intelligence and future multicentre studies. A brief discussion on the technique, other quantitative MR methods, and limitations of the MDME sequence is also presented.

CONCLUSION

The article intends to provide an explicit and comprehensive review of the applications of synthetic MR in neuroradiology, exploring its potential as a routine sequence in daily neuroimaging practice.

Improved gradient echo magnitude- and phase-based mapping of T 2 $$ {\mathrm{T}}_2 $$ using multiple RF spoiling increments at 3T and 7T

PURPOSE

The transverse relaxation time T  2 $$ {}_2 $$ holds significant relevance in clinical applications and research studies. Conventional T  2 $$ {}_2 $$ mapping approaches rely on spin-echo sequences, which require lengthy acquisition times and involve high radiofrequency (RF) power deposition. An alternative gradient echo (GRE) phase-based T  2 $$ {}_2 $$ mapping method, utilizing steady-state acquisitions at one small RF spoil phase increment, was recently demonstrated. Here, a modified magnitude- and phase-based T  2 $$ {}_2 $$ mapping approach is proposed, which improvesT 2 $$ {\mathrm{T}}_2 $$ estimations by simultaneous fitting ofT 1 $$ {\mathrm{T}}_1 $$ and signal amplitude (A ∝ P D $$ A\propto PD $$ ) at three or more RF spoiling phase increments, instead of assuming a fixedT 1 $$ {\mathrm{T}}_1 $$ value.

METHODS

The feasibility of the magnitude-phase-based method was assessed by simulations, in phantom and in vivo experiments using skipped-CAIPI three-dimensional-echo-planar imaging (3D-EPI) for rapid GRE imaging.T 2 $$ {\mathrm{T}}_2 $$ ,T 1 $$ {\mathrm{T}}_1 $$ and PD estimations obtained by our method were compared toT 2 $$ {\mathrm{T}}_2 $$ of the phase-based method andT 1 $$ {\mathrm{T}}_1 $$ and PD of spoiled GRE-based multi-parameter mapping using a multi-echo version of the same sequence.

RESULTS

The agreement of the proposedT 2 $$ {\mathrm{T}}_2 $$ with ground truth and referenceT 2 $$ {\mathrm{T}}_2 $$ values was higher than that of phase-basedT 2 $$ {\mathrm{T}}_2 $$ in simulations and in phantom data. While phase-basedT 2 $$ {\mathrm{T}}_2 $$ overestimation increases with actualT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ , the proposed method is accurate over a large range of physiologically meaningfulT 2 $$ {\mathrm{T}}_2 $$ andT 1 $$ {\mathrm{T}}_1 $$ values. At the same time, precision is improved. In vivo results were in line with these observations.

CONCLUSION

Accurate magnitude-phase-based T  2 $$ {}_2 $$ mapping is feasible in less than 5 min scan time for 1 mm nominal isotropic whole-head coverage at 3T and 7T.

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.