Feedback regulation in multistage cell lineages

Functional Classification of Fusion Proteins in Sarcoma

Sarcomas comprise a heterogeneous group of malignant tumors of mesenchymal origin. More than 80 entities are associated with different mesenchymal lineages. Sarcomas with fibroblastic, muscle, bone, vascular, adipocytic, and other characteristics are distinguished. Nearly half of all entities contain specific chromosomal translocations that give rise to fusion proteins. These are mostly pathognomonic, and their detection by various molecular techniques supports histopathologic classification. Moreover, the fusion proteins act as oncogenic drivers, and their blockade represents a promising therapeutic approach. This review summarizes the current knowledge on fusion proteins in sarcoma. We categorize the different fusion proteins into functional classes, including kinases, epigenetic regulators, and transcription factors, and describe their mechanisms of action. Interestingly, while fusion proteins acting as transcription factors are found in all mesenchymal lineages, the others have a more restricted pattern. Most kinase-driven sarcomas belong to the fibroblastic/myofibroblastic lineage. Fusion proteins with an epigenetic function are mainly associated with sarcomas of unclear differentiation, suggesting that epigenetic dysregulation leads to a major change in cell identity. Comparison of mechanisms of action reveals recurrent functional modes, including antagonism of Polycomb activity by fusion proteins with epigenetic activity and recruitment of histone acetyltransferases by fusion transcription factors of the myogenic lineage. Finally, based on their biology, we describe potential approaches to block the activity of fusion proteins for therapeutic intervention. Overall, our work highlights differences as well as similarities in the biology of fusion proteins from different sarcomas and provides the basis for a functional classification.

The inherent fragility of collective proliferative control

Tissues achieve and maintain their sizes through active feedback, whereby cells collectively regulate proliferation and differentiation so as to facilitate homeostasis and the ability to respond to disturbances. One of the best understood feedback mechanisms-renewal control-achieves remarkable feats of robustness in determining and maintaining desired sizes. Yet in a variety of biologically relevant situations, we show that stochastic effects should cause rare but catastrophic failures of renewal control. We define the circumstances under which this occurs and raise the possibility such events account for important non-genetic steps in the development of cancer. We further suggest that the spontaneous stochastic reversal of these events could explain cases of cancer normalization or dormancy following treatment. Indeed, we show that the kinetics of post-treatment recurrence for many cancers are often better fit by a model of stochastic re-emergence due to loss of collective proliferative control, than by deterministic models of cancer relapse.

Stochastic model for cell population dynamics quantifies homeostasis in colonic crypts and its disruption in early tumorigenesis

The questions of how healthy colonic crypts maintain their size, and how homeostasis is disrupted by driver mutations, are central to understanding colorectal tumorigenesis. We propose a three-type stochastic branching process, which accounts for stem, transit-amplifying (TA) and fully differentiated (FD) cells, to model the dynamics of cell populations residing in colonic crypts. Our model is simple in its formulation, allowing us to estimate all but one of the model parameters from the literature. Fitting the single remaining parameter, we find that model results agree well with data from healthy human colonic crypts, capturing the considerable variance in population sizes observed experimentally. Importantly, our model predicts a steady-state population in healthy colonic crypts for relevant parameter values. We show that APC and KRAS mutations, the most significant early alterations leading to colorectal cancer, result in increased steady-state populations in mutated crypts, in agreement with experimental results. Finally, our model predicts a simple condition for unbounded growth of cells in a crypt, corresponding to colorectal malignancy. This is predicted to occur when the division rate of TA cells exceeds their differentiation rate, with implications for therapeutic cancer prevention strategies.

A multiscale chemical-mechanical model predicts impact of morphogen spreading on tissue growth

The exact mechanism controlling cell growth remains a grand challenge in developmental biology and regenerative medicine. The Drosophila wing disc tissue serves as an ideal biological model to study mechanisms involved in growth regulation. Most existing computational models for studying tissue growth focus specifically on either chemical signals or mechanical forces. Here we developed a multiscale chemical-mechanical model to investigate the growth regulation mechanism based on the dynamics of a morphogen gradient. By comparing the spatial distribution of dividing cells and the overall tissue shape obtained in model simulations with experimental data of the wing disc, it is shown that the size of the domain of the Dpp morphogen is critical in determining tissue size and shape. A larger tissue size with a faster growth rate and more symmetric shape can be achieved if the Dpp gradient spreads in a larger domain. Together with Dpp absorbance at the peripheral zone, the feedback regulation that downregulates Dpp receptors on the cell membrane allows for further spreading of the morphogen away from its source region, resulting in prolonged tissue growth at a more spatially homogeneous growth rate.

Mathematics of neural stem cells: Linking data and processes

Adult stem cells are described as a discrete population of cells that stand at the top of a hierarchy of progressively differentiating cells. Through their unique ability to self-renew and differentiate, they regulate the number of end-differentiated cells that contribute to tissue physiology. The question of how discrete, continuous, or reversible the transitions through these hierarchies are and the precise parameters that determine the ultimate performance of stem cells in adulthood are the subject of intense research. In this review, we explain how mathematical modelling has improved the mechanistic understanding of stem cell dynamics in the adult brain. We also discuss how single-cell sequencing has influenced the understanding of cell states or cell types. Finally, we discuss how the combination of single-cell sequencing technologies and mathematical modelling provides a unique opportunity to answer some burning questions in the field of stem cell biology.

Modeling of ionizing radiation induced hair follicle regenerative dynamics

Hair follicles (HFs) are stem-cell-rich mammalian mini-organs that can undergo cyclic regenerations over the life span of the organism. The cycle of a HF consists of three consecutive phases: anagen-the active proliferation phase, catagen-the degeneration phase, and telogen-the resting phase. While HFs undergo irreversible degeneration during catagen, recent experimental research on mice shows that when anagen HFs are subject to ionizing radiation (IR), they undergo a transient degeneration, followed by a nearly full regeneration that makes the HFs return to homeostatic state. The mechanisms underlying these IR-induced HF regenerative dynamics and the catagen degenerative dynamics, remain unknown. In this work, we develop an ODE type cell differentiation population model to study the control mechanisms of HF regeneration. The model is built based on current theoretical knowledge in biology and mathematically formulated using feedback mechanisms. Model parameters are calibrated to IR experimental data, and we then provide modeling results with both deterministic ODE simulations and corresponding stochastic simulations. We perform stability and bifurcation analyses on the ODE model, which reveal that for anagen HFs, a low spontaneous apoptosis rate secures the stability of the HF homeostatic steady state, allowing the HF to regenerate even when subject to strong IR. On the other hand, the irreversible degeneration during catagen results from both strong spontaneous apoptosis rate and strong apoptosis feedback. Lastly, we perform sensitivity analysis to identify key parameters in the model to validate these hypotheses.

Regulation of stem cell dynamics through volume exclusion

The maintenance and regeneration of adult tissues rely on the self-renewal of stem cells. Regeneration without over-proliferation requires precise regulation of the stem cell proliferation and differentiation rates. The nature of such regulatory mechanisms in different tissues, and how to incorporate them in models of stem cell population dynamics, is incompletely understood. The critical birth-death (CBD) process is widely used to model stem cell populations, capturing key phenomena, such as scaling laws in clone size distributions. However, the CBD process neglects regulatory mechanisms. Here, we propose the birth-death process with volume exclusion (vBD), a variation of the birth-death process that considers crowding effects, such as may arise due to limited space in a stem cell niche. While the deterministic rate equations predict a single non-trivial attracting steady state, the master equation predicts extinction and transient distributions of stem cell numbers with three possible behaviours: long-lived quasi-steady state (QSS), and short-lived bimodal or unimodal distributions. In all cases, we approximate solutions to the vBD master equation using a renormalized system-size expansion, QSS approximation and the Wentzel–Kramers–Brillouin method. Our study suggests that the size distribution of a stem cell population bears signatures that are useful to detect negative feedback mediated via volume exclusion.

Microenvironmental sensing by fibroblasts controls macrophage population size

Animal tissues comprise diverse cell types. However, the mechanisms controlling the number of each cell type within tissue compartments remain poorly understood. Here, we report that different cell types utilize distinct strategies to control population numbers. Proliferation of fibroblasts, stromal cells important for tissue integrity, is limited by space availability. In contrast, proliferation of macrophages, innate immune cells involved in defense, repair, and homeostasis, is constrained by growth factor availability. Examination of density-dependent gene expression in fibroblasts revealed that Hippo and TGF-β target genes are both regulated by cell density. We found YAP1, the transcriptional coactivator of the Hippo signaling pathway, directly regulates expression of Csf1, the lineage-specific growth factor for macrophages, through an enhancer of Csf1 that is specifically active in fibroblasts. Activation of YAP1 in fibroblasts elevates Csf1 expression and is sufficient to increase the number of macrophages at steady state. Our data also suggest that expression programs in fibroblasts that change with density may result from sensing of mechanical force through actin-dependent mechanisms. Altogether, we demonstrate that two different modes of population control are connected and coordinated to regulate cell numbers of distinct cell types. Sensing of the tissue environment may serve as a general strategy to control tissue composition.

Effect of dedifferentiation on noise propagation in cellular hierarchy

Many fast renewing tissues have a hierarchical structure. Tissue-specific stem cells are at the root of this cellular hierarchy, which give rive to a whole range of specialized cells via cellular differentiation. However, increasing evidence shows that the hierarchical structure can be broken due to cellular dedifferentiation in which cells at differentiated stages can revert to the stem cell stage. Dedifferentiation has significant impacts on many aspects of hierarchical tissues. Here we investigate the effect of dedifferentiation on noise propagation by developing a stochastic model composed of different cell types. The moment equations are derived, via which we systematically investigate how the noise in the cell number is changed by dedifferentiation. Our results suggest that dedifferentiation have different effects on the noises in the numbers of stem cells and nonstem cells. Specifically, the noise in the number of stem cells is significantly reduced by increasing dedifferentiation probability. Due to the dual effect of dedifferentiation on nonstem cells, however, more complex changes could happen to the noise in the number of nonstem cells by increasing dedifferentiation probability. Furthermore, it is found that even though dedifferentiation could turn part of the noise propagation process into a noise-amplifying step, it is very unlikely to turn the entire process into a noise-amplifying cascade.

Spatial dynamics of feedback and feedforward regulation in cell lineages

Feedback mechanisms within cell lineages are thought to be important for maintaining tissue homeostasis. Mathematical models that assume well-mixed cell populations, together with experimental data, have suggested that negative feedback from differentiated cells on the stem cell self-renewal probability can maintain a stable equilibrium and hence homeostasis. Cell lineage dynamics, however, are characterized by spatial structure, which can lead to different properties. Here, we investigate these dynamics using spatially explicit computational models, including cell division, differentiation, death, and migration / diffusion processes. According to these models, the negative feedback loop on stem cell self-renewal fails to maintain homeostasis, both under the assumption of strong spatial restrictions and fast migration / diffusion. Although homeostasis cannot be maintained, this feedback can regulate cell density and promote the formation of spatial structures in the model. Tissue homeostasis, however, can be achieved if spatially restricted negative feedback on self-renewal is combined with an experimentally documented spatial feedforward loop, in which stem cells regulate the fate of transit amplifying cells. This indicates that the dynamics of feedback regulation in tissue cell lineages are more complex than previously thought, and that combinations of spatially explicit control mechanisms are likely instrumental.

Mathematical modelling identifies conditions for maintaining and escaping feedback control in the intestinal epithelium

The intestinal epithelium is one of the fastest renewing tissues in mammals. It shows a hierarchical organisation, where intestinal stem cells at the base of crypts give rise to rapidly dividing transit amplifying cells that in turn renew the pool of short-lived differentiated cells. Upon injury and stem-cell loss, cells can also de-differentiate. Tissue homeostasis requires a tightly regulated balance of differentiation and stem cell proliferation, and failure can lead to tissue extinction or to unbounded growth and cancerous lesions. Here, we present a two-compartment mathematical model of intestinal epithelium population dynamics that includes a known feedback inhibition of stem cell differentiation by differentiated cells. The model shows that feedback regulation stabilises the number of differentiated cells as these become invariant to changes in their apoptosis rate. Stability of the system is largely independent of feedback strength and shape, but specific thresholds exist which if bypassed cause unbounded growth. When dedifferentiation is added to the model, we find that the system can recover faster after certain external perturbations. However, dedifferentiation makes the system more prone to losing homeostasis. Taken together, our mathematical model shows how a feedback-controlled hierarchical tissue can maintain homeostasis and can be robust to many external perturbations.

Distinguishing between multiple mathematical models of neural stem cell quiescence and activation during age-related neural stem cell decline in neurogenesis

Stem cells are required for tissue maintenance and homeostasis during an organism's lifetime. Neural stem cells (NSCs) can be in an actively dividing state or in a quiescent state. The balance between stem cell quiescence and cycling activity determines the rate of neurogenesis. With age, more NSCs enter the quiescent state, while the total number of NSCs decreases. Here we reconsider an existing mathematical model of how neural stem cells switch between active and quiescent states from the point of view of control theory by considering the activation rate, self-renewal probability, and division rate as control parameters rather than as pre-defined functions. Our goal is to test whether those modifications to the basic model could explain the observed decline of neural stem cells with age better than Gomerzian time-dependent parameters, and compare the output from different model variants to experimental data from mice using AIC. We find that time-dependent activation rate provides the best fit to the activated cell fraction (ACF) of NSCs over time, but that other model variants with constant parameter values can better fit the total number of NSCs over time. We also consider an alternate model for NSCs with nonlinear feedback from progenitor cells that affect NSC parameters, and compare all models to experimental stem cell and progenitor data. However, all of the feedback models considered provide a worse fit to the experimental data. This suggests that when switching between active and quiescent stem cells is considered, a time-dependent linear model outperforms the integral feedback mechanism considered by other models of stem cell lineages. Fitting progenitor data for both the time varying models and feedback models indicates that four or five intermediate transit amplifying progenitor states are necessary. Our modeling suggests that in order to determine whether an increase in age-related neural stem cell quiescence is determined by by a decreasing stem cell activation rate or an increased stem cell depletion rate, additional experiments should be designed to explore whether or not depletion of the stem cell pool is occurring, and that a higher resolution time series for activated cell fraction (ACF) would be best to resolve this issue.

Regulatory feedback effects on tissue growth dynamics in a two-stage cell lineage model

Identifying the mechanism of intercellular feedback regulation is critical for the basic understanding of tissue growth control in organisms. In this paper, we analyze a tissue growth model consisting of a single lineage of two cell types regulated by negative feedback signaling molecules that undergo spatial diffusion. By deriving the fixed points for the uniform steady states and carrying out linear stability analysis, phase diagrams are obtained analytically for arbitrary parameters of the model. Two different generic growth modes are found: blow-up growth and final-state controlled growth which are governed by the nontrivial fixed point and the trivial fixed point, respectively, and can be sensitively switched by varying the negative feedback regulation on the proliferation of the stem cells. Analytic expressions for the characteristic timescales for these two growth modes are also derived. Remarkably, the trivial and nontrivial uniform steady states can coexist and a sharp transition occurs in the bistable regime as the relevant parameters are varied. Furthermore, the bistable growth properties allows for the external control to switch between these two growth modes. In addition, the condition for an early accelerated growth followed by a retarded growth can be derived. These analytical results are further verified by numerical simulations and provide insights on the growth behavior of the tissue. Our results are also discussed in the light of possible realistic biological experiments and tissue growth control strategy. Furthermore, by external feedback control of the concentration of regulatory molecules, it is possible to achieve a desired growth mode, as demonstrated with an analysis of boosted growth, catch-up growth and the design for the target of a linear growth dynamic.

Evolution of cancer stem cell lineage involving feedback regulation

Tumor emergence and progression is a complex phenomenon that assumes special molecular and cellular interactions. The hierarchical structuring and communication via feedback signaling of different cell types, which are categorized as the stem, progenitor, and differentiated cells in dependence of their maturity level, plays an important role. Under healthy conditions, these cells build a dynamical system that is responsible for facilitating the homeostatic regulation of the tissue. Generally, in this hierarchical setting, stem and progenitor cells are yet likely to undergo a mutation, when a cell divides into two daughter cells. This may lead to the development of abnormal characteristics, i.e. mutation in the cell, yielding an unrestrained number of cells. Therefore, the regulation of a stem cell’s proliferation and differentiation rate is crucial for maintaining the balance in the overall cell population. In this paper, a maturity based mathematical model with feedback regulation is formulated for healthy and mutated cell lineages. It is given in the form of coupled ordinary and partial differential equations. The focus is laid on the dynamical effects resulting from acquiring a mutation in the hierarchical structure of stem, progenitor and fully differentiated cells. Additionally, the effects of nonlinear feedback regulation from mature cells into both stem and progenitor cell populations have been inspected. The steady-state solutions of the model are derived analytically. Numerical simulations and results based on a finite volume scheme underpin various expected behavioral patterns of the homeostatic regulation and cancer evolution. For instance, it has been found that the mutated cells can experience significant growth even with a single somatic mutation, but under homeostatic regulation acquire a steady-state and thus, ensuing healthy cell population to either a steady-state or a lower cell concentration. Furthermore, the model behavior has been validated with different experimentally measured tumor values from the literature.

Multistage feedback-driven compartmental dynamics of hematopoiesis

Human hematopoiesis is surprisingly resilient to disruptions, providing suitable responses to severe bleeding, long-lasting immune activation, and even bone marrow transplants. Still, many blood disorders exist which push the system past its natural plasticity, resulting in abnormalities in the circulating blood. While proper treatment of such diseases can benefit from understanding the underlying cell dynamics, these are non-trivial to predict due to the hematopoietic system's hierarchical nature and complex feedback networks. To characterize the dynamics following different types of perturbations, we investigate a model representing hematopoiesis as a sequence of compartments covering all maturation stages-from stem to mature cells-where feedback regulates cell production to ongoing necessities. We find that a stable response to perturbations requires the simultaneous adaptation of cell differentiation and self-renewal rates, and show that under conditions of continuous disruption-as found in chronic hemolytic states-compartment cell numbers evolve to novel stable states.

CELLULAR FEEDBACK NETWORKS AND THEIR RESILIENCE AGAINST MUTATIONS

Many tissues undergo a steady turnover, where cell divisions are on average balanced with cell deaths. Cell fate decisions such as stem cell (SC) differentiations, proliferations, or differentiated cell (DC) deaths, may be controlled by cell populations through cell-to-cell signaling. Here, we examine a class of mathematical models of turnover in SC lineages to understand engineering design principles of control (feedback) loops, that may operate in such systems. By using ordinary differential equations that describe the co-dynamics of SCs and DCs, we study the effect of different types of mutations that interfere with feedback present within cellular networks. For instance, we find that mutants that do not participate in feedback are less dangerous in the sense that they will not rise from low numbers, whereas mutants that do not respond to feedback signals could rise and replace the wild-type population. Additionally, we asked if different feedback networks can have different degrees of resilience against such mutations. We found that all minimal networks, that is networks consisting of exactly one feedback loop that is sufficient for homeostatic stability of the wild-type population, are equally vulnerable. Mutants with a weakened/eliminated feedback parameter might expand from lower numbers and either enter unlimited growth or reach an equilibrium with an increased number of SCs and DCs. Therefore, from an evolutionary viewpoint, it appears advantageous to combine feedback loops, creating redundant feedback networks. Interestingly, from an engineering prospective, not all such redundant systems are equally resilient. For some of them, any mutation that weakens/eliminates one of the loops will lead to a population growth of SCs. For others, the population of SCs can actually shrink as a result of “cutting” one of the loops, thus slowing down further unwanted transformations.

Inference of Intercellular Communications and Multilayer Gene-Regulations of Epithelial-Mesenchymal Transition From Single-Cell Transcriptomic Data

Epithelial-to-mesenchymal transition (EMT) plays an important role in many biological processes during development and cancer. The advent of single-cell transcriptome sequencing techniques allows the dissection of dynamical details underlying EMT with unprecedented resolution. Despite several single-cell data analysis on EMT, how cell communicates and regulates dynamics along the EMT trajectory remains elusive. Using single-cell transcriptomic datasets, here we infer the cell-cell communications and the multilayer gene-gene regulation networks to analyze and visualize the complex cellular crosstalk and the underlying gene regulatory dynamics along EMT. Combining with trajectory analysis, our approach reveals the existence of multiple intermediate cell states (ICSs) with hybrid epithelial and mesenchymal features. Analyses on the time-series datasets from cancer cell lines with different inducing factors show that the induced EMTs are context-specific: the EMT induced by transforming growth factor B1 (TGFB1) is synchronous, whereas the EMTs induced by epidermal growth factor and tumor necrosis factor are asynchronous, and the responses of TGF-β pathway in terms of gene expression regulations are heterogeneous under different treatments or among various cell states. Meanwhile, network topology analysis suggests that the ICSs during EMT serve as the signaling in cellular communication under different conditions. Interestingly, our analysis of a mouse skin squamous cell carcinoma dataset also suggests regardless of the significant discrepancy in concrete genes between in vitro and in vivo EMT systems, the ICSs play dominant role in the TGF-β signaling crosstalk. Overall, our approach reveals the multiscale mechanisms coupling cell-cell communications and gene-gene regulations responsible for complex cell-state transitions.

Modelling stem cell ageing: a multi-compartment continuum approach

Stem cells are important to generate all specialized tissues at an early life stage, and in some systems, they also have repair functions to replenish the adult tissues. Repeated cell divisions lead to the accumulation of molecular damage in stem cells, which are commonly recognized as drivers of ageing. In this paper, a novel model is proposed to integrate stem cell proliferation and differentiation with damage accumulation in the stem cell ageing process. A system of two structured PDEs is used to model the population densities of stem cells (including all multiple progenitors) and terminally differentiated (TD) cells. In this system, cell cycle progression and damage accumulation are modelled by continuous dynamics, and damage segregation between daughter cells is considered at each division. Analysis and numerical simulations are conducted to study the steady-state populations and stem cell damage distributions under different damage segregation strategies. Our simulations suggest that equal distribution of the damaging substance between stem cells in a symmetric renewal and less damage retention in stem cells in the asymmetric division are favourable strategies, which reduce the death rate of the stem cells and increase the TD cell populations. Moreover, asymmetric damage segregation in stem cells leads to less concentrated damage distribution in the stem cell population, which may be more robust to the stochastic changes in the damage. The feedback regulation from stem cells can reduce oscillations and population overshoot in the process, and improve the fitness of stem cells by increasing the percentage of cells with less damage in the stem cell population.

STOCHASTIC DYNAMICS OF CELL LINEAGE IN TISSUE HOMEOSTASIS

During epithelium tissue maintenance, lineages of cells differentiate and proliferate in a coordinated way to provide the desirable size and spatial organization of different types of cells. While mathematical models through deterministic description have been used to dissect role of feedback regulations on tissue layer size and stratification, how the stochastic effects influence tissue maintenance remains largely unknown. Here we present a stochastic continuum model for cell lineages to investigate how both layer thickness and layer stratification are affected by noise. We find that the cell-intrinsic noise often causes reduction and oscillation of layer size whereas the cell-extrinsic noise increases the thickness, and sometimes, leads to uncontrollable growth of the tissue layer. The layer stratification usually deteriorates as the noise level increases in the cell lineage systems. Interestingly, the morphogen noise, which mixes both cell-intrinsic noise and cell-extrinsic noise, can lead to larger size of layer with little impact on the layer stratification. By investigating different combinations of the three types of noise, we find the layer thickness variability is reduced when cell-extrinsic noise level is high or morphogen noise level is low. Interestingly, there exists a tradeoff between low thickness variability and strong layer stratification due to competition among the three types of noise, suggesting robust layer homeostasis requires balanced levels of different types of noise in the cell lineage systems.

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.

Modelling acute myeloid leukaemia in a continuum of differentiation states

Here we present a mathematical model of movement in an abstract space representing states of cellular differentiation. We motivate this work with recent examples that demonstrate a continuum of cellular differentiation using single cell RNA sequencing data to characterize cellular states in a high-dimensional space, which is then mapped into ℝ 2 or ℝ 2 with dimension reduction techniques. We represent trajectories in the differentiation space as a graph, and model directed and random movement on the graph with partial differential equations. We hypothesize that flow in this space can be used to model normal and abnormal differentiation processes. We present a mathematical model of hematopoeisis parameterized with publicly available single cell RNA-Seq data and use it to simulate the pathogenesis of acute myeloid leukemia (AML). The model predicts the emergence of cells in novel intermediate states of differentiation consistent with immunophenotypic characterizations of a mouse model of AML.

Feedback regulation in a stem cell model with acute myeloid leukaemia

Background

The haematopoietic lineages with leukaemia lineages are considered in this paper. In particular, we mainly consider that haematopoietic lineages are tightly controlled by negative feedback inhibition of end-product. Actually, leukemia has been found 100 years ago. Up to now, the exact mechanism is still unknown, and many factors are thought to be associated with the pathogenesis of leukemia. Nevertheless, it is very necessary to continue the profound study of the pathogenesis of leukemia. Here, we propose a new mathematical model which include some negative feedback inhibition from the terminally differentiated cells of haematopoietic lineages to the haematopoietic stem cells and haematopoietic progenitor cells in order to describe the regulatory mechanisms mentioned above by a set of ordinary differential equations. Afterwards, we carried out detailed dynamical bifurcation analysis of the model, and obtained some meaningful results.

Results

In this work, we mainly perform the analysis of the mathematic model by bifurcation theory and numerical simulations. We have not only incorporated some new negative feedback mechanisms to the existing model, but also constructed our own model by using the modeling method of stem cell theory with probability method. Through a series of qualitative analysis and numerical simulations, we obtain that the weak negative feedback for differentiation probability is conducive to the cure of leukemia. However, with the strengthening of negative feedback, leukemia will be more difficult to be cured, and even induce death. In contrast, strong negative feedback for differentiation rate of progenitor cells can promote healthy haematopoiesis and suppress leukaemia.

Conclusions

These results demonstrate that healthy progenitor cells are bestowed a competitive advantage over leukaemia stem cells. Weak g₁, g₂, and h₁ enable the system stays in the healthy state. However, strong h₂ can promote healthy haematopoiesis and suppress leukaemia.

Multiscale modeling of layer formation in epidermis

The mammalian skin epidermis is a stratified epithelium composed of multiple layers of epithelial cells that exist in appropriate sizes and proportions, and with distinct boundaries separating each other. How the epidermis develops from a single layer of committed precursor cells to form a complex multilayered structure of multiple cell types remains elusive. Here, we construct stochastic, three-dimensional, and multiscale models consisting of a lineage of multiple cell types to study the control of epidermal development. Symmetric and asymmetric cell divisions, stochastic cell fate transitions within the lineage, extracellular morphogens, cell-to-cell adhesion forces, and cell signaling are included in model. A GPU algorithm was developed and implemented to accelerate the simulations. These simulations show that a balance between cell proliferation and differentiation during lineage progression is crucial for the development and maintenance of the epidermal tissue. We also find that selective intercellular adhesion is critical to sharpening the boundary between layers and to the formation of a highly ordered structure. The long-range action of a morphogen provides additional feedback regulations, enhancing the robustness of overall layer formation. Our model is built upon previous experimental findings revealing the role of Ovol transcription factors in regulating epidermal development. Direct comparisons of experimental and simulation perturbations show remarkable consistency. Taken together, our results highlight the major determinants of a well-stratified epidermis: balanced proliferation and differentiation, and a combination of both short- (symmetric/asymmetric division and selective cell adhesion) and long-range (morphogen) regulations. These underlying principles have broad implications for other developmental or regenerative processes leading to the formation of multilayered tissue structures, as well as for pathological processes such as epidermal wound healing.

Epidermal morphogenesis, which occurs during the second half of embryogenesis, is the developmental process that generates a skin permeability barrier essential for terrestrial survival. Defects with this barrier are associated with common skin disorders such as atopic dermatitis. Study of mechanisms that control epidermal development and differentiation is therefore highly relevant to human health. Motivated by recent experimental observations on the role of Ovol transcription factors in regulating epidermal development, we developed a multiscale model to investigate the underlying mechanisms responsible for epidermal layer formation and homeostasis. We report that regulation of proliferation and differentiation by Ovol plays an important role in epidermal development. In addition, our computational analysis shows that asymmetric cell division, selective cell adhesion, and morphogen regulation work in a synergetic manner to produce the well-stratified epidermal layers. Taken together, our results demonstrate that robust epidermal morphogenesis involves a balance between proliferation and differentiation, and an interplay between short- and long-range spatial control mechanisms. This principle may also be applicable to other complex systems of tissue development or regeneration.

Control of cell fraction and population recovery during tissue regeneration in stem cell lineages

Multicellular tissues are continually turning over, and homeostasis is maintained through regulated proliferation and differentiation of stem cells and progenitors. Following tissue injury, a dramatic increase in cell proliferation is commonly observed, resulting in rapid restoration of tissue size. This regulation is thought to occur via multiple feedback loops acting on cell self-renewal or differentiation. Models of ordinary differential equations have been widely used to study the cell lineage system. Prior modeling studies have suggested that loss of homeostasis and initiation of tumorigenesis can be contributed to the loss of control of these processes, and the rate of symmetric versus asymmetric division of the stem cells may also be altered. While most of the previous works focused on analysis of stability, existence and uniqueness of steady states of multistage cell lineage models, in this work we attempt to understand the cell lineage model from a different perspective. We compare three variants of hierarchical stem cell lineage tissue models with different combinations of negative feedbacks and use sensitivity analysis to examine the possible strategies for the cells to achieve certain performance objectives. Our results suggest that multiple negative feedback loops must be present in the stem cell lineage to keep the fractions of stem cells to differentiated cells in the total population as robust as possible to variations in cell division parameters, and to minimize the time for tissue recovery in a non-oscillatory manner.

Circuit Design Features of a Stable Two-Cell System

Cell communication within tissues is mediated by multiple paracrine signals including growth factors, which control cell survival and proliferation. Cells and the growth factors they produce and receive constitute a circuit with specific properties that ensure homeostasis. Here, we used computational and experimental approaches to characterize the features of cell circuits based on growth factor exchange between macrophages and fibroblasts, two cell types found in most mammalian tissues. We found that the macrophage-fibroblast cell circuit is stable and robust to perturbations. Analytical screening of all possible two-cell circuit topologies revealed the circuit features sufficient for stability, including environmental constraint and negative-feedback regulation. Moreover, we found that cell-cell contact is essential for the stability of the macrophage-fibroblast circuit. These findings illustrate principles of cell circuit design and provide a quantitative perspective on cell interactions.

Modeling the behavior of human induced pluripotent stem cells seeded on melt electrospun scaffolds

BACKGROUND

Human induced pluripotent stem cells (hiPSCs) can form any tissue found in the body, making them attractive for regenerative medicine applications. Seeding hiPSC aggregates into biomaterial scaffolds can control their differentiation into specific tissue types. Here we develop and analyze a mathematical model of hiPSC aggregate behavior when seeded on melt electrospun scaffolds with defined topography.

RESULTS

We used ordinary differential equations to model the different cellular populations (stem, progenitor, differentiated) present in our scaffolds based on experimental results and published literature. Our model successfully captures qualitative features of the cellular dynamics observed experimentally. We determined the optimal parameter sets to maximize specific cellular populations experimentally, showing that a physiologic oxygen level (∼ 5%) increases the number of neural progenitors and differentiated neurons compared to atmospheric oxygen levels (∼ 21%) and a scaffold porosity of ∼ 63% maximizes aggregate size.

CONCLUSIONS

Our mathematical model determined the key factors controlling hiPSC behavior on melt electrospun scaffolds, enabling optimization of experimental parameters.

Determining the control networks regulating stem cell lineages in colonic crypts

The question of stem cell control is at the center of our understanding of tissue functioning, both in healthy and cancerous conditions. It is well accepted that cellular fate decisions (such as divisions, differentiation, apoptosis) are orchestrated by a network of regulatory signals emitted by different cell populations in the lineage and the surrounding tissue. The exact regulatory network that governs stem cell lineages in a given tissue is usually unknown. Here we propose an algorithm to identify a set of candidate control networks that are compatible with (a) measured means and variances of cell populations in different compartments, (b) qualitative information on cell population dynamics, such as the existence of local controls and oscillatory reaction of the system to population size perturbations, and (c) statistics of correlations between cell numbers in different compartments. Using the example of human colon crypts, where lineages are comprised of stem cells, transit amplifying cells, and differentiated cells, we start with a theoretically known set of 32 smallest control networks compatible with tissue stability. Utilizing near-equilibrium stochastic calculus of stem cells developed earlier, we apply a series of tests, where we compare the networks' expected behavior with the observations. This allows us to exclude most of the networks, until only three, very similar, candidate networks remain, which are most compatible with the measurements. This work demonstrates how theoretical analysis of control networks combined with only static biological data can shed light onto the inner workings of stem cell lineages, in the absence of direct experimental assessment of regulatory signaling mechanisms. The resulting candidate networks are dominated by negative control loops and possess the following properties: (1) stem cell division decisions are negatively controlled by the stem cell population, (2) stem cell differentiation decisions are negatively controlled by the transit amplifying cell population.

Cellular Hierarchy as a Determinant of Tumor Sensitivity to Chemotherapy

Chemotherapy has been shown to enrich cancer stem cells in tumors. Recently, we demonstrated that administration of chemotherapy to human bladder cancer xenografts could trigger a wound-healing response that mobilizes quiescent tumor stem cells into active proliferation. This phenomenon leads to a loss of sensitivity to chemotherapy partly due to an increase in the number of tumor stem cells, which typically respond to chemotherapy-induced cell death less than more differentiated cells. Different bladder cancer xenografts, however, demonstrate differential sensitivities to chemotherapy, the basis of which is not understood. Using mathematical models, we show that characteristics of the tumor cell hierarchy can be crucial for determining the sensitivity of tumors to drug therapy, under the assumption that stem cell enrichment is the primary basis for drug resistance. Intriguingly, our model predicts a weaker response to therapy if there is negative feedback from differentiated tumor cells that inhibits the rate of tumor stem cell division. If this negative feedback is less pronounced, the treatment response is predicted to be enhanced. The reason is that negative feedback on the rate of tumor cell division promotes a permanent rise of the tumor stem cell population over time, both in the absence of treatment and even more so during drug therapy. Model application to data from chemotherapy-treated patient-derived xenografts indicates support for model predictions. These findings call for further research into feedback mechanisms that might remain active in cancers and potentially highlight the presence of feedback as an indication to combine chemotherapy with approaches that limit the process of tumor stem cell enrichment. Cancer Res; 77(9); 2231-41. ©2017 AACR.

Controlling Stochasticity in Epithelial-Mesenchymal Transition Through Multiple Intermediate Cellular States

Epithelial-mesenchymal transition (EMT) is an instance of cellular plasticity that plays critical roles in development, regeneration and cancer progression. Recent studies indicate that the transition between epithelial and mesenchymal states is a multi-step and reversible process in which several intermediate phenotypes might coexist. These intermediate states correspond to various forms of stem-like cells in the EMT system, but the function of the multi-step transition or the multiple stem cell phenotypes is unclear. Here, we use mathematical models to show that multiple intermediate phenotypes in the EMT system help to attenuate the overall fluctuations of the cell population in terms of phenotypic compositions, thereby stabilizing a heterogeneous cell population in the EMT spectrum. We found that the ability of the system to attenuate noise on the intermediate states depends on the number of intermediate states, indicating the stem-cell population is more stable when it has more sub-states. Our study reveals a novel advantage of multiple intermediate EMT phenotypes in terms of systems design, and it sheds light on the general design principle of heterogeneous stem cell population.

Near Equilibrium Calculus of Stem Cells in Application to the Airway Epithelium Lineage

Homeostatic maintenance of tissues is orchestrated by well tuned networks of cellular signaling. Such networks regulate, in a stochastic manner, fates of all cells within the respective lineages. Processes such as symmetric and asymmetric divisions, differentiation, de-differentiation, and death have to be controlled in a dynamic fashion, such that the cell population is maintained at a stable equilibrium, has a sufficiently low level of stochastic variation, and is capable of responding efficiently to external damage. Cellular lineages in real tissues may consist of a number of different cell types, connected by hierarchical relationships, albeit not necessarily linear, and engaged in a number of different processes. Here we develop a general mathematical methodology for near equilibrium studies of arbitrarily complex hierarchical cell populations, under regulation by a control network. This methodology allows us to (1) determine stability properties of the network, (2) calculate the stochastic variance, and (3) predict how different control mechanisms affect stability and robustness of the system. We demonstrate the versatility of this tool by using the example of the airway epithelium lineage. Recent research shows that airway epithelium stem cells divide mostly asymmetrically, while the so-called secretory cells divide predominantly symmetrically. It further provides quantitative data on the recovery dynamics of the airway epithelium, which can include secretory cell de-differentiation. Using our new methodology, we demonstrate that while a number of regulatory networks can be compatible with the observed recovery behavior, the observed division patterns of cells are the most optimal from the viewpoint of homeostatic lineage stability and minimizing the variation of the cell population size. This not only explains the observed yet poorly understood features of airway tissue architecture, but also helps to deduce the information on the still largely hypothetical regulatory mechanisms governing tissue turnover, and lends insight into how different control loops influence the stability and variance properties of cell populations.

Tissue stability is the basic property of healthy organs, and yet the mechanisms governing the stable, long-term maintenance of cell numbers in tissues are poorly understood. While more and more signaling pathways are being discovered, for the most part it remains unknown how they are being put together by different cell types into complex, nonlinear, hierarchical control networks that, on the one hand, reliably maintain constant cell numbers, and on the other hand, quickly adjust to oversee the robust response to tissue damage. Theoretical approaches can fill the gap by being able to reconstruct the underlying control network, based on the observations about the aspects of cellular dynamics. We argue that while many hypothetical networks may be capable of basic cell lineage maintenance, some are much more efficient from the viewpoint of variance minimization. Thus, we developed a new methodology that can test various control networks for stability, variance, and robustness. In the example of the airway epithelium that we highlight, it turns out that the evolutionary selected, actual architecture coincides with the mathematically optimal solution that minimizes the fluctuations of cell numbers at homeostasis.

Feedback, Lineages and Self-Organizing Morphogenesis

Feedback regulation of cell lineage progression plays an important role in tissue size homeostasis, but whether such feedback also plays an important role in tissue morphogenesis has yet to be explored. Here we use mathematical modeling to show that a particular feedback architecture in which both positive and negative diffusible signals act on stem and/or progenitor cells leads to the appearance of bistable or bi-modal growth behaviors, ultrasensitivity to external growth cues, local growth-driven budding, self-sustaining elongation, and the triggering of self-organization in the form of lamellar fingers. Such behaviors arise not through regulation of cell cycle speeds, but through the control of stem or progenitor self-renewal. Even though the spatial patterns that arise in this setting are the result of interactions between diffusible factors with antagonistic effects, morphogenesis is not the consequence of Turing-type instabilities.

Regulatory effects on the population dynamics and wave propagation in a cell lineage model

We consider the interplay of cell proliferation, cell differentiation (and de-differentiation), cell movement, and the effect of feedback regulations on the population and propagation dynamics of different cell types in a cell lineage model. Cells are assumed to secrete and respond to negative feedback molecules which act as a control on the cell lineage. The cell densities are described by coupled reaction-diffusion partial differential equations, and the propagating wave front solutions in one dimension are investigated analytically and by numerical solutions. In particular, wavefront propagation speeds are obtained analytically and verified by numerical solutions of the equations. The emphasis is on the effects of the feedback regulations on different stages in the cell lineage. It is found that when the progenitor cell is negatively regulated, the populations of the cell lineage are strongly down-regulated with the steady growth rate of the progenitor cell being driven to zero beyond a critical regulatory strength. An analytic expression for the critical regulation strength in terms of the model parameters is derived and verified by numerical solutions. On the other hand, if the inhibition is acting on the differentiated cells, the change in the population dynamics and wave propagation speed is small. In addition, it is found that only the propagating speed of the progenitor cells is affected by the regulation when the diffusion of the differentiated cells is large. In the presence of de-differentiation, the effect on down-regulating the progenitor population is weakened and there is no effect on the propagation speed due to regulation, suggesting that the effect of regulatory control is diminished by de-differentiation pathways.

The Role of Symmetric Stem Cell Divisions in Tissue Homeostasis

Successful maintenance of cellular lineages critically depends on the fate decision dynamics of stem cells (SCs) upon division. There are three possible strategies with respect to SC fate decision symmetry: (a) asymmetric mode, when each and every SC division produces one SC and one non-SC progeny; (b) symmetric mode, when 50% of all divisions produce two SCs and another 50%-two non-SC progeny; (c) mixed mode, when both the asymmetric and two types of symmetric SC divisions co-exist and are partitioned so that long-term net balance of the lineage output stays constant. Theoretically, either of these strategies can achieve lineage homeostasis. However, it remains unclear which strategy(s) are more advantageous and under what specific circumstances, and what minimal control mechanisms are required to operate them. Here we used stochastic modeling to analyze and quantify the ability of different types of divisions to maintain long-term lineage homeostasis, in the context of different control networks. Using the example of a two-component lineage, consisting of SCs and one type of non-SC progeny, we show that its tight homeostatic control is not necessarily associated with purely asymmetric divisions. Through stochastic analysis and simulations we show that asymmetric divisions can either stabilize or destabilize the lineage system, depending on the underlying control network. We further apply our computational model to biological observations in the context of a two-component lineage of mouse epidermis, where autonomous lineage control has been proposed and notable regional differences, in terms of symmetric division ratio, have been noted-higher in thickened epidermis of the paw skin as compared to ear and tail skin. By using our model we propose a possible explanation for the regional differences in epidermal lineage control strategies. We demonstrate how symmetric divisions can work to stabilize paw epidermis lineage, which experiences high level of micro-injuries and a lack of hair follicles as a back-up source of SCs.

Models, measurement and inference in epithelial tissue dynamics

The majority of solid tumours arise in epithelia and therefore much research effort has gone into investigating the growth, renewal and regulation of these tissues. Here we review different mathematical and computational approaches that have been used to model epithelia. We compare different models and describe future challenges that need to be overcome in order to fully exploit new data which present, for the first time, the real possibility for detailed model validation and comparison.

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.

Characterizing inhibited tumor growth in stem-cell-driven non-spatial cancers

Healthy human tissue is highly regulated to maintain homeostasis. Secreted negative feedback factors that inhibit stem cell division and stem cell self-renewal play a fundamental role in establishing this control. The appearance of abnormal cancerous growth requires an escape from these regulatory mechanisms. In a previous study we found that for non-solid tumors if feedback inhibition on stem cell self-renewal is lost, but the feedback on the division rate is still intact, then the tumor dynamics are characterized by a relatively slow sub-exponential growth that we called inhibited growth. Here we characterize the cell dynamics of inhibited cancer growth by modeling feedback inhibition using Hill equations. We find asymptotic approximations for the growth rates of the stem cell and differentiated cell populations in terms of the strength of the inhibitory signal: stem cells grow as a power law t(1/k+1),and the differentiated cells grow as t(1/k), where k is the Hill coefficient in the feedback law regulating cell divisions. It follows that as the tumor grows, undifferentiated cells take up an increasingly large fraction of the population. Implications of these results for specific cancers including CML are discussed. Understanding how the regulatory mechanisms that continue to operate in cancer affect the rate of disease progression can provide important insights relevant to chronic or other slow progressing types of cancer.

The Interplay between Wnt Mediated Expansion and Negative Regulation of Growth Promotes Robust Intestinal Crypt Structure and Homeostasis

The epithelium of the small intestinal crypt, which has a vital role in protecting the underlying tissue from the harsh intestinal environment, is completely renewed every 4–5 days by a small pool of stem cells at the base of each crypt. How is this renewal controlled and homeostasis maintained, particularly given the rapid nature of this process? Here, based on the recent observations from in vitro “mini gut” studies, we use a hybrid stochastic model of the crypt to investigate how exogenous niche signaling (from Wnt and BMP) combines with auto-regulation to promote homeostasis. This model builds on the sub-cellular element method to account for the three-dimensional structure of the crypt, external regulation by Wnt and BMP, internal regulation by Notch signaling, as well as regulation by internally generated diffusible signals. Results show that Paneth cell derived Wnt signals, which have been observed experimentally to sustain crypts in cultured organs, have a dramatically different influence on niche dynamics than does mesenchyme derived Wnt. While this signaling can indeed act as a redundant backup to the exogenous gradient, it introduces a positive feedback that destabilizes the niche and causes its uncontrolled expansion. We find that in this setting, BMP has a critical role in constraining this expansion, consistent with observations that its removal leads to crypt fission. Further results also point to a new hypothesis for the role of Ephrin mediated motility of Paneth cells, specifically that it is required to constrain niche expansion and maintain the crypt’s spatial structure. Combined, these provide an alternative view of crypt homeostasis where the niche is in a constant state of expansion and the spatial structure of the crypt arises as a balance between this expansion and the action of various sources of negative regulation that hold it in check.

The small intestinal epithelium, like our skin, is constantly being renewed. In the intestine however, this epithelium is exposed to the harsh digestive environment, necessitating much more rapid renewal. Remarkably, the entire epithelium is renewed every 4–5 days. This raises the question, how can the size and structure of this tissue be maintained given this pace. Motivated by recent experimental observations, we construct a three-dimensional, hybrid stochastic model to investigate the mechanisms responsible for homeostasis of this tissue. We find that there are redundant signals created by both the epithelium itself and surrounding tissues that act in parallel to maintain epithelial structure. This redundancy comes at a price however: it introduces the possibility of explosive stem cell population growth. Additional results suggest that other signals along with choreographed motion of cells are responsible for repressing this expansion. Taken together, our results provide a novel hypothesis for how robust but fast renewal of the crypt is achieved: as a balance between expansion, which drives fast renewal and repression, which holds that expansion in check to maintain the crypt’s structure.

Deciphering principles of morphogenesis from temporal and spatial patterns on the integument

BACKGROUND

How tissue patterns form in development and regeneration is a fundamental issue remaining to be fully understood. The integument often forms repetitive units in space (periodic patterning) and time (cyclic renewal), such as feathers and hairs. Integument patterns are visible and experimentally manipulatable, helping us reveal pattern formative processes. Variability is seen in regional phenotypic specificities and temporal cycling at different physiological stages.

RESULTS

Here we show some cellular/molecular bases revealed by analyzing integument patterns. (1) Localized cellular activity (proliferation, rearrangement, apoptosis, differentiation) transforms prototypic organ primordia into specific shapes. Combinatorial positioning of different localized activity zones generates diverse and complex organ forms. (2) Competitive equilibrium between activators and inhibitors regulates stem cells through cyclic quiescence and activation.

CONCLUSIONS

Dynamic interactions between stem cells and their adjacent niche regulate regenerative behavior, modulated by multi-layers of macro-environmental factors (dermis, body hormone status, and external environment). Genomics studies may reveal how positional information of localized cellular activity is stored. In vivo skin imaging and lineage tracing unveils new insights into stem cell plasticity. Principles of self-assembly obtained from the integumentary organ model can be applied to help restore damaged patterns during regenerative wound healing and for tissue engineering to rebuild tissues. Developmental Dynamics 244:905-920, 2015. © 2015 Wiley Periodicals, Inc.

Modeling spatial population dynamics of stem cell lineage in wound healing and cancerogenesis

Modeling the dynamics of cell population in tissues involving stem cell niches allows insight into the control mechanisms of the important wound healing process. It is well known that growth and divisions of stem cells are mainly repressed by niche cells, but can also be activated by signals released from wound. In addition, the proliferation and differentiation among three different types of cell: stem cells (SCs), intermediate progenitor cells (IPCs), and fully differentiated cells (FDCs) in stem cell lineage are under different activation and inhibition controls. We have developed a novel stochastic spatial dynamic model of cells. We can characterize not only overall cell population dynamics, but also details of temporal-spatial relationship of individual cells within a tissue. In our model, the shape, growth, and division of each cell are modeled using a realistic geometric model. Furthermore, the inhibited growth rate, proliferation and differentiation probabilities of individual cells are modeled through feedback loops controlled by secreted factors and wound signals from neighboring cells. With specific proliferation and differentiation probabilities, the actual division type that each cell will take is chosen by a Monte Carlo sampling process. With simulations, we study the effects of different strengths of wound signals to wound healing behaviors. We also study the correlations between chronic wound and cancerogenesis.

Analysis of stochastic stem cell models with control

Understanding the dynamics of stem cell lineages is of central importance both for healthy and cancerous tissues. We study stochastic population dynamics of stem cells and differentiated cells, where cell decisions, such as proliferation vs. differentiation decisions, or division and death decisions, are under regulation from surrounding cells. The goal is to understand how different types of control mechanisms affect the means and variances of cell numbers. We use the assumption of weak dependencies of the regulatory functions (the controls) on the cell populations near the equilibrium to formulate moment equations. We then study three different methods of closure, showing that they all lead to the same results for the highest order terms in the expressions for the moments. We derive simple explicit expressions for the means and the variances of stem cell and differentiated cell numbers. It turns out that the variance is expressed as an algebraic function of partial derivatives of the controls with respect to the population sizes at the equilibrium. We demonstrate that these findings are consistent with the results previously obtained in the context of particular systems, and also present two novel examples with negative and positive control of division and differentiation decisions. This methodology is formulated without any specific assumptions on the functional form of the controls, and thus can be used for any biological system.

A versatile mathematical work-flow to explore how Cancer Stem Cell fate influences tumor progression

Background

Nowadays multidisciplinary approaches combining mathematical models with experimental assays are becoming relevant for the study of biological systems. Indeed, in cancer research multidisciplinary approaches are successfully used to understand the crucial aspects implicated in tumor growth. In particular, the Cancer Stem Cell (CSC) biology represents an area particularly suited to be studied through multidisciplinary approaches, and modeling has significantly contributed to pinpoint the crucial aspects implicated in this theory.

More generally, to acquire new insights on a biological system it is necessary to have an accurate description of the phenomenon, such that making accurate predictions on its future behaviors becomes more likely. In this context, the identification of the parameters influencing model dynamics can be advantageous to increase model accuracy and to provide hints in designing wet experiments. Different techniques, ranging from statistical methods to analytical studies, have been developed. Their applications depend on case-specific aspects, such as the availability and quality of experimental data, and the dimension of the parameter space.

Results

The study of a new model on the CSC-based tumor progression has been the motivation to design a new work-flow that helps to characterize possible system dynamics and to identify those parameters influencing such behaviors. In detail, we extended our recent model on CSC-dynamics creating a new system capable of describing tumor growth during the different stages of cancer progression. Indeed, tumor cells appear to progress through lineage stages like those of normal tissues, being their division auto-regulated by internal feedback mechanisms. These new features have introduced some non-linearities in the model, making it more difficult to be studied by solely analytical techniques. Our new work-flow, based on statistical methods, was used to identify the parameters which influence the tumor growth. The effectiveness of the presented work-flow was firstly verified on two well known models and then applied to investigate our extended CSC model.

Conclusions

We propose a new work-flow to study in a practical and informative way complex systems, allowing an easy identification, interpretation, and visualization of the key model parameters. Our methodology is useful to investigate possible model behaviors and to establish factors driving model dynamics.

Analyzing our new CSC model guided by the proposed work-flow, we found that the deregulation of CSC asymmetric proliferation contributes to cancer initiation, in accordance with several experimental evidences. Specifically, model results indicated that the probability of CSC symmetric proliferation is responsible of a switching-like behavior which discriminates between tumorigenesis and unsustainable tumor growth.

Mechanical model of geometric cell and topological algorithm for cell dynamics from single-cell to formation of monolayered tissues with pattern

Geometric and mechanical properties of individual cells and interactions among neighboring cells are the basis of formation of tissue patterns. Understanding the complex interplay of cells is essential for gaining insight into embryogenesis, tissue development, and other emerging behavior. Here we describe a cell model and an efficient geometric algorithm for studying the dynamic process of tissue formation in 2D (e.g. epithelial tissues). Our approach improves upon previous methods by incorporating properties of individual cells as well as detailed description of the dynamic growth process, with all topological changes accounted for. Cell size, shape, and division plane orientation are modeled realistically. In addition, cell birth, cell growth, cell shrinkage, cell death, cell division, cell collision, and cell rearrangements are now fully accounted for. Different models of cell-cell interactions, such as lateral inhibition during the process of growth, can be studied in detail. Cellular pattern formation for monolayered tissues from arbitrary initial conditions, including that of a single cell, can also be studied in detail. Computational efficiency is achieved through the employment of a special data structure that ensures access to neighboring cells in constant time, without additional space requirement. We have successfully generated tissues consisting of more than 20,000 cells starting from 2 cells within 1 hour. We show that our model can be used to study embryogenesis, tissue fusion, and cell apoptosis. We give detailed study of the classical developmental process of bristle formation on the epidermis of D. melanogaster and the fundamental problem of homeostatic size control in epithelial tissues. Simulation results reveal significant roles of solubility of secreted factors in both the bristle formation and the homeostatic control of tissue size. Our method can be used to study broad problems in monolayered tissue formation. Our software is publicly available.

Cell lineage branching as a strategy for proliferative control

Background

How tissue and organ sizes are specified is one of the great unsolved mysteries in biology. Experiments and mathematical modeling implicate feedback control of cell lineage progression, but a broad understanding of what lineage feedback accomplishes is lacking.

Results

By exploring the possible effects of various biologically relevant disturbances on the dynamic and steady state behaviors of stem cell lineages, we find that the simplest and most frequently studied form of lineage feedback - which we term renewal control - suffers from several serious drawbacks. These reflect fundamental performance limits dictated by universal conservation-type laws, and are independent of parameter choice. Here we show that introducing lineage branches can circumvent all such limitations, permitting effective attenuation of a wide range of perturbations. The type of feedback that achieves such performance - which we term fate control - involves promotion of lineage branching at the expense of both renewal and (primary) differentiation. We discuss the evidence that feedback of just this type occurs in vivo, and plays a role in tissue growth control.

Conclusions

Regulated lineage branching is an effective strategy for dealing with disturbances in stem cell systems. The existence of this strategy provides a dynamics-based justification for feedback control of cell fate in vivo.

See commentary article: http://dx.doi.org/10.1186/s12915-015-0123-7.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0122-8) contains supplementary material, which is available to authorized users.

Engineering control into medicine

The human body is a tightly controlled engineering miracle. However, medical training generally does not cover "control" (in the engineering sense) in physiology, pathophysiology, and therapeutics. A better understanding of how evolved controls maintain normal homeostasis is critical for understanding the failure mode of controlled systems, that is, disease. We believe that teaching and research must incorporate an understanding of the control systems in physiology and take advantage of the quantitative tools used by engineering to understand complex systems. Control systems are ubiquitous in physiology, although often unrecognized. Here we provide selected examples of the role of control in physiology (heart rate variability, immunity), pathophysiology (inflammation in sepsis), and therapeutic devices (diabetes and the artificial pancreas). We also present a high-level background to the concept of robustly controlled systems and examples of clinical insights using the controls framework.

Mathematical Modelling as a Tool to Understand Cell Self-renewal and Differentiation

Mathematical modeling is a powerful technique to address key questions and paradigms in a variety of complex biological systems and can provide quantitative insights into cell kinetics, fate determination and development of cell populations. The chapter is devoted to a review of modeling of the dynamics of stem cell-initiated systems using mathematical methods of ordinary differential equations. Some basic concepts and tools for cell population dynamics are summarized and presented as a gentle introduction to non-mathematicians. The models take into account different plausible mechanisms regulating homeostasis. Two mathematical frameworks are proposed reflecting, respectively, a discrete (punctuated by division events) and a continuous character of transitions between differentiation stages. Advantages and constraints of the mathematical approaches are presented on examples of models of blood systems and compared to patients data on healthy hematopoiesis.

Stochastic control of proliferation and differentiation in stem cell dynamics

In self-renewing tissues, cell lineages consisting of stem cell and classes of daughter cells are the basic units which are responsible for the correct functioning of the organ. Cell proliferation and differentiation in lineages is thought to be mediated by feedback signals. In the simplest case a lineage is comprised of stem cells and differentiated cells. We create a model where stem cell proliferation and differentiation are controlled by the size of cell populations by means of a negative feedback loop. This two-dimensional Markov process allows for an analytical solution for the mean numbers and variances of stem and daughter cells. The mean values and the amounts of variation in cell numbers can be tightly regulated by the parameters of the control loop.

Interactions and tradeoffs between cell recruitment, proliferation, and differentiation affect CNS regeneration

Regeneration of central nervous system (CNS) lesions requires movement of progenitor cells and production of their differentiated progeny. Although damage to the CNS clearly promotes these two processes, the interplay between these complex events and how it affects a response remains elusive. Here, we use spatial stochastic modeling to show that tradeoffs arise between production and recruitment during regeneration. Proper spatial control of cell cycle timing can mitigate these tradeoffs, maximizing recruitment, improving infiltration into the lesion, and reducing wasteful production outside of it. Feedback regulation of cell lineage dynamics alone however leads to spatial defects in cell recruitment, suggesting a novel, to our knowledge, hypothesis for the aggregation of cells to the periphery of a lesion in multiple sclerosis. Interestingly, stronger chemotaxis does not correct this aggregation and instead, substantial random cell motions near the site of the lesion are required to improve CNS regeneration.