Drug-resilient cancer cell phenotype is acquired via polyploidization associated with early stress response coupled to HIF-2α transcriptional regulation

Therapeutic resistance and recurrence remain core challenges in cancer therapy. How therapy resistance arises is currently not fully understood with tumors surviving via multiple alternative routes. Here, we demonstrate that a subset of cancer cells survives therapeutic stress by entering a transient state characterized by whole genome doubling. At the onset of the polyploidization program, we identified an upregulation of key transcriptional regulators, including the early stress-response protein AP-1 and normoxic stabilization of HIF-2α. We found altered chromatin accessibility, ablated expression of RB1, and enrichment of AP-1 motif accessibility. We demonstrate that AP-1 and HIF-2α regulate a therapy resilient and survivor phenotype in cancer cells. Consistent with this, genetic or pharmacologic targeting of AP-1 and HIF-2α reduced the number of surviving cells following chemotherapy treatment. The role of AP-1 and HIF-2α in stress-response by polyploidy suggest a novel avenue for tackling chemotherapy-induced resistance in cancer

Evolutionary Game Theory: Darwinian Dynamics and the G Function Approach

Classical evolutionary game theory allows one to analyze the population dynamics of interacting individuals playing different strategies (broadly defined) in a population. To expand the scope of this framework to allow us to examine the evolution of these individuals’ strategies over time, we present the idea of a fitness-generating (G) function. Under this model, we can simultaneously consider population (ecological) and strategy (evolutionary) dynamics. In this paper, we briefly outline the differences between game theory and classical evolutionary game theory. We then introduce the G function framework, deriving the model from fundamental biological principles. We introduce the concept of a G-function species, explain the process of modeling with G functions, and define the conditions for evolutionary stable strategies (ESS). We conclude by presenting expository examples of G function model construction and simulations in the context of predator–prey dynamics and the evolution of drug resistance in cancer

Evolutionary Game Theory: Darwinian Dynamics and the <i>G</i> Function Approach

Classical evolutionary game theory allows one to analyze the population dynamics of interacting individuals playing different strategies (broadly defined) in a population. To expand the scope of this framework to allow us to examine the evolution of these individuals’ strategies over time, we present the idea of a fitness-generating (G) function. Under this model, we can simultaneously consider population (ecological) and strategy (evolutionary) dynamics. In this paper, we briefly outline the differences between game theory and classical evolutionary game theory. We then introduce the G function framework, deriving the model from fundamental biological principles. We introduce the concept of a G-function species, explain the process of modeling with G functions, and define the conditions for evolutionary stable strategies (ESS). We conclude by presenting expository examples of G function model construction and simulations in the context of predator–prey dynamics and the evolution of drug resistance in cancer

Integrating eco‐evolutionary dynamics into matrix population models for structured populations: Discrete and continuous frameworks

Abstract State‐structured populations are ubiquitous in biology, from the age‐structure of animal societies to the life cycles of parasitic species. Understanding how this structure contributes to eco‐evolutionary dynamics is critical not only for fundamental understanding but also for conservation and treatment purposes. Although some methods have been developed in the literature for modelling eco‐evolutionary dynamics in structured population, such methods are wholly lacking in the G function evolutionary game theoretic framework. In this paper, we integrate standard matrix population modelling into the G function framework to create a theoretical framework to probe eco‐evolutionary dynamics in structured populations. This framework encompasses age‐ and stage‐structured matrix models with basic density‐ and frequency‐dependent transition rates and probabilities. For both discrete and continuous time models, we define and characterize asymptotic properties of the system such as eco‐evolutionary equilibria (including ESSs) and the convergence stability of these equilibria. For multistate structured populations, we introduce an ergodic flow preserving folding method for analysing such models. The methods developed in this paper for state‐structured populations and their extensions to multistate‐structured populations provide a simple way to create, analyse and simulate eco‐evolutionary dynamics in structured populations. Furthermore, their generality allows these techniques to be applied to a variety of problems in ecology and evolution

Stochastic models of Mendelian and reverse transcriptional inheritance in state-structured cancer populations

Recent evidence suggests that a polyaneuploid cancer cell (PACC) state may play a key role in the adaptation of cancer cells to stressful environments and in promoting therapeutic resistance. The PACC state allows cancer cells to pause cell division and to avoid DNA damage and programmed cell death. Transition to the PACC state may also lead to an increase in the cancer cell’s ability to generate heritable variation (evolvability). One way this can occur is through evolutionary triage. Under this framework, cells gradually gain resistance by scaling hills on a fitness landscape through a process of mutation and selection. Another way this can happen is through self-genetic modification whereby cells in the PACC state find a viable solution to the stressor and then undergo depolyploidization, passing it on to their heritably resistant progeny. Here, we develop a stochastic model to simulate both of these evolutionary frameworks. We examine the impact of treatment dosage and extent of self-genetic modification on eco-evolutionary dynamics of cancer cells with aneuploid and PACC states. We find that under low doses of therapy, evolutionary triage performs better whereas under high doses of therapy, self-genetic modification is favored. This study generates predictions for teasing apart these biological hypotheses, examines the implications of each in the context of cancer, and provides a modeling framework to compare Mendelian and non-traditional forms of inheritance

A life history model of the ecological and evolutionary dynamics of polyaneuploid cancer cells

Therapeutic resistance is one of the main reasons for treatment failure in cancer patients. The polyaneuploid cancer cell (PACC) state has been shown to promote resistance by providing a refuge for cancer cells from the effects of therapy and by helping them adapt to a variety of environmental stressors. This state is the result of aneuploid cancer cells undergoing whole genome doubling and skipping mitosis, cytokinesis, or both. In this paper, we create a novel mathematical framework for modeling the eco-evolutionary dynamics of state-structured populations and use this framework to construct a model of cancer populations with an aneuploid and a PACC state. Using in silico simulations, we explore how the PACC state allows cancer cells to (1) survive extreme environmental conditions by exiting the cell cycle after S phase and protecting genomic material and (2) aid in adaptation to environmental stressors by increasing the cancer cell’s ability to generate heritable variation (evolvability) through the increase in genomic content that accompanies polyploidization. In doing so, we demonstrate the ability of the PACC state to allow cancer cells to persist under therapy and evolve therapeutic resistance. By eliminating cells in the PACC state through appropriately-timed PACC-targeted therapies, we show how we can prevent the emergence of resistance and promote cancer eradication

A mathematical investigation of polyaneuploid cancer cell memory and cross-resistance in state-structured cancer populations

Abstract The polyaneuploid cancer cell (PACC) state promotes cancer lethality by contributing to survival in extreme conditions and metastasis. Recent experimental evidence suggests that post-therapy PACC-derived recurrent populations display cross-resistance to classes of therapies with independent mechanisms of action. We hypothesize that this can occur through PACC memory, whereby cancer cells that have undergone a polyaneuploid transition (PAT) reenter the PACC state more quickly or have higher levels of innate resistance. In this paper, we build on our prior mathematical models of the eco-evolutionary dynamics of cells in the 2N+ and PACC states to investigate these two hypotheses. We show that although an increase in innate resistance is more effective at promoting cross-resistance, this trend can also be produced via PACC memory. We also find that resensitization of cells that acquire increased innate resistance through the PAT have a considerable impact on eco-evolutionary dynamics and extinction probabilities. This study, though theoretical in nature, can help inspire future experimentation to tease apart hypotheses surrounding how cross-resistance in structured cancer populations arises