Vascular phenotyping of the invasive front in breast cancer using a 3D angiogenesis atlas

Vascular remodeling at the invasive tumor front (ITF) plays a critical role in progression and metastasis of triple negative breast cancer (TNBC). Therefore, there is a crucial need to characterize the vascular phenotype (i.e. changes in the structure and function of vasculature) of the ITF and tumor core (TC) in TNBC. This requires high-resolution, 3D structural and functional microvascular data that spans the ITF and TC (i.e. ∼4–5 mm from the tumor's edge). Since such data are often challenging to obtain with most conventional imaging approaches, we employed a unique "3D whole-tumor angiogenesis atlas" derived from orthotopic xenografts to characterize the vascular phenotype of the ITF and TC in TNBC. First, high-resolution (8 μm) computed tomography (CT) images of "whole-tumor" microvasculature were acquired from eight orthotopic TNBC xenografts, of which three tumors were excised at post-inoculation day 21 (i.e. early-stage) and five tumors were excised at post-inoculation day 35 (i.e. advanced-stage). These 3D morphological CT data were combined with soft tissue contrast from MRI as well as functional data generated in silico using image-based hemodynamic modeling to generate a multi-layered "angiogenesis atlas". Employing this atlas, blood vessels were first spatially stratified within the ITF (i.e. ≤1 mm from the tumor's edge) and TC (i.e. >1 mm from the tumor's edge) of each tumor xenograft. Then, a novel method was developed to visualize and characterize microvascular remodeling and perfusion changes in terms of distance from the tumor's edge. The angiogenesis atlas enabled the 3D visualization of changes in tumor vessel growth patterns, morphology and perfusion within the ITF and TC. Early and advanced stage tumors demonstrated significant differences in terms of their edge-to-center distributions for vascular surface area density, vascular length density, intervessel distance and simulated perfusion density (p ≪ 0.01). Elevated vascular length density, vascular surface area density and perfusion density along the circumference of the ITF was suggestive of a preferential spatial pattern of angiogenic growth in this tumor cohort. Finally, we demonstrated the feasibility of differentiating the vascular phenotypes of ITF and TC in these TNBC xenografts. The combination of a 3D angiogenesis atlas and image-based hemodynamic modeling heralds a new approach for characterizing the role of vascular remodeling in cancer and other diseases.