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

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.

A modelling framework for cancer ecology and evolution

Cancer incidence increases rapidly with age, typically as a polynomial. The somatic mutation theory explains this increase through the waiting time for enough mutations to build up to generate cells with the full set of traits needed to grow without control. However, lines of evidence ranging from tumour reversion and dormancy to the prevalence of presumed cancer mutations in non-cancerous tissues argue that this is not the whole story, and that cancer is also an ecological process, and that mutations only lead to cancer when the systems of control within and across cells have broken down. Aging thus has two effects: the build-up of mutations and the breakdown of control. This paper presents a mathematical modelling framework to unify these theories with novel approaches to model the mutation and diversification of cell lineages and of the breakdown of the layers of control both within and between cells. These models correctly predict the polynomial increase of cancer with age, show how germline defects in control accelerate cancer initiation, and compute how the positive feedback between cell replication, ecology and layers of control leads to a doubly exponential growth of cell populations.

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.

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.

Stochastic stem cell models with mutation: A comparison of asymmetric and symmetric divisions

In order to fulfill cell proliferation and differentiation through cellular hierarchy, stem cells can undergo either asymmetric or symmetric divisions. Recent studies pay special attention to the effect of different modes of stem cell division on the lifetime risk of cancer, and report that symmetric division is more beneficial to delay the onset of cancer. The fate uncertainty of symmetric division is considered to be the reason for the cancer-delaying effect. In this paper we compare asymmetric and symmetric divisions of stem cells via studying stochastic stem cell models with mutation. Specially, by using rigorous mathematical analysis we find that both the asymmetric and symmetric models show the same statistical average, but the symmetric model shows higher fluctuation than the asymmetric model. We further show that the difference between the two models would be more remarkable for lower mutation rates. Our work quantifies the uncertainty of cell division and highlights the significance of stochasticity for distinguishing between different modes of stem cell division.

Dynamics of Phenotypic Heterogeneity Associated with EMT and Stemness during Cancer Progression

Genetic and phenotypic heterogeneity contribute to the generation of diverse tumor cell populations, thus enhancing cancer aggressiveness and therapy resistance. Compared to genetic heterogeneity, a consequence of mutational events, phenotypic heterogeneity arises from dynamic, reversible cell state transitions in response to varying intracellular/extracellular signals. Such phenotypic plasticity enables rapid adaptive responses to various stressful conditions and can have a strong impact on cancer progression. Herein, we have reviewed relevant literature on mechanisms associated with dynamic phenotypic changes and cellular plasticity, such as epithelial-mesenchymal transition (EMT) and cancer stemness, which have been reported to facilitate cancer metastasis. We also discuss how non-cell-autonomous mechanisms such as cell-cell communication can lead to an emergent population-level response in tumors. The molecular mechanisms underlying the complexity of tumor systems are crucial for comprehending cancer progression, and may provide new avenues for designing therapeutic strategies.

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.

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.

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.