Turing Patterning in Stratified Domains

Spatiotemporal patterns in coupled reaction-diffusion systems with nonidentical kinetics

Understanding of the effect of coupling interaction is at the heart of nonlinear science since some nonequilibrium systems are composed of different layers or units. In this paper, we demonstrate various spatio-temporal patterns in a nonlinearly coupled two-layer Turing system with nonidentical reaction kinetics. Both the type of Turing mode and coupling form play an important role in the pattern formation and pattern selection. Two kinds of Turing mode interactions, namely supercritical-subcritical and supercritical-supercritical Turing mode interaction, have been investigated. Stationary resonant superlattice patterns arise spontaneously in both cases, while dynamic patterns can also be formed in the latter case. The destabilization of spike solutions induced by spatial heterogeneity may be responsible for these dynamic patterns. In contrast to linear coupling, the nonlinear coupling not only increases the complexity of spatio-temporal patterns, but also reduces the requirements of spatial resonance conditions. The simulation results are in good agreement with the experimental observations in dielectric barrier discharge systems.

Linear coupling of patterning systems can have nonlinear effects

Isolated patterning systems have been repeatedly investigated. However, biological systems rarely work on their own. This paper presents a theoretical and quantitative analysis of a two-domain interconnected geometry, or bilayer, coupling two two-species reaction-diffusion systems mimicking interlayer communication, such as in mammary organoids. Each layer has identical kinetics and parameters, but differing diffusion coefficients. Critically, we show that despite a linear coupling between the layers, the model demonstrates nonlinear behavior; the coupling can lead to pattern suppression or pattern enhancement. Using the Routh-Hurwitz stability criterion multiple times, we investigate the pattern-forming capabilities of the uncoupled system, the weakly coupled system, and the strongly coupled system, using numerical simulations to back up the analysis. We show that although the dispersion relation of the entire system is a nontrivial octic polynomial, the patterning wave modes in the strongly coupled case can be approximated by a quartic polynomial, whose features are easier to understand.

A three-node Turing gene circuit forms periodic spatial patterns in bacteria

Turing patterns are self-organizing systems that can form spots, stripes, or labyrinths. Proposed examples in tissue organization include zebrafish pigmentation, digit spacing, and many others. The theory of Turing patterns in biology has been debated because of their stringent fine-tuning requirements, where patterns only occur within a small subset of parameters. This has complicated the engineering of synthetic Turing gene circuits from first principles, although natural genetic Turing networks have been identified. Here, we engineered a synthetic genetic reaction-diffusion system where three nodes interact according to a non-classical Turing network with improved parametric robustness. The system reproducibly generated stationary, periodic, concentric stripe patterns in growing E. coli colonies. A partial differential equation model reproduced the patterns, with a Turing parameter regime obtained by fitting to experimental data. Our synthetic Turing system can contribute to nanotechnologies, such as patterned biomaterial deposition, and provide insights into developmental patterning programs. A record of this paper's transparent peer review process is included in the supplemental information.

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.

Boundary Conditions Cause Different Generic Bifurcation Structures in Turing Systems

Turing’s theory of morphogenesis is a generic mechanism to produce spatial patterning from near homogeneity. Although widely studied, we are still able to generate new results by returning to common dogmas. One such widely reported belief is that the Turing bifurcation occurs through a pitchfork bifurcation, which is true under zero-flux boundary conditions. However, under fixed boundary conditions, the Turing bifurcation becomes generically transcritical. We derive these algebraic results through weakly nonlinear analysis and apply them to the Schnakenberg kinetics. We observe that the combination of kinetics and boundary conditions produce their own uncommon boundary complexities that we explore numerically. Overall, this work demonstrates that it is not enough to only consider parameter perturbations in a sensitivity analysis of a specific application. Variations in boundary conditions should also be considered.

Fixed and Distributed Gene Expression Time Delays in Reaction-Diffusion Systems

Time delays, modelling the process of intracellular gene expression, have been shown to have important impacts on the dynamics of pattern formation in reaction-diffusion systems. In particular, past work has shown that such time delays can shrink the Turing space, thereby inhibiting patterns from forming across large ranges of parameters. Such delays can also increase the time taken for pattern formation even when Turing instabilities occur. Here, we consider reaction-diffusion models incorporating fixed or distributed time delays, modelling the underlying stochastic nature of gene expression dynamics, and analyse these through a systematic linear instability analysis and numerical simulations for several sets of different reaction kinetics. We find that even complicated distribution kernels (skewed Gaussian probability density functions) have little impact on the reaction-diffusion dynamics compared to fixed delays with the same mean delay. We show that the location of the delay terms in the model can lead to changes in the size of the Turing space (increasing or decreasing) as the mean time delay, [Formula: see text], is increased. We show that the time to pattern formation from a perturbation of the homogeneous steady state scales linearly with [Formula: see text], and conjecture that this is a general impact of time delay on reaction-diffusion dynamics, independent of the form of the kinetics or location of the delayed terms. Finally, we show that while initial and boundary conditions can influence these dynamics, particularly the time-to-pattern, the effects of delay appear robust under variations of initial and boundary data. Overall, our results help clarify the role of gene expression time delays in reaction-diffusion patterning, and suggest clear directions for further work in studying more realistic models of pattern formation.

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'.

Reaction-diffusion models for morphological patterning of hESCs

In this paper we consider mathematical modeling of the dynamics of self-organized patterning of spatially confined human embryonic stem cells (hESCs) treated with BMP4 (gastruloids) described in recent experimental works (Warmflash in Nat Methods 11:847-854, 2014; Chhabra in PloS Biol 17: 3000498, 2019). In the first part of the paper we use the activator-inhibitor equations of Gierer and Meinhardt to identify 3 reaction-diffusion regimes for each of the three morphogenic proteins, BMP4, Wnt and Nodal, based on the characteristic features of the dynamic patterning. We identify appropriate boundary conditions which correspond to the experimental setup and perform numerical simulations of the reaction-diffusion (RD) systems, using the finite element approximation, to confirm that the RD systems in these regimes produce realistic dynamics of the protein concentrations. In the second part of the paper we use analytic tools to address the questions of the existence and stability of non-homogeneous steady states for the reaction-diffusion systems of the type considered in the first part of the paper.

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.