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.

A mathematical framework for the emergence of winners and losers in cell competition.

Cell competition is a process in multicellular organisms where cells interact with their neighbours to determine a "winner" or "loser" status. The loser cells are eliminated through programmed cell death, leaving only the winner cells to populate the tissue. Cell competition is context-dependent; the same cell type can win or lose depending on the cell type it is competing against. Hence, winner/loser status is an emergent property. A key question in cell competition is: how do cells acquire their winner/loser status? In this paper, we propose a mathematical framework for studying the emergence of winner/loser status based on a set of quantitative criteria that distinguishes competitive from non-competitive outcomes. We apply this framework in a cell-based modelling context, to both highlight the crucial role of active cell death in cell competition and identify the factors that drive cell competition.

Optimal Control of Collective Electrotaxis in Epithelial Monolayers.

Epithelial monolayers are some of the best-studied models for collective cell migration due to their abundance in multicellular systems and their tractability. Experimentally, the collective migration of epithelial monolayers can be robustly steered e.g. using electric fields, via a process termed electrotaxis. Theoretically, however, the question of how to design an electric field to achieve a desired spatiotemporal movement pattern is underexplored. In this work, we construct and calibrate an ordinary differential equation model to predict the average velocity of the centre of mass of a cellular monolayer in response to stimulation with an electric field. We use this model, in conjunction with optimal control theory, to derive physically realistic optimal electric field designs to achieve a variety of aims, including maximising the total distance travelled by the monolayer, maximising the monolayer velocity, and keeping the monolayer velocity constant during stimulation. Together, this work is the first to present a unified framework for optimal control of collective monolayer electrotaxis and provides a blueprint to optimally steer collective migration using other external cues.

Predicting Radiotherapy Patient Outcomes with Real-Time Clinical Data Using Mathematical Modelling.

Longitudinal tumour volume data from head-and-neck cancer patients show that tumours of comparable pre-treatment size and stage may respond very differently to the same radiotherapy fractionation protocol. Mathematical models are often proposed to predict treatment outcome in this context, and have the potential to guide clinical decision-making and inform personalised fractionation protocols. Hindering effective use of models in this context is the sparsity of clinical measurements juxtaposed with the model complexity required to produce the full range of possible patient responses. In this work, we present a compartment model of tumour volume and tumour composition, which, despite relative simplicity, is capable of producing a wide range of patient responses. We then develop novel statistical methodology and leverage a cohort of existing clinical data to produce a predictive model of both tumour volume progression and the associated level of uncertainty that evolves throughout a patient's course of treatment. To capture inter-patient variability, all model parameters are patient specific, with a bootstrap particle filter-like Bayesian approach developed to model a set of training data as prior knowledge. We validate our approach against a subset of unseen data, and demonstrate both the predictive ability of our trained model and its limitations.

Colec12 and Trail signaling confine cranial neural crest cell trajectories and promote collective cell migration.

BACKGROUND: Collective and discrete neural crest cell (NCC) migratory streams are crucial to vertebrate head patterning. However, the factors that confine NCC trajectories and promote collective cell migration remain unclear. RESULTS: Computational simulations predicted that confinement is required only along the initial one-third of the cranial NCC migratory pathway. This guided our study of Colec12 (Collectin-12, a transmembrane scavenger receptor C-type lectin) and Trail (tumor necrosis factor-related apoptosis-inducing ligand, CD253) which we show expressed in chick cranial NCC-free zones. NCC trajectories are confined by Colec12 or Trail protein stripes in vitro and show significant and distinct changes in cell morphology and dynamic migratory characteristics when cocultured with either protein. Gain- or loss-of-function of either factor or in combination enhanced NCC confinement or diverted cell trajectories as observed in vivo with three-dimensional confocal microscopy, respectively, resulting in disrupted collective migration. CONCLUSIONS: These data provide evidence for Colec12 and Trail as novel NCC microenvironmental factors playing a role to confine cranial NCC trajectories and promote collective cell migration.

Neural crest cells bulldoze through the microenvironment using Aquaporin 1 to stabilize filopodia.

Neural crest migration requires cells to move through an environment filled with dense extracellular matrix and mesoderm to reach targets throughout the vertebrate embryo. Here, we use high-resolution microscopy, computational modeling, and in vitro and in vivo cell invasion assays to investigate the function of Aquaporin 1 (AQP-1) signaling. We find that migrating lead cranial neural crest cells express AQP-1 mRNA and protein, implicating a biological role for water channel protein function during invasion. Differential AQP-1 levels affect neural crest cell speed and direction, as well as the length and stability of cell filopodia. Furthermore, AQP-1 enhances matrix metalloprotease activity and colocalizes with phosphorylated focal adhesion kinases. Colocalization of AQP-1 with EphB guidance receptors in the same migrating neural crest cells has novel implications for the concept of guided bulldozing by lead cells during migration.

Efficient Bayesian inference for mechanistic modelling with high-throughput data.

Bayesian methods are routinely used to combine experimental data with detailed mathematical models to obtain insights into physical phenomena. However, the computational cost of Bayesian computation with detailed models has been a notorious problem. Moreover, while high-throughput data presents opportunities to calibrate sophisticated models, comparing large amounts of data with model simulations quickly becomes computationally prohibitive. Inspired by the method of Stochastic Gradient Descent, we propose a minibatch approach to approximate Bayesian computation. Through a case study of a high-throughput imaging scratch assay experiment, we show that reliable inference can be performed at a fraction of the computational cost of a traditional Bayesian inference scheme. By applying a detailed mathematical model of single cell motility, proliferation and death to a data set of 118 gene knockdowns, we characterise functional subgroups of gene knockdowns, each displaying its own typical combination of local cell density-dependent and -independent motility and proliferation patterns. By comparing these patterns to experimental measurements of cell counts and wound closure, we find that density-dependent interactions play a crucial role in the process of wound healing.

Reliable and efficient parameter estimation using approximate continuum limit descriptions of stochastic models.

Stochastic individual-based mathematical models are attractive for modelling biological phenomena because they naturally capture the stochasticity and variability that is often evident in biological data. Such models also allow us to track the motion of individuals within the population of interest. Unfortunately, capturing this microscopic detail means that simulation and parameter inference can become computationally expensive. One approach for overcoming this computational limitation is to coarse-grain the stochastic model to provide an approximate continuum model that can be solved using far less computational effort. However, coarse-grained continuum models can be biased or inaccurate, particularly for certain parameter regimes. In this work, we combine stochastic and continuum mathematical models in the context of lattice-based models of two-dimensional cell biology experiments by demonstrating how to simulate two commonly used experiments: cell proliferation assays and barrier assays. Our approach involves building a simple statistical model of the discrepancy between the expensive stochastic model and the associated computationally inexpensive coarse-grained continuum model. We form this statistical model based on a limited number of expensive stochastic model evaluations at design points sampled from a user-chosen distribution, corresponding to a computer experiment design problem. With straightforward design point selection schemes, we show that using the statistical model of the discrepancy in tandem with the computationally inexpensive continuum model allows us to carry out prediction and inference while correcting for biases and inaccuracies due to the continuum approximation. We demonstrate this approach by simulating a proliferation assay, where the continuum limit model is the well-known logistic ordinary differential equation, as well as a barrier assay where the continuum limit model is closely related to the well-known Fisher-KPP partial differential equation. We construct an approximate likelihood function for parameter inference, both with and without discrepancy correction terms. Using maximum likelihood estimation, we provide point estimates of the unknown parameters, and use the profile likelihood to characterise the uncertainty in these estimates and form approximate confidence intervals. For the range of inference problems considered, working with the continuum limit model alone leads to biased parameter estimation and confidence intervals with poor coverage. In contrast, incorporating correction terms arising from the statistical model of the model discrepancy allows us to recover the parameters accurately with minimal computational overhead. The main tradeoff is that the associated confidence intervals are typically broader, reflecting the additional uncertainty introduced by the approximation process. All algorithms required to replicate the results in this work are written in the open source Julia language and are available at GitHub.

Parameter identifiability and model selection for sigmoid population growth models.

Sigmoid growth models, such as the logistic, Gompertz and Richards' models, are widely used to study population dynamics ranging from microscopic populations of cancer cells, to continental-scale human populations. Fundamental questions about model selection and parameter estimation are critical if these models are to be used to make practical inferences. However, the question of parameter identifiability - whether a data set contains sufficient information to give unique or sufficiently precise parameter estimates - is often overlooked. We use a profile-likelihood approach to explore practical parameter identifiability using data describing the re-growth of hard coral. With this approach, we explore the relationship between parameter identifiability and model misspecification, finding that the logistic growth model does not suffer identifiability issues for the type of data we consider whereas the Gompertz and Richards' models encounter practical non-identifiability issues. This analysis of parameter identifiability and model selection is important because different growth models are in biological modelling without necessarily considering whether parameters are identifiable. Standard practices that do not consider parameter identifiability can lead to unreliable or imprecise parameter estimates and potentially misleading mechanistic interpretations. For example, using the Gompertz model, the estimate of the time scale of coral re-growth is 625 days when we estimate the initial density from the data, whereas it is 1429 days using a more standard approach where variability in the initial density is ignored. While tools developed here focus on three standard sigmoid growth models only, our theoretical developments are applicable to any sigmoid growth model and any appropriate data set. MATLAB implementations of all software are available on GitHub.

Mathematical modelling of autoimmune myocarditis and the effects of immune checkpoint inhibitors.

Autoimmune myocarditis is a rare, but frequently fatal, side effect of immune checkpoint inhibitors (ICIs), a class of cancer therapies. Despite extensive experimental work on the causes, development and progression of this disease, much still remains unknown about the importance of the different immunological pathways involved. We present a mathematical model of autoimmune myocarditis and the effects of ICIs on its development and progression to either resolution or chronic inflammation. From this, we gain a better understanding of the role of immune cells, cytokines and other components of the immune system in driving the cardiotoxicity of ICIs. We parameterise the model using existing data from the literature, and show that qualitative model behaviour is consistent with disease characteristics seen in patients in an ICI-free context. The bifurcation structures of the model show how the presence of ICIs increases the risk of developing autoimmune myocarditis. This predictive modelling approach is a first step towards determining treatment regimens that balance the benefits of treating cancer with the risk of developing autoimmune myocarditis.

Self-assembly of biological networks via adaptive patterning revealed by avian intradermal muscle network formation.

Networked structures integrate numerous elements into one functional unit, while providing a balance between efficiency, robustness, and flexibility. Understanding how biological networks self-assemble will provide insights into how these features arise. Here, we demonstrate how nature forms exquisite muscle networks that can repair, regenerate, and adapt to external perturbations using the feather muscle network in chicken embryos as a paradigm. The self-assembled muscle networks arise through the implementation of a few simple rules. Muscle fibers extend outward from feather buds in every direction, but only those muscle fibers able to connect to neighboring buds are eventually stabilized. After forming such a nearest-neighbor configuration, the network can be reconfigured, adapting to perturbed bud arrangement or mechanical cues. Our computational model provides a bioinspired algorithm for network self-assembly, with intrinsic or extrinsic cues necessary and sufficient to guide the formation of these regenerative networks. These robust principles may serve as a useful guide for assembling adaptive networks in other contexts.

Quantifying the impact of electric fields on single-cell motility.

Cell motility in response to environmental cues forms the basis of many developmental processes in multicellular organisms. One such environmental cue is an electric field (EF), which induces a form of motility known as electrotaxis. Electrotaxis has evolved in a number of cell types to guide wound healing and has been associated with different cellular processes, suggesting that observed electrotactic behavior is likely a combination of multiple distinct effects arising from the presence of an EF. To determine the different mechanisms by which observed electrotactic behavior emerges, and thus to design EFs that can be applied to direct and control electrotaxis, researchers require accurate quantitative predictions of cellular responses to externally applied fields. Here, we use mathematical modeling to formulate and parameterize a variety of hypothetical descriptions of how cell motility may change in response to an EF. We calibrate our model to observed data using synthetic likelihoods and Bayesian sequential learning techniques and demonstrate that EFs bias cellular motility through only one of a selection of hypothetical mechanisms. We also demonstrate how the model allows us to make predictions about cellular motility under different EFs. The resulting model and calibration methodology will thus form the basis for future data-driven and model-based feedback control strategies based on electric actuation.

Visualizing mesoderm and neural crest cell dynamics during chick head morphogenesis.

Vertebrate head morphogenesis involves carefully-orchestrated tissue growth and cell movements of the mesoderm and neural crest to form the distinct craniofacial pattern. To better understand structural birth defects, it is important that we characterize the dynamics of these processes and learn how they rely on each other. Here we examine this question during chick head morphogenesis using time-lapse imaging, computational modeling, and experiments. We find that head mesodermal cells in culture move in random directions as individuals and move faster in the presence of neural crest cells. In vivo, mesodermal cells migrate in a directed manner and maintain neighbor relationships; neural crest cells travel through the mesoderm at a faster speed. The mesoderm grows with a non-uniform spatio-temporal profile determined by BrdU labeling during the period of faster and more-directed neural crest collective migration through this domain. We use computer simulations to probe the robustness of neural crest stream formation by varying the spatio-temporal growth profile of the mesoderm. We follow this with experimental manipulations that either stop mesoderm growth or prevent neural crest migration and observe changes in the non-manipulated cell population, implying a dynamic feedback between tissue growth and neural crest cell signaling to confer robustness to the system. Overall, we present a novel descriptive analysis of mesoderm and neural crest cell dynamics that reveals the coordination and co-dependence of these two cell populations during head morphogenesis.

The role of mechanical interactions in EMT.

The detachment of cells from the boundary of an epithelial tissue and the subsequent invasion of these cells into surrounding tissues is important for cancer development and wound healing, and is strongly associated with the epithelial-mesenchymal transition (EMT). Chemical signals, such as TGF-ß, produced by surrounding tissue can be uptaken by cells and induce EMT. In this work, we present a novel cell-based discrete mathematical model of mechanical cellular relaxation, cell proliferation, and cell detachment driven by chemically-dependent EMT in an epithelial tissue. A continuum description of the model is then derived in the form of a novel nonlinear free boundary problem. Using the discrete and continuum models we explore how the coupling of chemical transport and mechanical interactions influences EMT, and postulate how this could be used to help control EMT in pathological situations.

Biologically-informed neural networks guide mechanistic modeling from sparse experimental data.

Biologically-informed neural networks (BINNs), an extension of physics-informed neural networks [1], are introduced and used to discover the underlying dynamics of biological systems from sparse experimental data. In the present work, BINNs are trained in a supervised learning framework to approximate in vitro cell biology assay experiments while respecting a generalized form of the governing reaction-diffusion partial differential equation (PDE). By allowing the diffusion and reaction terms to be multilayer perceptrons (MLPs), the nonlinear forms of these terms can be learned while simultaneously converging to the solution of the governing PDE. Further, the trained MLPs are used to guide the selection of biologically interpretable mechanistic forms of the PDE terms which provides new insights into the biological and physical mechanisms that govern the dynamics of the observed system. The method is evaluated on sparse real-world data from wound healing assays with varying initial cell densities [2].

An interdisciplinary approach to investigate collective cell migration in neural crest.

The neural crest serves as a powerful and tractable model paradigm for understanding collective cell migration. The neural crest cell populations are well-known for their long-distance collective migration and contribution to diverse cell lineages during vertebrate development. If neural crest cells fail to reach a target or populate an incorrect location, then improper cell differentiation or uncontrolled cell proliferation can result. A wide range of interdisciplinary studies has been carried out to understand the response of neural crest cells to different stimuli and their ability to migrate to distant targets. In this critical commentary, we illustrate how an interdisciplinary collaboration involving experimental and mathematical modeling has led to a deeper understanding of cranial neural crest cell migration. We identify open questions and propose possible ways to start answering some of the challenges arising.

A practical guide to pseudo-marginal methods for computational inference in systems biology.

For many stochastic models of interest in systems biology, such as those describing biochemical reaction networks, exact quantification of parameter uncertainty through statistical inference is intractable. Likelihood-free computational inference techniques enable parameter inference when the likelihood function for the model is intractable but the generation of many sample paths is feasible through stochastic simulation of the forward problem. The most common likelihood-free method in systems biology is approximate Bayesian computation that accepts parameters that result in low discrepancy between stochastic simulations and measured data. However, it can be difficult to assess how the accuracy of the resulting inferences are affected by the choice of acceptance threshold and discrepancy function. The pseudo-marginal approach is an alternative likelihood-free inference method that utilises a Monte Carlo estimate of the likelihood function. This approach has several advantages, particularly in the context of noisy, partially observed, time-course data typical in biochemical reaction network studies. Specifically, the pseudo-marginal approach facilitates exact inference and uncertainty quantification, and may be efficiently combined with particle filters for low variance, high-accuracy likelihood estimation. In this review, we provide a practical introduction to the pseudo-marginal approach using inference for biochemical reaction networks as a series of case studies. Implementations of key algorithms and examples are provided using the Julia programming language; a high performance, open source programming language for scientific computing (https://github.com/davidwarne/Warne2019_GuideToPseudoMarginal).

Mechanical Cell Competition in Heterogeneous Epithelial Tissues.

Mechanical cell competition is important during tissue development, cancer invasion, and tissue ageing. Heterogeneity plays a key role in practical applications since cancer cells can have different cell stiffness and different proliferation rates than normal cells. To study this phenomenon, we propose a one-dimensional mechanical model of heterogeneous epithelial tissue dynamics that includes cell-length-dependent proliferation and death mechanisms. Proliferation and death are incorporated into the discrete model stochastically and arise as source/sink terms in the corresponding continuum model that we derive. Using the new discrete model and continuum description, we explore several applications including the evolution of homogeneous tissues experiencing proliferation and death, and competition in a heterogeneous setting with a cancerous tissue competing for space with an adjacent normal tissue. This framework allows us to postulate new mechanisms that explain the ability of cancer cells to outcompete healthy cells through mechanical differences rather than an intrinsic proliferative advantage. We advise when the continuum model is beneficial and demonstrate why naively adding source/sink terms to a continuum model without considering the underlying discrete model may lead to incorrect results.

Modelling collective cell migration: neural crest as a model paradigm.

A huge variety of mathematical models have been used to investigate collective cell migration. The aim of this brief review is twofold: to present a number of modelling approaches that incorporate the key factors affecting cell migration, including cell-cell and cell-tissue interactions, as well as domain growth, and to showcase their application to model the migration of neural crest cells. We discuss the complementary strengths of microscale and macroscale models, and identify why it can be important to understand how these modelling approaches are related. We consider neural crest cell migration as a model paradigm to illustrate how the application of different mathematical modelling techniques, combined with experimental results, can provide new biological insights. We conclude by highlighting a number of future challenges for the mathematical modelling of neural crest cell migration.

Self-organization process in newborn skin organoid formation inspires strategy to restore hair regeneration of adult cells.

Organoids made from dissociated progenitor cells undergo tissue-like organization. This in vitro self-organization process is not identical to embryonic organ formation, but it achieves a similar phenotype in vivo. This implies genetic codes do not specify morphology directly; instead, complex tissue architectures may be achieved through several intermediate layers of cross talk between genetic information and biophysical processes. Here we use newborn and adult skin organoids for analyses. Dissociated cells from newborn mouse skin form hair primordia-bearing organoids that grow hairs robustly in vivo after transplantation to nude mice. Detailed time-lapse imaging of 3D cultures revealed unexpected morphological transitions between six distinct phases: dissociated cells, cell aggregates, polarized cysts, cyst coalescence, planar skin, and hair-bearing skin. Transcriptome profiling reveals the sequential expression of adhesion molecules, growth factors, Wnts, and matrix metalloproteinases (MMPs). Functional perturbations at different times discern their roles in regulating the switch from one phase to another. In contrast, adult cells form small aggregates, but then development stalls in vitro. Comparative transcriptome analyses suggest suppressing epidermal differentiation in adult cells is critical. These results inspire a strategy that can restore morphological transitions and rescue the hair-forming ability of adult organoids: (i) continuous PKC inhibition and (ii) timely supply of growth factors (IGF, VEGF), Wnts, and MMPs. This comprehensive study demonstrates that alternating molecular events and physical processes are in action during organoid morphogenesis and that the self-organizing processes can be restored via environmental reprogramming. This tissue-level phase transition could drive self-organization behavior in organoid morphogenies beyond the skin.