Validation of Interstitial Volume Fraction Quantification Performed with Dynamic Contrast-Enhanced Magnetic Resonance Imaging in Skeletal Swine Muscle
Stefan Hindel1, Anika Söhner1, Marc Maaß2, and Lutz Lüdemann1

1University Hospital Essen, Essen, Germany, 2Wesel Protestant Hospital, Wesel, Germany

Synopsis

We assessed the accuracy of interstitial volume fraction v(e) measurements in low-perfused tissue performed using dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) with a gadolinium-based contrast agent. A 3D gradient echo sequence with k-space-sharing was used to determine v(e) in muscle tissue of twelve pigs. The evaluation was performed with the simple and extended Tofts model and the 2-compartment exchange model using different acquisition durations (ADs). The v(e) values determined by MRI were compared with the histologic analysis of muscle tissue sections. There was good agreement between histology and DCE-MRI modeling but also a strong dependence on AD with the Tofts models.

PURPOSE

The aim of our study was to assess the accuracy of interstitial volume fraction measurement in low-perfused and low-vascularized tissue performed with dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

METHODS

The interstitial tissue volume fraction, v(e), was determined in the medial thigh muscles of twelve female pigs using a 3D gradient echo sequence with k-space-sharing and with a gadolinium-based contrast agent (gadoterate meglumine) 1. Determination was performed using three pharmacokinetic models2,3, the simple (TM) and the extended Tofts model (ETM) and the two-compartment exchange model (2CXM). We investigated the effect of varying acquisition durations (ADs) on the model parameter estimates of the three models and compared the v(e) values with the results of histologic examinations of muscle of the medial thigh muscle.

RESULTS

Histology (illustrated in Fig. 1) yielded a median (25%-75% quartile) v(e) of 4.8(3.7-6.2)%. DCE-MRI measurement (illustrated by curve fits in Fig. 2) yielded similar interstitial tissue volume fractions depending on the model and AD. The lowest v(e) determined by DCE-MRI was 5.1(3.3-6.0)% and 5.2(3.3-6.1)% for the TM and ETM at 6 min AD, respectively. The highest v(e) was 7.7(4.4-8.9)% for the TM and 7.7(4.5-9.0)% for the ETM and 15 min AD. The variation of v(e) with AD was much smaller when the 2CXM was used: v(e)=6.2 (3.1-9.2)% for 6 min AD and v(e)=6.3 (4.3-9.8)% for 15 min. Fig. 3 summarizes median values of v(e), interstitium-to-plasma constant k(ep), and volume transfer constant K(trans) obtained with the three models and at the different ADs. The best fit was found for the 2CXM at 10 min AD (v(e)=6.6(3.7-8.2)%). With increasing AD, the values of K(trans) decreased while v(e) increased for both Tofts models. Figure 4 shows the behavior of k(ep) against AD for the three models used in our experiments. For the Tofts models, k(ep) decreased asymptotically with increasing AD and converged toward a value close to k(ep) determined with the 2CXM. The interplay of increasing v(e) values and declining values for K(trans) led to a decrease in the interstitium-to-plasma rate constant, k(ep)= K(trans)/v(e).

DISCUSSION

The 2CXM uses more fitting parameters and thus yields better fits. As a result, v(e) is less dependent on AD; however, the uncertainty expressed by the 25%-75% quartile range is larger. The underfitting with the Tofts models results in an uncertainty in parameter determination, which has a major impact even for analysis of low-vascularized and low-perfused tissue, where the estimated v(e) values depend on AD. Artzi et al. already reported dependence of k(ep) (determined with the ETM) on AD in human glioblastoma4. In a simulation study5, Luypeart et al. found that, compared to the 2CXM, the ETM underestimated v(e) for short ADs and increasingly overestimated it for increasing AD, and that, conversely, K(trans) was overestimated for short ADs and increasingly underestimated for longer ADs using the ETM.

CONCLUSION

Our results reveal good agreement between histologic estimates of interstitial volume and the estimates gained by DCE-MRI modeling. Due to its good fitting accuracy and independence of AD, the 2CXM appears to be suitable for use in clinical practice. In contrast to the 2CXM, estimation of v(e), K(trans), and k(ep) using both Tofts models significantly depend on AD. Thus, use of the Tofts models is more complex as AD must be taken into account when analyzing the fitting parameters.

Acknowledgements

The authors thank the Deutsche Forschungsgemeinschaft (DFG) for supporting this research.

References

1 Sauerbrey A, Hindel S, Maaß M, et al. Establishment of a swine model for validation of perfusion measurement by dynamic contrast-enhanced magnetic resonance imaging. BioMed research international. 2014 Feb;2014:390506.

2 Tofts PS, Brix G, Buckley DL, et al. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. Journal of magnetic resonance imaging : JMRI. 1999 Sep;10(3):223–232.

3 Sourbron SP, Buckley DL. Tracer kinetic modelling in MRI: estimating perfusion and capillary permeability. Physics in medicine and biology. 2012 Jan;57(2):R1–33.

4 Artzi M, Libermann G, Nadav G, et al. The Effect of Dynamic Contrast Enhanced Acquisition Duration on Estimated Pharmacokinetic Parameters: Study of Simulated and Real Data. Proc 23rd Annual Meeting of the International Society for Magnetic Resonance in Medicine (ISMRM, Toronto, Canada, May 2015).

5 Luypaert R, Sourbron S, de Mey J. Validity of perfusion parameters obtained using the modified Tofts model: a simulation study. Magnetic resonance in medicine. 2011 May;65(5):1491–1497.

Figures

Figure 1. Frozen section stained with modified van-Gieson staining. The intercellular space between the yellowish muscle cells shows connective tissue appearing in magenta and the white intercellular fluid space (A). The interstitial fluid space was labeled semi-automatically using a morphometry software to determine its area (B). The connective tissue fibers appear in magenta and were labeled in the same way (C).

Figure 2. Results of the curve fits of two exemplary experiments for an acquisition period of approx. 10 min, for the TM (blue), ETM (green) and 2CXM (red). The 2CXM yields a good fit for both the inflow region and the wash-out region and therefore is reasonably independent of AD. In contrast, the fit quality of both Tofts models in these two regions strongly depends on AD, resulting in a decrease of k(ep) towards longer ADs.

Figure 3. The table presents median values of interstitial volume, volume-transfer constant, and interstitium-to-plasma constant. The results are presented for the three models and the different acquisition durations. Strikingly, K(trans) and k(ep) drop considerably with increasing AD for the Tofts models and v(e) rises, whereas all quantities remain relatively constant for the 2CXM.

Figure 4. The median interstitium-to-plasma rate constant k(ep) obtained with the three different models is plotted against the duration of data acquisition. A convergent behavior of k(ep) can be assumed for large ADs when employing the Tofts models. Thereby, the TM and ETM convergence values seem to be in the range of the 2CXM k(ep) already determined for small ADs.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2460