Emergent structures in reaction-advection-diffusion systems on a sphere

Modelling of surface reactions and diffusion mediated by bulk diffusion

We develop a continuum framework applicable to solid-state hydrogen storage, cell biology and other scenarios where the diffusion of a single constituent within a bulk region is coupled via adsorption/desorption to reactions and diffusion on the boundary of the region. We formulate content balances for all relevant constituents and develop thermodynamically consistent constitutive equations. The latter encompass two classes of kinetics for adsorption/desorption and chemical reactions-fast and Marcelin-De Donder, and the second class includes mass action kinetics as a special case. We apply the framework to derive a system consisting of the standard diffusion equation in bulk and FitzHugh-Nagumo type surface reaction-diffusion system of equations on the boundary. We also study the linear stability of a homogeneous steady state in a spherical region and establish sufficient conditions for the occurrence of instabilities driven by surface diffusion. These findings are verified through numerical simulations which reveal that instabilities driven by diffusion lead to the emergence of steady-state spatial patterns from random initial conditions and that bulk diffusion can suppress spatial patterns, in which case temporal oscillations can ensue. We include an extension of our framework that accounts for mechanochemical coupling when the bulk region is occupied by a deformable solid. This article is part of the theme issue 'Foundational issues, analysis and geometry in continuum mechanics'.

A two-galectin network establishes mesenchymal condensation phenotype in limb development

Previous work showed that Gal-1A and Gal-8, two proteins belonging to the galactoside-binding galectin family, are the earliest determinants of the patterning of the skeletal elements of embryonic chicken limbs, and further, that their experimentally determined interactions in the embryonic limb bud can be interpreted via a reaction-diffusion-adhesion (2GL: two galectin plus ligands) model. Here, we use an ordinary differential equation-based approach to analyze the intrinsic switching modality of the 2GL network and characterize the network behavior independent of the diffusive and adhesive arms of the patterning mechanism. We identify two states: where the concentrations of both the galectins are respectively, negligible, and very high. This bistable switch-like system arises via a saddle-node bifurcation from a monostable state. For the case of mass-action production terms, we provide an explicit Lyapunov function for the system, which shows that it has no periodic solutions. Our model therefore predicts that the galectin network may exist in low expression and high expression states separated in space or time, without any intermediate states. We test these predictions in experiments performed with high density cultures of chick limb mesenchymal cells and observe that cells inside precartilage protocondensations express Gal-1A at a much higher rate than those outside, for which it was negligible. The Gal-1A and -8-based patterning network is therefore sufficient to partition the mesenchymal cell population into two discrete cell states with different developmental (chondrogenic vs. non-chondrogenic) fates. When incorporated into an adhesion and diffusion-enabled framework this system can generate a spatially patterned limb skeleton.

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

Bespoke Turing Systems

Reaction-diffusion systems are an intensively studied form of partial differential equation, frequently used to produce spatially heterogeneous patterned states from homogeneous symmetry breaking via the Turing instability. Although there are many prototypical "Turing systems" available, determining their parameters, functional forms, and general appropriateness for a given application is often difficult. Here, we consider the reverse problem. Namely, suppose we know the parameter region associated with the reaction kinetics in which patterning is required-we present a constructive framework for identifying systems that will exhibit the Turing instability within this region, whilst in addition often allowing selection of desired patterning features, such as spots, or stripes. In particular, we show how to build a system of two populations governed by polynomial morphogen kinetics such that the: patterning parameter domain (in any spatial dimension), morphogen phases (in any spatial dimension), and even type of resulting pattern (in up to two spatial dimensions) can all be determined. Finally, by employing spatial and temporal heterogeneity, we demonstrate that mixed mode patterns (spots, stripes, and complex prepatterns) are also possible, allowing one to build arbitrarily complicated patterning landscapes. Such a framework can be employed pedagogically, or in a variety of contemporary applications in designing synthetic chemical and biological patterning systems. We also discuss the implications that this freedom of design has on using reaction-diffusion systems in biological modelling and suggest that stronger constraints are needed when linking theory and experiment, as many simple patterns can be easily generated given freedom to choose reaction kinetics.

Turing conditions for pattern forming systems on evolving manifolds

The study of pattern-forming instabilities in reaction-diffusion systems on growing or otherwise time-dependent domains arises in a variety of settings, including applications in developmental biology, spatial ecology, and experimental chemistry. Analyzing such instabilities is complicated, as there is a strong dependence of any spatially homogeneous base states on time, and the resulting structure of the linearized perturbations used to determine the onset of instability is inherently non-autonomous. We obtain general conditions for the onset and structure of diffusion driven instabilities in reaction-diffusion systems on domains which evolve in time, in terms of the time-evolution of the Laplace-Beltrami spectrum for the domain and functions which specify the domain evolution. Our results give sufficient conditions for diffusive instabilities phrased in terms of differential inequalities which are both versatile and straightforward to implement, despite the generality of the studied problem. These conditions generalize a large number of results known in the literature, such as the algebraic inequalities commonly used as a sufficient criterion for the Turing instability on static domains, and approximate asymptotic results valid for specific types of growth, or specific domains. We demonstrate our general Turing conditions on a variety of domains with different evolution laws, and in particular show how insight can be gained even when the domain changes rapidly in time, or when the homogeneous state is oscillatory, such as in the case of Turing-Hopf instabilities. Extensions to higher-order spatial systems are also included as a way of demonstrating the generality of the approach.

A Non-local Cross-Diffusion Model of Population Dynamics I: Emergent Spatial and Spatiotemporal Patterns

We extend a spatially non-local cross-diffusion model of aggregation between multiple species with directed motion toward resource gradients to include many species and more general kinds of dispersal. We first consider diffusive instabilities, determining that for directed motion along fecundity gradients, the model permits the Turing instability leading to colony formation and persistence provided there are three or more interacting species. We also prove that such patterning is not possible in the model under the Turing mechanism for two species under directed motion along fecundity gradients, confirming earlier findings in the literature. However, when the directed motion is not along fecundity gradients, for instance, if foraging or migration is sub-optimal relative to fecundity gradients, we find that very different colony structures can emerge. This generalization also permits colony formation for two interacting species. In the advection-dominated case, aggregation patterns are more broad and global in nature, due to the inherent non-local nature of the advection which permits directed motion over greater distances, whereas in the diffusion-dominated case, more highly localized patterns and colonies develop, owing to the localized nature of random diffusion. We also consider the interplay between Turing patterning and spatial heterogeneity in resources. We find that for small spatial variations, there will be a combination of Turing patterns and patterning due to spatial forcing from the resources, whereas for large resource variations, spatial or spatiotemporal patterning can be modified greatly from what is predicted on homogeneous domains. For each of these emergent behaviors, we outline the theoretical mechanism leading to colony formation and then provide numerical simulations to illustrate the results. We also discuss implications this model has for studies of directed motion in different ecological settings.

Influence of Curvature, Growth, and Anisotropy on the Evolution of Turing Patterns on Growing Manifolds

We study two-species reaction–diffusion systems on growing manifolds, including situations where the growth is anisotropic yet dilational in nature. In contrast to the literature on linear instabilities in such systems, we study how growth and anisotropy impact the qualitative properties of nonlinear patterned states which have formed before growth is initiated. We produce numerical solutions to numerous reaction–diffusion systems with varying reaction kinetics, manner of growth (both isotropic and anisotropic), and timescales of growth on both planar elliptical and curved ellipsoidal domains. We find that in some parameter regimes, some of these factors have a negligible effect on the long-time patterned state. On the other hand, we find that some of these factors play a role in determining the patterns formed on surfaces and that anisotropic growth can produce qualitatively different patterns to those formed under isotropic growth.

Turing-Hopf patterns on growing domains: The torus and the sphere

This paper deals with the study of spatial and spatio-temporal patterns in the reaction-diffusion FitzHugh-Nagumo model on growing curved domains. This is carried out on two exemplar cases: a torus and a sphere. We compute bifurcation boundaries for the homogeneous steady state when the homogeneous system is monostable. We exhibit Turing and Turing-Hopf bifurcations, as well as additional patterning outside of these bifurcation regimes due to the multistability of patterned states. We consider static and growing domains, where the growth is slow, isotropic, and exponential in time, allowing for a simple analytical calculation of these bifurcations in terms of model parameters. Numerical simulations allow us to discuss the role played by the growth and the curvature of the domains on the pattern selection on the torus and the sphere. We demonstrate parameter regimes where the linear theory can successfully predict the kind of pattern (homogeneous and heterogeneous oscillations and stationary spatial patterns) but not their detailed nonlinear structure. We also find parameter regimes where the linear theory fails, such as Hopf regimes which give rise to spatial patterning (depending on geometric details), where we suspect that multistability plays a key role in the departure from homogeneity. Finally we also demonstrate effects due to the evolution of nonuniform patterns under growth, suggesting important roles for growth in reaction-diffusion systems beyond modifying instability regimes.

Heterogeneity induces spatiotemporal oscillations in reaction-diffusion systems

We report on an instability arising in activator-inhibitor reaction-diffusion (RD) systems with a simple spatial heterogeneity. This instability gives rise to periodic creation, translation, and destruction of spike solutions that are commonly formed due to Turing instabilities. While this behavior is oscillatory in nature, it occurs purely within the Turing space such that no region of the domain would give rise to a Hopf bifurcation for the homogeneous equilibrium. We use the shadow limit of the Gierer-Meinhardt system to show that the speed of spike movement can be predicted from well-known asymptotic theory, but that this theory is unable to explain the emergence of these spatiotemporal oscillations. Instead, we numerically explore this system and show that the oscillatory behavior is caused by the destabilization of a steady spike pattern due to the creation of a new spike arising from endogeneous activator production. We demonstrate that on the edge of this instability, the period of the oscillations goes to infinity, although it does not fit the profile of any well-known bifurcation of a limit cycle. We show that nearby stationary states are either Turing unstable or undergo saddle-node bifurcations near the onset of the oscillatory instability, suggesting that the periodic motion does not emerge from a local equilibrium. We demonstrate the robustness of this spatiotemporal oscillation by exploring small localized heterogeneity and showing that this behavior also occurs in the Schnakenberg RD model. Our results suggest that this phenomenon is ubiquitous in spatially heterogeneous RD systems, but that current tools, such as stability of spike solutions and shadow-limit asymptotics, do not elucidate understanding. This opens several avenues for further mathematical analysis and highlights difficulties in explaining how robust patterning emerges from Turing's mechanism in the presence of even small spatial heterogeneity.