The evolution of carrying capacity in constrained and expanding tumour cell populations

Modeling tumors as complex ecosystems

Many cancers resist therapeutic intervention. This is fundamentally related to intratumor heterogeneity: multiple cell populations, each with different phenotypic signatures, coexist within a tumor and its metastases. Like species in an ecosystem, cancer populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity or predict its consequences. Here, we propose that the generalized Lotka-Volterra model (GLV), a standard tool to describe species-rich ecological communities, provides a suitable framework to model the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties provide a new understanding of the disease. We discuss potential extensions of the model and their application to phenotypic plasticity, cancer-immune interactions, and metastatic growth. Our work outlines a set of questions and a road map for further research in cancer ecology.

Modeling tumors as species-rich ecological communities

Many advanced cancers resist therapeutic intervention. This process is fundamentally related to intra-tumor heterogeneity: multiple cell populations, each with different mutational and phenotypic signatures, coexist within a tumor and its metastatic nodes. Like species in an ecosystem, many cancer cell populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity nor are able to predict its consequences. Here we propose that the Generalized Lotka-Volterra model (GLV), a standard tool to describe complex, species-rich ecological communities, provides a suitable framework to describe the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties, such as outgrowth and multistability, provide a new understanding of the disease. Additionally, we discuss potential extensions of the model and their application to three active areas of cancer research, namely phenotypic plasticity, the cancer-immune interplay and the resistance of metastatic tumors to treatment. Our work outlines a set of questions and a tentative road map for further research in cancer ecology.

Mathematical Modeling of Clonal Interference by Density-Dependent Selection in Heterogeneous Cancer Cell Lines

Many cancer cell lines are aneuploid and heterogeneous, with multiple karyotypes co-existing within the same cell line. Karyotype heterogeneity has been shown to manifest phenotypically, thus affecting how cells respond to drugs or to minor differences in culture media. Knowing how to interpret karyotype heterogeneity phenotypically would give insights into cellular phenotypes before they unfold temporally. Here, we re-analyzed single cell RNA (scRNA) and scDNA sequencing data from eight stomach cancer cell lines by placing gene expression programs into a phenotypic context. Using live cell imaging, we quantified differences in the growth rate and contact inhibition between the eight cell lines and used these differences to prioritize the transcriptomic biomarkers of the growth rate and carrying capacity. Using these biomarkers, we found significant differences in the predicted growth rate or carrying capacity between multiple karyotypes detected within the same cell line. We used these predictions to simulate how the clonal composition of a cell line would change depending on density conditions during in-vitro experiments. Once validated, these models can aid in the design of experiments that steer evolution with density-dependent selection.

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.

Different costs of therapeutic resistance in cancer: Short- and long-term impact of population heterogeneity

Therapeutic resistance continues to undercut long-term success of many promising cancer treatments. At times, development of therapeutic resistance can come at a fitness cost for the cancer cell population, which could potentially be leveraged to the patient's advantage. A mathematical formulation of such a situation was proposed by Pressley et al. (2020), who discussed two scenarios, namely, when developing therapeutic resistance can come at a cost to proliferative capacity (such as when a drug targets a growth receptor), or to the total tumor carrying capacity (such as when a drug targets neovascularization). Here we expand the analysis of the two models and evaluate both short- and long-term dynamics of a population heterogeneous with respect to resistance. We analyze the four initial distributions with respect to resistance at the time of treatment initiation: uniform, bell-shaped, exponential, and U-shaped. We show that final population composition is invariant to the initial distribution, with a single clone eventually dominating within the population; the value of the resistance parameter of the final clone depends on other system parameters but not on the initial distribution. Transitional behaviors, however, which may have more significant implications for immediate treatment decisions, depend critically on the initial distribution. Furthermore, we show that depending on the mechanism for the cost of resistance (i.e., proliferation vs carrying capacity), increase in natural cell death rate has opposite effects, with higher natural death rate selecting for less resistant cell clones in the long term for proliferation-dependent model, and selecting for more resistant cell clones for carrying capacity-dependent model, a prediction that may have implications for combination therapy with cytotoxic agents. We conclude with a discussion of strengths and limitations of using modeling for understanding treatment trajectory, as well as the promise of model-informed evolutionary steering for improved long-term therapeutic outcomes.

The impact of the spatial heterogeneity of resistant cells and fibroblasts on treatment response

A long-standing practice in the treatment of cancer is that of hitting hard with the maximum tolerated dose to eradicate tumors. This continuous therapy, however, selects for resistant cells, leading to the failure of the treatment. A different type of treatment strategy, adaptive therapy, has recently been shown to have a degree of success in both preclinical xenograft experiments and clinical trials. Adaptive therapy is used to maintain a tumor's volume by exploiting the competition between drug-sensitive and drug-resistant cells with minimum effective drug doses or timed drug holidays. To further understand the role of competition in the outcomes of adaptive therapy, we developed a 2D on-lattice agent-based model. Our simulations show that the superiority of the adaptive strategy over continuous therapy depends on the local competition shaped by the spatial distribution of resistant cells. Intratumor competition can also be affected by fibroblasts, which produce microenvironmental factors that promote cancer cell growth. To this end, we simulated the impact of different fibroblast distributions on treatment outcomes. As a proof of principle, we focused on five types of distribution of fibroblasts characterized by different locations, shapes, and orientations of the fibroblast region with respect to the resistant cells. Our simulation shows that the spatial architecture of fibroblasts modulates tumor progression in both continuous and adaptive therapy. Finally, as a proof of concept, we simulated the outcomes of adaptive therapy of a virtual patient with four metastatic sites composed of different spatial distributions of fibroblasts and drug-resistant cell populations. Our simulation highlights the importance of undetected metastatic lesions on adaptive therapy outcomes.

The ecology of cancer differentiation therapy

A promising, yet still under development approach to cancer treatment is based on the idea of differentiation therapy (DTH). Most tumours are characterized by poorly differentiated cell populations exhibiting a marked loss of traits associated to communication and tissue homeostasis. DTH has been suggested as an alternative (or complement) to cytotoxic-based approaches, and has proven successful in some specific types of cancer such as acute promyelocytic leukemia (APL). While novel drugs favouring the activation of differentiation therapies are being tested, several open problems emerge in relation to its effectiveness on solid tumors. Here we present a mathematical framework to DTH based on a well-known ecological model used to describe habitat loss. The models presented here account for some of the observed clinical and in vitro outcomes of DTH, providing relevant insight into potential therapy design. Furthermore, the same ecological approach is tested in a hierarchical model that accounts for cancer stem cells, highlighting the role of niche specificity in CSC therapy resistance. We show that the lessons learnt from metapopulation ecology can help guide future developments and potential difficulties of DTH.

Eco-oncology: Applying ecological principles to understand and manage cancer

Cancer is a disease of single cells that expresses itself at the population level. The striking similarities between initiation and growth of tumors and dynamics of biological populations, and between metastasis and ecological invasion and community dynamics suggest that oncology can benefit from an ecological perspective to improve our understanding of cancer biology. Tumors can be viewed as complex, adaptive, and evolving systems as they are spatially and temporally heterogeneous, continually interacting with each other and with the microenvironment and evolving to increase the fitness of the cancer cells. We argue that an eco-evolutionary perspective is essential to understand cancer biology better. Furthermore, we suggest that ecologically informed therapeutic approaches that combine standard of care treatments with strategies aimed at decreasing the evolutionary potential and fitness of neoplastic cells, such as disrupting cell-to-cell communication and cooperation, and preventing successful colonization of distant organs by migrating cancer cells, may be effective in managing cancer as a chronic condition.

Adaptive dynamics of unstable cancer populations: The canonical equation

In most instances of tumour development, genetic instability plays a role in allowing cancer cell populations to respond to selection barriers, such as physical constraints or immune responses, and rapidly adapt to an always changing environment. Modelling instability is a nontrivial task, since by definition evolving instability leads to changes in the underlying landscape. In this article, we explore mathematically a simple version of unstable tumour progression using the formalism of adaptive dynamics (AD) where selection and mutation are explicitly coupled. Using a set of basic fitness landscapes, the so-called canonical equation for the evolution of genetic instability on a minimal scenario associated with a population of unstable cells is derived. We obtain explicit expressions for the evolution of mutation probabilities, and the implications of the model on further experimental studies and potential mutagenic therapies are discussed.

Bayesian inference of agent-based models: a tool for studying kidney branching morphogenesis

The adult mammalian kidney has a complex, highly-branched collecting duct epithelium that arises as a ureteric bud sidebranch from an epithelial tube known as the nephric duct. Subsequent branching of the ureteric bud to form the collecting duct tree is regulated by subcellular interactions between the epithelium and a population of mesenchymal cells that surround the tips of outgrowing branches. The mesenchymal cells produce glial cell-line derived neurotrophic factor (GDNF), that binds with RET receptors on the surface of the epithelial cells to stimulate several subcellular pathways in the epithelium. Such interactions are known to be a prerequisite for normal branching development, although competing theories exist for their role in morphogenesis. Here we introduce the first agent-based model of ex vivo kidney uretic branching. Through comparison with experimental data, we show that growth factor-regulated growth mechanisms can explain early epithelial cell branching, but only if epithelial cell division depends in a switch-like way on the local growth factor concentration; cell division occurring only if the driving growth factor level exceeds a threshold. We also show how a recently-developed method, "Approximate Approximate Bayesian Computation", can be used to infer key model parameters, and reveal the dependency between the parameters controlling a growth factor-dependent growth switch. These results are consistent with a requirement for signals controlling proliferation and chemotaxis, both of which are previously identified roles for GDNF.

Competition and niche construction in a model of cancer metastasis

Niche construction theory states that not only does the environment act on populations to generate Darwinian selection, but organisms reciprocally modify the environment and the sources of natural selection. Cancer cells participate in niche construction as they alter their microenvironments and create pre-metastatic niches; in fact, metastasis is a product of niche construction. Here, we present a mathematical model of niche construction and metastasis. Our model contains producers, which pay a cost to contribute to niche construction that benefits all tumor cells, and cheaters, which reap the benefits without paying the cost. We derive expressions for the conditions necessary for metastasis, showing that the establishment of a mutant lineage that promotes metastasis depends on niche construction specificity and strength of interclonal competition. We identify a tension between the arrival and invasion of metastasis-promoting mutants, where tumors composed only of cheaters remain small but are susceptible to invasion whereas larger tumors containing producers may be unable to facilitate metastasis depending on the level of niche construction specificity. Our results indicate that even if metastatic subclones arise through mutation, metastasis may be hindered by interclonal competition, providing a potential explanation for recent surprising findings that most metastases are derived from early mutants in primary tumors.

Incorporating drug delivery into an imaging-driven, mechanics-coupled reaction diffusion model for predicting the response of breast cancer to neoadjuvant chemotherapy: theory and preliminary clinical results

Clinical methods for assessing tumor response to therapy are largely rudimentary, monitoring only temporal changes in tumor size. Our goal is to predict the response of breast tumors to therapy using a mathematical model that utilizes magnetic resonance imaging (MRI) data obtained non-invasively from individual patients. We extended a previously established, mechanically coupled, reaction-diffusion model for predicting tumor response initialized with patient-specific diffusion weighted MRI (DW-MRI) data by including the effects of chemotherapy drug delivery, which is estimated using dynamic contrast-enhanced (DCE-) MRI data. The extended, drug incorporated, model is initialized using patient-specific DW-MRI and DCE-MRI data. Data sets from five breast cancer patients were used-obtained before, after one cycle, and at mid-point of neoadjuvant chemotherapy. The DCE-MRI data was used to estimate spatiotemporal variations in tumor perfusion with the extended Kety-Tofts model. The physiological parameters derived from DCE-MRI were used to model changes in delivery of therapy drugs within the tumor for incorporation in the extended model. We simulated the original model and the extended model in both 2D and 3D and compare the results for this five-patient cohort. Preliminary results show reductions in the error of model predicted tumor cellularity and size compared to the experimentally-measured results for the third MRI scan when therapy was incorporated. Comparing the two models for agreement between the predicted total cellularity and the calculated total cellularity (from the DW-MRI data) reveals an increased concordance correlation coefficient from 0.81 to 0.98 for the 2D analysis and 0.85 to 0.99 for the 3D analysis (p  <  0.01 for each) when the extended model was used in place of the original model. This study demonstrates the plausibility of using DCE-MRI data as a means to estimate drug delivery on a patient-specific basis in predictive models and represents a step toward the goal of achieving individualized prediction of tumor response to therapy.

Metabolism is the tie: The Bertalanffy-type cancer growth model as common denominator of various modelling approaches

Cancer or tumour growth has been addressed from a variety of mathematical modelling perspectives in the past. Examples are single variable growth models, reaction diffusion models, compartment models, individual cell-based models, clonal competition models, to name only a few. In this paper, we show that the so called Bertalanffy-type growth model is a macroscopic model variant that can be conceived as an optimal condensed modelling approach that to a high degree preserves complexity with respect to the aforementioned more complex modelling variants. The derivation of the Bertalanffy-type model is crucially based on features of metabolism. Therefore, this model contains a shape parameter that can be interpreted as a resource utilisation efficiency. This shape parameter reflects features that are usually captured in much more complex models. To be specific, the shape parameter is related to morphological structures of tumours, which in turn depend on metabolic conditions. We, furthermore, show that a single variable variant of the Bertalanffy-type model can straightforwardly be extended to a multiclonal competition model. Since competition is crucially based on available shared or clone-specific resources, the metabolism-based approach is an obvious candidate to capture clonal competition. Depending on the specific context, metabolic reprogramming or other oncogene driven changes either lead to a suppression of cancer cells or to an improved competition resulting in outgrowth of tumours. The parametrisation of the Bertalanffy-type growth model allows to account for this observed variety of cancer characteristics. The shape parameter, conceived as a classifier for healthy and oncogenic phenotypes, supplies a link to survival and evolutionary stability concepts discussed in demographic studies, such as opportunistic versus equilibrium strategies.

Extinction rates in tumour public goods games

Cancer evolution and progression are shaped by cellular interactions and Darwinian selection. Evolutionary game theory incorporates both of these principles, and has been proposed as a framework to understand tumour cell population dynamics. A cornerstone of evolutionary dynamics is the replicator equation, which describes changes in the relative abundance of different cell types, and is able to predict evolutionary equilibria. Typically, the replicator equation focuses on differences in relative fitness. We here show that this framework might not be sufficient under all circumstances, as it neglects important aspects of population growth. Standard replicator dynamics might miss critical differences in the time it takes to reach an equilibrium, as this time also depends on cellular turnover in growing but bounded populations. As the system reaches a stable manifold, the time to reach equilibrium depends on cellular death and birth rates. These rates shape the time scales, in particular, in coevolutionary dynamics of growth factor producers and free-riders. Replicator dynamics might be an appropriate framework only when birth and death rates are of similar magnitude. Otherwise, population growth effects cannot be neglected when predicting the time to reach an equilibrium, and cell-type-specific rates have to be accounted for explicitly.

Fighting Cancer with Mathematics and Viruses

After decades of research, oncolytic virotherapy has recently advanced to clinical application, and currently a multitude of novel agents and combination treatments are being evaluated for cancer therapy. Oncolytic agents preferentially replicate in tumor cells, inducing tumor cell lysis and complex antitumor effects, such as innate and adaptive immune responses and the destruction of tumor vasculature. With the availability of different vector platforms and the potential of both genetic engineering and combination regimens to enhance particular aspects of safety and efficacy, the identification of optimal treatments for patient subpopulations or even individual patients becomes a top priority. Mathematical modeling can provide support in this arena by making use of experimental and clinical data to generate hypotheses about the mechanisms underlying complex biology and, ultimately, predict optimal treatment protocols. Increasingly complex models can be applied to account for therapeutically relevant parameters such as components of the immune system. In this review, we describe current developments in oncolytic virotherapy and mathematical modeling to discuss the benefit of integrating different modeling approaches into biological and clinical experimentation. Conclusively, we propose a mutual combination of these research fields to increase the value of the preclinical development and the therapeutic efficacy of the resulting treatments.

Predicting Patient-Specific Radiotherapy Protocols Based on Mathematical Model Choice for Proliferation Saturation Index

Radiation is commonly used in cancer treatment. Over 50% of all cancer patients will undergo radiotherapy (RT) as part of cancer care. Scientific advances in RT have primarily focused on the physical characteristics of treatment including beam quality and delivery. Only recently have inroads been made into utilizing tumor biology and radiobiology to design more appropriate RT protocols. Tumors are composites of proliferating and growth-arrested cells, and overall response depends on their respective proportions at irradiation. Prokopiou et al. (Radiat Oncol 10:159, 2015) developed the concept of the proliferation saturation index (PSI) to augment the clinical decision process associated with RT. This framework is based on the application of the logistic equation to pre-treatment imaging data in order to estimate a patient-specific tumor carrying capacity, which is then used to recommend a specific RT protocol. It is unclear, however, how dependent clinical recommendations are on the underlying tumor growth law. We discuss a PSI framework with a generalized logistic equation that can capture kinetics of different well-known growth laws including logistic and Gompertzian growth. Estimation of model parameters on the basis of clinical data revealed that the generalized logistic model can describe data equally well for a wide range of the generalized logistic exponent value. Clinical recommendations based on the calculated PSI, however, are strongly dependent on the specific growth law assumed. Our analysis suggests that the PSI framework may best be utilized in clinical practice when the underlying tumor growth law is known, or when sufficiently many tumor growth models suggest similar fractionation protocols.

Effects of size refuge specificity on a predator-prey model

Predation is a major cause of early-stage mortality for prey individuals, which are often forced to use refuges in order to reduce the risk of being consumed. The ability of certain genotypes in a prey population to reach a size refuge from predation may contribute significantly to the preservation of community diversity. We investigate how the specificity of this behavior affects the evolution of a given population by using a modified Lotka-Volterra model, in which the proportion of each genotype available for predation consists of two components: an intrinsic part and a combination from all genotypes present in the population. The trade-off of these components is characterized by a specificity parameter. From the viewpoint of population dynamics, we observe that the ability of the mutant to invade the resident population strongly depends on the values of this parameter. Finally, we describe the possible evolutionary outcomes, analytically and numerically.

Stochastic simulation tools and continuum models for describing two-dimensional collective cell spreading with universal growth functions

Two-dimensional collective cell migration assays are used to study cancer and tissue repair. These assays involve combined cell migration and cell proliferation processes, both of which are modulated by cell-to-cell crowding. Previous discrete models of collective cell migration assays involve a nearest-neighbour proliferation mechanism where crowding effects are incorporated by aborting potential proliferation events if the randomly chosen target site is occupied. There are two limitations of this traditional approach: (i) it seems unreasonable to abort a potential proliferation event based on the occupancy of a single, randomly chosen target site; and, (ii) the continuum limit description of this mechanism leads to the standard logistic growth function, but some experimental evidence suggests that cells do not always proliferate logistically. Motivated by these observations, we introduce a generalised proliferation mechanism which allows non-nearest neighbour proliferation events to take place over a template of [Formula: see text] concentric rings of lattice sites. Further, the decision to abort potential proliferation events is made using a crowding function, f(C), which accounts for the density of agents within a group of sites rather than dealing with the occupancy of a single randomly chosen site. Analysing the continuum limit description of the stochastic model shows that the standard logistic source term, [Formula: see text], where λ is the proliferation rate, is generalised to a universal growth function, [Formula: see text]. Comparing the solution of the continuum description with averaged simulation data indicates that the continuum model performs well for many choices of f(C) and r. For nonlinear f(C), the quality of the continuum-discrete match increases with r.

Environmentally induced, occupational diseases with emphasis on chronic kidney disease of multifactorial origin affecting tropical countries

BACKGROUND

Environmentally induced, occupational diseases are increasing worldwide, especially in rural agricultural communities. Poverty-associated malnutrition, environmental hazards and pollution, and lack of access to clean water, safe sanitation, and modern healthcare facilities are often associated with these chronic illnesses.

METHOD

The authors systematically reviewed occupational public health issues that have been related to the environment. General interpretations of results were included as per the guidelines of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses. Pertinent publications from research databases were reviewed on (A) the risk-benefits, (B) the prevalence of risk factors for various diseases, (C) the benefits of not ignoring the risk factors (i.e., broader evidence), and (D) the risks, effects, and outcomes of different types of interventions. The authors used chronic kidney disease of multifactorial origin (CKDmfo) as an example to explore the theme. Emphasis was given to the regions with emerging economies and developing countries located in the vicinity of the equator.

FINDINGS

Geographical, socio-economic and aetiological similarities exist for many chronic non-communicable diseases that are affecting tropical countries around the equator. The authors identified manufacturing, mining, and agriculture as the biggest polluters of the environment. In addition, deforestation and associated soil erosion, overuse of agrochemicals, and irresponsible factory discharge (e.g., chemicals and paint, from rubber and textile factories, etc.), all contribute to pollution. To decrease the escalating incidences of environmentally induced diseases, governments should work proactively to protect the environment, especially watersheds, and take steps to minimise harmful occupational exposures and strictly enforce environmental regulations.

CONCLUSION

Creating public awareness of environmental issues and their relationship to public health is essential. This includes regular monitoring and periodic publication of the quality of water, air and soil; preventing deforestation and man-made soil erosion, increasing forest and ground cover, preventing occupational injuries, judicious and safe use of agrochemicals, sustainable agriculture and development programs, and implementing legislation to protect and conserve water heriage and the environment. These actions are essential both for a healthier environment and for the health of the people who live in that environment. Such measures would also decrease public health threats from such, including global-warming-related erratic environmental changes and the occurrence and the spread of non-communicable diseases, such as CKDmfo.