Turing Instabilities are Not Enough to Ensure Pattern Formation

Andrew L. Krause,Eamonn A. Gaffney,Thomas Jun Jewell,Václav Klika,Benjamin J. Walker
DOI: https://doi.org/10.1007/s11538-023-01250-4
2024-01-24
Bulletin of Mathematical Biology
Abstract:Symmetry-breaking instabilities play an important role in understanding the mechanisms underlying the diversity of patterns observed in nature, such as in Turing's reaction–diffusion theory, which connects cellular signalling and transport with the development of growth and form. Extensive literature focuses on the linear stability analysis of homogeneous equilibria in these systems, culminating in a set of conditions for transport-driven instabilities that are commonly presumed to initiate self-organisation. We demonstrate that a selection of simple, canonical transport models with only mild multistable non-linearities can satisfy the Turing instability conditions while also robustly exhibiting only transient patterns. Hence, a Turing-like instability is insufficient for the existence of a patterned state. While it is known that linear theory can fail to predict the formation of patterns, we demonstrate that such failures can appear robustly in systems with multiple stable homogeneous equilibria. Given that biological systems such as gene regulatory networks and spatially distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this raises important questions of how to analyse prospective mechanisms for self-organisation.
mathematical & computational biology,biology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is whether **Turing instability is sufficient to ensure pattern formation**. Specifically, the authors explored that in some simple transport models, even when the conditions for Turing instability are met, the system may only exhibit transient patterns rather than a persistent pattern state. Therefore, the core problem of the paper is to challenge the traditional linear stability analysis method, pointing out that Turing instability itself is not sufficient to guarantee the long - term pattern formation of the system. ### Background Turing instability is one of the key theories for understanding various pattern - formation mechanisms in nature, especially in reaction - diffusion systems. These systems explain the development of growth and morphology through cell signaling and mass transport. However, although a large amount of literature focuses on the linear stability analysis of these systems, these analyses usually can only predict the formation of initial patterns, but cannot guarantee the long - term stable existence of patterns. ### Main Findings 1. **Transient Patterns**: The authors showed several simple transport models that, even when the conditions for Turing instability are met, will only exhibit transient patterns rather than a persistent pattern state. 2. **The Influence of Multistability**: The existence of multistability can disrupt typical pattern predictions, causing the system to return to a stable, homogeneous equilibrium state after the appearance of transient patterns. 3. **The Complexity of Biological Systems**: Since biological systems such as gene regulatory networks and spatially - distributed ecosystems often exhibit a high degree of multistability and nonlinearity, this finding raises important questions for analyzing self - organizing mechanisms in these systems. ### Conclusion The paper emphasizes that although Turing instability is a necessary condition for pattern formation, it is not a sufficient condition. In order to more accurately understand the mechanisms of pattern formation, new tools and methods need to be developed, going beyond traditional linear stability analysis, especially when dealing with complex multi - component systems. ### Formula Examples Some of the key formulas involved in the paper are as follows: 1. **The Linearized Equation of the Reaction - Diffusion System**: \[ \lambda_k=\frac{1 - acd-\rho_k(1 + D)+\sqrt{(1 - acd-\rho_k(1 + D))^2-4(\rho_k^2D-\rho_k(D - acd)-acd + b)}}{2} \] 2. **The Linearized Equation of the Keller - Segel Model**: \[ \lambda_k=\frac{T_{\text{KS}}+\sqrt{T_{\text{KS}}^2-4Q_{\text{KS}}}}{2} \] where \(T_{\text{KS}}=-b(b - d)-a-\rho_k(1 + D)\), \(Q_{\text{KS}}=\rho_k^2D-\rho_k(cb - a - Db(b - d))+b(b - d)a\). 3. **The Linearized Equation of the Fourth - Order Model**: \[ \lambda_k = D\rho_k-\rho_k^2+ac(b - c) \] 4. **The Linearized Equation of the Non - local Convection Model**: \[ \lambda_k=-ac(c - b)-D\rho_k+\frac{c(1 - c)d\rho_k}{\pi^{2 - N}(1 + \rho_k)^{N + 1}} \] These formulas are used to describe the growth rate of perturbations in different models, thereby determining whether the system meets the conditions for Turing instability.