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.

Quantitative Susceptibility Mapping Using a Multispectral Autoregressive Moving Average Model to Assess Hepatic Iron Overload.

BACKGROUND: R2*-MRI is clinically used to noninvasively assess hepatic iron content (HIC) to guide potential iron chelation therapy. However, coexisting pathologies, such as fibrosis and steatosis, affect R2* measurements and may thus confound HIC estimations. PURPOSE: To evaluate whether a multispectral auto regressive moving average (ARMA) model can be used in conjunction with quantitative susceptibility mapping (QSM) to measure magnetic susceptibility as a confounder-free predictor of HIC. STUDY TYPE: Phantom study and in vivo cohort. SUBJECTS: Nine iron phantoms covering clinically relevant R2* range (20-1200/second) and 48 patients (22 male, 26 female, median age 18 years). FIELD STRENGTH/SEQUENCE: Three-dimensional (3D) and two-dimensional (2D) multi-echo gradient echo (GRE) at 1.5 T. ASSESSMENT: ARMA-QSM modeling was performed on the complex 3D GRE signal to estimate R2*, fat fraction (FF), and susceptibility measurements. R2*-based dry clinical HIC values were calculated from the 2D GRE acquisition using a published R2*-HIC calibration curve as reference standard. STATISTICAL TESTS: Linear regression analysis was performed to compare ARMA R2* and susceptibility-based estimates to iron concentrations and dry clinical HIC values in phantoms and patients, respectively. RESULTS: In phantoms, the ARMA R2* and susceptibility values strongly correlated with iron concentrations (R2 ≥ 0.9). In patients, the ARMA R2* values highly correlated (R2  = 0.97) with clinical HIC values with slope = 0.026, and the susceptibility values showed good correlation (R2  = 0.82) with clinical dry HIC values with slope = 3.3 and produced a dry-to-wet HIC ratio of 4.8. DATA CONCLUSION: This study shows the feasibility that ARMA-QSM can simultaneously estimate susceptibility-based wet HIC, R2*-based dry HIC and FFs from a single multi-echo GRE acquisition. Our results demonstrate that both, R2* and susceptibility-based wet HIC values estimated with ARMA-QSM showed good association with clinical dry HIC values with slopes similar to published R2*-biopsy HIC calibration and dry-to-wet tissue weight ratio, respectively. Hence, our study shows that ARMA-QSM can provide potentially confounder-free assessment of hepatic iron overload. LEVEL OF EVIDENCE: 3 TECHNICAL EFFICACY: Stage 2.

Autoregressive moving average modeling for hepatic iron quantification in the presence of fat.

BACKGROUND: Measuring hepatic R2* by fitting a monoexponential model to the signal decay of a multigradient-echo (mGRE) sequence noninvasively determines hepatic iron content (HIC). Concurrent hepatic steatosis introduces signal oscillations and confounds R2* quantification with standard monoexponential models. PURPOSE: To evaluate an autoregressive moving average (ARMA) model for accurate quantification of HIC in the presence of fat using biopsy as the reference. STUDY TYPE: Phantom study and in vivo cohort. POPULATION: Twenty iron-fat phantoms covering clinically relevant R2* (30-800 s-1 ) and fat fraction (FF) ranges (0-40%), and 10 patients (four male, six female, mean age 18.8 years). FIELD STRENGTH/SEQUENCE: 2D mGRE acquisitions at 1.5 T and 3 T. ASSESSMENT: Phantoms were scanned at both field strengths. In vivo data were analyzed using the ARMA model to determine R2* and FF values, and compared with biopsy results. STATISTICAL TESTS: Linear regression analysis was used to compare ARMA R2* and FF results with those obtained using a conventional monoexponential model, complex-domain nonlinear least squares (NLSQ) fat-water model, and biopsy. RESULTS: In phantoms and in vivo, all models produced R2* and FF values consistent with expected values in low iron and low/high fat conditions. For high iron and no fat phantoms, monoexponential and ARMA models performed excellently (slopes: 0.89-1.07), but NLSQ overestimated R2* (slopes: 1.14-1.36) and produced false FFs (12-17%) at 1.5 T; in high iron and fat phantoms, NLSQ (slopes: 1.02-1.16) outperformed monoexponential and ARMA models (slopes: 1.23-1.88). The results with NLSQ and ARMA improved in phantoms at 3 T (slopes: 0.96-1.04). In patients, mean R2*-HIC estimates for monoexponential and ARMA models were close to biopsy-HIC values (slopes: 0.90-0.95), whereas NLSQ substantially overestimated HIC (slope 1.4) and produced false FF values (4-28%) with very high SDs (15-222%) in patients with high iron overload and no steatosis. DATA CONCLUSION: ARMA is superior in quantifying R2* and FF under high iron and no fat conditions, whereas NLSQ is superior for high iron and concurrent fat at 1.5 T. Both models give improved R2* and FF results at 3 T. LEVEL OF EVIDENCE: 2 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2019;50:1620-1632.

Ultrashort echo time imaging for quantification of hepatic iron overload: Comparison of acquisition and fitting methods via simulations, phantoms, and in vivo data.

BACKGROUND: Current R2*-MRI techniques for measuring hepatic iron content (HIC) use various acquisition types and fitting models. PURPOSE: To evaluate the accuracy and precision of R2*-HIC acquisition and fitting methods. STUDY TYPE: Signal simulations, phantom study, and prospective in vivo cohort. POPULATION: In all, 132 patients (58/74 male/female, mean age 17.7 years). FIELD STRENGTH/SEQUENCE: 2D-multiecho gradient-echo (GRE) and ultrashort echo time (UTE) acquisitions at 1.5T. ASSESSMENT: Synthetic MR signals were created to mimic published GRE and UTE methods, using different R2* values (25-2000 s-1 ) and signal-to-noise ratios (SNR). Phantoms with varying iron concentrations were scanned at 1.5T. In vivo data were analyzed from 132 patients acquired at 1.5T. R2* was estimated by fitting using three signal models. Accuracy and precision of R2* measurements for UTE acquisition parameters (SNR, echo spacing [ΔTE], maximum echo time [TEmax ]) and fitting methods were compared for simulated, phantom, and in vivo datasets. STATISTICAL TESTS: R2* accuracy was determined from the relative error and by linear regression analysis. Precision was evaluated using coefficient of variation (CoV) analysis. RESULTS: In simulations, all models had high R2* accuracy (error <5%) and precision (CoV <10%) for all SNRs, shorter ΔTE (≤0.5 msec), and longer TEmax (≥10.1 msec); except the constant offset model overestimated R2* at the lowest SNR. In phantoms and in vivo, all models produced similar R2* values for different SNRs and shorter ΔTEs (slopes: 0.99-1.06, R2 > 0.99, P < 0.001). In all experiments, R2* results degraded for high R2* values with longer ΔTE (≥1 msec). In vivo, shorter and longer TEmax gave similar R2* results (slopes: 1.02-1.06, R2 > 0.99, P < 0.001) for the noise subtraction model for 25≤R2*≤2000 s-1 . However, both quadratic and constant offset models, using shorter TEmax (≤4.7 msec) overestimated R2* and yielded high CoVs up to ∼170% for low R2* (<250 s-1 ). DATA CONCLUSION: UTE with TEmax ≥ 10.1 msec and ΔTE ≤ 0.5 msec yields accurate R2* estimates over the entire clinical HIC range. Monoexponential fitting with noise subtraction is the most robust signal model to changes in UTE parameters and achieves the highest R2* accuracy and precision. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:1475-1488.

Fast quantitative parameter maps without fitting: Integration yields accurate mono-exponential signal decay rates.

PURPOSE: To develop a computationally fast and accurate algorithm for mono-exponential signal modelling and validate the new technique in the context of R2* mapping for iron overload assessment. METHODS: An algorithm is introduced that directly calculates R2* values from a series of images based on integration of the mono-exponential signal decay curve. The algorithm is fast, because fitting is avoided and only arithmetic computations without iterations are applied. Precision and accuracy of the method is determined in comparison to the conventional log-linear (LL), nonlinear least-squares-based Levenberg-Marquardt (NLM), and squared nonlinear Levenberg-Marquardt (SQNLM) methods, which rely on iterative curve fitting. RESULTS: In simulations, the signal integration based method consistently had the same or better accuracy than the LL, NLM, and SQNLM algorithms for R2* values ranging from 50 s-1 to 1200 s-1 . In phantoms and in vivo (12 participants), this method was robust over a wide range of R2* values and signal-to-noise ratios. Computation times were approximately 100, 1460, and 930 times faster than those of the LL, NLM, and SQNLM methods, respectively. CONCLUSIONS: The fast signal integration method accurately calculates R2* maps. It has the potential to replace conventional, mono-exponential fitting methods for quantitative MRI such as R2* parameter mapping. Magn Reson Med 79:2978-2985, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Automated vessel exclusion technique for quantitative assessment of hepatic iron overload by R2*-MRI.

BACKGROUND: Extraction of liver parenchyma is an important step in the evaluation of R2*-based hepatic iron content (HIC). Traditionally, this is performed by radiologists via whole-liver contouring and T2*-thresholding to exclude hepatic vessels. However, the vessel exclusion process is iterative, time-consuming, and susceptible to interreviewer variability. PURPOSE: To implement and evaluate an automatic hepatic vessel exclusion and parenchyma extraction technique for accurate assessment of R2*-based HIC. STUDY TYPE: Retrospective analysis of clinical data. SUBJECTS: Data from 511 MRI exams performed on 257 patients were analyzed. FIELD STRENGTH/SEQUENCE: All patients were scanned on a 1.5T scanner using a multiecho gradient echo sequence for clinical monitoring of HIC. ASSESSMENT: An automated method based on a multiscale vessel enhancement filter was investigated for three input data types-contrast-optimized composite image, T2* map, and R2* map-to segment blood vessels and extract liver tissue for R2*-based HIC assessment. Segmentation and R2* results obtained using this automated technique were compared with those from a reference T2*-thresholding technique performed by a radiologist. STATISTICAL TESTS: The Dice similarity coefficient was used to compare the segmentation results between the extracted parenchymas, and linear regression and Bland-Altman analyses were performed to compare the R2* results, obtained with the automated and reference techniques. RESULTS: Mean liver R2* values estimated from all three filter-based methods showed excellent agreement with the reference method (slopes 1.04-1.05, R2 > 0.99, P < 0.001). Parenchyma areas extracted using the reference and automated methods had an average overlap area of 87-88%. The T2*-thresholding technique included small vessels and pixels at the vessel/tissue boundaries as parenchymal area, potentially causing a small bias (<5%) in R2* values compared to the automated method. DATA CONCLUSION: The excellent agreement between reference and automated hepatic vessel segmentation methods confirms the accuracy and robustness of the proposed method. This automated approach might improve the radiologist's workflow by reducing the interpretation time and operator dependence for assessing HIC, an important clinical parameter that guides iron overload management. LEVEL OF EVIDENCE: 3 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;47:1542-1551.

Quantitative ultrashort echo time imaging for assessment of massive iron overload at 1.5 and 3 Tesla.

PURPOSE: Hepatic iron content (HIC) quantification via transverse relaxation rate (R2*)-MRI using multi-gradient echo (mGRE) imaging is compromised toward high HIC or at higher fields due to the rapid signal decay. Our study aims at presenting an optimized 2D ultrashort echo time (UTE) sequence for R2* quantification to overcome these limitations. METHODS: Two-dimensional UTE imaging was realized via half-pulse excitation and radial center-out sampling. The sequence includes chemically selective saturation pulses to reduce streaking artifacts from subcutaneous fat, and spatial saturation (sSAT) bands to suppress out-of-slice signals. The sequence employs interleaved multi-echo readout trains to achieve dense temporal sampling of rapid signal decays. Evaluation was done at 1.5 Tesla (T) and 3T in phantoms, and clinical applicability was demonstrated in five patients with biopsy-confirmed massively high HIC levels (>25 mg Fe/g dry weight liver tissue). RESULTS: In phantoms, the sSAT pulses were found to remove out-of-slice contamination, and R2* results were in excellent agreement to reference mGRE R2* results (slope of linear regression: 1.02/1.00 for 1.5/3T). UTE-based R2* quantification in patients with massive iron overload proved successful at both field strengths and was consistent with biopsy HIC values. CONCLUSION: The UTE sequence provides a means to measure R2* in patients with massive iron overload, both at 1.5T and 3T. Magn Reson Med 78:1839-1851, 2017. © 2017 Wiley Periodicals, Inc.

Radial Ultrashort TE Imaging Removes the Need for Breath-Holding in Hepatic Iron Overload Quantification by R2* MRI.

OBJECTIVE: The objective of this study is to evaluate radial free-breathing (FB) multiecho ultrashort TE (UTE) imaging as an alternative to Cartesian FB multiecho gradient-recalled echo (GRE) imaging for quantitative assessment of hepatic iron content (HIC) in sedated patients and subjects unable to perform breath-hold (BH) maneuvers. MATERIALS AND METHODS: FB multiecho GRE imaging and FB multiecho UTE imaging were conducted for 46 test group patients with iron overload who could not complete BH maneuvers (38 patients were sedated, and eight were not sedated) and 16 control patients who could complete BH maneuvers. Control patients also underwent standard BH multiecho GRE imaging. Quantitative R2* maps were calculated, and mean liver R2* values and coefficients of variation (CVs) for different acquisitions and patient groups were compared using statistical analysis. RESULTS: FB multiecho GRE images displayed motion artifacts and significantly lower R2* values, compared with standard BH multiecho GRE images and FB multiecho UTE images in the control cohort and FB multiecho UTE images in the test cohort. In contrast, FB multiecho UTE images produced artifact-free R2* maps, and mean R2* values were not significantly different from those measured by BH multiecho GRE imaging. Motion artifacts on FB multiecho GRE images resulted in an R2* CV that was approximately twofold higher than the R2* CV from BH multiecho GRE imaging and FB multiecho UTE imaging. The R2* CV was relatively constant over the range of R2* values for FB multiecho UTE, but it increased with increases in R2* for FB multiecho GRE imaging, reflecting that motion artifacts had a stronger impact on R2* estimation with increasing iron burden. CONCLUSION: FB multiecho UTE imaging was less motion sensitive because of radial sampling, produced excellent image quality, and yielded accurate R2* estimates within the same acquisition time used for multiaveraged FB multiecho GRE imaging. Thus, FB multiecho UTE imaging is a viable alternative for accurate HIC assessment in sedated children and patients who cannot complete BH maneuvers.

Can multi-slice or navigator-gated R2* MRI replace single-slice breath-hold acquisition for hepatic iron quantification?

BACKGROUND: Liver R2* values calculated from multi-gradient echo (mGRE) magnetic resonance images (MRI) are strongly correlated with hepatic iron concentration (HIC) as shown in several independently derived biopsy calibration studies. These calibrations were established for axial single-slice breath-hold imaging at the location of the portal vein. Scanning in multi-slice mode makes the exam more efficient, since whole-liver coverage can be achieved with two breath-holds and the optimal slice can be selected afterward. Navigator echoes remove the need for breath-holds and allow use in sedated patients. OBJECTIVE: To evaluate if the existing biopsy calibrations can be applied to multi-slice and navigator-controlled mGRE imaging in children with hepatic iron overload, by testing if there is a bias-free correlation between single-slice R2* and multi-slice or multi-slice navigator controlled R2*. MATERIALS AND METHODS: This study included MRI data from 71 patients with transfusional iron overload, who received an MRI exam to estimate HIC using gradient echo sequences. Patient scans contained 2 or 3 of the following imaging methods used for analysis: single-slice images (n = 71), multi-slice images (n = 69) and navigator-controlled images (n = 17). Small and large blood corrected region of interests were selected on axial images of the liver to obtain R2* values for all data sets. Bland-Altman and linear regression analysis were used to compare R2* values from single-slice images to those of multi-slice images and navigator-controlled images. RESULTS: Bland-Altman analysis showed that all imaging method comparisons were strongly associated with each other and had high correlation coefficients (0.98 ≤ r ≤ 1.00) with P-values ≤0.0001. Linear regression yielded slopes that were close to 1. CONCLUSION: We found that navigator-gated or breath-held multi-slice R2* MRI for HIC determination measures R2* values comparable to the biopsy-validated single-slice, single breath-hold scan. We conclude that these three R2* methods can be interchangeably used in existing R2*-HIC calibrations.