Spatial vs. non-spatial eco-evolutionary dynamics in a tumor growth model

Spatial interactions modulate tumor growth and immune infiltration

Direct observation of tumor-immune interactions is unlikely in tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen evolve as a mechanism of immune escape, but the exact nature of immune-collagen interactions is poorly understood. Spatial data quantifying collagen fiber alignment in squamous cell carcinomas indicates that late-stage disease is associated with highly aligned fibers. Our computational modeling framework discriminates between two hypotheses: immune cell migration that moves (1) parallel or (2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-extracellular matrix interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. Here, computational modeling provides important mechanistic insights by defining a kernel cell-cell interaction function that considers a spectrum of local (cell-scale) to global (tumor-scale) spatial interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interaction kernels drives tumor growth and infiltration.

The effect of tumor composition on the success of adaptive therapy: The case of metastatic Castrate-Resistant Prostate Cancer

Prostate-specific antigen (PSA) is the most commonly used serum marker for prostate cancer. It plays a role in cancer detection, treatment monitoring, and more recently, in guiding adaptive therapy protocols, where treatment is alternated based on PSA levels. However, the relationship between PSA levels and tumor volume remains poorly understood. Empirical evidence suggests that different cancer cell types produce varying amounts of PSA. Despite this, current mathematical cancer models often assume either that all cell types contribute equally to PSA levels or that only certain subpopulations produce PSA at fixed rates. In this study, we compare Zhang et al.'s classical adaptive therapy protocol with the standard of care, which involves continuous maximum tolerable dose treatment, under different assumptions regarding PSA production. Specifically, we explore the possibility that testosterone-dependent, testosterone-producing, and testosterone-independent cells contribute to PSA production to varying degrees. We use the time to competitive release as a proxy for the time to disease progression. Our findings indicate that adaptive therapy consistently results in a longer time to competitive release compared to the standard of care, regardless of the assumptions about PSA production. However, when testosterone-independent cells are the sole PSA producers, Zhang et al.'s adaptive therapy protocol becomes inapplicable, as PSA levels never fall to half of their initial value, preventing therapy discontinuation. Additionally, we observe that the number and duration of treatment cycles in adaptive therapy are highly sensitive to assumptions about how much each cell type contributes to PSA production. Overall, our results emphasize the need for a deeper understanding of patient-specific PSA dynamics, which could enhance the effectiveness of adaptive therapy in prostate cancer treatment.

Spatial interactions modulate tumor growth and immune infiltration

Direct observation of immune cell trafficking patterns and tumor-immune interactions is unlikely in human tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen formation evolve as a mechanism of immune escape, but the exact nature of the interaction between immune cells and collagen is poorly understood. Spatial data quantifying the degree of collagen fiber alignment in squamous cell carcinomas indicates that late stage disease is associated with highly aligned fibers. Here, we introduce a computational modeling framework (called Lenia) to discriminate between two hypotheses: immune cell migration that moves 1) parallel or 2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-ECM interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. We also illustrate the capabilities of Lenia to model the evolution of tumor progression and immune predation. Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining a kernel cell-cell interaction function that governs tumor growth dynamics under immune predation with immune cell migration. Mathematical modeling provides important mechanistic insights into cell interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interactions drives tumor growth and infiltration.

Intermetastatic and Intrametastatic Heterogeneity Shapes Adaptive Therapy Cycling Dynamics

Adaptive therapies that alternate between drug applications and drug-free vacations can exploit competition between sensitive and resistant cells to maximize the time to progression. However, optimal dosing schedules depend on the properties of metastases, which are often not directly measurable in clinical practice. Here, we proposed a framework for estimating features of metastases through tumor response dynamics during the first adaptive therapy treatment cycle. Longitudinal prostate-specific antigen (PSA) levels in 16 patients with metastatic castration-resistant prostate cancer undergoing adaptive androgen deprivation treatment were analyzed to investigate relationships between cycle dynamics and clinical variables such as Gleason score, the change in the number of metastases over a cycle, and the total number of cycles over the course of treatment. The first cycle of adaptive therapy, which consists of a response period (applying therapy until 50% PSA reduction), and a regrowth period (removing treatment until reaching initial PSA levels), delineated several features of the computational metastatic system: larger metastases had longer cycles; a higher proportion of drug-resistant cells slowed the cycles; and a faster cell turnover rate sped up drug response time and slowed regrowth time. The number of metastases did not affect cycle times, as response dynamics were dominated by the largest tumors rather than the aggregate. In addition, systems with higher intermetastasis heterogeneity responded better to continuous therapy and correlated with dynamics from patients with high or low Gleason scores. Conversely, systems with higher intrametastasis heterogeneity responded better to adaptive therapy and correlated with dynamics from patients with intermediate Gleason scores.

SIGNIFICANCE

Multiscale mathematical modeling combined with biomarker dynamics during adaptive therapy helps identify underlying features of metastatic cancer to inform treatment decisions.

State-dependent evolutionary models reveal modes of solid tumour growth

Spatial properties of tumour growth have profound implications for cancer progression, therapeutic resistance and metastasis. Yet, how spatial position governs tumour cell division remains difficult to evaluate in clinical tumours. Here, we demonstrate that faster division on the tumour periphery leaves characteristic genetic patterns, which become evident when a phylogenetic tree is reconstructed from spatially sampled cells. Namely, rapidly dividing peripheral lineages branch more extensively and acquire more mutations than slower-dividing centre lineages. We develop a Bayesian state-dependent evolutionary phylodynamic model (SDevo) that quantifies these patterns to infer the differential division rates between peripheral and central cells. We demonstrate that this approach accurately infers spatially varying birth rates of simulated tumours across a range of growth conditions and sampling strategies. We then show that SDevo outperforms state-of-the-art, non-cancer multi-state phylodynamic methods that ignore differential sequence evolution. Finally, we apply SDevo to single-time-point, multi-region sequencing data from clinical hepatocellular carcinomas and find evidence of a three- to six-times-higher division rate on the tumour edge. With the increasing availability of high-resolution, multi-region sequencing, we anticipate that SDevo will be useful in interrogating spatial growth restrictions and could be extended to model non-spatial factors that influence tumour progression.

A survey of open questions in adaptive therapy: Bridging mathematics and clinical translation

Adaptive therapy is a dynamic cancer treatment protocol that updates (or 'adapts') treatment decisions in anticipation of evolving tumor dynamics. This broad term encompasses many possible dynamic treatment protocols of patient-specific dose modulation or dose timing. Adaptive therapy maintains high levels of tumor burden to benefit from the competitive suppression of treatment-sensitive subpopulations on treatment-resistant subpopulations. This evolution-based approach to cancer treatment has been integrated into several ongoing or planned clinical trials, including treatment of metastatic castrate resistant prostate cancer, ovarian cancer, and BRAF-mutant melanoma. In the previous few decades, experimental and clinical investigation of adaptive therapy has progressed synergistically with mathematical and computational modeling. In this work, we discuss 11 open questions in cancer adaptive therapy mathematical modeling. The questions are split into three sections: (1) integrating the appropriate components into mathematical models (2) design and validation of dosing protocols, and (3) challenges and opportunities in clinical translation.

Modeling cancer's ecological and evolutionary dynamics

In this didactic paper, we present a theoretical modeling framework, called the G-function, that integrates both the ecology and evolution of cancer to understand oncogenesis. The G-function has been used in evolutionary ecology, but has not been widely applied to problems in cancer. Here, we build the G-function framework from fundamental Darwinian principles and discuss how cancer can be seen through the lens of ecology, evolution, and game theory. We begin with a simple model of cancer growth and add on components of cancer cell competition and drug resistance. To aid in exploration of eco-evolutionary modeling with this approach, we also present a user-friendly software tool. By the end of this paper, we hope that readers will be able to construct basic G function models and grasp the usefulness of the framework to understand the games cancer plays in a biologically mechanistic fashion.

Abiraterone-Docetaxel scheduling for metastatic castration-resistant prostate cancer based on evolutionary dynamics

Patients with metastatic castration-resistant prostate cancer (mCRPC) are divided into three groups based on their response to Abiraterone treatment: best responder, responder, and non-responder. In the latter two groups, successful outcomes may not be achieved due to the development of drug-resistant cells in the tumor environment during treatment. To overcome this challenge, a secondary drug can be used to control the population of drug-resistant cells, potentially leading to a longer period of disease inhibition. This paper proposes using a combination of Docetaxel and Abiraterone in some polytherapy methods to control both the overall cancer cell population and the drug-resistant subpopulation. To investigate the competition and evolution of mCRPC cancer phenotypes, as in previous studies, the Evolutionary Game Theory (EGT) has been used as a mathematical modeling of evolutionary biology concepts.

The role of evolutionary game theory in spatial and non-spatial models of the survival of cooperation in cancer: a review

Evolutionary game theory (EGT) is a branch of mathematics which considers populations of individuals interacting with each other to receive pay-offs. An individual’s pay-off is dependent on the strategy of its opponent(s) as well as on its own, and the higher its pay-off, the higher its reproductive fitness. Its offspring generally inherit its interaction strategy, subject to random mutation. Over time, the composition of the population shifts as different strategies spread or are driven extinct. In the last 25 years there has been a flood of interest in applying EGT to cancer modelling, with the aim of explaining how cancerous mutations spread through healthy tissue and how intercellular cooperation persists in tumour-cell populations. This review traces this body of work from theoretical analyses of well-mixed infinite populations through to more realistic spatial models of the development of cooperation between epithelial cells. We also consider work in which EGT has been used to make experimental predictions about the evolution of cancer, and discuss work that remains to be done before EGT can make large-scale contributions to clinical treatment and patient outcomes.

Classical mathematical models for prediction of response to chemotherapy and immunotherapy

Classical mathematical models of tumor growth have shaped our understanding of cancer and have broad practical implications for treatment scheduling and dosage. However, even the simplest textbook models have been barely validated in real world-data of human patients. In this study, we fitted a range of differential equation models to tumor volume measurements of patients undergoing chemotherapy or cancer immunotherapy for solid tumors. We used a large dataset of 1472 patients with three or more measurements per target lesion, of which 652 patients had six or more data points. We show that the early treatment response shows only moderate correlation with the final treatment response, demonstrating the need for nuanced models. We then perform a head-to-head comparison of six classical models which are widely used in the field: the Exponential, Logistic, Classic Bertalanffy, General Bertalanffy, Classic Gompertz and General Gompertz model. Several models provide a good fit to tumor volume measurements, with the Gompertz model providing the best balance between goodness of fit and number of parameters. Similarly, when fitting to early treatment data, the general Bertalanffy and Gompertz models yield the lowest mean absolute error to forecasted data, indicating that these models could potentially be effective at predicting treatment outcome. In summary, we provide a quantitative benchmark for classical textbook models and state-of-the art models of human tumor growth. We publicly release an anonymized version of our original data, providing the first benchmark set of human tumor growth data for evaluation of mathematical models.

Mathematical oncology uses quantitative models for prediction of tumor growth and treatment response. The theoretical foundation of mathematical oncology is provided by six classical mathematical models: the Exponential, Logistic, Classic Bertalanffy, General Bertalanffy, Classic Gompertz and General Gompertz model. These models have been introduced decades ago, have been used in thousands of scientific articles and are part of textbooks and curricula in mathematical oncology. However, these models have not been systematically tested in clinical data from actual patients. In this study, we have collected quantitative tumor volume measurements from thousands of patients in five large clinical trials of cancer immunotherapy. We use this dataset to systematically investigate how accurately mathematical models can describe tumor growth, showing that there are pronounced differences between models. In addition, we show that two of these models can predict tumor response to immunotherapy and chemotherapy at later time points when trained on early tumor growth dynamics. Thus, our article closes a conceptual gap in the literature and at the same time provides a simple tool to predict response to chemotherapy and immunotherapy on the level of individual patients.

Coordination games in cancer

We propose a model of cancer initiation and progression where tumor growth is modulated by an evolutionary coordination game. Evolutionary games of cancer are widely used to model frequency-dependent cell interactions with the most studied games being the Prisoner's Dilemma and public goods games. Coordination games, by their more obscure and less evocative nature, are left understudied, despite the fact that, as we argue, they offer great potential in understanding and treating cancer. In this paper we present the conditions under which coordination games between cancer cells evolve, we propose aspects of cancer that can be modeled as results of coordination games, and explore the ways through which coordination games of cancer can be exploited for therapy.

Coordination games in cancer

We propose a model of cancer initiation and progression where tumor growth is modulated by an evolutionary coordination game. Evolutionary games of cancer are widely used to model frequency-dependent cell interactions with the most studied games being the Prisoner's Dilemma and public goods games. Coordination games, by their more obscure and less evocative nature, are left understudied, despite the fact that, as we argue, they offer great potential in understanding and treating cancer. In this paper we present the conditions under which coordination games between cancer cells evolve, we propose aspects of cancer that can be modeled as results of coordination games, and explore the ways through which coordination games of cancer can be exploited for therapy.

Dynamic Phenotypic Switching and Group Behavior Help Non-Small Cell Lung Cancer Cells Evade Chemotherapy

Drug resistance, a major challenge in cancer therapy, is typically attributed to mutations and genetic heterogeneity. Emerging evidence suggests that dynamic cellular interactions and group behavior also contribute to drug resistance. However, the underlying mechanisms remain poorly understood. Here, we present a new mathematical approach with game theoretical underpinnings that we developed to model real-time growth data of non-small cell lung cancer (NSCLC) cells and discern patterns in response to treatment with cisplatin. We show that the cisplatin-sensitive and cisplatin-tolerant NSCLC cells, when co-cultured in the absence or presence of the drug, display dynamic group behavior strategies. Tolerant cells exhibit a 'persister-like' behavior and are attenuated by sensitive cells; they also appear to 'educate' sensitive cells to evade chemotherapy. Further, tolerant cells can switch phenotypes to become sensitive, especially at low cisplatin concentrations. Finally, switching treatment from continuous to an intermittent regimen can attenuate the emergence of tolerant cells, suggesting that intermittent chemotherapy may improve outcomes in lung cancer.

Darwinian Approaches for Cancer Treatment: Benefits of Mathematical Modeling

One of the major problems of traditional anti-cancer treatments is that they lead to the emergence of treatment-resistant cells, which results in treatment failure. To avoid or delay this phenomenon, it is relevant to take into account the eco-evolutionary dynamics of tumors. Designing evolution-based treatment strategies may help overcoming the problem of drug resistance. In particular, a promising candidate is adaptive therapy, a containment strategy which adjusts treatment cycles to the evolution of the tumors in order to keep the population of treatment-resistant cells under control. Mathematical modeling is a crucial tool to understand the dynamics of cancer in response to treatments, and to make predictions about the outcomes of these treatments. In this review, we highlight the benefits of in silico modeling to design adaptive therapy strategies, and to assess whether they could effectively improve treatment outcomes. Specifically, we review how two main types of models (i.e., mathematical models based on Lotka-Volterra equations and agent-based models) have been used to model tumor dynamics in response to adaptive therapy. We give examples of the advances they permitted in the field of adaptive therapy and discuss about how these models can be integrated in experimental approaches and clinical trial design.

The Contribution of Evolutionary Game Theory to Understanding and Treating Cancer

Evolutionary game theory mathematically conceptualizes and analyzes biological interactions where one's fitness not only depends on one's own traits, but also on the traits of others. Typically, the individuals are not overtly rational and do not select, but rather inherit their traits. Cancer can be framed as such an evolutionary game, as it is composed of cells of heterogeneous types undergoing frequency-dependent selection. In this article, we first summarize existing works where evolutionary game theory has been employed in modeling cancer and improving its treatment. Some of these game-theoretic models suggest how one could anticipate and steer cancer's eco-evolutionary dynamics into states more desirable for the patient via evolutionary therapies. Such therapies offer great promise for increasing patient survival and decreasing drug toxicity, as demonstrated by some recent studies and clinical trials. We discuss clinical relevance of the existing game-theoretic models of cancer and its treatment, and opportunities for future applications. Moreover, we discuss the developments in cancer biology that are needed to better utilize the full potential of game-theoretic models. Ultimately, we demonstrate that viewing tumors with evolutionary game theory has medically useful implications that can inform and create a lockstep between empirical findings and mathematical modeling. We suggest that cancer progression is an evolutionary competition between different cell types and therefore needs to be viewed as an evolutionary game.

Evolutionary Dynamics of Treatment-Induced Resistance in Cancer Informs Understanding of Rapid Evolution in Natural Systems

Rapid evolution is ubiquitous in nature. We briefly review some of this quite broadly, particularly in the context of response to anthropogenic disturbances. Nowhere is this more evident, replicated and accessible to study than in cancer. Curiously cancer has been late - relative to fisheries, antibiotic resistance, pest management and evolution in human dominated landscapes - in recognizing the need for evolutionarily informed management strategies. The speed of evolution matters. Here, we employ game-theoretic modeling to compare time to progression with continuous maximum tolerable dose to that of adaptive therapy where treatment is discontinued when the population of cancer cells gets below half of its initial size and re-administered when the cancer cells recover, forming cycles with and without treatment. We show that the success of adaptive therapy relative to continuous maximum tolerable dose therapy is much higher if the population of cancer cells is defined by two cell types (sensitive vs. resistant in a polymorphic population). Additionally, the relative increase in time to progression increases with the speed of evolution. These results hold with and without a cost of resistance in cancer cells. On the other hand, treatment-induced resistance can be modeled as a quantitative trait in a monomorphic population of cancer cells. In that case, when evolution is rapid, there is no advantage to adaptive therapy. Initial responses to therapy are blunted by the cancer cells evolving too quickly. Our study emphasizes how cancer provides a unique system for studying rapid evolutionary changes within tumor ecosystems in response to human interventions; and allows us to contrast and compare this system to other human managed or dominated systems in nature.

Cooperative success in epithelial public goods games

Cancer cells obtain mutations which rely on the production of diffusible growth factors to confer a fitness benefit. These mutations can be considered cooperative, and studied as public goods games within the framework of evolutionary game theory. The population structure, benefit function and update rule all influence the evolutionary success of cooperators. We model the evolution of cooperation in epithelial cells using the Voronoi tessellation model. Unlike traditional evolutionary graph theory, this allows us to implement global updating, for which birth and death events are spatially decoupled. We compare, for a sigmoid benefit function, the conditions for cooperation to be favoured and/or beneficial for well-mixed and structured populations. We find that when population structure is combined with global updating, cooperation is more successful than if there were local updating or the population were well-mixed. Interestingly, the qualitative behaviour for the well-mixed population and the Voronoi tessellation model is remarkably similar, but the latter case requires significantly lower incentives to ensure cooperation.

IsoMaTrix: a framework to visualize the isoclines of matrix games and quantify uncertainty in structured populations

 

Evolutionary game theory describes frequency-dependent selection for fixed, heritable strategies in a population of competing individuals using a payoff matrix. We present a software package to aid in the construction, analysis, and visualization of three-strategy matrix games. The IsoMaTrix package computes the isoclines (lines of zero growth) of matrix games, and facilitates direct comparison of well-mixed dynamics to structured populations on a lattice grid. IsoMaTrix computes fixed points, phase flow, trajectories, (sub)velocities, and uncertainty quantification for stochastic effects in spatial matrix games. We describe a result obtained via IsoMaTrix's spatial games functionality, which shows that the timing of competitive release in a cancer model (under continuous treatment) critically depends on the initial spatial configuration of the tumor.

AVAILABILITY

The code is available at: https://github.com/mathonco/isomatrix.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Optimal control to reach eco-evolutionary stability in metastatic castrate-resistant prostate cancer

In the absence of curative therapies, treatment of metastatic castrate-resistant prostate cancer (mCRPC) using currently available drugs can be improved by integrating evolutionary principles that govern proliferation of resistant subpopulations into current treatment protocols. Here we develop what is coined as an 'evolutionary stable therapy', within the context of the mathematical model that has been used to inform the first adaptive therapy clinical trial of mCRPC. The objective of this therapy is to maintain a stable polymorphic tumor heterogeneity of sensitive and resistant cells to therapy in order to prolong treatment efficacy and progression free survival. Optimal control analysis shows that an increasing dose titration protocol, a very common clinical dosing process, can achieve tumor stabilization for a wide range of potential initial tumor compositions and volumes. Furthermore, larger tumor volumes may counter intuitively be more likely to be stabilized if sensitive cells dominate the tumor composition at time of initial treatment, suggesting a delay of initial treatment could prove beneficial. While it remains uncertain if metastatic disease in humans has the properties that allow it to be truly stabilized, the benefits of a dose titration protocol warrant additional pre-clinical and clinical investigations.

Integrating evolutionary dynamics into cancer therapy

Many effective drugs for metastatic and/or advanced-stage cancers have been developed over the past decade, although the evolution of resistance remains the major barrier to disease control or cure. In large, diverse populations such as the cells that compose metastatic cancers, the emergence of cells that are resistant or that can quickly develop resistance is virtually inevitable and most likely cannot be prevented. However, clinically significant resistance occurs only when the pre-existing resistant phenotypes are able to proliferate extensively, a process governed by eco-evolutionary dynamics. Attempts to disrupt the molecular mechanisms of resistance have generally been unsuccessful in clinical practice. In this Review, we focus on the Darwinian processes driving the eco-evolutionary dynamics of treatment-resistant cancer populations. We describe a variety of evolutionarily informed strategies designed to increase the probability of disease control or cure by anticipating and steering the evolutionary dynamics of acquired resistance.

Can Environmental Manipulation Help Suppress Cancer? Non-Linear Competition Among Tumor Cells in Periodically Changing Conditions

It has been shown that the tumor population growth dynamics in a periodically varying environment can drastically differ from the one in a fixed environment. Thus, the environment of a tumor can potentially be manipulated to suppress cancer progression. Diverse evolutionary processes play vital roles in cancer progression and accordingly, understanding the interplay between these processes is essential in optimizing the treatment strategy. Somatic evolution and genetic instability result in intra-tumor cell heterogeneity. Various models have been developed to analyze the interactions between different types of tumor cells. Here, models of density-dependent interaction between different types of tumor cells under fast periodical environmental changes are examined. It is illustrated that tumor population densities, which vary on a slow time scale, are affected by fast environmental variations. Finally, the intriguing density-dependent interactions in metastatic castration-resistant prostate cancer (mCRPC) in which the different types of tumor cells are defined with respect to the production of and dependence on testosterone are discussed.

Optimizing adaptive cancer therapy: dynamic programming and evolutionary game theory

Recent clinical trials have shown that adaptive drug therapies can be more efficient than a standard cancer treatment based on a continuous use of maximum tolerated doses (MTD). The adaptive therapy paradigm is not based on a preset schedule; instead, the doses are administered based on the current state of tumour. But the adaptive treatment policies examined so far have been largely ad hoc. We propose a method for systematically optimizing adaptive policies based on an evolutionary game theory model of cancer dynamics. Given a set of treatment objectives, we use the framework of dynamic programming to find the optimal treatment strategies. In particular, we optimize the total drug usage and time to recovery by solving a Hamilton-Jacobi-Bellman equation. We compare MTD-based treatment strategy with optimal adaptive treatment policies and show that the latter can significantly decrease the total amount of drugs prescribed while also increasing the fraction of initial tumour states from which the recovery is possible. We conclude that the use of optimal control theory to improve adaptive policies is a promising concept in cancer treatment and should be integrated into clinical trial design.

Macronutrient intakes and the lifespan-fecundity trade-off: a geometric framework agent-based model

Lifespan and fecundity, the main components in evolutionary fitness, are both strongly affected by nutritional state. Geometric framework of nutrition (GFN) experiments has shown that lifespan and fecundity are separated in nutrient space leading to a functional trade-off between the two traits. Here we develop a spatially explicit agent-based model (ABM) using the GFN to explore how ecological factors may cause selection on macronutrient appetites to optimally balance these life-history traits. We show that increasing the risk of extrinsic mortality favours intake of a mixture of nutrients that is associated with maximal fecundity at the expense of reduced longevity and that this result is robust across spatial and nutritional environments. These model behaviours are consistent with what has been observed in studies that quantify changes in life history in response to environmental manipulations. Previous GFN-derived ABMs have treated fitness as a single value. This is the first such model to instead decompose fitness into its primary component traits, longevity and fecundity, allowing evolutionary fitness to be an emergent property of the two. Our model demonstrates that selection on macronutrient appetites may affect life-history trade-offs and makes predictions that can be directly tested in artificial selection experiments.

Modeling Tumor Evolutionary Dynamics to Predict Clinical Outcomes for Patients with Metastatic Colorectal Cancer: A Retrospective Analysis

Over 50% of colorectal cancer patients develop resistance after a transient response to therapy. Understanding tumor resistance from an evolutionary perspective leads to better predictions of treatment outcomes. The objectives of this study were to develop a computational framework to analyze tumor longitudinal measurements and recapitulate the individual evolutionary dynamics in metastatic colorectal cancer (mCRC) patients. A stochastic modeling framework was developed to depict the whole spectrum of tumor evolution prior to diagnosis and during and after therapy. The evolutionary model was optimized using a nonlinear mixed effect (NLME) method based on the longitudinal measurements of liver metastatic lesions from 599 mCRC patients. The deterministic limits in the NLME model were applied to optimize the stochastic model for each patient. Cox proportional hazards models coupled with the least absolute shrinkage and selection operator (LASSO) algorithm were applied to predict patients' progression-free survival (PFS) and overall survival (OS). The stochastic evolutionary model well described the longitudinal profiles of tumor sizes. The evolutionary parameters optimized for each patient indicated substantial interpatient variability. The number of resistant subclones at diagnosis was found to be a significant predictor to survival, and the hazard ratios with 95% CI were 1.09 (0.79-1.49) and 1.54 (1.01-2.34) for patients with three or more resistant subclones. Coupled with several patient characteristics, evolutionary parameters strongly predict patients' PFS and OS. A stochastic computational framework was successfully developed to recapitulate individual patient evolutionary dynamics, which could predict clinical survival outcomes in mCRC patients. SIGNIFICANCE: A data analysis framework depicts the individual evolutionary dynamics of mCRC patients and can be generalized to project patient survival outcomes.

Towards Multidrug Adaptive Therapy

A new ecologically inspired paradigm in cancer treatment known as "adaptive therapy" capitalizes on competitive interactions between drug-sensitive and drug-resistant subclones. The goal of adaptive therapy is to maintain a controllable stable tumor burden by allowing a significant population of treatment-sensitive cells to survive. These, in turn, suppress proliferation of the less-fit resistant populations. However, there remain several open challenges in designing adaptive therapies, particularly in extending these therapeutic concepts to multiple treatments. We present a cancer treatment case study (metastatic castrate-resistant prostate cancer) as a point of departure to illustrate three novel concepts to aid the design of multidrug adaptive therapies. First, frequency-dependent "cycles" of tumor evolution can trap tumor evolution in a periodic, controllable loop. Second, the availability and selection of treatments may limit the evolutionary "absorbing region" reachable by the tumor. Third, the velocity of evolution significantly influences the optimal timing of drug sequences. These three conceptual advances provide a path forward for multidrug adaptive therapy. SIGNIFICANCE: Driving tumor evolution into periodic, repeatable treatment cycles provides a path forward for multidrug adaptive therapy.

Multidrug Cancer Therapy in Metastatic Castrate-Resistant Prostate Cancer: An Evolution-Based Strategy

PURPOSE

Integration of evolutionary dynamics into systemic therapy for metastatic cancers can prolong tumor control compared with standard maximum tolerated dose (MTD) strategies. Prior investigations have focused on monotherapy, but many clinical cancer treatments combine two or more drugs. Optimizing the evolutionary dynamics in multidrug therapy is challenging because of the complex cellular interactions and the large parameter space of potential variations in drugs, doses, and treatment schedules. However, multidrug therapy also represents an opportunity to further improve outcomes using evolution-based strategies.

EXPERIMENTAL DESIGN

We examine evolution-based strategies for two-drug therapy and identify an approach that divides the treatment drugs into primary and secondary roles. The primary drug has the greatest efficacy and/or lowest toxicity. The secondary drug is applied solely to reduce the resistant population to the primary drug.

RESULTS

Simulations from the mathematical model demonstrate that the primary-secondary approach increases time to progression (TTP) compared with conventional strategies in which drugs are administered without regard to evolutionary dynamics. We apply our model to an ongoing adaptive therapy clinical trial of evolution-based administration of abiraterone to treat metastatic castrate-resistant prostate cancer. Model simulations, parameterized with data from individual patients who progressed, demonstrate that strategic application of docetaxel during abiraterone therapy would have significantly increased their TTP.

CONCLUSIONS

Mathematical models can integrate evolutionary dynamics into multidrug cancer clinical trials. This has the potential to improve outcomes and to develop clinical trials in which these mathematical models are also used to estimate the mechanism(s) of treatment failure and explore alternative strategies to improve outcomes in future trials.

Including Blood Vasculature into a Game-Theoretic Model of Cancer Dynamics

For cancer, we develop a 2-D agent-based continuous-space game-theoretical model that considers cancer cells’ proximity to a blood vessel. Based on castrate resistant metastatic prostate cancer (mCRPC), the model considers the density and frequency (eco-evolutionary) dynamics of three cancer cell types: those that require exogenous testosterone ( T + ), those producing testosterone ( T P ), and those independent of testosterone ( T - ). We model proximity to a blood vessel by imagining four zones around the vessel. Zone 0 is the blood vessel. As rings, zones 1–3 are successively farther from the blood vessel and have successively lower carrying capacities. Zone 4 represents the space too far from the blood vessel and too poor in nutrients for cancer cell proliferation. Within the other three zones that are closer to the blood vessel, the cells’ proliferation probabilities are determined by zone-specific payoff matrices. We analyzed how zone width, dispersal, interactions across zone boundaries, and blood vessel dynamics influence the eco-evolutionary dynamics of cell types within zones and across the entire cancer cell population. At equilibrium, zone 3’s composition deviates from its evolutionary stable strategy (ESS) towards that of zone 2. Zone 2 sees deviations from its ESS because of dispersal from zones 1 and 3; however, its composition begins to resemble zone 1’s more so than zone 3’s. Frequency-dependent interactions between cells across zone boundaries have little effect on zone 2’s and zone 3’s composition but have decisive effects on zone 1. The composition of zone 1 diverges dramatically from both its own ESS, but also that of zone 2. That is because T + cells (highest frequency in zone 1) benefit from interacting with T P cells (highest frequency in zone 2). Zone 1 T + cells interacting with cells in zone 2 experience a higher likelihood of encountering a T P cell than when restricted to their own zone. As expected, increasing the width of zones decreases these impacts of cross-boundary dispersal and interactions. Increasing zone widths increases the persistence likelihood of the cancer subpopulation in the face of blood vessel dynamics, where the vessel may die or become occluded resulting in the “birth” of another blood vessel elsewhere in the space. With small zone widths, the cancer cell subpopulations cannot persist. With large zone widths, blood vessel dynamics create cancer cell subpopulations that resemble the ESS of zone 3 as the larger area of zone 3 and its contribution to cells within the necrotic zone 4 mean that zones 3 and 4 provide the likeliest colonizers for the new blood vessel. In conclusion, our model provides an alternative modeling approach for considering density-dependent, frequency-dependent, and dispersal dynamics into cancer models with spatial gradients around blood vessels. Additionally, our model can consider the occurrence of circulating tumor cells (cells that disperse into the blood vessel from zone 1) and the presence of live cancer cells within the necrotic regions of a tumor.

Cooperation among cancer cells: applying game theory to cancer

Cell cooperation promotes many of the hallmarks of cancer via the secretion of diffusible factors that can affect cancer cells or stromal cells in the tumour microenvironment. This cooperation cannot be explained simply as the collective action of cells for the benefit of the tumour because non-cooperative subclones can constantly invade and free-ride on the diffusible factors produced by the cooperative cells. A full understanding of cooperation among the cells of a tumour requires methods and concepts from evolutionary game theory, which has been used successfully in other areas of biology to understand similar problems but has been underutilized in cancer research. Game theory can provide insights into the stability of cooperation among cells in a tumour and into the design of potentially evolution-proof therapies that disrupt this cooperation.

Nonlinear adaptive control of competitive release and chemotherapeutic resistance

We use a three-component replicator system with healthy cells, sensitive cells, and resistant cells, with a prisoner's dilemma payoff matrix from evolutionary game theory, to model and control the nonlinear dynamical system governing the ecological mechanism of competitive release by which tumors develop chemotherapeutic resistance. The control method we describe is based on nonlinear trajectory design and energy transfer methods first introduced in the orbital mechanics literature for Hamiltonian systems. For continuous therapy, the basin boundaries of attraction associated with the chemo-sensitive population and the chemo-resistant population for increasing values of chemo-concentrations have an intertwined spiral structure with extreme sensitivity to changes in chemo-concentration level as well as sensitivity with respect to resistant mutations. For time-dependent therapies, we introduce an orbit transfer method to construct continuous families of periodic (closed) orbits by switching the chemo-dose at carefully chosen times and appropriate levels to design schedules that are superior to both maximum tolerated dose (MTD) schedules and low-dose metronomic (LDM) schedules, both of which ultimately lead to fixation of sensitive cells or resistant cells. By keeping the three subpopulations of cells in competition with each other indefinitely, we avoid fixation of the cancer cell population and regrowth of a resistant tumor. The method can be viewed as a way to dynamically shape the average population fitness landscape of a tumor to steer the chemotherapeutic response curve. We show that the method is remarkably insensitive to initial conditions and small changes in chemo-dosages, an important criterion for turning the method into an actionable strategy.

Structured environments fundamentally alter dynamics and stability of ecological communities

The dynamics and stability of ecological communities are intimately linked with the specific interactions-like cooperation or predation-between constituent species. In microbial communities, like those found in soils or the mammalian gut, physical anisotropies produced by fluid flow and chemical gradients impact community structure and ecological dynamics, even in structurally isotropic environments. Although natural communities existing in physically unstructured environments are rare, the role of environmental structure in determining community dynamics and stability remains poorly studied. To address this gap, we used modified Lotka-Volterra simulations of competitive microbial communities to characterize the effects of surface structure on community dynamics. We find that environmental structure has profound effects on communities, in a manner dependent on the specific pattern of interactions between community members. For two mutually competing species, eventual extinction of one competitor is effectively guaranteed in isotropic environments. However, addition of environmental structure enables long-term coexistence of both species via local "pinning" of competition interfaces, even when one species has a significant competitive advantage. In contrast, while three species competing in an intransitive loop (as in a game of rock-paper-scissors) coexist stably in isotropic environments, structural anisotropy disrupts the spatial patterns on which coexistence depends, causing chaotic population fluctuations and subsequent extinction cascades. These results indicate that the stability of microbial communities strongly depends on the structural environment in which they reside. Therefore, a more complete ecological understanding, including effective manipulation and interventions in natural communities of interest, must account for the physical structure of the environment.

Moving the systemic evolutionary approach to cancer forward: Therapeutic implications

We have previously presented a new Systemic Evolutionary Theory of Cancer (SETOC) based on the failure of proper endosymbiosis in eukaryotic cells. Here, we propose that the progressive uncoupling of two endosymbiotic subsystems (information and energy) inside the cell, as a consequence of long-term injuries, gives rise to alterations (i) in tissue interactions and (ii) in cell organization. In the first case, we argue that the impairment of both the coherent state and the synergy between intercellular communications underpins the onset of tissue dysplasia, that usually evolves towards cancer development. In the second case, we suggest that the rupture of endosymbiosis drives a sort of cell regression towards a protist-like entity represented by the concept of "de-emergence" postulated in our systemic evolutionary approach to carcinogenesis. This conceptual association of the cancer cell with a protist-like organism could support the development of novel cancer therapeutic approaches. To this end, we propose a paradigm shift in cancer pharmacology since: i) our knowledge of cancer pathophysiology as a complex system is insufficient, despite a vast knowledge of molecular mechanisms underlying cancer; ii) current cancer pharmacology deals only with microvariables (e.g. gene or protein targets), which do not account for the integrated pathophysiology of cancer, rather than with macrovariables (e.g. pH, membrane potential, electromagnetic fields, cell communications and so on) and mesovariables (between micro and macro), such as the interaction between various cellular components including cellular organelles. This paradigm shift should allow cancer pharmacology to move forward from molecular treatments (focusing on single targets) to modular treatments that consider cancer-related processes (i.e. inflammation, coagulation, etc.) or even to a sort of ecosystemic treatment addressing the whole functioning of the "cancer ecosystem". Examples of ecosystems treatment may be natural plant derivatives that act synergistically or pulsed electromagnetic fields which can act on particular biological processes in cancer cells. In addition, we need different working theoretical models on which to base new anticancer pharmacological approaches. Finally, we examine what value our systemic evolutionary approach could add to cancer treatments, in particular in liver cancer as a paradigm for developing potential applications.

Eco-evolutionary causes and consequences of temporal changes in intratumoural blood flow

Temporal changes in blood flow are commonly observed in malignant tumours, but the evolutionary causes and consequences are rarely considered. We propose that stochastic temporal variations in blood flow and microenvironmental conditions arise from the eco-evolutionary dynamics of tumour angiogenesis in which cancer cells, as individual units of selection, can influence and respond only to local environmental conditions. This leads to new vessels arising from the closest available vascular structure regardless of the size or capacity of this parental vessel. These dynamics produce unstable vascular networks with unpredictable spatial and temporal variations in blood flow and microenvironmental conditions. Adaptations of evolving populations to temporally varying environments in nature include increased diversity, greater motility and invasiveness, and highly plastic phenotypes, allowing for broad metabolic adaptability and rapid shifts to high rates of proliferation and profound quiescence. These adaptive strategies, when adopted in cancer cells, promote many commonly observed phenotypic properties including those found in the stem phenotype and in epithelial-to-mesenchymal transition. Temporal variations in intratumoural blood flow, which occur through the promotion of cancer cell phenotypes that facilitate both metastatic spread and resistance to therapy, may have substantial clinical consequences.

Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer

Abiraterone treats metastatic castrate-resistant prostate cancer by inhibiting CYP17A, an enzyme for testosterone auto-production. With standard dosing, evolution of resistance with treatment failure (radiographic progression) occurs at a median of ~16.5 months. We hypothesize time to progression (TTP) could be increased by integrating evolutionary dynamics into therapy. We developed an evolutionary game theory model using Lotka–Volterra equations with three competing cancer “species”: androgen dependent, androgen producing, and androgen independent. Simulations with standard abiraterone dosing demonstrate strong selection for androgen-independent cells and rapid treatment failure. Adaptive therapy, using patient-specific tumor dynamics to inform on/off treatment cycles, suppresses proliferation of androgen-independent cells and lowers cumulative drug dose. In a pilot clinical trial, 10 of 11 patients maintained stable oscillations of tumor burdens; median TTP is at least 27 months with reduced cumulative drug use of 47% of standard dosing. The outcomes show significant improvement over published studies and a contemporaneous population.

Evolution of resistance is a common cause of cancer treatment failure and tumor progression. Here, the authors present a method for integrating evolutionary principles based on adaptive therapy into abiraterone therapy for metastatic castrate-resistant prostate cancer and show the positive results of an interim analysis of a trial cohort.