Spiral MR fingerprinting at 7T with simultaneous B1 estimation

QRAGE-Simultaneous multiparametric quantitative MRI of water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength

PURPOSE

To introduce quantitative rapid gradient-echo (QRAGE), a novel approach for the simultaneous mapping of multiple quantitative MRI parameters, including water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength.

METHODS

QRAGE leverages a newly developed multi-echo MPnRAGE sequence, facilitating the acquisition of 171 distinct contrast images across a range of inversion and TE points. To maintain a short acquisition time, we introduce MIRAGE2, a novel model-based reconstruction method that exploits prior knowledge of temporal signal evolution, represented as damped complex exponentials. MIRAGE2 minimizes local Block-Hankel and Casorati matrices. Parameter maps are derived from the reconstructed contrast images through postprocessing steps. We validate QRAGE through extensive simulations, phantom studies, and in vivo experiments, demonstrating its capability for high-precision imaging.

RESULTS

In vivo brain measurements show the promising performance of QRAGE, with test-retest SDs and deviations from reference methods of < 0.8% for water content, < 17 ms for T1, and < 0.7 ms for T2*. QRAGE achieves whole-brain coverage at a 1-mm isotropic resolution in just 7 min and 15 s, comparable to the acquisition time of an MP2RAGE scan. In addition, QRAGE generates a contrast image akin to the UNI image produced by MP2RAGE.

CONCLUSION

QRAGE is a new, successful approach for simultaneously mapping multiple MR parameters at ultrahigh field.

High SNR full brain relaxometry at 7T by accelerated MR-STAT

PURPOSE

To demonstrate the feasibility and robustness of the Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) framework for fast, high SNR relaxometry at 7T.

METHODS

To deploy MR-STAT on 7T-systems, we designed optimized flip-angles using the BLAKJac-framework that incorporates the SAR-constraints. Transmit RF-inhomogeneities were mitigated by including a measuredB 1 + $$ {B}_1^{+} $$ -map in the reconstruction. Experiments were performed on a gel-phantom and on five volunteers to explore the robustness of the sequence and its sensitivity toB 1 + $$ {B}_1^{+} $$ inhomogeneities. The SNR-gain at 7T was explored by comparing phantom and in vivo results to MR-STAT at 3T in terms of SNR-efficiency.

RESULTS

The higher SNR at 7T enabled two-fold acceleration with respect to current 2D MR-STAT protocols at lower field strengths. The resulting scan had whole-brain coverage, with 1 x 1 x 3 mm3 resolution (1.5 mm slice-gap) and was acquired within 3 min including theB 1 + $$ {B}_1^{+} $$ -mapping. AfterB 1 + $$ {B}_1^{+} $$ -correction, the estimated T1 and T2 in a phantom showed a mean relative error of, respectively, 1.7% and 4.4%. In vivo, the estimated T1 and T2 in gray and white matter corresponded to the range of values reported in literature with a variation over the subjects of 1.0%-2.1% (WM-GM) for T1 and 4.3%-5.3% (WM-GM) for T2. We measured a higher SNR-efficiency at 7T (R = 2) than at 3T for both T1 and T2 with, respectively, a 4.1 and 2.3 times increase in SNR-efficiency.

CONCLUSION

We presented an accelerated version of MR-STAT tailored to high field (7T) MRI using a low-SAR flip-angle train and showed high quality parameter maps with an increased SNR-efficiency compared to MR-STAT at 3T.

Quantitative brain mapping using magnetic resonance fingerprinting on a 50-mT portable MRI scanner

Ultralow-field magnetic resonance imaging (ULF-MRI) has broad application prospects because of its portable hardware system and low cost. However, the low B0 magnitude of ULF-MRI results in a reduced signal-to-noise ratio in qualitative images compared with that of commercial high-field MRI, which can affect the visibility and delineation of tissues and lesions. In this work, a magnetic resonance fingerprinting (MRF) approach is applied to a homemade 50-mT ULF-MRI scanner to achieve efficient quantitative brain imaging, which is an original and promising disease-diagnosis approach for portable MRI systems. An inversion recovery fast imaging with steady-state precession-based sequence is utilized for MRF through Cartesian acquisition. A microdictionary analysis method is proposed to select the optimal repetition time and flip angle variation schedule and ensure the best possible tissue discriminative ability of MRF. The T1 and T2 relaxation properties and the B1 + distribution are considered for estimation, and the results are compared with those of gold standard (GS) quantitative imaging or qualitative imaging methods. The phantom experiment indicates that the quantitative values obtained by schedule-optimized MRF show good agreement, and the bias from the GS results is acceptable. The in vivo experiment shows that the relaxation times of white and gray matter estimated by MRF are slightly lower than the reference data, and the relaxation times of lipid are within the range of the reference data. Compared with qualitative MRI under ULF, MRF can intuitively reflect various items of brain tissue information in a single scan, so it is a valuable addition to point-of-care imaging approaches.

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.

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

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

Fast high-resolution electric properties tomography using three-dimensional quantitative transient-state imaging-based water fraction estimation

In this study, we aimed to develop a fast and robust high-resolution technique for clinically feasible electrical properties tomography based on water content maps (wEPT) using Quantitative Transient-state Imaging (QTI), a multiparametric transient state-based method that is similar to MR fingerprinting. Compared with the original wEPT implementation based on standard spin-echo acquisition, QTI provides robust electrical properties quantification towards B1 + inhomogeneities and full quantitative relaxometry data. To validate the proposed approach, 3D QTI data of 12 healthy volunteers were acquired on a 1.5 T scanner. QTI-provided T1 maps were used to compute water content maps of the tissues using an empirical relationship based on literature ex-vivo measurements. Assuming that electrical properties are modulated mainly by tissue water content, the water content maps were used to derive electrical conductivity and relative permittivity maps. The proposed technique was compared with a conventional phase-only Helmholtz EPT (HH-EPT) acquisition both within whole white matter, gray matter, and cerebrospinal fluid masks, and within different white and gray matter subregions. In addition, QTI-based wEPT was retrospectively applied to four multiple sclerosis adolescent and adult patients, compared with conventional contrast-weighted imaging in terms of lesion delineation, and quantitatively assessed by measuring the variation of electrical properties in lesions. Results obtained with the proposed approach agreed well with theoretical predictions and previous in vivo findings in both white and gray matter. The reconstructed maps showed greater anatomical detail and lower variability compared with standard phase-only HH-EPT. The technique can potentially improve delineation of pathology when compared with conventional contrast-weighted imaging and was able to detect significant variations in lesions with respect to normal-appearing tissues. In conclusion, QTI can reliably measure conductivity and relative permittivity of brain tissues within a short scan time, opening the way to the study of electric properties in clinical settings.

Cardiac MR Fingerprinting: Overview, Technical Developments, and Applications

Cardiovascular magnetic resonance (CMR) is an established imaging modality with proven utility in assessing cardiovascular diseases. The ability of CMR to characterize myocardial tissue using T1 - and T2 -weighted imaging, parametric mapping, and late gadolinium enhancement has allowed for the non-invasive identification of specific pathologies not previously possible with modalities like echocardiography. However, CMR examinations are lengthy and technically complex, requiring multiple pulse sequences and different anatomical planes to comprehensively assess myocardial structure, function, and tissue composition. To increase the overall impact of this modality, there is a need to simplify and shorten CMR exams to improve access and efficiency, while also providing reproducible quantitative measurements. Multiparametric MRI techniques that measure multiple tissue properties offer one potential solution to this problem. This review provides an in-depth look at one such multiparametric approach, cardiac magnetic resonance fingerprinting (MRF). The article is structured as follows. First, a brief review of single-parametric and (non-Fingerprinting) multiparametric CMR mapping techniques is presented. Second, a general overview of cardiac MRF is provided covering pulse sequence implementation, dictionary generation, fast k-space sampling methods, and pattern recognition. Third, recent technical advances in cardiac MRF are covered spanning a variety of topics, including simultaneous multislice and 3D sampling, motion correction algorithms, cine MRF, synthetic multicontrast imaging, extensions to measure additional clinically important tissue properties (proton density fat fraction, T2 *, and T1ρ ), and deep learning methods for image reconstruction and parameter estimation. The last section will discuss potential clinical applications, concluding with a perspective on how multiparametric techniques like MRF may enable streamlined CMR protocols. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

Fast in vivo assay of creatine kinase activity in the human brain by 31 P magnetic resonance fingerprinting

A new and efficient magnetisation transfer 31 P magnetic resonance fingerprinting (MT-31 P-MRF) approach is introduced to measure the creatine kinase metabolic ratek CK between phosphocreatine (PCr) and adenosine triphosphate (ATP) in human brain. The MRF framework is extended to overcome challenges in conventional 31 P measurement methods in the human brain, enabling reduced acquisition time and specific absorption rate (SAR). To address the challenge of creating and matching large multiparametric dictionaries in an MRF scheme, a nested iteration interpolation method (NIIM) is introduced. As the number of parameters to estimate increases, the size of the dictionary grows exponentially. NIIM can reduce the computational load by breaking dictionary matching into subsolutions of linear computational order. MT-31 P-MRF combined with NIIM providesT 1 PCr ,T 1 ATP andk CK estimates in good agreement with those obtained by the exchange kinetics by band inversion transfer (EBIT) method and literature values. Furthermore, the test-retest reproducibility results showed that MT-31 P-MRF achieves a similar or better coefficient of variation (<12%) forT 1 ATP andk CK measurements in 4 min 15 s, than EBIT with 17 min 4 s scan time, enabling a fourfold reduction in scan time. We conclude that MT-31 P-MRF in combination with NIIM is a fast, accurate, and reproducible approach for in vivok CK assays in the human brain, which enables the potential to investigate energy metabolism in a clinical setting.

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.

A Self-Supervised Deep Learning Reconstruction for Shortening the Breathhold and Acquisition Window in Cardiac Magnetic Resonance Fingerprinting

The aim of this study is to shorten the breathhold and diastolic acquisition window in cardiac magnetic resonance fingerprinting (MRF) for simultaneous T1, T2, and proton spin density (M0) mapping to improve scan efficiency and reduce motion artifacts. To this end, a novel reconstruction was developed that combines low-rank subspace modeling with a deep image prior, termed DIP-MRF. A system of neural networks is used to generate spatial basis images and quantitative tissue property maps, with training performed using only the undersampled k-space measurements from the current scan. This approach avoids difficulties with obtaining in vivo MRF training data, as training is performed de novo for each acquisition. Calculation of the forward model during training is accelerated by using GRAPPA operator gridding to shift spiral k-space data to Cartesian grid points, and by using a neural network to rapidly generate fingerprints in place of a Bloch equation simulation. DIP-MRF was evaluated in simulations and at 1.5 T in a standardized phantom, 18 healthy subjects, and 10 patients with suspected cardiomyopathy. In addition to conventional mapping, two cardiac MRF sequences were acquired, one with a 15-heartbeat(HB) breathhold and 254 ms acquisition window, and one with a 5HB breathhold and 150 ms acquisition window. In simulations, DIP-MRF yielded decreased nRMSE compared to dictionary matching and a sparse and locally low rank (SLLR-MRF) reconstruction. Strong correlation (R2 &gt; 0.999) with T1 and T2 reference values was observed in the phantom using the 5HB/150 ms scan with DIP-MRF. DIP-MRF provided better suppression of noise and aliasing artifacts in vivo, especially for the 5HB/150 ms scan, and lower intersubject and intrasubject variability compared to dictionary matching and SLLR-MRF. Furthermore, it yielded a better agreement between myocardial T1 and T2 from 15HB/254 ms and 5HB/150 ms MRF scans, with a bias of −9 ms for T1 and 2 ms for T2. In summary, this study introduces an extension of the deep image prior framework for cardiac MRF tissue property mapping, which does not require pre-training with in vivo scans, and has the potential to reduce motion artifacts by enabling a shortened breathhold and acquisition window.

Optimized multi-axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole-brain high-isotropic-resolution quantitative imaging

PURPOSE

To improve image quality and accelerate the acquisition of 3D MR fingerprinting (MRF).

METHODS

Building on the multi-axis spiral-projection MRF technique, a subspace reconstruction with locally low-rank constraint and a modified spiral-projection spatiotemporal encoding scheme called tiny golden-angle shuffling were implemented for rapid whole-brain high-resolution quantitative mapping. Reconstruction parameters such as the locally low-rank regularization parameter and the subspace rank were tuned using retrospective in vivo data and simulated examinations. B₀ inhomogeneity correction using multifrequency interpolation was incorporated into the subspace reconstruction to further improve the image quality by mitigating blurring caused by off-resonance effect.

RESULTS

The proposed MRF acquisition and reconstruction framework yields high-quality 1-mm isotropic whole-brain quantitative maps in 2 min at better quality compared with 6-min acquisitions of prior approaches. The proposed method was validated to not induce bias in T₁ and T₂ mapping. High-quality whole-brain MRF data were also obtained at 0.66-mm isotropic resolution in 4 min using the proposed technique, where the increased resolution was shown to improve visualization of subtle brain structures.

CONCLUSIONS

The proposed tiny golden-angle shuffling, MRF with optimized spiral-projection trajectory and subspace reconstruction enables high-resolution quantitative mapping in ultrafast acquisition time.

Cartesian dictionary-based native T1 and T2 mapping of the myocardium

PURPOSE

To implement and evaluate a new dictionary-based technique for native myocardial T₁ and T₂ mapping using Cartesian sampling.

METHODS

The proposed technique (Multimapping) consisted of single-shot Cartesian image acquisitions in 10 consecutive cardiac cycles, with inversion pulses in cycle 1 and 5, and T₂ preparation (TE: 30 ms, 50 ms, and 70 ms) in cycles 8-10. Multimapping was simulated for different T₁ and T₂ , where entries corresponding to the k-space centers were matched to acquired data. Experiments were performed in a phantom, 16 healthy subjects, and 3 patients with cardiovascular disease.

RESULTS

Multimapping phantom measurements showed good agreement with reference values for both T₁ and T₂ , with no discernable heart-rate dependency for T₁ and T₂ within the range of myocardium. In vivo mean T₁ in healthy subjects was significantly higher using Multimapping (T₁ = 1114 ± 14 ms) compared to the reference (T₁ = 991 ± 26 ms) (p < 0.01). Mean Multimapping T₂ (47.1 ± 1.3 ms) and T₂ spatial variability (5.8 ± 1.0 ms) was significantly lower compared to the reference (T₂ = 54.7 ± 2.2 ms, p < 0.001; spatial variability = 8.4 ± 2.0 ms, p < 0.01). Increased T₁ and T₂ was detected in all patients using Multimapping.

CONCLUSIONS

Multimapping allows for simultaneous native myocardial T₁ and T₂ mapping with a conventional Cartesian trajectory, demonstrating promising in vivo image quality and parameter quantification results.

Quantitative MR relaxation using MR fingerprinting with fractional-order signal evolution

The fractional-order Bloch equations have been shown to describe a wider range of experimental situations involving heterogeneous, porous, or composite materials. This paper introduces a novel dictionary of quantitative MR fingerprinting generated by signal evolution model with fractional-order Bloch equations to describe magnetic resonance (MR) relaxation. Here, the fractional-order relaxation models are implemented into Bloch equations through phase transitions using EPG simulation. In the phantom experiments, the fractional-order analysis showed smaller root mean squared error (T1: RMSE = 5.21%, T2: RMSE=3.75%) using the proposed method compared to using conventional method. Among the in vivo experiments of human brains, the estimated T1 and T2 values (mean ± SD) were 843 ± 46.3 ms and 70 ± 4.7 ms in white matter, 1323 ± 28.5 ms and 95 ± 3.8 ms in gray matter. So the proposed method can provide well extensions of current MR fingerprinting and has shown potential to apply into the phantom experiments and the in vivo applications to approach the standard methods for quantitative imaging.

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

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

Rapid B1-Insensitive MR Fingerprinting for Quantitative Kidney Imaging

Background MR fingerprinting (MRF) provides rapid and simultaneous quantification of multiple tissue parameters in a single scan. Purpose To evaluate a rapid kidney MRF technique at 3.0 T in phantoms, healthy volunteers, and patients. Materials and Methods A 15-second kidney MRF acquisition was designed with 12 acquisition segments, a range of low flip angles (5°-12°), multiple magnetization preparation schema (T1, T2, and fat suppression), and an undersampled spiral trajectory. This technique was first validated in vitro using standardized T1 and T2 phantoms. Kidney T1 and T2 maps were then obtained for 10 healthy adult volunteers (mean age ± standard deviation, 35 years ± 13; six men) and three pediatric patients with autosomal recessive polycystic kidney disease (ARPKD) (mean age, 10 years ± 3; two boys) between August 2019 and October 2020 to evaluate the method in vivo. Results Results in nine phantoms showed good agreement with spin-echo-based T1 and T2 values (R² > 0.99). In vivo MRF kidney T1 and T2 assessments in healthy adult volunteers (cortex: T1, 1362 msec ± 5; T2, 64 msec ± 5; medulla: T1, 1827 msec ± 94; T2, 69 msec ± 3) were consistent with values in the literature but with improved precision in comparison with prior MRF implementations. In vivo MRF-based kidney T1 and T2 values with and without B₁ correction were in good agreement (R² > 0.96, P < .001), demonstrating limited sensitivity to B₁ field inhomogeneities. Additional MRF reconstructions using the first nine segments of the MRF profiles (11-second acquisition time) were in good agreement with the reconstructions using 12 segments (15-second acquisition time) (R² > 0.87, P < .001). Repeat kidney MRF scans for the three patients with ARPKD on successive days also demonstrated good reproducibility (T1 and T2: <3% difference). Conclusion A kidney MR fingerprinting method provided in vivo kidney T1 and T2 maps at 3.0 T in a single breath hold with improved precision and no need for B₁ correction. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Laustsen in this issue.

Three dimensional MRF obtains highly repeatable and reproducible multi-parametric estimations in the healthy human brain at 1.5T and 3T

Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories. Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated a high repeatability (CVs 0.7-1.3% for T1, 2.0-7.8% for T2, 1.4-2.5% for normalized PD) and reproducibility (CVs of 2.0-5.8% for T1, 7.4-10.2% for T2, 5.2-9.2% for normalized PD) in gray and white matter. Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.

Confounding factors in breast magnetic resonance fingerprinting: B 1 + , slice profile, and diffusion effects

PURPOSE

Magnetic resonance fingerprinting (MRF) offers rapid quantitative imaging but may be subject to confounding effects (CE) if these are not included in the model-based reconstruction. This study characterizes the influence of in-plane B 1 + , slice profile and diffusion effects on T₁ and T₂ estimation in the female breast at 1.5T.

METHODS

Simulations were used to predict the influence of each CE on the accuracy of MRF and to investigate the influence of electronic noise and spiral aliasing artefacts. The experimentally observed bias in regions of fibroglandular tissue (FGT) and fatty tissue (FT) was analyzed for undersampled spiral breast MRF data of 6 healthy volunteers by performing MRF reconstruction with and without a CE.

RESULTS

Theoretic analysis predicts T₁ under-/T₂ overestimation if the nominal flip angles are underestimated and inversely, T₁ under-/T₂ overestimation if omitting slice profile correction, and T₁ under-/T₂ underestimation if omitting diffusion in the signal model. Averaged over repeated signal simulations, including spiral aliasing artefacts affected precision more than accuracy. Strong in-plane B 1 + effects occurred in vivo, causing T₂ left-right inhomogeneity between both breasts. Their correction decreased the T₂ difference from 29 to 5 ms in FGT and from 29 to 9 ms in FT. Slice profile correction affected FGT T₂ most strongly, resulting in -22% smaller values. For the employed spoiler gradient strengths, diffusion did not affect the parameter maps, corresponding well with theoretic predictions.

CONCLUSION

Understanding CEs and their relative significance for an MRF sequence is important when defining an MRF signal model for accurate parameter mapping.

Deep learning reconstruction for cardiac magnetic resonance fingerprinting T1 and T2 mapping

PURPOSE

To develop a deep learning method for rapidly reconstructing T₁ and T₂ maps from undersampled electrocardiogram (ECG) triggered cardiac magnetic resonance fingerprinting (cMRF) images.

METHODS

A neural network was developed that outputs T₁ and T₂ values when given a measured cMRF signal time course and cardiac RR interval times recorded by an ECG. Over 8 million cMRF signals, corresponding to 4000 random cardiac rhythms, were simulated for training. The training signals were corrupted by simulated k-space undersampling artifacts and random phase shifts to promote robust learning. The deep learning reconstruction was evaluated in Monte Carlo simulations for a variety of cardiac rhythms and compared with dictionary-based pattern matching in 58 healthy subjects at 1.5T.

RESULTS

In simulations, the normalized root-mean-square error (nRMSE) for T₁ was below 1% in myocardium, blood, and liver for all tested heart rates. For T₂ , the nRMSE was below 4% for myocardium and liver and below 6% for blood for all heart rates. The difference in the mean myocardial T₁ or T₂ observed in vivo between dictionary matching and deep learning was 3.6 ms for T₁ and -0.2 ms for T₂ . Whereas dictionary generation and pattern matching required more than 4 min per slice, the deep learning reconstruction only required 336 ms.

CONCLUSION

A neural network is introduced for reconstructing cMRF T₁ and T₂ maps directly from undersampled spiral images in under 400 ms and is robust to arbitrary cardiac rhythms, which paves the way for rapid online display of cMRF maps.

Musculoskeletal MRI at 7 T: do we need more or is it more than enough?

Ultra-high field magnetic resonance imaging (UHF-MRI) provides important diagnostic improvements in musculoskeletal imaging. The higher signal-to-noise ratio leads to higher spatial and temporal resolution which results in improved anatomic detail and higher diagnostic confidence. Several methods, such as T2, T2*, T1rho mapping, delayed gadolinium-enhanced, diffusion, chemical exchange saturation transfer, and magnetisation transfer techniques, permit a better tissue characterisation. Furthermore, UHF-MRI enables in vivo measurements by low-γ nuclei (²³Na, ³¹P, ¹³C, and ³⁹K) and the evaluation of different tissue metabolic pathways. European Union and Food and Drug Administration approvals for clinical imaging at UHF have been the first step towards a more routinely use of this technology, but some drawbacks are still present limiting its widespread clinical application. This review aims to provide a clinically oriented overview about the application of UHF-MRI in the different anatomical districts and tissues of musculoskeletal system and its pros and cons. Further studies are needed to consolidate the added value of the use of UHF-MRI in the routine clinical practice and promising efforts in technology development are already in progress.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

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.

Targeted Biopsy Validation of Peripheral Zone Prostate Cancer Characterization With Magnetic Resonance Fingerprinting and Diffusion Mapping

OBJECTIVE

This study aims for targeted biopsy validation of magnetic resonance fingerprinting (MRF) and diffusion mapping for characterizing peripheral zone (PZ) prostate cancer and noncancers.

MATERIALS AND METHODS

One hundred four PZ lesions in 85 patients who underwent magnetic resonance imaging were retrospectively analyzed with apparent diffusion coefficient (ADC) mapping, MRF, and targeted biopsy (cognitive or in-gantry). A radiologist blinded to pathology drew regions of interest on targeted lesions and visually normal peripheral zone on MRF and ADC maps. Mean T1, T2, and ADC were analyzed using linear mixed models. Generalized estimating equations logistic regression analyses were used to evaluate T1 and T2 relaxometry combined with ADC in differentiating pathologic groups.

RESULTS

Targeted biopsy revealed 63 cancers (low-grade cancer/Gleason score 6 = 10, clinically significant cancer/Gleason score ≥7 = 53), 15 prostatitis, and 26 negative biopsies. Prostate cancer T1, T2, and ADC (mean ± SD, 1660 ± 270 milliseconds, 56 ± 20 milliseconds, 0.70 × 10 ± 0.24 × 10 mm/s) were significantly lower than prostatitis (mean ± SD, 1730 ± 350 milliseconds, 77 ± 36 milliseconds, 1.00 × 10 ± 0.30 × 10 mm/s) and negative biopsies (mean ± SD, 1810 ± 250 milliseconds, 71 ± 37 milliseconds, 1.00 × 10 ± 0.33 × 10 mm/s). For cancer versus prostatitis, ADC was sensitive and T2 specific with comparable area under curve (AUC; (AUCT2 = 0.71, AUCADC = 0.79, difference between AUCs not significant P = 0.37). T1 + ADC (AUCT1 + ADC = 0.83) provided the best separation between cancer and negative biopsies. Low-grade cancer T2 and ADC (mean ± SD, 75 ± 29 milliseconds, 0.96 × 10 ± 0.34 × 10 mm/s) were significantly higher than clinically significant cancers (mean ± SD, 52 ± 16 milliseconds, 0.65 ± 0.18 × 10 mm/s), and T2 + ADC (AUCT2 + ADC = 0.91) provided the best separation.

CONCLUSIONS

T1 and T2 relaxometry combined with ADC mapping may be useful for quantitative characterization of prostate cancer grades and differentiating cancer from noncancers for PZ lesions seen on T2-weighted images.

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.

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

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

MR fingerprinting with simultaneous T1, T2, and fat signal fraction estimation with integrated B0 correction reduces bias in water T1 and T2 estimates

PURPOSE

MR fingerprinting (MRF) sequences permit efficient T₁ and T₂ estimation in cranial and extracranial regions, but these areas may include substantial fat signals that bias T₁ and T₂ estimates. MRI fat signal fraction estimation is also a topic of active research in itself, but may be complicated by B₀ heterogeneity and blurring during spiral k-space acquisitions, which are commonly used for MRF. An MRF method is proposed that separates fat and water signals, estimates water T₁ and T₂, and accounts for B₀ effects with spiral blurring correction, in a single sequence.

THEORY AND METHODS

A k-space-based fat-water separation method is further extended to unbalanced steady-state free precession MRF with swept echo time. Repeated application of this k-space fat-water separation to demodulated forms of the measured data allows a B₀ map and correction to be approximated. The method is compared with MRF without fat separation across a broad range of fat signal fractions (FSFs), water T₁s and T₂s, and under heterogeneous static fields in simulations, phantoms, and in vivo.

RESULTS

The proposed method's FSF estimates had a concordance correlation coefficient of 0.990 with conventional measurements, and reduced biases in the T₁ and T₂ estimates due to fat signal relative to other MRF sequences by several hundred ms. The B₀ correction improved the FSF, T₁, and T₂ estimation compared to those estimates without correction.

CONCLUSION

The proposed method improves MRF water T₁ and T₂ estimation in the presence of fat and provides accurate FSF estimation with inline B₀ correction.

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.

Effect of spiral undersampling patterns on FISP MRF parameter maps

PURPOSE

Artifacts arising from undersampling are not always treatable as incoherent noise for the pattern matching process in Magnetic Resonance Fingerprinting (MRF). To estimate the effect of undersampling artifacts on MRF quantitative results, spiral sampling trajectories and their temporal variation is examined.

METHODS

The effect of sampling trajectories and their variation during the MRF experiment was assessed by characterizing aliasing artifacts. Temporal rearrangements of sampling trajectories were tested and evaluated in simulations and scans of phantoms and in a volunteer brain.

RESULTS

Results show that some temporal variations of sampling patterns can lead to spatial biases in MRF parameter maps. Observed effects are consistent with derived performance indicators for different interleaving schemes, leading to substantially improved MRF sampling patterns.

CONCLUSION

With the help of the presented simulation framework, MRF implementations can be investigated and improved. This was demonstrated for a spiral FISP (Fast imaging with steady-state free precession) MRF implementation, where a significantly improved interleaving scheme was identified, and confirmed by experiment.

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 field fingerprinting

PURPOSE

To develop and evaluate the magnetic resonance field fingerprinting method that simultaneously generates T₁ , T₂ , B₀ , and B 1 + maps from a single continuous measurement.

METHODS

An encoding pattern was designed to integrate true fast imaging with steady-state precession (TrueFISP), fast imaging with steady-state precession (FISP), and fast low-angle shot (FLASH) sequence segments with varying flip angles, radio frequency (RF) phases, TEs, and gradient moments in a continuous acquisition. A multistep matching process was introduced that includes steps for integrated spiral deblurring and the correction of intravoxel phase dispersion. The method was evaluated in phantoms as well as in vivo studies in brain and lower abdomen.

RESULTS

Simultaneous measurement of T₁ , T₂ , B₀ , and B 1 + is achieved with T₁ and T₂ subsequently being less afflicted by B₀ and B 1 + variations. Phantom results demonstrate the stability of generated parameter maps. Higher undersampling factors and spatial resolution can be achieved with the proposed method as compared with solely FISP-based magnetic resonance fingerprinting. High-resolution B₀ maps can potentially be further used as diagnostic information.

CONCLUSION

The proposed magnetic resonance field fingerprinting method can estimate T₁ , T₂ , B₀ , and B 1 + maps accurately in phantoms, in the brain, and in the lower abdomen.

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.

Fast 3D magnetic resonance fingerprinting for a whole-brain coverage

PURPOSE

The purpose of this study was to accelerate the acquisition and reconstruction time of 3D magnetic resonance fingerprinting scans.

METHODS

A 3D magnetic resonance fingerprinting scan was accelerated by using a single-shot spiral trajectory with an undersampling factor of 48 in the x-y plane, and an interleaved sampling pattern with an undersampling factor of 3 through plane. Further acceleration came from reducing the waiting time between neighboring partitions. The reconstruction time was accelerated by applying singular value decomposition compression in k-space. Finally, a 3D premeasured B1 map was used to correct for the B1 inhomogeneity.

RESULTS

The T1 and T2 values of the International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology MRI phantom showed a good agreement with the standard values, with an average concordance correlation coefficient of 0.99, and coefficient of variation of 7% in the repeatability scans. The results from in vivo scans also showed high image quality in both transverse and coronal views.

CONCLUSIONS

This study applied a fast acquisition scheme for a fully quantitative 3D magnetic resonance fingerprinting scan with a total acceleration factor of 144 as compared with the Nyquist rate, such that 3D T1 , T2 , and proton density maps can be acquired with whole-brain coverage at clinical resolution in less than 5 min. Magn Reson Med 79:2190-2197, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

3D MR fingerprinting with accelerated stack-of-spirals and hybrid sliding-window and GRAPPA reconstruction

PURPOSE

Whole-brain high-resolution quantitative imaging is extremely encoding intensive, and its rapid and robust acquisition remains a challenge. Here we present a 3D MR fingerprinting (MRF) acquisition with a hybrid sliding-window (SW) and GRAPPA reconstruction strategy to obtain high-resolution T1, T2 and proton density (PD) maps with whole brain coverage in a clinically feasible timeframe.

METHODS

3D MRF data were acquired using a highly under-sampled stack-of-spirals trajectory with a steady-state precession (FISP) sequence. For data reconstruction, kx-ky under-sampling was mitigated using SW combination along the temporal axis. Non-uniform fast Fourier transform (NUFFT) was then applied to create Cartesian k-space data that are fully-sampled in the in-plane direction, and Cartesian GRAPPA was performed to resolve kz under-sampling to create an alias-free SW dataset. T1, T2 and PD maps were then obtained using dictionary matching.

RESULTS

Phantom study demonstrated that the proposed 3D-MRF acquisition/reconstruction method is able to produce quantitative maps that are consistent with conventional quantification techniques. Retrospectively under-sampled in vivo acquisition revealed that SW + GRAPPA substantially improves quantification accuracy over the current state-of-the-art accelerated 3D MRF. Prospectively under-sampled in vivo study showed that whole brain T1, T2 and PD maps with 1 mm3 resolution could be obtained in 7.5 min.

CONCLUSIONS

3D MRF stack-of-spirals acquisition with hybrid SW + GRAPPA reconstruction may provide a feasible approach for rapid, high-resolution quantitative whole-brain imaging.

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.