Free-Breathing Liver Fat, R₂* and B₀ Field Mapping Using Multi-Echo Radial FLASH and Regularized Model-Based Reconstruction

5D image reconstruction exploiting space-motion-echo sparsity for accelerated free-breathing quantitative liver MRI

Recent advances in 3D non-Cartesian multi-echo gradient-echo (mGRE) imaging and compressed sensing (CS)-based 4D (3D image space + 1D respiratory motion) motion-resolved image reconstruction, which applies temporal total variation to the respiratory motion dimension, have enabled free-breathing liver tissue MR parameter mapping. This technology now allows for robust reconstruction of high-resolution proton density fat fraction (PDFF), R2∗, and quantitative susceptibility mapping (QSM), previously unattainable with conventional Cartesian mGRE imaging. However, long scan times remain a persistent challenge in free-breathing 3D non-Cartesian mGRE imaging. Recognizing that the underlying dimension of the imaging data is essentially 5D (4D + 1D echo signal evolution), we propose a CS-based 5D motion-resolved mGRE image reconstruction method to further accelerate the acquisition. Our approach integrates discrete wavelet transforms along the echo and spatial dimensions into a CS-based reconstruction model and devises a solution algorithm capable of handling such a 5D complex-valued array. Through phantom and in vivo human subject studies, we evaluated the effectiveness of leveraging unexplored correlations by comparing the proposed 5D reconstruction with the 4D reconstruction (i.e., motion-resolved reconstruction with temporal total variation) across a wide range of acceleration factors. The 5D reconstruction produced more reliable and consistent measurements of PDFF, R2∗, and QSM compared to the 4D reconstruction. In conclusion, the proposed 5D motion-resolved image reconstruction demonstrates the feasibility of achieving accelerated, reliable, and free-breathing liver mGRE imaging for the measurement of PDFF, R2∗, and QSM.

Accelerated model-based T1, T2* and proton density mapping using a Bayesian approach with automatic hyperparameter estimation

PURPOSE

To achieve automatic hyperparameter estimation for the model-based recovery of quantitative MR maps from undersampled data, we propose a Bayesian formulation that incorporates the signal model and sparse priors among multiple image contrasts.

THEORY

We introduce a novel approximate message passing framework "AMP-PE" that enables the automatic and simultaneous recovery of hyperparameters and quantitative maps.

METHODS

We employed the variable-flip-angle method to acquire multi-echo measurements using gradient echo sequence. We explored undersampling schemes to incorporate complementary sampling patterns across different flip angles and echo times. We further compared AMP-PE with conventional compressed sensing approaches such as thel 1 $$ {l}_1 $$ -norm minimization, PICS and other model-based approaches such as GraSP, MOBA.

RESULTS

Compared to conventional compressed sensing approaches such as thel 1 $$ {l}_1 $$ -norm minimization and PICS, AMP-PE achieved superior reconstruction performance with lower errors inT 2 ∗ $$ {\mathrm{T}}_2^{\ast } $$ mapping and comparable performance inT 1 $$ {\mathrm{T}}_1 $$ and proton density mappings. When compared to other model-based approaches including GraSP and MOBA, AMP-PE exhibited greater robustness and outperformed GraSP in reconstruction error. AMP-PE offers faster speed than MOBA. AMP-PE performed better than MOBA at higher sampling rates and worse than MOBA at a lower sampling rate. Notably, AMP-PE eliminates the need for hyperparameter tuning, which is a requisite for all the other approaches.

CONCLUSION

AMP-PE offers the benefits of model-based recovery with the additional key advantage of automatic hyperparameter estimation. It works adeptly in situations where ground-truth is difficult to obtain and in clinical environments where it is desirable to automatically adapt hyperparameters to individual protocol, scanner and patient.

Rapid, High-resolution and Distortion-free R 2 * Mapping of Fetal Brain using Multi-echo Radial FLASH and Model-based Reconstruction

Purpose

To develop a rapid, high-resolution and distortion-free quantitativeR 2 * mapping technique for fetal brain at 3 T.

Methods

A 2D multi-echo radial FLASH sequence with blip gradients is adapted for fetal brain data acquisition during maternal free breathing at 3 T. A calibrationless model-based reconstruction with sparsity constraints is developed to jointly estimate water, fat,R 2 * andB 0 field maps directly from the acquired k-space data. Validations have been performed on numerical and NIST phantoms and five fetal subjects ranging from 27 weeks to 36 weeks gestation age.

Results

Both numerical and experimental phantom studies confirm good accuracy and precision of the proposed method. In fetal studies, both the parallel imaging compressed sensing (PICS) technique with a Graph Cut algorithm and the model-based approach proved effective for parameter quantification, with the latter providing enhanced image details. Compared to commonly used multi-echo EPI approaches, the proposed radial technique shows improved spatial resolution (1.1 × 1.1 × 3 mm3 vs. 2-3 × 2-3 × 3 mm3) and reduced distortion. QuantitativeR 2 * results confirm good agreement between the two acquisition strategies. Additionally, high-resolution, distortion-freeR 2 * -weighted images can be synthesized, offering complementary information to HASTE.

Conclusion

This work demonstrates the feasibility of radial acquisition for motion-robust quantitativeR 2 * mapping of the fetal brain. This proposed multi-echo radial FLASH, combined with calibrationless model-based reconstruction, achieves accurate, distortion-free fetal brainR 2 * mapping at a nominal resolution of 1.1 × 1.1 × 3 mm3 within 2 seconds.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Accelerated free-breathing liver fat and R 2 * quantification using multi-echo stack-of-radial MRI with motion-resolved multidimensional regularized reconstruction: Initial retrospective evaluation

PURPOSE

To improve image quality, mitigate quantification biases and variations for free-breathing liver proton density fat fraction (PDFF) andR 2 * $$ {\mathrm{R}}_2^{\ast } $$ quantification accelerated by radial k-space undersampling.

METHODS

A free-breathing multi-echo stack-of-radial MRI method was developed with compressed sensing with multidimensional regularization. It was validated in motion phantoms with reference acquisitions without motion and in 11 subjects (6 patients with nonalcoholic fatty liver disease) with reference breath-hold Cartesian acquisitions. Images, PDFF, andR 2 * $$ {\mathrm{R}}_2^{\ast } $$ maps were reconstructed using different radial view k-space sampling factors and reconstruction settings. Results were compared with reference-standard results using Bland-Altman analysis. Using linear mixed-effects model fitting (p < 0.05 considered significant), mean and SD were evaluated for biases and variations of PDFF andR 2 * $$ {\mathrm{R}}_2^{\ast } $$ , respectively, and coefficient of variation on the first echo image was evaluated as a surrogate for image quality.

RESULTS

Using the empirically determined optimal sampling factor of 0.25 in the accelerated in vivo protocols, mean differences and limits of agreement for the proposed method were [-0.5; -33.6, 32.7] s-1 forR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and [-1.0%; -5.8%, 3.8%] for PDFF, close to those of a previous self-gating method using fully sampled radial views: [-0.1; -27.1, 27.0] s-1 forR 2 * $$ {\mathrm{R}}_2^{\ast } $$ and [-0.4%; -4.5%, 3.7%] for PDFF. The proposed method had significantly lower coefficient of variation than other methods (p < 0.001). Effective acquisition time of 64 s or 59 s was achieved, compared with 171 s or 153 s for two baseline protocols with different radial views corresponding to sampling factor of 1.0.

CONCLUSION

This proposed method may allow accelerated free-breathing liver PDFF andR 2 * $$ {\mathrm{R}}_2^{\ast } $$ mapping with reduced biases and variations.

Free-breathing MRI techniques for fat and R2* quantification in the liver

OBJECTIVE

To review the recent advancements in free-breathing MRI techniques for proton-density fat fraction (PDFF) and R2* quantification in the liver, and discuss the current challenges and future opportunities.

MATERIALS AND METHODS

This work focused on recent developments of different MRI pulse sequences, motion management strategies, and reconstruction approaches that enable free-breathing liver PDFF and R2* quantification.

RESULTS

Different free-breathing liver PDFF and R2* quantification techniques have been evaluated in various cohorts, including healthy volunteers and patients with liver diseases, both in adults and children. Initial results demonstrate promising performance with respect to reference measurements. These techniques have a high potential impact on providing a solution to the clinical need of accurate liver fat and iron quantification in populations with limited breath-holding capacity.

DISCUSSION

As these free-breathing techniques progress toward clinical translation, studies of the linearity, bias, and repeatability of free-breathing PDFF and R2* quantification in a larger cohort are important. Scan acceleration and improved motion management also hold potential for further enhancement.

An Improved Postprocessing Method to Mitigate the Macroscopic Cross-Slice B0 Field Effect on R2* Measurements in the Mouse Brain at 7T

The MR transverse relaxation rate, R2*, has been widely used to detect iron and myelin content in tissue. However, it is also sensitive to macroscopic B0 inhomogeneities. One approach to correct for the B0 effect is to fit gradient-echo signals with the three-parameter model, a sinc function-weighted monoexponential decay. However, such three-parameter models are subject to increased noise sensitivity. To address this issue, this study presents a two-stage fitting procedure based on the three-parameter model to mitigate the B0 effect and reduce the noise sensitivity of R2* measurement in the mouse brain at 7T. MRI scans were performed on eight healthy mice. The gradient-echo signals were fitted with the two-stage fitting procedure to generate R2corr_t*. The signals were also fitted with the monoexponential and three-parameter models to generate R2nocorr* and R2corr*, respectively. Regions of interest (ROIs), including the corpus callosum, internal capsule, somatosensory cortex, caudo-putamen, thalamus, and lateral ventricle, were selected to evaluate the within-ROI mean and standard deviation (SD) of the R2* measurements. The results showed that the Akaike information criterion of the monoexponential model was significantly reduced by using the three-parameter model in the selected ROIs (p = 0.0039-0.0078). However, the within-ROI SD of R2corr* using the three-parameter model was significantly higher than that of the R2nocorr* in the internal capsule, caudo-putamen, and thalamus regions (p = 0.0039), a consequence partially due to the increased noise sensitivity of the three-parameter model. With the two-stage fitting procedure, the within-ROI SD of R2corr* was significantly reduced by 7.7-30.2% in all ROIs, except for the somatosensory cortex region with a fast in-plane variation of the B0 gradient field (p = 0.0039-0.0078). These results support the utilization of the two-stage fitting procedure to mitigate the B0 effect and reduce noise sensitivity for R2* measurement in the mouse brain.

Highly accelerated parameter mapping using model-based alternating reconstruction coupling fitting

Objective.A model-based alternating reconstruction coupling fitting, termed Model-based Alternating Reconstruction COupling fitting (MARCO), is proposed for accurate and fast magnetic resonance parameter mapping.Approach.MARCO utilizes the signal model as a regularization by minimizing the bias between the image series and the signal produced by the suitable signal model based on iteratively updated parameter maps when reconstructing. The technique can incorporate prior knowledge of both image series and parameters by adding sparsity constraints. The optimization problem is decomposed into three subproblems and solved through three alternating steps involving reconstruction and nonlinear least-square fitting, which can produce both contrast-weighted images and parameter maps simultaneously.Main results.The algorithm is applied toT2mapping with extended phase graph algorithm integrated and validated on undersampled multi-echo spin-echo data from both phantom and in vivo sources. Compared with traditional compressed sensing and model-based methods, the proposed approach yields more accurateT2maps with more details at high acceleration factors.Significance.The proposed method provides a basic framework for quantitative MR relaxometry, theoretically applicable to all quantitative MR relaxometry. It has the potential to improve the diagnostic utility of quantitative imaging techniques.

Joint MAPLE: Accelerated joint T1 and T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ mapping with scan-specific self-supervised networks

PURPOSE

Quantitative MRI finds important applications in clinical and research studies. However, it is encoding intensive and may suffer from prohibitively long scan times. Accelerated MR parameter mapping techniques have been developed to help address these challenges. Here, an accelerated joint T1,T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ , frequency and proton density mapping technique with scan-specific self-supervised network reconstruction is proposed to synergistically combine parallel imaging, model-based, and deep learning approaches to speed up parameter mapping.

METHODS

Proposed framework, Joint MAPLE, includes parallel imaging, signal modeling, and data consistency blocks which are optimized jointly in a combined loss function. A scan-specific self-supervised reconstruction is embedded into the framework, which takes advantage of multi-contrast data from a multi-echo, multi-flip angle, gradient echo acquisition.

RESULTS

In comparison with parallel reconstruction techniques powered by low-rank methods, emerging scan specific networks, and model-basedT 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ estimation approaches, the proposed framework reduces the reconstruction error in parameter maps by approximately two-fold on average at acceleration rates as high as R = 16 with uniform sampling. It can outperform evaluated parallel reconstruction techniques up to four-fold on average in the presence of challenging sub-sampling masks. It is observed that Joint MAPLE performs well at extreme acceleration rates of R = 25 and R = 36 with error values less than 20%.

CONCLUSION

Joint MAPLE enables higher fidelity parameter estimation at high acceleration rates by synergistically combining parallel imaging and model-based parameter mapping and exploiting multi-echo, multi-flip angle datasets. Utilizing a scan-specific self-supervised reconstruction obviates the need for large data sets for training while improving the parameter estimation ability.