Three-dimensional MR Fingerprinting for Quantitative Breast Imaging

Efficient pulse sequence design framework for high-dimensional MR fingerprinting scans using systematic error index

PURPOSE

For effective optimization of MR fingerprinting (MRF) pulse sequences, estimating and minimizing errors from actual scan conditions are crucial. Although virtual-scan simulations offer an approximation to these errors, their computational demands become expensive for high-dimensional MRF frameworks, where interactions between more than two tissue properties are considered. This complexity makes sequence optimization impractical. We introduce a new mathematical model, the systematic error index (SEI), to address the scalability challenges for high-dimensional MRF sequence design.

METHODS

By eliminating the need to perform dictionary matching, the SEI model approximates quantification errors with low computational costs. The SEI model was validated in comparison with virtual-scan simulations. The SEI model was further applied to optimize three high-dimensional MRF sequences that quantify two to four tissue properties. The optimized scans were examined in simulations and healthy subjects.

RESULTS

The proposed SEI model closely approximated the virtual-scan simulation outcomes while achieving hundred- to thousand-times acceleration in the computational speed. In both simulation and in vivo experiments, the optimized MRF sequences yield higher measurement accuracy with fewer undersampling artifacts at shorter scan times than the heuristically designed sequences.

CONCLUSION

We developed an efficient method for estimating real-world errors in MRF scans with high computational efficiency. Our results illustrate that the SEI model could approximate errors both qualitatively and quantitatively. We also proved the practicality of the SEI model of optimizing sequences for high-dimensional MRF frameworks with manageable computational power. The optimized high-dimensional MRF scans exhibited enhanced robustness against undersampling and system imperfections with faster scan times.

The diagnostic performance of T1 mapping in the assessment of breast lesions: A preliminary study

PURPOSE

To assess T1 mapping performance in distinguishing between benign and malignant breast lesions and to explore its correlation with histopathologic features in breast cancer.

METHODS

This study prospectively enrolled 103 participants with a total of 108 lesions, including 25 benign and 83 malignant lesions. T1 mapping, diffusion-weighted imaging (DWI), and dynamic contrast-enhanced (DCE) were performed. Two radiologists independently outlined the ROIs and analyzed T1 and apparent diffusion coefficient (ADC) values for each lesion, assessing interobserver reliability with the intraclass correlation coefficient (ICC). T1 and ADC values were compared between benign and malignant lesions, across different histopathological characteristics (histological grades, estrogen, progesterone and HER2 receptors expression, Ki67, N status). Receiver operating characteristic (ROC) analysis and Pearson correlation coefficient (ρ) were performed.

RESULTS

T1 values showed statistically significant differences between benign and malignant groups (P < 0.001), with higher values in the malignant (1817.08 ms ± 126.64) compared to the benign group (1429.31 ms ± 167.66). In addition, T1 values significantly increased in the ER (-) group (P = 0.001). No significant differences were found in T1 values among HER2, Ki67, N status, and histological grades groups. Furthermore, T1 values exhibited a significant correlation (ρ) with ER (P < 0.01) and PR (P = 0.03). The AUC for T1 value in distinguishing benign from malignant lesions was 0.69 (95 % CI: 0.55 - 0.82, P = 0.005), and for evaluating ER status, it was 0.75 (95 % CI: 0.62 - 0.87, P = 0.002).

CONCLUSIONS

T1 mapping holds the potential as an imaging biomarker to assist in the discrimination of benign and malignant breast lesions and assessing the ER expression status in breast cancer.

Multiparametric quantification of T1 and T2 relaxation time of bone metastasis in comparison with red or fatty bone marrow using magnetic resonance fingerprinting

OBJECTIVES

To assess the T1 and T2 values of bone marrow lesions in spine and pelvis derived from magnetic resonance fingerprinting (MRF) and to evaluate the differences in values among bone metastasis, red marrow and fatty marrow.

METHODS

Sixty patients who underwent lumbar spine and pelvic MRI with magnetic resonance fingerprinting were retrospectively included. Among eligible patients, those with bone metastasis, benign red marrow deposition and normal fatty marrow were identified. Two radiologists independently measured the T1 and T2 values from metastatic bone lesions, fatty marrow, and red marrow deposition on three-dimensional-magnetic resonance fingerprinting. Intergroup comparison and interobserver agreement were analyzed.

RESULTS

T1 relaxation time was significantly higher in osteoblastic metastasis than in red marrow (1674.6 ± 436.3 vs 858.7 ± 319.5, p < .001). Intraclass correlation coefficients for T1 and T2 values were 0.96 (p < 0.001) and 0.83 (p < 0.001), respectively. T2 relaxation time of osteoblastic metastasis and red marrow deposition had no evidence of a difference (osteoblastic metastasis, 57.9 ± 25.0 vs red marrow, 58.0 ± 34.4, p = 0.45), as were the average T2 values of osteolytic metastasis and red marrow deposition (osteolytic metastasis, 45.3 ± 15.1 vs red marrow, 58.0 ± 34.4, p = 0.63).

CONCLUSIONS

We report the feasibility of three-dimensional-magnetic resonance fingerprinting based quantification of bone marrow to differentiate bone metastasis from red marrow. Simultaneous T1 and T2 quantification of metastasis and red marrow deposition was possible in spine and pelvis and showed significant different values with excellent inter-reader agreement.

ADVANCE IN KNOWLEDGE

T1 values from three-dimensional-magnetic resonance fingerprinting might be a useful quantifier for evaluating bone marrow lesions.

A Non-Contrast Multi-Parametric MRI Biomarker for Assessment of MR-Guided Focused Ultrasound Thermal Therapies

OBJECTIVE

We present the development of a non-contrast multi-parametric magnetic resonance (MPMR) imaging biomarker to assess treatment outcomes for magnetic resonance-guided focused ultrasound (MRgFUS) ablations of localized tumors. Images obtained immediately following MRgFUS ablation were inputs for voxel-wise supervised learning classifiers, trained using registered histology as a label for thermal necrosis.

METHODS

VX2 tumors in New Zealand white rabbits quadriceps were thermally ablated using an MRgFUS system under 3 T MRI guidance. Animals were re-imaged three days post-ablation and euthanized. Histological necrosis labels were created by 3D registration between MR images and digitized H&E segmentations of thermal necrosis to enable voxel-wise classification of necrosis. Supervised MPMR classifier inputs included maximum temperature rise, cumulative thermal dose (CTD), post-FUS differences in T2-weighted images, and apparent diffusion coefficient, or ADC, maps. A logistic regression, support vector machine, and random forest classifier were trained in red a leave-one-out strategy in test data from four subjects.

RESULTS

In the validation dataset, the MPMR classifiers achieved higher recall and Dice than a clinically adopted 240 cumulative equivalent minutes at 43 °C (CEM 43) threshold (0.43) in all subjects. The average Dice scores of overlap with the registered histological label for the logistic regression (0.63) and support vector machine (0.63) MPMR classifiers were within 6% of the acute contrast-enhanced non-perfused volume (0.67).

CONCLUSIONS

Voxel-wise registration of MPMR data to histological outcomes facilitated supervised learning of an accurate non-contrast MR biomarker for MRgFUS ablations in a rabbit VX2 tumor model.

Synthetic MRI in breast cancer: differentiating benign from malignant lesions and predicting immunohistochemical expression status

To evaluate and compare the performance of synthetic magnetic resonance imaging (SyMRI) in classifying benign and malignant breast lesions and predicting the expression status of immunohistochemistry (IHC) markers. We retrospectively analysed 121 patients with breast lesions who underwent dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and SyMRI before surgery in our hospital. DCE-MRI was used to assess the lesions, and then regions of interest (ROIs) were outlined on SyMRI (before and after enhancement), and apparent diffusion coefficient (ADC) maps to obtain quantitative values. After being grouped according to benign and malignant status, the malignant lesions were divided into high and low expression groups according to the expression status of IHC markers. Logistic regression was used to analyse the differences in independent variables between groups. The performance of the modalities in classification and prediction was evaluated by receiver operating characteristic (ROC) curves. In total, 57 of 121 lesions were benign, the other 64 were malignant, and 56 malignant lesions performed immunohistochemical staining. Quantitative values from proton density-weighted imaging prior to an injection of the contrast agent (PD-Pre) and T2-weighted imaging (T2WI) after the injection (T2-Gd), as well as its standard deviation (SD of T2-Gd), were valuable SyMRI parameters for the classification of benign and malignant breast lesions, but the performance of SyMRI (area under the curve, AUC = 0.716) was not as good as that of ADC values (AUC = 0.853). However, ADC values could not predict the expression status of breast cancer markers, for which SyMRI had excellent performance. The AUCs of androgen receptor (AR), estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER-2), p53 and Ki-67 were 0.687, 0.890, 0.852, 0.746, 0.813 and 0.774, respectively. SyMRI had certain value in distinguishing between benign and malignant breast lesions, and ADC values were still the ideal method. However, to predict the expression status of IHC markers, SyMRI had an incomparable value compared with ADC values.

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.

Only-train-once MR fingerprinting for B0 and B1 inhomogeneity correction in quantitative magnetization-transfer contrast

PURPOSE

To develop a fast, deep-learning approach for quantitative magnetization-transfer contrast (MTC)-MR fingerprinting (MRF) that simultaneously estimates multiple tissue parameters and corrects the effects of B0 and B1 variations.

METHODS

An only-train-once recurrent neural network was designed to perform the fast tissue-parameter quantification for a large range of different MRF acquisition schedules. It enabled a dynamic scan-wise linear calibration of the scan parameters using the measured B0 and B1 maps, which allowed accurate, multiple-tissue parameter mapping. MRF images were acquired from 8 healthy volunteers at 3 T. Estimated parameter maps from the MRF images were used to synthesize the MTC reference signal (Zref ) through Bloch equations at multiple saturation power levels.

RESULTS

The B0 and B1 errors in MR fingerprints, if not corrected, would impair the tissue quantification and subsequently corrupt the synthesized MTC reference images. Bloch equation-based numerical phantom studies and synthetic MRI analysis demonstrated that the proposed approach could correctly estimate water and semisolid macromolecule parameters, even with severe B0 and B1 inhomogeneities.

CONCLUSION

The only-train-once deep-learning framework can improve the reconstruction accuracy of brain-tissue parameter maps and be further combined with any conventional MRF or CEST-MRF method.

Assessment of MRF for simultaneous T1 and T2 quantification and water-fat separation in the liver at 0.55 T

OBJECTIVE

The goal of this work was to assess the feasibility of performing MRF in the liver on a 0.55 T scanner and to examine the feasibility of water-fat separation using rosette MRF at 0.55 T.

MATERIALS AND METHODS

Spiral and rosette MRF sequences were implemented on a commercial 0.55 T scanner. The accuracy of both sequences in T1 and T2 quantification was validated in the ISMRM/NIST system phantom. The efficacy of rosette MRF in water-fat separation was evaluated in simulations and water/oil phantoms. Both spiral and rosette MRF were performed in the liver of healthy subjects.

RESULTS

In the ISMRM/NIST phantom, both spiral and rosette MRF achieved good agreement with reference values in T1 and T2 measurements. In addition, rosette MRF enables water-fat separation and can generate water- and fat- specific T1 maps, T2 maps, and proton density images from the same dataset for a spatial resolution of 1.56 × 1.56 × 5mm3 within the acquisition time of 15 s.

CONCLUSION

It is feasible to measure T1 and T2 simultaneously in the liver using MRF on a 0.55 T system with lower performance gradients compared to state-of-the-art 1.5 T and 3 T systems within an acquisition time of 15 s. In addition, rosette MRF enables water-fat separation along with T1 and T2 quantification with no time penalty.

Optimized three-dimensional ultrashort echo time: Magnetic resonance fingerprinting for myelin tissue fraction mapping

Quantitative assessment of brain myelination has gained attention for both research and diagnosis of neurological diseases. However, conventional pulse sequences cannot directly acquire the myelin-proton signals due to its extremely short T2 and T2* values. To obtain the myelin-proton signals, dedicated short T2 acquisition techniques, such as ultrashort echo time (UTE) imaging, have been introduced. However, it remains challenging to isolate the myelin-proton signals from tissues with longer T2. In this article, we extended our previous two-dimensional ultrashort echo time magnetic resonance fingerprinting (UTE-MRF) with dual-echo acquisition to three dimensional (3D). Given a relatively low proton density (PD) of myelin-proton, we utilized Cramér-Rao Lower Bound to encode myelin-proton with the maximal SNR efficiency for optimizing the MR fingerprinting design, in order to improve the sensitivity of the sequence to myelin-proton. In addition, with a second echo of approximately 3 ms, myelin-water component can be also captured. A myelin-tissue (myelin-proton and myelin-water) fraction mapping can be thus calculated. The optimized 3D UTE-MRF with dual-echo acquisition is tested in simulations, physical phantom and in vivo studies of both healthy subjects and multiple sclerosis patients. The results suggest that the rapidly decayed myelin-proton and myelin-water signal can be depicted with UTE signals of our method at clinically relevant resolution (1.8 mm isotropic) in 15 min. With its good sensitivity to myelin loss in multiple sclerosis patients demonstrated, our method for the whole brain myelin-tissue fraction mapping in clinical friendly scan time has the potential for routine clinical imaging.

Probing lipids relaxation times in breast cancer using magnetic resonance spectroscopic fingerprinting

OBJECTIVES

To investigate the clinical relevance of the relaxation times of lipids within breast cancer and normal fibroglandular tissue in vivo, using magnetic resonance spectroscopic fingerprinting (MRSF).

METHODS

Twelve patients with biopsy-confirmed breast cancer and 14 healthy controls were prospectively scanned at 3 T using a protocol consisting of diffusion tensor imaging (DTI), MRSF, and dynamic contrast-enhanced (DCE) MRI. Single-voxel MRSF data was recorded from the tumor (patients) - identified using DTI - or normal fibroglandular tissue (controls), in under 20 s. MRSF data was analyzed using in-house software. Linear mixed model analysis was used to compare the relaxation times of lipids in breast cancer VOIs vs. normal fibroglandular tissue.

RESULTS

Seven distinguished lipid metabolite peaks were identified and their relaxation times were recorded. Of them, several exhibited statistically significant changes between controls and patients, with strong significance (p < 10^-3) recorded for several of the lipid resonances at 1.3 ppm (T₁ = 355 ± 17 ms vs. 389 ± 27 ms), 4.1 ppm (T₁ = 255 ± 86 ms vs. 127 ± 33 ms), 5.22 ppm (T₁ = 724 ± 81 ms vs. 516 ± 62 ms), and 5.31 ppm (T₂ = 56 ± 5 ms vs. 44 ± 3.5 ms, respectively).

CONCLUSIONS

The application of MRSF to breast cancer imaging is feasible and achievable in clinically relevant scan time. Further studies are required to verify and comprehend the underling biological mechanism behind the differences in lipid relaxation times in cancer and normal fibroglandular tissue.

KEY POINTS

•The relaxation times of lipids in breast tissue are potential markers for quantitative characterization of the normal fibroglandular tissue and cancer. •Lipid relaxation times can be acquired rapidly in a clinically relevant manner using a single-voxel technique, termed MRSF. •Relaxation times of T₁ at 1.3 ppm, 4.1 ppm, and 5.22 ppm, as well as of T₂ at 5.31 ppm, were significantly different between measurements within breast cancer and the normal fibroglandular tissue.

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.

Feasibility study of 2D Dixon-Magnetic Resonance Fingerprinting (MRF) of breast cancer

Purpose

Application of MRF to evaluate the feasibility of 2D Dixon blurring-corrected MRF (2DDb-cMRF) to differentiate breast cancer (BC) from normal fibroglandular tissue (FGT).

Methods

Prospective study on 14 patients with unilateral BC on 1.5 T system/axial T2w-TSE sequence, 2DDb-cMRF, B1 map, dynamic contrast-enhanced (DCE) T1-w GE-series. Mean T1 and T2 values and standard deviations were computed in the BC-/FGT-ROI on pre-/post-contrast MRF-maps and their differences were tested by two-tailed student t-test.Accuracy and repeatability of MRF were evaluated in a phantom experiment with gelatin with Primovist surrounded by fat.The T1 reduction between pre-/post-contrast MRF-maps was correlated to DCE signal enhancement in the last image post-contrast through the Pearson´s correlation coefficient (r) and for the phantom validation experiment through the Lin's concordance correlation coefficient (CCC).Visual evaluation of cancers on MRF-Maps was performed by rating each MRF-Map by 3 radiologists.

Results

T1- and T2-MRF values of BC vs. FGT were for T1 and T2 pre-contrast respectively: 1147 ± 1 ms vs. 1052 ± 9 ms (p = 0.007) and 83 ± 1 ms vs. 73 ± 1 ms (p = 0.03); post-contrast respectively: 367.3 ± 121.5 ms vs. 690.3 ± 200.3 ms (p = 0.0005) and 76.9 ± 11.5 ms vs. 69.8 ± 15.2 ms (p = 0.12). r was positive (FGT r = 0.7; BC r = 0.6). CCC was 0.999 for T1 and 0.994 for T2. In the T1- and T2-MRF-Maps before contrast respectively (7,7,8)/14 and (5,9,8)/14 cancers were visible to the readers; afterwards, (11,12,12)/14 and (5,6,11)/14.

Conclusions

MRF is promising for distinction between BC and FGT as well as for analyzing pre-/post-contrast T1 changes. However, its potential for differential diagnosis warrants further studies.

MR fingerprinting of the prostate

Multiparametric magnetic resonance imaging (mpMRI) has been adopted as the key tool for detection, localization, characterization, and risk stratification of patients suspected to have prostate cancer. Despite advantages over systematic biopsy, the interpretation of prostate mpMRI has limitations including a steep learning curve, leading to considerable interobserver variation. There is growing interest in clinical translation of quantitative imaging techniques for more objective lesion assessment. However, traditional mapping techniques are slow, precluding their use in the clinic. Magnetic resonance fingerprinting (MRF) is an efficient approach for quantitative maps of multiple tissue properties simultaneously. The T1 and T2 values obtained with MRF have been validated with phantom studies as well as in normal volunteers and patients. Studies have shown that MRF-derived T1 and T2 along with ADC values are all significant independent predictors in the differentiation between normal prostate tissue and prostate cancer, and hold promise in differentiating low and intermediate/high-grade cancers. This review seeks to introduce the basics of the prostate MRF technique, discuss the potential applications of prostate MRF for the characterization of prostate cancer, and describes ongoing areas of research.

The Impact of Dense Breasts on the Stage of Breast Cancer at Diagnosis: A Review and Options for Supplemental Screening

The purpose of breast cancer screening is to find cancers early to reduce mortality and to allow successful treatment with less aggressive therapy. Mammography is the gold standard for breast cancer screening. Its efficacy in reducing mortality from breast cancer was proven in randomized controlled trials (RCTs) conducted from the early 1960s to the mid 1990s. Panels that recommend breast cancer screening guidelines have traditionally relied on the old RCTs, which did not include considerations of breast density, race/ethnicity, current hormone therapy, and other risk factors. Women do not all benefit equally from mammography. Mortality reduction is significantly lower in women with dense breasts because normal dense tissue can mask cancers on mammograms. Moreover, women with dense breasts are known to be at increased risk. To provide equity, breast cancer screening guidelines should be created with the goal of maximizing mortality reduction and allowing less aggressive therapy, which may include decreasing the interval between screening mammograms and recommending consideration of supplemental screening for women with dense breasts. This review will address the issue of dense breasts and the impact on the stage of breast cancer at the time of diagnosis, and discuss options for supplemental screening.

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.

3D MR fingerprinting using Seiffert spirals

PURPOSE

Seiffert spirals were recently explored as an efficient way to traverse 3D k-space compared to traditional 3D techniques. Several studies have shown the ability of 3D MR fingerprinting (MRF) techniques to acquire T₁ and T₂ relaxation maps in a short period of time. However, these sequences do not sample across a large region of 3D k-space every TR, especially in the way that Seiffert trajectories can.

METHODS

A 3D MRF sequence was designed using 8 Seiffert spirals rotated in 3D k-space, with flip angle modulation for T₁ and T₂ sensitivity. The sequence was compared to an MRF sequence using a 2D spiral rotated in 3D k-space using the tiny golden angle acquisition with similar resolution/readout duration. Both sequences were evaluated using simulations, phantom validation, and in vivo imaging.

RESULTS

In all experiments, the Seiffert spiral MRF sequence performed similar to if not better than the multi-axis 2D spiral MRF sequence. Strong intraclass correlation coefficients (> 0.9) were found between conventional and MRF sequences in phantoms, whereas the in vivo results showed slightly less aliasing artifact with the Seiffert trajectory.

CONCLUSION

In this study, Seiffert spirals were used within the MRF framework to acquire high-resolution T₁ and T₂ relaxation time maps in less than 2.5 min. The reduced aliasing artifacts seen with the Seiffert sequence suggests that sampling over 3D k-space evenly each TR can improve quantification or shorten scan times.

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.

Screening Breast MRI and Gadolinium Deposition: Cause for Concern?

Gadolinium-based contrast agents (GBCAs) have been used worldwide for over 30 years and have enabled lifesaving diagnoses. Contrast-enhanced breast MRI is frequently used as supplemental screening for women with an elevated lifetime risk of breast cancer. Data have emerged that indicate a fractional amount of administered gadolinium is retained in the bone, skin, solid organs, and brain tissues of patients with normal renal function, although there are currently no reliable data regarding the clinical or biological significance of this retention. Linear GBCAs are associated with a higher risk of gadolinium retention than macrocyclic agents. Over the course of their lives, screened women may receive high cumulative doses of GBCA. Therefore, as breast MRI screening utilization increases, thoughtful use of GBCA is indicated in this patient population.

Characterization of the diffusion signal of breast tissues using multi-exponential models

PURPOSE

Restriction spectrum imaging (RSI) decomposes the diffusion-weighted MRI signal into separate components of known apparent diffusion coefficients (ADCs). The number of diffusion components and optimal ADCs for RSI are organ-specific and determined empirically. The purpose of this work was to determine the RSI model for breast tissues.

METHODS

The diffusion-weighted MRI signal was described using a linear combination of multiple exponential components. A set of ADC values was estimated to fit voxels in cancer and control ROIs. Later, the signal contributions of each diffusion component were estimated using these fixed ADC values. Relative-fitting residuals and Bayesian information criterion were assessed. Contrast-to-noise ratio between cancer and fibroglandular tissue in RSI-derived signal contribution maps was compared to DCE imaging.

RESULTS

A total of 74 women with breast cancer were scanned at 3.0 Tesla MRI. The fitting residuals of conventional ADC and Bayesian information criterion suggest that a 3-component model improves the characterization of the diffusion signal over a biexponential model. Estimated ADCs of triexponential model were D_1,3 = 0, D_2,3 = 1.5 × 10^-3 , and D_3,3 = 10.8 × 10^-3 mm² /s. The RSI-derived signal contributions of the slower diffusion components were larger in tumors than in fibroglandular tissues. Further, the contrast-to-noise and specificity at 80% sensitivity of DCE and a subset of RSI-derived maps were equivalent.

CONCLUSION

Breast diffusion-weighted MRI signal was best described using a triexponential model. Tumor conspicuity in breast RSI model is comparable to that of DCE without the use of exogenous contrast. These data may be used as differential features between healthy and malignant breast tissues.

Magnetic Resonance Fingerprinting of the Pediatric Brain

Magnetic resonance fingerprinting (MRF) is increasingly being used to evaluate brain development and differentiate normal and pathologic tissues in children. MRF can provide reliable and accurate intrinsic tissue properties, such as T1 and T2 relaxation times. MRF is a powerful tool in evaluating brain disease in pediatric population. MRF is a new quantitative MR imaging technique for rapid and simultaneous quantification of multiple tissue properties.

Technical overview of magnetic resonance fingerprinting and its applications in radiation therapy

Magnetic resonance fingerprinting (MRF) is an emerging imaging technique for rapid and simultaneous quantification of multiple tissue properties. The technique has been developed for quantitative imaging of different organs. The obtained quantitative measures have the potential to improve multiple steps of a typical radiotherapy workflow and potentially further improve integration of magnetic resonance imaging guided clinical decision making. In this review paper, we first provide a technical overview of the MRF method from data acquisition to postprocessing, along with recent development in advanced reconstruction methods. We further discuss critical aspects that could influence its usage in radiation therapy, such as accuracy and precision, repeatability and reproducibility, geometric distortion, and motion robustness. Finally, future directions for MRF application in radiation therapy are discussed.

Repeatability of MR fingerprinting in normal cervix and utility in cervical carcinoma

Background

Magnetic resonance fingerprinting (MRF) is a fast-imaging acquisition technique that generates quantitative and co-registered parametric maps. The aim of this feasibility study was to evaluate the agreement between MRF and phantom reference values, scan-rescan repeatability of MRF in normal cervix, and its ability to distinguish cervical carcinoma (CC) from normal cervical tissues.

Methods

An International Society of Magnetic Resonance in Medicine/National Institute of Standards and Technology (ISMRM/NIST) phantom was scanned using MRF 15 times over 65 days. Agreement between MRF and phantom reference T1 and T2 values was assessed by linear regression. Healthy volunteers and patients with suspected CC were prospectively recruited. MRF was repeated twice for healthy volunteers (MRF1 and MRF2). Volumes of interest of normal cervical tissues and CC were delineated on T1 and T2 maps. MRF scan-rescan repeatability was evaluated by Bland-Altman plots, within-subject coefficients of variation (wCV), and intraclass correlation coefficients (ICC). T1 and T2 values were compared between CC and normal cervical tissues using Mann-Whitney U test. Receiver operating characteristic (ROC) analysis was performed to evaluate diagnostic efficiency.

Results

Strong correlations were observed between MRF and phantom (R²=0.999 for T1, 0.981 for T2). Twelve healthy volunteers (28.7±5.1 years) and 28 patients with CC (54.6±15.2 years) were recruited for the in-vivo experiments. Repeatability of MRF parameters were wCV <3% for T1, <5% for T2 and ICC ≥0.92 for T1, ≥0.94 for T2. T1 value of CC (1,529±112 ms) was higher than normal mucosa [MRF1: 1,430±129 ms, MRF2: 1,440±130 ms; P=0.031, area under the curve (AUC) ≥0.717] and normal stroma (MRF1: 1,258±101 ms, MRF2: 1,276±105 ms; P<0.001, AUC ≥0.946). T2 value of CC (69±9 ms) was lower than normal mucosa (MRF1: 88±16 ms, MRF2: 87±13 ms; P<0.001, AUC ≥0.854), but was not different from normal stroma (P=0.919).

Conclusions

Excellent agreement was observed between MRF and phantom reference values. MRF exhibited excellent scan-rescan repeatability in normal cervix with potential value in differentiating CC from normal cervical tissues.

Workload of diagnostic radiologists in the foreseeable future based on recent scientific advances: growth expectations and role of artificial intelligence

Objective

To determine the anticipated contribution of recently published medical imaging literature, including artificial intelligence (AI), on the workload of diagnostic radiologists.

Methods

This study included a random sample of 440 medical imaging studies published in 2019. The direct contribution of each study to patient care and its effect on the workload of diagnostic radiologists (i.e., number of examinations performed per time unit) was assessed. Separate analyses were done for an academic tertiary care center and a non-academic general teaching hospital.

Results

In the academic tertiary care center setting, 65.0% (286/440) of studies could directly contribute to patient care, of which 48.3% (138/286) would increase workload, 46.2% (132/286) would not change workload, 4.5% (13/286) would decrease workload, and 1.0% (3/286) had an unclear effect on workload. In the non-academic general teaching hospital setting, 63.0% (277/240) of studies could directly contribute to patient care, of which 48.7% (135/277) would increase workload, 46.2% (128/277) would not change workload, 4.3% (12/277) would decrease workload, and 0.7% (2/277) had an unclear effect on workload. Studies with AI as primary research area were significantly associated with an increased workload (p < 0.001), with an odds ratio (OR) of 10.64 (95% confidence interval (CI) 3.25–34.80) in the academic tertiary care center setting and an OR of 10.45 (95% CI 3.19–34.21) in the non-academic general teaching hospital setting.

Conclusions

Recently published medical imaging studies often add value to radiological patient care. However, they likely increase the overall workload of diagnostic radiologists, and this particularly applies to AI studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13244-021-01031-4.

A multi-inversion-recovery magnetic resonance fingerprinting for multi-compartment water mapping

PURPOSE

This study aimed at introducing short-T1/T2 compartment to MR fingerprinting (MRF) at 3 T. Water that is bound to myelin macromolecules have significantly shorter T1 and T2 than free water and can be distinguished from free water by multi-compartment analysis.

METHODS

We developed a new multi-inversion-recovery (mIR) water mapping-MRF based on an unbalanced steady-state coherent sequence (FISP). mIR pulses with an interval of 400 or 500 repetition times (TRs) were inserted into the conventional FISP MRF sequence. Data from our proposed mIR MRF was used to quantify different compartments, including myelin water, gray matter free water, and white matter free water, of brain water by virtue of the iterative non-negative least square (NNLS) with reweighting. Three healthy volunteers were scanned with mIR MRF on a clinical 3 T MRI.

RESULTS

Using an extended phase graph simulation, we found that our proposed mIR scheme with four IR pulses allowed differentiation between short and long T1/T2 components. For in vivo experiments, we achieved the quantification of myelin water, gray matter water, and white matter water at an image resolution of 1.17 × 1.17 × 5 mm³/pixel. As compared to the conventional MRF technique with single IR, our proposed mIR improved the detection of myelin water content. In addition, mIR MRF using spiral-in/out trajectory provided a higher signal level compared with that with spiral-out trajectory. Myelin water quantification using mIR MRF with 4 IR and 5 IR pulses were qualitatively similar. Meanwhile, 5 IR MRF showed fewer artifacts in myelin water detection.

CONCLUSION

We developed a new mIR MRF sequence for the rapid quantification of brain water compartments.

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.

Magnetic resonance fingerprinting: an overview

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

Free-Breathing Abdominal Magnetic Resonance Fingerprinting Using a Pilot Tone Navigator

BACKGROUND

Quantitative T1 and T2 mapping in the abdomen provides valuable information in tissue characterization but is technically challenging due to respiratory motions. The proposed technique integrates magnetic resonance fingerprinting (MRF) and pilot tone (PT) navigator with retrospective gating to provide simultaneous quantification of multiple tissue properties in a single acquisition without breath-holding or patient set-up.

PURPOSE

To develop a free-breathing abdominal MRF technique for quantitative mapping in the abdomen.

STUDY TYPE

Prospective.

POPULATION

Twelve healthy volunteers.

FIELD STRENGTH/SEQUENCE

A 3 T, two-dimensional (2D) and three-dimensional (3D) spiral MRF sequence with fast imaging with steady-state free precession (FISP) readout.

ASSESSMENT

The PT navigator was compared to standard respiratory belt performance. The T1 and T2 values acquired using 2D and 3D MRF with and without PT were obtained in a phantom and compared to reference values. Digital phantom simulation was performed to evaluate PT MRF reconstruction with varying breathing patterns. In the in vivo studies, T1 and T2 values derived from PT 2D MRF were compared to 2D breath-hold MRF. T1 and T2 values derived from PT 3D MRF were compared to published values.

STATISTICAL TESTS

Principal component analysis (PCA), linear regression, relative error, Pearson correlation, paired Student's t-test, Bland-Altman Analysis.

RESULTS

The phantom study showed PT MRF T1 values had a mean difference of 0.2% ± 0.1%, and T2 values had a mean difference of 0.1% ± 0.4% when compared to no-PT MRF values. The digital phantom experiment suggested the T1 and T2 maps at both end-exhalation and end-inhalation states resemble the corresponding ground-truth maps.

DATA CONCLUSION

The phantom study showed good agreement between MRF T1 and T2 values and with reference values. In vivo studies demonstrated that 2D and 3D quantitative imaging in the abdomen could be achieved with integration of PT navigation with MRF reconstruction using retrospective gating of respiratory motion. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

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

Purpose

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

Materials and methods

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

Results

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

Conclusions

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

Diffusion-weighted double-echo steady-state with a three-dimensional cones trajectory for non-contrast-enhanced breast MRI

The image quality limitations of echo-planar diffusion-weighted imaging (DWI) are an obstacle to its widespread adoption in the breast. Steady-state DWI is an alternative DWI method with more robust image quality but its contrast for imaging breast cancer is not well-understood. The aim of this study was to develop and evaluate diffusion-weighted double-echo steady-state imaging with a three-dimensional cones trajectory (DW-DESS-Cones) as an alternative to conventional DWI for non-contrast-enhanced MRI in the breast. This prospective study included 28 women undergoing clinically indicated breast MRI and six asymptomatic volunteers. In vivo studies were performed at 3 T and included DW-DESS-Cones, DW-DESS-Cartesian, DWI, and CE-MRI acquisitions. Phantom experiments (diffusion phantom, High Precision Devices) and simulations were performed to establish framework for contrast of DW-DESS-Cones in comparison to DWI in the breast. Motion artifacts of DW-DESS-Cones were measured with artifact-to-noise ratio in volunteers and patients. Lesion-to-fibroglandular tissue signal ratios were measured, lesions were categorized as hyperintense or hypointense, and an image quality observer study was performed in DW-DESS-Cones and DWI in patients. Effect of DW-DESS-Cones method on motion artifacts was tested by mixed-effects generalized linear model. Effect of DW-DESS-Cones on signal in phantom was tested by quadratic regression. Correlation was calculated between DW-DESS-Cones and DWI lesion-to-fibroglandular tissue signal ratios. Inter-observer agreement was assessed with Gwet's AC. Simulations predicted hyperintensity of lesions with DW-DESS-Cones but at a 3% to 67% lower degree than with DWI. Motion artifacts were reduced with DW-DESS-Cones versus DW-DESS-Cartesian (p < 0.05). Lesion-to-fibroglandular tissue signal ratios were not correlated between DW-DESS-Cones and DWI (r = 0.25, p = 0.38). Concordant hyperintensity/hypointensity was observed between DW-DESS-Cones and DWI in 11/14 lesions. DW-DESS-Cones improved sharpness, distortion, and overall image quality versus DWI. DW-DESS-Cones may be able to eliminate motion artifacts in the breast allowing for investigation of higher degrees of steady-state diffusion weighting. Malignant breast lesions in DW-DESS-Cones demonstrated hyperintensity with respect to surrounding tissue without an injection of contrast. LEVEL OF EVIDENCE: 2. TECHNICAL EFFICACY STAGE: 1.

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.

Five-dimensional quantitative low-dose Multitasking dynamic contrast- enhanced MRI: Preliminary study on breast cancer

PURPOSE

To develop a low-dose Multitasking DCE technique (LD-MT-DCE) for breast imaging, enabling dynamic T₁ mapping-based quantitative characterization of tumor blood flow and vascular properties with whole-breast coverage, a spatial resolution of 0.9 × 0.9 × 1.1 mm³ , and a temporal resolution of 1.4 seconds using a 20% gadolinium dose (0.02 mmol/kg).

METHODS

Magnetic resonance Multitasking was used to reconstruct 5D images with three spatial dimensions, one T₁ recovery dimension for dynamic T₁ quantification, and one DCE dimension for contrast kinetics. Kinetic parameters F p , v p , K trans , and v e were estimated from dynamic T₁ maps using the two-compartment exchange model. The LD-MT-DCE repeatability and agreement against standard-dose MT-DCE were evaluated in 20 healthy subjects. In 7 patients with triple-negative breast cancer, LD-MT-DCE image quality and diagnostic results were compared with that of standard-dose clinical DCE in the same imaging session. One-way unbalanced analysis of variance with Tukey test was performed to evaluate the statistical significance of the kinetic parameters between control and patient groups.

RESULTS

The LD-MT-DCE technique was repeatable, agreed with standard-dose MT-DCE, and showed excellent image quality. The diagnosis using LD-MT-DCE matched well with clinical results. The values of F p , v p , and K trans were significantly different between malignant tumors and normal breast tissue (P < .001, < .001, and < .001, respectively), and between malignant and benign tumors (P = .020, .003, and < .001, respectively).

CONCLUSION

The LD-MT-DCE technique was repeatable and showed excellent image quality and equivalent diagnosis compared with standard-dose clinical DCE. The estimated kinetic parameters were capable of differentiating between normal breast tissue and benign and malignant tumors.

The diagnostic performance of quantitative mapping in breast cancer patients: a preliminary study using synthetic MRI

Background

Previous studies have indicated that quantitative MRI (qMR) is beneficial for diagnosis of breast cancer. As a novel qMR technology, synthetic MRI (syMRI) may be advantageous by offering simultaneous generation of T1 and T2 mapping in one scan within a few minutes and without concern to the deposition of the gadolinium contrast agent in cell nucleus. In this study, the potential of quantitative mapping derived from Synthetic MRI (SyMRI) to diagnose breast cancer was investigated.

Methods

From April 2018 to May 2019, a total of 87 patients with suspicious breast lesions underwent both conventional and SyMRI before treatment. The quantitative metrics derived from SyMRI, including T1 and T2 values, were measured in breast lesions. The diagnostic performance of SyMRI was evaluated with unpaired Student’s t-tests, receiver operating characteristic curve analysis and multivariate logistic regression analysis. The AUCs of quantitative values were compared using Delong test.

Results

Among 77 patients who met the inclusion criteria, 48 were diagnosed with histopathological confirmed breast cancers, and the rest had benign lesions. The breast cancers showed significantly higher T1 (1611.61 ± 215.88 ms) values and lower T2 (80.93 ± 7.51 ms) values than benign lesions. The area under the ROC curve (AUC) values were 0.931 (95% CI: 0.874–0.989) and 0.883 (95% CI: 0.810–0.956) for T1 and T2 maps, respectively, in diagnostic discrimination between breast cancers and benign lesions. A slightly increased AUC of 0.978 (95% CI: 0.915–0.993) was achieved by combining those two relaxation-based quantitative metrics.

Conclusion

In conclusion, our preliminary study showed that the quantitative T1 and T2 values obtained by SyMRI could distinguish effectively between benign and malignant breast lesions, and T1 relaxation time showed the highest diagnostic efficiency. Furthermore, combining the two quantitative relaxation metrics further improved their diagnostic performance.

Magnetic resonance fingerprinting: from evolution to clinical applications

In 2013, Magnetic Resonance Fingerprinting (MRF) emerged as a method for fast, quantitative Magnetic Resonance Imaging. This paper reviews the current status of MRF up to early 2020 and aims to highlight the advantages MRF can offer medical imaging professionals. By acquiring scan data as pseudorandom samples, MRF elicits a unique signal evolution, or 'fingerprint', from each tissue type. It matches 'randomised' free induction decay acquisitions against pre-computed simulated tissue responses to generate a set of quantitative images of T1 , T2 and proton density (PD) with co-registered voxels, rather than as traditional relative T1 - and T2 -weighted images. MRF numeric pixel values retain accuracy and reproducibility between 2% and 8%. MRF acquisition is robust to strong undersampling of k-space. Scan sequences have been optimised to suppress sub-sampling artefacts, while artificial intelligence and machine learning techniques have been employed to increase matching speed and precision. MRF promises improved patient comfort with reduced scan times and fewer image artefacts. Quantitative MRF data could be used to define population-wide numeric biomarkers that classify normal versus diseased tissue. Certification of clinical centres for MRF scan repeatability would permit numeric comparison of sequential images for any individual patient and the pooling of multiple patient images across large, cross-site imaging studies. MRF has to date shown promising results in early clinical trials, demonstrating reliable differentiation between malignant and benign prostate conditions, and normal and sclerotic hippocampal tissue. MRF is now undergoing small-scale trials at several sites across the world; moving it closer to routine clinical application.

Investigation of Synthetic Relaxometry and Diffusion Measures in the Differentiation of Benign and Malignant Breast Lesions as Compared to BI-RADS

BACKGROUND

Breast cancer is the most common malignant tumor in women and a quantitative contrast-free method is highly desirable for its diagnosis.

PURPOSE

To investigate the performance of quantitative MRI in differentiating malignant from benign breast lesions and to compare with the Breast Imaging Reporting and Data System (BI-RADS).

STUDY TYPE

Retrospective.

SUBJECTS

Eighty patients (56 with malignant lesions and 24 with benign lesions).

FIELD STRENGTH/SEQUENCE

Diffusion-weighted imaging (DWI) with a single-shot echo planar sequence and synthetic MRI with magnetic resonance image compilation (MAGiC) were performed at 3T.

ASSESSMENT

T₁ relaxation time (T₁ ), T₂ relaxation time (T₂ ), and proton density (PD) from synthetic MRI and apparent diffusion coefficient (ADC) from DWI were analyzed by two radiologists (Reader A, Reader B). Univariable and multivariable models were developed to optimize differentiation between malignant and benign lesions and their performances compared to BI-RADS.

STATISTICAL TESTS

The diagnostic performance was evaluated using multivariate logistic regression analysis and area under the receiver operating characteristic (ROC) curves (AUC).

RESULTS

T₂ , PD, and ADC values for malignant lesions were significantly lower than those in benign breast lesions for both radiologists (all P < 0.05). The combined T₂ , PD, and ADC model had the best performance for differentiating malignant and benign lesions with AUC, sensitivity, specificity, positive predictive value, and negative predictive values of 0.904, 94.6%, 87.5%, 94.6%, and 87.5%, respectively. The corresponding results for BI-RADS were no AUC, 94.6%, 75.0%, 89.8%, and 85.7%, respectively.

DATA CONCLUSION

The approach that combined synthetic MRI and DWI outperformed BI-RADS in the differential diagnosis of malignant and benign breast lesions and was achieved without contrast agents. This approach may serve as an alternative and effective strategy for the improvement of breast lesion differentiation.

LEVEL OF EVIDENCE

3.

TECHNICAL EFFICACY STAGE

3.

Technical feasibility of magnetic resonance fingerprinting on a 1.5T MRI-linac

Hybrid MRI-linac (MRL) systems enable daily multiparametric quantitative MRI to assess tumor response to radiotherapy. Magnetic resonance fingerprinting (MRF) may provide time efficient means of rapid multiparametric quantitative MRI. The accuracy of MRF, however, relies on adequate control over system imperfections, such as eddy currents and [Formula: see text], which are different and not as well established on MRL systems compared to diagnostic systems. In this study we investigate the technical feasibility of gradient spoiled 2D MRF on a 1.5T MRL. We show with phantom experiments that the MRL generates reliable MRF signals that are temporally stable during the day and have good agreement with spin-echo reference measurements. Subsequent in-vivo MRF scans in healthy volunteers and a patient with a colorectal liver metastasis showed good image quality, where the quantitative values of selected organs corresponded with the values reported in literature. Therefore we conclude that gradient spoiled 2D MRF is feasible on a 1.5T MRL with similar performance as on a diagnostic system. The precision and accuracy of the parametric maps are sufficient for further investigation of the clinical utility of MRF for online quantitatively MRI-guided radiotherapy.

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.

Three-dimensional whole-brain simultaneous T1, T2, and T1ρ quantification using MR Multitasking: Method and initial clinical experience in tissue characterization of multiple sclerosis

PURPOSE

To develop a 3D whole-brain simultaneous T1/T2/T1ρ quantification method with MR Multitasking that provides high quality, co-registered multiparametric maps in 9 min.

METHODS

MR Multitasking conceptualizes T1/T2/T1ρ relaxations as different time dimensions, simultaneously resolving all three dimensions with a low-rank tensor image model. The proposed method was validated on a phantom and in healthy volunteers, comparing quantitative measurements against corresponding reference methods and evaluating the scan-rescan repeatability. Initial clinical validation was performed in age-matched relapsing-remitting multiple sclerosis (RRMS) patients to examine the feasibility of quantitative tissue characterization and to compare with the healthy control cohort. The feasibility of synthesizing six contrast-weighted images was also examined.

RESULTS

Our framework produced high quality, co-registered T1/T2/T1ρ maps that closely resemble the reference maps. Multitasking T1/T2/T1ρ measurements showed substantial agreement with reference measurements on the phantom and in healthy controls. Bland-Altman analysis indicated good in vivo repeatability of all three parameters. In RRMS patients, lesions were conspicuously delineated on all three maps and on four synthetic weighted images (T2-weighted, T2-FLAIR, double inversion recovery, and a novel "T1ρ-FLAIR" contrast). T1 and T2 showed significant differences for normal appearing white matter between patients and controls, while T1ρ showed significant differences for normal appearing white matter, cortical gray matter, and deep gray matter. The combination of three parameters significantly improved the differentiation between RRMS patients and healthy controls, compared to using any single parameter alone.

CONCLUSION

MR Multitasking simultaneously quantifies whole-brain T1/T2/T1ρ and is clinically promising for quantitative tissue characterization of neurological diseases, such as MS.

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.

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.

Current and Emerging Magnetic Resonance-Based Techniques for Breast Cancer

Breast cancer is the most commonly diagnosed cancer among women worldwide, and early detection remains a principal factor for improved patient outcomes and reduced mortality. Clinically, magnetic resonance imaging (MRI) techniques are routinely used in determining benign and malignant tumor phenotypes and for monitoring treatment outcomes. Static MRI techniques enable superior structural contrast between adipose and fibroglandular tissues, while dynamic MRI techniques can elucidate functional characteristics of malignant tumors. The preferred clinical procedure-dynamic contrast-enhanced MRI-illuminates the hypervascularity of breast tumors through a gadolinium-based contrast agent; however, accumulation of the potentially toxic contrast agent remains a major limitation of the technique, propelling MRI research toward finding an alternative, noninvasive method. Three such techniques are magnetic resonance spectroscopy, chemical exchange saturation transfer, and non-contrast diffusion weighted imaging. These methods shed light on underlying chemical composition, provide snapshots of tissue metabolism, and more pronouncedly characterize microstructural heterogeneity. This review article outlines the present state of clinical MRI for breast cancer and examines several research techniques that demonstrate capacity for clinical translation. Ultimately, multi-parametric MRI-incorporating one or more of these emerging methods-presently holds the best potential to afford improved specificity and deliver excellent accuracy to clinics for the prediction, detection, and monitoring of breast cancer.

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.

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.

Initial assessment of 3D magnetic resonance fingerprinting (MRF) towards quantitative brain imaging for radiation therapy

PURPOSE

Magnetic resonance fingerprinting (MRF) provides quantitative T1/T2 maps, enabling applications in clinical radiotherapy such as large-scale, multi-center clinical trials for longitudinal assessment of therapy response. We evaluated the feasibility of a quantitative three-dimensional-MRF (3D-MRF) towards its radiotherapy applications of primary brain tumors.

METHODS

A fast whole-brain 3D-MRF sequence initially developed for diagnostic radiology was optimized using flexible body coils, which is the typical MR imaging setup for radiotherapy treatment planning and for MR imaging (MRI)-guided treatment delivery. Optimization criteria included the accuracy and the precision of T1/T2 quantifications of polyvinylpyrrolidone (PVP) solutions, compared to those from the 3D-MRF using a 32-channel head coil. The accuracy of T1/T2 quantifications from the optimized MRF was first examined in healthy volunteers with two different coil setups. The intra- and inter-scanner variations of image intensity from the optimized sequence were quantified by longitudinal scans of the PVP solutions on two 3T scanners. Using a 3D-printed MRI geometry phantom, susceptibility-induced distortion with the optimized 3D-MRF was quantified as the Dice coefficient of phantom contours, compared to those from CT images. By introducing intentional head motion during 10% of the scan, the robustness of the optimized 3D-MRF towards motion was evaluated through visual inspection of motion artifacts and through quantitative analysis of image sharpness in brain MRF maps.

RESULTS

The optimized sequence acquired whole-brain T1, T2 and proton density maps and with a resolution of 1.2 × 1.2 × 3 mm³ in 10 min, similar to the total acquisition time of 3D T1- and T2-weighted images of the same resolution. In vivo T1 and T2 values of the white and gray matter were consistent with literature. The intra- and inter-scanner variability of the intensity-normalized MRF T1 was 1.0% ± 0.7% and 2.3% ± 1.0% respectively, in contrast to 5.3% ± 3.8% and 3.2% ± 1.6% from the normalized T1-weighted MRI. Repeatability and reproducibility of MRF T1 were independent of intensity normalization. Both phantom and human data demonstrated that the optimized 3D-MRF is more robust to subject motion and artifacts from subject-specific susceptibility difference. Compared to CT contours, the Dice coefficient of phantom contours from 3D-MRF was 0.93, improved from 0.87 from the T1-weighted MRI.

CONCLUSION

Compared to conventional MRI, the optimized 3D-MRF demonstrated improved repeatability across time points and reproducibility across scanners for better tissue quantification, as well as improved robustness to subject-specific susceptibility and motion artifacts under a typical MR imaging setup for radiotherapy. More importantly, quantitative MRF T1/T2 measurements lead to promising potentials towards longitudinal quantitative assessment of treatment response for better adaptive therapy and for large-scale, multi-center clinical trials.

Three-dimensional simultaneous brain T1 , T2 , and ADC mapping with MR Multitasking

PURPOSE

To develop a simultaneous T₁ , T₂ , and ADC mapping method that provides co-registered, distortion-free images and enables multiparametric quantification of 3D brain coverage in a clinically feasible scan time with the MR Multitasking framework.

METHODS

The T₁ /T₂ /diffusion weighting was generated by a series of T₂ preparations and diffusion preparations. The underlying multidimensional image containing 3 spatial dimensions, 1 T₁ weighting dimension, 1 T₂ -preparation duration dimension, 1 b-value dimension, and 1 diffusion direction dimension was modeled as a 5-way low-rank tensor. A separate real-time low-rank model incorporating time-resolved phase correction was also used to compensate for both inter-shot and intra-shot phase inconsistency induced by physiological motion. The proposed method was validated on both phantom and 16 healthy subjects. The quantification of T₁ /T₂ /ADC was evaluated for each case. Three post-surgery brain tumor patients were scanned for demonstration of clinical feasibility.

RESULTS

Multitasking T₁ /T₂ /ADC maps were perfectly co-registered and free from image distortion. Phantom studies showed substantial quantitative agreement ( R 2 = 0.999 ) with reference protocols for T₁ /T₂ /ADC. In vivo studies showed nonsignificant T₁ (P = .248), T₂ (P = .97), ADC (P = .328) differences among the frontal, parietal, and occipital regions. Although Multitasking showed significant differences of T₁ (P = .03), T₂ (P < .001), and ADC (P = .001) biases against the references, the mean bias estimates were small (ΔT₁ % < 5%, ΔT₂ % < 7%, ΔADC% < 5%), with all intraclass correlation coefficients greater than 0.82 indicating "excellent" agreement. Patient studies showed that Multitasking T₁ /T₂ /ADC maps were consistent with the clinical qualitative images.

CONCLUSION

The Multitasking approach simultaneously quantifies T₁ /T₂ /ADC with substantial agreement with the references and is promising for clinical applications.

High-resolution 3D MR Fingerprinting using parallel imaging and deep learning

MR Fingerprinting (MRF) is a relatively new imaging framework capable of providing accurate and simultaneous quantification of multiple tissue properties for improved tissue characterization and disease diagnosis. While 2D MRF has been widely available, extending the method to 3D MRF has been an actively pursued area of research as a 3D approach can provide a higher spatial resolution and better tissue characterization with an inherently higher signal-to-noise ratio. However, 3D MRF with a high spatial resolution requires lengthy acquisition times, especially for a large volume, making it impractical for most clinical applications. In this study, a high-resolution 3D MR Fingerprinting technique, combining parallel imaging and deep learning, was developed for rapid and simultaneous quantification of T₁ and T₂ relaxation times. Parallel imaging was first applied along the partition-encoding direction to reduce the amount of acquired data. An advanced convolutional neural network was then integrated with the MRF framework to extract features from the MRF signal evolution for improved tissue characterization and accelerated mapping. A modified 3D-MRF sequence was also developed in the study to acquire data to train the deep learning model that can be directly applied to prospectively accelerate 3D MRF scans. Our results of quantitative T₁ and T₂ maps demonstrate that improved tissue characterization can be achieved using the proposed method as compared to prior methods. With the integration of parallel imaging and deep learning techniques, whole-brain (26 × 26 × 18 cm³) quantitative T₁ and T₂ mapping with 1-mm isotropic resolution were achieved in ~7 min. In addition, a ~7-fold improvement in processing time to extract tissue properties was also accomplished with the deep learning approach as compared to the standard template matching method. All of these improvements make high-resolution whole-brain quantitative MR imaging feasible for clinical applications.