2122

Echo time dependence of the biexponential and triexponential IVIM parameters in the liver
Tobit Führes1, Andreas Julian Riexinger1, Jan Martin2, Martin Gerhard Loh1, Bernhard Hensel3, Michael Uder1, and Frederik Bernd Laun1
1University Hospital Erlangen, Erlangen, Germany, 2Lund University, Lund, Sweden, 3Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany

Synopsis

In this study, we acquired diffusion-weighted image data of 15 healthy volunteers at four different TE. A ROI-based IVIM analysis showed that the pseudodiffusion coefficients did not change with TE, whereas the perfusion fractions rose significantly for bi- and triexponential IVIM curves at longer TE. A quantitative analysis indicates that the triexponential perfusion fractions are unlikely attributable to venous and arterial blood compartments. Our results indicate moreover that the pseudodiffusion coefficients can be compared without TE-correction among studies performed with different TE.

Introduction

The intravoxel incoherent motion (IVIM) effect allows one to distinguish between diffusion and perfusion by considering the diffusion-weighted signal at various b-values. Recent studies showed that a triexponential signal curve describes the data better in several organs than the classically used biexponential signal curve1-4. In our study, the TE dependence of the IVIM parameters in the liver was measured with a twofold purpose. First, to investigate whether the well-documented TE-dependence of the biexponential perfusion fraction f also occurs for the biexponential pseudodiffusion coefficient and for the triexponential IVIM parameters; and would thus require a correction strategy similar to the f-correction approach by Jerome et al.5 Second, to elucidate whether the triexponential perfusion fractions f1 and f2 might be linked to arterial and venous compartments. This might be the case because arterial and venous blood exhibit different T2 times.

Methods

At 3 T, diffusion-weighted images of the liver of 15 healthy volunteers were acquired at TE = 45 ms, 60 ms, 75 ms, 90 ms.6 For 24 b-values, which were spaced between 0.2 s/mm2 and 800 s/mm2, data for six diffusion directions and of six transversal slices were acquired. The diffusion time and the diffusion gradient length were kept constant at 21640 µs and 12240 µs, respectively, for all diffusion weightings and echo times. The b-values used in the evaluation were calculated numerically in order to take additional diffusion weightings generated e.g. by imaging gradients into account.
The evaluation of the images was done with MATLAB R2017b. For each slice, a single region of interest (ROI) under exclusion of large vessels was drawn and evaluated using the median of all included voxel signals. The signals were normalized to the signal at the lowest b-value of the respective echo time. The triexponential IVIM curve was fitted to all data points; data from each direction and slice was treated as equally-weighted data points. The biexponential IVIM curve was fitted to all b-values except for those in the interval [0.3 s/mm2, 6 s/mm2] to ensure consistency with previous biexponential IVIM studies. The fits were performed in a two-step multistart approach. The formula for the biexponential IVIM signal curve reads:
$$S(b) = S_0 ((1-f) \exp(-bD) + f \exp(-bD^*))$$ where f denotes the biexponential perfusion fraction, D the diffusion coefficient and D* the pseudodiffusion coefficient. The formula for the triexponential IVIM signal curve reads:
$$S(b) = S_0 ((1-f_1-f_2) \exp(-bD) + f_1 \exp(-bD_1^*) + f_2 \exp(-bD_2^*))$$ where f1 and f2 denote the triexponential perfusion fractions and D1* and D2* denote the triexponential pseudodiffusion coefficients. Statistical analysis was performed with the Kruskal-Wallis test and the Dunn-Sidak test as post-hoc test.

Results

Representative data and fits are shown in Figure 1. The resulting parameters are shown in boxplots in Figure 2. The median diffusion coefficient D was in a range between 9.45 x 10-4 mm2/s and 9.75 x 10-4 mm2/s and not echo-time dependent (p = 0.848).
Biexponential parameters: The median biexponential perfusion fraction f was in the range [19.76 %, 35.94 %] and was significantly dependent on the echo time (p < 0.001). Significant changes occured between (45 ms, 75 ms), (45 ms, 90 ms), and (60 ms, 90 ms). The median biexponential pseudodiffusion coefficient D* assumed values in a range of [5.26 x 10-2 mm2/s, 6.13 x 10-2 mm2/s] and was not echo time dependent (p = 0.360).
Triexponential parameters: The median perfusion fractions f1 and f2 were in the range [11.25 %, 20.49 %] and [11.48 %, 19.38 %], respectively. Both were echo time dependent (p = 0.001 and p < 0.001, respectively). For f1, significant changes occured between the echo time pairs (45 ms, 90 ms) and (60 ms, 90 ms), for f2 between (45 ms, 75 ms) and (45 ms, 90 ms). The median pseudodiffusion coefficients D1* and D2* were in the range [1.73 x 10-2 mm2/s, 2.91 x 10-2 mm2/s] and [0.478 mm2/s, 1.385 mm2/s], respectively. Both of them did not show a significant echo time dependence (p = 0.360 and p = 0.053, respectively).

Discussion

The values of the biexponential IVIM parameters are well in line with other studies.7 The fact that the biexponential perfusion fraction rises with the echo time was also described in previous studies and can be explained by the different T2 times of blood and tissue.8 The behavior of the triexponential perfusion fractions is in accordance with the observations made for biexponential IVIM.8 The observation that the pseudodiffusion coefficients are not significantly echo time dependent is consistent with previous reports.9

Conclusion

The almost identical monotonic rise of the triexponential perfusion fractions f1 and f2 with echo time indicates that their T2 times are similar and that corresponding compartments may not be linked straightforwardly to arterial and venous blood compartment, which is in agreement with the study by Riexinger et al.1 The independence of bi- and triexponential pseudodiffusion coefficients indicates that they can be compared without TE-correction among studies.

Acknowledgements

Funding by the Deutsche Forschungsgemeinschaft is gratefully acknowledged (LA 2804/12-1, LA 2804/13-1). We thank the Imaging Science Institute (Erlangen, Germany) for providing us with measurement time.

References

1. Riexinger AJ, Martin J, Rauh S, et al. On the Field Strength Dependence of Bi- and Triexponential Intravoxel Incoherent Motion (IVIM) Parameters in the Liver. Journal of magnetic resonance imaging : JMRI. 2019.

2. van Baalen S, Leemans A, Dik P, Lilien MR, Ten Haken B, Froeling M. Intravoxel incoherent motion modeling in the kidneys: Comparison of mono-, bi-, and triexponential fit. Journal of magnetic resonance imaging : JMRI. 2017;46(1):228-239.

3. Cercueil JP, Petit JM, Nougaret S, et al. Intravoxel incoherent motion diffusion-weighted imaging in the liver: comparison of mono-, bi- and tri-exponential modelling at 3.0-T. Eur Radiol. 2015;25(6):1541-1550.

4. Wurnig MC, Germann M, Boss A. Is there evidence for more than two diffusion components in abdominal organs? - A magnetic resonance imaging study in healthy volunteers. NMR in biomedicine. 2017.

5. Jerome NP, d'Arcy JA, Feiweier T, et al. Extended T2-IVIM model for correction of TE dependence of pseudo-diffusion volume fraction in clinical diffusion-weighted magnetic resonance imaging. Physics in medicine and biology. 2016;61(24):N667-N680.

6. Hanspach J, Nagel AM, Hensel B, Uder M, Koros L, Laun FB. Sample size estimation: Current practice and considerations for original investigations in MRI technical development studies. Magnetic resonance in medicine. 2020.

7. Li YT, Cercueil JP, Yuan J, Chen W, Loffroy R, Wang YX. Liver intravoxel incoherent motion (IVIM) magnetic resonance imaging: a comprehensive review of published data on normal values and applications for fibrosis and tumor evaluation. Quantitative imaging in medicine and surgery. 2017;7(1):59-78.

8. Lemke A, Laun FB, Simon D, Stieltjes B, Schad LR. An in vivo verification of the intravoxel incoherent motion effect in diffusion-weighted imaging of the abdomen. Magnetic resonance in medicine. 2010;64(6):1580-1585.

9. Barbieri S, Donati OF, Froehlich JM, Thoeny HC. Comparison of Intravoxel Incoherent Motion Parameters across MR Imagers and Field Strengths: Evaluation in Upper Abdominal Organs. Radiology. 2016;279(3):784-794.

Figures

Figure 1: Signal curves. For better readability, the displayed signal values represent the values averaged over diffusion directions and slices . (a) Representative normalized signals of one volunteer at TE = 45 ms. Both the biexponential and triexponential fit curves are shown. Data points with lower brightness were used for the triexponential fit only. (b) Representative absolute signal of one volunteer at four echo times. Lines represent the triexponential IVIM curves.

Figure 2: Box plots of the bi- and triexponential IVIM parameters. Each point indicates the data of one volunteer. The central mark of each box plot indicates the respective median, the bottom and top lines indicate the 25th and 75th percentile. The whiskers reach out to ±2.7 $$$\sigma$$$, outliers are shown as “+”. Dashed horizontal lines show the threshold where the scale of the y-axis is changed in order to display all points in a reasonable way.

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
2122