Understanding temporal and spatial hemodynamic heterogeneity as a function of tumor growth or therapy impacts the development of novel therapeutic strategies. In this study, we employed eigenspectra multispectral optoacoustic tomography (eMSOT) as a next-generation optoacoustic method to impart high accuracy in resolving tumor hemodynamics during bevacizumab therapy in two types of breast cancer xenografts (KPL-4, MDA-MB-468). Patterns of tumor total hemoglobin concentration (THb) and oxygen saturation (sO2) were imaged in control and bevacizumab-treated tumors over the course of 58d (KPL-4) and 16d (MDA-MB-468), and the evolution of functional vasculature “normalization” was resolved macroscopically. An initial sharp drop in tumor sO2 and THb content shortly after the initiation of bevacizumab treatment was followed by a recovery in oxygenation levels. Rim-core subregion analysis revealed steep spatial oxygenation gradients in growing tumors that were reduced after bevacizumab treatment. Critically, eMSOT imaging findings were validated directly by histopathologic assessment of hypoxia (pimonidazole) and vascularity (CD31). These data demonstrate how eMSOT brings new abilities for accurate observation of entire tumor responses to challenges at spatial and temporal dimensions not available by other techniques today.
