Assessment of tumor perfusion, oxygenation, and metabolism using DCE, BOLD, and hyperpolarized 13C MRI in a mouse model of breast cancer
Erin B Adamson1, Roberta M Strigel1,2,3, David J Niles1, Kai D Ludwig1, Ben L Cox1,4,5, Amy R Moser2,6, and Sean B Fain1,3,7

1Medical Physics, University of Wisconsin-Madison, Madison, WI, United States, 2Carbone Cancer Center, University of Wisconsin-Madison, Madison, WI, United States, 3Radiology, University of Wisconsin-Madison, Madison, WI, United States, 4Morgridge Institute for Research, Madison, WI, United States, 5Laboratory for Optical and Computational Instrumentation, University of Wisconsin-Madison, Madison, WI, United States, 6Human Oncology, University of Wisconsin-Madison, Madison, WI, United States, 7Biomedical Engineering, University of Wisconsin-Madison, Madison, WI, United States

Synopsis

Hyperpolarized (HP) 13C MRSI, dynamic contrast-enhanced (DCE) MRI, and blood-oxygen-level dependent (BOLD) MRI have the potential to non-invasively characterize tumor metabolism, perfusion, and oxygenation, respectively, and aid in the development of individualized treatment plans for cancer patients. However, a regional comparison of these non-invasive techniques for probing the tumor microenvironment has not been explored. This work aims to test the feasibility of performing quantitative, spatial analysis and comparison of HP 13C MRSI and BOLD and DCE MRI in a murine breast cancer model.

Purpose

Developing individualized treatment plans for breast cancer patients necessitates improved characterization of tumor microenvironment to predict tumor progression and response to treatment. Factors implicated in tumor progression include perfusion, oxygenation, and metabolism. Due to the poorly structured vasculature associated with oncologic angiogenesis, it is hypothesized that dynamic contrast-enhanced (DCE) and blood-oxygen-level dependent (BOLD) MRI can be used to characterize associations between regional vascular supply and oxygen availability in breast cancer. DCE parameters have been successfully correlated with exogenous hypoxia markers in cancerous tumors1 and it is suggested that BOLD parameters provide additional tissue oxygenation information not solely controlled by perfusion2. Additionally, aerobic glycolysis (Warburg Effect) is postulated to be upregulated in cancerous tumors3 and has the potential for in vivo interrogation using hyperpolarized (HP) 13C magnetic resonance spectroscopic imaging (MRSI) of [1-13C]pyruvate and its downstream metabolites4. In this study, murine mammary adenocarcinomas were imaged using DCE, BOLD, and HP 13C MRI to assess their viability to spatially and quantitatively evaluate tumor vasculature, oxygenation, and metabolism.

Methods

A pilot study complying with institutional animal care and use committee regulations was conducted using a murine breast cancer model. Cells from ethylnitrosourea-induced mammary adenocarcinomas developed in FVB.B6-ApcMin/+ mice were injected into axillary fat pads of syngeneic FVB/Tac mice. Mice with both moderately aggressive (fast-growing) and markedly aggressive (very fast-growing) tumor lines were imaged. Imaging was performed on a 4.7T small animal scanner (Agilent, Palo Alto, CA) with a 1H/13C dual-tuned volume coil and 13C surface coil (Doty Scientific, Columbia, SC). A high-resolution (0.25×0.25×2mm3) T2-weighted FSE sequence was acquired for anatomical reference (TR/TEĀ­eff= 3500/66ms) followed by a multi-echo SPGR T2*-weighted BOLD sequence (TR/TE/ΔTE= 350/2/2.8ms, 32 echoes) with identical resolution. [1-13C]pyruvate was polarized (10-20%) via dynamic nuclear polarization (HyperSense, Oxford Instruments, UK) and 10μL/g was injected into the tail vein for imaging. A single-shot spiral acquisition (3×3×5mm3, TR/TE/ΔTE=90/1.00/1.19ms, echoes=5, FA=10°) acquired dynamic 13C images interleaved with slice-selective spectra (FA=5°) at ~5s temporal resolution. Spiral image reconstruction used an iterative, least-squares estimation technique5. Actual flip-angle imaging (AFI) and T1 maps were generated from 3D SPGR sequences (TR/TE= 6.1/1.2ms and 5.9/1.7ms, respectively) with 0.5×0.5×0.5mm3 resolution. Finally, a T1-weighted DCE SPGR sequence was acquired (TR/TE=20.5/2.9ms, 0.25×0.25×2mm3) with ~5s temporal resolution. 10 frames were acquired prior to IV injection of 0.15mmol/kg of gadodiamide followed by a 10min acquisition. R2* maps were generated by fitting BOLD data to a linearized exponential decay model. Ratio maps of the volume transfer constant (Ktrans) to extravascular-extracellular volume fraction (ve) were developed by converting DCE signal to gadodiamide concentration using AFI-corrected T1-maps6, then fitting to a linearized reference region model7,8. Voxel-wise ratio maps of the area-under-the-curve (AUC) of HP lactate-to-pyruvate signal were also generated. R2* and Ktrans/ve were corrected for major outliers with removal of values greater than three times the inter-quartile range. Mean R2* and Ktrans/ve values were calculated for the total volume of tumor present in each slice.

Results

Three mice each hosting markedly and moderately aggressive tumor lines underwent HP 13C, DCE, and BOLD MRI. Ktrans/ve and R2* maps displayed heterogeneity in both tumor strains, with a tendency of Ktrans/ve to be larger towards tumor peripheries and R2* to be largest towards the tumor core or air-tissue interfaces secondary to susceptibility (Figure 1). Mean R2* and Ktrans/ve values were similar for both tumor lines (Figure 2). The 13C AUC ratio map (Figure 3) displayed heterogeneous lactate/pyruvate signal within the tumors.

Discussion

R2* indicates tissue deoxyhemoglobin concentration, with higher R2* values indicating relative hypoxia, while Ktrans/ve indicates perfusion/permeability of the vasculature. Larger Ktrans/ve values were observed around the tumor periphery, suggesting more leaky, perfused vasculature in this region compared with the core. High R2* values near the tumor core may indicate the local oxygenated hemoglobin content is low, likely due to a combination of rapid proliferation inhibiting oxygen diffusion and necrosis. Follow-up histology will reveal the degree of necrosis and predominance of tissue oxygen stress markers. 13C AUC maps demonstrated higher pyruvate and lactate values along the medial tumor, possibly due to more HP substrate delivery via native vasculature from the abdominal region. The 13C AUC ratio map demonstrates heterogeneous distribution of lactate-to-pyruvate. Ongoing studies will allow comparisons of patterns of angiogenesis, oxygen availability, and metabolism with tissue histology to better understand the role of these processes in tumor models of breast cancer.

Conclusion

This study demonstrates the technical feasibility of using BOLD, DCE, and HP 13C MRI to characterize mammary carcinomas. Future studies are planned to investigate regional associations of elevated R2*, perfusion/permeability, and glycolytic flux with tissue histology for different breast cancer models.

Acknowledgements

This project was supported by the RSNA Research & Education Foundation, the Department of Radiology Research and Development fund, and the Department of Medical Physics at the author’s institution. This project also received support from an AAPM Graduate Fellowship and GE Healthcare.

References

1. Newbold K, Castellano I, Charles-Edwards E, Mears D, Sohaib A, Leach M, Rhys-Evans P, Clarke P, Fisher C, Harrington K, Nutting C. An Exploratory Study Into the Role of Dynamic Contrast-Enhanced Magnetic Resonance Imaging or Perfusion Computed Tomography for Detection of Intratumoral Hypoxia in Head-and-Neck Cancer. Int. J. Radiation Oncology Biol. Phys. 2009;74(1):29-37.

2. Jiang L, Weatherall P, McColl R, Tripathy D, Mason R. Blood oxygenation level-dependent (BOLD) contrast magnetic resonance imaging (MRI) for prediction of breast cancer chemotherapy response: A pilot study. J. Magn. Reson. Imaging. 2013;37(5):1083-1092.

3. Sotgia F, Martinez-Outschoorn U, Pavlides S, Howell A, Pestell R, Lisanti M. Understanding the Warburg effect and the prognostic value of stromal caveolin-1 as a marker of a lethal tumor microenvironment. Breast Cancer Research. 2011;13(4):213.

4. Brindle K. Watching tumors gasp and die with MRI : the promise of hyperpolarized 13C MR spectroscopic imaging. The British Journal of Radiology. 2012;85:697-708.

5. Reeder S, Brittain J, Grist T, Yen Y. Least-Squares Chemical Shift Separation for 13C Metabolic Imaging. J Magn. Reson. Imaging. 2007;26:1145-1152.

6. Li J, Yu Y, Zhang Y, Bao S. A clinically feasible method to estimate pharmacokinetic parameters in breast cancer. Med. Phys. 2009;36(8):3786-3794.

7. Murase K. Efficient Method for Calculating Kinetic Parameters Using T1-Weighted Dynamic Contrast-Enhanced Magnetic Resonance Imaging. Magn. Reson.in Medicine. 2004;51:858-862.

8. Yankeelov T, Luci J, Lepage M, Li R, Debusk L, Lin P, Price R, Gore J. Quantitative pharmacokinetic analysis of DCE-MRI data without an arterial input function: a reference region model. Magn. Reson. Imaging. 2005;23:519-529.

Figures

Figure 1: Overlaid on anatomical T2-weighted images (grayscale) are maps R2* (A, D), Ktrans/ve (B, E), and median thresholding of R2* and Ktrans/ve (C, F) for a mouse with the markedly aggressive tumor strain (top row) and the moderately aggressive tumor strain (bottom row).

Figure 2: Table of quantitative DCE (Ktrans/ve) and BOLD (R2*) measurements.

Figure 3: Overlaid on anatomical T2-weighted images (grayscale) are maps of lactate AUC (A), pyruvate AUC (B), ratio of lactate AUC to pyruvate AUC (C), R2* (D), Ktrans/ve (E), and median thresholding of R2* and Ktrans/ve (F) for a mouse with the most aggressive tumor strain.



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