Cancer Foraging Ecology: Diet Choice, Patch Use, and Habitat Selection of Cancer Cells

A gene for all seasons: The evolutionary consequences of HIF-1 in carcinogenesis, tumor growth and metastasis

Oxygen played a pivotal role in the evolution of multicellularity during the Cambrian Explosion. Not surprisingly, responses to fluctuating oxygen concentrations are integral to the evolution of cancer-a disease characterized by the breakdown of multicellularity. Poorly organized tumor vasculature results in chaotic patterns of blood flow characterized by large spatial and temporal variations in intra-tumoral oxygen concentrations. Hypoxia-inducible growth factor (HIF-1) plays a pivotal role in enabling cells to adapt, metabolize, and proliferate in low oxygen conditions. HIF-1 is often constitutively activated in cancers, underscoring its importance in cancer progression. Here, we argue that the phenotypic changes mediated by HIF-1, in addition to adapting the cancer cells to their local environment, also "pre-adapt" them for proliferation at distant, metastatic sites. HIF-1-mediated adaptations include a metabolic shift towards anaerobic respiration or glycolysis, activation of cell survival mechanisms like phenotypic plasticity and epigenetic reprogramming, and formation of tumor vasculature through angiogenesis. Hypoxia induced epigenetic reprogramming can trigger epithelial to mesenchymal transition in cancer cells-the first step in the metastatic cascade. Highly glycolytic cells facilitate local invasion by acidifying the tumor microenvironment. New blood vessels, formed due to angiogenesis, provide cancer cells a conduit to the circulatory system. Moreover, survival mechanisms acquired by cancer cells in the primary site allow them to remodel tissue at the metastatic site generating tumor promoting microenvironment. Thus, hypoxia in the primary tumor promoted adaptations conducive to all stages of the metastatic cascade from the initial escape entry into a blood vessel, intravascular survival, extravasation into distant tissues, and establishment of secondary tumors.

Metastatic cells exploit their stoichiometric niche in the network of cancer ecosystems

Metastasis is a nonrandom process with varying degrees of organotropism-specific source-acceptor seeding. Understanding how patterns between source and acceptor tumors emerge remains a challenge in oncology. We hypothesize that organotropism results from the macronutrient niche of cells in source and acceptor organs. To test this, we constructed and analyzed a metastatic network based on 9303 records across 28 tissue types. We found that the topology of the network is nested and modular with scale-free degree distributions, reflecting organotropism along a specificity/generality continuum. The variation in topology is significantly explained by the matching of metastatic cells to their stoichiometric niche. Specifically, successful metastases are associated with higher phosphorus content in the acceptor compared to the source organ, due to metabolic constraints in proliferation crucial to the invasion of new tissues. We conclude that metastases are codetermined by processes at source and acceptor organs, where phosphorus content is a limiting factor orchestrating tumor ecology.

GLUT1 production in cancer cells: a tragedy of the commons

The tragedy of the commons occurs when competition among individual members of a group leads to overexploitation of a shared resource to the detriment of the overall population. We hypothesize that cancer cells may engage in a tragedy of the commons when competing for a shared resource such as glucose. To formalize this notion, we create a game theoretic model of glucose uptake based on a cell’s investment in transporters relative to that of its neighboring cells. We show that production of transporters per cell increases as the number of competing cells in a microenvironment increases and nutrient uptake per cell decreases. Furthermore, the greater the resource availability, the more intense the tragedy of the commons at the ESS. Based on our simulations, cancer cells produce 2.2–2.7 times more glucose transporters than would produce optimal fitness for all group members. A tragedy of the commons affords novel therapeutic strategies. By simulating GLUT1 inhibitor and glucose deprivation treatments, we demonstrate a synergistic combination with standard-of-care therapies, while also displaying the existence of a trade-off between competition among cancer cells and depression of their gain function. Assuming cancer cell transporter production is heritable, we then show the potential for a sucker’s gambit therapy by exploiting this trade-off. By strategically changing environmental conditions, we can take advantage of cellular competition and gain function depression.

Cancer cell foraging to explain bone-specific metastatic progression

Metastatic cancer is lethal and patients who suffer bone metastases fare especially poorly. Bone-specific metastatic progression in prostate and breast cancers is a highly observed example of organ-specific metastasis, or organotropism. Though research has delineated the sequential steps of the metastatic cascade, the determinants of bone-specific metastasis have remained elusive for decades. Applying fundamental ecological principles to cancer biology models of metastasis provides novel insights into metastatic organotropism. We use critical concepts from foraging theory and movement ecology to propose that observed bone-specific metastasis is the result of habitat selection by foraging cancer cells. Furthermore, we posit that cancer cells can only perform habitat selection if and when they employ a reversible motile foraging strategy. Only a very small percentage of cells in a primary tumor harbor this ability. Therefore, our habitat selection model emphasizes the importance of identifying the rare subset of cancer cells that might exhibit habitat selection, ergo achieve bone-specific metastatic colonization.

Microeconomics of Metabolism: The Warburg Effect as Giffen Behaviour

Metabolic behaviours of proliferating cells are often explained as a consequence of rational optimization of cellular growth rate, whereas microeconomics formulates consumption behaviours as optimization problems. Here, we pushed beyond the analogy to precisely map metabolism onto the theory of consumer choice. We thereby revealed the correspondence between long-standing mysteries in both fields: the Warburg effect, a seemingly wasteful but ubiquitous strategy where cells favour aerobic glycolysis over more energetically efficient oxidative phosphorylation, and Giffen behaviour, the unexpected consumer behaviour where a good is demanded more as its price rises. We identified the minimal, universal requirements for the Warburg effect: a trade-off between oxidative phosphorylation and aerobic glycolysis and complementarity, i.e. impossibility of substitution for different metabolites. Thus, various hypotheses for the Warburg effect are integrated into an identical optimization problem with the same universal structure. Besides, the correspondence between the Warburg effect and Giffen behaviour implies that oxidative phosphorylation is counter-intuitively stimulated when its efficiency is decreased by metabolic perturbations such as drug administration or mitochondrial dysfunction; the concept of Giffen behaviour bridges the Warburg effect and the reverse Warburg effect. This highlights that the application of microeconomics to metabolism can offer new predictions and paradigms for both biology and economics.

Coexistence of “Cream Skimmer” and “Crumb Picker” Phenotypes in Nature and in Cancer

Over 40 years ago, seminal papers by Armstrong and McGehee and by Levins showed that temporal fluctuations in resource availability could permit coexistence of two species on a single resource. Such coexistence results from non-linearities or non-additivities in the way resource supply translates into fitness. These reflect trade-offs where one species benefits more than the other during good periods and suffers more (or does less well) than the other during less good periods, be the periods stochastic, unstable population dynamics, or seasonal. Since, coexistence based on fluctuating conditions has been explored under the guises of “grazers” and “diggers,” variance partitioning, relative non-linearity, “opportunists” and “gleaners,” and as the storage effect. Here we focus on two phenotypes, “cream skimmers” and “crumb pickers,” the former having the advantage in richer times and the latter in less rich times. In nature, richer and poorer times, with regular or stochastic appearances, are the norm and occur on many time scales. Fluctuations among richer and poorer times also appear to be the norm in cancer ecosystems. Within tumors, nutrient availability, oxygen, and pH can fluctuate stochastically or periodically, with swings occurring over seconds to minutes to hours. Despite interest in tumor heterogeneity and how it promotes the coexistence of different cancer cell types, the effects of fluctuating resource availability have not been explored for cancer. Here, in the context of pulsed resources, we (1) develop models of foraging consumers who experience pulsed resources to examine four types of trade-offs that can promote coexistence of phenotypes that do relatively better in richer versus in poorer times, (2) establish that conditions in tumors are conducive for this mechanism, (3) propose and empirically explore biomarkers indicative of the two phenotypes (HIF-1, GLUT-1, CA IX, CA XII), and (4) and compare cream skimmer and crumb picker biology and ecology in nature and cancer to provide cross-disciplinary insights into this interesting, and, we argue, likely very common, mechanism of coexistence.

Predator-Prey in Tumor-Immune Interactions: A Wrong Model or Just an Incomplete One?

Tumor-immune interactions are often framed as predator-prey. This imperfect analogy describes how immune cells (the predators) hunt and kill immunogenic tumor cells (the prey). It allows for evaluation of tumor cell populations that change over time during immunoediting and it also considers how the immune system changes in response to these alterations. However, two aspects of predator-prey type models are not typically observed in immuno-oncology. The first concerns the conversion of prey killed into predator biomass. In standard predator-prey models, the predator relies on the prey for nutrients, while in the tumor microenvironment the predator and prey compete for resources (e.g. glucose). The second concerns oscillatory dynamics. Standard predator-prey models can show a perpetual cycling in both prey and predator population sizes, while in oncology we see increases in tumor volume and decreases in infiltrating immune cell populations. Here we discuss the applicability of predator-prey models in the context of cancer immunology and evaluate possible causes for discrepancies. Key processes include "safety in numbers", resource availability, time delays, interference competition, and immunoediting. Finally, we propose a way forward to reconcile differences between model predictions and empirical observations. The immune system is not just predator-prey. Like natural food webs, the immune-tumor community of cell types forms an immune-web of different and identifiable interactions.

Estrogen as an Essential Resource and the Coexistence of ER+ and ER– Cancer Cells

Diagnosis of estrogen sensitivity in breast cancer is largely predicated on the ratio of ER+ and ER– cancer cells obtained from biopsies. Estrogen is a growth factor necessary for cell survival and division. It can also be thought of as an essential resource that can act in association with other nutrients, glucose, glutamine, fatty acids, amino acids, etc. All of these nutrients, collectively or individually, may limit the growth of the cancer cells (Liebig’s Law of the Minimum). Here we model estrogen susceptibility in breast cancer as a consumer-resource interaction: ER+ cells require both estrogen and glucose as essential resources, whereas ER– only require the general resource. The model predicts that when estrogen is the limiting factor, other nutrients may go unconsumed and available at higher levels, thus permitting the invasion of ER– cells. Conversely, when ER– cells are less efficient on glucose than ER+ cells, then ER– cells limited by glucose may be susceptible to invasion by ER+ cells, provided that sufficient levels of estrogen are available. ER+ cells will outcompete ER– cells when estrogen is abundant, resulting in low concentrations of interstitial glucose within the tumor. In the absence of estrogen, ER– cells will outcompete ER+ cells, leaving a higher concentration of interstitial glucose. At intermediate delivery rates of estrogen and glucose, ER+ and ER– cells are predicted to coexist. In modeling the dynamics of cells in the same tumor with different resource requirements, we can apply concepts and terms familiar to many ecologists. These include: resource supply points, R∗, ZNGI (zero net growth isoclines), resource depletion, and resource uptake rates. Based on the circumstances favoring ER+ vs. ER– breast cancer, we use the model to explore the consequences of therapeutic regimens that may include hormonal therapies, possible roles of diet in changing cancer cell composition, and potential for evolutionarily informed therapies. More generally, the model invites the viewpoint that cancer’s eco-evolutionary dynamics are a consumer-resource interaction, and that other growth factors such as EGFR or androgens may be best viewed as essential resources within these dynamics.

Cancer Community Ecology

Here we advocate Cancer Community Ecology as a valuable focus of study in Cancer Biology. We hypothesize that the heterogeneity and characteristics of cancer cells within tumors should vary systematically in space and time and that cancer cells form local ecological communities within tumors. These communities possess limited numbers of species determined by local conditions, with each species in a community possessing predictable traits that enable them to cope with their particular environment and coexist with each other. We start with a discussion of concepts and assumptions that ecologists use to study closely related species. We then discuss the competitive exclusion principle as a means for knowing when two species should not coexist, and as an opening towards understanding how they can. We present the five major categories of mechanisms of coexistence that operate in nature and suggest that the same mechanisms apply towards understanding the diversification and coexistence of cancer cell species. They are: Food-Safety Tradeoffs, Diet Choice, Habitat Selection, Variance Partitioning, and Competition-Colonization Tradeoffs. For each mechanism, we discuss how it works in nature, how it might work in cancers, and its implications for therapy.

Frequency-dependent interactions determine outcome of competition between two breast cancer cell lines

Tumors are highly dynamic ecosystems in which diverse cancer cell subpopulations compete for space and resources. These complex, often non-linear interactions govern continuous spatial and temporal changes in the size and phenotypic properties of these subpopulations. Because intra-tumoral blood flow is often chaotic, competition for resources may be a critical selection factor in progression and prognosis. Here, we quantify resource competition using 3D spheroid cultures with MDA-MB-231 and MCF-7 breast cancer cells. We hypothesized that MCF-7 cells, which primarily rely on efficient aerobic glucose metabolism, would dominate the population under normal pH and low glucose conditions; and MDA-MB-231 cells, which exhibit high levels of glycolytic metabolism, would dominate under low pH and high glucose conditions. In spheroids with single populations, MCF-7 cells exhibited equal or superior intrinsic growth rates (density-independent measure of success) and carrying capacities (density-dependent measure of success) when compared to MDA-MB-231 cells under all pH and nutrient conditions. Despite these advantages, when grown together, MCF-7 cells do not always outcompete MDA-MB-231 cells. MDA-MB-231 cells outcompete MCF-7 cells in low glucose conditions and coexistence is achieved in low pH conditions. Under all conditions, MDA-MB-231 has a stronger competitive effect (frequency-dependent interaction) on MCF-7 cells than vice-versa. This, and the inability of growth rate or carrying capacity when grown individually to predict the outcome of competition, suggests a reliance on frequency-dependent interactions and the need for competition assays. We frame these results in a game-theoretic (frequency-dependent) model of cancer cell interactions and conclude that competition assays can demonstrate critical density-independent, density-dependent and frequency-dependent interactions that likely contribute to in vivo outcomes.

A Mathematical Dissection of the Adaptation of Cell Populations to Fluctuating Oxygen Levels

The disordered network of blood vessels that arises from tumour angiogenesis results in variations in the delivery of oxygen into the tumour tissue. This brings about regions of chronic hypoxia (i.e. sustained low oxygen levels) and regions with alternating periods of low and relatively higher oxygen levels, and makes it necessary for cancer cells to adapt to fluctuating environmental conditions. We use a phenotype-structured model to dissect the evolutionary dynamics of cell populations exposed to fluctuating oxygen levels. In this model, the phenotypic state of every cell is described by a continuous variable that provides a simple representation of its metabolic phenotype, ranging from fully oxidative to fully glycolytic, and cells are grouped into two competing populations that undergo heritable, spontaneous phenotypic variations at different rates. Model simulations indicate that, depending on the rate at which oxygen is consumed by the cells, dynamic nonlinear interactions between cells and oxygen can stimulate chronic hypoxia and cycling hypoxia. Moreover, the model supports the idea that under chronic-hypoxic conditions lower rates of phenotypic variation lead to a competitive advantage, whereas higher rates of phenotypic variation can confer a competitive advantage under cycling-hypoxic conditions. In the latter case, the numerical results obtained show that bet-hedging evolutionary strategies, whereby cells switch between oxidative and glycolytic phenotypes, can spontaneously emerge. We explain how these results can shed light on the evolutionary process that may underpin the emergence of phenotypic heterogeneity in vascularised tumours.

Cancer cell lines show high heritability for motility but not generation time

Tumour evolution depends on heritable differences between cells in traits affecting cell survival or replication. It is well established that cancer cells are genetically and phenotypically heterogeneous; however, the extent to which this phenotypic variation is heritable is far less well explored. Here, we estimate the broad-sense heritability (H ²) of two cell traits related to cancer hallmarks--cell motility and generation time--within populations of four cancer cell lines in vitro and find that motility is strongly heritable. This heritability is stable across multiple cell generations, with heritability values at the high end of those measured for a range of traits in natural populations of animals or plants. These findings confirm a central assumption of cancer evolution, provide a first quantification of the evolvability of key traits in cancer cells and indicate that there is ample raw material for experimental evolution in cancer cell lines. Generation time, a trait directly affecting cell fitness, shows substantially lower values of heritability than cell speed, consistent with its having been under directional selection removing heritable variation.