Magnetic resonance fingerprinting: a technical review

Bringing MRI to low- and middle-income countries: Directions, challenges and potential solutions

The global disparity of magnetic resonance imaging (MRI) is a major challenge, with many low- and middle-income countries (LMICs) experiencing limited access to MRI. The reasons for limited access are technological, economic and social. With the advancement of MRI technology, we explore why these challenges still prevail, highlighting the importance of MRI as the epidemiology of disease changes in LMICs. In this paper, we establish a framework to develop MRI with these challenges in mind and discuss the different aspects of MRI development, including maximising image quality using cost-effective components, integrating local technology and infrastructure and implementing sustainable practices. We also highlight the current solutions-including teleradiology, artificial intelligence and doctor and patient education strategies-and how these might be further improved to achieve greater access to MRI.

Magnetic resonance relaxometry in quantitative imaging of brain gliomas: A literature review

Magnetic resonance (MR) relaxometry is a quantitative imaging method that measures tissue relaxation properties. This review discusses the state of the art of clinical proton MR relaxometry for glial brain tumors. Current MR relaxometry technology also includes MR fingerprinting and synthetic MRI, which solve the inefficiencies and challenges of earlier techniques. Despite mixed results regarding its capability for brain tumor differential diagnosis, there is growing evidence that MR relaxometry can differentiate between gliomas and metastases and between glioma grades. Studies of the peritumoral zones have demonstrated their heterogeneity and possible directions of tumor infiltration. In addition, relaxometry offers T2* mapping that can define areas of tissue hypoxia not discriminated by perfusion assessment. Studies of tumor therapy response have demonstrated an association between survival and progression terms and dynamics of native and contrast-enhanced tumor relaxometric profiles. In conclusion, MR relaxometry is a promising technique for glial tumor diagnosis, particularly in association with neuropathological studies and other imaging techniques.

Emerging Trends in Magnetic Resonance Fingerprinting for Quantitative Biomedical Imaging Applications: A Review

Magnetic resonance imaging (MRI) stands as a vital medical imaging technique, renowned for its ability to offer high-resolution images of the human body with remarkable soft-tissue contrast. This enables healthcare professionals to gain valuable insights into various aspects of the human body, including morphology, structural integrity, and physiological processes. Quantitative imaging provides compositional measurements of the human body, but, currently, either it takes a long scan time or is limited to low spatial resolutions. Undersampled k-space data acquisitions have significantly helped to reduce MRI scan time, while compressed sensing (CS) and deep learning (DL) reconstructions have mitigated the associated undersampling artifacts. Alternatively, magnetic resonance fingerprinting (MRF) provides an efficient and versatile framework to acquire and quantify multiple tissue properties simultaneously from a single fast MRI scan. The MRF framework involves four key aspects: (1) pulse sequence design; (2) rapid (undersampled) data acquisition; (3) encoding of tissue properties in MR signal evolutions or fingerprints; and (4) simultaneous recovery of multiple quantitative spatial maps. This paper provides an extensive literature review of the MRF framework, addressing the trends associated with these four key aspects. There are specific challenges in MRF for all ranges of magnetic field strengths and all body parts, which can present opportunities for further investigation. We aim to review the best practices in each key aspect of MRF, as well as for different applications, such as cardiac, brain, and musculoskeletal imaging, among others. A comprehensive review of these applications will enable us to assess future trends and their implications for the translation of MRF into these biomedical imaging applications.

Reconstruction of multi-phase parametric maps in 4D-magnetic resonance fingerprinting (4D-MRF) by optimization of local T1 and T2 sensitivities

BACKGROUND

Time-resolved magnetic resonance fingerprinting (MRF), or 4D-MRF, has been demonstrated its feasibility in motion management in radiotherapy (RT). However, the prohibitive long acquisition time is one of challenges of the clinical implementation of 4D-MRF. The shortening of acquisition time causes data insufficiency in each respiratory phase, leading to poor accuracies and consistencies of the predicted tissues' properties of each phase.

PURPOSE

To develop a technique for the reconstruction of multi-phase parametric maps in four-dimensional magnetic resonance fingerprinting (4D-MRF) through the optimization of local T1 and T2 sensitivities.

METHODS

The proposed technique employed an iterative optimization to tailor the data arrangement of each phase by manipulation of inter-phase frames, such that the T1 and T2 sensitivities, which were quantified by the modified Minkowski distance, of the truncated signal evolution curve was maximized. The multi-phase signal evolution curves were modified by sliding window reconstruction and inter-phase frame sharing (SWIFS). Motion correction (MC) and dot product matching were sequentially performed on the modified signal evolution and dictionary to reconstruct the multi-parametric maps. The proposed technique was evaluated by numerical simulations using the extended cardiac-torso (XCAT) phantom with regular and irregular breathing patterns, and by in vivo MRF data of three health volunteers and six liver cancer patients acquired at a 3.0 T scanner.

RESULTS

In simulation study, the proposed SWIFS approach achieved the overall mean absolute percentage error (MAPE) of 8.62% ± 1.59% and 16.2% ± 3.88% for the eight-phases T1 and T2 maps, respectively, in the sagittal view with irregular breathing patterns. In contrast, the overall MAPE of T1 and T2 maps generated by the conventional approach with multiple MRF repetitions were 22.1% ± 11.0% and 30.8% ± 14.9%, respectively. For in-vivo study, the predicted mean T1 and T2 of liver by the proposed SWIFS approach were 795 ms ± 38.9 ms and 58.3 ms ± 11.7 ms, respectively.

CONCLUSIONS

Both simulation and in vivo results showed that the approach empowered by T1 and T2 sensitivities optimization and sliding window under the shortened acquisition of MRF had superior performance in the estimation of multi-phase T1 and T2 maps as compared to the conventional approach with oversampling of MRF data.

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.

Free-breathing simultaneous native myocardial T1, T2 and T1ρ mapping with Cartesian acquisition and dictionary matching

BACKGROUND

T1, T2 and T1ρ are well-recognized parameters for quantitative cardiac MRI. Simultaneous estimation of these parameters allows for comprehensive myocardial tissue characterization, such as myocardial fibrosis and edema. However, conventional techniques either quantify the parameters individually with separate breath-hold acquisitions, which may result in unregistered parameter maps, or estimate multiple parameters in a prolonged breath-hold acquisition, which may be intolerable to patients. We propose a free-breathing multi-parametric mapping (FB-MultiMap) technique that provides co-registered myocardial T1, T2 and T1ρ maps in a single efficient acquisition.

METHODS

The proposed FB-MultiMap performs electrocardiogram-triggered single-shot Cartesian acquisition over 16 consecutive cardiac cycles, where inversion, T2 and T1ρ preparations are introduced for varying contrasts. A diaphragmatic navigator was used for prospective through-plane motion correction and the in-plane motion was corrected retrospectively with a group-wise image registration method. Quantitative mapping was conducted through dictionary matching of the motion corrected images, where the subject-specific dictionary was created using Bloch simulations for a range of T1, T2 and T1ρ values, as well as B1 factors to account for B1 inhomogeneities. The FB-MultiMap was optimized and validated in numerical simulations, phantom experiments, and in vivo imaging of 15 healthy subjects and six patients with suspected cardiac diseases.

RESULTS

The phantom T1, T2 and T1ρ values estimated with FB-MultiMap agreed well with reference measurements with no dependency on heart rate. In healthy subjects, FB-MultiMap T1 was higher than MOLLI T1 mapping (1218 ± 50 ms vs. 1166 ± 38 ms, p < 0.001). The myocardial T2 and T1ρ estimated with FB-MultiMap were lower compared to the mapping with T2- or T1ρ-prepared 2D balanced steady-state free precession (T2: 41.2 ± 2.8 ms vs. 42.5 ± 3.1 ms, p = 0.06; T1ρ: 45.3 ± 4.4 ms vs. 50.2 ± 4.0, p < 0.001). The pathological changes in myocardial parameters measured with FB-MultiMap were consistent with conventional techniques in all patients.

CONCLUSION

The proposed free-breathing multi-parametric mapping technique provides co-registered myocardial T1, T2 and T1ρ maps in 16 heartbeats, achieving similar mapping quality to conventional breath-hold mapping methods.

Magnetic Resonance Fingerprinting: A Review of Clinical Applications

ABSTRACT

Magnetic resonance fingerprinting (MRF) is an approach to quantitative magnetic resonance imaging that allows for efficient simultaneous measurements of multiple tissue properties, which are then used to create accurate and reproducible quantitative maps of these properties. As the technique has gained popularity, the extent of preclinical and clinical applications has vastly increased. The goal of this review is to provide an overview of currently investigated preclinical and clinical applications of MRF, as well as future directions. Topics covered include MRF in neuroimaging, neurovascular, prostate, liver, kidney, breast, abdominal quantitative imaging, cardiac, and musculoskeletal applications.

Tailored magnetic resonance fingerprinting

Neuroimaging of certain pathologies requires both multi-parametric qualitative and quantitative imaging. The role of the quantitative MRI (qMRI) is well accepted but suffers from long acquisition times leading to patient discomfort, especially in geriatric and pediatric patients. Previous studies show that synthetic MRI can be used in order to reduce the scan time and provide qMRI as well as multi-contrast data. However, this approach suffers from artifacts such as partial volume and flow. In order to increase the scan efficiency (the number of contrasts and quantitative maps acquired per unit time), we designed, simulated, and demonstrated rapid, simultaneous, multi-contrast qualitative (T₁ weighted, T₁ fluid attenuated inversion recovery (FLAIR), T₂ weighted, water, and fat), and quantitative imaging (T₁ and T₂ maps) through the approach of tailored MR fingerprinting (TMRF) to cover whole-brain in approximately four minutes. We performed TMRF on in vivo four healthy human brains and in vitro ISMRM/NIST phantom and compared with vendor supplied gold standard (GS) and MRF sequences. All scans were performed on a 3 T GE Premier system and images were reconstructed offline using MATLAB. The reconstructed qualitative images were then subjected to custom DL denoising and gradient anisotropic diffusion denoising. The quantitative tissue parametric maps were reconstructed using a dense neural network to gain computational speed compared to dictionary matching. The grey matter and white matter tissues in qualitative and quantitative data for the in vivo datasets were segmented semi-automatically. The SNR and mean contrasts were plotted and compared across all three methods. The GS images show better SNR in all four subjects compared to MRF and TMRF (GS > TMRF>MRF). The T₁ and T₂ values of MRF are relatively overestimated as compared to GS and TMRF. The scan efficiency for TMRF is 1.72 min^-1 which is higher compared to GS (0.32 min^-1) and MRF (0.90 min^-1).

Adaptive model-based Magnetic Resonance

PURPOSE

Conventional sequences are static in nature, fixing measurement parameters in advance in anticipation of a wide range of expected tissue parameter values. We set out to design and benchmark a new, personalized approach-termed adaptive MR-in which incoming subject data is used to update and fine-tune the pulse sequence parameters in real time.

METHODS

We implemented an adaptive, real-time multi-echo (MTE) experiment for estimating T₂ s. Our approach combined a Bayesian framework with model-based reconstruction. It maintained and continuously updated a prior distribution of the desired tissue parameters, including T₂ , which was used to guide the selection of sequence parameters in real time.

RESULTS

Computer simulations predicted accelerations between 1.7- and 3.3-fold for adaptive multi-echo sequences relative to static ones. These predictions were corroborated in phantom experiments. In healthy volunteers, our adaptive framework accelerated the measurement of T₂ for n-acetyl-aspartate by a factor of 2.5.

CONCLUSION

Adaptive pulse sequences that alter their excitations in real time could provide substantial reductions in acquisition times. Given the generality of our proposed framework, our results motivate further research into other adaptive model-based approaches to MRI and MRS.

Accuracy, repeatability, and reproducibility of T1 and T2 relaxation times measurement by 3D magnetic resonance fingerprinting with different dictionary resolutions

OBJECTIVES

To assess the accuracy, repeatability, and reproducibility of T₁ and T₂ relaxation time measurements by three-dimensional magnetic resonance fingerprinting (3D MRF) using various dictionary resolutions.

METHODS

The ISMRM/NIST phantom was scanned daily for 10 days in two 3 T MR scanners using a 3D MRF sequence reconstructed using four dictionaries with varying step sizes and one dictionary with wider ranges. Thirty-nine healthy volunteers were enrolled: 20 subjects underwent whole-brain MRF scans in both scanners and the rest in one scanner. ROI/VOI analyses were performed on phantom and brain MRF maps. Accuracy, repeatability, and reproducibility metrics were calculated.

RESULTS

In the phantom study, all dictionaries showed high T₁ linearity to the reference values (R² > 0.99), repeatability (CV < 3%), and reproducibility (CV < 3%) with lower linearity (R² > 0.98), repeatability (CV < 6%), and reproducibility (CV ≤ 4%) for T₂ measurement. The volunteer study demonstrated high T₁ reproducibility of within-subject CV (wCV) < 4% by all dictionaries with the same ranges, both in the brain parenchyma and CSF. Yet, reproducibility was moderate for T₂ measurement (wCV < 8%). In CSF measurement, dictionaries with a smaller range showed a seemingly better reproducibility (T₁, wCV 3%; T₂, wCV 8%) than the much wider range dictionary (T₁, wCV 5%; T₂, wCV 13%). Truncated CSF relaxometry values were evident in smaller range dictionaries.

CONCLUSIONS

The accuracy, repeatability, and reproducibility of 3D MRF across various dictionary resolutions were high for T₁ and moderate for T₂ measurements. A lower-resolution dictionary with a well-defined range may be adequate, thus significantly reducing the computational load.

KEY POINTS

• A lower-resolution dictionary with a well-defined range may be sufficient for 3D MRF reconstruction. • CSF relaxation times might be underestimated due to truncation by the upper dictionary range. • Dictionary with a higher upper range might be advisable, especially for CSF evaluation and elderly subjects whose perivascular spaces are more prominent.

GPU accelerated grouped magnetic resonance fingerprinting using clustering techniques

Magnetic Resonance Fingerprinting (MRF) is a new quantitative technique of Magnetic Resonance Imaging (MRI). Conventionally, MRF requires sequential correlation of the acquired MRF signals with all the signals of (a large sized) MRF dictionary. This is a computationally intensive matching process and is a major challenge in MRF image reconstruction. This paper introduces the use of clustering techniques (to reduce the effective size of MRF dictionary) by splitting MRF dictionary into multiple small sized MRF dictionary components called MRF signal groups. The proposed method has been further optimized for parallel processing to reduce the computation time of MRF pattern matching. A multi-core GPU based parallel framework has been developed that enables the MRF algorithm to process multiple MRF signals simultaneously. Experiments have been performed on human head and phantom datasets. The results show that the proposed method accelerates the conventional MRF (MATLAB based) reconstruction time up to 25× with single-core CPU implementation, 300× with multi- core CPU implementation and 1035× with the proposed multi-core GPU based framework by keeping the SNR of the resulting images in a clinically acceptable range. Furthermore, experimental results show that the memory requirements of MRF dictionary get significantly reduced (due to efficient memory utilization) in the proposed method.

High-Resolution 3D Magnetic Resonance Fingerprinting With a Graph Convolutional Network

Magnetic resonance fingerprinting (MRF) is a novel quantitative imaging framework for rapid and simultaneous quantification of multiple tissue properties. 3D MRF allows higher through-plane resolution, but the acquisition process is slow when whole-brain coverage is needed. Existing methods for acceleration mainly rely on GRAPPA for k-space interpolation in the partition-encoding direction, limiting the acceleration factor to 2 or 3. In this work, we replace GRAPPA with a deep learning approach for accurate tissue quantification with greater acceleration. Specifically, a graph convolution network (GCN) is developed to cater to the non-Cartesian spiral sampling trajectories typical in MRF acquisition. The GCN maintains high quantification accuracy with up to 6-fold acceleration and allows 1mm isotropic resolution whole-brain 3D MRF data to be acquired in 3min and submillimeter 3D MRF (0.8mm) in 5min, greatly improving the feasibility of MRF in clinical settings.

Approximate B 1 + scaling of the SSFP steady state

PURPOSE

It is shown that the steady state of rapid, TR-periodic steady-state free precession (SSFP) sequences at small to moderate flip angles exhibits a universal, approximate scaling law with respect to variations of B 1 + $$ {B}_1^{+} $$ . Implications for the accuracy and precision of relaxometry experiments are discussed.

METHODS

The approximate scaling law is derived from and numerically tested against known analytical solutions. To assess the attainable estimator precision in a typical relaxometry experiment, we calculate the Cramér-Rao bound (CRB) and perform Monte Carlo (MC) simulations.

RESULTS

The approximate universal scaling holds well up to moderate flip angles. For pure steady state relaxometry, we observe a significant precision penalty for simultaneous estimation of R 1 $$ {R}_1 $$ and B 1 + $$ {B}_1^{+} $$ , whereas good R 2 $$ {R}_2 $$ estimates can be obtained without even knowing the correct actual flip angle.

CONCLUSION

Simultaneous estimation of R 1 $$ {R}_1 $$ and B 1 + $$ {B}_1^{+} $$ from a set of SSFP steady states alone is not advisable. Apart from separate B 1 + $$ {B}_1^{+} $$ measurements, the problem can be addressed by adding transient state information, but, depending on the situation, residual effects due to the scaling may still require some attention.

Optimisation of data acquisition towards continuous cardiac Magnetic Resonance Fingerprinting applications

PURPOSE

Assess and optimise acquisition parameters for continuous cardiac Magnetic Resonance Fingerprinting (MRF).

METHODS

Different acquisition schemes (flip angle amplitude, lobe size, T2-preparation pulses) for cardiac MRF were assessed in simulations and phantom and demonstrated in one healthy volunteer. Three different experimental designs were evaluated using central composite and fractional factorial designs. Relative errors for T1 and T2 were calculated for a wide range of realistic T1 and T2 value combinations. The effect of different designs on the accuracy of T1 and T2 was assessed using response surface modelling and Cohen's f calculations.

RESULTS

Larger flip angle amplitudes lead to an improvement of T2 accuracy and precision for simulations and phantom experiments. Similar effects could also be shown qualitatively in in-vivo scans. Accuracy and precision of T1 were robust to different design parameters with improved values for faster flip angle variation. Cohen's f showed that T2-preparation pulses influence the accuracy of T2. The number of pulses used is the most important parameter. Without T2-preparation pulses, RMSE were 3.0 ± 8.09 % for T1 and 16.24 ± 14.47 % for T2. Using those pulses reduced the RMSE to 2.3 ± 8.4 % for T1 and 14.11 ± 13.46 % for T2. Nonetheless, even if the improvement is significant, RMSE are still too high for reliable quantification.

CONCLUSION

In contrast to previous study using triggered MRF sequences using < 30° flip angles, large flip angle amplitudes led to better results for continuous cardiac MRF sequences. T2-preparation pulse can improve the accuracy of T2 estimation but lead to longer scan times.

Liver T1 Relaxation Quantification Using a 3-Dimensional Interleaved Look-Locker Acquisition With T2 Preparation Pulse Sequence (3D-QALAS): Comparison With Conventional 2-Dimensional MOLLI

BACKGROUND

Changes in liver magnetic resonance imaging T1 relaxation times are associated with histologic inflammation and fibrosis.

OBJECTIVE

To compare liver T1 measurements obtained using a novel single-breath-hold 3-dimensional (3D) whole-liver T1 estimation method (3D-QALAS) to standard-of-care 2-dimensional (2D) modified Look-Locker (2D-MOLLI) measurements.

METHODS

With institutional review board approval, research magnetic resonance imaging examinations were performed in 19 participants at 1.5 T. T1 relaxometry of the liver was performed using a novel 3D whole-liver T1 estimation method (3D-QALAS) as well as a 2D modified Look-Locker (2D-MOLLI) method. The 3D method covered the entire liver in a single breath hold, whereas 2D imaging was performed at 4 anatomic levels in 4 consecutive breath holds. T1 measurements from parametric maps were obtained by a single operator, and region-of-interest area-weighted mean T1 values were calculated. Pearson correlation ( r ) was used to assess correlation between T1 estimation methods, and the paired t test and Bland-Altman analysis were used to compare agreement in T1 measurements.

RESULTS

In 18 participants (1 participant was excluded from analysis because of respiratory motion artifacts on 3D-QALAS images), 2D-MOLLI and 3D-QALAS mean T1 measurements were strongly correlated ( r = 0.95, [95% CI: 0.87-0.98]; P < 0.0001). 2D-MOLLI T1 values were significantly longer than 3D-QALAS values (647.2 ± 87.3 milliseconds vs. 554.7 ± 75.8 milliseconds; P < 0.0001) with mean bias = 92.5 milliseconds (95% limits of agreement, 36.8, 148.2 milliseconds).

CONCLUSION

Whole-liver T1 measurements obtained using a novel single-breath-hold 3D T1 estimation method correlate with a standard-of-care multiple consecutive-breath-hold 2D single-slice method but demonstrate systematic bias that should be considered or corrected when used in a clinical or research setting.

Deep-Learning-Based Contrast Synthesis From MRF Parameter Maps in the Knee Joint

BACKGROUND

Magnetic resonance fingerprinting (MRF) is a method to speed up acquisition of quantitative MRI data. However, MRF does not usually produce contrast-weighted images that are required by radiologists, limiting reachable total scan time improvement. Contrast synthesis from MRF could significantly decrease the imaging time.

PURPOSE

To improve clinical utility of MRF by synthesizing contrast-weighted MR images from the quantitative data provided by MRF, using U-nets that were trained for the synthesis task utilizing L1- and perceptual loss functions, and their combinations.

STUDY TYPE

Retrospective.

POPULATION

Knee joint MRI data from 184 subjects from Northern Finland 1986 Birth Cohort (ages 33-35, gender distribution not available).

FIELD STRENGTH AND SEQUENCE

A 3 T, multislice-MRF, proton density (PD)-weighted 3D-SPACE (sampling perfection with application optimized contrasts using different flip angle evolution), fat-saturated T2-weighted 3D-space, water-excited double echo steady state (DESS).

ASSESSMENT

Data were divided into training, validation, test, and radiologist's assessment sets in the following way: 136 subjects to training, 3 for validation, 3 for testing, and 42 for radiologist's assessment. The synthetic and target images were evaluated using 5-point Likert scale by two musculoskeletal radiologists blinded and with quantitative error metrics.

STATISTICAL TESTS

Friedman's test accompanied with post hoc Wilcoxon signed-rank test and intraclass correlation coefficient. The statistical cutoff P <0.05 adjusted by Bonferroni correction as necessary was utilized.

RESULTS

The networks trained in the study could synthesize conventional images with high image quality (Likert scores 3-4 on a 5-point scale). Qualitatively, the best synthetic images were produced with combination of L1- and perceptual loss functions and perceptual loss alone, while L1-loss alone led to significantly poorer image quality (Likert scores below 3). The interreader and intrareader agreement were high (0.80 and 0.92, respectively) and significant. However, quantitative image quality metrics indicated best performance for the pure L1-loss.

DATA CONCLUSION

Synthesizing high-quality contrast-weighted images from MRF data using deep learning is feasible. However, more studies are needed to validate the diagnostic accuracy of these synthetic images.

EVIDENCE LEVEL

4.

TECHNICAL EFFICACY

Stage 1.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

[Applications of Artificial Intelligence in MR Image Acquisition and Reconstruction]

Recently, artificial intelligence (AI) technology has shown potential clinical utility in a wide range of MRI fields. In particular, AI models for improving the efficiency of the image acquisition process and the quality of reconstructed images are being actively developed by the MR research community. AI is expected to further reduce acquisition times in various MRI protocols used in clinical practice when compared to current parallel imaging techniques. Additionally, AI can help with tasks such as planning, parameter optimization, artifact reduction, and quality assessment. Furthermore, AI is being actively applied to automate MR image analysis such as image registration, segmentation, and object detection. For this reason, it is important to consider the effects of protocols or devices in MR image analysis. In this review article, we briefly introduced issues related to AI application of MR image acquisition and reconstruction.

High-Quality MR Fingerprinting Reconstruction Using Structured Low-Rank Matrix Completion and Subspace Projection

Due to the capability of fast multiparametric quantitative imaging, magnetic resonance fingerprinting (MRF) is becoming a promising quantitative magnetic resonance imaging approach. However, the artifacts caused by the highly undersampled data acquisition lead to inaccurate estimation of the tissue parameter maps. Based on the assumption that the 3-D MRF data can be modeled as a piecewise smooth signal, with the discontinuities localized to the zero sets of a bandlimited function, we exploit the low-rank property of the structured Toeplitz matrix constructed from the Fourier measurements. In addition, we adopt the subspace projection scheme to improve the accuracy of parameter estimation. In order to efficiently solve the regularized problem, we propose an iterative two-stage algorithm, which alternately updates the k -space data and projects the space-time matrix into the dictionary space. Numerical experiments demonstrate that the proposed algorithm shows significant improvement in MRF time-series images reconstruction and can provide more accurate parameter maps over the state-of-the-art algorithms.

Circumventing the curse of dimensionality in magnetic resonance fingerprinting through a deep learning approach

Magnetic resonance fingerprinting (MRF) is a rapidly developing approach for fast quantitative MRI. A typical drawback of dictionary-based MRF is an explosion of the dictionary size as a function of the number of reconstructed parameters, according to the "curse of dimensionality", which determines an explosion of resource requirements. Neural networks (NNs) have been proposed as a feasible alternative, but this approach is still in its infancy. In this work, we design a deep learning approach to MRF using a fully connected network (FCN). In the first part we investigate, by means of simulations, how the NN performance scales with the number of parameters to be retrieved in comparison with the standard dictionary approach. Four MRF sequences were considered: IR-FISP, bSSFP, IR-FISP-B₁ , and IR-bSSFP-B₁ , the latter two designed to be more specific for B1+ parameter encoding. Estimation accuracy, memory usage, and computational time required to perform the estimation task were considered to compare the scalability capabilities of the dictionary-based and the NN approaches. In the second part we study optimal training procedures by including different data augmentation and preprocessing strategies during training to achieve better accuracy and robustness to noise and undersampling artifacts. The study is conducted using the IR-FISP MRF sequence exploiting both simulations and in vivo acquisitions. Results demonstrate that the NN approach outperforms the dictionary-based approach in terms of scalability capabilities. Results also allow us to heuristically determine the optimal training strategy to make an FCN able to predict T₁ ,  T₂ , and M₀ maps that are in good agreement with those obtained with the original dictionary approach. k-SVD denoising is proposed and found to be critical as a preprocessing step to handle undersampled data.

Radiomics with 3-dimensional magnetic resonance fingerprinting: influence of dictionary design on repeatability and reproducibility of radiomic features

Objectives

We aimed to investigate the influence of magnetic resonance fingerprinting (MRF) dictionary design on radiomic features using in vivo human brain scans.

Methods

Scan-rescans of three-dimensional MRF and conventional T1-weighted imaging were performed on 21 healthy volunteers (9 males and 12 females; mean age, 41.3 ± 14.6 years; age range, 22–72 years). Five patients with multiple sclerosis (3 males and 2 females; mean age, 41.2 ± 7.3 years; age range, 32–53 years) were also included. MRF data were reconstructed using various dictionaries with different step sizes. First- and second-order radiomic features were extracted from each dataset. Intra-dictionary repeatability and inter-dictionary reproducibility were evaluated using intraclass correlation coefficients (ICCs). Features with ICCs > 0.90 were considered acceptable. Relative changes were calculated to assess inter-dictionary biases.

Results

The overall scan-rescan ICCs of MRF-based radiomics ranged from 0.86 to 0.95, depending on dictionary step size. No significant differences were observed in the overall scan-rescan repeatability of MRF-based radiomic features and conventional T1-weighted imaging (p = 1.00). Intra-dictionary repeatability was insensitive to dictionary step size differences. MRF-based radiomic features varied among dictionaries (overall ICC for inter-dictionary reproducibility, 0.62–0.99), especially when step sizes were large. First-order and gray level co-occurrence matrix features were the most reproducible feature classes among different step size dictionaries. T1 map-derived radiomic features provided higher repeatability and reproducibility among dictionaries than those obtained with T2 maps.

Conclusion

MRF-based radiomic features are highly repeatable in various dictionary step sizes. Caution is warranted when performing MRF-based radiomics using datasets containing maps generated from different dictionaries.

Key Points

• MRF-based radiomic features are highly repeatable in various dictionary step sizes.

• Use of different MRF dictionaries may result in variable radiomic features, even when the same MRF acquisition data are used.

• Caution is needed when performing radiomic analysis using data reconstructed from different dictionaries.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00330-022-08555-3.

A simultaneous multi-slice T2 mapping framework based on overlapping-echo detachment planar imaging and deep learning reconstruction

PURPOSE

Quantitative MRI (qMRI) is of great importance to clinical medicine and scientific research. However, most qMRI techniques are time-consuming and sensitive to motion, especially when a large 3D volume is imaged. To accelerate the acquisition, a framework is proposed to realize reliable simultaneous multi-slice T₂ mapping.

METHODS

The simultaneous multi-slice T₂ mapping framework is based on overlapping-echo detachment (OLED) planar imaging (dubbed SMS-OLED). Multi-slice overlapping-echo signals were generated by multiple excitation pulses together with echo-shifting gradients. The signals were excited and acquired with a single-channel coil. U-Net was used to reconstruct T₂ maps from the acquired overlapping-echo image.

RESULTS

Single-shot double-slice and two-shot triple-slice SMS-OLED scan schemes were designed according to the framework for evaluation. Simulations, water phantom, and in vivo rat brain experiments were carried out. Overlapping-echo signals were acquired, and T₂ maps were reconstructed and compared with references. The results demonstrate the superior performance of our method.

CONCLUSION

Two slices of T₂ maps can be obtained in a single shot within hundreds of milliseconds. Higher quality multi-slice T₂ maps can be obtained via multiple shots. SMS-OLED provides a lower specific absorption rate scheme compared with conventional SMS methods with a coil with only a single receiver channel. The new method is of potential in dynamic qMRI and functional qMRI where temporal resolution is vital.

An MR fingerprinting approach for quantitative inhomogeneous magnetization transfer imaging

PURPOSE

Magnetization transfer (MT) and inhomogeneous MT (ihMT) contrasts are used in MRI to provide information about macromolecular tissue content. In particular, MT is sensitive to macromolecules, and ihMT appears to be specific to myelinated tissue. This study proposes a technique to characterize MT and ihMT properties from a single acquisition, producing both semiquantitative contrast ratios and quantitative parameter maps.

THEORY AND METHODS

Building on previous work that uses multiband RF pulses to efficiently generate ihMT contrast, we propose a cyclic steady-state approach that cycles between multiband and single-band pulses to boost the achieved contrast. Resultant time-variable signals are reminiscent of an MR fingerprinting acquisition, except that the signal fluctuations are entirely mediated by MT effects. A dictionary-based low-rank inversion method is used to reconstruct the resulting images and to produce both semiquantitative MT ratio and ihMT ratio maps, as well as quantitative parameter estimates corresponding to an ihMT tissue model.

RESULTS

Phantom and in vivo brain data acquired at 1.5 Tesla demonstrate the expected contrast trends, with ihMT ratio maps showing contrast more specific to white matter, as has been reported by others. Quantitative estimation of semisolid fraction and dipolar T₁ was also possible and yielded measurements consistent with literature values in the brain.

CONCLUSION

By cycling between multiband and single-band pulses, an entirely MT-mediated fingerprinting method was demonstrated. This proof-of-concept approach can be used to generate semiquantitative maps and quantitatively estimate some macromolecular-specific tissue parameters.

Accelerated reconstruction of dictionary-based T2 relaxation maps based on dictionary compression and gradient descent search algorithms

Background Quantitative T₂-relaxation-based contrast maps have shown to be highly beneficial for clinical diagnosis and follow-up. The generation of quantitative maps, however, is impaired by long acquisition times, and time-consuming post-processing schemes. The EMC platform is a dictionary-based technique, which involves simulating theoretical signal curves for different physical and experimental values, followed by matching the experimentally acquired signals to the set simulated ones. Purpose Although the EMC technique has shown to produce accurate T₂ maps, it involves computationally intensive post-processing procedures. In this work we present an approach for accelerating the reconstruction of T₂ relaxation maps. Methods This work presents two alternative post-processing approaches for accelerating the reconstruction of EMC-based T₂ relaxation maps. These are (a) Dictionary compression using principal component analysis (PCA) and (b) gradient-descent search algorithm. Additional acceleration was achieved by finding the optimal MATLAB C++ compiler. The utility of the two suggested approaches was examined by calculating the relative error, produced by each technique. Results Gradient descent method was in perfect agreement with the ground truth exhaustive search matching process. PCA based acceleration produced root mean square error (RMSE) of up to 4% compared to exhaustive matching process. Overall acceleration of x16 was achieved using gradient descent in addition to x7 acceleration by choosing the optimal MATLAB C++ compiler. Conclusions Postprocessing of EMC-based T₂ relaxation maps can be accelerated without impairing the accuracy of the ensuing T₂ values.

Data-driven algorithm for myelin water imaging: Probing subvoxel compartmentation based on identification of spatially global tissue features

PURPOSE

Multicomponent analysis of MRI T₂ relaxation time (mcT₂ ) is commonly used for estimating myelin content by separating the signal at each voxel into its underlying distribution of T₂ values. This voxel-based approach is challenging due to the large ambiguity in the multi-T₂ space and the low SNR of MRI signals. Herein, we present a data-driven mcT₂ analysis, which utilizes the statistical strength of identifying spatially global mcT₂ motifs in white matter segments before deconvolving the local signal at each voxel.

METHODS

Deconvolution is done using a tailored optimization scheme, which incorporates the global mcT₂ motifs without additional prior assumptions regarding the number of microscopic components. The end results of this process are voxel-wise myelin water fraction maps.

RESULTS

Validations are shown for computer-generated signals, uniquely designed subvoxel mcT₂ phantoms, and in vivo human brain. Results demonstrated excellent fitting accuracy, both for the numerical and the physical mcT₂ phantoms, exhibiting excellent agreement between calculated myelin water fraction and ground truth. Proof-of-concept in vivo validation is done by calculating myelin water fraction maps for white matter segments of the human brain. Interscan stability of myelin water fraction values was also estimated, showing good correlation between scans.

CONCLUSION

We conclude that studying global tissue motifs prior to performing voxel-wise mcT₂ analysis stabilizes the optimization scheme and efficiently overcomes the ambiguity in the T₂ space. This new approach can improve myelin water imaging and the investigation of microstructural compartmentation in general.

T2 quantification in brain using 3D fast spin-echo imaging with long echo trains

Purpose

Three‐dimensional fast spin‐echo (FSE) sequences commonly use very long echo trains (>64 echoes) and severely reduced refocusing angles. They are increasingly used in brain exams due to high, isotropic resolution and reasonable scan time when using long trains and short interecho spacing. In this study, T₂ quantification in 3D FSE is investigated to achieve increased resolution when comparing with established 2D (proton‐density dual‐echo and multi‐echo spin‐echo) methods.

Methods

The FSE sequence design was explored to use long echo trains while minimizing T₂ fitting error and maintaining typical proton density and T₂‐weighted contrasts. Constant and variable flip angle trains were investigated using extended phase graph and Bloch equation simulations. Optimized parameters were analyzed in phantom experiments and validated in vivo in comparison to 2D methods for eight regions of interest in brain, including deep gray‐matter structures and white‐matter tracts.

Results

Phantom and healthy in vivo brain T₂ measurements showed that optimized variable echo‐train 3D FSE performs similarly to previous 2D methods, while achieving three‐fold‐higher slice resolution, evident visually in the 3D T₂ maps. Optimization resulted in better T₂ fitting and compared well with standard multi‐echo spin echo (within the 8‐ms confidence limits defined based on Bland‐Altman analysis).

Conclusion

T₂ mapping using 3D FSE with long echo trains and variable refocusing angles provides T₂ accuracy in agreement with 2D methods with additional high‐resolution benefits, allowing isotropic views while avoiding incidental magnetization transfer effects. Consequently, optimized 3D sequences should be considered when choosing T₂ mapping methods for high anatomic detail.

A Spatial Off-Resonance Correction in Spirals for Magnetic Resonance Fingerprinting

In MR Fingerprinting (MRF), balanced Steady-State Free Precession (bSSFP) has advantages over unbalanced SSFP because it retains the spin history achieving a higher signal-to-noise ratio (SNR) and scan efficiency. However, bSSFP-MRF is not frequently used because it is sensitive to off-resonance, producing artifacts and blurring, and affecting the parametric map quality. Here we propose a novel Spatial Off-resonance Correction (SOC) approach for reducing these artifacts in bSSFP-MRF with spiral trajectories. SOC-MRF uses each pixel's Point Spread Function to create system matrices that encode both off-resonance and gridding effects. We iteratively compute the inverse of these matrices to reduce the artifacts. We evaluated the proposed method using brain simulations and actual MRF acquisitions of a standardized T1/T2 phantom and five healthy subjects. The results show that the off-resonance distortions in T1/T2 maps were considerably reduced using SOC-MRF. For T2, the Normalized Root Mean Square Error (NRMSE) was reduced from 17.3 to 8.3% (simulations) and from 35.1 to 14.9% (phantom). For T1, the NRMS was reduced from 14.7 to 7.7% (simulations) and from 17.7 to 6.7% (phantom). For in-vivo, the mean and standard deviation in different ROI in white and gray matter were significantly improved. For example, SOC-MRF estimated an average T2 for white matter of 77ms (the ground truth was 74ms) versus 50 ms of MRF. For the same example the standard deviation was reduced from 18 ms to 6ms. The corrections achieved with the proposed SOC-MRF may expand the potential applications of bSSFP-MRF, taking advantage of its better SNR property.

Small-Molecule Prodrug Nanoassemblies: An Emerging Nanoplatform for Anticancer Drug Delivery

The antitumor efficiency and clinical translation of traditional nanomedicines is mainly restricted by low drug loading, complex preparation technology, and potential toxicity caused by the overused carrier materials. In recent decades, small-molecule prodrug nanoassemblies (SMP-NAs), which are formed by the self-assembly of prodrugs themselves, have been widely investigated with distinct advantages of ultrahigh drug-loading and negligible excipients-trigged adverse reaction. Benefited from the simple preparation process, SMP-NAs are widely used for chemotherapy, phototherapy, immunotherapy, and tumor diagnosis. In addition, combination therapy based on the accurate co-delivery behavior of SMP-NAs can effectively address the challenges of tumor heterogeneity and multidrug resistance. Recent trends in SMP-NAs are outlined, and the corresponding self-assembly mechanisms are discussed in detail. Besides, the smart stimuli-responsive SMP-NAs and the combination therapy based on SMP-NAs are summarized, with special emphasis on the structure-function relationships. Finally, the outlooks and potential challenges of SMP-NAs in cancer therapy are highlighted.

Temperature dependence, accuracy, and repeatability of T1 and T2 relaxation times for the ISMRM/NIST system phantom measured using MR fingerprinting

PURPOSE

Before MR fingerprinting (MRF) can be adopted clinically, the derived quantitative values must be proven accurate and repeatable over a range of T₁ and T₂ values and temperatures. Correct assessment of accuracy and precision as well as comparison between measurements can only be performed when temperature is either controlled or corrected for. The purpose of this study was to investigate the temperature dependence of T₁ and T₂ MRF values and evaluate the accuracy and repeatability of temperature-corrected relaxation values derived from a B₁ -corrected MRF-fast imaging with steady-state precession implementation using 2 different dictionary sizes.

METHODS

The International Society of MR in Medicine/National Institute of Standards and Technology phantom was scanned using an MRF sequence of 2 different lengths, a variable flip angle T₁ , and a multi-echo spin echo T₂ at 14 temperatures ranging from 15°C to 28°C and investigated with a linear regression model. Temperature-corrected accuracy was evaluated by correlating T₁ and T₂ times from each MRF dictionary with reference values. Repeatability was assessed using the coefficient of variation, with measurements taken over 30 separate sessions.

RESULTS

There was a statistically significant fit of the model for MRF-derived T₁ and T₂ and temperature (p < 0.05) for all the spheres with a T₁ > 500 ms. Both MRF methods showed a strong linear correlation with reference values for T₁ (R² = 0.996) and T₂ (R² = 0.982). MRF repeatability for T₁ values was ≤1.4% and for T₂ values was ≤3.4%.

CONCLUSION

MRF demonstrated relaxation times with a temperature dependence similar to that of conventional mapping methods. Temperature-corrected T₁ and T₂ values from both dictionaries showed adequate accuracy and excellent repeatability in this phantom study.

Current Applications and Future Development of Magnetic Resonance Fingerprinting in Diagnosis, Characterization, and Response Monitoring in Cancer

Magnetic resonance imaging (MRI) has enabled non-invasive cancer diagnosis, monitoring, and management in common clinical settings. However, inadequate quantitative analyses in MRI continue to limit its full potential and these often have an impact on clinicians' judgments. Magnetic resonance fingerprinting (MRF) has recently been introduced to acquire multiple quantitative parameters simultaneously in a reasonable timeframe. Initial retrospective studies have demonstrated the feasibility of using MRF for various cancer characterizations. Further trials with larger cohorts are still needed to explore the repeatability and reproducibility of the data acquired by MRF. At the moment, technical difficulties such as undesirable processing time or lack of motion robustness are limiting further implementations of MRF in clinical oncology. This review summarises the latest findings and technology developments for the use of MRF in cancer management and suggests possible future implications of MRF in characterizing tumour heterogeneity and response assessment.

Rapid 2D variable flip angle method for accurate and precise T1 measurements over a wide range of T1 values

PURPOSE

To perform dynamic T₁ mapping using a 2D variable flip angle (VFA) method, a correction for the slice profile effect is needed. In this work we investigated the impact of flip angle selection and excitation RF pulse profile on the performance of slice profile correction when applied to T₁ mapping over a range of T₁ values.

METHODS

A correction of the slice profile effect is proposed, based on Bloch simulation of steady-state signals. With this correction, Monte Carlo simulations were performed to assess the accuracy and precision of 2D VFA T₁ mapping in the presence of noise, for RF pulses with time-bandwidth products of 2, 3 and 10 and with flip angle pairs in the range [1°-90°]. To evaluate its performance over a wide range of T₁ , maximum errors were calculated for six T₁ values between 50 ms and 1250 ms. The method was demonstrated using in vitro and in vivo experiments.

RESULTS

Without corrections, 2D VFA severely underestimates T₁ . Slice profile errors were effectively reduced with the correction based on simulations, both in vitro and in vivo. The precision and accuracy of the method depend on the nominal T₁ values, the FA pair, and the RF pulse shape. FA pairs leading to <5% errors in T₁ can be identified for the common RF shapes, for T₁ values between 50 ms and 1250 ms.

CONCLUSIONS

2D VFA T₁ mapping with Bloch-simulation-based correction can deliver T₁ estimates that are accurate and precise to within 5% over a wide T₁ range.

Bayesian Inverse Regression for Vascular Magnetic Resonance Fingerprinting

Standard parameter estimation from vascular magnetic resonance fingerprinting (MRF) data is based on matching the MRF signals to their best counterparts in a grid of coupled simulated signals and parameters, referred to as a dictionary. To reach a good accuracy, the matching requires an informative dictionary whose cost, in terms of design, storage and exploration, is rapidly prohibitive for even moderate numbers of parameters. In this work, we propose an alternative dictionary-based statistical learning (DB-SL) approach made of three steps: 1) a quasi-random sampling strategy to produce efficiently an informative dictionary, 2) an inverse statistical regression model to learn from the dictionary a correspondence between fingerprints and parameters, and 3) the use of this mapping to provide both parameter estimates and their confidence indices. The proposed DB-SL approach is compared to both the standard dictionary-based matching (DBM) method and to a dictionary-based deep learning (DB-DL) method. Performance is illustrated first on synthetic signals including scalable and standard MRF signals with spatial undersampling noise. Then, vascular MRF signals are considered both through simulations and real data acquired in tumor bearing rats. Overall, the two learning methods yield more accurate parameter estimates than matching and to a range not limited to the dictionary boundaries. DB-SL in particular resists to higher noise levels and provides in addition confidence indices on the estimates at no additional cost. DB-SL appears as a promising method to reduce simulation needs and computational requirements, while modeling sources of uncertainty and providing both accurate and interpretable results.

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.

Quantification of T1, T2 relaxation times from Magnetic Resonance Fingerprinting radially undersampled data using analytical transformations

Quantitative magnetic resonance imaging (MRI) estimates magnetic parameters related to tissue, such as T1, T2 relaxation times and proton density. MR fingerprinting (MRF) is a new concept that uses pseudo-random, incoherent measurements to create a unique fingerprint for each tissue type to quantify magnet parameters. This paper aims to enhance MRF performance by investigating (i) the most suitable acquisition trajectory, and (ii) analytical transformations, suitable for radial acquisitions. Highly undersampled MRF brain (k, t)-space data have been simulated and non-linearly reconstructed to exploit the low-rank property of dynamic imaging. Based on our findings, the radial trajectory is the most suitable for MRF compared to Cartesian and spiral acquisitions. Perhaps this is due to the fact that its aliasing artifacts are more noise-like, and that unlike spiral trajectories, it can use analytical transformations that do not require re-gridding. One such analytical algorithm is the spline reconstruction technique (SRT) that is based on a novel numerical implementation of an analytic representation of the inverse Radon transform. Here, for the first time, this algorithm is applied to MR radial data. Reconstructions using SRT were compared to the ones using filtered back-projection. SRT provided images of higher contrast, lower bias, which resulted in more accurate T1, T2 values.

Quantitative imaging metrics derived from magnetic resonance fingerprinting using ISMRM/NIST MRI system phantom: An international multicenter repeatability and reproducibility study

PURPOSE

To compare the bias and inherent reliability of the quantitative (T₁ and T₂ ) imaging metrics generated from the magnetic resonance fingerprinting (MRF) technique using the ISMRM/NIST system phantom in an international multicenter setting.

METHOD

ISMRM/NIST MRI system phantom provides standard reference T₁ and T₂ relaxation values (vendor-provided) for each of the 14 vials in T₁ and T₂ arrays. MRF-SSFP scans repeated over 30 days on GE 1.5 and 3.0 T scanners at three collaborative centers. MRF estimated T_1, and T₂ values averaged over 30 days were compared with the phantom vendor-provided and spin-echo (SE) based convention gold standard (GS) method. Repeatability and reproducibility were characterized by the within-case coefficient of variation (wCV) of the MRF data acquired over 30 days, along with the biases.

RESULT

For the wide ranges of MRF estimated T₁ values, vials #1-8 (T₁ relaxation time between 2033 and 184 ms) exhibited a wCV less than 3% and 4%, respectively, on 3.0 and 1.5 T scanners. T₂ values in vials #1-8 (T₂ relaxation, 1044-45 ms) have shown wCV to be <7% on both 3.0 and 1.5 T scanners. A stronger linear correlation overall for T₁ (R²  = 0.9960 and 0.9963 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9951 and 0.9988 at center-1 and center-3 on 1.5 T scanner) compared to T₂ (R²  = 0.9971 and 0.9972 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9815 and 0.9754 at center-1 and center-3 on 1.5 T scanner). Bland-Altman (BA) analysis showed MRF based T₁ and T₂ values were within the limit of agreement (LOA) except for one data point. The mean difference or bias and 95% lower bound (LB) and upper bound (UB) LOA are reported in the format; mean bias: 95% LB LOA: 95% UB LOA. The biases for T₁ values were 21.34: -50.00: 92.69, 21.32: -47.29: 89.94 ms, and for T₂ values were -19.88: -42.37: 2.61, -19.06: -43.58: 5.45 ms on 3.0 T scanner at center-1 and center-2, respectively. Similarly, on 1.5 T scanner biases for T₁ values were 26.54: -53.41: 106.50, 9.997: -51.94: 71.94 ms, and for T₂ values were -23.84: -135.40: 87.76, -37.30: 134.30: 59.73 ms at center-1 and center-3, respectively. Additionally, the correlation between the SE based GS and MRF estimated T₁ and T₂ values (R²  = 0.9969 and 0.9977) showed a similar trend as we observed between vendor-provided and MRF estimated T₁ and T₂ values (R²  = 0.9963 and 0.9972). In addition to correlation, BA analysis showed that all the vials are within the LOA between the GS and vendor-provided for the T₁ values and except one vial for T₂ . All the vials are within the LOA between GS and MRF except one vial in T₁ and T₂ array. The wCV for reproducibility was <3% for both T₁ and T₂ values in vials #1-8, for all the 14 vials, wCV calculated for reproducibility was <4% for T₁ values and <5% for T₂ .

CONCLUSION

This study shows that MRF is highly repeatable (wCV <4% for T₁ and <7% for T₂ ) and reproducible (wCV < 3% for both T₁ and T₂ ) in certain vials (vials #1-8).

Simultaneous T1 , T2 , and T1ρ relaxation mapping of the lower leg muscle with MR fingerprinting

PURPOSE

To develop a novel MR-fingerprinting (MRF) pulse sequence that is insensitive to B 1 + and B₀ imperfections for simultaneous T₁ , T₂ , and T_1ρ relaxation mapping.

METHODS

We implemented a totally balanced spin-lock (TB-SL) module to encode T_1ρ relaxation into an existing MRF framework that encoded T₁ and T₂ . The spin-lock module used two 180° pulses with compensatory phases to reduce T_1ρ sensitivity to B₁ and B₀ inhomogeneities. We compared T_1ρ measured using TB-SL MRF in Bloch simulations, model agar phantoms, and in vivo experiments to those with a self-compensated spin-lock preparation module (SC-SL). The TB-SL MRF repeatability was evaluated in maps acquired in the lower leg skeletal muscle of 12 diabetic peripheral neuropathy patients, scanned two times each during visits separated by about 30 days.

RESULTS

The phantom relaxation times measured with TB-SL and SC-SL MRF were in good agreement with reference values in regions with low B₁ inhomogeneities. Compared with SC-SL, TB-SL MRF showed in experiments greater robustness against severe B₁ inhomogeneities and in Bloch simulations greater robustness against B₁ and B₀ . We measured with TB-SL MRF an average T₁ = 950.1 ± 28.7 ms, T₂ = 26.0 ± 1.2 ms, and T_1ρ = 31.7 ± 3.2 ms in skeletal muscle across patients. Bland-Altman analysis demonstrated low bias between TB-SL and SC-SL MRF and between TB-SL MRF maps acquired in two visits. The coefficient of variation was less than 3% for all measurements.

CONCLUSION

The proposed TB-SL MRF sequence is fast and insensitive to B 1 + and B₀ imperfections. It can simultaneously map T₁ , T₂ , T_1ρ , and B 1 + in a single scan and can potentially be used to study muscle composition.

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.

Reproducibility of magnetic resonance fingerprinting-based T1 mapping of the healthy prostate at 1.5 and 3.0 T: A proof-of-concept study

Facilitating clinical translation of quantitative imaging techniques has been suggested as means of improving interobserver agreement and diagnostic accuracy of multiparametric magnetic resonance imaging (mpMRI) of the prostate. One such technique, magnetic resonance fingerprinting (MRF), has significant competitive advantages over conventional mapping techniques in terms of its multi-site reproducibility, short scanning time and inherent robustness to motion. It has also been shown to improve the detection of clinically significant prostate cancer when added to standard mpMRI sequences, however, the existing studies have all been conducted on 3.0 T MRI systems, limiting the technique's use on 1.5 T MRI scanners that are still more widely used for prostate imaging across the globe. The aim of this proof-of-concept study was, therefore, to evaluate the cross-system reproducibility of prostate MRF T1 in healthy volunteers (HVs) using 1.5 and 3.0 T MRI systems. The initial validation of MRF T1 against gold standard inversion recovery fast spin echo (IR-FSE) T1 in the ISMRM/NIST MRI system revealed a strong linear correlation between phantom-derived MRF and IR-FSE T1 values was observed at both field strengths (R2 = 0.998 at 1.5T and R2 = 0.993 at 3T; p = < 0.0001 for both). In young HVs, inter-scanner CVs demonstrated marginal differences across all tissues with the highest difference of 3% observed in fat (2% at 1.5T vs 5% at 3T). At both field strengths, MRF T1 could confidently differentiate prostate peripheral zone from transition zone, which highlights the high quantitative potential of the technique given the known difficulty of tissue differentiation in this age group. The high cross-system reproducibility of MRF T1 relaxometry of the healthy prostate observed in this preliminary study, therefore, supports the technique's prospective clinical validation as part of larger trials employing 1.5 T MRI systems, which are still widely used clinically for routine mpMRI of the prostate.

Sensitive detection of extremely small iron oxide nanoparticles in living mice using MP2RAGE with advanced image co-registration

Magnetic resonance imaging (MRI) is a widely used non-invasive methodology for both preclinical and clinical studies. However, MRI lacks molecular specificity. Molecular contrast agents for MRI would be highly beneficial for detecting specific pathological lesions and quantitatively evaluating therapeutic efficacy in vivo. In this study, an optimized Magnetization Prepared—RApid Gradient Echo (MP-RAGE) with 2 inversion times called MP2RAGE combined with advanced image co-registration is presented as an effective non-invasive methodology to quantitatively detect T1 MR contrast agents. The optimized MP2RAGE produced high quality in vivo mouse brain T1 (or R1 = 1/T1) map with high spatial resolution, 160 × 160 × 160 µm³ voxel at 9.4 T. Test–retest signal to noise was > 20 for most voxels. Extremely small iron oxide nanoparticles (ESIONPs) having 3 nm core size and 11 nm hydrodynamic radius after polyethylene glycol (PEG) coating were intracranially injected into mouse brain and detected as a proof-of-concept. Two independent MP2RAGE MR scans were performed pre- and post-injection of ESIONPs followed by advanced image co-registration. The comparison of two T1 (or R1) maps after image co-registration provided precise and quantitative assessment of the effects of the injected ESIONPs at each voxel. The proposed MR protocol has potential for future use in the detection of T1 molecular contrast agents.

The effect of gadolinium-based contrast agent administration on magnetic resonance fingerprinting-based T1 relaxometry in patients with prostate cancer

Magnetic resonance fingerprinting (MRF) is a rapidly developing fast quantitative mapping technique able to produce multiple property maps with reduced sensitivity to motion. MRF has shown promise in improving the diagnosis of clinically significant prostate cancer but requires further validation as part of a prostate multiparametric (mp) MRI protocol. mpMRI protocol mandates the inclusion of dynamic contrast enhanced (DCE) imaging, known for its significant T1 shortening effect. MRF could be used to measure both pre- and post-contrast T1 values, but its utility must be assessed. In this proof-of-concept study, we sought to evaluate the variation in MRF T1 measurements post gadolinium-based contrast agent (GBCA) injection and the utility of such T1 measurements to differentiate peripheral and transition zone tumours from normal prostatic tissue. We found that the T1 variation in all tissues increased considerably post-GBCA following the expected significant T1 shortening effect, compromising the ability of MRF T1 to identify transition zone lesions. We, therefore, recommend performing MRF T1 prior to DCE imaging to maintain its benefit for improving detection of both peripheral and transition zone lesions while reducing additional scanning time. Demonstrating the effect of GBCA on MRF T1 relaxometry in patients also paves the way for future clinical studies investigating the added value of post-GBCA MRF in PCa, including its dynamic analysis as in DCE-MRF.

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 assessment of the cerebral alterations associated with obesity development

Obesity is a current threat to health care systems, affecting approximately 13% of the world's adult population, and over 18% children and adolescents. The rise of obesity is fuelled by inadequate life style habits, as consumption of diets rich in fats and sugars which promote, additionally, the development of associated comorbidities. Obesity results from a neuroendocrine imbalance in the cerebral mechanisms controlling food intake and energy expenditure, including the hypothalamus and the reward and motivational centres. Specifically, high-fat diets are known to trigger an early inflammatory response in the hypothalamus that precedes weight gain, is time-dependent, and eventually extends to the remaining appetite regulating regions in the brain. Multiple magnetic resonance imaging (MRI) and spectroscopy (MRS) methods are currently available to characterize different features of cerebral obesity, including diffusion weighted, T2 and volumetric imaging and 1H and 13C spectroscopic evaluations. In particular, consistent evidences have revealed increased water diffusivity and T2 values, decreased grey matter volumes, and altered metabolic profiles and fluxes, in the brain of animal models and in obese humans. This review provides an integrative interpretation of the physio-pathological processes associated with obesity development in the brain, and the MRI and MRS methods implemented to characterize them.

Repeatability and reproducibility of human brain morphometry using three-dimensional magnetic resonance fingerprinting

Three-dimensional (3D) Magnetic resonance fingerprinting (MRF) permits whole-brain volumetric quantification of T1 and T2 relaxation values, potentially replacing conventional T1-weighted structural imaging for common brain imaging analysis. The aim of this study was to evaluate the repeatability and reproducibility of 3D MRF in evaluating brain cortical thickness and subcortical volumetric analysis in healthy volunteers using conventional 3D T1-weighted images as a reference standard. Scan-rescan tests of both 3D MRF and conventional 3D fast spoiled gradient recalled echo (FSPGR) were performed. For each sequence, the regional cortical thickness and volume of the subcortical structures were measured using standard automatic brain segmentation software. Repeatability and reproducibility were assessed using the within-subject coefficient of variation (wCV), intraclass correlation coefficient (ICC), and mean percent difference and ICC, respectively. The wCV and ICC of cortical thickness were similar across all regions with both 3D MRF and FSPGR. The percent relative difference in cortical thickness between 3D MRF and FSPGR across all regions was 8.0 ± 3.2%. The wCV and ICC of the volume of subcortical structures across all structures were similar between 3D MRF and FSPGR. The percent relative difference in the volume of subcortical structures between 3D MRF and FSPGR across all structures was 7.1 ± 3.6%. 3D MRF measurements of human brain cortical thickness and subcortical volumes are highly repeatable, and consistent with measurements taken on conventional 3D T1-weighted images. A slight, consistent bias was evident between the two, and thus careful attention is required when combining data from MRF and conventional acquisitions.

Contemporary and emerging magnetic resonance imaging methods for evaluation of moyamoya disease

Numerous recent technological advances offer the potential to substantially enhance the MRI evaluation of moyamoya disease (MMD). These include high-resolution volumetric imaging, high-resolution vessel wall characterization, improved cerebral angiographic and perfusion techniques, high-field imaging, fast scanning methods, and artificial intelligence. This review discusses the current state-of-the-art MRI applications in these realms, emphasizing key imaging findings, clinical utility, and areas that will benefit from further investigation. Although these techniques may apply to imaging of a wide array of neurovascular or other neurological conditions, consideration of their application to MMD is useful given the comprehensive multidimensional MRI assessment used to evaluate MMD. These MRI techniques span from basic cross-sectional to advanced functional sequences, both qualitative and quantitative.The aim of this review was to provide a comprehensive summary and analysis of current key relevant literature of advanced MRI techniques for the evaluation of MMD with image-rich case examples. These imaging methods can aid clinical characterization, help direct treatment, assist in the evaluation of treatment response, and potentially improve the understanding of the pathophysiology of MMD.

Acceleration of High-Resolution 3D MR Fingerprinting via a Graph Convolutional Network

Magnetic resonance fingerprinting (MRF) is a novel imaging framework for fast and simultaneous quantification of multiple tissue properties. Recently, 3D MRF methods have been developed, but the acquisition speed needs to be improved before they can be adopted for clinical use. The purpose of this study is to develop a novel deep learning approach to accelerate 3D MRF acquisition along the slice-encoding direction in k-space. We introduce a graph-based convolutional neural network that caters to non-Cartesian spiral trajectories commonly used for MRF acquisition. We improve tissue quantification accuracy compared with the state of the art. Our method enables fast 3D MRF with high spatial resolution, allowing whole-brain coverage within 5min, making MRF more feasible in clinical settings.

Variability and Standardization of Quantitative Imaging: Monoparametric to Multiparametric Quantification, Radiomics, and Artificial Intelligence

Radiological images have been assessed qualitatively in most clinical settings by the expert eyes of radiologists and other clinicians. On the other hand, quantification of radiological images has the potential to detect early disease that may be difficult to detect with human eyes, complement or replace biopsy, and provide clear differentiation of disease stage. Further, objective assessment by quantification is a prerequisite of personalized/precision medicine. This review article aims to summarize and discuss how the variability of quantitative values derived from radiological images are induced by a number of factors and how these variabilities are mitigated and standardization of the quantitative values are achieved. We discuss the variabilities of specific biomarkers derived from magnetic resonance imaging and computed tomography, and focus on diffusion-weighted imaging, relaxometry, lung density evaluation, and computer-aided computed tomography volumetry. We also review the sources of variability and current efforts of standardization of the rapidly evolving techniques, which include radiomics and artificial intelligence.

A Perspective on MR Fingerprinting

This article reviews the basic concept of MR fingerprinting (MRF) with the goal of highlighting MRF's key contributions, putting them in the context of other quantitative MRI literature, and refining MRF's terminology. The article discusses the robustness and flexibility of MRF's signature dictionary-matching reconstruction along with more advanced MRF reconstructions. A key feature of MRF is the lack of assumptions about the signal evolution, which gives scientists the flexibility to tailor sequences for their needs. The article argues that the concept of unique fingerprints does not capture the requirements for successful parameter mapping and that an analysis of the signal's derivatives with respect to biophysical parameters, such as relaxation times, is more informative, as it allows one to evaluate the efficiency of a pulse sequence. The article points at the source of MRF's efficiency, namely, flip angle variations at the time scale of the relaxation times, and reveals that MRF's advantages are strongest at long scan times, as required for 3D imaging. Further, it outlines how MRF's flexibility can be used to design mutually tailored pulse sequences and biophysical models with the goal of improving the reproducibility of parameter mapping biological tissue. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 1.

Magnetic resonance fingerprinting review part 2: Technique and directions

Magnetic resonance fingerprinting (MRF) is a general framework to quantify multiple MR-sensitive tissue properties with a single acquisition. There have been numerous advances in MRF in the years since its inception. In this work we highlight some of the recent technical developments in MRF, focusing on sequence optimization, modifications for reconstruction and pattern matching, new methods for partial volume analysis, and applications of machine and deep learning. Level of Evidence: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:993-1007.