4280

Time-dependent IVIM and non-Gaussian parameters in breast cancer and melanoma xenograft models; correlation with histological markers
Yuko Someya1, Mami Iima1,2, Hirohiko Imai3, Tomomi Nobashi1, Akihiko Yoshizawa4, Masako Kataoka1, Hiroyoshi Isoda1, Denis Le Bihan5,6,7, and Yuji Nakamoto1
1Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University, Kyoto, Japan, 2Institute for advancement of clinical and translational science, Kyoto University Hospital, Kyoto, Japan, 3Department of Systems science, Graduate School of Informatics, Kyoto University, Kyoto, Japan, 4Department of Diagnostic Pathology, Kyoto University, Kyoto, Japan, 5NeuroSpin/Joliot, CEA-Saclay Center, Paris-Saclay University, Gif-sur-Yvette, France, 6Human Brain Research Center, Kyoto University Graduate School of Medicine, Kyoto, Japan, 7National Institute for Physiological Sciences, Okazaki, Japan

Synopsis

Time-dependent IVIM and non-Gaussian diffusion parameters might provide additional information over conventional DWI to evaluate features linked to tissue microstructure. Time-dependent IVIM and non-Gaussian diffusion parameters were evaluated both in-vivo and post-mortem in MDA-MB-231 (breast cancer) and B16 (melanoma) models. As expected, fIVIM values significantly dropped from in-vivo to post-mortem conditions with both models at two different diffusion times, confirming the association between fIVIM and tumor perfusion, but did not reach 0, suggesting limitations in the IVIM/non-Gaussian Kurtosis model. ADCo values significantly decreased at long diffusion times and upon sacrification while Kurtosis moderately increased.

Introduction

IVIM MRI has the potential to provide quantitative information on microcirculation (perfusion) non-invasively, and the increase of fIVIM (perfusion-related fraction) has been reported as associated with the proliferation of neovessels in some cancer tumors1. Thus, IVIM parameters are expected to play an important role as imaging biomarkers to reveal and monitor tumor angiogenesis without the need for a contrast agent. Evaluating the degree of tumor angiogenesis is important for diagnosis, characterization, or staging of malignant tumors. However, fitting data to the IVIM model sometimes fails, revealing some issues, such as effects of non-Gaussian diffusion and noise floor issues2. While a decrease in fIVIM values after sacrification is expected with the association between IVIM parameters and tumor perfusion, as reported previously3, fIVIM does not always reach 0, which has been ascribed to limitations to the IVIM/Kurtosis model.
Here, we investigated IVIM and non-Gaussian diffusion parameters in-vivo and post-mortem conditions using two different tumor cell lines (breast cancer and melanoma) and evaluated the association between these parameters and pathological information, such as microvessel density and vessel size.

Materials & Methods

Seven ICR nu/nu mice and seven BALB mice were injected subcutaneously with 1×106 MDA-MB-231 cells and B16-F0-Lac cells, respectively in the hind limb. MRI was conducted on a 7T MRI scanner (Bruker, Germany) using a 1H quadrature transmit/receive volume coil. The SE-EPI acquisition parameters were set as follows; resolution 250x250μm², FOV 25x25 mm², slice-thickness 1.5 mm, TE=57ms, TR=2500ms, 8averages, 4segments. IVIM/Diffusion-weighted (DW) images were acquired using 2 different diffusion times (Td, 9 and 27.6ms) and 19 b values (7-4105 s/mm²). MRI images were obtained both in-vivo and post-mortem. The total acquisition time was 22min 40sec.
Shifted ADC (sADC) values were calculated using b values of 200 and 1500 sec/mm2 4. Non-Gaussian diffusion and fIVIM values were estimated using the combined IVIM/non-Gaussian diffusion kurtosis model corrected for Noise Floor effects5. Xenograft tumors were harvested and stained immunohistochemically for CD31. CD31-positive microvessel density (MVD) and diameters of microvessels (minor axis) were measured using QuPath (version 0.3.0, open-source software). MRI data analysis was performed using a code developed in Matlab (Mathworks, Natick, MA). Diffusion parameters were compared using the paired-t test and Wilcoxon test, and the correlation of MR parameters with histopathological parameters were investigated using the Spearman's rank correlation test.

Results

The fIVIM values significantly dropped after sacrification in B16 tumors (Mean±SD; 13.0±6.1% vs. 5.0±2.1%; P<0.05 at long-Td and 9.8±4.0% vs. 4.8±2.7%; P<0.05 at short-Td) and tended to decrease in post-mortem condition in MDA-MB-231 tumors (7.7±4.2% vs. 2.9±1.2%; P=0.058 at long-Td and 6.6±5.7% vs. 2.6±1.2%; P=0.124 at short-Td). However, fIVIM values after sacrification did not reach 0 despite the lack of perfusion, as previously reported3. fIVIM in B16 tended to be larger than MDA-MB-231 at both Td and in both in-vivo and post-mortem conditions (Fig. 1). CD31-positive MVD tended to be larger in MDA-MB-231 compared to B16 tumors (Table 1). In MDA-MB-231, numerous blood vessels with small round lumens were observed while B16 tumors had blood vessels with irregular, expanded lumens (Fig. 2). There was no significant positive association between fIVIM and MVD, but there was a correlation between fIVIM and the square of the vessel's diameters (vessel lumen area) (R=0.75, P=0.012). In both tumors, sADC values significantly decreased both when increasing the Td and after sacrification (Table 2, P<0.05, respectively). ADCo values had almost the same behavior as sADC, however, the short-Td ADCo decrease after sacrification, did not reach significance (P=0.06). K after sacrification showed a mild increase for both tumors and both Td.

Discussion

The absence of correlation in-vivo of fIVIM with MVD, while they were correlated with the vessels lumen area (more vascularization and larger capillary lumen space for B16 compared with MDA-MB-231 model6), confirms that fIVIM corresponds more to blood volume than vessel density. A significant decrease of fIVIM values after sacrification was observed both in breast cancer and melanoma xenograft models, confirming the strong association between fIVIM and tumor microperfusion, but not reaching 0. Residual post-mortem fIVIM in B16 tended to be slightly higher (5%) than MDA-MB-231 (3%) independently of Td. Those residual fIVIM effects may result from imperfect fitting of the data with the IVIM/Kurtosis model which has known limitations. However, as ADCo (and sADC) values were also higher for the B16 model we may hypothesize that the larger capillary volume might represent a pool of fast diffusing (quasi-free) water molecules in tumor capillaries, larger for B16 than MDA. This pool would artifactually contribute to the IVIM compartment, as D* in the IVIM model must be taken as the sum of the perfusion-driven pseudo-diffusion coefficient and the water diffusion coefficient in the vessel which remains present even in the absence of blood flow. This hypothesis was supported by simulations using data from this study, however, further work is in progress to investigate this hypothesis.

Conclusion

Our findings confirm that fIVIM values reflect blood microcirculation in tissues more in terms of volume than vessel density. Residual post-mortem fIVIM values may suggest the imperfections in the IVIM/Kurtosis model or the sensitivity of the model to water pools with high diffusion, close to the D* values expected by the perfusion-driven IVIM model, as, thus contaminating fIVIM values.

Acknowledgements

The authors thank the Center for Anatomical, Pathological and Forensic Medical Research, Kyoto University Graduate School of Medicine, for preparing microscope slides.

References

1. Iima, M. Perfusion-driven Intravoxel Incoherent Motion (IVIM) MRI in Oncology: Applications, Challenges, and Future Trends. Magn. Reson. Med. Sci. 20:125-138 (2021).

2. Iima, M. & Le Bihan, D. Clinical Intravoxel Incoherent Motion and Diffusion MR Imaging: Past, Present, and Future. Radiology 278, 13–32 (2016).

3. Yuko Someya, Mami Iima, Hirohiko Imai, Hiroyoshi Isoda, Masako Kataoka, Denis LeBihan, Yuji Nakamoto. Time-dependent IVIM/Non-Gaussian parameters between in vivo and post-mortem breast cancer xenograft models. ISMRM & SMRT Annual Meeting & Exhibition #1150 (2021).

4. Iima, M. et al. Effects of diffusion time on non-Gaussian diffusion and intravoxel incoherent motion (IVIM) MRI parameters in breast cancer and hepatocellular carcinoma xenograft models. Acta Radiol Open 11;7 (2018).

5. Someya, Y. et al. Investigation of Breast Cancer Microstructure and Microvasculature from Time-Dependent DWI and CEST in Correlation with Histological Biomarkers. in submission (2021).

6. Britto, D. D. et al. Macrophages enhance Vegfa-driven angiogenesis in an embryonic zebrafish tumour xenograft model. Dis. Model. Mech. 11(12) (2018).

Figures

Fig.1 Box-and-whisker plots of fIVIM values. The fIVIM value dropped post-mortem conditions in both tumors (B16 tumors revealed statistically significant decreases at both diffusion times*). After sacrification, fIVIM values did not reach 0 despite lack of perfusion. Time-dependency of fIVIM value was not observed. B16 tumors tended to have larger fIVIM values compared with MDA-MB-231 tumors, but not statistically significant.

Fig.2 Representative case of MDA-MB-231 tumor (A-C) and B16 tumor (D-F) with CD31 immunostaining. In fig. B, C, E and F, the brownish area (CD-31 positive part) was extracted and adjusted using QuPath. In MDA-MB-231 tumor (B), numerous blood vessels with small round lumens were shown. B16 tumors (E) had blood vessels with irregular, expanded lumens. MVD was 16.8 in MDA-MB-231 (A-C) and 6.3 in B16 tumors (D-F).

Fig.3 A representative case of sADC map and fIVIM map of the same B16 tumor shown in Fig. 2D-F. These maps were obtained at long diffusion time. The sADC values were clearly decreased after sacrification. fIVIM values dropped post-mortem as expected, due to the lack of perfusion, but did not reach 0. The lower sADC area (deep blue area) was correlated well with the higher fIVIM area (orange area) in vivo.

Table 1. IVIM parameters and pathological information associated with microvessel. IVIM parameters were obtained using two different diffusion times and in-vivo/postmortem conditions. * indicates a significant difference.

Table 2. Diffusion and non-Gaussian diffusion parameters obtained using two different diffusion times and in-vivo/postmortem conditions. * indicates a significant difference.

Proc. Intl. Soc. Mag. Reson. Med. 30 (2022)
4280
DOI: https://doi.org/10.58530/2022/4280