A mathematical model for the dynamics of cancer hepatocytes under therapeutic actions

Evolutionary dynamics of competing phenotype-structured populations in periodically fluctuating environments

Living species, ranging from bacteria to animals, exist in environmental conditions that exhibit spatial and temporal heterogeneity which requires them to adapt. Risk-spreading through spontaneous phenotypic variations is a known concept in ecology, which is used to explain how species may survive when faced with the evolutionary risks associated with temporally varying environments. In order to support a deeper understanding of the adaptive role of spontaneous phenotypic variations in fluctuating environments, we consider a system of non-local partial differential equations modelling the evolutionary dynamics of two competing phenotype-structured populations in the presence of periodically oscillating nutrient levels. The two populations undergo heritable, spontaneous phenotypic variations at different rates. The phenotypic state of each individual is represented by a continuous variable, and the phenotypic landscape of the populations evolves in time due to variations in the nutrient level. Exploiting the analytical tractability of our model, we study the long-time behaviour of the solutions to obtain a detailed mathematical depiction of the evolutionary dynamics. The results suggest that when nutrient levels undergo small and slow oscillations, it is evolutionarily more convenient to rarely undergo spontaneous phenotypic variations. Conversely, under relatively large and fast periodic oscillations in the nutrient levels, which bring about alternating cycles of starvation and nutrient abundance, higher rates of spontaneous phenotypic variations confer a competitive advantage. We discuss the implications of our results in the context of cancer metabolism.

The role of spatial variations of abiotic factors in mediating intratumour phenotypic heterogeneity

We present here a space- and phenotype-structured model of selection dynamics between cancer cells within a solid tumour. In the framework of this model, we combine formal analyses with numerical simulations to investigate in silico the role played by the spatial distribution of abiotic components of the tumour microenvironment in mediating phenotypic selection of cancer cells. Numerical simulations are performed both on the 3D geometry of an in silico multicellular tumour spheroid and on the 3D geometry of an in vivo human hepatic tumour, which was imaged using computerised tomography. The results obtained show that inhomogeneities in the spatial distribution of oxygen, currently observed in solid tumours, can promote the creation of distinct local niches and lead to the selection of different phenotypic variants within the same tumour. This process fosters the emergence of stable phenotypic heterogeneity and supports the presence of hypoxic cells resistant to cytotoxic therapy prior to treatment. Our theoretical results demonstrate the importance of integrating spatial data with ecological principles when evaluating the therapeutic response of solid tumours to cytotoxic therapy.

Emergence of spatial patterns in a mathematical model for the co-culture dynamics of epithelial-like and mesenchymal-like cells

Accumulating evidence indicates that the interaction between epithelial and mesenchymal cells plays a pivotal role in cancer development and metastasis formation. Here we propose an integro-differential model for the co-culture dynamics of epithelial-like and mesenchymal-like cells. Our model takes into account the effects of chemotaxis, adhesive interactions between epithelial-like cells, proliferation and competition for nutrients. We present a sample of numerical results which display the emergence of spots, stripes and honeycomb patterns, depending on parameters and initial data. These simulations also suggest that epithelial-like and mesenchymal-like cells can segregate when there is little competition for nutrients. Furthermore, our computational results provide a possible explanation for how the concerted action between epithelial-cell adhesion and mesenchymal-cell spreading can precipitate the formation of ring-like structures, which resemble the fibrous capsules frequently observed in hepatic tumours.

Emergence of heterogeneity in acute leukemias

BACKGROUND

Leukemias are malignant proliferative disorders of the blood forming system. Sequencing studies demonstrate that the leukemic cell population consists of multiple clones. The genetic relationship between the different clones, referred to as the clonal hierarchy, shows high interindividual variability. So far, the source of this heterogeneity and its clinical relevance remain unknown. We propose a mathematical model to study the emergence and evolution of clonal heterogeneity in acute leukemias. The model allows linking properties of leukemic clones in terms of self-renewal and proliferation rates to the structure of the clonal hierarchy.

RESULTS

Computer simulations imply that the self-renewal potential of the first emerging leukemic clone has a major impact on the total number of leukemic clones and on the structure of their hierarchy. With increasing depth of the clonal hierarchy the self-renewal of leukemic clones increases, whereas the proliferation rates do not change significantly. The emergence of deep clonal hierarchies is a complex process that is facilitated by a cooperativity of different mutations.

CONCLUSION

Comparison of patient data and simulation results suggests that the self-renewal of leukemic clones increases with the emergence of clonal heterogeneity. The structure of the clonal hierarchy may serve as a marker for patient prognosis.

REVIEWERS

This article was reviewed by Marek Kimmel, Tommaso Lorenzi and Tomasz Lipniacki.

Tracking the evolution of cancer cell populations through the mathematical lens of phenotype-structured equations

Background

A thorough understanding of the ecological and evolutionary mechanisms that drive the phenotypic evolution of neoplastic cells is a timely and key challenge for the cancer research community. In this respect, mathematical modelling can complement experimental cancer research by offering alternative means of understanding the results of in vitro and in vivo experiments, and by allowing for a quick and easy exploration of a variety of biological scenarios through in silico studies.

Results

To elucidate the roles of phenotypic plasticity and selection pressures in tumour relapse, we present here a phenotype-structured model of evolutionary dynamics in a cancer cell population which is exposed to the action of a cytotoxic drug. The analytical tractability of our model allows us to investigate how the phenotype distribution, the level of phenotypic heterogeneity, and the size of the cell population are shaped by the strength of natural selection, the rate of random epimutations, the intensity of the competition for limited resources between cells, and the drug dose in use.

Conclusions

Our analytical results clarify the conditions for the successful adaptation of cancer cells faced with environmental changes. Furthermore, the results of our analyses demonstrate that the same cell population exposed to different concentrations of the same cytotoxic drug can take different evolutionary trajectories, which culminate in the selection of phenotypic variants characterised by different levels of drug tolerance. This suggests that the response of cancer cells to cytotoxic agents is more complex than a simple binary outcome, i.e., extinction of sensitive cells and selection of highly resistant cells. Also, our mathematical results formalise the idea that the use of cytotoxic agents at high doses can act as a double-edged sword by promoting the outgrowth of drug resistant cellular clones. Overall, our theoretical work offers a formal basis for the development of anti-cancer therapeutic protocols that go beyond the ‘maximum-tolerated-dose paradigm’, as they may be more effective than traditional protocols at keeping the size of cancer cell populations under control while avoiding the expansion of drug tolerant clones.

Reviewers

This article was reviewed by Angela Pisco, Sébastien Benzekry and Heiko Enderling.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0143-4) contains supplementary material, which is available to authorized users.

Dissecting the dynamics of epigenetic changes in phenotype-structured populations exposed to fluctuating environments

An enduring puzzle in evolutionary biology is to understand how individuals and populations adapt to fluctuating environments. Here we present an integro-differential model of adaptive dynamics in a phenotype-structured population whose fitness landscape evolves in time due to periodic environmental oscillations. The analytical tractability of our model allows for a systematic investigation of the relative contributions of heritable variations in gene expression, environmental changes and natural selection as drivers of phenotypic adaptation. We show that environmental fluctuations can induce the population to enter an unstable and fluctuation-driven epigenetic state. We demonstrate that this can trigger the emergence of oscillations in the size of the population, and we establish a full characterisation of such oscillations. Moreover, the results of our analyses provide a formal basis for the claim that higher rates of epimutations can bring about higher levels of intrapopulation heterogeneity, whilst intense selection pressures can deplete variation in the phenotypic pool of asexual populations. Finally, our work illustrates how the dynamics of the population size is led by a strong synergism between the rate of phenotypic variation and the frequency of environmental oscillations, and identifies possible ecological conditions that promote the maximisation of the population size in fluctuating environments.

Spatial heterogeneity in drug concentrations can facilitate the emergence of resistance to cancer therapy

Acquired resistance is one of the major barriers to successful cancer therapy. The development of resistance is commonly attributed to genetic heterogeneity. However, heterogeneity of drug penetration of the tumor microenvironment both on the microscopic level within solid tumors as well as on the macroscopic level across metastases may also contribute to acquired drug resistance. Here we use mathematical models to investigate the effect of drug heterogeneity on the probability of escape from treatment and the time to resistance. Specifically we address scenarios with sufficiently potent therapies that suppress growth of all preexisting genetic variants in the compartment with the highest possible drug concentration. To study the joint effect of drug heterogeneity, growth rate, and evolution of resistance, we analyze a multi-type stochastic branching process describing growth of cancer cells in multiple compartments with different drug concentrations and limited migration between compartments. We show that resistance is likely to arise first in the sanctuary compartment with poor drug penetrations and from there populate non-sanctuary compartments with high drug concentrations. Moreover, we show that only below a threshold rate of cell migration does spatial heterogeneity accelerate resistance evolution, otherwise deterring drug resistance with excessively high migration rates. Our results provide new insights into understanding why cancers tend to quickly become resistant, and that cell migration and the presence of sanctuary sites with little drug exposure are essential to this end.

Failure of cancer therapy is commonly attributed to the outgrowth of pre-existing resistant mutants already present prior to treatment, yet there is increasing evidence that the tumor microenvironment influences cell sensitivity to drugs and thus mediates the evolution of resistance during treatment. Here, we take into consideration important aspects of the tumor microenvironment, including spatial drug gradients and differential rates of cell proliferation. We show that the dependence of fitness on space together with cell migration facilitates the emergence of acquired resistance. Our analysis indicates that resistant cells that are selected for in compartments with high concentrations are likely to disseminate from sanctuary sites where they first acquire resistance preceding migration. The results suggest that it would be helpful to improve clinical outcomes by combining targeted therapy with anti-metastatic treatment aimed at constraining cell motility as well as by enhancing drug transportation and distribution throughout all metastatic compartments.

Emergence of drug tolerance in cancer cell populations: an evolutionary outcome of selection, nongenetic instability, and stress-induced adaptation

In recent experiments on isogenetic cancer cell lines, it was observed that exposure to high doses of anticancer drugs can induce the emergence of a subpopulation of weakly proliferative and drug-tolerant cells, which display markers associated with stem cell-like cancer cells. After a period of time, some of the surviving cells were observed to change their phenotype to resume normal proliferation and eventually repopulate the sample. Furthermore, the drug-tolerant cells could be drug resensitized following drug washout. Here, we propose a theoretical mechanism for the transient emergence of such drug tolerance. In this framework, we formulate an individual-based model and an integro-differential equation model of reversible phenotypic evolution in a cell population exposed to cytotoxic drugs. The outcomes of both models suggest that nongenetic instability, stress-induced adaptation, selection, and the interplay between these mechanisms can push an actively proliferating cell population to transition into a weakly proliferative and drug-tolerant state. Hence, the cell population experiences much less stress in the presence of the drugs and, in the long run, reacquires a proliferative phenotype, due to phenotypic fluctuations and selection pressure. These mechanisms can also reverse epigenetic drug tolerance following drug washout. Our study highlights how the transient appearance of the weakly proliferative and drug-tolerant cells is related to the use of high-dose therapy. Furthermore, we show how stem-like characteristics can act to stabilize the transient, weakly proliferative, and drug-tolerant subpopulation for a longer time window. Finally, using our models as in silico laboratories, we propose new testable hypotheses that could help uncover general principles underlying the emergence of cancer drug tolerance.

Hepatocellular carcinoma: a systems biology perspective

Hepatocellular carcinomas (HCCs) have different etiology and heterogenic genomic alterations lead to high complexity. The molecular features of HCC have largely been studied by gene expression and proteome profiling focusing on the correlations between the expression of specific markers and clinical data. Integration of the increasing amounts of data in databases has facilitated the link of genomic and proteomic profiles of HCC to disease state and clinical outcome. Despite the current knowledge, specific molecular markers remain to be identified and new strategies are required to establish novel-targeted therapies. In the last years, mathematical models reconstructing gene and protein networks based on experimental data of HCC have been developed providing powerful tools to predict candidate interactions and potential targets for therapy. Furthermore, the combination of dynamic and logical mathematical models with quantitative data allows detailed mechanistic insights into system properties. To address effects at the organ level, mathematical models reconstructing the three-dimensional organization of liver lobules were developed. In the future, integration of different modeling approaches capturing the effects at the cellular up to the organ level is required to address the complex properties of HCC and to enable the discovery of new targets for HCC prevention or treatment.