A mathematical investigation of the interplay between vasculature and intratumoral cellular heterogeneity during tumor progression

Tumors consist of heterogeneous phenotypic cells, such as normoxic cells, which are highly proliferative, and hypoxic cells, which are less proliferative. Their phenotypic switching depends on tumor microenvironmental factors, such as oxygen and nutrient concentrations supplied by local blood vessels. However, during ongoing angiogenesis, the process of sprouting new blood vessels at the tumor site from pre-existing blood vessels, and how this phenotypic switching affects and impacts tumor growth, remains poorly understood. In this article, we formulate a mathematical model to elucidate the crosstalk between vasculature and tumor cellular heterogeneity during tumor progression. The model results show a strong agreement with the experimental data. Our simulation results demonstrate that ongoing angiogenesis increases tumor growth rate. In addition, we observe that the influence of hypoxic cells on phenotypic switching from normoxic to hypoxic is more pronounced than their influence on the transition from hypoxic to normoxic. Furthermore, we perform a global sensitivity analysis using the Sobol’s method to assess the importance of the model’s parameters. It highlights that the volume at which blood vessels attain half-maximal rate has the maximum effect on the model.