A multi-phenotypic cancer model with cell plasticity

Inferring absolute cell numbers from relative proportion in stochastic models with cell plasticity

Quantifying dynamic changes in cell populations is crucial for a comprehensive understanding of biological processes such as cell proliferation, injury repair, and disease progression. However, compared to directly measuring the absolute cell numbers of specific subpopulations, relative proportion data demonstrate greater reproducibility and yield more stable, reliable outcomes. Therefore, inferring absolute cell numbers from relative proportion data may present a novel approach for effectively predicting changes in cell population sizes. To address this, we establish two mathematical mappings between cell proportions and population sizes using moment equations derived from stochastic cell-plasticity models. Notably, our findings indicate that one of these mappings does not require prior knowledge of the initial population size, highlighting the value of incorporating variance information into cell proportion data. We evaluated the robustness of our methods from multiple perspectives and extended their application to various biological mechanisms within the context of cell plasticity models. These methods help mitigate the limitations associated with the direct measurement of absolute cell counts through experimental techniques. Moreover, they provide new insights into leveraging the stochastic dynamics of cell populations to quantify interactions between different biomasses within the system.

Negligible Long-Term Impact of Nonlinear Growth Dynamics on Heterogeneity in Models of Cancer Cell Populations

Linear compartmental models are often employed to capture the change in cell type composition of cancer cell populations. Yet, these populations usually grow in a nonlinear fashion. This begs the question of how linear compartmental models can successfully describe the dynamics of cell types. Here, we propose a general modeling framework with a nonlinear part capturing growth dynamics and a linear part capturing cell type transitions. We prove that dynamics in this general model are asymptotically equivalent to those governed only by its linear part under a wide range of assumptions for nonlinear growth.

Algorithms for the Uniqueness of the Longest Common Subsequence

Given several number sequences, determining the longest common subsequence is a classical problem in computer science. This problem has applications in bioinformatics, especially determining transposable genes. Nevertheless, related works only consider how to find one longest common subsequence. In this paper, we consider how to determine the uniqueness of the longest common subsequence. If there are multiple longest common subsequences, we also determine which number appears in all/some/none of the longest common subsequences. We focus on four scenarios: (1) linear sequences without duplicated numbers; (2) circular sequences without duplicated numbers; (3) linear sequences with duplicated numbers; (4) circular sequences with duplicated numbers. We develop corresponding algorithms and apply them to gene sequencing data.

Cell population growth kinetics in the presence of stochastic heterogeneity of cell phenotype

Recent studies at individual cell resolution have revealed phenotypic heterogeneity in nominally clonal tumor cell populations. The heterogeneity affects cell growth behaviors, which can result in departure from the idealized uniform exponential growth of the cell population. Here we measured the stochastic time courses of growth of an ensemble of populations of HL60 leukemia cells in cultures, starting with distinct initial cell numbers to capture a departure from the uniform exponential growth model for the initial growth ("take-off"). Despite being derived from the same cell clone, we observed significant variations in the early growth patterns of individual cultures with statistically significant differences in growth dynamics, which could be explained by the presence of inter-converting subpopulations with different growth rates, and which could last for many generations. Based on the hypothesis of existence of multiple subpopulations, we developed a branching process model that was consistent with the experimental observations.

Effects of a differentiating therapy on cancer-stem-cell-driven tumors

The growth of many solid tumors has been found to be driven by chemo- and radiotherapy-resistant cancer stem cells (CSCs). A suitable therapeutic avenue in these cases may involve the use of a differentiating agent (DA) to force the differentiation of the CSCs and of conventional therapies to eliminate the remaining differentiated cancer cells (DCCs). To describe the effects of a DA that reprograms CSCs into DCCs, we adapt a differential equation model developed to investigate tumorspheres, which are assumed to consist of jointly evolving CSC and DCC populations. We analyze the mathematical properties of the model, finding the equilibria and their stability. We also present numerical solutions and phase diagrams to describe the system evolution and the therapy effects, denoting the DA strength by a parameter adif. To obtain realistic predictions, we choose the other model parameters to be those determined previously from fits to various experimental datasets. These datasets characterize the progression of the tumor under various culture conditions. Typically, for small values of adif the tumor evolves towards a final state that contains a CSC fraction, but a strong therapy leads to the suppression of this phenotype. Nonetheless, different external conditions lead to very diverse behaviors. For microchamber-grown tumorspheres, there is a threshold in therapy strength below which both subpopulations survive, while high values of adif lead to the complete elimination of the CSC phenotype. For tumorspheres grown on hard and soft agar and in the presence of growth factors, the model predicts a threshold not only in the therapy strength, but also in its starting time, an early beginning being potentially crucial. In summary, our model shows how the effects of a DA depend critically not only on the dosage and timing of the drug application, but also on the tumor nature and its environment.

Order-of-mutation effects on cancer progression: models for myeloproliferative neoplasm

In some patients with myeloproliferative neoplasms (MPN), two genetic mutations are often found, JAK2 V617F and one in the TET2 gene. Whether or not one mutation is present will influence how the other subsequent mutation affects the regulation of gene expression. When both mutations are present, the order of their occurrence has been shown to influence disease progression and prognosis. We propose a nonlinear ordinary differential equation (ODE), Moran process, and Markov chain models to explain the non-additive and non-commutative mutation effects on recent clinical observations of gene expression patterns, proportions of cells with different mutations, and ages at diagnosis of MPN. These observations consistently shape our modeling framework. Our key proposal is that bistability in gene expression provides a natural explanation for many observed order-of-mutation effects. We also propose potential experimental measurements that can be used to confirm or refute predictions of our models.

Inference on autoregulation in gene expression with variance-to-mean ratio

Some genes can promote or repress their own expressions, which is called autoregulation. Although gene regulation is a central topic in biology, autoregulation is much less studied. In general, it is extremely difficult to determine the existence of autoregulation with direct biochemical approaches. Nevertheless, some papers have observed that certain types of autoregulations are linked to noise levels in gene expression. We generalize these results by two propositions on discrete-state continuous-time Markov chains. These two propositions form a simple but robust method to infer the existence of autoregulation from gene expression data. This method only needs to compare the mean and variance of the gene expression level. Compared to other methods for inferring autoregulation, our method only requires non-interventional one-time data, and does not need to estimate parameters. Besides, our method has few restrictions on the model. We apply this method to four groups of experimental data and find some genes that might have autoregulation. Some inferred autoregulations have been verified by experiments or other theoretical works.

Cell Population Growth Kinetics in the Presence of Stochastic Heterogeneity of Cell Phenotype

Recent studies at individual cell resolution have revealed phenotypic heterogeneity in nominally clonal tumor cell populations. The heterogeneity affects cell growth behaviors, which can result in departure from the idealized exponential growth. Here we measured the stochastic time courses of growth of an ensemble of populations of HL60 leukemia cells in cultures, starting with distinct initial cell numbers to capture the departure from the exponential growth model in the initial growth phase. Despite being derived from the same cell clone, we observed significant variations in the early growth patterns of individual cultures with statistically significant differences in growth kinetics and the presence of subpopulations with different growth rates that endured for many generations. Based on the hypothesis of existence of multiple inter-converting subpopulations, we developed a branching process model that captures the experimental observations.

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.

Growth dynamics of breast cancer stem cells: effects of self-feedback and EMT mechanisms

Breast cancer stem cells (BCSCs) with the ability to self-renew and differentiate have been identified in primary breast cancer tissues and cell lines. The BCSCs are often resistant to traditional radiation and/or chemotherapies. Previous studies have also shown that successful therapy must eradicate cancer stem cells. The purpose of this paper is to develop a mathematical model with self-feedback mechanism to illustrate the issues regarding the difficulties of absolutely eliminating a breast cancer. In addition, we introduce the mechanism of the epithelial-mesenchymal transition (EMT) to investigate the influence of EMT on the effects of breast cancer growth and treatment. Results indicate that the EMT mechanism facilitates the growth of breast cancer and makes breast cancer more difficult to be cured. Therefore, targeting the signals involved in EMT can halt tumor progression in breast cancer. Finally, we apply the experimental data to carry out numerical simulations and validate our theoretical conclusions.

Inference on the structure of gene regulatory networks

In this paper, we conduct theoretical analyses on inferring the structure of gene regulatory networks. Depending on the experimental method and data type, the inference problem is classified into 20 different scenarios. For each scenario, we discuss the problem that with enough data, under what assumptions, what can be inferred about the structure. For scenarios that have been covered in the literature, we provide a brief review. For scenarios that have not been covered in literature, if the structure can be inferred, we propose new mathematical inference methods and evaluate them on simulated data. Otherwise, we prove that the structure cannot be inferred.

Alternative models of cancer stem cells: The stemness phenotype model, 10 years later

The classical cancer stem cell (CSCs) theory proposed the existence of a rare but constant subpopulation of CSCs. In this model cancer cells are organized hierarchically and are responsible for tumor resistance and tumor relapse. Thus, eliminating CSCs will eventually lead to cure of cancer. This simplistic model has been challenged by experimental data. In 2010 we proposed a novel and controversial alternative model of CSC biology (the Stemness Phenotype Model, SPM). The SPM proposed a non-hierarchical model of cancer biology in which there is no specific subpopulation of CSCs in tumors. Instead, cancer cells are highly plastic in term of stemness and CSCs and non-CSCs can interconvert into each other depending on the microenvironment. This model predicts the existence of cancer cells ranging from a pure CSC phenotype to pure non-CSC phenotype and that survival of a single cell can originate a new tumor. During the past 10 years, a plethora of experimental evidence in a variety of cancer types has shown that cancer cells are indeed extremely plastic and able to interconvert into cells with different stemness phenotype. In this review we will (1) briefly describe the cumulative evidence from our laboratory and others supporting the SPM; (2) the implications of the SPM in translational oncology; and (3) discuss potential strategies to develop more effective therapeutic regimens for cancer treatment.

The invasion of de-differentiating cancer cells into hierarchical tissues

Many fast renewing tissues are characterized by a hierarchical cellular architecture, with tissue specific stem cells at the root of the cellular hierarchy, differentiating into a whole range of specialized cells. There is increasing evidence that tumors are structured in a very similar way, mirroring the hierarchical structure of the host tissue. In some tissues, differentiated cells can also revert to the stem cell phenotype, which increases the risk that mutant cells lead to long lasting clones in the tissue. However, it is unclear under which circumstances de-differentiating cells will invade a tissue. To address this, we developed mathematical models to investigate how de-differentiation is selected as an adaptive mechanism in the context of cellular hierarchies. We derive thresholds for which de-differentiation is expected to emerge, and it is shown that the selection of de-differentiation is a result of the combination of the properties of cellular hierarchy and de-differentiation patterns. Our results suggest that de-differentiation is most likely to be favored provided stem cells having the largest effective self-renewal rate. Moreover, jumpwise de-differentiation provides a wider range of favorable conditions than stepwise de-differentiation. Finally, the effect of de-differentiation on the redistribution of self-renewal and differentiation probabilities also greatly influences the selection for de-differentiation.

How can a tissue such as the blood system or the skin, which constantly produces a huge number of cells, avoids that errors accumulate in the cells over time? Such tissues are typically organized in cellular hierarchies, which induce a directional relation between different stages of cellular differentiation, minimizing the risk of retention of mutations. However, recent evidence also shows that some differentiated cells can de-differentiate into the stem cell phenotype. Why does de-differentiation arise in some tumors, but not in others? We developed a mathematical model to study the growth competition between de-differentiating mutant cell populations and non de-differentiating resident cell population. Our results suggest that the invasion of de-differentiation is jointly influenced by the cellular hierarchy (e.g. number of cell compartments, inherent cell division pattern) and the de-differentiation pattern, i.e. how exactly cells acquire their stem-cell like properties.

Interplay of Darwinian Selection, Lamarckian Induction and Microvesicle Transfer on Drug Resistance in Cancer

Development of drug resistance in cancer has major implications for patients’ outcome. It is related to processes involved in the decrease of drug efficacy, which are strongly influenced by intratumor heterogeneity and changes in the microenvironment. Heterogeneity arises, to a large extent, from genetic mutations analogously to Darwinian evolution, when selection of tumor cells results from the adaptation to the microenvironment, but could also emerge as a consequence of epigenetic mutations driven by stochastic events. An important exogenous source of alterations is the action of chemotherapeutic agents, which not only affects the signalling pathways but also the interactions among cells. In this work we provide experimental evidence from in vitro assays and put forward a mathematical kinetic transport model to describe the dynamics displayed by a system of non-small-cell lung carcinoma cells (NCI-H460) which, depending on the effect of a chemotherapeutic agent (doxorubicin), exhibits a complex interplay between Darwinian selection, Lamarckian induction and the nonlocal transfer of extracellular microvesicles. The role played by all of these processes to multidrug resistance in cancer is elucidated and quantified.

Modeling differentiation-state transitions linked to therapeutic escape in triple-negative breast cancer

Drug resistance in breast cancer cell populations has been shown to arise through phenotypic transition of cancer cells to a drug-tolerant state, for example through epithelial-to-mesenchymal transition or transition to a cancer stem cell state. However, many breast tumors are a heterogeneous mixture of cell types with numerous epigenetic states in addition to stem-like and mesenchymal phenotypes, and the dynamic behavior of this heterogeneous mixture in response to drug treatment is not well-understood. Recently, we showed that plasticity between differentiation states, as identified with intracellular markers such as cytokeratins, is linked to resistance to specific targeted therapeutics. Understanding the dynamics of differentiation-state transitions in this context could facilitate the development of more effective treatments for cancers that exhibit phenotypic heterogeneity and plasticity. In this work, we develop computational models of a drug-treated, phenotypically heterogeneous triple-negative breast cancer (TNBC) cell line to elucidate the feasibility of differentiation-state transition as a mechanism for therapeutic escape in this tumor subtype. Specifically, we use modeling to predict the changes in differentiation-state transitions that underlie specific therapy-induced changes in differentiation-state marker expression that we recently observed in the HCC1143 cell line. We report several statistically significant therapy-induced changes in transition rates between basal, luminal, mesenchymal, and non-basal/non-luminal/non-mesenchymal differentiation states in HCC1143 cell populations. Moreover, we validate model predictions on cell division and cell death empirically, and we test our models on an independent data set. Overall, we demonstrate that changes in differentiation-state transition rates induced by targeted therapy can provoke distinct differentiation-state aggregations of drug-resistant cells, which may be fundamental to the design of improved therapeutic regimens for cancers with phenotypic heterogeneity.

Phenotypic equilibrium as probabilistic convergence in multi-phenotype cell population dynamics

We consider the cell population dynamics with n different phenotypes. Both the Markovian branching process model (stochastic model) and the ordinary differential equation (ODE) system model (deterministic model) are presented, and exploited to investigate the dynamics of the phenotypic proportions. We will prove that in both models, these proportions will tend to constants regardless of initial population states ("phenotypic equilibrium") under weak conditions, which explains the experimental phenomenon in Gupta et al.'s paper. We also prove that Gupta et al.'s explanation is the ODE model under a special assumption. As an application, we will give sufficient and necessary conditions under which the proportion of one phenotype tends to 0 (die out) or 1 (dominate). We also extend our results to non-Markovian cases.

Mathematical modelling of plasticity and phenotype switching in cancer cell populations

The cancer stem cell (CSC) hypothesis suggests that cancer stem cells proliferate via a hierarchical model of unidirectional differentiation. However, growing experimental evidence has advanced this hypothesis by introducing a bidirectional hierarchy, in which non-CSCs may dedifferentiate into CSCs. Various models have been developed enabling the incorporation of this plasticity within cancer cell populations, focusing on behaviour in the limit of a large number of cells. However, stochastic effects predominate in the limit of small numbers of cells, which correlates with biologically relevant assays such as the mammosphere formation assay (MFA). Here, we consider two mathematical models incorporating cellular plasticity, namely a two-compartment model and a hierarchical model, and by parameterizing these models with experimental data, we show this behavioural difference in the limits of large and small numbers of cells. Additionally, we analyse the effects of varying cellular plasticity on the survival of the cancer cell population, and show that interestingly, increased plasticity, in certain cases, may be advantageous in reducing the survival probability. Thus, this analysis highlights the necessity of experimentally studying both small and large populations of cancer cells concurrently to obtain valid model predictions, potentially aiding the design of novel therapeutics.

Variation is function: Are single cell differences functionally important?: Testing the hypothesis that single cell variation is required for aggregate function

There is a growing appreciation of the extent of transcriptome variation across individual cells of the same cell type. While expression variation may be a byproduct of, for example, dynamic or homeostatic processes, here we consider whether single-cell molecular variation per se might be crucial for population-level function. Under this hypothesis, molecular variation indicates a diversity of hidden functional capacities within an ensemble of identical cells, and this functional diversity facilitates collective behavior that would be inaccessible to a homogenous population. In reviewing this topic, we explore possible functions that might be carried by a heterogeneous ensemble of cells; however, this question has proven difficult to test, both because methods to manipulate molecular variation are limited and because it is complicated to define, and measure, population-level function. We consider several possible methods to further pursue the hypothesis that variation is function through the use of comparative analysis and novel experimental techniques.

The overshoot and phenotypic equilibrium in characterizing cancer dynamics of reversible phenotypic plasticity

The paradigm of phenotypic plasticity indicates reversible relations of different cancer cell phenotypes, which extends the cellular hierarchy proposed by the classical cancer stem cell (CSC) theory. Since it is still questionable if the phenotypic plasticity is a crucial improvement to the hierarchical model or just a minor extension to it, it is worthwhile to explore the dynamic behavior characterizing the reversible phenotypic plasticity. In this study we compare the hierarchical model and the reversible model in predicting the cell-state dynamics observed in biological experiments. Our results show that the hierarchical model shows significant disadvantages over the reversible model in describing both long-term stability (phenotypic equilibrium) and short-term transient dynamics (overshoot) in cancer cell populations. In a very specific case in which the total growth of population due to each cell type is identical, the hierarchical model predicts neither phenotypic equilibrium nor overshoot, whereas the reversible model succeeds in predicting both of them. Even though the performance of the hierarchical model can be improved by relaxing the specific assumption, its prediction to the phenotypic equilibrium strongly depends on a precondition that may be unrealistic in biological experiments. Moreover, it still does not show as rich dynamics as the reversible model in capturing the overshoots of both CSCs and non-CSCs. By comparison, it is more likely for the reversible model to correctly predict the stability of the phenotypic mixture and various types of overshoot behavior.

The phenotypic equilibrium of cancer cells: From average-level stability to path-wise convergence

The phenotypic equilibrium, i.e. heterogeneous population of cancer cells tending to a fixed equilibrium of phenotypic proportions, has received much attention in cancer biology very recently. In the previous literature, some theoretical models were used to predict the experimental phenomena of the phenotypic equilibrium, which were often explained by different concepts of stabilities of the models. Here we present a stochastic multi-phenotype branching model by integrating conventional cellular hierarchy with phenotypic plasticity mechanisms of cancer cells. Based on our model, it is shown that: (i) our model can serve as a framework to unify the previous models for the phenotypic equilibrium, and then harmonizes the different kinds of average-level stabilities proposed in these models; and (ii) path-wise convergence of our model provides a deeper understanding to the phenotypic equilibrium from stochastic point of view. That is, the emergence of the phenotypic equilibrium is rooted in the stochastic nature of (almost) every sample path, the average-level stability just follows from it by averaging stochastic samples.

MCL-1, BCL-XL and MITF Are Diversely Employed in Adaptive Response of Melanoma Cells to Changes in Microenvironment

Melanoma cells can switch their phenotypes in response to microenvironmental insults. Heterogeneous melanoma populations characterized by long-term growth and a high self-renewal capacity can be obtained in vitro in EGF(+)bFGF(+) medium whilst invasive potential of melanoma cells is increased in serum-containing cultures. In the present study, we have shown that originally these patient-derived melanoma populations exhibit variable expression of pro-survival genes from the BCL-2 family and inhibitors of apoptosis (IAPs), and differ in the baseline MCL-1 transcript stability as well. While being transferred to serum-containing medium, melanoma cells are well protected from death. Immediate adaptive response of melanoma cells selectively involves a temporary MCL-1 increase, both at mRNA and protein levels, and BCL-X_L can complement MCL-1, especially in MITF_low populations. Thus, the extent of MCL-1 and BCL-XL contributions seems to be cell context-dependent. An increase in MCL-1 level results from a transiently enhanced stability of its transcript, but not from altered protein turnover. Inhibition of MCL-1 preceding transfer to serum-containing medium caused the induction of cell death in a subset of melanoma cells, which confirms the involvement of MCL-1 in melanoma cell survival during the rapid alteration of growth conditions. Additionally, immediate response to serum involves the transient increase in MITF expression and inhibition of ERK-1/2 activity. Uncovering the mechanisms of adaptive response to rapid changes in microenvironment may extend our knowledge on melanoma biology, especially at the stage of dissemination.

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.