The invasion speed of cell migration models with realistic cell cycle time distributions

Random walk models in the life sciences: including births, deaths and local interactions

Random walks and related spatial stochastic models have been used in a range of application areas, including animal and plant ecology, infectious disease epidemiology, developmental biology, wound healing and oncology. Classical random walk models assume that all individuals in a population behave independently, ignoring local physical and biological interactions. This assumption simplifies the mathematical description of the population considerably, enabling continuum-limit descriptions to be derived and used in model analysis and fitting. However, interactions between individuals can have a crucial impact on population-level behaviour. In recent decades, research has increasingly been directed towards models that include interactions, including physical crowding effects and local biological processes such as adhesion, competition, dispersal, predation and adaptive directional bias. In this article, we review the progress that has been made with models of interacting individuals. We aim to provide an overview that is accessible to researchers in application areas, as well as to specialist modellers. We focus particularly on derivation of asymptotically exact or approximate continuum-limit descriptions and simplified deterministic models of mean-field behaviour and resulting spatial patterns. We provide worked examples and illustrative results of selected models. We conclude with a discussion of current areas of focus and future challenges.

Quantifying cell cycle regulation by tissue crowding

The spatiotemporal coordination and regulation of cell proliferation is fundamental in many aspects of development and tissue maintenance. Cells have the ability to adapt their division rates in response to mechanical constraints, yet we do not fully understand how cell proliferation regulation impacts cell migration phenomena. Here, we present a minimal continuum model of cell migration with cell cycle dynamics, which includes density-dependent effects and hence can account for cell proliferation regulation. By combining minimal mathematical modeling, Bayesian inference, and recent experimental data, we quantify the impact of tissue crowding across different cell cycle stages in epithelial tissue expansion experiments. Our model suggests that cells sense local density and adapt cell cycle progression in response, during G1 and the combined S/G2/M phases, providing an explicit relationship between each cell-cycle-stage duration and local tissue density, which is consistent with several experimental observations. Finally, we compare our mathematical model's predictions to different experiments studying cell cycle regulation and present a quantitative analysis on the impact of density-dependent regulation on cell migration patterns. Our work presents a systematic approach for investigating and analyzing cell cycle data, providing mechanistic insights into how individual cells regulate proliferation, based on population-based experimental measurements.

Solving the stochastic dynamics of population growth

Population growth is a fundamental process in ecology and evolution. The population size dynamics during growth are often described by deterministic equations derived from kinetic models. Here, we simulate several population growth models and compare the size averaged over many stochastic realizations with the deterministic predictions. We show that these deterministic equations are generically bad predictors of the average stochastic population dynamics. Specifically, deterministic predictions overestimate the simulated population sizes, especially those of populations starting with a small number of individuals. Describing population growth as a stochastic birth process, we prove that the discrepancy between deterministic predictions and simulated data is due to unclosed-moment dynamics. In other words, the deterministic approach does not consider the variability of birth times, which is particularly important with small population sizes. We show that some moment-closure approximations describe the growth dynamics better than the deterministic prediction. However, they do not reduce the error satisfactorily and only apply to some population growth models. We explicitly solve the stochastic growth dynamics, and our solution applies to any population growth model. We show that our solution exactly quantifies the dynamics of a community composed of different strains and correctly predicts the fixation probability of a strain in a serial dilution experiment. Our work sets the foundations for a more faithful modeling of community and population dynamics. It will allow the development of new tools for a more accurate analysis of experimental and empirical results, including the inference of important growth parameters.

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.

Equivalence framework for an age-structured multistage representation of the cell cycle

We develop theoretical equivalences between stochastic and deterministic models for populations of individual cells stratified by age. Specifically, we develop a hierarchical system of equations describing the full dynamics of an age-structured multistage Markov process for approximating cell cycle time distributions. We further demonstrate that the resulting mean behavior is equivalent, over large timescales, to the classical McKendrick-von Foerster integropartial differential equation. We conclude by extending this framework to a spatial context, facilitating the modeling of traveling wave phenomena and cell-mediated pattern formation. More generally, this methodology may be extended to myriad reaction-diffusion processes for which the age of individuals is relevant to the dynamics.

A novel mathematical model of heterogeneous cell proliferation

We present a novel mathematical model of heterogeneous cell proliferation where the total population consists of a subpopulation of slow-proliferating cells and a subpopulation of fast-proliferating cells. The model incorporates two cellular processes, asymmetric cell division and induced switching between proliferative states, which are important determinants for the heterogeneity of a cell population. As motivation for our model we provide experimental data that illustrate the induced-switching process. Our model consists of a system of two coupled delay differential equations with distributed time delays and the cell densities as functions of time. The distributed delays are bounded and allow for the choice of delay kernel. We analyse the model and prove the nonnegativity and boundedness of solutions, the existence and uniqueness of solutions, and the local stability characteristics of the equilibrium points. We find that the parameters for induced switching are bifurcation parameters and therefore determine the long-term behaviour of the model. Numerical simulations illustrate and support the theoretical findings, and demonstrate the primary importance of transient dynamics for understanding the evolution of many experimental cell populations.

Synchronized oscillations in growing cell populations are explained by demographic noise

Understanding synchrony in growing populations is important for applications as diverse as epidemiology and cancer treatment. Recent experiments employing fluorescent reporters in melanoma cell lines have uncovered growing subpopulations exhibiting sustained oscillations, with nearby cells appearing to synchronize their cycles. In this study, we demonstrate that the behavior observed is consistent with long-lasting transient phenomenon initiated and amplified by the finite-sample effects and demographic noise. We present a novel mathematical analysis of a multistage model of cell growth, which accurately reproduces the synchronized oscillations. As part of the analysis, we elucidate the transient and asymptotic phases of the dynamics and derive an analytical formula to quantify the effect of demographic noise in the appearance of the oscillations. The implications of these findings are broad, such as providing insight into experimental protocols that are used to study the growth of asynchronous populations and, in particular, those investigations relating to anticancer drug discovery.

Coupled differentiation and division of embryonic stem cells inferred from clonal snapshots

The deluge of single-cell data obtained by sequencing, imaging and epigenetic markers has led to an increasingly detailed description of cell state. However, it remains challenging to identify how cells transition between different states, in part because data are typically limited to snapshots in time. A prerequisite for inferring cell state transitions from such snapshots is to distinguish whether transitions are coupled to cell divisions. To address this, we present two minimal branching process models of cell division and differentiation in a well-mixed population. These models describe dynamics where differentiation and division are coupled or uncoupled. For each model, we derive analytic expressions for each subpopulation's mean and variance and for the likelihood, allowing exact Bayesian parameter inference and model selection in the idealised case of fully observed trajectories of differentiation and division events. In the case of snapshots, we present a sample path algorithm and use this to predict optimal temporal spacing of measurements for experimental design. We then apply this methodology to an in vitro dataset assaying the clonal growth of epiblast stem cells in culture conditions promoting self-renewal or differentiation. Here, the larger number of cell states necessitates approximate Bayesian computation. For both culture conditions, our inference supports the model where cell state transitions are coupled to division. For culture conditions promoting differentiation, our analysis indicates a possible shift in dynamics, with these processes becoming more coupled over time.

NetLogo agent-based models as tools for understanding the self-organization of cell fate, morphogenesis and collective migration of the zebrafish posterior Lateral Line primordium

Interactions between primordium cells and their environment determines the self-organization of the zebrafish posterior Lateral Line primordium as it migrates under the skin from the ear to the tip of the tail forming and depositing neuromasts to spearhead formation of the posterior Lateral Line sensory system. In this review we describe how the NetLogo agent-based programming environment has been used in our lab to visualize and explore how self-generated chemokine gradients determine collective migration, how the dynamics of Wnt signaling can be used to predict patterns of neuromast deposition, and how previously defined interactions between Wnt and Fgf signaling systems have the potential to determine the periodic formation of center-biased Fgf signaling centers in the wake of a shrinking Wnt system. We also describe how NetLogo was used as a database for storing and visualizing the results of in toto lineage analysis of all cells in the migrating primordium. Together, the models illustrate how this programming environment can be used in diverse ways to integrate what has been learnt from biological experiments about the nature of interactions between cells and their environment, and explore how these interactions could potentially determine emergent patterns of cell fate specification, morphogenesis and collective migration of the zebrafish posterior Lateral Line primordium.