Magnetic Resonance Fingerprinting with short relaxation intervals

Accelerated liver water T1 mapping using single-shot continuous inversion-recovery spiral imaging

PURPOSE

Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency.

METHODS

The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients.

RESULTS

The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility.

CONCLUSION

A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.

Nonuniform sliding-window reconstruction for accelerated dual contrast agent quantification with MR fingerprinting

OBJECTIVE

MR fingerprinting (MRF) can enable preclinical studies of cell tracking by quantifying multiple contrast agents simultaneously, but faster scan times are required for in vivo applications. Sliding window (SW)-MRF is one option for accelerating MRF, but standard implementations are not sufficient to preserve the accuracy of T2*, which is critical for tracking iron-labelled cells in vivo.

PURPOSE

To develop a SW approach to MRF which preserves the T2* accuracy required for accelerated concentration mapping of iron-labelled cells on single-channel preclinical systems.

METHODS

A nonuniform SW was applied to the MRF sequence and dictionary. Segments of the sequence most sensitive to T2* were subject to a shorter window length, preserving the T2* sensitivity. Phantoms containing iron-labelled CD8+ T cells and gadolinium were used to compare 24× undersampled uniform and nonuniform SW-MRF parameter maps. Dual concentration maps were generated for both uniform and nonuniform MRF and compared.

RESULTS

Lin's concordance correlation coefficient, compared to gold standard parameter values, was much greater for nonuniform SW-MRF than for uniform SW-MRF. A Wilcoxon signed-rank test showed no significant difference between nonuniform SW-MRF and gold standards. Nonuniform SW-MRF outperformed the uniform SW-MRF concentration maps for all parameters, providing a balance between T2* sensitivity of short window lengths, and SNR of longer window lengths.

CONCLUSIONS

Nonuniform SW-MRF improves the accuracy of matching compared to uniform SW-MRF, allowing higher accelerated concentration mapping for preclinical systems.

A three-dimensional Magnetic Resonance Spin Tomography in Time-domain protocol for high-resolution multiparametric quantitative magnetic resonance imaging

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) is a multiparametric quantitative MR framework, which allows for simultaneously acquiring quantitative tissue parameters such as T1, T2, and proton density from one single short scan. A typical two-dimensional (2D) MR-STAT acquisition uses a gradient-spoiled, gradient-echo sequence with a slowly varying RF flip-angle train and Cartesian readouts, and the quantitative tissue maps are reconstructed by an iterative, model-based optimization algorithm. In this work, we design a three-dimensional (3D) MR-STAT framework based on previous 2D work, in order to achieve better image signal-to-noise ratio, higher though-plane resolution, and better tissue characterization. Specifically, we design a 7-min, high-resolution 3D MR-STAT sequence, and the corresponding two-step reconstruction algorithm for the large-scale dataset. To reduce the long acquisition time, Cartesian undersampling strategies such as SENSE are adopted in our transient-state quantitative framework. To reduce the computational burden, a data-splitting scheme is designed for decoupling the 3D reconstruction problem into independent 2D reconstructions. The proposed 3D framework is validated by numerical simulations, phantom experiments, and in vivo experiments. High-quality knee quantitative maps with 0.8 × 0.8 × 1.5 mm3 resolution and bilateral lower leg maps with 1.6 mm isotropic resolution can be acquired using the proposed 7-min acquisition sequence and the 3-min-per-slice decoupled reconstruction algorithm. The proposed 3D MR-STAT framework could have wide clinical applications in the future.

3D magnetic resonance fingerprinting on a low-field 50 mT point-of-care system prototype: evaluation of muscle and lipid relaxation time mapping and comparison with standard techniques

OBJECTIVE

To implement magnetic resonance fingerprinting (MRF) on a permanent magnet 50 mT low-field system deployable as a future point-of-care (POC) unit and explore the quality of the parameter maps.

MATERIALS AND METHODS

3D MRF was implemented on a custom-built Halbach array using a slab-selective spoiled steady-state free precession sequence with 3D Cartesian readout. Undersampled scans were acquired with different MRF flip angle patterns and reconstructed using matrix completion and matched to the simulated dictionary, taking excitation profile and coil ringing into account. MRF relaxation times were compared to that of inversion recovery (IR) and multi-echo spin echo (MESE) experiments in phantom and in vivo. Furthermore, B₀ inhomogeneities were encoded in the MRF sequence using an alternating TE pattern, and the estimated map was used to correct for image distortions in the MRF images using a model-based reconstruction.

RESULTS

Phantom relaxation times measured with an optimized MRF sequence for low field were in better agreement with reference techniques than for a standard MRF sequence. In vivo muscle relaxation times measured with MRF were longer than those obtained with an IR sequence (T₁: 182 ± 21.5 vs 168 ± 9.89 ms) and with an MESE sequence (T₂: 69.8 ± 19.7 vs 46.1 ± 9.65 ms). In vivo lipid MRF relaxation times were also longer compared with IR (T₁: 165 ± 15.1 ms vs 127 ± 8.28 ms) and with MESE (T₂: 160 ± 15.0 ms vs 124 ± 4.27 ms). Integrated ΔB₀ estimation and correction resulted in parameter maps with reduced distortions.

DISCUSSION

It is possible to measure volumetric relaxation times with MRF at 2.5 × 2.5 × 3.0 mm³ resolution in a 13 min scan time on a 50 mT permanent magnet system. The measured MRF relaxation times are longer compared to those measured with reference techniques, especially for T₂. This discrepancy can potentially be addressed by hardware, reconstruction and sequence design, but long-term reproducibility needs to be further improved.

Fast Quantitative Low-Field Magnetic Resonance Imaging With OPTIMUM-Optimized Magnetic Resonance Fingerprinting Using a Stationary Steady-State Cartesian Approach and Accelerated Acquisition Schedules

OBJECTIVE

The aim of the proposed work is to develop model-based, fast multiparametric magnetic resonance imaging (MRI) in field regimes where signal-to-noise ratio is poor, such as encountered at low-field and in low γ nuclei.

MATERIALS AND METHODS

A custom, optimized MRI pipeline was developed at low field (0.1 T) that relies on the magnetic resonance fingerprinting framework, called OPTIMUM. An optimization algorithm was used to select a short acquisition schedule (n = 18 images) that favors maximal discrimination across varying magnetic properties (T1, T2) and off-resonance effects while maintaining high transverse magnetization at the steady state. In the presented study, a stationary balanced steady-state approach was investigated that allows for Cartesian (used here) and non-Cartesian acquisition schemes. Images were collected in calibrated samples containing different concentrations of manganese(II) chloride (MnCl2) in deionized water and compared with gold standard techniques (ie, inversion recovery for T1, Carr-Purcell-Meiboom-Gill for T2). Images were then collected in vivo in the human hand and wrist.

RESULTS

OPTIMUM successfully provided sets of quantified maps (T1, T2, T2*, M0, ΔB0, B1+) in calibrated samples and in vivo in the human hand and wrist in 3 dimensions, in ~8.5 minutes, with a voxel resolution of [1.5 ×1.5 × 6.5] mm3. Relaxation parameters (T1, T2) scale linearly with [MnCl2] and are in good agreement with the calibrations performed for T1, with a consistent trend to underestimate T2.

CONCLUSION

We show that low-field MRI can benefit from innovative multiparametric approaches to gain speed and become realistic in clinical environments. For the first time, we report simultaneous, multiparametric imaging (6 quantitative maps) in 3 dimensions, in vivo in the human hand and wrist, obtained in just 8.5 minutes. It is sometimes overlooked that low magnetic fields provide higher dispersion of nuclear spin relaxation rates. Rapid quantification such as offered by OPTIMUM could be an enabling technology to explore new metrics and contrasts in point-of-care MRI diagnosis, making it an important step toward broad democratization.

Free-breathing abdominal T1 mapping using an optimized MR fingerprinting sequence

In this work, we propose a free-breathing magnetic resonance fingerprinting (MRF) method that can be used to obtain B₁⁺ -robust quantitative T₁ maps of the abdomen in a clinically acceptable time. A three-dimensional MRF sequence with a radial stack-of-stars trajectory was implemented, and its k-space acquisition ordering was adjusted to improve motion-robustness in the context of MRF. The flip angle pattern was optimized using the Cramér-Rao Lower Bound, and the encoding efficiency of sequences with 300, 600, 900 and 1800 flip angles was evaluated. To validate the sequence, a movable multicompartment phantom was developed. Reference multiparametric maps were acquired under stationary conditions using a previously validated MRF method. Periodic motion of the phantom was used to investigate the motion-robustness of the proposed sequence. The best performing sequence length (600 flip angles) was used to image the abdomen during a free-breathing volunteer scan. When using a series of 600 or more flip angles, the estimated T₁ values in the stationary phantom showed good agreement with the reference scan. Phantom experiments revealed that motion-related artifacts can appear in the quantitative maps and confirmed that a motion-robust k-space ordering is essential. The in vivo scan demonstrated that the proposed sequence can produce clean parameter maps while the subject breathes freely. Using this sequence, it is possible to generate B₁⁺ -robust quantitative maps of T₁ and B₁⁺ next to M₀ -weighted images under free-breathing conditions at a clinically usable resolution within 5 min.

Whole-brain 3D MR fingerprinting brain imaging: clinical validation and feasibility to patients with meningioma

Purpose

MR fingerprinting (MRF) is a MR technique that allows assessment of tissue relaxation times. The purpose of this study is to evaluate the clinical application of this technique in patients with meningioma.

Materials and methods

A whole-brain 3D isotropic 1mm³ acquisition under a 3.0T field strength was used to obtain MRF T₁ and T₂-based relaxometry values in 4:38 s. The accuracy of values was quantified by scanning a quantitative MR relaxometry phantom. In vivo evaluation was performed by applying the sequence to 20 subjects with 25 meningiomas. Regions of interest included the meningioma, caudate head, centrum semiovale, contralateral white matter and thalamus. For both phantom and subjects, mean values of both T₁ and T₂ estimates were obtained. Statistical significance of differences in mean values between the meningioma and other brain structures was tested using a Friedman’s ANOVA test.

Results

MR fingerprinting phantom data demonstrated a linear relationship between measured and reference relaxometry estimates for both T₁ (r² = 0.99) and T₂ (r² = 0.97). MRF T₁ relaxation times were longer in meningioma (mean ± SD 1429 ± 202 ms) compared to thalamus (mean ± SD 1054 ± 58 ms; p = 0.004), centrum semiovale (mean ± SD 825 ± 42 ms; p < 0.001) and contralateral white matter (mean ± SD 799 ± 40 ms; p < 0.001). MRF T₂ relaxation times were longer for meningioma (mean ± SD 69 ± 27 ms) as compared to thalamus (mean ± SD 27 ± 3 ms; p < 0.001), caudate head (mean ± SD 39 ± 5 ms; p < 0.001) and contralateral white matter (mean ± SD 35 ± 4 ms; p < 0.001)

Conclusions

Phantom measurements indicate that the proposed 3D-MRF sequence relaxometry estimations are valid and reproducible. For in vivo, entire brain coverage was obtained in clinically feasible time and allows quantitative assessment of meningioma in clinical practice.

Investigation of the effect of acquisition schemes on time-resolved magnetic resonance fingerprinting

Purpose.This study aims to investigate the feasibility of different acquisition methods for time-resolved magnetic resonance fingerprinting (TR-MRF) in computer simulation.Methods.An extended cardiac-torso (XCAT) phantom is used to generate abdominal T1, T2, and proton density maps for MRF simulation. The simulated MRF technique consists of an IR-FISP MRF sequence with spiral trajectory acquisition. MRF maps were simulated with different numbers of repetitions from 1 to 15. Three different methods were used to generate TR-MRF maps: (1) continuous acquisition without delay between MRF repetitions; (2) continuous acquisition with 5 s delay between MRF repetitions; (3) triggered acquisition with variable delay between MRF repetitions to allow the next acquisition to start at different respiration phase. After the generation of TR-MRF maps, the image quality indexes including the absolute T1 and T2 values, signal-to-noise-ratio (SNR), tumor-to-liver contrast-to-noise ratio, error in the amplitude of diaphragm motion and tumor volume error were used to evaluate the reconstructed parameter maps. Three volunteers were recruited to test the feasibility of the selected acquisition method.Results.Dynamic MR parametric maps using three different acquisition methods were estimated. The overall and liver T1 value error, liver SNR in T1 and T2 maps, and tumor SNR from T1 maps from triggered method is statistically significantly better than the other two methods (p-value < 0.05). The other image quality indexes have no significant difference between the triggered method and the other two continuous acquisition methods. All image quality indexes exhibit no significant difference between the acquisition methods with 0 s and 5 s delay. The triggered method was successfully performed in three healthy volunteers.Conclusion.TR-MRF technique was investigated using three different acquisition methods in computer simulation where the triggered method showed better performance than the other two methods. The triggered method has been tested successfully in healthy volunteers.

Simultaneous quantification of SPIO and gadolinium contrast agents using MR fingerprinting

PURPOSE

Develop a magnetic resonance fingerprinting (MRF) methodology with R₂^∗ quantification, intended for use with simultaneous contrast agent concentration mapping, particularly gadolinium (Gd) and iron labelled CD8+ T cells.

METHODS

Variable-density spiral SSFP MRF was used, modified to allow variable TE, and with an exp.(-TE·R₂^∗) dictionary modulation. In vitro phantoms containing SPIO labelled cells and/or gadolinium were used to validate parameter maps, probe undersampling capacity, and verify dual quantification capabilities. A C57BL/6 mouse was imaged using MRF to demonstrate acceptable in vivo resolution and signal at 8× undersampling necessary for a 25-min scan.

RESULTS

Strong agreement was found between conventional and MRF-derived values for R₁, R₂, and R₂^∗. Expanded MRF allowed quantification of iron-loaded CD8+ T cells. Results were robust to 8× undersampling and enabled recreation of relaxation profiles for both a Gd agent and iron labelled cells simultaneously. In vivo data demonstrated sufficient SNR in undersampled data for parameter mapping to visualise key features.

CONCLUSION

MRF can be expanded to include R₁, R₂, and R₂^∗ mapping required for simultaneous quantification of gadolinium and SPIO in vitro, allowing for potential implementation of a variety of future in vivo studies using dual MR contrast agents, including molecular imaging of labelled cells.

Water-fat separation in spiral magnetic resonance fingerprinting for high temporal resolution tissue relaxation time quantification in muscle

Purpose

To minimize the known biases introduced by fat in rapid T₁ and T₂ quantification in muscle using a single‐run magnetic resonance fingerprinting (MRF) water–fat separation sequence.

Methods

The single‐run MRF acquisition uses an alternating in‐phase/out‐of‐phase TE pattern to achieve water–fat separation based on a 2‐point DIXON method. Conjugate phase reconstruction and fat deblurring were applied to correct for B₀ inhomogeneities and chemical shift blurring. Water and fat signals were matched to the on‐resonance MRF dictionary. The method was first tested in a multicompartment phantom. To test whether the approach is capable of measuring small in vivo dynamic changes in relaxation times, experiments were run in 9 healthy volunteers; parameter values were compared with and without water–fat separation during muscle recovery after plantar flexion exercise.

Results

Phantom results show the robustness of the water–fat resolving MRF approach to undersampling. Parameter maps in volunteers show a significant (P < .01) increase in T₁ (105 ± 94 ms) and decrease in T₂ (14 ± 6 ms) when using water–fat‐separated MRF, suggesting improved parameter quantification by reducing the well‐known biases introduced by fat. Exercise results showed smooth T₁ and T₂ recovery curves.

Conclusion

Water–fat separation using conjugate phase reconstruction is possible within a single‐run MRF scan. This technique can be used to rapidly map relaxation times in studies requiring dynamic scanning, in which the presence of fat is problematic.

Unipolar MR elastography: Theory, numerical analysis and implementation

In MR elastography (MRE), zeroth moment balanced motion-encoding gradients (MEGs) are incorporated into MRI sequences to induce a phase shift proportional to the local displacement caused by external actuation. To maximize the signal-to-noise ratio (SNR), fractional encoding is employed, i.e., the MEG duration is reduced below the wave period. Here, gradients encode primarily the velocity of the motion-reducing encoding efficiency. Thus, in GRE-MRE, T₂ * decay and motion sensitivity have to be balanced, imposing a lower limit on repetition times (TRs). We propose to use a single trapezoidal gradient, a "unipolar gradient", to directly encode spin displacement. Such gradients cannot be used in conventional sequences as they exhibit a large zeroth moment and dephase magnetization. By time-reversing a spoiled SSFP sequence, the spoiling gradient becomes an efficient unipolar MEG. The proposed "unipolar MRE" technique benefits from this approach in three ways: first, displacement encoding is split over multiple TRs increasing motion sensitivity; second, spoiler and MEG coincide, allowing a reduction in TR; third, motion sensitivity of a typical unipolar lobe is of an order of magnitude higher than a bipolar MEG of equal duration. In this work, motion encoding using unipolar MRE is analyzed using the extended phase graph (EPG) formalism with a periodic motion propagator. As an approximation, the two-transverse TR approximation for diffusion-weighted SSFP is extended to incorporate cyclic motion. A complex encoding efficiency metric is introduced to compare the displacement fields of unipolar and conventional GRE-MRE sequences in both magnitude and phase. The derived theoretical encoding equations are used to characterize the proposed sequence using an extensive parameter study. Unipolar MRE is validated against conventional GRE-MRE in a phantom study showing excellent agreement between measured displacement fields. In addition, unipolar MRE yields significantly increased octahedral shear strain-SNR relative to conventional GRE-MRE and allows for the recovery of high stiffness inclusions, where conventional GRE-MRE fails.

Magnetic Resonance Fingerprinting: Implications and Opportunities for PET/MR

Magnetic Resonance Imaging (MRI) can be used to assess anatomical structure, and its sensitivity to a variety of tissue properties enables superb contrast between tissues as well as the ability to characterize these tissues. However, despite vast potential for quantitative and functional evaluation, MRI is typically used qualitatively, in which the underlying tissue properties are not measured, and thus the brightness of each pixel is not quantitatively meaningful. Positron Emission Tomography (PET) is an inherently quantitative imaging modality that interrogates functional activity within a tissue, probed by a molecule of interest coupled with an appropriate tracer. These modalities can complement one another to provide clinical information regarding both structure and function, but there are still technical and practical hurdles in the way of the integrated use of both modalities. Recent advances in MRI have moved the field in an increasingly quantitative direction, which is complementary to PET, and could also potentially help solve some of the challenges in PET/MR. Magnetic Resonance Fingerprinting (MRF) is a recently described MRI-based technique which can efficiently and simultaneously quantitatively map several tissue properties in a single exam. Here, the basic principles behind the quantitative approach of MRF are laid out, and the potential implications for combined PET/MR are discussed.

Cartesian MR fingerprinting in the eye at 7T using compressed sensing and matrix completion-based reconstructions

PURPOSE

To explore the feasibility of MR Fingerprinting (MRF) to rapidly quantify relaxation times in the human eye at 7T, and to provide a data acquisition and processing framework for future tissue characterization in eye tumor patients.

METHODS

In this single-element receive coil MRF approach with Cartesian sampling, undersampling is used to shorten scan time and, therefore, to reduce the degree of motion artifacts. For reconstruction, approaches based on compressed sensing (CS) and matrix completion (MC) were used, while their effects on the quality of the MRF parameter maps were studied in simulations and experiments. Average relaxation times in the eye were measured in 6 healthy volunteers. One uveal melanoma patient was included to show the feasibility of MRF in a clinical context.

RESULTS

Simulation results showed that an MC-based reconstruction enables large undersampling factors and also results in more accurate parameter maps compared with using CS. Experiments in 6 healthy volunteers used a reduction in scan time from 7:02 to 1:16 min, producing images without visible loss of detail in the parameter maps when using the MC-based reconstruction. Relaxation times from 6 healthy volunteers are in agreement with values obtained from fully sampled scans and values in literature, and parameter maps in a uveal melanoma patient show clear difference in relaxation times between tumor and healthy tissue.

CONCLUSION

Cartesian-based MRF is feasible in the eye at 7T. High undersampling factors can be achieved by means of MC, significantly shortening scan time and increasing patient comfort, while also mitigating the risk of motion artifacts.

Time efficient whole-brain coverage with MR Fingerprinting using slice-interleaved echo-planar-imaging

Magnetic resonance fingerprinting (MRF) is a promising method for fast simultaneous quantification of multiple tissue parameters. The objective of this study is to improve the coverage of MRF based on echo-planar imaging (MRF-EPI) by using a slice-interleaved acquisition scheme. For this, the MRF-EPI is modified to acquire several slices in a randomized interleaved manner, increasing the effective repetition time of the spoiled gradient echo readout acquisition in each slice. Per-slice matching of the signal-trace to a precomputed dictionary allows the generation of T1 and T2* maps with integrated B1+ correction. Subsequent compensation for the coil sensitivity profile and normalization to the cerebrospinal fluid additionally allows for quantitative proton density (PD) mapping. Numerical simulations are performed to optimize the number of interleaved slices. Quantification accuracy is validated in phantom scans and feasibility is demonstrated in-vivo. Numerical simulations suggest the acquisition of four slices as a trade-off between quantification precision and scan-time. Phantom results indicate good agreement with reference measurements (Difference T1: -2.4 ± 1.1%, T2*: -0.5 ± 2.5%, PD: -0.5 ± 7.2%). In-vivo whole-brain coverage of T1, T2* and PD with 32 slices was acquired within 3:36 minutes, resulting in parameter maps of high visual quality and comparable performance with single-slice MRF-EPI at 4-fold scan-time reduction.