Limits of Fat Quantification in the Presence of Iron Overload

Confidence maps for reliable estimation of proton density fat fraction and R 2 * in the liver

PURPOSE

The objective was to develop a fully automated algorithm that generates confidence maps to identify regions valid for analysis of quantitative proton density fat fraction (PDFF) andR 2 * $$ {R}_2^{\ast } $$ maps of the liver, generated with chemical shift-encoded MRI (CSE-MRI). Confidence maps are urgently needed for automated quality assurance, particularly with the emergence of automated segmentation and analysis algorithms.

METHODS

Confidence maps for both PDFF andR 2 * $$ {R}_2^{\ast } $$ maps are generated based on goodness of fit, measured by normalized RMS error between measured complex signals and the CSE-MRI signal model. Based on Cramér-Rao lower bound and Monte-Carlo simulations, normalized RMS error threshold criteria were developed to identify unreliable regions in quantitative maps. Simulation, phantom, and in vivo clinical studies were included. To analyze the clinical data, a board-certified radiologist delineated regions of interest (ROIs) in each of the nine liver segments for PDFF andR 2 * $$ {R}_2^{\ast } $$ analysis in consecutive clinical CSE-MRI data sets. The percent area of ROIs in areas deemed unreliable by confidence maps was calculated to assess the impact of confidence maps on real-world clinical PDFF andR 2 * $$ {R}_2^{\ast } $$ measurements.

RESULTS

Simulations and phantom studies demonstrated that the proposed algorithm successfully excluded regions with unreliable PDFF andR 2 * $$ {R}_2^{\ast } $$ measurements. ROI analysis by the radiologist revealed that 2.6% and 15% of the ROIs were placed in unreliable areas of PDFF andR 2 * $$ {R}_2^{\ast } $$ maps, as identified by confidence maps.

CONCLUSION

A proposed confidence map algorithm that identifies reliable areas of PDFF andR 2 * $$ {R}_2^{\ast } $$ measurements from CSE-MRI acquisitions was successfully developed. It demonstrated technical and clinical feasibility.

Simulation of a virtual liver iron overload model and R2 * estimation using multispectral fat-water models for GRE and UTE acquisitions at 1.5 T and 3 T

R2 *-MRI has emerged as a noninvasive alternative to liver biopsy for assessment of hepatic iron content (HIC). Multispectral fat-water R2 * modeling techniques such as the nonlinear least squares (NLSQ) fitting and autoregressive moving average (ARMA) models have been proposed for the accurate assessment of iron overload by also considering fat, which can otherwise confound R2 *-based HIC measurements in conditions of coexisting iron overload and steatosis. However, the R2 * estimation by these multispectral models has not been systematically investigated for various acquisition methods in iron overload only conditions and across the full clinically relevant range of HICs (0-40 mg Fe/g dry liver weight). The purpose of this study is to evaluate the R2 * accuracy and precision of multispectral models for various multiecho gradient echo (GRE) and ultrashort echo time (UTE) imaging acquisitions by constructing virtual iron overload models based on true histology and synthesizing MRI signals via Monte Carlo simulations at 1.5 T and 3 T, and comparing their results with monoexponential model and published in vivo R2 *-HIC calibrations. The signals were synthesized with TE1  = 1.0 ms for GRE and TE1  = 0.1 ms for UTE acquisition for varying echo spacing, ΔTE (0.1, 0.5, 1, 2 ms), and maximum echo time, TEmax (2, 4, 6, 10 ms). An iron-doped phantom study is also conducted to validate the simulation results in experimental GRE (TE1  = 1.2 ms, ΔTE  = 0.72 ms, TEmax  = 6.24 ms) and UTE (TE1  = 0.1 ms, ΔTE  = 0.5 ms, TEmax  = 6.1 ms) acquisitions. For GRE acquisitions, the multispectral ARMA and NLSQ models produced higher slopes (0.032-0.035) compared with the monoexponential model and published in vivo R2 *-HIC calibrations (0.025-0.028). However, for UTE acquisition for shorter echo spacing (≤0.5 ms) and longer maximum echo time, TEmax (≥6 ms), the multispectral and monoexponential signal models produced similar R2 *-HIC slopes (1.5 T, 0.028-0.032; 3 T, 0.014-0.016) and precision values (coefficient of variation < 25%) across the full clinical spectrum of HICs at both 1.5 T and 3 T. The phantom analysis also showed that all signal models demonstrated a significant improvement in R2 * estimation for UTE acquisition compared with GRE, confirming our simulation findings. Future work should investigate the performance of multispectral fat-water models by simulating liver models in coexisting conditions of iron overload and steatosis for accurate R2 * and fat quantification.

Two-point Dixon and six-point Dixon magnetic resonance techniques in the detection, quantification and grading of hepatic steatosis

BACKGROUND

Hepatic steatosis is a very common problem worldwide.

AIM

To assess the performance of two- and six-point Dixon magnetic resonance (MR) techniques in the detection, quantification and grading of hepatic steatosis.

METHODS

A single-center retrospective study was performed in 62 patients with suspected parenchymal liver disease. MR sequences included two-point Dixon, six-point Dixon, MR spectroscopy (MRS) and MR elastography. Fat fraction (FF) estimates on the Dixon techniques were compared to the MRS-proton density FF (PDFF). Statistical tests used included Pearson's correlation and receiver operating characteristic.

RESULTS

FF estimates on the Dixon techniques showed excellent correlation (≥ 0.95) with MRS-PDFF, and excellent accuracy [area under the receiver operating characteristic (AUROC) ≥ 0.95] in: (1) Detecting steatosis; and (2) Grading severe steatosis, (P < 0.001). In iron overload, two-point Dixon was not evaluable due to confounding T2* effects. FF estimates on six-point Dixon vs MRS-PDFF showed a moderate correlation (0.82) in iron overload vs an excellent correlation (0.97) without iron overload, (P < 0.03). The accuracy of six-point Dixon in grading mild steatosis improved (AUROC: 0.59 to 0.99) when iron overload cases were excluded. The excellent correlation (> 0.9) between the Dixon techniques vs MRS-PDFF did not change in the presence of liver fibrosis (P < 0.01).

CONCLUSION

Dixon techniques performed satisfactorily for the evaluation of hepatic steatosis but with exceptions.

Liver PDFF estimation using a multi-decoder water-fat separation neural network with a reduced number of echoes

OBJECTIVE

To accurately estimate liver PDFF from chemical shift-encoded (CSE) MRI using a deep learning (DL)-based Multi-Decoder Water-Fat separation Network (MDWF-Net), that operates over complex-valued CSE-MR images with only 3 echoes.

METHODS

The proposed MDWF-Net and a U-Net model were independently trained using the first 3 echoes of MRI data from 134 subjects, acquired with conventional 6-echoes abdomen protocol at 1.5 T. Resulting models were then evaluated using unseen CSE-MR images obtained from 14 subjects that were acquired with a 3-echoes CSE-MR pulse sequence with a shorter duration compared to the standard protocol. Resulting PDFF maps were qualitatively assessed by two radiologists, and quantitatively assessed at two corresponding liver ROIs, using Bland Altman and regression analysis for mean values, and ANOVA testing for standard deviation (STD) (significance level: .05). A 6-echo graph cut was considered ground truth.

RESULTS

Assessment of radiologists demonstrated that, unlike U-Net, MDWF-Net had a similar quality to the ground truth, despite it considered half of the information. Regarding PDFF mean values at ROIs, MDWF-Net showed a better agreement with ground truth (regression slope = 0.94, R2 = 0.97) than U-Net (regression slope = 0.86, R2 = 0.93). Moreover, ANOVA post hoc analysis of STDs showed a statistical difference between graph cuts and U-Net (p < .05), unlike MDWF-Net (p = .53).

CONCLUSION

MDWF-Net showed a liver PDFF accuracy comparable to the reference graph cut method, using only 3 echoes and thus allowing a reduction in the acquisition times.

CLINICAL RELEVANCE STATEMENT

We have prospectively validated that the use of a multi-decoder convolutional neural network to estimate liver proton density fat fraction allows a significant reduction in MR scan time by reducing the number of echoes required by 50%.

KEY POINTS

• Novel water-fat separation neural network allows for liver PDFF estimation by using multi-echo MR images with a reduced number of echoes. • Prospective single-center validation demonstrated that echo reduction leads to a significant shortening of the scan time, compared to standard 6-echo acquisition. • Qualitative and quantitative performance of the proposed method showed no significant differences in PDFF estimation with respect to the reference technique.

Quantitative MRI of diffuse liver diseases: techniques and tissue-mimicking phantoms

Quantitative magnetic resonance imaging (MRI) techniques are emerging as non-invasive alternatives to biopsy for assessment of diffuse liver diseases of iron overload, steatosis and fibrosis. For testing and validating the accuracy of these techniques, phantoms are often used as stand-ins to human tissue to mimic diffuse liver pathologies. However, currently, there is no standardization in the preparation of MRI-based liver phantoms for mimicking iron overload, steatosis, fibrosis or a combination of these pathologies as various sizes and types of materials are used to mimic the same liver disease. Liver phantoms that mimic specific MR features of diffuse liver diseases observed in vivo are important for testing and calibrating new MRI techniques and for evaluating signal models to accurately quantify these features. In this study, we review the liver morphology associated with these diffuse diseases, discuss the quantitative MR techniques for assessing these liver pathologies, and comprehensively examine published liver phantom studies and discuss their benefits and limitations.

Intra-examination agreement between multi-echo gradient echo and confounder-corrected chemical shift-encoded MR sequences for R2* estimation as a biomarker of liver iron content in patients with a wide range of T2*/R2* and proton density fat fraction values

PURPOSE

To investigate the intra-examination agreement between multi-echo gradient echo (MEGE) and confounder-corrected chemical shift-encoded (CSE) sequences for liver T2*/R2* estimations in a wide range of T2*/R2* and proton density fat fraction (PDFF) values. Exploratorily, to search for the T2*/R2* value where the agreement line breaks and examine differences between regions of low and high agreement.

METHODS

Consecutive patients at risk for liver iron overload who underwent MEGE and CSE sequences within the same exam at 1.5 T were retrospectively selected. Regions of interest were drawn in the right and one in the left liver lobes on post-processed images for R2*(sec-1) and PDFF (%) estimation. Agreement between MEGE-R2* and CSE-R2* was evaluated using intra-class correlation coefficient (ICC) and Bland-Altman analysis. 95% confidence intervals (CI) were computed. Segment-and-regression analysis was performed to find the point where the agreement between sequences is interrupted. Regions of low and high agreement were examined using tree-based partitioning analyses.

RESULTS

49 patients were included. Mean MEGE-R2* was 94.2 s-1 (range: 31.0-737.1) and mean CSE-R2* 87.7 (29.7-748.1). Mean CSE-PDFF was 9.12% (0.1-43.3). Agreement was strong for R2* estimations (ICC: 0.992,95%CI 0.987,0.996), but the relation was nonlinear and possibly heteroskedastic. Lower agreement occurred when MEGE-R2* > 235 s-1, with MEGE-R2* values consistently lower than CSE-R2*. Higher agreement was observed when PDFF < 14%.

CONCLUSION

MEGE-R2* and CSE-R2* strongly agree, though at higher iron content, MEGE-R2* is consistently lower than CSE-R2*. In this preliminary dataset, a breaking point for agreement was found at R2* > 235. Lower agreement was observed in patients with moderate to severe liver steatosis.

Quantification of Liver Iron Overload with MRI: Review and Guidelines from the ESGAR and SAR

Accumulation of excess iron in the body, or systemic iron overload, results from a variety of causes. The concentration of iron in the liver is linearly related to the total body iron stores and, for this reason, quantification of liver iron concentration (LIC) is widely regarded as the best surrogate to assess total body iron. Historically assessed using biopsy, there is a clear need for noninvasive quantitative imaging biomarkers of LIC. MRI is highly sensitive to the presence of tissue iron and has been increasingly adopted as a noninvasive alternative to biopsy for detection, severity grading, and treatment monitoring in patients with known or suspected iron overload. Multiple MRI strategies have been developed in the past 2 decades, based on both gradient-echo and spin-echo imaging, including signal intensity ratio and relaxometry strategies. However, there is a general lack of consensus regarding the appropriate use of these methods. The overall goal of this article is to summarize the current state of the art in the clinical use of MRI to quantify liver iron content and to assess the overall level of evidence of these various methods. Based on this summary, expert consensus panel recommendations on best practices for MRI-based quantification of liver iron are provided.

Ultrasound-Based Machine Learning Approach for Detection of Nonalcoholic Fatty Liver Disease

OBJECTIVES

Current diagnosis of nonalcoholic fatty liver disease (NAFLD) relies on biopsy or MR-based fat quantification. This prospective study explored the use of ultrasound with artificial intelligence for the detection of NAFLD.

METHODS

One hundred and twenty subjects with clinical suspicion of NAFLD and 10 healthy volunteers consented to participate in this institutional review board-approved study. Subjects were categorized as NAFLD and non-NAFLD according to MR proton density fat fraction (PDFF) findings. Ultrasound images from 10 different locations in the right and left hepatic lobes were collected following a standard protocol. MRI-based liver fat quantification was used as the reference standard with >6.4% indicative of NAFLD. A supervised machine learning model was developed for assessment of NAFLD. To validate model performance, a balanced testing dataset of 24 subjects was used. Sensitivity, specificity, positive predictive value, negative predictive value, and overall accuracy with 95% confidence interval were calculated.

RESULTS

A total of 1119 images from 106 participants was used for model development. The internal evaluation achieved an average precision of 0.941, recall of 88.2%, and precision of 89.0%. In the testing set AutoML achieved a sensitivity of 72.2% (63.1%-80.1%), specificity of 94.6% (88.7%-98.0%), positive predictive value (PPV) of 93.1% (86.0%-96.7%), negative predictive value of 77.3% (71.6%-82.1%), and accuracy of 83.4% (77.9%-88.0%). The average agreement for an individual subject was 92%.

CONCLUSIONS

An ultrasound-based machine learning model for identification of NAFLD showed high specificity and PPV in this prospective trial. This approach may in the future be used as an inexpensive and noninvasive screening tool for identifying NAFLD in high-risk patients.

Multicenter Reproducibility of Liver Iron Quantification with 1.5-T and 3.0-T MRI

Background MRI is a standard of care tool to measure liver iron concentration (LIC). Compared with regulatory-approved R2 MRI, R2* MRI has superior speed and is available in most MRI scanners; however, the cross-vendor reproducibility of R2*-based LIC estimation remains unknown. Purpose To evaluate the reproducibility of LIC via single-breath-hold R2* MRI at both 1.5 T and 3.0 T with use of a multicenter, multivendor study. Materials and Methods Four academic medical centers using MRI scanners from three different vendors (three 1.5-T scanners, one 2.89-T scanner, and two 3.0-T scanners) participated in this prospective cross-sectional study. Participants with known or suspected liver iron overload were recruited to undergo multiecho gradient-echo MRI for R2* mapping at 1.5 T and 3.0 T (2.89 T or 3.0 T) on the same day. R2* maps were reconstructed from the multiecho images and analyzed at a single center. Reference LIC measurements were obtained with a commercial R2 MRI method performed using standardized 1.5-T spin-echo imaging. R2*-versus-LIC calibrations were generated across centers and field strengths using linear regression and compared using F tests. Receiver operating characteristic (ROC) curve analysis was used to determine the diagnostic performance of R2* MRI in the detection of clinically relevant LIC thresholds. Results A total of 207 participants (mean age, 38 years ± 20 [SD]; 117 male participants) were evaluated between March 2015 and September 2019. A linear relationship was confirmed between R2* and LIC. All calibrations within the same field strength were highly reproducible, showing no evidence of statistically significant center-specific differences (P > .43 across all comparisons). Calibrations for 1.5 T and 3.0 T were generated, as follows: for 1.5 T, LIC (in milligrams per gram [dry weight]) = -0.16 + 2.603 × 10-2 R2* (in seconds-1); for 2.89 T, LIC (in milligrams per gram) = -0.03 + 1.400 × 10-2 R2* (in seconds-1); for 3.0 T, LIC (in milligrams per gram) = -0.03 + 1.349 × 10-2 R2* (in seconds-1). Liver R2* had high diagnostic performance in the detection of clinically relevant LIC thresholds (area under the ROC curve, >0.98). Conclusion R2* MRI enabled accurate and reproducible quantification of liver iron overload over clinically relevant ranges of liver iron concentration (LIC). The data generated in this study provide the necessary calibrations for broad clinical dissemination of R2*-based LIC quantification. ClinicalTrials.gov registration no.: NCT02025543 © RSNA, 2022 Online supplemental material is available for this article.

Hepatic Iron Quantification Using a Free-Breathing 3D Radial Gradient Echo Technique and Validation With a 2D Biopsy-Calibrated R2* Relaxometry Method

BACKGROUND

Hepatic iron content (HIC) is an important parameter for the management of iron overload. Non-invasive HIC assessment is often performed using biopsy-calibrated two-dimensional breath-hold Cartesian gradient echo (2D BH GRE) R2* -MRI. However, breath-holding is not possible in most pediatric patients or those with respiratory problems, and three-dimensional free-breathing radial GRE (3D FB rGRE) has emerged as a viable alternative.

PURPOSE

To evaluate the performance of a 3D FB rGRE and validate its R2* and fat fraction (FF) quantification with 3D breath-hold Cartesian GRE (3D BH cGRE) and biopsy-calibrated 2D BH GRE across a wide range of HICs.

STUDY TYPE

Retrospective.

SUBJECTS

Twenty-nine patients with hepatic iron overload (22 females, median age: 15 [5-25] years).

FIELD STRENGTH/SEQUENCE

Three-dimensional radial and 2D and 3D Cartesian multi-echo GRE at 1.5 T.

ASSESSMENT

R2* and FF maps were computed for 3D GREs using a multi-spectral fat model and 2D GRE R2* maps were calculated using a mono-exponential model. Mean R2* and FF values were calculated via whole-liver contouring and T2* -thresholding by three operators.

STATISTICAL TESTS

Inter- and intra-observer reproducibility was assessed using Bland-Altman and intraclass correlation coefficient (ICC). Linear regression and Bland-Altman analysis were performed to compare R2* and FF values among the three acquisitions. One-way repeated-measures ANOVA and Wilcoxon signed-rank tests, respectively, were used to test for significant differences between R2* and FF values obtained with different acquisitions. Statistical significance was assumed at P < 0.05.

RESULTS

The mean biases and ICC for inter- and intra-observer reproducibility were close to 0% and >0.99, respectively for both R2* and FF. The 3D FB rGRE R2* and FF values were not significantly different (P > 0.44) and highly correlated (R2  ≥ 0.98) with breath-hold Cartesian GREs, with mean biases ≤ ±2.5% and slopes 0.90-1.12. In non-breath-holding patients, Cartesian GREs showed motion artifacts, whereas 3D FB rGRE exhibited only minimal streaking artifacts.

DATA CONCLUSION

Free-breathing 3D radial GRE is a viable alternative in non-breath-hold patients for accurate HIC estimation.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 2.