Mathematical modelling of plasticity and phenotype switching in cancer cell populations

Modeling stress-induced responses: plasticity in continuous state space and gradual clonal evolution

Mathematical models of cancer and bacterial evolution have generally stemmed from a gene-centric framework, assuming clonal evolution via acquisition of resistance-conferring mutations and selection of their corresponding subpopulations. More recently, the role of phenotypic plasticity has been recognized and models accounting for phenotypic switching between discrete cell states (e.g., epithelial and mesenchymal) have been developed. However, seldom do models incorporate both plasticity and mutationally driven resistance, particularly when the state space is continuous and resistance evolves in a continuous fashion. In this paper, we develop a framework to model plastic and mutational mechanisms of acquiring resistance in a continuous gradual fashion. We use this framework to examine ways in which cancer and bacterial populations can respond to stress and consider implications for therapeutic strategies. Although we primarily discuss our framework in the context of cancer and bacteria, it applies broadly to any system capable of evolving via plasticity and genetic evolution.

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

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

On tumoural growth and treatment under cellular dedifferentiation

Differentiated cancer cells may regain stem cell characteristics; however, the effects of such a cellular dedifferentiation on tumoural growth and treatment are currently understudied. Thus, we here extend a mathematical model of cancer stem cell (CSC) driven tumour growth to also include dedifferentiation. We show that dedifferentiation increases the likelihood of tumorigenesis and the speed of tumoural growth, both modulated by the proliferative potential of the non-stem cancer cells (NSCCs). We demonstrate that dedifferentiation also may lead to treatment evasion, especially when a treatment solely targets CSCs. Conversely, targeting both CSCs and NSCCs in parallel is shown to be more robust to dedifferentiation. Despite dedifferentiation, perturbing CSC-related parameters continues to exert the largest relative effect on tumoural growth; however, we show the existence of synergies between specific CSC- and NSCC-directed treatments which cause superadditive reductions of tumoural growth. Overall, our study demonstrates various effects of dedifferentiation on growth and treatment of tumoural lesions, and we anticipate our results to be helpful in guiding future molecular and clinical research on limiting tumoural growth in vivo.

Optimal chemotherapy counteracts cancer adaptive resistance in a cell-based, spatially-extended, evolutionary model

Most aggressive cancers are incurable due to their fast evolution of drug resistance. We model cancer growth and adaptive response in a simplified cell-based (CB) setting, assuming a genetic resistance to two chemotherapeutic drugs. We show that optimal administration protocols can steer cells resistance and turned it into a weakness for the disease. Our work extends the population-based (PB) model proposed by Orlando et al. (Physical Biology, 2012), in which a homogeneous population of cancer cells evolves according to a fitness landscape. The landscape models three types of trade-offs, differing on whether the cells are more, less, or equal effective when generalizing resistance to two drugs as opposed to specializing to a single one. The CB framework allows us to include genetic heterogeneity, spatial competition, and drugs diffusion, as well as realistic administration protocols. By calibrating our model on Orlando et al.'s assumptions, we show that dynamical protocols that alternate the two drugs minimize the cancer size at the end of (or at mid-points during) treatment. These results significantly differ from those obtained with the homogeneous model---suggesting static protocols under the pro-generalizing and neutral allocation trade-offs---highlighting the important role of spatial and genetic heterogeneities. Our work is the first attempt to search for optimal treatments in a CB setting, a step forward toward realistic clinical applications.

MicroRNA-222 Regulates Melanoma Plasticity

Melanoma is one of the most aggressive and highly resistant tumors. Cell plasticity in melanoma is one of the main culprits behind its metastatic capabilities. The detailed molecular mechanisms controlling melanoma plasticity are still not completely understood. Here we combine mathematical models of phenotypic switching with experiments on IgR39 human melanoma cells to identify possible key targets to impair phenotypic switching. Our mathematical model shows that a cancer stem cell subpopulation within the tumor prevents phenotypic switching of the other cancer cells. Experiments reveal that hsa-mir-222 is a key factor enabling this process. Our results shed new light on melanoma plasticity, providing a potential target and guidance for therapeutic studies.

Effect of cellular de-differentiation on the dynamics and evolution of tissue and tumor cells in mathematical models with feedback regulation

Tissues are maintained by adult stem cells that self-renew and also differentiate into functioning tissue cells. Homeostasis is achieved by a set of complex mechanisms that involve regulatory feedback loops. Similarly, tumors are believed to be maintained by a minority population of cancer stem cells, while the bulk of the tumor is made up of more differentiated cells, and there is indication that some of the feedback loops that operate in tissues continue to be functional in tumors. Mathematical models of such tissue hierarchies, including feedback loops, have been analyzed in a variety of different contexts. Apart from stem cells giving rise to differentiated cells, it has also been observed that more differentiated cells can de-differentiate into stem cells, both in healthy tissue and tumors, aspects of which have also been investigated mathematically. This paper analyses the effect of de-differentiation on the basic and evolutionary dynamics of cells in the context of tissue hierarchy models that include negative feedback regulation of the cell populations. The models predict that in the presence of de-differentiation, the fixation probability of a neutral mutant is lower than in its absence. Therefore, if de-differentiation occurs, a mutant with identical parameters compared to the wild-type cell population behaves like a disadvantageous mutant. Similarly, the process of de-differentiation is found to lower the fixation probability of an advantageous mutant. These results indicate that the presence of de-differentiation can lower the rates of tumor initiation and progression in the context of the models considered here.