Parameter-free methods distinguish Wnt pathway models and guide design of experiments

Expressive rule-based modeling and fast simulation for dynamic compartments

Compartmentalization is vital for cell biological processes. The field of rule-based stochastic simulation has acknowledged this, and many tools and methods have capabilities for compartmentalization. However, mostly, this is limited to a static compartmental hierarchy and does not integrate compartmental changes. Integrating compartmental dynamics is challenging for the design of the modeling language and the simulation engine. The language should support a concise yet flexible modeling of compartmental dynamics. Our work is based on ML-Rules, a rule-based language for multi-level cell biological modeling that supports a wide variety of compartmental dynamics, whose syntax we slightly adapt. To develop an efficient simulation engine for compartmental dynamics, we combine specific data structures and new and existing algorithms and implement them in the Rust programming language. We evaluate the concept and implementation using two case studies from existing cell-biological models. The execution of these models outperforms previous simulations of ML-Rules by two orders of magnitude. Finally, we present a prototype of a WebAssembly-based implementation to allow for a low barrier of entry when exploring the language and associated models without the need for local installation.

Comparative analysis of kinetic realizations of insulin signaling

Several studies have developed dynamical models to understand the underlying mechanisms of insulin signaling, a signaling cascade that leads to the translocation of glucose, the human body's main source of energy. Fortunately, reaction network analysis allows us to extract properties of dynamical systems without depending on their model parameter values. This study focuses on the comparison of insulin signaling in healthy state (INSMS or INSulin Metabolic Signaling) and in type 2 diabetes (INRES or INsulin RESistance) using reaction network analysis. The analysis uses network decomposition to identify the different subsystems involved in insulin signaling (e.g., insulin receptor binding and recycling, GLUT4 translocation, and ERK signaling pathway, among others). Furthermore, results show that INSMS and INRES are similar with respect to some network, structo-kinetic, and kinetic properties. Their differences, however, provide insights into what happens when insulin resistance occurs. First, the variation in the number of species involved in INSMS and INRES suggests that when irregularities occur in the insulin signaling pathway, other complexes (and, hence, other processes) get involved, characterizing insulin resistance. Second, the loss of concordance exhibited by INRES suggests less restrictive interplay between the species involved in insulin signaling, leading to unusual activities in the signaling cascade. Lastly, GLUT4 losing its absolute concentration robustness in INRES may signify that the transporter has lost its reliability in shuttling glucose to the cell, inhibiting efficient cellular energy production. This study also suggests possible applications of the equilibria parametrization and network decomposition, resulting from the analysis, to potentially establish absolute concentration robustness in a species.

Julia for biologists

Major computational challenges exist in relation to the collection, curation, processing and analysis of large genomic and imaging datasets, as well as the simulation of larger and more realistic models in systems biology. Here we discuss how a relative newcomer among programming languages-Julia-is poised to meet the current and emerging demands in the computational biosciences and beyond. Speed, flexibility, a thriving package ecosystem and readability are major factors that make high-performance computing and data analysis available to an unprecedented degree. We highlight how Julia's design is already enabling new ways of analyzing biological data and systems, and we provide a list of resources that can facilitate the transition into Julian computing.

Algebra, Geometry and Topology of ERK Kinetics

The MEK/ERK signalling pathway is involved in cell division, cell specialisation, survival and cell death (Shaul and Seger in Biochim Biophys Acta (BBA)-Mol Cell Res 1773(8):1213-1226, 2007). Here we study a polynomial dynamical system describing the dynamics of MEK/ERK proposed by Yeung et al. (Curr Biol 2019, https://doi.org/10.1016/j.cub.2019.12.052 ) with their experimental setup, data and known biological information. The experimental dataset is a time-course of ERK measurements in different phosphorylation states following activation of either wild-type MEK or MEK mutations associated with cancer or developmental defects. We demonstrate how methods from computational algebraic geometry, differential algebra, Bayesian statistics and computational algebraic topology can inform the model reduction, identification and parameter inference of MEK variants, respectively. Throughout, we show how this algebraic viewpoint offers a rigorous and systematic analysis of such models.

Relating simulation studies by provenance-Developing a family of Wnt signaling models

For many biological systems, a variety of simulation models exist. A new simulation model is rarely developed from scratch, but rather revises and extends an existing one. A key challenge, however, is to decide which model might be an appropriate starting point for a particular problem and why. To answer this question, we need to identify entities and activities that contributed to the development of a simulation model. Therefore, we exploit the provenance data model, PROV-DM, of the World Wide Web Consortium and, building on previous work, continue developing a PROV ontology for simulation studies. Based on a case study of 19 Wnt/β-catenin signaling models, we identify crucial entities and activities as well as useful metadata to both capture the provenance information from individual simulation studies and relate these forming a family of models. The approach is implemented in WebProv, a web application for inserting and querying provenance information. Our specialization of PROV-DM contains the entities Research Question, Assumption, Requirement, Qualitative Model, Simulation Model, Simulation Experiment, Simulation Data, and Wet-lab Data as well as activities referring to building, calibrating, validating, and analyzing a simulation model. We show that most Wnt simulation models are connected to other Wnt models by using (parts of) these models. However, the overlap, especially regarding the Wet-lab Data used for calibration or validation of the models is small. Making these aspects of developing a model explicit and queryable is an important step for assessing and reusing simulation models more effectively. Exposing this information helps to integrate a new simulation model within a family of existing ones and may lead to the development of more robust and valid simulation models. We hope that our approach becomes part of a standardization effort and that modelers adopt the benefits of provenance when considering or creating simulation models.

We revise a provenance ontology for simulation studies of cellular biochemical models. Provenance information is useful for understanding the creation of a simulation model because it not only contains information about the entities and activities that have led to a simulation model but also their relations, all of which can be visualized. It provides additional structure by explicitly recording research questions, assumptions, and requirements and relating them along with data, qualitative models, simulation models, and simulation experiments through a small set of predefined but extensible activities.

We have applied our concept to a family of 19 Wnt signaling models and implemented a web-based tool (WebProv) to store the provenance information from these studies. The resulting provenance graph visualizes the story line of simulation studies and demonstrates the creation and calibration of simulation models, the successive attempts of validation and extension, and shows, beyond an individual simulation study, how the Wnt models are related. Thereby, the steps and sources that contributed to a simulation model are made explicit.

Our approach complements other approaches aimed at facilitating the reuse and assessment of simulation products in systems biology such as model repositories as well as annotation and documentation guidelines.

Joining and decomposing reaction networks

In systems and synthetic biology, much research has focused on the behavior and design of single pathways, while, more recently, experimental efforts have focused on how cross-talk (coupling two or more pathways) or inhibiting molecular function (isolating one part of the pathway) affects systems-level behavior. However, the theory for tackling these larger systems in general has lagged behind. Here, we analyze how joining networks (e.g., cross-talk) or decomposing networks (e.g., inhibition or knock-outs) affects three properties that reaction networks may possess-identifiability (recoverability of parameter values from data), steady-state invariants (relationships among species concentrations at steady state, used in model selection), and multistationarity (capacity for multiple steady states, which correspond to multiple cell decisions). Specifically, we prove results that clarify, for a network obtained by joining two smaller networks, how properties of the smaller networks can be inferred from or can imply similar properties of the original network. Our proofs use techniques from computational algebraic geometry, including elimination theory and differential algebra.

Modelling the effect of subcellular mutations on the migration of cells in the colorectal crypt

Background

Many cancers arise from mutations in cells within epithelial tissues. Mutations manifesting at the subcellular level influence the structure and function of the tissue resulting in cancer. Previous work has proposed how cell level properties can lead to mutant cell invasion, but has not incorporated detailed subcellular modelling

Results

We present a framework that allows the straightforward integration and simulation of SBML representations of subcellular dynamics within multiscale models of epithelial tissues. This allows us to investigate the effect of mutations in subcellular pathways on the migration of cells within the colorectal crypt. Using multiple models we find that mutations in APC, a key component in the Wnt signalling pathway, can bias neutral drift and can also cause downward invasion of mutant cells in the crypt.

Conclusions

Our framework allows us to investigate how subcellular mutations, i.e. knockouts and knockdowns, affect cell-level properties and the resultant migration of cells within epithelial tissues. In the context of the colorectal crypt, we see that mutations in APC can lead directly to mutant cell invasion.

Cross-talk between Hippo and Wnt signalling pathways in intestinal crypts: Insights from an agent-based model

Intestinal crypts are responsible for the total cell renewal of the lining of the intestines; this turnover is governed by the interplay between signalling pathways and the cell cycle. The role of Wnt signalling in cell proliferation and differentiation in the intestinal crypt has been extensively studied, with increased signalling found towards the lower regions of the crypt. Recent studies have shown that the Wnt signalling gradient found within the crypt may arise as a result of division-based spreading from a Wnt ‘reservoir’ at the crypt base. The discovery of the Hippo pathway’s involvement in maintaining crypt homeostasis is more recent; a mechanistic understanding of Hippo pathway dynamics, and its possible cross-talk with the Wnt pathway, remains lacking. To explore how the interplay between these pathways may control crypt homeostasis, we extended an ordinary differential equation model of the Wnt signalling pathway to include a phenomenological description of Hippo signalling in single cells, and then coupled it to a cell-based description of cell movement, proliferation and contact inhibition in agent-based simulations. Furthermore, we compared an imposed Wnt gradient with a division-based Wnt gradient model. Our results suggest that Hippo signalling affects the Wnt pathway by reducing the presence of free cytoplasmic β-catenin, causing cell cycle arrest. We also show that a division-based spreading of Wnt can form a Wnt gradient, resulting in proliferative dynamics comparable to imposed-gradient models. Finally, a simulated APC double mutant, with misregulated Wnt and Hippo signalling activity, is predicted to cause monoclonal conversion of the crypt.

Dynamically characterizing individual clinical change by the steady state of disease-associated pathway

Background

Along with the development of precision medicine, individual heterogeneity is attracting more and more attentions in clinical research and application. Although the biomolecular reaction seems to be some various when different individuals suffer a same disease (e.g. virus infection), the final pathogen outcomes of individuals always can be mainly described by two categories in clinics, i.e. symptomatic and asymptomatic. Thus, it is still a great challenge to characterize the individual specific intrinsic regulatory convergence during dynamic gene regulation and expression. Except for individual heterogeneity, the sampling time also increase the expression diversity, so that, the capture of similar steady biological state is a key to characterize individual dynamic biological processes.

Results

Assuming the similar biological functions (e.g. pathways) should be suitable to detect consistent functions rather than chaotic genes, we design and implement a new computational framework (ABP: Attractor analysis of Boolean network of Pathway). ABP aims to identify the dynamic phenotype associated pathways in a state-transition manner, using the network attractor to model and quantify the steady pathway states characterizing the final steady biological sate of individuals (e.g. normal or disease). By analyzing multiple temporal gene expression datasets of virus infections, ABP has shown its effectiveness on identifying key pathways associated with phenotype change; inferring the consensus functional cascade among key pathways; and grouping pathway activity states corresponding to disease states.

Conclusions

Collectively, ABP can detect key pathways and infer their consensus functional cascade during dynamical process (e.g. virus infection), and can also categorize individuals with disease state well, which is helpful for disease classification and prediction.

Src- and confinement-dependent FAK activation causes E-cadherin relaxation and β-catenin activity

β-Catenin is a transcription cofactor proposed to be released from E-cadherin upon mechanically induced phosphorylation. However, evidence for this mechanism is lacking. Gayrard et al. show instead that during epithelial-to-mesenchymal transition, Src- and multicellular confinement–dependent FAK-induced cytoskeleton remodeling causes E-cadherin tension relaxation and phosphorylation-independent β-catenin nuclear translocation from the membrane.

In epithelia, E-cadherin cytoplasmic tail is under cytoskeleton-generated tension via a link that contains β-catenin. A cotranscription factor, β-catenin, is also active in morphogenetic processes associated with epithelial-to-mesenchymal transition. β-Catenin signaling appears mechanically inducible and was proposed to follow phosphorylation-induced β-catenin release from E-cadherin. Evidence for this mechanism is lacking, and whether E-cadherin tension is involved is unknown. To test this, we combined quantitative fluorescence microscopies with genetic and pharmacological perturbations of epithelial-to-mesenchymal transition–induced cells in culture. We showed that β-catenin nuclear activity follows a substantial release from the membrane specific to migrating cells and requires multicellular deconfinement and Src activity. Selective nuclear translocation occurs downstream of focal adhesion kinase activation, which targets E-cadherin tension relaxation through actomyosin remodeling. In contrast, phosphorylations of the cadherin/catenin complex are not substantially required. These data demonstrate that E-cadherin acts as a sensor of intracellular mechanics in a crosstalk with cell-substrate adhesions that target β-catenin signaling.

Hilbert-Schmidt and Sobol sensitivity indices for static and time series Wnt signaling measurements in colorectal cancer - part A

Background

Ever since the accidental discovery of Wingless [Sharma R.P., Drosophila information service, 1973, 50, p 134], research in the field of Wnt signaling pathway has taken significant strides in wet lab experiments and various cancer clinical trials, augmented by recent developments in advanced computational modeling of the pathway. Information rich gene expression profiles reveal various aspects of the signaling pathway and help in studying different issues simultaneously. Hitherto, not many computational studies exist which incorporate the simultaneous study of these issues.

Results

This manuscript ∙ explores the strength of contributing factors in the signaling pathway, ∙ analyzes the existing causal relations among the inter/extracellular factors effecting the pathway based on prior biological knowledge and ∙ investigates the deviations in fold changes in the recently found prevalence of psychophysical laws working in the pathway. To achieve this goal, local and global sensitivity analysis is conducted on the (non)linear responses between the factors obtained from static and time series expression profiles using the density (Hilbert-Schmidt Information Criterion) and variance (Sobol) based sensitivity indices.

Conclusion

The results show the advantage of using density based indices over variance based indices mainly due to the former’s employment of distance measures & the kernel trick via Reproducing kernel Hilbert space (RKHS) that capture nonlinear relations among various intra/extracellular factors of the pathway in a higher dimensional space. In time series data, using these indices it is now possible to observe where in time, which factors get influenced & contribute to the pathway, as changes in concentration of the other factors are made. This synergy of prior biological knowledge, sensitivity analysis & representations in higher dimensional spaces can facilitate in time based administration of target therapeutic drugs & reveal hidden biological information within colorectal cancer samples.

Numerical algebraic geometry for model selection and its application to the life sciences

Researchers working with mathematical models are often confronted by the related problems of parameter estimation, model validation and model selection. These are all optimization problems, well known to be challenging due to nonlinearity, non-convexity and multiple local optima. Furthermore, the challenges are compounded when only partial data are available. Here, we consider polynomial models (e.g. mass-action chemical reaction networks at steady state) and describe a framework for their analysis based on optimization using numerical algebraic geometry. Specifically, we use probability-one polynomial homotopy continuation methods to compute all critical points of the objective function, then filter to recover the global optima. Our approach exploits the geometrical structures relating models and data, and we demonstrate its utility on examples from cell signalling, synthetic biology and epidemiology.

Paracrine and autocrine regulation of gene expression by Wnt-inhibitor Dickkopf in wild-type and mutant hepatocytes

Background

Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/ β-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/ β-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells.

Results

Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/ β-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength.

Conclusions

Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter- and intracellular processes, such as Dkk feedback, diffusion and Wnt/ β-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0470-9) contains supplementary material, which is available to authorized users.

A Theoretical Model of the Wnt Signaling Pathway in the Epithelial Mesenchymal Transition

BACKGROUND

Following the formation of a primary carcinoma, neoplastic cells metastasize by undergoing the epithelial mesenchymal transition (EMT), which is triggered by cues from inflammatory and stromal cells in the microenvironment. EMT allows epithelial cells to lose their highly adhesive nature and instead adopt the spindle-like appearance, as well as the invasive and migratory behavior, of mesenchymal cells. We hypothesize that a bistable switch between the epithelial and mesenchymal phenotypes governs EMT, allowing the cell to maintain its mesenchymal phenotype even after it leaves the primary tumor microenvironment and EMT-inducing extracellular signal.

RESULTS

This work presents a simple mathematical model of EMT, specifically the roles played by four key proteins in the Wnt signaling pathway: Dishevelled (Dvl), E-cadherin, β-catenin, and Slug. The model predicts that following activation of the Wnt pathway, an epithelial cell in the primary carcinoma must attain a threshold level of membrane-bound Dvl to convert to the mesenchymal-like phenotype and maintain that phenotype once it has migrated away from the primary tumor. Furthermore, sensitivity analysis of the model suggests that in both the epithelial and the mesenchymal states, the steady state behavior of E-cadherin and the transcription factor Slug are sensitive to changes in the degradation rate of Slug, while E-cadherin is also sensitive to the IC₅₀ (half-maximal) concentration of Slug necessary to inhibit E-cadherin production. The steady state behavior of Slug exhibits sensitivity to changes in the rate at which it is induced by β-catenin upon activation of the Wnt pathway. In the presence of sufficient amount of Wnt ligand, E-cadherin levels are sensitive to the ratio of the rate of Slug activation via β-catenin to the IC₅₀ concentration of Slug necessary to inhibit E-cadherin production.

CONCLUSIONS

The sensitivity of E-cadherin to the degradation rate of Slug, as well as the IC₅₀ concentration of Slug necessary to inhibit E-cadherin production, shows how the adhesive nature of the cell depends on finely-tuned regulation of Slug. By highlighting the role of β-catenin in the activation of EMT and the relationship between E-cadherin and Slug, this model identifies critical parameters of therapeutic concern, such as the threshold level of Dvl necessary to inactivate the GSK-3β complex mediating β-catenin degradation, the rate at which β-catenin translocates to the nucleus, and the IC₅₀ concentration of Slug needed to inhibit E-cadherin production.

Chemical Reaction Network Theory elucidates sources of multistability in interferon signaling

Bistability has important implications in signaling pathways, since it indicates a potential cell decision between alternative outcomes. We present two approaches developed in the framework of the Chemical Reaction Network Theory for easy and efficient search of multiple steady state behavior in signaling networks (both with and without mass conservation), and apply them to search for sources of bistability at different levels of the interferon signaling pathway. Different type I interferon subtypes and/or doses are known to elicit differential bioactivities (ranging from antiviral, antiproliferative to immunomodulatory activities). How different signaling outcomes can be generated through the same receptor and activating the same JAK/STAT pathway is still an open question. Here, we detect bistability at the level of early STAT signaling, showing how two different cell outcomes are achieved under or above a threshold in ligand dose or ligand-receptor affinity. This finding could contribute to explain the differential signaling (antiviral vs apoptotic) depending on interferon dose and subtype (α vs β) observed in type I interferons.

The role of the Hes1 crosstalk hub in Notch-Wnt interactions of the intestinal crypt

The Notch pathway plays a vital role in determining whether cells in the intestinal epithelium adopt a secretory or an absorptive phenotype. Cell fate specification is coordinated via Notch's interaction with the canonical Wnt pathway. Here, we propose a new mathematical model of the Notch and Wnt pathways, in which the Hes1 promoter acts as a hub for pathway crosstalk. Computational simulations of the model can assist in understanding how healthy intestinal tissue is maintained, and predict the likely consequences of biochemical knockouts upon cell fate selection processes. Chemical reaction network theory (CRNT) is a powerful, generalised framework which assesses the capacity of our model for monostability or multistability, by analysing properties of the underlying network structure without recourse to specific parameter values or functional forms for reaction rates. CRNT highlights the role of β-catenin in stabilising the Notch pathway and damping oscillations, demonstrating that Wnt-mediated actions on the Hes1 promoter can induce dynamic transitions in the Notch system, from multistability to monostability. Time-dependent model simulations of cell pairs reveal the stabilising influence of Wnt upon the Notch pathway, in which β-catenin- and Dsh-mediated action on the Hes1 promoter are key in shaping the subcellular dynamics. Where Notch-mediated transcription of Hes1 dominates, there is Notch oscillation and maintenance of fate flexibility; Wnt-mediated transcription of Hes1 favours bistability akin to cell fate selection. Cells could therefore regulate the proportion of Wnt- and Notch-mediated control of the Hes1 promoter to coordinate the timing of cell fate selection as they migrate through the intestinal epithelium and are subject to reduced Wnt stimuli. Furthermore, mutant cells characterised by hyperstimulation of the Wnt pathway may, through coupling with Notch, invert cell fate in neighbouring healthy cells, enabling an aberrant cell to maintain its neighbours in mitotically active states.

Intermediates, catalysts, persistence, and boundary steady states

For dynamical systems arising from chemical reaction networks, persistence is the property that each species concentration remains positively bounded away from zero, as long as species concentrations were all positive in the beginning. We describe two graphical procedures for simplifying reaction networks without breaking known necessary or sufficient conditions for persistence, by iteratively removing so-called intermediates and catalysts from the network. The procedures are easy to apply and, in many cases, lead to highly simplified network structures, such as monomolecular networks. For specific classes of reaction networks, we show that these conditions for persistence are equivalent to one another. Furthermore, they can also be characterized by easily checkable strong connectivity properties of a related graph. In particular, this is the case for (conservative) monomolecular networks, as well as cascades of a large class of post-translational modification systems (of which the MAPK cascade and the n-site futile cycle are prominent examples). Since one of the aforementioned sufficient conditions for persistence precludes the existence of boundary steady states, our method also provides a graphical tool to check for that.

Mathematical and Statistical Techniques for Systems Medicine: The Wnt Signaling Pathway as a Case Study

The last decade has seen an explosion in models that describe phenomena in systems medicine. Such models are especially useful for studying signaling pathways, such as the Wnt pathway. In this chapter we use the Wnt pathway to showcase current mathematical and statistical techniques that enable modelers to gain insight into (models of) gene regulation and generate testable predictions. We introduce a range of modeling frameworks, but focus on ordinary differential equation (ODE) models since they remain the most widely used approach in systems biology and medicine and continue to offer great potential. We present methods for the analysis of a single model, comprising applications of standard dynamical systems approaches such as nondimensionalization, steady state, asymptotic and sensitivity analysis, and more recent statistical and algebraic approaches to compare models with data. We present parameter estimation and model comparison techniques, focusing on Bayesian analysis and coplanarity via algebraic geometry. Our intention is that this (non-exhaustive) review may serve as a useful starting point for the analysis of models in systems medicine.

Algebraic Systems Biology: A Case Study for the Wnt Pathway

Steady-state analysis of dynamical systems for biological networks gives rise to algebraic varieties in high-dimensional spaces whose study is of interest in their own right. We demonstrate this for the shuttle model of the Wnt signaling pathway. Here, the variety is described by a polynomial system in 19 unknowns and 36 parameters. It has degree 9 over the parameter space. This case study explores multistationarity, model comparison, dynamics within regions of the state space, identifiability, and parameter estimation, from a geometric point of view. We employ current methods from computational algebraic geometry, polyhedral geometry, and combinatorics.

Cellular population dynamics control the robustness of the stem cell niche

Within populations of cells, fate decisions are controlled by an indeterminate combination of cell-intrinsic and cell-extrinsic factors. In the case of stem cells, the stem cell niche is believed to maintain ‘stemness’ through communication and interactions between the stem cells and one or more other cell-types that contribute to the niche conditions. To investigate the robustness of cell fate decisions in the stem cell hierarchy and the role that the niche plays, we introduce simple mathematical models of stem and progenitor cells, their progeny and their interplay in the niche. These models capture the fundamental processes of proliferation and differentiation and allow us to consider alternative possibilities regarding how niche-mediated signalling feedback regulates the niche dynamics. Generalised stability analysis of these stem cell niche systems enables us to describe the stability properties of each model. We find that although the number of feasible states depends on the model, their probabilities of stability in general do not: stem cell–niche models are stable across a wide range of parameters. We demonstrate that niche-mediated feedback increases the number of stable steady states, and show how distinct cell states have distinct branching characteristics. The ecological feedback and interactions mediated by the stem cell niche thus lend (surprisingly) high levels of robustness to the stem and progenitor cell population dynamics. Furthermore, cell–cell interactions are sufficient for populations of stem cells and their progeny to achieve stability and maintain homeostasis. We show that the robustness of the niche – and hence of the stem cell pool in the niche – depends only weakly, if at all, on the complexity of the niche make-up: simple as well as complicated niche systems are capable of supporting robust and stable stem cell dynamics.

Summary: Stem cell niche dynamics are very robust to external and physiological perturbations because proliferation and differentiation are naturally balanced and controlled by the reliance on a shared niche environment.