Magnetic resonance fingerprinting using echo-planar imaging: Joint quantification of T1 and T2∗ relaxation times

Rapid simultaneous estimation of relaxation rates using multi-echo, multi-contrast MRI

PURPOSE

Multi-echo, multi-contrast methods are increasingly used in dynamic imaging studies to simultaneously quantify R2∗ and R2. To overcome the computational challenges associated with nonlinear least squares (NLSQ) fitting, we propose a generalized linear least squares (LLSQ) solution to rapidly fit R2∗ and R2.

METHODS

Spin- and gradient-echo (SAGE) data were simulated across T2∗ and T2 values at high (200) and low (20) SNR. Full (four-parameter) and reduced (three-parameter) parameter fits were implemented and compared with both LLSQ and NLSQ fitting. Fit data were compared to ground truth using concordance correlation coefficient (CCC) and coefficient of variation (CV). In vivo SAGE perfusion data were acquired in 20 subjects with relapsing-remitting multiple sclerosis. LLSQ R2∗ and R2, as well as cerebral blood volume (CBV), were compared with the standard NLSQ approach.

RESULTS

Across all fitting methods, T2∗ was well-fit at high (CCC = 1, CV = 0) and low (CCC ≥ 0.87, CV ≤ 0.08) SNR. Except for short T2∗ values (5-15 ms), T2 was well-fit at high (CCC = 1, CV = 0) and low (CCC ≥ 0.99, CV ≤ 0.03) SNR. In vivo, LLSQ R2∗ and R2 estimates were similar to NLSQ, and there were no differences in R2∗ across fitting methods at high SNR. However, there were some differences at low SNR and for R2 at high and low SNR. In vivo NLSQ and LLSQ three parameter fits performed similarly, as did NLSQ and LLSQ four-parameter fits. LLSQ CBV nearly matched the standard NLSQ method for R2∗- (0.97 ratio) and R2-CBV (0.98 ratio). Voxel-wise whole-brain fitting was faster for LLSQ (3-4 min) than NLSQ (16-18 h).

CONCLUSIONS

LLSQ reliably fit for R2∗ and R2 in simulated and in vivo data. Use of LLSQ methods reduced the computational demand, enabling rapid estimation of R2∗ and R2.

Phase-sensitive deep reconstruction method for rapid multiparametric MR fingerprinting and quantitative susceptibility mapping in the brain

INTRODUCTION

This study explores the potential of Magnetic Resonance Fingerprinting (MRF) with a novel Phase-Sensitivity Deep Reconstruction Network (PS-DRONE) for simultaneous quantification of T1, T2, Proton Density, B1+, phase and quantitative susceptibility mapping (QSM).

METHODS

Data were acquired at 3 T in vitro and in vivo using an optimized EPI-based MRF sequence. Phantom experiments were conducted using a standardized phantom for T1 and T2 maps and a custom-made agar-based gadolinium phantom for B1 and QSM maps. In vivo experiments included five healthy volunteers and one patient diagnosed with brain metastasis. PSDRONE maps were compared to reference maps obtained through standard imaging sequences.

RESULTS

Total scan time was 2 min for 32 slices and a resolution of [1 mm, 1 mm, 4.5 mm]. The reconstruction of T1, T2, Proton Density, B1+ and phase maps were reconstructed within 1 s. In the phantoms, PS-DRONE analysis presented accurate and strongly correlated T1 and T2 maps (r = 0.99) compared to the reference maps. B1 maps from PS-DRONE showed slightly higher values, though still correlated (r = 0.6) with the reference. QSM values showed a small bias but were strongly correlated (r = 0.99) with reference data. In the in vivo analysis, PS-DRONE-derived T1 and T2 values for gray and white matter matched reference values in healthy volunteers. PS-DRONE B1 and QSM maps showed strong correlations with reference values.

CONCLUSION

The PS-DRONE network enables concurrent acquisition of T1, T2, PD, B1+, phase and QSM maps, within 2 min of acquisition time and 1 s of reconstruction time.

Fast deep learning reconstruction techniques for preclinical magnetic resonance fingerprinting

We propose a deep learning (DL) model and a hyperparameter optimization strategy to reconstruct T1 and T2 maps acquired with the magnetic resonance fingerprinting (MRF) methodology. We applied two different MRF sequence routines to acquire images of ex vivo rat brain phantoms using a 7-T preclinical scanner. Subsequently, the DL model was trained using experimental data, completely excluding the use of any theoretical MRI signal simulator. The best combination of the DL parameters was implemented by an automatic hyperparameter optimization strategy, whose key aspect is to include all the parameters to the fit, allowing the simultaneous optimization of the neural network architecture, the structure of the DL model, and the supervised learning algorithm. By comparing the reconstruction performances of the DL technique with those achieved from the traditional dictionary-based method on an independent dataset, the DL approach was shown to reduce the mean percentage relative error by a factor of 3 for T1 and by a factor of 2 for T2 , and to improve the computational time by at least a factor of 37. Furthermore, the proposed DL method enables maintaining comparable reconstruction performance, even with a lower number of MRF images and a reduced k-space sampling percentage, with respect to the dictionary-based method. Our results suggest that the proposed DL methodology may offer an improvement in reconstruction accuracy, as well as speeding up MRF for preclinical, and in prospective clinical, investigations.

Rapid Whole-Knee Quantification of Cartilage Using T1, T2*, and TRAFF2 Mapping With Magnetic Resonance Fingerprinting

OBJECTIVE

Quantitative Magnetic Resonance Imaging (MRI) holds great promise for the early detection of cartilage deterioration. Here, a Magnetic Resonance Fingerprinting (MRF) framework is proposed for comprehensive and rapid quantification of T1, T2*, and TRAFF2 with whole-knee coverage.

METHODS

A MRF framework was developed to achieve quantification of Relaxation Along a Fictitious Field in the 2nd rotating frame of reference ( TRAFF2) along with T1 and T2*. The proposed sequence acquires 65 measurements of 25 high-resolution slices, interleaved with 7 inversion pulses and 40 RAFF2 trains, for whole-knee quantification in a total acquisition time of 3:25 min. Comparison with reference T1, T2*, and TRAFF2 methods was performed in phantom and in seven healthy subjects at 3 T. Repeatability (test-retest) with and without repositioning was also assessed.

RESULTS

Phantom measurements resulted in good agreement between MRF and the reference with mean biases of -54, 2, and 5 ms for T1, T2*, and TRAFF2, respectively. Complete characterization of the whole-knee cartilage was achieved for all subjects, and, for the femoral and tibial compartments, a good agreement between MRF and reference measurements was obtained. Across all subjects, the proposed MRF method yielded acceptable repeatability without repositioning ( R2 ≥ 0.94) and with repositioning ( R2 ≥ 0.57) for T1, T2*, and TRAFF2.

SIGNIFICANCE

The short scan time combined with the whole-knee coverage makes the proposed MRF framework a promising candidate for the early assessment of cartilage degeneration with quantitative MRI, but further research may be warranted to improve repeatability after repositioning and assess clinical value in patients.

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.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Quantification of the Intrinsic T1 and T2 of Heschl's Gyri with MR Fingerprinting

PURPOSE

The human primary auditory cortex is located in the Heschl's gyrus (HG). To assess the intrinsic MR property in the gray matter of the HG (GM-HG) with T1 and T2 values using a commercially available MR fingerprinting (MRF) technique.

METHODS

The subjects were 10 healthy volunteers (with 20 HGs; mean age, 31.5 years old; range, 25-53 years old). Coronal T1 and T2 maps were obtained with commercially available MRF using a 3-Tesla MR system. Two radiologists measured the T1 and T2 values of the GM-HG, the GM in the superior temporal gyrus (GM-STG), and the GM in the middle temporal gyrus (GM-MTG) by drawing a ROI on coronal maps.

RESULTS

For both radiologists, the mean T1 and T2 values of the GM-HG were significantly lower than those in the GM-STG or GM-MTG (P < 0.01). The interobserver reliability using the intraclass correlation coefficients (ICC) (2,1) showed strong agreement for the measurement of the T1 and T2 values (ICCs =⃥ 0.80 and 0.78 for T1 and T2 values, respectively).

CONCLUSION

The T1 and T2 values on MRF for the GM-HG were lower than those for the GM-STG and GM-MTG, likely reflecting a higher myelin content and iron deposition in the GM-HG. Quantitative measurements using the MRF can clarify cortical properties with high reliability, which may indicate that MRF mapping provides new insights into the structure of the human cortical GM.

Acceleration of Magnetic Resonance Fingerprinting Reconstruction Using Denoising and Self-Attention Pyramidal Convolutional Neural Network

Magnetic resonance fingerprinting (MRF) based on echo-planar imaging (EPI) enables whole-brain imaging to rapidly obtain T1 and T2* relaxation time maps. Reconstructing parametric maps from the MRF scanned baselines by the inner-product method is computationally expensive. We aimed to accelerate the reconstruction of parametric maps for MRF-EPI by using a deep learning model. The proposed approach uses a two-stage model that first eliminates noise and then regresses the parametric maps. Parametric maps obtained by dictionary matching were used as a reference and compared with the prediction results of the two-stage model. MRF-EPI scans were collected from 32 subjects. The signal-to-noise ratio increased significantly after the noise removal by the denoising model. For prediction with scans in the testing dataset, the mean absolute percentage errors between the standard and the final two-stage model were 3.1%, 3.2%, and 1.9% for T1, and 2.6%, 2.3%, and 2.8% for T2* in gray matter, white matter, and lesion locations, respectively. Our proposed two-stage deep learning model can effectively remove noise and accurately reconstruct MRF-EPI parametric maps, increasing the speed of reconstruction and reducing the storage space required by dictionaries.

Magnetic Resonance Imaging compatible Elastic Loading Mechanism (MELM): A minimal footprint device for MR imaging under load

Quantitative Magnetic Resonance Imaging (MRI) can enable early diagnosis of knee cartilage damage if imaging is performed during the application of load. Mechanical loading via ropes, pulleys and suspended weights can be obstructive and require adaptations to the patient table. In this paper, a new lightweight MRI-compatible elastic loading mechanism is introduced. The new device showed sufficient linearity (|α/β| = 0.42 ± 0.25), reproducibility (CoV = 5 ± 2%), and stability (CoV = 0.5 ± 0.1%). In vivo and ex vivo scans confirmed the ability of the device to exert sufficient force to study the knee cartilage under loading conditions, inducing up to a 29% decrease in $T_2^{\ast}$ of the central medial cartilage. With this device mechanical loading can become more accessible for researchers and clinicians, thus facilitating the translational use of MRI biomarkers for the detection of cartilage deterioration.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

3D MR fingerprinting (MRF) for simultaneous T1 and T2 quantification of the bone metastasis: Initial validation in prostate cancer patients

PURPOSE

To investigate the feasibility of using 3-dimensional MRF for bone marrow evaluation in the field of view of prostate MRI for T1 and T2 quantification of prostate cancer bone metastases, as well as comparing it to the ADC value.

METHODS

In this retrospective study, 30 prostate MRIs were included: 14 cases with prostate cancer bone metastasis and 16 cases without prostate cancer (control). MRF was obtained twice before (nonenhanced [NE] MRF) and after contrast injection (contrast-enhanced [CE] MRF), and T1 and T2 maps were generated from each MRF. Two radiologists independently drew regions of interest (ROIs) on the MRF maps and the ADC maps. Mann-Whitney U tests and the area under the receiver operating characteristic curve (AUROC) evaluated the two-reader means of T1, T2 and ADC values between bone metastasis and normal bone.

RESULTS

There were 83 ROIs, including 39 bone metastases and 44 normal bone. The two-reader average ADC, NE T2 and CE T2 values were significantly lower and NE T1 and CE T1 values were significantly higher in metastatic bone compared with normal bone (P < 0.001). The AUROC of the ADC was lowest (0.685), which was significantly lower than those of NE T1 (1.0, P = 0.001), NE T2 (0.932, P = 0.004), and CE T2 (0.876, P = 0.031).

CONCLUSION

MRF to assess the pelvic bone during a prostate gland evaluation provides a reliable parametric map for skeletal work-up. With higher diagnostic performance than the ADC value, NE MRF is a potential alternative for quantifying bone marrow metastases in prostate cancer patients.

A deep learning approach for synthetic MRI based on two routine sequences and training with synthetic data

BACKGROUND AND OBJECTIVE

Synthetic magnetic resonance imaging (MRI) is a low cost procedure that serves as a bridge between qualitative and quantitative MRI. However, the proposed methods require very specific sequences or private protocols which have scarcely found integration in clinical scanners. We propose a learning-based approach to compute T1, T2, and PD parametric maps from only a pair of T1- and T2-weighted images customarily acquired in the clinical routine.

METHODS

Our approach is based on a convolutional neural network (CNN) trained with synthetic data; specifically, a synthetic dataset with 120 volumes was constructed from the anatomical brain model of the BrainWeb tool and served as the training set. The CNN learns an end-to-end mapping function to transform the input T1- and T2-weighted images to their underlying T1, T2, and PD parametric maps. Then, conventional weighted images unseen by the network are analytically synthesized from the parametric maps. The network can be fine tuned with a small database of actual weighted images and maps for better performance.

RESULTS

This approach is able to accurately compute parametric maps from synthetic brain data achieving normalized squared error values predominantly below 1%. It also yields realistic parametric maps from actual MR brain acquisitions with T1, T2, and PD values in the range of the literature and with correlation values above 0.95 compared to the T1 and T2 maps obtained from relaxometry sequences. Further, the synthesized weighted images are visually realistic; the mean square error values are always below 9% and the structural similarity index is usually above 0.90. Network fine tuning with actual maps improves performance, while training exclusively with a small database of actual maps shows a performance degradation.

CONCLUSIONS

These results show that our approach is able to provide realistic parametric maps and weighted images out of a CNN that (a) is trained with a synthetic dataset and (b) needs only two inputs, which are in turn obtained from a common full-brain acquisition that takes less than 8 min of scan time. Although a fine tuning with actual maps improves performance, synthetic data is crucial to reach acceptable performance levels. Hence, we show the utility of our approach for both quantitative MRI in clinical viable times and for the synthesis of additional weighted images to those actually acquired.

Efficiency analysis for quantitative MRI of T1 and T2 relaxometry methods

This study presents a comparison of quantitative MRI methods based on an efficiency metric that quantifies their intrinsic ability to extract information about tissue parameters. Under a regime of unbiased parameter estimates, an intrinsic efficiency metricηwas derived for fully-sampled experiments which can be used to both optimize and compare sequences. Here we optimize and compare several steady-state and transient gradient-echo based qMRI methods, such as magnetic resonance fingerprinting (MRF), for jointT1andT2mapping. The impact of undersampling was also evaluated, assuming incoherent aliasing that is treated as noise by parameter estimation.In vivovalidation of the efficiency metric was also performed. Transient methods such as MRF can be up to 3.5 times more efficient than steady-state methods, when spatial undersampling is ignored. If incoherent aliasing is treated as noise during least-squares parameter estimation, the efficiency is reduced in proportion to the SNR of the data, with reduction factors of 5 often seen for practical SNR levels.In vivovalidation showed a very good agreement between the theoretical and experimentally predicted efficiency. This work presents and validates an efficiency metric to optimize and compare the performance of qMRI methods. Transient methods were found to be intrinsically more efficient than steady-state methods, however the effect of spatial undersampling can significantly erode this advantage.

Lesion probability mapping in MS patients using a regression network on MR fingerprinting

BACKGROUND

To develop a regression neural network for the reconstruction of lesion probability maps on Magnetic Resonance Fingerprinting using echo-planar imaging (MRF-EPI) in addition to [Formula: see text], [Formula: see text], NAWM, and GM- probability maps.

METHODS

We performed MRF-EPI measurements in 42 patients with multiple sclerosis and 6 healthy volunteers along two sites. A U-net was trained to reconstruct the denoised and distortion corrected [Formula: see text] and [Formula: see text] maps, and to additionally generate NAWM-, GM-, and WM lesion probability maps.

RESULTS

WM lesions were predicted with a dice coefficient of [Formula: see text] and a lesion detection rate of [Formula: see text] for a threshold of 33%. The network jointly enabled accurate [Formula: see text] and [Formula: see text] times with relative deviations of 5.2% and 5.1% and average dice coefficients of [Formula: see text] and [Formula: see text] for NAWM and GM after binarizing with a threshold of 80%.

CONCLUSION

DL is a promising tool for the prediction of lesion probability maps in a fraction of time. These might be of clinical interest for the WM lesion analysis in MS patients.

Streamlined magnetic resonance fingerprinting: Fast whole-brain coverage with deep-learning based parameter estimation

Magnetic resonance fingerprinting (MRF) is a quantitative MRI (qMRI) framework that provides simultaneous estimates of multiple relaxation parameters as well as metrics of field inhomogeneity in a single acquisition. However, current challenges exist in the forms of (1) scan time; (2) need for custom image reconstruction; (3) large dictionary sizes; (4) long dictionary-matching time. This study aims to introduce a novel streamlined magnetic-resonance fingerprinting (sMRF) framework based on a single-shot echo-planar imaging (EPI) sequence to simultaneously estimate tissue T1, T2, and T2* with integrated B1⁺ correction. Encouraged by recent work on EPI-based MRF, we developed a method that combines spin-echo EPI with gradient-echo EPI to achieve T2 in addition to T1 and T2* quantification. To this design, we add simultaneous multi-slice (SMS) acceleration to enable full-brain coverage in a few minutes. Moreover, in the parameter-estimation step, we use deep learning to train a deep neural network (DNN) to accelerate the estimation process by orders of magnitude. Notably, due to the high image quality of the EPI scans, the training process can rely simply on Bloch-simulated data. The DNN also removes the need for storing large dictionaries. Phantom scans along with in-vivo multi-slice scans from seven healthy volunteers were acquired with resolutions of 1.1×1.1×3 mm³ and 1.7×1.7×3 mm³, and the results were validated against ground truth measurements. Excellent correspondence was found between our T1, T2, and T2* estimates and results obtained from standard approaches. In the phantom scan, a strong linear relationship (R = 1-1.04, R²>0.96) was found for all parameter estimates, with a particularly high agreement for T2 estimation (R²>0.99). Similar findings are reported for the in-vivo human data for all of our parameter estimates. Incorporation of DNN results in a reduction of parameter estimation time on the order of 1000 x and a reduction in storage requirements on the order of 2500 x while achieving highly similar results as conventional dictionary matching (%differences of 7.4 ± 0.4%, 3.6 ± 0.3% and 6.0 ± 0.4% error in T1, T2, and T2* estimation). Thus, sMRF has the potential to be the method of choice for future MRF studies by providing ease of implementation, fast whole-brain coverage, and ultra-fast T1/T2/T2* estimation.

Magnetic resonance fingerprinting: an overview

Magnetic resonance fingerprinting (MRF) is an evolving quantitative MRI framework consisting of unique data acquisition, processing, visualization, and interpretation steps. MRF is capable of simultaneously producing multiple high-resolution property maps including T1, T2, M0, ADC, and T2* measurements. While a relatively new technology, MRF has undergone rapid development for a variety of clinical applications from brain tumor characterization and epilepsy imaging to characterization of prostate cancer, cardiac imaging, among others. This paper will provide a brief overview of current state of MRF technology including highlights of technical and clinical advances. We will conclude with a brief discussion of the challenges that need to be overcome to establish MRF as a quantitative imaging biomarker.

Accelerated white matter lesion analysis based on simultaneous T 1 and T 2 ∗ quantification using magnetic resonance fingerprinting and deep learning

PURPOSE

To develop an accelerated postprocessing pipeline for reproducible and efficient assessment of white matter lesions using quantitative magnetic resonance fingerprinting (MRF) and deep learning.

METHODS

MRF using echo-planar imaging (EPI) scans with varying repetition and echo times were acquired for whole brain quantification of T 1 and T 2 ∗ in 50 subjects with multiple sclerosis (MS) and 10 healthy volunteers along 2 centers. MRF T 1 and T 2 ∗ parametric maps were distortion corrected and denoised. A CNN was trained to reconstruct the T 1 and T 2 ∗ parametric maps, and the WM and GM probability maps.

RESULTS

Deep learning-based postprocessing reduced reconstruction and image processing times from hours to a few seconds while maintaining high accuracy, reliability, and precision. Mean absolute error performed the best for T 1 (deviations 5.6%) and the logarithmic hyperbolic cosinus loss the best for T 2 ∗ (deviations 6.0%).

CONCLUSIONS

MRF is a fast and robust tool for quantitative T 1 and T 2 ∗ mapping. Its long reconstruction and several postprocessing steps can be facilitated and accelerated using deep learning.

Compressive MRI quantification using convex spatiotemporal priors and deep encoder-decoder networks

We propose a dictionary-matching-free pipeline for multi-parametric quantitative MRI image computing. Our approach has two stages based on compressed sensing reconstruction and deep learned quantitative inference. The reconstruction phase is convex and incorporates efficient spatiotemporal regularisations within an accelerated iterative shrinkage algorithm. This minimises the under-sampling (aliasing) artefacts from aggressively short scan times. The learned quantitative inference phase is purely trained on physical simulations (Bloch equations) that are flexible for producing rich training samples. We propose a deep and compact encoder-decoder network with residual blocks in order to embed Bloch manifold projections through multi-scale piecewise affine approximations, and to replace the non-scalable dictionary-matching baseline. Tested on a number of datasets we demonstrate effectiveness of the proposed scheme for recovering accurate and consistent quantitative information from novel and aggressively subsampled 2D/3D quantitative MRI acquisition protocols.

Synthetic MRI: Technologies and Applications in Neuroradiology

Synthetic MRI is a technique that synthesizes contrast-weighted images from multicontrast MRI data. There have been advances in synthetic MRI since the technique was introduced. Although a number of synthetic MRI methods have been developed for quantifying one or more relaxometric parameters and for generating multiple contrast-weighted images, this review focuses on several methods that quantify all three relaxometric parameters (T₁ , T₂ , and proton density) and produce multiple contrast-weighted images. Acquisition, quantification, and image synthesis techniques are discussed for each method. We discuss the image quality and diagnostic accuracy of synthetic MRI methods and their clinical applications in neuroradiology. Based on this analysis, we highlight areas that need to be addressed for synthetic MRI to be widely implemented in the clinic. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

3D magnetic resonance fingerprinting with quadratic RF phase

PURPOSE

To implement 3D magnetic resonance fingerprinting (MRF) with quadratic RF phase (qRF-MRF) for simultaneous quantification of T₁ , T₂ , ΔB₀ , and T 2 ∗ .

METHODS

3D MRF data with effective undersampling factor of 3 in the slice direction were acquired with quadratic RF phase patterns for T₁ , T₂ , and T 2 ∗ sensitivity. Quadratic RF phase encodes the off-resonance by modulating the on-resonance frequency linearly in time. Transition to 3D brings practical limitations for reconstruction and dictionary matching because of increased data and dictionary sizes. Randomized singular value decomposition (rSVD)-based compression in time and reduction in dictionary size with a quadratic interpolation method are combined to be able to process prohibitively large data sets in feasible reconstruction and matching times.

RESULTS

Accuracy of 3D qRF-MRF maps in various resolutions and orientations are compared to 3D fast imaging with steady-state precession (FISP) for T₁ and T₂ contrast and to 2D qRF-MRF for T 2 ∗ contrast and ΔB₀ . The precision of 3D qRF-MRF was 1.5-2 times higher than routine clinical scans. 3D qRF-MRF ΔB₀ maps were further processed to highlight the susceptibility contrast.

CONCLUSION

Natively co-registered 3D whole brain T₁ , T₂ , T 2 ∗ , ΔB₀ , and QSM maps can be acquired in as short as 5 min with 3D qRF-MRF.

Magnetic resonance fingerprinting for simultaneous renal T1 and T 2 * mapping in a single breath-hold

PURPOSE

To evaluate the use of magnetic resonance fingerprinting (MRF) for simultaneous quantification of T 1 and T 2 ∗ in a single breath-hold in the kidneys.

METHODS

The proposed kidney MRF sequence was based on MRF echo-planar imaging. Thirty-five measurements per slice and overall 4 slices were measured in 15.4 seconds. Group matching was performed for in-line quantification of T 1 and T 2 ∗ . Images were acquired in a phantom and 8 healthy volunteers in coronal orientation. To evaluate our approach, region of interests were drawn in the kidneys to calculate mean values and standard deviations of the T 1 and T 2 ∗ times. Precision was calculated across multiple repeated MRF scans. Gaussian filtering is applied on baseline images to improve SNR and match stability.

RESULTS

T 1 and T 2 ∗ times acquired with MRF in the phantom showed good agreement with reference measurements and conventional mapping methods with deviations of less than 5% for T 1 and less than 10% for T 2 ∗ . Baseline images in vivo were free of artifacts and relaxation times yielded good agreement with conventional methods and literature (deviation T 1 : 7 ± 4 % , T 2 ∗ : 6 ± 3 % ).

CONCLUSIONS

In this feasibility study, the proposed renal MRF sequence resulted in accurate T 1 and T 2 ∗ quantification in a single breath-hold.

Machine Learning for Rapid Magnetic Resonance Fingerprinting Tissue Property Quantification

Magnetic Resonance Fingerprinting (MRF) is an MRI-based method that can provide quantitative maps of multiple tissue properties simultaneously from a single rapid acquisition. Tissue property maps are generated by matching the complex signal evolutions collected at the scanner to a dictionary of signals derived using Bloch equation simulations. However, in some circumstances, the process of dictionary generation and signal matching can be time-consuming, reducing the utility of this technique. Recently, several groups have proposed using machine learning to accelerate the extraction of quantitative maps from MRF data. This article will provide an overview of current research that combines MRF and machine learning, as well as present original research demonstrating how machine learning can speed up dictionary generation for cardiac MRF.

Dynamic, Simultaneous Concentration Mapping of Multiple MRI Contrast Agents with Dual Contrast - Magnetic Resonance Fingerprinting

Synchronous assessment of multiple MRI contrast agents in a single scanning session would provide a new "multi-color" imaging capability similar to fluorescence imaging but with high spatiotemporal resolution and unlimited imaging depth. This multi-agent MRI technology would enable a whole new class of basic science and clinical MRI experiments that simultaneously explore multiple physiologic/molecular events in vivo. Unfortunately, conventional MRI acquisition techniques are only capable of detecting and quantifying one paramagnetic MRI contrast agent at a time. Herein, the Dual Contrast - Magnetic Resonance Fingerprinting (DC-MRF) methodology was extended for in vivo application and evaluated by simultaneously and dynamically mapping the intra-tumoral concentration of two MRI contrast agents (Gd-BOPTA and Dy-DOTA-azide) in a mouse glioma model. Co-registered gadolinium and dysprosium concentration maps were generated with sub-millimeter spatial resolution and acquired dynamically with just over 2-minute temporal resolution. Mean tumor Gd and Dy concentration measurements from both single agent and dual agent DC-MRF studies demonstrated significant correlations with ex vivo mass spectrometry elemental analyses. This initial in vivo study demonstrates the potential for DC-MRF to provide a useful dual-agent MRI platform.

Flexible and efficient optimization of quantitative sequences using automatic differentiation of Bloch simulations

PURPOSE

To investigate a computationally efficient method for optimizing the Cramér-Rao Lower Bound (CRLB) of quantitative sequences without using approximations or an analytical expression of the signal.

METHODS

Automatic differentiation was applied to Bloch simulations and used to optimize several quantitative sequences without the need for approximations or an analytical expression. The results were validated with in vivo measurements and comparisons to prior art. Multi-echo spin echo and DESPO T 1 were used as benchmarks to verify the CRLB implementation. The CRLB of the Magnetic Resonance Fingerprinting (MRF) sequence, which has a complicated analytical formulation, was also optimized using automatic differentiation.

RESULTS

The sequence parameters obtained for multi-echo spin echo and DESPO T 1 matched results obtained using conventional methods. In vivo, MRF scans demonstrate that the CRLB optimization obtained with automatic differentiation can improve performance in presence of white noise. For MRF, the CRLB optimization converges in 1.1 CPU hours for N TR = 400 and has O ( N TR ) asymptotic runtime scaling for the calculation of the CRLB objective and gradient.

CONCLUSIONS

Automatic differentiation can be used to optimize the CRLB of quantitative sequences without using approximations or analytical expressions. For MRF, the runtime is computationally efficient and can be used to investigate confounding factors as well as MRF sequences with a greater number of repetitions.

Magnetic resonance fingerprinting Part 1: Potential uses, current challenges, and recommendations

Magnetic resonance fingerprinting (MRF) is a powerful quantitative MRI technique capable of acquiring multiple property maps simultaneously in a short timeframe. The MRF framework has been adapted to a wide variety of clinical applications, but faces challenges in technical development, and to date has only demonstrated repeatability and reproducibility in small studies. In this review, we discuss the current implementations of MRF and their use in a clinical setting. Based on this analysis, we highlight areas of need that must be addressed before MRF can be fully adopted into the clinic and make recommendations to the MRF community on standardization and validation strategies of MRF techniques. Level of Evidence: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:675-692.

Comprehensive Evaluation of B1+-corrected FISP-based Magnetic Resonance Fingerprinting: Accuracy, Repeatability and Reproducibility of T1 and T2 Relaxation Times for ISMRM/NIST System Phantom and Volunteers

Purpose:

This study aimed to evaluate comprehensively; accuracy, repeatability and reproducibility of T₁ and T₂ relaxation times measured by magnetic resonance fingerprinting using B1+-corrected fast imaging with steady-state precession (FISP–MRF).

Methods:

The International Society of Magnetic Resonance in Medicine/National Institute of Standards and Technology (ISMRM/NIST) phantom was scanned for 100 days, and six healthy volunteers for 5 days using a FISP–MRF prototype sequence. Accuracy was evaluated on the phantom by comparing relaxation times measured by FISP–MRF with the reference values provided by the phantom manufacturer. Daily repeatability was characterized as the coefficient of variation (CV) of the measurements over 100 days for the phantom and over 5 days for volunteers. In addition, the cross-scanner reproducibility was evaluated in volunteers.

Results:

In the phantom study, T₁ and T₂ values from FISP–MRF showed a strong linear correlation with the reference values of the phantom (R² = 0.9963 for T₁; R² = 0.9966 for T₂). CVs were <1.0% for T₁ values larger than 300 ms, and <3.0% for T₂ values across a wide range. In the volunteer study, CVs for both T₁ and T₂ values were <5.0%, except for one subject. In addition, all T₂ values estimated by FISP–MRF in vivo were lower than those measured with conventional mapping sequences reported in previous studies. The cross-scanner variation of T₁ and T₂ showed good agreement between two different scanners in the volunteers.

Conclusion:

B1+-corrected FISP-MRF showed an acceptable accuracy, repeatability and reproducibility in the phantom and volunteer studies.

Ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) for simultaneous quantification of long and ultrashort T2 tissues

PURPOSE

To demonstrate an ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) method that allows quantifying relaxation times for muscle and bone in the musculoskeletal system and generating bone enhanced images that mimic CT scans.

METHODS

A fast imaging steady-state free precession MRF sequence with half pulse excitation and half projection readout was designed to sample fast T₂ decay signals. Varying echo time (TE) of a sinusoidal pattern was applied to enhance sensitivity for tissues with short and ultrashort T₂ values. The performance of UTE-MRF was evaluated via simulations, phantom, and in vivo experiments.

RESULTS

A minimal TE of 0.05 ms was achieved. Simulations indicated the sinusoidal TE sampling increased T₂ quantification accuracy in the cortical bone and tendon but had little impact on long T₂ muscle quantifications. For the rubber phantom, the averaged relaxometries from UTE-MRF (T₁ = 162 ms and T₂ = 1.07 ms) compared well with the gold standard (T₁ = 190 ms and T 2 ∗ = 1.03 ms). For the long T₂ agarose phantom, the linear regression slope between UTE-MRF and gold standard was 1.07 (R² = 0.991) for T₁ and 1.04 (R² = 0.994) for T₂ . In vivo experiments showed the detection of the cortical bone (averaged T₂ = 1.0 ms) and Achilles tendon (averaged T₂ = 15 ms). Scalp structures from the bone enhanced image show high similarity with CT.

CONCLUSION

The UTE-MRF with sinusoidal TEs can simultaneously quantify T₁ , T₂ , proton density, and B₀ in long, short, even ultrashort T₂ musculoskeletal structures. Bone enhanced images can be achieved in the brain with UTE-MRF.

Optimal Experiment Design for Magnetic Resonance Fingerprinting: Cramér-Rao Bound Meets Spin Dynamics

Magnetic resonance (MR) fingerprinting is a new quantitative imaging paradigm, which simultaneously acquires multiple MR tissue parameter maps in a single experiment. In this paper, we present an estimation-theoretic framework to perform experiment design for MR fingerprinting. Specifically, we describe a discrete-time dynamic system to model spin dynamics, and derive an estimation-theoretic bound, i.e., the Cramér-Rao bound, to characterize the signal-to-noise ratio (SNR) efficiency of an MR fingerprinting experiment. We then formulate an optimal experiment design problem, which determines a sequence of acquisition parameters to encode MR tissue parameters with the maximal SNR efficiency, while respecting the physical constraints and other constraints from the image decoding/reconstruction process. We evaluate the performance of the proposed approach with numerical simulations, phantom experiments, and in vivo experiments. We demonstrate that the optimized experiments substantially reduce data acquisition time and/or improve parameter estimation. For example, the optimized experiments achieve about a factor of two improvement in the accuracy of T2 maps, while keeping similar or slightly better accuracy of T1 maps. Finally, as a remarkable observation, we find that the sequence of optimized acquisition parameters appears to be highly structured rather than randomly/pseudo-randomly varying as is prescribed in the conventional MR fingerprinting experiments.

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.

Magnetic resonance fingerprinting with dictionary-based fat and water separation (DBFW MRF): A multi-component approach

Purpose

To obtain a fast and robust fat‐water separation with simultaneous estimation of water T₁, fat T₁, and fat fraction maps.

Methods

We modified an MR fingerprinting (MRF) framework to use a single dictionary combination of a water and fat dictionary. A variable TE acquisition pattern with maximum TE = 4.8 ms was used to increase the fat–water separability. Radiofrequency (RF) spoiling was used to reduce the size of the dictionary by reducing T₂ sensitivity. The technique was compared both in vitro and in vivo to an MRF method that incorporated 3‐point Dixon (DIXON MRF), as well as Cartesian IDEAL with different acquisition parameters.

Results

The proposed dictionary‐based fat–water separation technique (DBFW MRF) successfully provided fat fraction, water, and fat T₁, B₀, and B₁₊ maps both in vitro and in vivo. The fat fraction and water T₁ values obtained with DBFW MRF show excellent agreement with DIXON MRF as well as with the reference values obtained using a Cartesian IDEAL with a long TR (concordance correlation coefficient: 0.97/0.99 for fat fraction–water T₁). Whereas fat fraction values with Cartesian IDEAL were degraded in the presence of T₁ saturation, MRF methods successfully estimated and accounted for T₁ in the fat fraction estimates.

Conclusion

The DBFW MRF technique can successfully provide T₁ and fat fraction quantification in under 20 s per slice, intrinsically correcting T₁ biases typical of fast Dixon techniques. These features could improve the diagnostic quality and use of images in presence of fat.

Magnetic resonance fingerprinting with quadratic RF phase for measurement of T2 * simultaneously with δf , T1 , and T2

PURPOSE

This study explores the possibility of using a gradient moment balanced sequence with a quadratically varied RF excitation phase in the magnetic resonance fingerprinting (MRF) framework to quantify T₂ ^* in addition to δ f , T₁ , and T₂ tissue properties.

METHODS

The proposed quadratic RF phase-based MRF method (qRF-MRF) combined a varied RF excitation phase with the existing balanced SSFP (bSSFP)-based MRF method to generate signals that were uniquely sensitive to δ f , T₁ , T₂ , as well as the distribution width of intravoxel frequency dispersion, Γ . A dictionary, generated through Bloch simulation, containing possible signal evolutions within the physiological range of δ f , T₁ , T₂ , and Γ , was used to perform parameter estimation. The estimated T₂ and Γ were subsequently used to estimate T₂ ^* . The proposed method was evaluated in phantom experiments and healthy volunteers (N = 5).

RESULTS

The T₁ and T₂ values from the phantom by qRF-MRF demonstrated good agreement with values obtained by traditional gold standard methods (r² = 0.995 and 0.997, respectively; concordance correlation coefficient = 0.978 and 0.995, respectively). The T₂ ^* values from the phantom demonstrated good agreement with values obtained through the multi-echo gradient-echo method (r² = 0.972, concordance correlation coefficient = 0.983). In vivo qRF-MRF-measured T₁ , T₂ , and T₂ ^* values were compared with measurements by existing methods and literature values.

CONCLUSION

The proposed qRF-MRF method demonstrated the potential for simultaneous quantification of δ f , T₁ , T₂ , and T₂ ^* tissue properties.

Magnetic resonance fingerprinting: a technical review

Multiparametric quantitative imaging is gaining increasing interest due to its widespread advantages in clinical applications. Magnetic resonance fingerprinting is a recently introduced approach of fast multiparametric quantitative imaging. In this article, magnetic resonance fingerprinting acquisition, dictionary generation, reconstruction, and validation are reviewed.

Multi-parametric T2 * magnetic resonance fingerprinting using variable echo times

The use of quantitative imaging biomarkers in the imaging of various disease states, including cancer and neurodegenerative disease, has increased in recent years. T₁ , T₂ , and T₂ * relaxation time constants have been shown to be affected by tissue structure or contrast infusion. Acquiring these biomarkers simultaneously in a multi-parametric acquisition could provide more robust detection of tissue changes in various disease states including neurodegeneration and cancer. Traditional magnetic resonance fingerprinting (MRF) has been shown to provide quick, quantitative mapping of T₁ and T₂ relaxation time constants. In this study, T₂ * relaxation is added to the MRF framework using variable echo times (TE). To demonstrate the feasibility of the method and compare incremental and golden angle spiral rotations, simulated phantom data was fit using the proposed method. Additionally, T₁ /T₂ /T₂ */δf MRF as well as conventional T₁ , T₂ , and T₂ * acquisitions were acquired in agar phantoms and the brains of three healthy volunteers. Golden angle spiral rotation was found to reduce inaccuracy resulting from off resonance effects. Strong correlations were found between conventional and MRF values in the T₁ , T₂ , and T₂ * relaxation time constants of the agar phantoms and healthy volunteers. In this study, T₂ * relaxation has been incorporated into the MRF framework by using variable echo times, while still fitting for T₁ and T₂ relaxation time constants. In addition to fitting these relaxation time constants, a novel method for fitting and correcting off resonance effects has been developed.

MR fingerprinting Deep RecOnstruction NEtwork (DRONE)

PURPOSE

Demonstrate a novel fast method for reconstruction of multi-dimensional MR Fingerprinting (MRF) data using Deep Learning methods.

METHODS

A neural network (NN) is defined using the TensorFlow framework and trained on simulated MRF data computed with the Extended Phase Graph formalism. The NN reconstruction accuracy for noiseless and noisy data is compared to conventional MRF template matching as a function of training data size, and quantified in simulated numerical brain phantom data and ISMRM/NIST phantom data measured on 1.5T and 3T scanners with an optimized MRF EPI and MRF FISP sequences with spiral readout. The utility of the method is demonstrated in a healthy subject in vivo at 1.5 T.

RESULTS

Network training required 10 to 74 minutes and once trained, data reconstruction required approximately 10 ms for the MRF EPI and 76 ms for the MRF FISP sequence. Reconstruction of simulated, noiseless brain data using the NN resulted in a root-mean-square error (RMSE) of 2.6 ms for T1 and 1.9 ms for T2. The reconstruction error in the presence of noise was less than 10% for both T1 and T2 for signal-to-noise greater than 25 dB. Phantom measurements yielded good agreement (R2=0.99/0.99 for MRF EPI T1/T2 and 0.94/0.98 for MRF FISP T1/T2) between the T1 and T2 estimated by the NN and reference values from the ISMRM/NIST phantom.

CONCLUSION

Reconstruction of MRF data with a NN is accurate, 300–5000 fold faster and more robust to noise and undersampling than conventional MRF dictionary matching.

Accelerated magnetic resonance fingerprinting using soft-weighted key-hole (MRF-SOHO)

Object

To develop a novel approach for highly accelerated Magnetic Resonance Fingerprinting (MRF) acquisition.

Materials and methods

The proposed method combines parallel imaging, soft-gating and key-hole approaches to highly accelerate MRF acquisition. Slowly varying flip angles (FA), commonly used during MRF acquisition, lead to a smooth change in the signal contrast of consecutive time-point images. This assumption enables sharing of high frequency data between different time-points, similar to what is done in some dynamic MR imaging methods such as key-hole. The proposed approach exploits this information using a SOft-weighted key-HOle (MRF-SOHO) reconstruction to achieve high acceleration factors and/or increased resolution without compromising image quality or increasing scan time. MRF-SOHO was validated on a standard T₁/T₂ phantom and in in-vivo brain acquisitions reconstructing T₁, T₂ and proton density parametric maps.

Results

Accelerated MRF-SOHO using less data per time-point and less time-point images enabled a considerable reduction in scan time (up to 4.6x), while obtaining similar T₁ and T₂ accuracy and precision when compared to zero-filled MRF reconstruction. For the same number of spokes and time-points, the proposed method yielded an enhanced performance in quantifying parameters than the zero-filled MRF reconstruction, which was verified with 2, 1 and 0.7 (sub-millimetre) resolutions.

Conclusion

The proposed MRF-SOHO enabled a 4.6x scan time reduction for an in-plane spatial resolution of 2x2 mm² when compared to zero-filled MRF and enabled sub-millimetric (0.7x0.7 mm²) resolution MRF.

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.

The Combined Quantification and Interpretation of Multiple Quantitative Magnetic Resonance Imaging Metrics Enlightens Longitudinal Changes Compatible with Brain Repair in Relapsing-Remitting Multiple Sclerosis Patients

Objective

Quantitative and semi-quantitative MRI (qMRI) metrics provide complementary specificity and differential sensitivity to pathological brain changes compatible with brain inflammation, degeneration, and repair. Moreover, advanced magnetic resonance imaging (MRI) metrics with overlapping elements amplify the true tissue-related information and limit measurement noise. In this work, we combined multiple advanced MRI parameters to assess focal and diffuse brain changes over 2 years in a group of early-stage relapsing-remitting MS patients.

Methods

Thirty relapsing-remitting MS patients with less than 5 years disease duration and nine healthy subjects underwent 3T MRI at baseline and after 2 years including T1, T2, T2* relaxometry, and magnetization transfer imaging. To assess longitudinal changes in normal-appearing (NA) tissue and lesions, we used analyses of variance and Bonferroni correction for multiple comparisons. Multivariate linear regression was used to assess the correlation between clinical outcome and multiparametric MRI changes in lesions and NA tissue.

Results

In patients, we measured a significant longitudinal decrease of mean T2 relaxation times in NA white matter (p = 0.005) and a decrease of T1 relaxation times in the pallidum (p < 0.05), which are compatible with edema reabsorption and/or iron deposition. No longitudinal changes in qMRI metrics were observed in controls. In MS lesions, we measured a decrease in T1 relaxation time (p-value < 2.2e−16) and a significant increase in MTR (p-value < 1e−6), suggesting repair mechanisms, such as remyelination, increased axonal density, and/or a gliosis. Last, the evolution of advanced MRI metrics—and not changes in lesions or brain volume—were correlated to motor and cognitive tests scores evolution (Adj-R² > 0.4, p < 0.05). In summary, the combination of multiple advanced MRI provided evidence of changes compatible with focal and diffuse brain repair at early MS stages as suggested by histopathological studies.

Magnetic Resonance Fingerprinting-An Overview

Magnetic Resonance Fingerprinting (MRF) is a new approach to quantitative magnetic resonance imaging that allows simultaneous measurement of multiple tissue properties in a single, time-efficient acquisition. The ability to reproducibly and quantitatively measure tissue properties could enable more objective tissue diagnosis, comparisons of scans acquired at different locations and time points, longitudinal follow-up of individual patients and development of imaging biomarkers. This review provides a general overview of MRF technology, current preclinical and clinical applications and potential future directions. MRF has been initially evaluated in brain, prostate, liver, cardiac, musculoskeletal imaging, and measurement of perfusion and microvascular properties through MR vascular fingerprinting.

Regularly incremented phase encoding - MR fingerprinting (RIPE-MRF) for enhanced motion artifact suppression in preclinical cartesian MR fingerprinting

PURPOSE

The regularly incremented phase encoding-magnetic resonance fingerprinting (RIPE-MRF) method is introduced to limit the sensitivity of preclinical MRF assessments to pulsatile and respiratory motion artifacts.

METHODS

As compared to previously reported standard Cartesian-MRF methods (SC-MRF), the proposed RIPE-MRF method uses a modified Cartesian trajectory that varies the acquired phase-encoding line within each dynamic MRF dataset. Phantoms and mice were scanned without gating or triggering on a 7T preclinical MRI scanner using the RIPE-MRF and SC-MRF methods. In vitro phantom longitudinal relaxation time (T₁ ) and transverse relaxation time (T₂ ) measurements, as well as in vivo liver assessments of artifact-to-noise ratio (ANR) and MRF-based T₁ and T₂ mean and standard deviation, were compared between the two methods (n = 5).

RESULTS

RIPE-MRF showed significant ANR reductions in regions of pulsatility (P < 0.005) and respiratory motion (P < 0.0005). RIPE-MRF also exhibited improved precision in T₁ and T₂ measurements in comparison to the SC-MRF method (P <  0.05). The RIPE-MRF and SC-MRF methods displayed similar mean T₁ and T₂ estimates (difference in mean values < 10%).

CONCLUSION

These results show that the RIPE-MRF method can provide effective motion artifact suppression with minimal impact on T₁ and T₂ accuracy for in vivo small animal MRI studies. Magn Reson Med 79:2176-2182, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Spiral MR fingerprinting at 7T with simultaneous B1 estimation

Magnetic resonance fingerprinting is an efficient, new approach for quantitative imaging with MR. We aimed to extend this technique to cases with B1+ inhomogeneities within the imaging volume. Previous approaches have used abrupt changes in flip angles to estimate the B1+ field simultaneously with T1 and T2, using a Cartesian approach in a small-animal scanner at 4.7T. Here, we evaluated whether a similar approach would be suitable for imaging human brains using spiral readouts with a 7T scanner. We found that our modified scheme could significantly reduce the adverse effects of B1+ inhomogeneities even in extreme cases, reducing both the bias and the variance in T2 estimations by an order of magnitude when compared to literature methods. Acquisitions used less than 1.5W/kg SAR and could be performed in 12s per slice. In conclusion, our approach can be used to perform quantitative imaging of the brain at 7T in a short time, simultaneously estimating the B1+ profile.