Instability of turing patterns in reaction-diffusion-ODE systems

Nonlocal models in biology and life sciences: Sources, developments, and applications

Mathematical modeling is one of the fundamental techniques for understanding biophysical mechanisms in developmental biology. It helps researchers to analyze complex physiological processes and connect like a bridge between theoretical and experimental observations. Various groups of mathematical models have been studied to analyze these processes, and the nonlocal models are one of them. Nonlocality is important in realistic mathematical models of physical and biological systems when local models fail to capture the essential dynamics and interactions that occur over a range of distances (e.g., cell-cell, cell-tissue adhesions, neural networks, the spread of diseases, intra-specific competition, nanobeams, etc.). This review illustrates different nonlocal mathematical models applied to biology and life sciences. The major focus has been given to sources, developments, and applications of such models. Among other things, a systematic discussion has been provided for the conditions of pattern formations in biological systems of population dynamics. Special attention has also been given to nonlocal interactions on networks, network coupling and integration, including brain dynamics models that provide an important tool to understand neurodegenerative diseases better. In addition, we have discussed nonlocal modeling approaches for cancer stem cells and tumor cells that are widely applied in the cell migration processes, growth, and avascular tumors in any organ. Furthermore, the discussed nonlocal continuum models can go sufficiently smaller scales, including nanotechnology, where classical local models often fail to capture the complexities of nanoscale interactions, applied to build biosensors to sense biomaterial and its concentration. Piezoelectric and other smart materials are among them, and these devices are becoming increasingly important in the digital and physical world that is intrinsically interconnected with biological systems. Additionally, we have reviewed a nonlocal theory of peridynamics, which deals with continuous and discrete media and applies to model the relationship between fracture and healing in cortical bone, tissue growth and shrinkage, and other areas increasingly important in biomedical and bioengineering applications. Finally, we provided a comprehensive summary of emerging trends and highlighted future directions in this rapidly expanding field.

Linking discrete and continuous models of cell birth and migration

Self-organization of individuals within large collectives occurs throughout biology. Mathematical models can help elucidate the individual-level mechanisms behind these dynamics, but analytical tractability often comes at the cost of biological intuition. Discrete models provide straightforward interpretations by tracking each individual yet can be computationally expensive. Alternatively, continuous models supply a large-scale perspective by representing the 'effective' dynamics of infinite agents, but their results are often difficult to translate into experimentally relevant insights. We address this challenge by quantitatively linking spatio-temporal dynamics of continuous models and individual-based data in settings with biologically realistic, time-varying cell numbers. Specifically, we introduce and fit scaling parameters in continuous models to account for discrepancies that can arise from low cell numbers and localized interactions. We illustrate our approach on an example motivated by zebrafish-skin pattern formation, in which we create a continuous framework describing the movement and proliferation of a single cell population by upscaling rules from a discrete model. Our resulting continuous models accurately depict ensemble average agent-based solutions when migration or proliferation act alone. Interestingly, the same parameters are not optimal when both processes act simultaneously, highlighting a rich difference in how combining migration and proliferation affects discrete and continuous dynamics.

Travelling waves due to negative plant-soil feedbacks in a model including tree life-stages

The emergence and maintenance of tree species diversity in tropical forests is commonly attributed to the Janzen-Connell (JC) hypothesis, which states that growth of seedlings is suppressed in the proximity of conspecific adult trees. As a result, a JC distribution due to a density-dependent negative feedback emerges in the form of a (transient) pattern where conspecific seedling density is highest at intermediate distances away from parent trees. Several studies suggest that the required density-dependent feedbacks behind this pattern could result from interactions between trees and soil-borne pathogens. However, negative plant-soil feedback may involve additional mechanisms, including the accumulation of autotoxic compounds generated through tree litter decomposition. An essential task therefore consists in constructing mathematical models incorporating both effects showing the ability to support the emergence of JC distributions. In this work, we develop and analyse a novel reaction-diffusion-ODE model, describing the interactions within tropical tree species across different life stages (seeds, seedlings, and adults) as driven by negative plant-soil feedback. In particular, we show that under strong negative plant-soil feedback travelling wave solutions exist, creating transient distributions of adult trees and seedlings that are in agreement with the Janzen-Connell hypothesis. Moreover, we show that these travelling wave solutions are pulled fronts and a robust feature as they occur over a broad parameter range. Finally, we calculate their linear spreading speed and show its (in)dependence on relevant nondimensional parameters.

Patterns governed by chemotaxis and time delay

In this paper, sufficient conditions of Turing instability are established for general delayed reaction-diffusion-chemotaxis models with no-flux boundary conditions. In particular, we address the difficulty brought about by the time delay in investigating the Turing instability. These models involve the time delay parameter τ and the general density (concentration) function in chemotaxis terms with the chemotaxis parameter χ. Theoretical results reveal that the time delay parameter τ could determine the stability of positive equilibrium for ordinary differential equations, while the chemotaxis parameter χ could describe the stability of positive equilibrium for partial differential equations. In this fashion, the general conditions of Turing instability are presented. To confirm their validity, the delayed chemotaxis-type predator-prey model and the phytoplankton-zooplankton model are considered. It is found that these two models admit Turing instability, and the numerical results are in good agreement with the theoretical analysis. The obtained results are helpful in the application of the Turing pattern in models with time delay and chemotaxis.

Turing patterns in simplicial complexes

The spontaneous emergence of patterns in nature, such as stripes and spots, can be mathematically explained by reaction-diffusion systems. These patterns are often referred as Turing patterns to honor the seminal work of Alan Turing in the early 1950s. With the coming of age of network science, and with its related departure from diffusive nearest-neighbor interactions to long-range links between nodes, additional layers of complexity behind pattern formation have been discovered, including irregular spatiotemporal patterns. Here we investigate the formation of Turing patterns in simplicial complexes, where links no longer connect just pairs of nodes but can connect three or more nodes. Such higher-order interactions are emerging as a new frontier in network science, in particular describing group interaction in various sociological and biological systems, so understanding pattern formation under these conditions is of the utmost importance. We show that a canonical reaction-diffusion system defined over a simplicial complex yields Turing patterns that fundamentally differ from patterns observed in traditional networks. For example, we observe a stable distribution of Turing patterns where the fraction of nodes with reactant concentrations above the equilibrium point is exponentially related to the average degree of 2-simplexes, and we uncover parameter regions where Turing patterns will emerge only under higher-order interactions, but not under pairwise interactions.

Mathematical models of nitrogen-fixing cell patterns in filamentous cyanobacteria

The Anabaena genus is a model organism of filamentous cyanobacteria whose vegetative cells can differentiate under nitrogen-limited conditions into a type of cell called a heterocyst. These heterocysts lose the possibility to divide and are necessary for the filament because they can fix and share environmental nitrogen. In order to distribute the nitrogen efficiently, heterocysts are arranged to form a quasi-regular pattern whose features are maintained as the filament grows. Recent efforts have allowed advances in the understanding of the interactions and genetic mechanisms underlying this dynamic pattern. Here, we present a systematic review of the existing theoretical models of nitrogen-fixing cell differentiation in filamentous cyanobacteria. These filaments constitute one of the simplest forms of multicellular organization, and this allows for several modeling scales of this emergent pattern. The system has been approached at three different levels. From bigger to smaller scale, the system has been considered as follows: at the population level, by defining a mean-field simplified system to study the ratio of heterocysts and vegetative cells; at the filament level, with a continuous simplification as a reaction-diffusion system; and at the cellular level, by studying the genetic regulation that produces the patterning for each cell. In this review, we compare these different approaches noting both the virtues and shortcomings of each one of them.

Pattern Formation in a Spatially Extended Model of Pacemaker Dynamics in Smooth Muscle Cells

Spatiotemporal patterns are common in biological systems. For electrically coupled cells, previous studies of pattern formation have mainly used applied current as the primary bifurcation parameter. The purpose of this paper is to show that applied current is not needed to generate spatiotemporal patterns for smooth muscle cells. The patterns can be generated solely by external mechanical stimulation (transmural pressure). To do this we study a reaction-diffusion system involving the Morris-Lecar equations and observe a wide range of spatiotemporal patterns for different values of the model parameters. Some aspects of these patterns are explained via a bifurcation analysis of the system without coupling - in particular Type I and Type II excitability both occur. We show the patterns are not due to a Turing instability and that the spatially extended model exhibits spatiotemporal chaos. We also use travelling wave coordinates to analyse travelling waves.

Optimal control of networked reaction-diffusion systems

Patterns in nature are fascinating both aesthetically and scientifically. Alan Turing's celebrated reaction-diffusion model of pattern formation from the 1950s has been extended to an astounding diversity of applications: from cancer medicine, via nanoparticle fabrication, to computer architecture. Recently, several authors have studied pattern formation in underlying networks, but thus far, controlling a reaction-diffusion system in a network to obtain a particular pattern has remained elusive. We present a solution to this problem in the form of an analytical framework and numerical algorithm for optimal control of Turing patterns in networks. We demonstrate our method's effectiveness and discuss factors that affect its performance. We also pave the way for multidisciplinary applications of our framework beyond reaction-diffusion models.

Robust controlled formation of Turing patterns in three-component systems

Over the past few decades, formation of Turing patterns in reaction-diffusion systems has been shown to be the underlying process in several examples of biological morphogenesis, confirming Alan Turing's hypothesis, put forward in 1952. However, theoretical studies suggest that Turing patterns formation via classical "short-range activation and long-range inhibition" concept in general can happen within only narrow parameter ranges. This feature seemingly contradicts the accuracy and reproducibility of biological morphogenesis given the stochasticity of biochemical processes and the influence of environmental perturbations. Moreover, it represents a major hurdle to synthetic engineering of Turing patterns. In this work it is shown that this problem can be overcome in some systems under certain sets of interactions between their elements, one of which is immobile and therefore corresponding to a cell-autonomous factor. In such systems Turing patterns formation can be guaranteed by a simple universal control under any values of kinetic parameters and diffusion coefficients of mobile elements. This concept is illustrated by analysis and simulations of a specific three-component system, characterized in absence of diffusion by a presence of codimension two pitchfork-Hopf bifurcation.

Beyond Turing: far-from-equilibrium patterns and mechano-chemical feedback

Turing patterns are commonly understood as specific instabilities of a spatially homogeneous steady state, resulting from activator-inhibitor interaction destabilized by diffusion. We argue that this view is restrictive and its agreement with biological observations is problematic. We present two alternatives to the classical Turing analysis of patterns. First, we employ the abstract framework of evolution equations to enable the study of far-from-equilibrium patterns. Second, we introduce a mechano-chemical model, with the surface on which the pattern forms being dynamic and playing an active role in the pattern formation, effectively replacing the inhibitor. We highlight the advantages of these two alternatives vis-à-vis the classical Turing analysis, and give an overview of recent results and future challenges for both approaches. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

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

Signaling gradients in surface dynamics as basis for planarian regeneration

Based on experimental data, we introduce and analyze a system of reaction-diffusion equations for the regeneration of planarian flatworms. We model dynamics of head and tail cells expressing positional control genes that translate into localized signals which in turn guide stem cell differentiation. Tissue orientation and positional information are encoded in a long range wnt-related signaling gradient. Our system correctly reproduces typical cut and graft experiments, and improves on previous models by preserving polarity in regeneration over orders of magnitude in body size during growth phases. Key to polarity preservation in our model flatworm is the sensitivity of cell differentiation to gradients of wnt-related signals relative to the tissue surface. This process is particularly relevant in small tissue layers close to cuts during their healing, and modeled in a robust fashion through dynamic boundary conditions.

Widening the criteria for emergence of Turing patterns

The classical concept for emergence of Turing patterns in reaction-diffusion systems requires that a system should be composed of complementary subsystems, one of which is unstable and diffuses sufficiently slowly while the other one is stable and diffuses sufficiently rapidly. In this work, the phenomena of emergence of Turing patterns are studied and do not fit into this concept, yielding the following results. (1) The criteria are derived, under which a reaction-diffusion system with immobile species should spontaneously produce Turing patterns under any diffusion coefficients of its mobile species. It is shown for such systems that under certain sets of types of interactions between their species, Turing patterns should be produced under any parameter values, at least provided that the corresponding spatially non-distributed system is stable. (2) It is demonstrated that in a reaction-diffusion system, which contains more than two species and is stable in absence of diffusion, the presence of a sufficiently slowly diffusing unstable subsystem is already sufficient for diffusion instability (i.e., Turing or wave instability), while its complementary subsystem can also be unstable. (3) It is shown that the presence of an immobile unstable subsystem, which leads to destabilization of waves within an infinite range of wavenumbers, in a spatially discrete case can result in the generation of large-scale stationary or oscillatory patterns. (4) It is demonstrated that under the presence of subcritical Turing and supercritical wave bifurcations, the interaction of two diffusion instabilities can result in the spontaneous formation of Turing structures outside the region of Turing instability.

Turing patterns in a network-reduced FitzHugh-Nagumo model

Reduction of a two-component FitzHugh-Nagumo model to a single-component model with long-range connection is considered on general networks. The reduced model describes a single chemical species reacting on the nodes and diffusing across the links with weighted long-range connections, which can be interpreted as a class of networked dynamical systems on a multigraph with local and nonlocal Laplace matrices that self-consistently emerge from the adiabatic elimination. We study the conditions for the instability of homogeneous states in the original and reduced models and show that Turing patterns can emerge in both models. We also consider generality of the adiabatic elimination for a wider class of slow-fast systems and discuss the peculiarity of the FitzHugh-Nagumo model.

Post-Turing tissue pattern formation: Advent of mechanochemistry

Chemical and mechanical pattern formation is fundamental during embryogenesis and tissue development. Yet, the underlying molecular and cellular mechanisms are still elusive in many cases. Most current theories assume that tissue development is driven by chemical processes: either as a sequence of chemical patterns each depending on the previous one, or by patterns spontaneously arising from specific chemical interactions (such as “Turing-patterns”). Within both theories, mechanical patterns are usually regarded as passive by-products of chemical pre-patters. However, several experiments question these theories, and an increasing number of studies shows that tissue mechanics can actively influence chemical patterns during development. In this study, we thus focus on the interplay between chemical and mechanical processes during tissue development. On one hand, based on recent experimental data, we develop new mechanochemical simulation models of evolving tissues, in which the full 3D representation of the tissue appears to be critical for obtaining a realistic mechanochemical behaviour. The presented modelling approach is flexible and numerically studied using state of the art finite element methods. Thus, it may serve as a basis to combine simulations with new experimental methods in tissue development. On the other hand, we apply the developed approach and demonstrate that even simple interactions between tissue mechanics and chemistry spontaneously lead to robust and complex mechanochemical patterns. Especially, we demonstrate that the main contradictions arising in the framework of purely chemical theories are naturally and automatically resolved using the mechanochemical patterning theory.

During embryogenesis, biological tissues gradually increase their complexity by self-organised creation of diverse chemical and mechanical patterns. Detailed mechanisms driving and controlling these patterns are not well understood. Previous theories mostly assume that these patterns are driven by chemical processes. Based on these theories, mechanical patterns are usually considered being mainly determined by chemical pre-patterns. However, experimental evidence for these theories is sparse, and several inconsistencies have been discovered. Furthermore, an increasing amount of data shows that tissue mechanics plays an important role in pattern formation. In this study, we present 3D computer simulations of evolving tissues to investigate the capacity of mechanochemical interactions for pattern formation. We show that even simple interactions between tissue mechanics and tissue chemistry spontaneously lead to robust chemical and mechanical pattern formation. We additionally demonstrate that main contradictions arising in the framework of purely chemical theories are naturally and automatically resolved using the mechanochemical patterning theory. The presented modelling approach can be used to combine simulations with recent experimental developments, to help unravel one of the big mysteries in biology: The mechanisms of self-organised pattern formation during embryogenesis.

Delay-induced wave instabilities in single-species reaction-diffusion systems

The Turing (wave) instability is only possible in reaction-diffusion systems with more than one (two) components. Motivated by the fact that a time delay increases the dimension of a system, we investigate the presence of diffusion-driven instabilities in single-species reaction-diffusion systems with delay. The stability of arbitrary one-component systems with a single discrete delay, with distributed delay, or with a variable delay is systematically analyzed. We show that a wave instability can appear from an equilibrium of single-species reaction-diffusion systems with fluctuating or distributed delay, which is not possible in similar systems with constant discrete delay or without delay. More precisely, we show by basic analytic arguments and by numerical simulations that fast asymmetric delay fluctuations or asymmetrically distributed delays can lead to wave instabilities in these systems. Examples, for the resulting traveling waves are shown for a Fisher-KPP equation with distributed delay in the reaction term. In addition, we have studied diffusion-induced instabilities from homogeneous periodic orbits in the same systems with variable delay, where the homogeneous periodic orbits are attracting resonant periodic solutions of the system without diffusion, i.e., periodic orbits of the Hutchinson equation with time-varying delay. If diffusion is introduced, standing waves can emerge whose temporal period is equal to the period of the variable delay.