Abstract
Histopathologic knowledge that extensive heterogeneity exists between and within tumors has been confirmed and deepened recently by molecular studies. However, the impact of tumor heterogeneity on prognosis and treatment remains as poorly understood as ever. Using a hybrid multiscale mathematical model of tumor growth in vascularized tissue, we investigated the selection pressures exerted by spatial and temporal variations in tumor microenvironment and the resulting phenotypic adaptations. A key component of this model is normal and tumor metabolism and its interaction with microenvironmental factors. The metabolic phenotype of tumor cells is plastic, and microenvironmental selection leads to increased tumor glycolysis and decreased pH. Once this phenotype emerges, the tumor dramatically changes its behavior due to acid-mediated invasion, an effect that depends on both variations in the tumor cell phenotypes and their spatial distribution within the tumor. In early stages of growth, tumors are stratified, with the most aggressive cells developing within the interior of the tumor. These cells then grow to the edge of the tumor and invade into the normal tissue using acidosis. Simulations suggest that diffusible cytotoxic treatments, such as chemotherapy, may increase the metabolic aggressiveness of a tumor due to drug-mediated selection. Chemotherapy removes the metabolic stratification of the tumor and allows more aggressive cells to grow toward blood vessels and normal tissue. Antiangiogenic therapy also selects for aggressive phenotypes due to degradation of the tumor microenvironment, ultimately resulting in a more invasive tumor. In contrast, pH buffer therapy slows down the development of aggressive tumors, but only if administered when the tumor is still stratified. Overall, findings from this model highlight the risks of cytotoxic and antiangiogenic treatments in the context of tumor heterogeneity resulting from a selection for more aggressive behaviors.