Extended phase graphs: Dephasing, RF pulses, and echoes - pure and simple: Extended Phase Graphs

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

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

A New Method for Cartilage Evaluation in Femoroacetabular Impingement Using Quantitative T2 Magnetic Resonance Imaging: Preliminary Validation against Arthroscopic Findings

OBJECTIVE

The outcome of arthroscopic treatment for femoroacetabular impingement (FAI) depends on the preoperative status of the hip cartilage. Quantitative T2 can detect early biochemical cartilage changes, but its routine implementation is challenging. Furthermore, intrinsic T2 variability between patients makes it difficult to define a threshold to identify cartilage lesions. To address this, we propose a normalized T2-index as a new method to evaluate cartilage in FAI.

DESIGN

We retrospectively analyzed magnetic resonance imaging (MRI) data of 18 FAI patients with arthroscopically confirmed cartilage defects. Cartilage T2 maps were reconstructed from multi-spin-echo 3-T data using the echo-modulation-curve (EMC) model-based technique. The central femoral cartilage, assumed healthy in early-stage FAI, was used as the normalization reference to define a T2-index. We investigated the ability of the T2-index to detect surgically confirmed cartilage lesions.

RESULTS

The average T2-index was 1.14 ± 0.1 and 1.13 ± 0.1 for 2 separated segmentations. Using T2-index >1 as the threshold for damaged cartilage, accuracy was 88% and 100% for the 2 segmentations. We found moderate intraobserver repeatability, although separate segmentations yielded comparable accuracy. Damaged cartilage could not be identified using nonnormalized average T2 values.

CONCLUSIONS

This preliminary study confirms the importance of normalizing T2 values to account for interpatient variability and suggests that the T2-index is a promising biomarker for the detection of cartilage lesions in FAI. Future work is needed to confirm that combining T2-index with morphologic MRI and other quantitative biomarkers could improve cartilage assessment in FAI.

B1 inhomogeneity correction of RARE MRI at low SNR: Quantitative in vivo 19 F MRI of mouse neuroinflammation with a cryogenically-cooled transceive surface radiofrequency probe

PURPOSE

Low SNR in fluorine-19 (¹⁹ F) MRI benefits from cryogenically-cooled transceive surface RF probes (CRPs), but strong B₁ inhomogeneities hinder quantification. Rapid acquisition with refocused echoes (RARE) is an SNR-efficient method for MRI of neuroinflammation with perfluorinated compounds but lacks an analytical signal intensity equation to retrospectively correct B₁ inhomogeneity. Here, a workflow was proposed and validated to correct and quantify ¹⁹ F-MR signals from the inflamed mouse brain using a ¹⁹ F-CRP.

METHODS

In vivo ¹⁹ F-MR images were acquired in a neuroinflammation mouse model with a quadrature ¹⁹ F-CRP using an imaging setup including 3D-printed components to acquire co-localized anatomical and ¹⁹ F images. Model-based corrections were validated on a uniform ¹⁹ F phantom and in the neuroinflammatory model. Corrected ¹⁹ F-MR images were benchmarked against reference images and overlaid on in vivo ¹ H-MR images. Computed concentration uncertainty maps using Monte Carlo simulations served as a measure of performance of the B₁ corrections.

RESULTS

Our study reports on the first quantitative in vivo ¹⁹ F-MR images of an inflamed mouse brain using a ¹⁹ F-CRP, including in vivo T₁ calculations for ¹⁹ F-nanoparticles during pathology and B₁ corrections for ¹⁹ F-signal quantification. Model-based corrections markedly improved ¹⁹ F-signal quantification from errors > 50% to < 10% in a uniform phantom (p < 0.001). Concentration uncertainty maps ex vivo and in vivo yielded uncertainties that were generally < 25%. Monte Carlo simulations prescribed SNR ≥ 10.1 to reduce uncertainties < 10%, and SNR ≥ 4.25 to achieve uncertainties < 25%.

CONCLUSION

Our model-based correction method facilitated ¹⁹ F signal quantification in the inflamed mouse brain when using the SNR-boosting ¹⁹ F-CRP technology, paving the way for future low-SNR ¹⁹ F-MRI applications in vivo.

Simultaneous T1 , T2 , and T1ρ cardiac magnetic resonance fingerprinting for contrast agent-free myocardial tissue characterization

PURPOSE

To develop a simultaneous T₁ , T₂ , and T_1ρ cardiac magnetic resonance fingerprinting (MRF) approach to enable comprehensive contrast agent-free myocardial tissue characterization in a single breath-hold scan.

METHODS

A 2D gradient-echo electrocardiogram-triggered cardiac MRF sequence with low flip angles, varying magnetization preparation, and spiral trajectory was acquired at 1.5 T to encode T₁ , T₂ , and T_1⍴ simultaneously. The MRF images were reconstructed using low-rank inversion, regularized with a multicontrast patch-based higher-order reconstruction. Parametric maps were generated and matched in the singular value domain to extended phase graph-based dictionaries. The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional methods in terms of coefficients of determination and best fits for the phantom study, and in terms of Bland-Altman agreement, average values and coefficient of variation of T₁ , T₂ , and T_1⍴ for the healthy subjects study.

RESULTS

The T₁ , T₂ , and T_1⍴ MRF values showed excellent correlation with conventional spin-echo and clinical mapping methods in phantom studies (r² > 0.97). Measured MRF values in myocardial tissue (mean ± SD) were 1133 ± 33 ms, 38.8 ± 3.5 ms, and 52.0 ± 4.0 ms for T₁ , T₂ and T1_⍴ , respectively, against 1053 ± 47 ms, 50.4 ± 3.9 ms, and 55.9 ± 3.3 ms for T₁ modified Look-Locker inversion imaging, T₂ gradient and spin echo, and T_1⍴ turbo field echo, respectively.

CONCLUSION

A cardiac MRF approach for simultaneous quantification of myocardial T₁ , T₂ , and T_1ρ in a single breath-hold MR scan of about 16 seconds has been proposed. The approach has been investigated in phantoms and healthy subjects showing good agreement with reference spin echo measurements and conventional clinical maps.

Simultaneous comprehensive liver T1 , T2 , T 2 ∗ , T1ρ , and fat fraction characterization with MR fingerprinting

PURPOSE

To develop a novel simultaneous co-registered T₁ , T₂ , T 2 ∗ , T_1ρ , and fat fraction abdominal MR fingerprinting (MRF) approach for fully comprehensive liver-tissue characterization in a single breath-hold scan.

METHODS

A gradient-echo liver MRF sequence with low fixed flip angle, multi-echo radial readout, and varying magnetization preparation pulses for multiparametric encoding is performed at 1.5 T. The T 2 ∗ and fat fraction are estimated from a graph/cut water/fat separation method using a six-peak fat model. Water/fat singular images obtained are then matched to an MRF dictionary, estimating water-specific T₁ , T₂ , and T_1ρ . The proposed approach was tested in phantoms and 10 healthy subjects and compared against conventional sequences.

RESULTS

For the phantom studies, linear fits show excellent coefficients of determination (r² > 0.9) for every parametric map. For in vivo studies, the average values measured within regions of interest drawn on liver, spleen, muscle, and fat are statistically different from the reference scans (p < 0.05) for T₁ , T₂ , and T_1⍴ but not for T 2 ∗ and fat fraction, whereas correlation between MRF and reference scans is excellent for each parameter (r² > 0.92 for every parameter).

CONCLUSION

The proposed multi-echo inversion-recovery, T₂ , and T_1⍴ prepared liver MRF sequence presented in this work allows for quantitative T₁ , T₂ , T 2 ∗ , T_1⍴ , and fat fraction liver-tissue characterization in a single breath-hold scan of 18 seconds. The approach showed good agreement and correlation with respect to reference clinical maps.

A strategy to reduce the sensitivity of inhomogeneous magnetization transfer (ihMT) imaging to radiofrequency transmit field variations at 3 T

PURPOSE

To minimize the sensitivity of inhomogeneous magnetization transfer gradient-echo (ihMT-GRE) imaging to radiofrequency (RF) transmit field ( B 1 + ) inhomogeneities at 3 T.

METHODS

The ihMT-GRE sequence was optimized by varying the concentration of the RF saturation energy over time, obtained by increasing the saturation pulse power while extending the sequence repetition time (TR). Different protocols were tested using numerical simulations and human in vivo experiments in the brain white matter (WM) of healthy subjects at 3 T. The sensitivity of the ihMT ratio (ihMTR) to B 1 + variations was investigated by comparing measurements obtained at nominal transmitter adjustments and following a 20% global B 1 + drop. The resulting relative variations (δ_ihMTR ) were evaluated voxelwise as a function of the local B 1 + distribution. The reproducibility of the protocol providing minimal B 1 + bias was assessed in a test-retest experiment.

RESULTS

In line with simulations, ihMT-GRE experiments conducted at high concentration of the RF energy over time demonstrated strong reduction of the B 1 + inhomogeneity effects in the human WM. Under the optimal conditions of 350-ms TR and 3-µT root mean square (RMS) saturation power, 73% of all WM voxels presented δ_ihMTR below 10%. Reproducibility analysis yielded a close-to-zero systematic bias (ΔihMTR = -0.081%) and a high correlation (ρ² = 0.977) between test and retest experiments.

CONCLUSION

Concentrating RF saturation energy in ihMT-GRE sequences mitigates the sensitivity of the ihMTR to B 1 + variations and allows for clinical-ready ihMT imaging at 3 T. This feature is of particular interest for high and ultra-high field applications.

Complex B1 + mapping with Carr-Purcell spin echoes and its application to electrical properties tomography

Purpose

To present a new complex‐valued B₁⁺ mapping method for electrical properties tomography using Carr‐Purcell spin echoes.

Methods

A Carr‐Purcell (CP) echo train generates pronounced flip‐angle dependent oscillations that can be used to estimate the magnitude of B₁⁺. To this end, a dictionary is used that takes into account the slice profile as well as T₂ relaxation along the echo train. For validation, the retrieved B₁⁺ map is compared with the actual flip angle imaging (AFI) method in a phantom (79 ε₀, 0.34 S/m). Moreover, the phase of the first echo reflects the transceive phase. Overall, the CP echo train yields an estimate of the complex‐valued B₁⁺, allowing electrical properties tomography with both permittivity and conductivity. The presented method is evaluated in phantom scans as well as for in vivo brain at 3 T.

Results

In the phantom, the obtained magnitude B₁⁺ maps retrieved from the CP echo train and the AFI method show excellent agreement, and both the reconstructed estimated permittivity (79 ± 3) ε₀ and conductivity (0.35 ± 0.04) S/m values are in accordance with expectations. In the brain, the obtained electrical properties are also close to expectations. In addition to the retrieved complex B₁⁺ information, the decay of the CP echo trains also yields an estimate for T₂.

Conclusion

The CP sequence can be used to simultaneously provide both B₁⁺ magnitude and phase estimations, and therefore allows for full reconstruction of the electrical properties.

A unifying view on extended phase graphs and Bloch simulations for quantitative MRI

Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.

Explicit phenomenological solutions for magnetization exposed to an arbitrary NMR diffusion steady state pulse sequence

Explicit phenomenological solutions to recurrence relations for the bulk transverse and longitudinal magnetization found using the Torrey-Bloch equations with relaxation effects are used to investigate nuclear magnetic resonance (NMR) diffusion measurements. Of particular interest are steady state NMR (self-)diffusion measurements that reduce experimental time that can extend the techniques to quickly reacting systems. The solutions for bulk transverse and longitudinal magnetization presented here are used to investigate the average behavior of the transverse and longitudinal magnetization in forming a steady state and are used to derive new expressions for the steady state longitudinal magnetization. These solutions can be applied to a noninteracting spin 1/2 ensemble undergoing free diffusion exposed to an arbitrary NMR pulse sequence containing arbitrary magnetic field gradient waveforms. The closed algebraic form method presented here has an advantage over iterative procedures for calculating transverse and longitudinal magnetization for the analysis and development of steady state pulse sequences. Previous theoretical results for steady state diffusion measurements are also reproduced. The Mathematica code for these solutions is provided in the supplementary material.

Generalized low-rank nonrigid motion-corrected reconstruction for MR fingerprinting

PURPOSE

Develop a novel low-rank motion-corrected (LRMC) reconstruction for nonrigid motion-corrected MR fingerprinting (MRF).

METHODS

Generalized motion-corrected (MC) reconstructions have been developed for steady-state imaging. Here we extend this framework to enable nonrigid MC for transient imaging applications with varying contrast, such as MRF. This is achieved by integrating low-rank dictionary-based compression into the generalized MC model to reconstruct MC singular images, reducing motion artifacts in the resulting parametric maps. The proposed LRMC reconstruction was applied for cardiac motion correction in 2D myocardial MRF (T₁ and T₂ ) with extended cardiac acquisition window (~450 ms) and for respiratory MC in free-breathing 3D myocardial and 3D liver MRF. Experiments were performed in phantom and 22 healthy subjects. The proposed approach was compared with reference spin echo (phantom) and with 2D electrocardiogram-triggered/breath-hold MOLLI and T₂ gradient-and-spin echo conventional maps (in vivo 2D and 3D myocardial MRF).

RESULTS

Phantom results were in general agreement with reference spin-echo measurements, presenting relative errors of approximately 5.4% and 5.5% for T₁ and short T₂ (<100 ms), respectively. The proposed LRMC MRF reduced residual blurring artifacts with respect to no MC for cardiac or respiratory motion in all cases (2D and 3D myocardial, 3D abdominal). In 2D myocardial MRF, left-ventricle T₁ values were 1150 ± 41 ms for LRMC MRF and 1010 ± 56 ms for MOLLI; T₂ values were 43.8 ± 2.3 ms for LRMC MRF and 49.5 ± 4.5 ms for T₂ gradient and spin echo. Corresponding measurements for 3D myocardial MRF were 1085 ± 30 ms and 1062 ± 29 ms for T₁ , and 43.5 ± 1.9 ms and 51.7 ± 1.7 ms for T₂ . For 3D liver, LRMC MRF measured liver T₁ at 565 ± 44 ms and liver T₂ at 35.4 ± 2.4 ms.

CONCLUSION

The proposed LRMC reconstruction enabled generalized (nonrigid) MC for 2D and 3D MRF, both for cardiac and respiratory motion. The proposed approach reduced motion artifacts in the MRF maps with respect to no motion compensation and achieved good agreement with reference measurements.

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

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

Quantitative T1 and T2 mapping by magnetic resonance fingerprinting (MRF) of the placenta before and after maternal hyperoxia

INTRODUCTION

MR relaxometry has been used to assess placental exchange function, but methods to date are not sufficiently fast to be robust to placental motion. Magnetic resonance fingerprinting (MRF) permits rapid, voxel-wise, intrinsically co-registered T₁ and T₂ mapping. After characterizing measurement error, we scanned pregnant women during air and oxygen breathing to demonstrate MRF's ability to detect placental oxygenation changes.

METHODS

The accuracy of FISP-based, sliding-window reconstructed MRF was tested on phantoms. MRF scans in 9-s breath holds were acquired at 3T in 31 pregnant women during air and oxygen breathing. A mixed effects model was used to test for changes in placenta relaxation times between physiological states, to assess the dependency on gestational age (GA), and the impact of placental motion.

RESULTS

MRF estimates of known phantom relaxation times resulted in mean absolute errors for T₁ of 92 ms (4.8%), but T₂ was less accurate at 16 ms (13.6%). During normoxia, placental T₁ = 1825 ± 141 ms (avg ± standard deviation) and T₂ = 60 ± 16 ms (gestational age range 24.3-36.7, median 32.6 weeks). In the statistical model, placental T₂ rose and T₁ remained contant after hyperoxia, and no GA dependency was observed for T₁ or T₂.

DISCUSSION

Well-characterized, motion-robust MRF was used to acquire T₁ and T₂ maps of the placenta. Changes with hyperoxia are consistent with a net increase in oxygen saturation. Toward the goal of whole-placenta quantitative oxygenation imaging over time, we aim to implement 3D MRF with integrated motion correction to improve T₂ accuracy.

Accuracy and repeatability of QRAPMASTER and MRF-vFA

Our purpose is to evaluate bias and repeatability of the quantitative MRI sequences QRAPMASTER, based on steady-state imaging, and variable Flip Angle MRF (MRF-VFA), based on the transient response. Both techniques are assessed with a standardized phantom and five volunteers on 1.5 T and 3 T clinical scanners. All scans were repeated eight times in consecutive weeks. In the phantom, the mean bias±95% confidence interval for T1 values with QRAPMASTER was 10 ± 10% on 1.5 T and 4 ± 13% on 3.0 T. The mean bias for T1 values with MRF-vFA was 21 ± 17% on 1.5 T and 9 ± 9% on 3.0 T. For T2 values the mean bias with QRAPMASTER was 12 ± 3% on 1.5 T and 23 ± 1% on 3.0 T. For T2 values the mean bias with MRF-vFA was 17 ± 1% on 1.5 T and 19 ± 2% on 3.0 T. QRAPMASTER estimated lower T1 and T2 values than MRF-vFA. Repeatability was good with low coefficients of variation (CoV). Mean CoV ± 95% confidence interval for T1 were 3.2 ± 0.4% on 1.5 T and 4.5 ± 0.8% on 3.0 T with QRAPMASTER and 2.7% ± 0.2% on 1.5 T and 2.5 ± 0.2% on 3.0 T with MRF-vFA. For T2 were 3.3 ± 1.9% on 1.5 T and 3.2 ± 0.6% on 3.0 T with QRAPMASTER and 2.0 ± 0.4% on 1.5 T and 5.7 ± 1.0% on 3.0 T with MRF-vFA. The in-vivo T1 and T2 are in the range of values previously reported by other authors. The in-vivo mean CoV ± 95% confidence interval in gray matter were for T1 1.7 ± 0.2% using QRAPMASTER and 0.7 ± 0.5% using MRF-vFA and for T2 were 0.9 ± 0.4% using QRAPMASTER and 2.4 ± 0.5% using MRF-vFA. In white matter were for T1 0.9 ± 0.3% using QRAPMASTER and 1.3 ± 1.1% using MRF-vFA and for T2 were 0.7 ± 0.4% using QRAPMASTER and 2.4 ± 0.4% using MRF-vFA. A GLM analysis showed that the variations in T1 and T2 mainly depend on the field strength and the subject, but not on the follow-up repetition in different days. This confirms the high repeatability of QRAPMASTER and MRF-vFA. In summary, QRAPMASTER and MRF-vFA on both systems were highly repeatable with moderate accuracy, providing results comparable to standard references. While repeatability was similar for both methods, QRAPMASTER was more accurate. QRAPMASTER is a tested commercial product but MRF-vFA is 4.77 times faster, which would ease the inclusion of quantitative relaxometry.

Toward magnetic resonance fingerprinting for low-field MR-guided radiation therapy

PURPOSE

The acquisition of multiparametric quantitative magnetic resonance imaging (qMRI) is becoming increasingly important for functional characterization of cancer prior to- and throughout the course of radiation therapy. The feasibility of a qMRI method known as magnetic resonance fingerprinting (MRF) for rapid T₁ and T₂ mapping was assessed on a low-field MR-linac system.

METHODS

A three-dimensional MRF sequence was implemented on a 0.35T MR-guided radiotherapy system. MRF-derived measurements of T₁ and T₂ were compared to those obtained with gold standard single spin echo methods, and the impacts of the radiofrequency field homogeneity and scan times ranging between 6 and 48 min were analyzed by acquiring between 1 and 8 spokes per time point in a standard quantitative system phantom. The short-term repeatability of MRF was assessed over three measurements taken over a 10-h period. To evaluate transferability, MRF measurements were acquired on two additional MR-guided radiotherapy systems. Preliminary human volunteer studies were performed.

RESULTS

The phantom benchmarking studies showed that MRF is capable of mapping T₁ and T₂ values within 8% and 10% of gold standard measures, respectively, at 0.35T. The coefficient of variation of T₁ and T₂ estimates over three repeated scans was < 5% over a broad range of relaxation times. The T₁ and T₂ times derived using a single-spoke MRF acquisition across three scanners were near unity and mean percent errors in T₁ and T₂ estimates using the same phantom were < 3%. The mean percent differences in T₁ and T₂ as a result of truncating the scan time to 6 min over the large range of relaxation times in the system phantom were 0.65% and 4.05%, respectively.

CONCLUSIONS

The technical feasibility and accuracy of MRF on a low-field MR-guided radiation therapy device has been demonstrated. MRF can be used to measure accurate T₁ and T₂ maps in three dimensions from a brief 6-min scan, offering strong potential for efficient and reproducible qMRI for future clinical trials in functional plan adaptation and tumor/normal tissue response assessment.

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.

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.

Improving B 1 + parametric estimation in the brain from multispin-echo sequences using a fusion bootstrap moves solver

PURPOSE

To simultaneously estimate the B 1 + field (along with the T₂ ) in the brain with multispin-echo (MSE) sequences and dictionary matching.

METHODS

T₂ mapping provides clinically relevant information such as in the assessment of brain degenerative diseases. It is commonly obtained with MSE sequences, and accuracy can be further improved by matching the MSE signal to a precomputed dictionary of echo-modulation curves. For additional T₁ quantification, transmit B 1 + field knowledge is also required. Preliminary work has shown that although simultaneous brain B 1 + estimation along with T₂ is possible, it presents a bimodal distribution with the main peak coinciding with the true value. By taking advantage of this, the B 1 + maps are expected to be spatially smooth by applying an iterative method that takes into account each pixel neighborhood known as the fusion bootstrap moves solver (FBMS). The effect of the FBMS on B 1 + accuracy and piecewise smoothness is investigated and different spatial regularization levels are compared. Total variation regularization was used for both B 1 + and T₂ simultaneous estimation because of its simplicity as an initial proof-of-concept; future work could explore non edge-preserving regularization independently for B 1 + .

RESULTS

Improvements in B 1 + accuracy (up to 45.37% and 16.81% B 1 + error decrease) and recovery of spatially homogeneous maps are shown in simulations and in vivo 3.0T brain data, respectively.

CONCLUSION

Accurate B 1 + estimated values can be obtained from widely available MSE sequences while jointly estimating T₂ maps with the use of echo-modulation curve matching and FBMS at no further cost.

Fast and accurate modeling of transient-state, gradient-spoiled sequences by recurrent neural networks

Fast and accurate modeling of MR signal responses are typically required for various quantitative MRI applications, such as MR fingerprinting. This work uses a new extended phase graph (EPG)-Bloch model for accurate simulation of transient-state, gradient-spoiled MR sequences, and proposes a recurrent neural network (RNN) as a fast surrogate of the EPG-Bloch model for computing large-scale MR signals and derivatives. The computational efficiency of the RNN model is demonstrated by comparisons with other existing models, showing one to three orders of acceleration compared with the latest GPU-accelerated, open-source EPG package. By using numerical and in vivo brain data, two used cases, namely, MRF dictionary generation and optimal experimental design, are also provided. Results show that the RNN surrogate model can be efficiently used for computing large-scale dictionaries of transient-state signals and derivatives within tens of seconds, resulting in several orders of magnitude acceleration with respect to state-of-the-art implementations. The practical application of transient-state quantitative techniques can therefore be substantially facilitated.

Synergistic multi-contrast cardiac magnetic resonance image reconstruction

Cardiac magnetic resonance imaging (CMR) is an important tool for the non-invasive diagnosis of a variety of cardiovascular diseases. Parametric mapping with multi-contrast CMR is able to quantify tissue alterations in myocardial disease and promises to improve patient care. However, magnetic resonance imaging is an inherently slow imaging modality, resulting in long acquisition times for parametric mapping which acquires a series of cardiac images with different contrasts for signal fitting or dictionary matching. Furthermore, extra efforts to deal with respiratory and cardiac motion by triggering and gating further increase the scan time. Several techniques have been developed to speed up CMR acquisitions, which usually acquire less data than that required by the Nyquist-Shannon sampling theorem, followed by regularized reconstruction to mitigate undersampling artefacts. Recent advances in CMR parametric mapping speed up CMR by synergistically exploiting spatial-temporal and contrast redundancies. In this article, we will review the recent developments in multi-contrast CMR image reconstruction for parametric mapping with special focus on low-rank and model-based reconstructions. Deep learning-based multi-contrast reconstruction has recently been proposed in other magnetic resonance applications. These developments will be covered to introduce the general methodology. Current technical limitations and potential future directions are discussed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

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.

Optimizing trajectory ordering for fast radial ultra-short TE (UTE) acquisitions

PURPOSE

Additional spoiler gradients are required in 3D UTE sequences with random view ordering to suppress magnetization refocusing. By leveraging the encoding gradient induced spoiling effect, the spoiler gradients could potentially be reduced or removed to shorten the TR and increase encoding efficiency. An analysis framework is built that models the gradient spoiling effects and a new ordering scheme is proposed for fast 3D UTE acquisition.

THEORY AND METHODS

UTE signal evolution and spatial encoding gradient induced spoiling effect are derived from the Bloch equations. And the concept is validated in 2D radial UTE simulation. Then an optimized ordering scheme, named reordered 2D golden angle (r2DGA) scheme, for 3D UTE acquisition is proposed. The r2DGA scheme is compared to the sequential and 3D golden angle schemes in both phantom and volunteer studies.

RESULTS

The proposed r2DGA ordering scheme was applied to two applications, single breath-holding and free breathing 3D lung MRI. With r2DGA ordering scheme, breath-holding lung MRI scan increased 60% scan efficiency by removing the spoiler gradients and the free breathing scan reduced 20% scan time compared to the 3D golden angle scheme by reducing the spoiler gradients.

CONCLUSIONS

The proposed r2DGA ordering scheme UTE acquisition reduces the need of spoiler gradients and increases the encoding efficiency, and shows improvements in both breath-holding and free breathing lung MRI applications.

Preserved thenar muscles in non-ambulant Duchenne muscular dystrophy patients

BACKGROUND

Clinical trials in Duchenne muscular dystrophy (DMD) focus primarily on ambulant patients. Results cannot be extrapolated to later disease stages due to a decline in targeted muscle tissue. In non-ambulant DMD patients, hand function is relatively preserved and crucial for daily-life activities. We used quantitative MRI (qMRI) to establish whether the thenar muscles could be valuable to monitor treatment effects in non-ambulant DMD patients.

METHODS

Seventeen non-ambulant DMD patients (range 10.2-24.1 years) and 13 healthy controls (range 9.5-25.4 years) underwent qMRI of the right hand at 3 T at baseline. Thenar fat fraction (FF), total volume (TV), and contractile volume (CV) were determined using 4-point Dixon, and T2_water was determined using multiecho spin-echo. Clinical assessments at baseline (n = 17) and 12 months (n = 13) included pinch strength (kg), performance of the upper limb (PUL) 2.0, DMD upper limb patient reported outcome measure (PROM), and playing a video game for 10 min using a game controller. Group differences and correlations were assessed with non-parametric tests.

RESULTS

Total volume was lower in patients compared with healthy controls (6.9 cm³ , 5.3-9.0 cm³ vs. 13.0 cm³ , 7.6-15.8 cm³ , P = 0.010). CV was also lower in patients (6.3 cm³ , 4.6-8.3 cm³ vs. 11.9 cm³ , 6.9-14.6 cm³ , P = 0.010). FF was slightly elevated (9.7%, 7.3-11.4% vs. 7.7%, 6.6-8.4%, P = 0.043), while T2_water was higher (31.5 ms, 30.0-32.6 ms vs. 28.1 ms, 27.8-29.4 ms, P < 0.001). Pinch strength and PUL decreased over 12 months (2.857 kg, 2.137-4.010 to 2.243 kg, 1.930-3.339 kg, and 29 points, 20-36 to 23 points, 17-30, both P < 0.001), while PROM did not (49 points, 36-57 to 44 points, 30-54, P = 0.041). All patients were able to play for 10 min at baseline or follow-up, but some did not comply with the study procedures regarding this endpoint. Pinch strength correlated with TV and CV in patients (rho = 0.72 and rho = 0.68) and controls (both rho = 0.89). PUL correlated with TV, CV, and T2_water (rho = 0.57, rho = 0.51, and rho = -0.59).

CONCLUSIONS

Low thenar FF, increased T2_water , correlation of muscle size with strength and function, and the decrease in strength and function over 1 year indicate that the thenar muscles are a valuable and quantifiable target for therapy in later stages of DMD. Further studies are needed to relate these data to the loss of a clinically meaningful milestone.

Non-contrast MRI of synovitis in the knee using quantitative DESS

OBJECTIVES

To determine whether synovitis graded by radiologists using hybrid quantitative double-echo in steady-state (qDESS) images can be utilized as a non-contrast approach to assess synovitis in the knee, compared against the reference standard of contrast-enhanced MRI (CE-MRI).

METHODS

Twenty-two knees (11 subjects) with moderate to severe osteoarthritis (OA) were scanned using CE-MRI, qDESS with a high diffusion weighting (qDESS_High), and qDESS with a low diffusion weighting (qDESS_Low). Four radiologists graded the overall impression of synovitis, their diagnostic confidence, and regional grading of synovitis severity at four sites (suprapatellar pouch, intercondylar notch, and medial and lateral peripatellar recesses) in the knee using a 4-point scale. Agreement between CE-MRI and qDESS, inter-rater agreement, and intra-rater agreement were assessed using a linearly weighted Gwet's AC2.

RESULTS

Good agreement was seen between CE-MRI and both qDESS_Low (AC2 = 0.74) and qDESS_High (AC2 = 0.66) for the overall impression of synovitis, but both qDESS sequences tended to underestimate the severity of synovitis compared to CE-MRI. Good inter-rater agreement was seen for both qDESS sequences (AC2 = 0.74 for qDESS_Low, AC2 = 0.64 for qDESS_High), and good intra-rater agreement was seen for both sequences as well (qDESS_Low AC2 = 0.78, qDESS_High AC2 = 0.80). Diagnostic confidence was moderate to high for qDESS_Low (mean = 2.36) and slightly less than moderate for qDESS_High (mean = 1.86), compared to mostly high confidence for CE-MRI (mean = 2.73).

CONCLUSIONS

qDESS shows potential as an alternative MRI technique for assessing the severity of synovitis without the use of a gadolinium-based contrast agent.

KEY POINTS

The use of the quantitative double-echo in steady-state (qDESS) sequence for synovitis assessment does not require the use of a gadolinium-based contrast agent. Preliminary results found that low diffusion-weighted qDESS (qDESS_Low) shows good agreement to contrast-enhanced MRI for characterization of the severity of synovitis, with a relative bias towards underestimation of severity. Preliminary results also found that qDESS_Low shows good inter- and intra-rater agreement for the depiction of synovitis, particularly for readers experienced with the sequence.

Ensuring both velocity and spatial responses robust to B0/B1+ field inhomogeneities for velocity-selective arterial spin labeling through dynamic phase-cycling

PURPOSE

To evaluate both velocity and spatial responses of velocity-selective arterial spin labeling (VS-ASL), using velocity-insensitive and velocity-compensated waveforms for control modules, as well as a novel dynamic phase-cycling approach, at different B0 / B 1 + field inhomogeneities.

METHODS

In the presence of imperfect refocusing, the mechanism of phase-cycling the refocusing pulses through four dynamics was first theoretically analyzed with the conventional velocity-selective saturation (VSS) pulse train. Numerical simulations were then deployed to compare the performance of the Fourier-transform based velocity-selective inversion (FT-VSI) with these three different schemes in terms of both velocity and spatial responses under various B0 / B 1 + conditions. Phantom and human brain scans were performed to evaluate the three methods at B 1 + scales of 0.8, 1.0, and 1.2.

RESULTS

The simulations of FT-VSI showed that, under nonuniform B0 / B 1 + conditions, the scheme with velocity-insensitive control was susceptible to DC bias of the static spins as systematic error, while the scheme with velocity-compensated control had deteriorated velocity-selective labeling profiles and, thus, reduced labeling efficiency. Through numerical simulation, phantom scans, and brain perfusion measurements, the dynamic phase-cycling method demonstrated considerable improvements over these issues.

CONCLUSION

The proposed dynamic phase-cycling approach was demonstrated for the velocity-selective label and control modules with both velocity and spatial responses robust to a wide range of B0 and B 1 + field inhomogeneities.

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

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

Accelerated 3D whole-brain T1, T2, and proton density mapping: feasibility for clinical glioma MR imaging

Purpose

Advanced MRI-based biomarkers offer comprehensive and quantitative information for the evaluation and characterization of brain tumors. In this study, we report initial clinical experience in routine glioma imaging with a novel, fully 3D multiparametric quantitative transient-state imaging (QTI) method for tissue characterization based on T1 and T2 values.

Methods

To demonstrate the viability of the proposed 3D QTI technique, nine glioma patients (grade II–IV), with a variety of disease states and treatment histories, were included in this study. First, we investigated the feasibility of 3D QTI (6:25 min scan time) for its use in clinical routine imaging, focusing on image reconstruction, parameter estimation, and contrast-weighted image synthesis. Second, for an initial assessment of 3D QTI-based quantitative MR biomarkers, we performed a ROI-based analysis to characterize T1 and T2 components in tumor and peritumoral tissue.

Results

The 3D acquisition combined with a compressed sensing reconstruction and neural network-based parameter inference produced parametric maps with high isotropic resolution (1.125 × 1.125 × 1.125 mm³ voxel size) and whole-brain coverage (22.5 × 22.5 × 22.5 cm³ FOV), enabling the synthesis of clinically relevant T1-weighted, T2-weighted, and FLAIR contrasts without any extra scan time. Our study revealed increased T1 and T2 values in tumor and peritumoral regions compared to contralateral white matter, good agreement with healthy volunteer data, and high inter-subject consistency.

Conclusion

3D QTI demonstrated comprehensive tissue assessment of tumor substructures captured in T1 and T2 parameters. Aiming for fast acquisition of quantitative MR biomarkers, 3D QTI has potential to improve disease characterization in brain tumor patients under tight clinical time-constraints.

Combining prospective and retrospective motion correction based on a model for fast continuous motion

PURPOSE

Prospective motion correction (PMC) and retrospective motion correction (RMC) have different advantages and limitations. The present work aims to combine the advantages of both for rigid body motion, aiming at correcting for faster motions than was previously achievable. Additionally, it provides insights into the effects of motion on pulse sequences and MR signals with a goal of further improving motion correction in the future.

METHODS

The effective encoding trajectory and a global phase offset in a moving object are calculated based on complete gradient waveforms of an arbitrary sequence and a continuous motion model. These data are used to feed the forward signal model, which is then used in iterative image reconstruction to suppress the artifacts still present after the PMC.

RESULTS

Verification experiments with a rotation phantom and in vivo were performed. Predictions of simulated motion artifacts for PMC based on sequence waveforms are very accurate. The performance at combined PMC+RMC is limited by Nyquist violations in the sampled k-space and can be compensated by oversampling.

CONCLUSION

The combined correction results in better images than pure PMC in the presence of fast motion. The predictions of artifacts are very accurate, allowing for comparing sequences or protocols in simulations. The observed artifacts due to Nyquist violations are expected to be corrected by utilizing parallel imaging.

Pilot Study on Quantitative Cervical Cord and Muscular MRI in Spinal Muscular Atrophy: Promising Biomarkers of Disease Evolution and Treatment?

Introduction: Nusinersen is a recent promising therapy approved for the treatment of spinal muscular atrophy (SMA), a rare disease characterized by the degeneration of alpha motor neurons (αMN) in the spinal cord (SC) leading to progressive muscle atrophy and dysfunction. Muscle and cervical SC quantitative magnetic resonance imaging (qMRI) has never been used to monitor drug treatment in SMA. The aim of this pilot study is to investigate whether qMRI can provide useful biomarkers for monitoring treatment efficacy in SMA. Methods: Three adult SMA 3a patients under treatment with nusinersen underwent longitudinal clinical and qMRI examinations every 4 months from baseline to 21-month follow-up. The qMRI protocol aimed to quantify thigh muscle fat fraction (FF) and water-T2 (w-T2) and to characterize SC volumes and microstructure. Eleven healthy controls underwent the same SC protocol (single time point). We evaluated clinical and imaging outcomes of SMA patients longitudinally and compared SC data between groups transversally. Results: Patient motor function was stable, with only Patient 2 showing moderate improvements. Average muscle FF was already high at baseline (50%) and progressed over time (57%). w-T2 was also slightly higher than previously published data at baseline and slightly decreased over time. Cross-sectional area of the whole SC, gray matter (GM), and ventral horns (VHs) of Patients 1 and 3 were reduced compared to controls and remained stable over time, while GM and VHs areas of Patient 2 slightly increased. We found altered diffusion and magnetization transfer parameters in SC structures of SMA patients compared to controls, thus suggesting changes in tissue microstructure and myelin content. Conclusion: In this pilot study, we found a progression of FF in thigh muscles of SMA 3a patients during nusinersen therapy and a concurrent slight reduction of w-T2 over time. The SC qMRI analysis confirmed previous imaging and histopathological studies suggesting degeneration of αMN of the VHs, resulting in GM atrophy and demyelination. Our longitudinal data suggest that qMRI could represent a feasible technique for capturing microstructural changes induced by SMA in vivo and a candidate methodology for monitoring the effects of treatment, once replicated on a larger cohort.

Simultaneous quantification of SPIO and gadolinium contrast agents using MR fingerprinting

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Strategies to improve intratrain prospective motion correction for turbo spin-echo sequences with constant flip angles

PURPOSE

To investigate the effects of prospective motion correction on turbo spin echo sequences and optimize motion correction approaches, mitigating signal dropout artifacts caused by the imperfections of motion tracking data.

METHODS

Signal dropout artifacts caused by undesired phase deviations introduced by tracking errors are analyzed theoretically. To reduce the adverse effect of such deviations, two approaches are proposed: (1) freezing the correction for example, for even-numbered or higher number of echoes and (2) shifting the correction event prior to the left crusher gradient preceding the refocusing pulse. A comprehensive analysis is presented, including both signal simulations and experimental verifications in phantoms and in vivo. Performance of the proposed approach is validated in two healthy volunteers imaged under two types of motion conditions simulating inadvertent fast motions associated with discomfort and continuous large motions.

RESULTS

The results show that the proposed optimization is able to efficiently correct for the motion artifacts and at the same time avoid signal dropout artifacts. Specifically, performing correction every 4^th echo prior to the left crusher gradient was shown to improve image quality.

CONCLUSION

An optimization approach is proposed to exploit the potential of external tracking for intra-echo-train motion artifact correction for turbo spin echo sequences.

Bayesian uncertainty quantification for magnetic resonance fingerprinting

Magnetic Resonance Fingerprinting (MRF) is a promising technique for fast quantitative imaging of human tissue. In general, MRF is based on a sequence of highly undersampled MR images which are analyzed with a pre-computed dictionary. MRF provides valuable diagnostic parameters such as the $T_1$ and $T_2$ MR relaxation times. However, uncertainty characterization of dictionary-based MRF estimates for $T_1$ and $T_2$ has not been achieved so far, which makes it challenging to assess if observed differences in these estimates are significant and may indicate pathological changes of the underlying tissue. We propose a Bayesian approach for the uncertainty quantification of dictionary-based MRF which leads to probability distributions for $T_1$ and $T_2$ in every voxel. The distributions can be used to make probability statements about the relaxation times, and to assign uncertainties to their dictionary-based MRF estimates. All uncertainty calculations are based on the pre-computed dictionary and the observed sequence of undersampled MR images, and they can be calculated in short time. The approach is explored by analyzing MRF measurements of a phantom consisting of several tubes across which MR relaxation times are constant. The proposed uncertainty quantification is quantitatively consistent with the observed within-tube variability of estimated relaxation times. Furthermore, calculated uncertainties are shown to characterize well observed differences between the MRF estimates and the results obtained from high-accurate reference measurements. These findings indicate that a reliable uncertainty quantification is achieved. We also present results for simulated MRF data and an uncertainty quantification for an in vivo MRF measurement. MATLAB$^{\scriptsize \text{\textregistered}}$ source code implementing the proposed approach is made available.

High-Spatial-Resolution 3D Whole-Heart MRI T2 Mapping for Assessment of Myocarditis

Background Clinical guidelines recommend the use of established T2 mapping sequences to detect and quantify myocarditis and edema, but T2 mapping is performed in two dimensions with limited coverage and repetitive breath holds. Purpose To assess the reproducibility of an accelerated free-breathing three-dimensional (3D) whole-heart T2 MRI mapping sequence in phantoms and participants without a history of cardiac disease and to investigate its clinical performance in participants with suspected myocarditis. Materials and Methods Eight participants (three women, mean age, 31 years ± 4 [standard deviation]; cohort 1) without a history of cardiac disease and 25 participants (nine women, mean age, 45 years ± 17; cohort 2) with clinically suspected myocarditis underwent accelerated free-breathing 3D whole-heart T2 mapping with 100% respiratory scanning efficiency at 1.5 T. The participants were enrolled from November 2018 to August 2020. Three repeated scans were performed on 2 separate days in cohort 1. Segmental variations in T2 relaxation times of the left ventricular myocardium were assessed, and intrasession and intersession reproducibility were measured. In cohort 2, segmental myocardial T2 values, detection of focal inflammation, and map quality were compared with those obtained from clinical breath-hold two-dimensional (2D) T2 mapping. Statistical differences were assessed using the nonparametric Mann-Whitney and Kruskal-Wallis tests, whereas the paired Wilcoxon signed-rank test was used to assess subjective scores. Results Whole-heart T2 maps were acquired in a mean time of 6 minutes 53 seconds ± 1 minute 5 seconds at 1.5 mm3 resolution. Breath-hold 2D and free-breathing 3D T2 mapping had similar intrasession (mean T2 change of 3.2% and 2.3% for 2D and 3D, respectively) and intersession (4.8% and 4.9%, respectively) reproducibility. The two T2 mapping sequences showed similar map quality (P = .23, cohort 2). Abnormal myocardial segments were identified with confidence (score 3) in 14 of 25 participants (56%) with 3D T2 mapping and only in 10 of 25 participants (40%) with 2D T2 mapping. Conclusion High-spatial-resolution three-dimensional (3D) whole-heart T2 mapping shows high intrasession and intersession reproducibility and helps provide T2 myocardial characterization in agreement with clinical two-dimensional reference, while enabling 3D assessment of focal disease with higher confidence. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Friedrich in this issue.

Fast Open-Source Toolkit for Water T2 Mapping in the Presence of Fat From Multi-Echo Spin-Echo Acquisitions for Muscle MRI

Imaging has become a valuable tool in the assessment of neuromuscular diseases, and, specifically, quantitative MR imaging provides robust biomarkers for the monitoring of disease progression. Quantitative evaluation of fat infiltration and quantification of the T2 values of the muscular tissue's water component (wT2) are two of the most essential indicators currently used. As each voxel of the image can contain both water and fat, a two-component model for the estimation of wT2 must be used. In this work, we present a fast method for reconstructing wT2 maps obtained from conventional multi-echo spin-echo (MESE) acquisitions and released as Free Open Source Software. The proposed software is capable of fast reconstruction thanks to extended phase graphs (EPG) simulations and dictionary matching implemented on a general-purpose graphic processing unit. The program can also perform more conventional biexponential least-squares fitting of the data and incorporate information from an external water-fat acquisition to increase the accuracy of the results. The method was applied to the scans of four healthy volunteers and five subjects suffering from facioscapulohumeral muscular dystrophy (FSHD). Conventional multi-slice MESE acquisitions were performed with 17 echoes, and additionally, a 6-echo multi-echo gradient-echo (MEGE) sequence was used for an independent fat fraction calculation. The proposed reconstruction software was applied on the full datasets, and additionally to reduced number of echoes, respectively, to 8, 5, and 3, using EPG and biexponential least-squares fitting, with and without incorporating information from the MEGE acquisition. The incorporation of external fat fraction maps increased the robustness of the fitting with a reduced number of echoes per datasets, whereas with unconstrained fitting, the total of 17 echoes was necessary to retain an independence of wT2 from the level of fat infiltration. In conclusion, the proposed software can successfully be used to calculate wT2 maps from conventional MESE acquisition, allowing the usage of an optimized protocol with similar precision and accuracy as a 17-echo acquisition. As it is freely released to the community, it can be used as a reference for more extensive cohort studies.

Three-dimensional spatially resolved phase graph framework

PURPOSE

An open-source spatially resolved phase graph framework is proposed for simulating arbitrary pulse sequences in the presence of piece-wise constant gradients with arbitrary orientations in three dimensions. It generalizes the extended phase graph algorithm for analysis of nonperiodic sequences while preserving its efficiency, and is able to estimate the signal modulation in the 3D spatial domain.

METHODS

The framework extends the recursive magnetization-evolution algorithm to account for anisotropic diffusion and exploits a novel 3D k-space grid-merging method to balance the computational effort and memory requirements against acceptable simulation errors. A new postsimulation module is proposed to track and visualize the signal evolution both in the k-space and in the image domain, which can be used for simulating image artifacts or finding frequency-response profiles. To illustrate the developed technique, three examples are presented: (1) fast off-resonance calculation for dictionary building in MR fingerprinting, (2) validation of a steady-state sequence with quasi-isotropic diffusion weighting, and (3) investigation of the magnetization evolution in PRESS-based spectroscopic imaging.

RESULTS

The grid-merging algorithm of the proposed framework demonstrates high calculation efficiency exemplified by frequency-response simulation of pseudo steady-state or diffusion-weighted steady-state sequences. It further helps to visualize the signal evolution in PRESS-based sequences.

CONCLUSIONS

The proposed simulation framework has been validated based on several different example applications for analyzing signal evolution in the frequency and spatial domain.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Utilizing shared information between gradient-spoiled and RF-spoiled steady-state MRI signals

This work presents an analytical relationship between gradient-spoiled and RF-spoiled steady-state signals. The two echoes acquired in double-echo in steady-state scans are shown to lie on a line in the signal plane, where the two axes represent the amplitudes of each echo. The location along the line depends on the amount of spoiling and the diffusivity. The line terminates in a point corresponding to an RF-spoiled signal. In addition to the main contribution of demonstrating this signal relationship, we also include the secondary contribution of preliminary results from an example application of the relationship, in the form of a heuristic denoising method when both types of scans are performed. This is investigated in simulations, phantom scans, and in vivo scans. For the signal model, the main topic of this study, simulations confirmed its accuracy and explored its dependency on signal parameters and image noise. For the secondary topic of its preliminary application to reduce noise, simulations demonstrated the denoising method giving a reduction in noise-induced standard deviation of about 30%. The relative effect of the method on the signals is shown to depend on the slope of the described line, which is demonstrated to be zero at the Ernst angle. The phantom scans show a similar effect as the simulations. In vivo scans showed a slightly lower average improvement of about 28%.

Toward Quantification: Microstructure and Magnetic Resonance Fingerprinting

Quantitative magnetic resonance imaging (MRI) is a long-standing challenge. We advocate that the origin of the problem is the simplification applied in commonly used models of the MRI signal relation to the target parameters of biological tissues. Two research fields are briefly reviewed as ways to respond to the challenge of quantitative MRI, both experiencing an exponential growth right now. Microstructure MRI strives to build physiology-based models from cells to signal and, given the signal, back to the cells again. Magnetic resonance fingerprinting aims at efficient simultaneous determination of multiple signal parameters. The synergy of these yet disjoined approaches promises truly quantitative MRI with specific target-oriented diagnostic tools rather than universal imaging methods.

The feasibility of quantitative MRI of extra-ocular muscles in myasthenia gravis and Graves' orbitopathy

Although quantitative MRI can be instrumental in the diagnosis and assessment of disease progression in orbital diseases involving the extra-ocular muscles (EOM), acquisition can be challenging as EOM are small and prone to eye-motion artefacts. We explored the feasibility of assessing fat fractions (FF), muscle volumes and water T2 (T2_water ) of EOM in healthy controls (HC), myasthenia gravis (MG) and Graves' orbitopathy (GO) patients. FF, EOM volumes and T2_water values were determined in 12 HC (aged 22-65 years), 11 MG (aged 28-71 years) and six GO (aged 28-64 years) patients at 7 T using Dixon and multi-echo spin-echo sequences. The EOM were semi-automatically 3D-segmented by two independent observers. MANOVA and t-tests were used to assess differences in FF, T2_water and volume of EOM between groups (P < .05). Bland-Altman limits of agreement (LoA) were used to assess the reproducibility of segmentations and Dixon scans. The scans were well tolerated by all subjects. The bias in FF between the repeated Dixon scans was -0.7% (LoA: ±2.1%) for the different observers; the bias in FF was -0.3% (LoA: ±2.8%) and 0.03 cm³ (LoA: ± 0.36 cm³ ) for volume. Mean FF of EOM in MG (14.1% ± 1.6%) was higher than in HC (10.4% ± 2.5%). Mean muscle volume was higher in both GO (1.2 ± 0.4 cm³ ) and MG (0.8 ± 0.2 cm³ ) compared with HC (0.6 ± 0.2 cm³ ). The average T2_water for all EOM was 24.6 ± 4.0 ms for HC, 24.0 ± 4.7 ms for MG patients and 27.4 ± 4.2 ms for the GO patient. Quantitative MRI at 7 T is feasible for measuring FF and muscle volumes of EOM in HC, MG and GO patients. The measured T2_water was on average comparable with skeletal muscle, although with higher variation between subjects. The increased FF in the EOM in MG patients suggests that EOM involvement in MG is accompanied by fat replacement. The unexpected EOM volume increase in MG may provide novel insights into underlying pathophysiological processes.

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

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

Estimating myelin-water content from anatomical and diffusion images using spatially undersampled myelin-water imaging through machine learning

Myelin is vital for healthy neuronal development, and can therefore provide valuable information regarding neuronal maturation. Anatomical and diffusion weighted images (DWI) possess information related to the myelin content and the current study investigates whether quantitative myelin markers can be extracted from anatomical and DWI using neural networks. Thirteen volunteers (mean age 29y) are included, and for each subject, a residual neural network was trained using spatially undersampled reference myelin-water markers. The network is trained on a voxel-by-voxel basis, resulting in a large amount of training data for each volunteer. The inputs used are the anatomical contrasts (cT1w, cT2w), the standardized T1w/T2w ratio, estimates of the relaxation times (T1, T2) and their ratio (T1/T2), and common DWI metrics (FA, RD, MD, λ₁, λ₂, λ₃). Furthermore, to estimate the added value of the DWI metrics, neural networks were trained using either the combined set (DWI, T1w and T2w) or only the anatomical (T1w and T2w) images. The reconstructed myelin-water maps are in good agreement with the reference myelin-water content in terms of the coefficient of variation (CoV) and the intraclass correlation coefficient (ICC). A 6-fold undersampling using both anatomical and DWI metrics resulted in ICC = 0.68 and CoV = 5.9%. Moreover, using twice the training data (3-fold undersampling) resulted in an ICC that is comparable to the reproducibility of the myelin-water imaging itself (CoV = 5.5% vs. CoV = 6.7% and ICC = 0.74 vs ICC = 0.80). To achieve this, beside the T1w, T2w images, DWI is required. This preliminary study shows the potential of machine learning approaches to extract specific myelin-content from anatomical and diffusion-weighted scans.

Texture analysis and machine learning to predict water T2 and fat fraction from non-quantitative MRI of thigh muscles in Facioscapulohumeral muscular dystrophy

PURPOSE

Quantitative MRI (qMRI) plays a crucial role for assessing disease progression and treatment response in neuromuscular disorders, but the required MRI sequences are not routinely available in every center. The aim of this study was to predict qMRI values of water T2 (wT2) and fat fraction (FF) from conventional MRI, using texture analysis and machine learning.

METHOD

Fourteen patients affected by Facioscapulohumeral muscular dystrophy were imaged at both thighs using conventional and quantitative MR sequences. Muscle FF and wT2 were calculated for each muscle of the thighs. Forty-seven texture features were extracted for each muscle on the images obtained with conventional MRI. Multiple machine learning regressors were trained to predict qMRI values from the texture analysis dataset.

RESULTS

Eight machine learning methods (linear, ridge and lasso regression, tree, random forest (RF), generalized additive model (GAM), k-nearest-neighbor (kNN) and support vector machine (SVM) provided mean absolute errors ranging from 0.110 to 0.133 for FF and 0.068 to 0.115 for wT2. The most accurate methods were RF, SVM and kNN to predict FF, and tree, RF and kNN to predict wT2.

CONCLUSION

This study demonstrates that it is possible to estimate with good accuracy qMRI parameters starting from texture analysis of conventional MRI.

Quantitative magnetization transfer imaging for non-contrast enhanced detection of myocardial fibrosis

PURPOSE

To develop a novel gadolinium-free model-based quantitative magnetization transfer (qMT) technique to assess macromolecular changes associated with myocardial fibrosis.

METHODS

The proposed sequence consists of a two-dimensional breath-held dual shot interleaved acquisition of five MT-weighted (MTw) spoiled gradient echo images, with variable MT flip angles (FAs) and off-resonance frequencies. A two-pool exchange model and dictionary matching were used to quantify the pool size ratio (PSR) and bound pool T2 relaxation ( T 2 B ). The signal model was developed and validated using 25 MTw images on a bovine serum albumin (BSA) phantom and in vivo human thigh muscle. A protocol with five MTw images was optimized for single breath-hold cardiac qMT imaging. The proposed sequence was tested in 10 healthy subjects and 5 patients with myocardial fibrosis and compared to late gadolinium enhancement (LGE).

RESULTS

PSR values in the BSA phantom were within the confidence interval of previously reported values (concentration 10% BSA = 5.9 ± 0.1%, 15% BSA = 9.4 ± 0.2%). PSR and T 2 B in thigh muscle were also in agreement with literature (PSR = 10.9 ± 0.3%, T 2 B = 6.4 ± 0.4 us). In 10 healthy subjects, global left ventricular PSR was 4.30 ± 0.65%. In patients, PSR was reduced in areas associated with LGE (remote: 4.68 ± 0.70% vs. fibrotic: 3.12 ± 0.78 %, n = 5, P < .002).

CONCLUSION

In vivo model-based qMT mapping of the heart was performed for the first time, with promising results for non-contrast enhanced assessment of myocardial fibrosis.

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.

Fast T2 mapping using multi-echo spin-echo MRI: A linear order approach

PURPOSE

Multi-echo spin-echo sequence is commonly used for T2 mapping. The estimated values using conventional exponential fit, however, are hampered by stimulated and indirect echoes leading to overestimation of T2 . Here, we present fast analysis of multi-echo spin-echo (FAMESE) as a novel approach to decrease the complexity of the search space, which leads to accelerated measurement of T2 .

METHODS

We developed FAMESE based on mathematical analysis of the Bloch equations in which the search space dimension decreased to only one. Then, we tested it in both phantom and human brain. Bland-Altman plot was used to assess the agreement between the estimated T2 values from FAMESE and the ones estimated from single-echo spin-echo sequence. The reliability of FAMESE was assessed by intraclass correlation coefficients. In addition, we investigated the noise stability of the method in synthetic and experimental data.

RESULTS

In both phantom and healthy participants, FAMESE provided accelerated and SNR-resistant T2 maps. The FAMESE had a very good agreement with the single-echo spin echo for the whole range of T2 values. The intraclass correlation coefficient values for FAMESE were excellent (ie, 0.9998 and 0.9860 < intraclass correlation coefficient < 0.9942 for the phantom and humans, respectively).

CONCLUSION

Our developed method FAMESE could be considered as a candidate for rapid T2 mapping with a clinically feasible scan time.

MR fingerprinting for rapid simultaneous T1 , T2 , and T1ρ relaxation mapping of the human articular cartilage at 3T

PURPOSE

To implement a novel technique for simultaneous, quantitative multiparametric mapping of the knee articular cartilage.

METHODS

A novel MR fingerprinting pulse sequence is proposed and implemented for simultaneous measurements of proton density, T1 , T2, and T1ρ relaxation times at 3T. The repeatability and reproducibility of the proposed technique were assessed in model phantoms. Institutional review board-approved MR fingerprinting imaging sequence was performed on healthy volunteers and patients with mild knee osteoarthritis. The Wilcoxon test was used to compare healthy controls and patients. The intra- and intersubject repeatability were assessed with coefficient of variation and the RMS coefficient of variation, respectively RESULTS: The Bland-Altman plots demonstrated an average difference of 4.67 ms, -0.09 ms, and 0.05 ms between 2 scans in the same scanner; and 9.68 ms, 0.29 ms, and -0.72 ms between the scans acquired on 2 different scanners for T1 , T2 , and T1ρ , respectively. The in vivo knee study showed excellent repeatability with RMS coefficient of variation less than 3%, 6%, and 5% for T1 , T2 , and T1ρ , respectively. The Wilcoxon test showed a significant difference between control and mild osteoarthritis patients for T1 (P = .04), T2 (P = .01), and T1ρ (P = .02) relaxation time in medial tibial cartilage, as well as for T2 relaxation time (P = .02) in medial femoral cartilage.

CONCLUSION

The proposed MRF sequence is fast and can simultaneously measure the T1 , T2 , T1ρ , and B 1 + maps in a single scan. It is able to discriminate between mild osteoarthritis patients and healthy volunteers.

Rapid Knee MRI Acquisition and Analysis Techniques for Imaging Osteoarthritis

Osteoarthritis (OA) of the knee is a major source of disability that has no known treatment or cure. Morphological and compositional MRI is commonly used for assessing the bone and soft tissues in the knee to enhance the understanding of OA pathophysiology. However, it is challenging to extend these imaging methods and their subsequent analysis techniques to study large population cohorts due to slow and inefficient imaging acquisition and postprocessing tools. This can create a bottleneck in assessing early OA changes and evaluating the responses of novel therapeutics. The purpose of this review article is to highlight recent developments in tools for enhancing the efficiency of knee MRI methods useful to study OA. Advances in efficient MRI data acquisition and reconstruction tools for morphological and compositional imaging, efficient automated image analysis tools, and hardware improvements to further drive efficient imaging are discussed in this review. For each topic, we discuss the current challenges as well as potential future opportunities to alleviate these challenges. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 3.

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.

T1, T2, and Fat Fraction Cardiac MR Fingerprinting: Preliminary Clinical Evaluation

BACKGROUND

Dixon cardiac magnetic resonance fingerprinting (MRF) has been recently introduced to simultaneously provide water T₁ , water T₂ , and fat fraction (FF) maps.

PURPOSE

To assess Dixon cardiac MRF repeatability in healthy subjects and its clinical feasibility in a cohort of patients with cardiovascular disease.

POPULATION

T1MES phantom, water-fat phantom, 11 healthy subjects and 19 patients with suspected cardiovascular disease.

STUDY TYPE

Prospective.

FIELD STRENGTH/SEQUENCE

1.5T, inversion recovery spin echo (IRSE), multiecho spin echo (MESE), modified Look-Locker inversion recovery (MOLLI), T₂ gradient spin echo (T₂ -GRASE), 6-echo gradient rewound echo (GRE), and Dixon cardiac MRF.

ASSESSMENT

Dixon cardiac MRF precision was assessed through repeated scans against conventional MOLLI, T₂ -GRASE, and PDFF in phantom and 11 healthy subjects. Dixon cardiac MRF native T₁ , T₂ , FF, postcontrast T₁ and synthetic extracellular volume (ECV) maps were assessed in 19 patients in comparison to conventional sequences. Measurements in patients were performed in the septum and in late gadolinium enhanced (LGE) areas and assessed using mean value distributions, correlation, and Bland-Altman plots. Image quality and diagnostic confidence were assessed by three experts using 5-point scoring scales.

STATISTICAL TESTS

Paired Wilcoxon rank signed test and paired t-tests were applied. Statistical significance was indicated by *(P < 0.05).

RESULTS

Dixon cardiac MRF showed good overall precision in phantom and in vivo. Septal average repeatability was ~23 msec for T₁ , ~2.2 msec for T₂ , and ~1% for FF. Biases in healthy subjects/patients were measured at +37 msec*/+60 msec* and -8.8 msec*/-8 msec* when compared to MOLLI and T₂ -GRASE, respectively. No statistically significant differences in postcontrast T₁ (P = 0.17) and synthetic ECV (P = 0.19) measurements were observed in patients.

DATA CONCLUSION

Dixon cardiac MRF attained good overall precision in phantom and healthy subjects, while providing coregistered T₁ , T₂ , and fat fraction maps in a single breath-hold scan with similar or better image quality than conventional methods in patients.

LEVEL OF EVIDENCE

2.

TECHNICAL EFFICACY STAGE

2.