On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver

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.

Effect of simultaneous multislice imaging, slice properties, and repetition time on the measured magnetic resonance biexponential intravoxel incoherent motion in the liver

OBJECTIVES

This study aims to investigate the previously reported dependency of intravoxel incoherent motion (IVIM) parameters on simultaneous multislice (SMS) acquisition and repetition time (TR). This includes the influence of slice thickness, slice gaps, and slice order on measured IVIM parameters.

MATERIALS AND METHODS

Diffusion-weighted imaging (DWI) of the liver was performed on 10 healthy volunteers (aged 20-30 years) at 3T with a slice thickness of 5 mm, a slice gap of 5 mm, and a linear slice order. Diffusion-weighted images were acquired with 19 b-values (0-800 s/mm2) using both conventional slice excitation with an acceleration factor of one (AF1) and SMS excitation with an acceleration factor of three (AF3). Each of these measurements were carried out with two repetition times (TRs)- 1,300 ms (prefix s) and 4,500 ms (prefix l)-resulting in four different combinations: sAF1, sAF3, lAF1, and lAF3. Five volunteers underwent additional measurements using a 10 mm slice thickness and with AF1. Median signal values in the liver were used to determine the biexponential IVIM parameters. Statistical significances were assessed using the Kruskal-Wallis test, Wilcoxon signed-rank test, and Student's t-test. In-silico investigations were also used to interpret the data.

RESULTS

There were no significant differences between the biexponential IVIM parameters acquired from sAF1, sAF3, lAF1, and lAF3. Median values of the perfusion fraction f were as follows: 29.9% (sAF1), 26.9% (sAF3), 28.1% (lAF1), and 27.5% (lAF3). In the 10 mm-thick slices, f decreased from 31.3% (lAF1) to 27.4% (sAF1) (p = 0.141).

CONCLUSION

The slice excitation mode did not appear to have any significant influence on the biexponential IVIM parameters. However, our simulations, as well as values reported from previous publications, show that slice thickness, slice gaps, and slice order are relevant and should thus be reported in IVIM studies.

Influence of Magnetic Field Strength on Intravoxel Incoherent Motion Parameters in Diffusion MRI of the Calf

Background: The purpose of this study was to investigate the dependence of Intravoxel Incoherent Motion (IVIM) parameters measured in the human calf on B0. Methods: Diffusion-weighted image data of eight healthy volunteers were acquired using five b-values (0-600 s/mm2) at rest and after muscle activation at 0.55 and 7 T. The musculus gastrocnemius mediale (GM, activated) was assessed. The perfusion fraction f and diffusion coefficient D were determined using segmented fits. The dependence on field strength was assessed using Student's t-test for paired samples and the Wilcoxon signed-rank test. A biophysical model built on the three non-exchanging compartments of muscle, venous blood, and arterial blood was used to interpret the data using literature relaxation times. Results: The measured perfusion fraction of the GM was significantly lower at 7 T, both for the baseline measurement and after muscle activation. For 0.55 and 7 T, the mean f values were 7.59% and 3.63% at rest, and 14.03% and 6.92% after activation, respectively. The biophysical model estimations for the mean proton-density-weighted perfusion fraction were 3.37% and 6.50% for the non-activated and activated states, respectively. Conclusions: B0 may have a significant effect on the measured IVIM parameters. The blood relaxation times suggest that 7 T IVIM may be arterial-weighted whereas 0.55 T IVIM may exhibit an approximately equal weighting of arterial and venous blood.

Image denoising and model-independent parameterization for IVIM MRI

Objective. To improve intravoxel incoherent motion imaging (IVIM) magnetic resonance Imaging quality using a new image denoising technique and model-independent parameterization of the signal versusb-value curve.Approach. IVIM images were acquired for 13 head-and-neck patients prior to radiotherapy. Post-radiotherapy scans were also acquired for five of these patients. Images were denoised prior to parameter fitting using neural blind deconvolution, a method of solving the ill-posed mathematical problem of blind deconvolution using neural networks. The signal decay curve was then quantified in terms of several area under the curve (AUC) parameters. Improvements in image quality were assessed using blind image quality metrics, total variation (TV), and the correlations between parameter changes in parotid glands with radiotherapy dose levels. The validity of blur kernel predictions was assessed by the testing the method's ability to recover artificial 'pseudokernels'. AUC parameters were compared with monoexponential, biexponential, and triexponential model parameters in terms of their correlations with dose, contrast-to-noise (CNR) around parotid glands, and relative importance via principal component analysis.Main results. Image denoising improved blind image quality metrics, smoothed the signal versusb-value curve, and strengthened correlations between IVIM parameters and dose levels. Image TV was reduced and parameter CNRs generally increased following denoising.AUCparameters were more correlated with dose and had higher relative importance than exponential model parameters.Significance. IVIM parameters have high variability in the literature and perfusion-related parameters are difficult to interpret. Describing the signal versusb-value curve with model-independent parameters like theAUCand preprocessing images with denoising techniques could potentially benefit IVIM image parameterization in terms of reproducibility and functional utility.

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.

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.

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.

b value and first-order motion moment optimized data acquisition for repeatable quantitative intravoxel incoherent motion DWI

PURPOSE

To design a b value and first-order motion moment (M1 ) optimized data acquisition for repeatable intravoxel incoherent motion (IVIM) quantification in the liver.

METHODS

Cramer-Rao lower bound optimization was performed to determine optimal monopolar and optimal 2D samplings of the b-M1 space based on noise performance. Monte Carlo simulations were used to evaluate the bias and variability in estimates obtained using the proposed optimal samplings and conventional monopolar sampling. Diffusion MRI of the liver was performed in 10 volunteers using 3 IVIM acquisitions: conventional monopolar, optimized monopolar, and b-M1 -optimized gradient waveforms (designed based on the optimal 2D sampling). IVIM parameter maps of diffusion coefficient, perfusion fraction, and blood velocity SD were obtained using nonlinear least squares fitting. Noise performance (SDs), stability (outlier percentage), and test-retest or scan-rescan repeatability (intraclass correlation coefficients) were evaluated and compared across acquisitions.

RESULTS

Cramer-Rao lower bound and Monte Carlo simulations demonstrated improved noise performance of the optimal 2D sampling in comparison to monopolar samplings. Evaluating the designed b-M1 -optimized waveforms in healthy volunteers, significant decreases (p < 0.05) in the SDs and outlier percentages were observed for measurements of diffusion coefficient, perfusion fraction, and blood velocity SD in comparison to measurements obtained using monopolar samplings. Good-to-excellent repeatability (intraclass correlation coefficients ≥ 0.77) was observed for all 3 parameters in both the right and left liver lobes using the b-M1 -optimized waveforms.

CONCLUSIONS

2D b-M1 -optimized data acquisition enables repeatable IVIM quantification with improved noise performance. 2D acquisitions may advance the establishment of IVIM quantitative biomarkers for liver diseases.

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.

Correlation Between IVIM-DWI Parameters and Pathological Classification of Idiopathic Orbital Inflammatory Pseudotumors: A Preliminary Study

Objective

To investigate the correlation between intravoxel incoherent motion diffusion-weighted imaging (IVIM-DWI) and the pathological classification of idiopathic orbital inflammatory pseudotumors (IOIPs).

Methods

Nineteen patients who were diagnosed with IOIPs (a total of 24 affected eyes) between November 2018 and December 2020 were included in the study. All the patients underwent magnetic resonance imaging orbital plain scans and IVIM-DWI multiparameter scans before an operation. The true diffusion coefficient (D), pseudodiffusion coefficient (D*), and perfusion fraction (f) values were obtained. Based on histopathology, the lesions were divided into three types: lymphocytic infiltration, fibrosclerotic, and mixed. The correlation between IVIM-DWI parameters and pathological classification was tested with the histopathological results as the gold standard. The data were analyzed using SPSS version 17.0, with P < 0.05 defined as significant.

Results

Among the 19 patients (24 eyes) affected by IOIP, there were no significant differences between IOIP pathological classification and gender or age (P > 0.05). There were statistically significant differences between the D and f values for different pathological types of IOIP and IVIM parameters (P < 0.05), and there was no significant difference in D* value between the different pathological types (P > 0.05).

Conclusion

The D and f values showed correlation with different types of IOIP, and the sensitivity of the D value was higher than that of the f value. The D* value showed no significant distinction between pathological types of IOIP.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Image Quality and Detection of Small Focal Liver Lesions in Diffusion-Weighted Imaging: Comparison of Navigator Tracking and Free-Breathing Acquisition

OBJECTIVES

The aim of this study was to compare intraindividual diffusion-weighted imaging (DWI) of the liver acquired with free breathing (FB) versus navigator triggering (NT) for assessing small focal liver lesions (FLLs) in noncirrhotic patients.

MATERIALS AND METHODS

Patients with known or suspected multiple FLLs were prospectively included, and spin-echo echo-planar DWI with NT and FB acquisition was performed (b-values, 50 and 800 s/mm2 [b50 and b800]). NT and FB DWI sequences with similar acquisitions times were used. Liver and lesion signal-to-noise ratios were measured at b800. The DWI scans were analyzed independently by 2 readers. Liver edge delineation, presence of stair-step artifacts, vessel sharpness, severity of cardiac motion artifacts, overall image quality, and lesion conspicuity were rated with 5-point Likert scales. Small and large FLLs (ie, <1 cm or ≥1 cm) were rated separately for lesion conspicuity. The FLL detectability was estimated by comparing the number of lesions visible with FB to those visible with NT.

RESULTS

Forty-three patients were included in the study. The FB acquisition performed better in terms of severity of cardiac motion artifacts. The NT performed better in terms of liver edge delineation and vessel sharpness. Little difference was found for stair-step artifact, overall image quality, and conspicuity of large FLL, whereas the conspicuity of small FLL was better for NT. For small FLL, both readers found more lesions with NT in 11 cases at b800. For large FLL, this effect was much less pronounced (1 case at b800 reported by 1 of the readers). The mean liver and lesion signal-to-noise ratios were 16.8/41.5 and 19.8/38.4 for NT/FB, respectively.

CONCLUSIONS

Small FLL detection is better with NT. Large FLL detection by FB and NT is similarly good. We conclude that NT should be used.

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.

Characterisation of microvessel blood velocity and segment length in the brain using multi-diffusion-time diffusion-weighted MRI

Multi-diffusion-time diffusion-weighted MRI can probe tissue microstructure, but the method has not been widely applied to the microvasculature. At long diffusion-times, blood flow in capillaries is in the diffusive regime, and signal attenuation is dependent on blood velocity (v) and capillary segment length (l). It is described by the pseudo-diffusion coefficient (D*=vl/6) of intravoxel incoherent motion (IVIM). At shorter diffusion-times, blood flow is in the ballistic regime, and signal attenuation depends on v, and not l. In theory, l could be estimated using D* and v. In this study, we compare the accuracy and repeatability of three approaches to estimating v, and therefore l: the IVIM ballistic model, the velocity autocorrelation model, and the ballistic approximation to the velocity autocorrelation model. Twenty-nine rat datasets from two strains were acquired at 7 T, with b-values between 0 and 1000 smm^-2 and diffusion times between 11.6 and 50 ms. Five rats were scanned twice to assess scan-rescan repeatability. Measurements of l were validated using corrosion casting and micro-CT imaging. The ballistic approximation of the velocity autocorrelation model had lowest bias relative to corrosion cast estimates of l, and had highest repeatability.

Age and gender dependence of liver diffusion parameters and the possibility that intravoxel incoherent motion modeling of the perfusion component is constrained by the diffusion component

The aim of this study was to establish reference values for middle-aged subjects and to investigate the age and gender dependence of liver diffusion MRI parameters. The IVIM type of liver diffusion scan was based on a single-shot spin-echo-type echo-planar sequence using a 1.5 T magnet with 16 b-values. Diffusion-derived vessel density (DDVD)(b0b2) or DDVD(b0b10) was the signal difference between b = 0 and b = 2 (or b = 10) s/mm² images after removing visible vessels. IVIM analysis was performed with full fitting and segmented fitting, and with a threshold b-value of 60 or 200 s/mm² , and fitting started from b = 2 s/mm² . Thirty-one men (age range: 25-71 years) and 26 men (age: 22-69 years) had DDVD and IVIM analysis, respectively, while 37 women (age: 20-71 years) and 36 women (age: 20-71 years) had DDVD and IVIM analysis, respectively. DDVD results showed a significant age-related reduction for women. IVIM results for full fitting showed excellent agreement with those for segmented fitting using a threshold b of 60 s/mm² , but this was less good for results with a threshold b of 200 s/mm² . As age increased, female subjects' D_slow measure showed a significant reduction, while their PF and D_fast measures showed a significant increase. For the age group of 40-55 years, DDVD(b0b2), DDVD(b0b10), D_slow , PF and D_fast were 12.26 ± 3.90 au/pixel, 16.95 ± 5.45 au/pixel, 1.072 ± 0.067 (10^-3 mm² /s), 0.141 ± 0.025 and 61.0 ± 14.0 (10^-3 mm² /s) for men, and 13.35 ± 3.6 au/pixel, 17.20 ± 3.62 au/pixel, 1.069 ± 0.074 (10^-3 mm² /s), 0.119 ± 0.014 and 57.1 ± 13.2 (10^-3 mm² /s) for women, respectively. DDVD measure of this study suggest that aging is associated with a reduction in liver perfusion. There is a possibility that a lower D_slow measure is associated with artificially higher PF and D_fast measures, and that IVIM modeling of the perfusion component is constrained by the diffusion component.

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

Diffusion weighted imaging (DWI) and intravoxel incoherent motion (IVIM) have been explored to assess liver tumors and diffused liver diseases. IVIM reflects the microscopic translational motions that occur in voxels in magnetic resonance (MR) DWI. In biologic tissues, molecular diffusion of water and microcirculation of blood in the capillary network can be assessed using IVIM DWI. The most commonly applied model to describe the DWI signal is a bi-exponential model, with a slow compartment of diffusion linked to pure molecular diffusion (represented by the coefficient D_slow), and a fast compartment of diffusion, related to microperfusion (represented by the coefficient D_fast). However, high variance in D_fast estimates has been consistently shown in literature for liver IVIM, restricting its application in clinical practice. This variation could be explained by the presence of another very fast compartment of diffusion in the liver. Therefore, a tri-exponential model would be more suitable to describe the DWI signal. This article reviews the published evidence of the existence of this additional very fast diffusion compartment and discusses the performance and limitations of the tri-exponential model for liver IVIM in current clinical settings.

An optimized b-value distribution for triexponential intravoxel incoherent motion (IVIM) in the liver

PURPOSE

To find an optimized b-value distribution for reproducible triexponential intravoxel incoherent motion (IVIM) exams in the liver.

METHODS

A numeric optimization of b-value distributions was performed using the triexponential IVIM equation and 27 different IVIM parameter sets. Starting with an initially optimized distribution of 6 b-values, the number of b-values was increased stepwise. Each new b-value was chosen from a set of 64 predefined b-values based on the computed summed relative mean error of the fitted triexponential IVIM parameters. This process was repeated for up to 100 b-values. In simulations and in vivo measurements, optimized b-value distributions were compared to 4 representative distributions found in literature.

RESULTS

The first 16 optimized b-values were 0, 0.3, 0.3, 70, 200, 800, 70, 1, 3.5, 5, 70, 1.2, 6, 45, 1.5, and 60 in units of s/mm² . Low b-values were much more frequent than high b-values. The optimized b-value distribution resulted in a higher fit stability compared to distributions used in literature in both, simulation and in vivo measurements. Using more than 6 b-values, ideally 16 or more, increased the fit stability considerably.

CONCLUSION

Using optimized b-values, the fit uncertainty in triexponential IVIM can be largely reduced. Ideally, 16 or more b-values should be acquired.

Value of intravoxel incoherent motion in detecting and staging liver fibrosis: A meta-analysis

BACKGROUND

Liver fibrosis (LF) is a common pathological feature of all chronic liver diseases. With the accumulation of extracellular matrix in the fibrotic liver, true molecular water diffusion and perfusion-related diffusion are restricted. Intravoxel incoherent motion (IVIM) can capture the information on tissue diffusivity and microcapillary perfusion separately and reflect the fibrotic severity with diffusion coefficients.

AIM

To investigate the diagnostic performance of IVIM in detecting and staging LF with histology as a reference standard.

METHODS

A comprehensive literature search was conducted to identify studies on the diagnostic accuracy of IVIM for assessment of histologically proven LF. The stages of LF were classified as F0 (no fibrosis), F1 (portal fibrosis without septa), F2 (periportal fibrosis with few septa), F3 (septal fibrosis), and F4 (cirrhosis) according to histopathological findings. Data were extracted to calculate the pooled sensitivity, specificity, positive and negative likelihood ratios, and diagnostic odds ratio, as well as the area under the summary receiver operating characteristic curve (AUC) in each group.

RESULTS

A total of 12 studies with 923 subjects were included in this meta-analysis with 5 studies (n = 465) for LF ≥ F1, 9 studies (n = 757) for LF ≥ F2, 4 studies (n = 413) for LF ≥ F3, and 6 studies (n = 562) for LF = F4. The pooled sensitivity and specificity were estimated to be 0.78 (95% confidence interval: 0.73-0.82) and 0.81 (0.74-0.86) for LF ≥ F1 detection with IVIM; 0.82 (0.79-0.86) and 0.80 (0.75-0.84) for staging F2 fibrosis; 0.85 (0.79-0.90) and 0.83 (0.77-0.87) for staging F3 fibrosis, and 0.90 (0.84-0.94) and 0.75 (0.70-0.79) for detecting F4 cirrhosis, respectively. The AUCs for LF ≥ F1, F2, F3, F4 detection were 0.862 (0.811-0.914), 0.883 (0.856-0.909), 0.886 (0.865-0.907), and 0.899 (0.866-0.932), respectively. Moderate to substantial heterogeneity was observed with inconsistency index (I²) ranging from 0% to 77.9%. No publication bias was detected.

CONCLUSION

IVIM is a noninvasive tool with good diagnostic performance in detecting and staging LF. Optimized and standardized IVIM protocols are needed to further improve its diagnostic accuracy in clinical practice.

Optimal acquisition scheme for flow-compensated intravoxel incoherent motion diffusion-weighted imaging in the abdomen: An accurate and precise clinically feasible protocol

PURPOSE

Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time.

METHODS

Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys.

RESULTS

Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively.

CONCLUSIONS

We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min).

A mixed waveform protocol for reduction of the cardiac motion artifact in black-blood diffusion-weighted imaging of the liver

OBJECTIVE

Diffusion-weighted imaging (DWI) in the liver suffers from signal loss due to the cardiac motion artifact, especially in the left liver lobe. The purpose of this work was to improve the image quality of liver DWI in terms of cardiac motion artifact reduction and achievement of black-blood images in low b-value images.

MATERIAL AND METHODS

Ten healthy volunteers (age 20-31 years) underwent MRI examinations at 1.5 T with a prototype DWI sequence provided by the vendor. Two diffusion encodings (i.e. waveforms), monopolar and flow-compensated, and the b-values 0, 20, 50, 100, 150, 600 and 800 s/mm² were used. Two Likert scales describing the severity of the pulsation artifact and the quality of the black-blood state were defined and evaluated by two experienced radiologists. Regions of interest (ROIs) were manually drawn in the right and left liver lobe in each slice and combined to a volume of interest (VOI). The mean and coefficient of variation were calculated for each normalized VOI-averaged signal to assess the severity of the cardiac motion artifact. The ADC was calculated using two b-values once for the monopolar data and once with mixed data, using the monopolar data for the small and the flow-compensated data for the high b-value. A Wilcoxon rank sum test was used to compare the Likert scores obtained for monopolar and flow-compensated data.

RESULTS

At b-values from 20 to 150 s/mm², unlike the flow-compensated diffusion encoding, the monopolar encoding yielded black blood in all images with a negligible signal loss due to the cardiac motion artifact. At the b-values 600 and 800 s/mm², the flow-compensated encoding resulted in a significantly reduced cardiac motion artifact, especially in the left liver lobe, and in a black-blood state. The ADC calculated with monopolar data was significantly higher in the left than in the right liver lobe.

CONCLUSION

It is recommendable to use the following mixed waveform protocol: Monopolar diffusion encodings at small b-values and flow-compensated diffusion encodings at high b-values.