Cancer stem cells are controlled by developmental networks that are often topologically indistinguishable from normal, healthy stem cells. The question is why cancer stem cells can be both phenotypically distinct and have morphological effects so different from normal stem cells. The difference between cancer stem cells and normal stem cells lies not in differences their network architecture, but rather in the spatial-temporal locality of their activation in the genome and the resulting expression in the body. The metastatic potential cancer stem cells is not based primarily on their network divergence from normal stem cells, but on non-network based genetic changes that enable the evolution of gene-based phenotypic properties of the cell that permit its escape and travel to other parts of the body. Stem cell network theory allows the precise prediction of stem cell behavioral dynamics and a mathematical description of stem cell proliferation for both normal and cancer stem cells. It indicates that the best therapeutic approach is to tackle the highest order stem cells first, otherwise spontaneous remission of so called cured cancers will always be a danger. Stem cell networks poin
Stem cells have the potential to produce lineages of non-stem cell populations (differentiated cells) via a ubiquitous hierarchal division scheme. Differentiation of a stem cell into (partially) differentiated cells can happen either symmetrically or asymmetrically. The selection dynamics of a mutant cancer stem cell should be investigated in the light of a stem cell proliferation hierarchy and presence of a non-stem cell population. By constructing a three-compartment Moran-type model composed of normal stem cells, mutant (cancer) stem cells and differentiated cells, we derive the replicator dynamics of stem cell frequencies where asymmetric differentiation and differentiated cell death rates are included in the model. We determine how these new factors change the conditions for a successful mutant invasion and discuss the variation on the steady state fraction of the population as different model parameters are changed. By including the phenotypic plasticity/dedifferentiation, in which a progenitor/differentiated cell can transform back into a cancer stem cell, we show that the effective fitness of mutant stem cells is not only determined by their proliferation and death rates bu
Induced pluripotent stem cells (iPSCs) provide a great model to study the process of reprogramming and differentiation of stem cells. Single-cell RNA sequencing (scRNA-seq) enables us to investigate the reprogramming process at single-cell level. Here, we introduce single-cell entropy (scEntropy) as a macroscopic variable to quantify the cellular transcriptome from scRNA-seq data during reprogramming and differentiation of iPSCs. scEntropy measures the relative order parameter of genomic transcriptions at single cell level during the cell fate change process, which shows increasing during differentiation, and decreasing upon reprogramming. Moreover, based on the scEntropy dynamics, we construct a phenomenological stochastic differential equation model and the corresponding Fokker-Plank equation for cell state transitions during iPSC differentiation, which provide insights to infer cell fates changes and stem cell differentiation. This study is the first to introduce the novel concept of scEntropy to the biological process of iPSC, and suggests that the scEntropy can provide a suitable quantify to describe cell fate transition in differentiation and reprogramming of stem cells.
Stem cells can precisely and robustly undergo cellular differentiation and lineage commitment, referred to as stemness. However, how the gene network underlying stemness regulation reliably specifies cell fates is not well understood. To address this question, we applied a recently developed computational method, Random Circuit Perturbation (RACIPE), to a nine-component gene regulatory network (GRN) governing stemness, from which we identified fifteen robust gene states. Among them, four out of the five most probable gene states exhibit gene expression patterns observed in single mouse embryonic cells at 32-cell and 64-cell stages. These gene states can be robustly predicted by the stemness GRN but not by randomized versions of the stemness GRN. Strikingly, we found a hierarchical structure of the GRN with the Oct4/Cdx2 motif functioning as the first decision-making module followed by Gata6/Nanog. We propose that stem cell populations, instead of being viewed as all having a specific cellular state, can be regarded as a heterogeneous mixture including cells in various states. Upon perturbations by external signals, stem cells lose the capacity to access certain cellular states, the
Cycling tissues such as the intestinal epithelium, germ line, and hair follicles, require a constant flux of differentiated cells. These tissues are maintained by a population of stem cells, which generate differentiated progenies and self-renew. Asymmetric division of each stem cell into one stem cell and one differentiated cell can accomplish both tasks. However, in mammalian cycling tissues, some stem cells divide symmetrically into two differentiated cells and are replaced by a neighbor that divides symmetrically into two stem cells. Besides this heterogeneity in fate (population asymmetry), stem cells also exhibit heterogenous proliferation-rates; in the long run, however, all stem cells proliferate at the same average rate (equipotency). We construct and simulate a mathematical model based on these experimental observations. We show that the complex steady-state dynamics of population-asymmetric stem cells reduces the rate of replicative aging of the tissue --potentially lowering the incidence of somatic mutations and genetics diseases such as cancer. Essentially, slow-dividing stem cells proliferate and purge the population of the fast-dividing --older-- cells which had unde
Cancer stem cells have been shown to be critical to the development of a variety of solid cancers. The precise interplay mechanisms between cancer stem cells and the rest of a tissue are still not elucidated. To shed light on the interactions between stem and non-stem cancer cell populations we develop a two-population mathematical model, which is suitable to describe tumorsphere growth. Both interspecific and intraspecific interactions, mediated by the microenvironment, are included. We show that there is a tipping point, characterized by a transcritical bifurcation, where a purely non-stem cell attractor is replaced by a new attractor that contains both stem and differentiated cancer cells. The model is then applied to describe the outcome of a recent experiment. This description reveals that, while the intraspecific interactions are inhibitory, the interspecific interactions stimulate growth. This can be understood in terms of stem cells needing differentiated cells to reinforce their niches, and phenotypic plasticity favoring the de-differentiation of differentiated cells into cancer stem cells. We posit that this is a consequence of the deregulation of the quorum sensing that
Statistical and mathematical modeling are crucial to describe, interpret, compare and predict the behavior of complex biological systems including the organization of hematopoietic stem and progenitor cells in the bone marrow environment. The current prominence of high-resolution and live-cell imaging data provides an unprecedented opportunity to study the spatiotemporal dynamics of these cells within their stem cell niche and learn more about aberrant, but also unperturbed, normal hematopoiesis. However, this requires careful quantitative statistical analysis of the spatial and temporal behavior of cells and the interaction with their microenvironment. Moreover, such quantification is a prerequisite for the construction of hypothesis-driven mathematical models that can provide mechanistic explanations by generating spatiotemporal dynamics that can be directly compared to experimental observations. Here, we provide a brief overview of statistical methods in analyzing spatial distribution of cells, cell motility, cell shapes and cellular genealogies. We also describe cell- based modeling formalisms that allow researchers to simulate emergent behavior in a multicellular system based
We present a general computational theory of stem cell networks and their developmental dynamics. Stem cell networks are special cases of developmental control networks. Our theory generates a natural classification of all possible stem cell networks based on their network architecture. Each stem cell network has a unique topology and semantics and developmental dynamics that result in distinct phenotypes. We show that the ideal growth dynamics of multicellular systems generated by stem cell networks have mathematical properties related to the coefficients of Pascal's Triangle. The relationship to cancer stem cells and their control networks is indicated. The theory lays the foundation for a new research paradigm for understanding and investigating stem cells. The theory of stem cell networks implies that new methods for generating and controlling stem cells will become possible.
The crawling motility of many eukaryotic cells is driven by filamentous actin (F-actin), and regulated by a network of signaling proteins and lipids (including small GTPases). The tangle of positive and negative feedback loops gives rise to various experimentally observed dynamic patterns (``actin waves''). Here we consider a recent prototypical model for actin waves in which F-actin exerts negative feedback onto a GTPase. Guided by recent numerical PDE bifurcation analysis in Hughes (2025) and Hughes et al (2026), we explore cell shapes and motility associated with polar, oscillatory, and traveling waves solutions of a mass-conserved partial differential equation (PDE) model. We use Morpheus (cellular Potts) simulations to investigate the implications of such regimes of behavior on the shapes and motion of cells, and on transitions between modes of behavior. The model demonstrates various cell states, including resting (spatially uniform GTPase), polar cells (static ``zones'' of GTPase), and traveling waves along the cell edge. In some parameter regimes, such states can coexist, so that cells can transition from one behavior to another in response to noisy stimuli.
Stem cell heterogeneity is essential for the homeostasis in tissue development. This paper established a general formulation for understanding the dynamics of stem cell regeneration with cell heterogeneity and random transitions of epigenetic states. The model generalizes the classical G0 cell cycle model, and incorporates the epigenetic states of stem cells that are represented by a continuous multidimensional variable and the kinetic rates of cell behaviors, including proliferation, differentiation, and apoptosis, that are dependent on their epigenetic states. Moreover, the random transition of epigenetic states is represented by an inheritance probability that can be described as a conditional beta distribution. This model can be extended to investigate gene mutation-induced tumor development. The proposed formula is a generalized formula that helps us to understand various dynamic processes of stem cell regeneration, including tissue development, degeneration, and abnormal growth.
Oncolytic viruses, which may be naturally occurring or genetically engineered, are a type of virus that infects and destroy cancer cells preferentially. Owing to their selectivity, they outperform conventional chemotherapy and radiotherapy, which both have a tendency to impact non-target cells and cause unwanted adverse side effects. Oncolytic virotherapy is a type of cancer treatment in which oncolytic viruses are deliberately introduced into patients affected with cancers in order for them to infect and destroy cancer cells locally or systemically, in a manner analogous to chemotherapy but with a greater degree of selectivity. Multiple studies indicate that oncolytic virotherapy is effective in vitro but in vivo findings remain ambiguous due to the approach's primary limitation: inefficient therapeutic agent delivery to its target, which is heavily influenced by the immune system. Here, we propose overcoming this limitation by exploiting a recent discovery in cancer research: a carrier cell. By exploiting their tumor-promoting activities, mesenchymal stem cells may be employed for cancer therapy by serving as a carrier for the oncolytic viruses toward their target. This approach
Stem cell regeneration is a vital biological process in self-renewing tissues, governing development and tissue homeostasis. Gene regulatory network dynamics are pivotal in controlling stem cell regeneration and cell type transitions. However, integrating the quantitative dynamics of gene regulatory networks at the single-cell level with stem cell regeneration at the population level poses significant challenges. This study presents a computational framework connecting gene regulatory network dynamics with stem cell regeneration through a data-driven formulation of the inheritance function. The inheritance function captures epigenetic state transitions during cell division in heterogeneous stem cell populations. Our scheme allows the derivation of the inheritance function based on a hybrid model of cross-cell-cycle gene regulation network dynamics. The proposed scheme enables us to derive the inheritance function based on the hybrid model of cross-cell-cycle gene regulation network dynamics. By explicitly incorporating gene regulatory network structure, it replicates cross-cell-cycling gene regulation dynamics through individual-cell-based modeling. The numerical scheme holds the p
Heterogeneity is ubiquitous in stem cells (SC), cancer cells (CS), and cancer stem cells (CSC). SC and CSC heterogeneity is manifested as diverse sub-populations with self-renewing and unique regeneration capacity. Moreover, the CSC progeny possesses multiple plasticity and cancerous characteristics. Many studies have demonstrated that cancer heterogeneity is one of the greatest obstacle for therapy. This leads to the incomplete anti-cancer therapies and transitory efficacy. Furthermore, numerous micro-metastasis leads to the wide spread of the tumor cells across the body which is the beginning of metastasis. The epigenetic processes (DNA methylation or histone remodification etc.) can provide a source for certain heterogeneity. In this study, we develop a mathematical model to quantify the heterogeneity of SC, CSC and cancer taking both genetic and epigenetic effects into consideration. We uncovered the roles and physical mechanisms of heterogeneity from the three aspects (SC, CSC and cancer). In the adiabatic regime (relatively fast regulatory binding and effective coupling among genes), seven native states (SC, CSC, Cancer, Premalignant, Normal, Lesion and Hyperplasia) emerge. I
In many adult tissues, stem cells and differentiated cells are not homogeneously distributed : stem cells are arranged in periodic "niches", and differentiated cells are constantly produced and migrate out of these niches. In this article, we provide a general theoretical framework to study mixtures of dividing and actively migrating particles, which we apply to biological tissues. We show in particular that the interplay between the stresses arising from active cell migration and stem cell division give rise to robust stem cell patterns. The instability of the tissue leads to spatial patterns which are either steady or oscillating in time. The wavelength of the instability has an order of magnitude consistent with the biological observations. We also discuss the implications of these results for future in vitro and in vivo experiments.
We generated a computational approach to analyze the biomechanics of epithelial cell aggregates, either island or stripes or entire monolayers, that combines both vertex and contact-inhibition-of-locomotion models to include both cell-cell and cell-substrate adhesion. Examination of the distribution of cell protrusions (adhesion to the substrate) in the model predicted high order profiles of cell organization that agree with those previously seen experimentally. Cells acquired an asymmetric distribution of basal protrusions, traction forces and apical aspect ratios that decreased when moving from the edge to the island center. Our in silico analysis also showed that tension on cell-cell junctions and apical stress is not homogeneous across the island. Instead, these parameters are higher at the island center and scales up with island size, which we confirmed experimentally using laser ablation assays and immunofluorescence. Without formally being a 3-dimensional model, our approach has the minimal elements necessary to reproduce the distribution of cellular forces and mechanical crosstalk as well as distribution of principal stress in cells within epithelial cell aggregates. By mak
Discoveries of endogenous human stem cell-associated retroviruses (SCARs) revealed consistent activation of specific endogenous retroviral elements in human preimplantation embryos and documented the essential role of the sustained retroviral activities in the maintenance of pluripotency, functional identity and integrity of naive-state embryonic stem cells, and anti-viral resistance of the early-stage human embryos. Activation of specific SCARs, namely LTR7.HERVH and LTR5Hs.HERVK, has been demonstrated in patients diagnosed with multiple types of cancer, autoimmune diseases, neurodegenerative disorders and it is likely associated with the emergence of clinically lethal therapy resistant death-from-cancer phenotypes in a sub-set of cancer patients diagnosed with different types of malignant tumors.
Numerous biological approaches are available to characterise the mechanisms which govern the formation of human embryonic stem cell (hESC) colonies. To understand how the kinematics of single and pairs of hESCs impact colony formation, we study their mobility characteristics using time-lapse imaging. We perform a detailed statistical analysis of their speed, survival, directionality, distance travelled and diffusivity. We confirm that single and pairs of cells migrate as a diffusive random walk. Moreover, we show that the presence of Cell Tracer significantly reduces hESC mobility. Our results open the path to employ the theoretical framework of the diffusive random walk for the prognostic modelling and optimisation of the growth of hESC colonies. Indeed, we employ this random walk model to estimate the seeding density required to minimise the occurrence of hESC colonies arising from more than one founder cell and the minimal cell number needed for successful colony formation. We believe that our prognostic model can be extended to investigate the kinematic behaviour of somatic cells emerging from hESC differentiation and to enable its wide application in phenotyping of pluripotent
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 \(a_{dif}\).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 \(a_{dif}\) the tumor evolves towards a final state that contains a CSC fraction, but a strong the
The entropy production rate of cancer cell is always higher than healthy cell under the case of no external field applied. Different entropy production between two kinds of cells determines the direction of entropy flow among cells. The entropy flow is the carrier of information flow. The entropy flow from cancer to healthy cell takes along the harmful information of cancerous cell, propagating its toxic action to healthy tissues. We demonstrate that a low-frequency and low-intensity electromagnetic field or ultrasound irradiation may increase the entropy production rate of a cell in normal tissue than that in cancer, consequently reverse the direction of entropy current between two kinds of cells. The modification of PH value of cells may also cause the reversal of the direction of entropy flow between healthy and cancerous cells. So, the biological tissue under the irradiation of electromagnetic field or ultrasound or under the appropriate change of cell acidity can avoid the propagation of harmful information from cancer cells. We suggest that this entropy mechanism possibly provides a basis for a novel approach to anticancer therapy.
To explain the differentiation of stem cells in terms of dynamical systems theory, models of interacting cells with intracellular protein expression dynamics are analyzed and simulated. Simulations were carried out for all possible protein expression networks consisting of two genes under cell--cell interactions mediated by the diffusion of a protein. Networks that show cell differentiation are extracted and two forms of symmetric differentiation based on Turing's mechanism and asymmetric differentiation are identified. In the latter network, the intracellular protein levels show oscillatory dynamics at a single-cell level, while cell-to-cell synchronicity of the oscillation is lost with an increase in the number of cells. Differentiation to a fixed-point type behavior follows with a further increase in the number of cells. The cell type with oscillatory dynamics corresponds to a stem cell that can both proliferate and differentiate, while the latter fixed-point type only proliferates. This differentiation is analyzed as a saddle-node bifurcation on an invariant circle, while the number ratio of each cell type is shown to be robust against perturbations due to self-consistent deter