Parsing patterns: Emerging roles of tissue self-organization in health and disease.

Patterned morphologies, such as segments, spirals, stripes, and spots, frequently emerge during embryogenesis through self-organized coordination between cells. Yet, complex patterns also emerge in adults, suggesting that the capacity for spontaneous self-organization is a ubiquitous property of biological tissues. We review current knowledge on the principles and mechanisms of self-organized patterning in embryonic tissues and explore how these principles and mechanisms apply to adult tissues that exhibit features of patterning. We discuss how and why spontaneous pattern generation is integral to homeostasis and healing of tissues, illustrating it with examples from regenerative biology. We examine how aberrant self-organization underlies diverse pathological states, including inflammatory skin disorders and tumors. Lastly, we posit that based on such blueprints, targeted engineering of pattern-driving molecular circuits can be leveraged for synthetic biology and the generation of organoids with intricate patterns.

Gap junctions in Turing-type periodic feather pattern formation.

Periodic patterning requires coordinated cell-cell interactions at the tissue level. Turing showed, using mathematical modeling, how spatial patterns could arise from the reactions of a diffusive activator-inhibitor pair in an initially homogeneous 2D field. Most activators and inhibitors studied in biological systems are proteins, and the roles of cell-cell interaction, ions, bioelectricity, etc. are only now being identified. Gap junctions (GJs) mediate direct exchanges of ions or small molecules between cells, enabling rapid long-distance communications in a cell collective. They are therefore good candidates for propagating nonprotein-based patterning signals that may act according to the Turing principles. Here, we explore the possible roles of GJs in Turing-type patterning using feather pattern formation as a model. We found 7 of the 12 investigated GJ isoforms are highly dynamically expressed in the developing chicken skin. In ovo functional perturbations of the GJ isoform, connexin 30, by siRNA and the dominant-negative mutant applied before placode development led to disrupted primary feather bud formation. Interestingly, inhibition of gap junctional intercellular communication (GJIC) in the ex vivo skin explant culture allowed the sequential emergence of new feather buds at specific spatial locations relative to the existing primary buds. The results suggest that GJIC may facilitate the propagation of long-distance inhibitory signals. Thus, inhibition of GJs may stimulate Turing-type periodic feather pattern formation during chick skin development, and the removal of GJ activity would enable the emergence of new feather buds if the local environment were competent and the threshold to form buds was reached. We further propose Turing-based computational simulations that can predict the sequential appearance of these ectopic buds. Our models demonstrate how a Turing activator-inhibitor system can continue to generate patterns in the competent morphogenetic field when the level of intercellular communication at the tissue scale is modulated.

Enhanced perfusion following exposure to radiotherapy: A theoretical investigation.

Tumour angiogenesis leads to the formation of blood vessels that are structurally and spatially heterogeneous. Poor blood perfusion, in conjunction with increased hypoxia and oxygen heterogeneity, impairs a tumour's response to radiotherapy. The optimal strategy for enhancing tumour perfusion remains unclear, preventing its regular deployment in combination therapies. In this work, we first identify vascular architectural features that correlate with enhanced perfusion following radiotherapy, using in vivo imaging data from vascular tumours. Then, we present a novel computational model to determine the relationship between these architectural features and blood perfusion in silico. If perfusion is defined to be the proportion of vessels that support blood flow, we find that vascular networks with small mean diameters and large numbers of angiogenic sprouts show the largest increases in perfusion post-irradiation for both biological and synthetic tumours. We also identify cases where perfusion increases due to the pruning of hypoperfused vessels, rather than blood being rerouted. These results indicate the importance of considering network composition when determining the optimal irradiation strategy. In the future, we aim to use our findings to identify tumours that are good candidates for perfusion enhancement and to improve the efficacy of combination therapies.

Using a probabilistic approach to derive a two-phase model of flow-induced cell migration.

Interstitial fluid flow is a feature of many solid tumors. In vitro experiments have shown that such fluid flow can direct tumor cell movement upstream or downstream depending on the balance between the competing mechanisms of tensotaxis (cell migration up stress gradients) and autologous chemotaxis (downstream cell movement in response to flow-induced gradients of self-secreted chemoattractants). In this work we develop a probabilistic-continuum, two-phase model for cell migration in response to interstitial flow. We use a kinetic description for the cell velocity probability density function, and model the flow-dependent mechanical and chemical stimuli as forcing terms that bias cell migration upstream and downstream. Using velocity-space averaging, we reformulate the model as a system of continuum equations for the spatiotemporal evolution of the cell volume fraction and flux in response to forcing terms that depend on the local direction and magnitude of the mechanochemical cues. We specialize our model to describe a one-dimensional cell layer subject to fluid flow. Using a combination of numerical simulations and asymptotic analysis, we delineate the parameter regime where transitions from downstream to upstream cell migration occur. As has been observed experimentally, the model predicts downstream-oriented chemotactic migration at low cell volume fractions, and upstream-oriented tensotactic migration at larger volume fractions. We show that the locus of the critical volume fraction, at which the system transitions from downstream to upstream migration, is dominated by the ratio of the rate of chemokine secretion and advection. Our model also predicts that, because the tensotactic stimulus depends strongly on the cell volume fraction, upstream, tensotaxis-dominated migration occurs only transiently when the cells are initially seeded, and transitions to downstream, chemotaxis-dominated migration occur at later times due to the dispersive effect of cell diffusion.

Turing Pattern Formation in Reaction-Cross-Diffusion Systems with a Bilayer Geometry.

Conditions for self-organisation via Turing's mechanism in biological systems represented by reaction-diffusion or reaction-cross-diffusion models have been extensively studied. Nonetheless, the impact of tissue stratification in such systems is under-explored, despite its ubiquity in the context of a thin epithelium overlying connective tissue, for instance the epidermis and underlying dermal mesenchyme of embryonic skin. In particular, each layer can be subject to extensively different biochemical reactions and transport processes, with chemotaxis - a special case of cross-diffusion - often present in the mesenchyme, contrasting the solely molecular transport typically found in the epidermal layer. We study Turing patterning conditions for a class of reaction-cross-diffusion systems in bilayered regions, with a thin upper layer and coupled by a linear transport law. In particular, the role of differential transport through the interface is explored together with the presence of asymmetry between the homogeneous equilibria of the two layers. A linear stability analysis is carried out around a spatially homogeneous equilibrium state in the asymptotic limit of weak and strong coupling strengths, where quantitative approximations of the bifurcation curve can be computed. Our theoretical findings, for an arbitrary number of reacting species, reveal quantitative Turing conditions, highlighting when the coupling mechanism between the layered regions can either trigger patterning or stabilize a spatially homogeneous equilibrium regardless of the independent patterning state of each layer. We support our theoretical results through direct numerical simulations, and provide an open source code to explore such systems further.

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.

Getting the most out of maths: How to coordinate mathematical modelling research to support a pandemic, lessons learnt from three initiatives that were part of the COVID-19 response in the UK.

In March 2020 mathematics became a key part of the scientific advice to the UK government on the pandemic response to COVID-19. Mathematical and statistical modelling provided critical information on the spread of the virus and the potential impact of different interventions. The unprecedented scale of the challenge led the epidemiological modelling community in the UK to be pushed to its limits. At the same time, mathematical modellers across the country were keen to use their knowledge and skills to support the COVID-19 modelling effort. However, this sudden great interest in epidemiological modelling needed to be coordinated to provide much-needed support, and to limit the burden on epidemiological modellers already very stretched for time. In this paper we describe three initiatives set up in the UK in spring 2020 to coordinate the mathematical sciences research community in supporting mathematical modelling of COVID-19. Each initiative had different primary aims and worked to maximise synergies between the various projects. We reflect on the lessons learnt, highlighting the key roles of pre-existing research collaborations and focal centres of coordination in contributing to the success of these initiatives. We conclude with recommendations about important ways in which the scientific research community could be better prepared for future pandemics. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Investigating the Influence of Growth Arrest Mechanisms on Tumour Responses to Radiotherapy.

Cancer is a heterogeneous disease and tumours of the same type can differ greatly at the genetic and phenotypic levels. Understanding how these differences impact sensitivity to treatment is an essential step towards patient-specific treatment design. In this paper, we investigate how two different mechanisms for growth control may affect tumour cell responses to fractionated radiotherapy (RT) by extending an existing ordinary differential equation model of tumour growth. In the absence of treatment, this model distinguishes between growth arrest due to nutrient insufficiency and competition for space and exhibits three growth regimes: nutrient limited, space limited (SL) and bistable (BS), where both mechanisms for growth arrest coexist. We study the effect of RT for tumours in each regime, finding that tumours in the SL regime typically respond best to RT, while tumours in the BS regime typically respond worst to RT. For tumours in each regime, we also identify the biological processes that may explain positive and negative treatment outcomes and the dosing regimen which maximises the reduction in tumour burden.

Abnormal morphology biases hematocrit distribution in tumor vasculature and contributes to heterogeneity in tissue oxygenation.

Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal vascular structure of the tumor, but the precise mechanisms linking abnormal structure and compromised oxygen transport are only partially understood. In this paper, we investigate the role that red blood cell (RBC) transport plays in establishing oxygen heterogeneity in tumor tissue. We focus on heterogeneity driven by network effects, which are challenging to observe experimentally due to the reduced fields of view typically considered. Motivated by our findings of abnormal vascular patterns linked to deviations from current RBC transport theory, we calculated average vessel lengths [Formula: see text] and diameters [Formula: see text] from tumor allografts of three cancer cell lines and observed a substantial reduction in the ratio [Formula: see text] compared to physiological conditions. Mathematical modeling reveals that small values of the ratio λ (i.e., [Formula: see text]) can bias hematocrit distribution in tumor vascular networks and drive heterogeneous oxygenation of tumor tissue. Finally, we show an increase in the value of λ in tumor vascular networks following treatment with the antiangiogenic cancer agent DC101. Based on our findings, we propose λ as an effective way of monitoring the efficacy of antiangiogenic agents and as a proxy measure of perfusion and oxygenation in tumor tissue undergoing antiangiogenic treatment.

A DNA-structured mathematical model of cell-cycle progression in cyclic hypoxia.

New experimental data have shown how the periodic exposure of cells to low oxygen levels (i.e., cyclic hypoxia) impacts their progress through the cell-cycle. Cyclic hypoxia has been detected in tumours and linked to poor prognosis and treatment failure. While fluctuating oxygen environments can be reproduced in vitro, the range of oxygen cycles that can be tested is limited. By contrast, mathematical models can be used to predict the response to a wide range of cyclic dynamics. Accordingly, in this paper we develop a mechanistic model of the cell-cycle that can be combined with in vitro experiments to better understand the link between cyclic hypoxia and cell-cycle dysregulation. A distinguishing feature of our model is the inclusion of impaired DNA synthesis and cell-cycle arrest due to periodic exposure to severely low oxygen levels. Our model decomposes the cell population into five compartments and a time-dependent delay accounts for the variability in the duration of the S phase which increases in severe hypoxia due to reduced rates of DNA synthesis. We calibrate our model against experimental data and show that it recapitulates the observed cell-cycle dynamics. We use the calibrated model to investigate the response of cells to oxygen cycles not yet tested experimentally. When the re-oxygenation phase is sufficiently long, our model predicts that cyclic hypoxia simply slows cell proliferation since cells spend more time in the S phase. On the contrary, cycles with short periods of re-oxygenation are predicted to lead to inhibition of proliferation, with cells arresting from the cell-cycle in the G2 phase. While model predictions on short time scales (about a day) are fairly accurate (i.e, confidence intervals are small), the predictions become more uncertain over longer periods. Hence, we use our model to inform experimental design that can lead to improved model parameter estimates and validate model predictions.

A multiscale model of complex endothelial cell dynamics in early angiogenesis.

We introduce a hybrid two-dimensional multiscale model of angiogenesis, the process by which endothelial cells (ECs) migrate from a pre-existing vascular bed in response to local environmental cues and cell-cell interactions, to create a new vascular network. Recent experimental studies have highlighted a central role of cell rearrangements in the formation of angiogenic networks. Our model accounts for this phenomenon via the heterogeneous response of ECs to their microenvironment. These cell rearrangements, in turn, dynamically remodel the local environment. The model reproduces characteristic features of angiogenic sprouting that include branching, chemotactic sensitivity, the brush border effect, and cell mixing. These properties, rather than being hardwired into the model, emerge naturally from the gene expression patterns of individual cells. After calibrating and validating our model against experimental data, we use it to predict how the structure of the vascular network changes as the baseline gene expression levels of the VEGF-Delta-Notch pathway, and the composition of the extracellular environment, vary. In order to investigate the impact of cell rearrangements on the vascular network structure, we introduce the mixing measure, a scalar metric that quantifies cell mixing as the vascular network grows. We calculate the mixing measure for the simulated vascular networks generated by ECs of different lineages (wild type cells and mutant cells with impaired expression of a specific receptor). Our results show that the time evolution of the mixing measure is directly correlated to the generic features of the vascular branching pattern, thus, supporting the hypothesis that cell rearrangements play an essential role in sprouting angiogenesis. Furthermore, we predict that lower cell rearrangement leads to an imbalance between branching and sprout elongation. Since the computation of this statistic requires only individual cell trajectories, it can be computed for networks generated in biological experiments, making it a potential biomarker for pathological angiogenesis.

Combining Mechanisms of Growth Arrest in Solid Tumours: A Mathematical Investigation.

The processes underpinning solid tumour growth involve the interactions between various healthy and tumour tissue components and the vasculature, and can be affected in different ways by cancer treatment. In particular, the growth-limiting mechanisms at play may influence tumour responses to treatment. In this paper, we propose a simple ordinary differential equation model of solid tumour growth to investigate how tumour-specific mechanisms of growth arrest may affect tumour response to different combination cancer therapies. We consider the interactions of tumour cells with the physical space in which they proliferate and a nutrient supplied by the tumour vasculature, with the aim of representing two distinct growth arrest mechanisms. More specifically, we wish to consider growth arrest due to (1) nutrient deficiency, which corresponds to balancing cell proliferation and death rates, and (2) competition for space, which corresponds to cessation of proliferation without cell death. We perform numerical simulations of the model and a steady-state analysis to determine the possible tumour growth scenarios described by the model. We find that there are three distinct growth regimes: the nutrient- and spatially limited regimes and a bi-stable regime, in which both growth arrest mechanisms are simultaneously active. Thus, the proposed model has the features required to investigate and distinguish tumour responses to different cancer treatments.

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.

Inferring Tumor Proliferative Organization from Phylogenetic Tree Measures in a Computational Model.

We use a computational modeling approach to explore whether it is possible to infer a solid tumor's cellular proliferative hierarchy under the assumptions of the cancer stem cell hypothesis and neutral evolution. We work towards inferring the symmetric division probability for cancer stem cells, since this is believed to be a key driver of progression and therapeutic response. Motivated by the advent of multiregion sampling and resulting opportunities to infer tumor evolutionary history, we focus on a suite of statistical measures of the phylogenetic trees resulting from the tumor's evolution in different regions of parameter space and through time. We find strikingly different patterns in these measures for changing symmetric division probability which hinge on the inclusion of spatial constraints. These results give us a starting point to begin stratifying tumors by this biological parameter and also generate a number of actionable clinical and biological hypotheses regarding changes during therapy, and through tumor evolutionary time. [Cancer; evolution; phylogenetics.].

Modern perspectives on near-equilibrium analysis of Turing systems.

In the nearly seven decades since the publication of Alan Turing's work on morphogenesis, enormous progress has been made in understanding both the mathematical and biological aspects of his proposed reaction-diffusion theory. Some of these developments were nascent in Turing's paper, and others have been due to new insights from modern mathematical techniques, advances in numerical simulations and extensive biological experiments. Despite such progress, there are still important gaps between theory and experiment, with many examples of biological patterning where the underlying mechanisms are still unclear. Here, we review modern developments in the mathematical theory pioneered by Turing, showing how his approach has been generalized to a range of settings beyond the classical two-species reaction-diffusion framework, including evolving and complex manifolds, systems heterogeneous in space and time, and more general reaction-transport equations. While substantial progress has been made in understanding these more complicated models, there are many remaining challenges that we highlight throughout. We focus on the mathematical theory, and in particular linear stability analysis of 'trivial' base states. We emphasize important open questions in developing this theory further, and discuss obstacles in using these techniques to understand biological reality. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Introduction to 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Elucidating pattern forming processes is an important problem in the physical, chemical and biological sciences. Turing's contribution, after being initially neglected, eventually catalysed a huge amount of work from mathematicians, physicists, chemists and biologists aimed towards understanding how steady spatial patterns can emerge from homogeneous chemical mixtures due to the reaction and diffusion of different chemical species. While this theory has been developed mathematically and investigated experimentally for over half a century, many questions still remain unresolved. This theme issue places Turing's theory of pattern formation in a modern context, discussing the current frontiers in foundational aspects of pattern formation in reaction-diffusion and related systems. It highlights ongoing work in chemical, synthetic and developmental settings which is helping to elucidate how important Turing's mechanism is for real morphogenesis, while highlighting gaps that remain in matching theory to reality. The theme issue also surveys a variety of recent mathematical research pushing the boundaries of Turing's original theory to more realistic and complicated settings, as well as discussing open theoretical challenges in the analysis of such models. It aims to consolidate current research frontiers and highlight some of the most promising future directions. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Isolating Patterns in Open Reaction-Diffusion Systems.

Realistic examples of reaction-diffusion phenomena governing spatial and spatiotemporal pattern formation are rarely isolated systems, either chemically or thermodynamically. However, even formulations of 'open' reaction-diffusion systems often neglect the role of domain boundaries. Most idealizations of closed reaction-diffusion systems employ no-flux boundary conditions, and often patterns will form up to, or along, these boundaries. Motivated by boundaries of patterning fields related to the emergence of spatial form in embryonic development, we propose a set of mixed boundary conditions for a two-species reaction-diffusion system which forms inhomogeneous solutions away from the boundary of the domain for a variety of different reaction kinetics, with a prescribed uniform state near the boundary. We show that these boundary conditions can be derived from a larger heterogeneous field, indicating that these conditions can arise naturally if cell signalling or other properties of the medium vary in space. We explain the basic mechanisms behind this pattern localization and demonstrate that it can capture a large range of localized patterning in one, two, and three dimensions and that this framework can be applied to systems involving more than two species. Furthermore, the boundary conditions proposed lead to more symmetrical patterns on the interior of the domain and plausibly capture more realistic boundaries in developmental systems. Finally, we show that these isolated patterns are more robust to fluctuations in initial conditions and that they allow intriguing possibilities of pattern selection via geometry, distinct from known selection mechanisms.

Comparative analysis of continuum angiogenesis models.

Although discrete approaches are increasingly employed to model biological phenomena, it remains unclear how complex, population-level behaviours in such frameworks arise from the rules used to represent interactions between individuals. Discrete-to-continuum approaches, which are used to derive systems of coarse-grained equations describing the mean-field dynamics of a microscopic model, can provide insight into such emergent behaviour. Coarse-grained models often contain nonlinear terms that depend on the microscopic rules of the discrete framework, however, and such nonlinearities can make a model difficult to mathematically analyse. By contrast, models developed using phenomenological approaches are typically easier to investigate but have a more obscure connection to the underlying microscopic system. To our knowledge, there has been little work done to compare solutions of phenomenological and coarse-grained models. Here we address this problem in the context of angiogenesis (the creation of new blood vessels from existing vasculature). We compare asymptotic solutions of a classical, phenomenological "snail-trail" model for angiogenesis to solutions of a nonlinear system of partial differential equations (PDEs) derived via a systematic coarse-graining procedure (Pillay et al. in Phys Rev E 95(1):012410, 2017. https://doi.org/10.1103/PhysRevE.95.012410 ). For distinguished parameter regimes corresponding to chemotaxis-dominated cell movement and low branching rates, both continuum models reduce at leading order to identical PDEs within the domain interior. Numerical and analytical results confirm that pointwise differences between solutions to the two continuum models are small if these conditions hold, and demonstrate how perturbation methods can be used to determine when a phenomenological model provides a good approximation to a more detailed coarse-grained system for the same biological process.