Evidence of Tri-Exponential Decay for Liver Intravoxel Incoherent Motion MRI: A Review of Published Results and Limitations

Reproducibility of intravoxel incoherent motion quantification in the liver across field strengths and gradient hardware

PURPOSE

To evaluate reproducibility and interlobar agreement of intravoxel incoherent motion (IVIM) quantification in the liver across field strengths and MR scanners with different gradient hardware.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimized monopolar and motion-robust 2D (b-value and first-order motion moment [M1]) IVIM-DWI acquisitions. Eleven healthy volunteers underwent diffusion MRI of the liver, where each optimized acquisition was obtained five times across three MRI scanners. For each data set, IVIM estimates (diffusion coefficient (D), pseudo-diffusion coefficients (d 1 * $$ {d}_1^{\ast } $$ andd 2 * $$ {d}_2^{\ast } $$ ), blood velocity SDs (Vb1 and Vb2), and perfusion fractions [f1 and f2]) were obtained in the right and left liver lobes using two signal models (pseudo-diffusion and M1-dependent physical) with and without T2 correction (fc1 and fc2) and three fitting techniques (tri-exponential region of interest-based full and segmented fitting and blood velocity SD distribution fitting). Reproducibility and interlobar agreement were compared across methods using within-subject and pairwise coefficients of variation (CVw and CVp), paired sample t-tests, and Bland-Altman analysis.

RESULTS

Using a combination of motion-robust 2D (b-M1) data acquisition, M1-dependent physical signal modeling with T2 correction, and blood velocity SD distribution fitting, multiscanner reproducibility with median CVw = 5.09%, 11.3%, 9.20%, 14.2%, and 12.6% for D, Vb1, Vb2, fc1, and fc2, respectively, and interlobar agreement with CVp = 8.14%, 11.9%, 8.50%, 49.9%, and 42.0%, respectively, was achieved.

CONCLUSION

Recently proposed advanced IVIM acquisition, signal modeling, and fitting techniques may facilitate reproducible IVIM quantification in the liver, as needed for establishment of IVIM-based quantitative biomarkers for detection, staging, and treatment monitoring of diseases.

Repeatability and reproducibility comparisons of liver IVIM imaging with free-breathing or respiratory-triggered sequences

For liver intravoxel incoherent motion (IVIM) data acquisition, respiratory-triggering (RT) MRI is commonly used, and there are strong motivations to shorten the scan duration. For the same scan duration, more b values or higher numbers of excitations can be allowed for free-breathing (FB) imaging than for RT. We studied whether FB can be used to replace RT when careful IVIM image acquisition and image processing are conducted. MRI data of 22 healthy participants were acquired using a 3.0 T scanner. Diffusion imaging was based on a single-shot spin-echo-type echo-planar sequence and 16 b values of 0, 2, 4, 7, 10, 15, 20, 30, 46, 60, 72, 100, 150, 200, 400, and 600 s/mm2 . Each subject attended two scan sessions with an interval of 10-20 days. For each scan session, a subject was scanned twice, first with RT and then with FB. The mean image acquisition time was 5.4 min for FB and 10.8 min for RT. IVIM parameters were calculated with bi-exponential model segmented fitting with a threshold b value of 60 s/mm2 , and fitting started from b = 2 s/mm2 . There was no statistically significant difference between IVIM parameters measured with FB imaging or RT imaging. Perfusion fraction ICC (intraclass correlation coefficient) for FB imaging and RT imaging in the same scan session was 0.824. For perfusion fraction, wSD (within-subject standard deviation), BA (Bland-Altman) difference, BA 95% limit, and ICC were 0.022, 0.0001, -0.0635~0.0637, and 0.687 for FB and 0.031, 0.0122, -0.0723~0.0967, and 0.611 for RT. For Dslow (×10-3  s/mm2 ), wSD, BA difference, BA 95% limit, and ICC were 0.057, 0.0268, -0.1258~0.1793, and 0.471 for FB and 0.073, -0.0078, -0.2170-0.2014, and <0.4 for RT. The Dfast coefficient of variation was 0.20 for FB imaging and 0.28 for RT imaging. All reproducibility indicators slightly favored FB imaging.

Reproducibility and Repeatability of Intravoxel Incoherent Motion MRI Acquisition Methods in Liver

BACKGROUND

Intravoxel incoherent motion (IVIM) diffusion weighted MRI (DWI) has potential for evaluating hepatic fibrosis but image acquisition technique influence on diffusion parameter estimation bears investigation.

PURPOSE

To minimize variability and maximize repeatably in abdominal DWI in terms of IVIM parameter estimates.

STUDY TYPE

Prospective test-retest and image quality comparison.

SUBJECTS

Healthy volunteers (3F/7M, 29.9 ± 12.9 years) and Family Study subjects (18F/12M, 51.7 ± 16.7 years), without and with liver steatosis.

FIELD STRENGTH/SEQUENCE

Abdominal single-shot echo-planar imaging (EPI) and simultaneous multi-slice (SMS) DWI sequences with respiratory triggering (RT), breath-holding (BH), and navigator echo (NE) at 3 Tesla.

ASSESSMENT

SMS-BH, EPI-NE, and SMS-RT data from twice-scanned healthy volunteers were analyzed using 6 × b-values (0-800 s⋅mm-2 ) and lower (LO) and higher (HI) b-value ranges. Family Study subjects were scanned using SMS and standard EPI sequences. The biexponential IVIM model was used to estimate fast-diffusion coefficient (Df ), fraction of fast diffusion (f), and slow-diffusion coefficient (Ds ). Scan time, estimated signal-to-noise ratio (eSNR), eSNR per acquisition, and distortion ratio were compared.

STATISTICAL TESTS

Coefficients of variation (CoV) and Bland Altman analyses were performed for test-retest repeatability. Interclass correlation coefficient (ICC) assessed interobserver agreement with P < 0.05 deemed significant.

RESULTS

Within-subject CoVs among volunteers (N = 10) for f and Ds were lowest in EPI-NE-LO (11.6%) and SMS-RT-HI (11.1%). Inter-observer ICCs for f and Ds were highest for EPI-NE-LO (0.63) and SMS-RT-LO (0.76). Df could not be estimated for most subjects. Estimated eSNR (EPI = 21.9, SMS = 4.7) and eSNR time (EPI = 6.7, SMS = 16.6) were greater for SMS, with less distortion in the liver region (DR-PE: EPI = 23.6, SMS = 13.1).

DATA CONCLUSION

Simultaneous multislice acquisitions had significantly less variability and higher ICCs of Ds , higher eSNR, less distortion, and reduced scan time compared to EPI.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

The Feasibility of Using Tri-Exponential Intra-Voxel Incoherent Motion DWI for Identifying the Microvascular Invasion in Hepatocellular Carcinoma

Purpose

To assess the effectiveness of tri-exponential Intra-Voxel Incoherent Motion (tri-IVIM) MRI in preoperatively identifying microvascular invasion (MVI) in hepatocellular carcinoma (HCC).

Patients and Methods

In this prospective study, 67 patients with HCC were included. Metrics from bi-exponential IVIM (bi-IVIM) and tri-IVIM were calculated. Subgroup comparisons were analyzed using the independent Student's t-test or Mann-Whitney U-test. Logistic regression was performed to explore clinical risk factors. Diagnostic performance was assessed using receiver operating characteristic (ROC) curves, calibration curves and decision curve analysis.

Results

MVI-positive HCCs exhibited significantly lower true diffusion coefficient (Dt) from bi-IVIM, as well as fast-diffusion coefficients (Df) and slow-diffusion coefficients (Ds) from tri-IVIM, compared to MVI-negative HCCs (p < 0.05). Tumor size and alpha-fetoprotein (AFP) were identified as risk factors. The combination of tri-IVIM-derived metrics (Ds and Df) yielded higher diagnostic accuracy (AUC = 0.808) compared to bi-IVIM (AUC = 0.741). A predictive model based on a nomogram was constructed using Ds, Df, tumor size, and AFP, resulting in the highest diagnostic accuracy (AUC = 0.859). Decision curve analysis indicated that the constructed model, provided the highest net benefit by accurately stratifying the risk of MVI, followed by tri-IVIM and bi-IVIM.

Conclusion

Tri-IVIM can provide information on perfusion and diffusion for evaluating MVI in HCC. Additionally, tri-IVIM outperformed bi-IVIM in identifying MVI-positive HCC. By integrating clinical risk factors and metrics from tri-IVIM, a predictive nomogram exhibited the highest diagnostic accuracy, potentially aiding in the noninvasive and preoperative assessment of MVI.

Usefulness of intravoxel incoherent motion MRI for visualizing slow cerebrospinal fluid motion

BACKGROUND

In the cerebrospinal fluid (CSF) dynamics, the pulsations of cerebral arteries and brain is considered the main driving force for the reciprocating bidirectional CSF movements. However, measuring these complex CSF movements on conventional flow-related MRI methods is difficult. We tried to visualize and quantify the CSF motion by using intravoxel incoherent motion (IVIM) MRI with low multi-b diffusion-weighted imaging.

METHODS

Diffusion-weighted sequence with six b values (0, 50, 100, 250, 500, and 1000 s/mm²) was performed on 132 healthy volunteers aged ≥ 20 years and 36 patients with idiopathic normal pressure hydrocephalus (iNPH). The healthy volunteers were divided into three age groups (< 40, 40 to < 60, and ≥ 60 years). In the IVIM analysis, the bi-exponential IVIM fitting method using the Levenberg-Marquardt algorithm was adapted. The average, maximum, and minimum values of ADC, D, D*, and fraction of incoherent perfusion (f) calculated by IVIM were quantitatively measured in 45 regions of interests in the whole ventricles and subarachnoid spaces.

RESULTS

Compared with healthy controls aged ≥ 60 years, the iNPH group had significantly lower mean f values in all the parts of the lateral and 3rd ventricles, whereas significantly higher mean f value in the bilateral foramina of Luschka. In the bilateral Sylvian fossa, which contain the middle cerebral bifurcation, the mean f values increased gradually with increasing age, whereas those were significantly lower in the iNPH group. In the 45 regions of interests, the f values in the bilateral foramina of Luschka were the most positively correlated with the ventricular size and indices specific to iNPH, whereas that in the anterior part of the 3rd ventricle was the most negatively correlated with the ventricular size and indices specific to iNPH. Other parameters of ADC, D, and D* were not significantly different between the two groups in any locations.

CONCLUSIONS

The f value on IVIM MRI is useful for evaluating small pulsatile complex motion of CSF throughout the intracranial CSF spaces. Patients with iNPH had significantly lower mean f values in the whole lateral ventricles and 3rd ventricles and significantly higher mean f value in the bilateral foramina of Luschka, compared with healthy controls aged ≥ 60 years.

Influence of saturation effects on biexponential liver intravoxel incoherent motion

PURPOSE

Studies on intravoxel incoherent motion (IVIM) imaging in the liver have been carried out with different acquisition protocols. The number of acquired slices and the distances between slices can influence IVIM measurements due to saturation effects, but these effects have often been disregarded. This study investigated differences in biexponential IVIM parameters between two slice settings.

METHODS

Fifteen healthy volunteers (21-30 years) were examined at a field strength of 3 T. Diffusion-weighted images of the abdomen were acquired with 16 b values (0-800 s/mm² ), with four slices for the few slices setting and 24-27 slices for the many slices setting. Regions of interest were manually drawn in the liver. The data were fitted with a monoexponential signal curve and a biexponential IVIM curve, and biexponential IVIM parameters were determined. The dependence on the slice setting was assessed with Student's t test for paired samples (normally distributed IVIM parameters) and the Wilcoxon signed-rank test (non-normally distributed parameters).

RESULTS

None of the parameters were significantly different between the settings. For few slices and many slices, respectively, the mean values (SDs) for D $$ D $$ were 1.21 μm 2 / ms $$ 1.21{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.19 μm 2 / ms $$ 0.19\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ) and 1.20 μm 2 / ms $$ 1.20{\upmu \mathrm{m}}^2/\mathrm{ms} $$ ( 0.11 μm 2 / ms $$ 0.11\kern0.3em {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ); for f $$ f $$ they were 29.7% (6.2%) and 27.7% (3.6%); and for D * $$ {D}^{\ast } $$ they were 8.76 ⋅ 10 - 2 mm 2 / s $$ 8.76\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.54 ⋅ 10 - 2 mm 2 / s $$ 4.54\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ) and 8.71 ⋅ 10 - 2 mm 2 / s $$ 8.71\cdot {10}^{-2}{\mathrm{mm}}^2/\mathrm{s} $$ ( 4.06 ⋅ 10 - 2 mm 2 / s $$ 4.06\cdot {10}^{-2}\kern0.3em {\mathrm{mm}}^2/\mathrm{s} $$ ).

CONCLUSION

Biexponential IVIM parameters in the liver are comparable among IVIM studies that use different slice settings, with mostly negligible saturation effects. However, this may not hold for studies that use much shorter TR.

Precision of region of interest-based tri-exponential intravoxel incoherent motion quantification and the role of the Intervoxel spatial distribution of flow velocities

PURPOSE

The purpose of this work was to obtain precise tri-exponential intravoxel incoherent motion (IVIM) quantification in the liver using 2D (b-value and first-order motion moment [M1 ]) IVIM-DWI acquisitions and region of interest (ROI)-based fitting techniques.

METHODS

Diffusion MRI of the liver was performed in 10 healthy volunteers using three IVIM-DWI acquisitions: conventional monopolar, optimized monopolar, and optimized 2D (b-M1 ). For each acquisition, bi-exponential and tri-exponential full, segmented, and over-segmented ROI-based fitting and a newly proposed blood velocity SDdistribution (BVD) fitting technique were performed to obtain IVIM estimates in the right and left liver lobes. Fitting quality was evaluated using corrected Akaike information criterion. Precision metrics (test-retest repeatability, inter-reader reproducibility, and inter-lobar agreement) were evaluated using Bland-Altman analysis, repeatability/reproducibility coefficients (RPCs), and paired sample t-tests. Precision was compared across acquisitions and fitting methods.

RESULTS

High repeatability and reproducibility was observed in the estimations of the diffusion coefficient (Dtri  = [1.03 ± 0.11] × 10-3  mm2 /s; RPCs ≤ 1.34 × 10-4  mm2 /s), perfusion fractions (F1  = 3.19 ± 1.89% and F2  = 16.4 ± 2.07%; RPCs ≤ 2.51%), and blood velocity SDs (Vb,1  = 1.44 ± 0.14 mm/s and Vb,2  = 3.62 ± 0.13 mm/s; RPCs ≤ 0.41 mm/s) in the right liver lobe using the 2D (b-M1 ) acquisition in conjunction with BVD fitting. Using these methods, significantly larger (p < 0.01) estimates of Dtri and F1 were observed in the left lobe in comparison to the right lobe, while estimates of Vb,1 and Vb,2 demonstrated high interlobar agreement (RPCs ≤ 0.45 mm/s).

CONCLUSIONS

The 2D (b-M1 ) IVIM-DWI data acquisition in conjunction with BVD fitting enables highly precise tri-exponential IVIM quantification in the right liver lobe.

Self-supervised neural network improves tri-exponential intravoxel incoherent motion model fitting compared to least-squares fitting in non-alcoholic fatty liver disease

Recent literature suggests that tri-exponential models may provide additional information and fit liver intravoxel incoherent motion (IVIM) data more accurately than conventional bi-exponential models. However, voxel-wise fitting of IVIM results in noisy and unreliable parameter maps. For bi-exponential IVIM, neural networks (NN) were able to produce superior parameter maps than conventional least-squares (LSQ) generated images. Hence, to improve parameter map quality of tri-exponential IVIM, we developed an unsupervised physics-informed deep neural network (IVIM3-NET). We assessed its performance in simulations and in patients with non-alcoholic fatty liver disease (NAFLD) and compared outcomes with bi-exponential LSQ and NN fits and tri-exponential LSQ fits. Scanning was performed using a 3.0T free-breathing multi-slice diffusion-weighted single-shot echo-planar imaging sequence with 18 b-values. Images were analysed for visual quality, comparing the bi- and tri-exponential IVIM models for LSQ fits and NN fits using parameter-map signal-to-noise ratios (SNR) and adjusted R2. IVIM parameters were compared to histological fibrosis, disease activity and steatosis grades. Parameter map quality improved with bi- and tri-exponential NN approaches, with a significant increase in average parameter-map SNR from 3.38 to 5.59 and 2.45 to 4.01 for bi- and tri-exponential LSQ and NN models respectively. In 33 out of 36 patients, the tri-exponential model exhibited higher adjusted R2 values than the bi-exponential model. Correlating IVIM data to liver histology showed that the bi- and tri-exponential NN outperformed both LSQ models for the majority of IVIM parameters (10 out of 15 significant correlations). Overall, our results support the use of a tri-exponential IVIM model in NAFLD. We show that the IVIM3-NET can be used to improve image quality compared to a tri-exponential LSQ fit and provides promising correlations with histopathology similar to the bi-exponential neural network fit, while generating potentially complementary additional parameters.

Bi-exponential fitting excluding b=0 data improves the scan-rescan stability of liver IVIM parameter measures and particularly so for the perfusion fraction

BACKGROUND

A prerequisite to translating intravoxel incoherent motion (IVIM) imaging into meaningful clinical applications is sufficient scan-rescan reproducibility. This study aims to confirm the hypothesis that IVIM data fitting by not using b=0 images will improve the stability of liver IVIM measurement.

METHODS

Healthy volunteers' liver IVIM images were prospectively acquired using a 1.5-T magnet or a 3.0 T with 16 b-values. Repeatability study subjects were scanned twice during the same session, resulted in 35 paired scans for 35 subjects (11 men, mean age: 41.82 years, range: 32-60 years; 24 women, mean age: 42.67 years, range: 20-71 years). IVIM analysis was performed with full-fitting and segmented-fitting with a threshold b-value of 60 s/mm2, and fitting started from b=0 s/mm2 or from b=2 s/mm2. Reproducibility study subjects were scanned and then rescanned with an interval of 5-18 days, resulted in 20 paired scans for 11 subjects (4 men, mean age: 26.25 years, range: 25-27 years; 7 women, mean age: 25.57 years, range: 24-27 years). IVIM analysis was performed with segmented-fitting with a threshold b-value of 50 s/mm2, and fitting started from b=0 s/mm2 or from b=3 s/mm2.

RESULTS

Fitting without b=0 data generally improved the repeatability and reproducibility for both PF and Dslow, and particularly so for PF. For with b=0 data segmented fitting repeatability, PF had within-subject standard deviation of 0.019, bland-Atman 75% agreement limit of -31.52% to 28.35%, and ICC of 0.647, while these values were 0.009, -20.78% to 16.86%, and 0.837 for without b=0 analysis. Though the repeatability and reproducibility for Dfast generally also improved, they remained suboptimal. Measurement stability was better for repeatability than for reproducibility.

CONCLUSIONS

Scan-rescan repeatability and reproducibility of liver IVIM parameters can be improved by fitting without b=0 data, which is particularly so for PF.

Usefulness of diffusion derived vessel density computed from a simplified IVIM imaging protocol: An experimental study with rat biliary duct blockage induced liver fibrosis

OBJECTIVES

Liver vessel density can be evaluated by DDVD (diffusion derived vessel density): DDVD_(b0b1) = Sb0/ROIarea0 - Sb1/ROIarea1, where Sb0 and Sb1 refer to the liver signal when b is 0 or 1 s/mm². Sb1 and ROIarea1 may be replaced by other b-values. With a rat biliary duct ligation (BDL) model, this study assesses the usefulness of liver DDVD computed from a simplified IVIM imaging protocol using b = 25 and b = 50 to replace b = 1 s/mm², alone and in combination with other IVIM parameters.

METHODS

Male Sprague-Dawley rats were used. The rat number was 5, 5, 5, and 3 respectively, for the timepoints of 7, 14, 21, 28 days post-BDL surgery. 12 rats had partial biliary duct recanalization performed after the rats had BDL for 7 days and then again followed-up for a mean of 14 days. Liver diffusion MRIs were acquired at 3.0 T with a b-value distribution of 0, 25, 50, 75, 100, 150, 300, 700, 1000 s/mm². DDVD_mean (control rats n = 6) was the mean of DDVD_(b0b25) and DDVD_(b0b50). IVIM fitting started from b = 0 s/mm² with segmented fitting and a threshold b of 50 s/mm² (n = 5 for control rats). Three 3-D spaces were constructed using a combination of the four diffusion parameters.

RESULTS

The control rats and BDL rats (n = 18) had a liver DDVD_mean of 84.0 ± 26.2 and 44.7 ± 14.4 au/pixel (p < 0.001). All 3-D spaces totally separated healthy livers and all fibrotic livers (n = 30, BDL rats and recanalization rats). The mean relative distance between healthy liver cluster and fibrotic liver cluster was 0.331 for PF, D_slow, and D_fast; 0.381 for PF, D_fast, and DDVD_mean; and 0.384 for PF, D_slow, and DDVD_mean.

CONCLUSION

A combination of PF, D_slow, and D_fast allows total separation of healthy livers and fibrotic livers and the integration of DDVD improved the separation.

Echo time dependence of biexponential and triexponential intravoxel incoherent motion parameters in the liver

PURPOSE

Intravoxel incoherent motion (IVIM) studies are performed with different acquisition protocols. Comparing them requires knowledge of echo time (TE) dependencies. The TE-dependence of the biexponential perfusion fraction f is well-documented, unlike that of its triexponential counterparts f₁ and f₂ and the biexponential and triexponential pseudodiffusion coefficients D^* , D 1 ∗ , and D 2 ∗ . The purpose was to investigate the TE-dependence of these parameters and to check whether the triexponential pseudodiffusion compartments are associated with arterial and venous blood.

METHODS

Fifteen healthy volunteers (19-58 y; mean: 24.7 y) underwent diffusion-weighted imaging of the abdomen with 24 b-values (0.2-800 s/mm² ) at TEs of 45, 60, 75, and 90 ms. Regions of interest (ROIs) were manually drawn in the liver. One set of bi- and triexponential IVIM parameters per volunteer and TE was determined. The TE-dependence was assessed with the Kruskal-Wallis test.

RESULTS

TE-dependence was observed for f (P < .001), f₁ (P = .001), and f₂ (P < .001). Their median values at the four measured TEs were: f: 0.198/0.240/0.274/0.359, f₁ : 0.113/0.139/0.146/0.205, f₂ : 0.115/0.155/0.182/0.194. D, D^* , D 1 ∗ , and D 2 ∗ showed no significant TE-dependence. Their values were: diffusion coefficient D (10^-4 mm² /s): 9.45/9.63/9.75/9.41, biexponential D^* (10^-2 mm² /s): 5.26/5.52/6.13/5.82, triexponential D 1 ∗ (10^-2 mm² /s): 1.73/2.91/2.25/2.51, triexponential D 2 ∗ (mm² /s): 0.478/1.385/0.616/0.846.

CONCLUSION

f₁ and f₂ show similar TE-dependence as f, ie, increase with rising TE; an effect that must be accounted for when comparing different studies. The diffusion and pseudodiffusion coefficients might be compared without TE correction. Because of the similar TE-dependence of f₁ and f₂ , the triexponential pseudodiffusion compartments are most probably not associated to venous and arterial blood.