1
|
Dúzs B, Holló G, Kitahata H, Ginder E, Suematsu NJ, Lagzi I, Szalai I. Appearance and suppression of Turing patterns under a periodically forced feed. Commun Chem 2023; 6:3. [PMID: 36697882 PMCID: PMC9814632 DOI: 10.1038/s42004-022-00800-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2022] [Accepted: 12/19/2022] [Indexed: 01/04/2023] Open
Abstract
Turing instability is a general and straightforward mechanism of pattern formation in reaction-diffusion systems, and its relevance has been demonstrated in different biological phenomena. Still, there are many open questions, especially on the robustness of the Turing mechanism. Robust patterns must survive some variation in the environmental conditions. Experiments on pattern formation using chemical systems have shown many reaction-diffusion patterns and serve as relatively simple test tools to study general aspects of these phenomena. Here, we present a study of sinusoidal variation of the input feed concentrations on chemical Turing patterns. Our experimental, numerical and theoretical analysis demonstrates that patterns may appear even at significant amplitude variation of the input feed concentrations. Furthermore, using time-dependent feeding opens a way to control pattern formation. The patterns settled at constant feed may disappear, or new patterns may appear from a homogeneous steady state due to the periodic forcing.
Collapse
Affiliation(s)
- Brigitta Dúzs
- grid.5591.80000 0001 2294 6276Laboratory of Nonlinear Chemical Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/A, H-1117 Budapest, Hungary ,grid.5802.f0000 0001 1941 7111Present Address: University of Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Gábor Holló
- ELKH-BME Condensed Matter Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - Hiroyuki Kitahata
- grid.136304.30000 0004 0370 1101Graduate School of Science, Chiba University, Yayoi-cho 1-33, Inage-ku, Chiba 263-8522 Japan
| | - Elliott Ginder
- grid.411764.10000 0001 2106 7990School of Interdisciplinary Mathematical Sciences, Graduate School of Advanced Mathematical Sciences, and Meiji Institute for Advanced Study of Mathematical Sciences (MIMS), Meiji University, 4-21-1, Nakano Tokyo, 164-8525 Japan
| | - Nobuhiko J. Suematsu
- grid.411764.10000 0001 2106 7990School of Interdisciplinary Mathematical Sciences, Graduate School of Advanced Mathematical Sciences, and Meiji Institute for Advanced Study of Mathematical Sciences (MIMS), Meiji University, 4-21-1, Nakano Tokyo, 164-8525 Japan
| | - István Lagzi
- ELKH-BME Condensed Matter Research Group, Műegyetem rkp. 3, H-1111 Budapest, Hungary ,grid.6759.d0000 0001 2180 0451Department of Physics, Institute of Physics, Budapest University of Technology and Economics, Műegyetem rkp. 3, H-1111 Budapest, Hungary
| | - István Szalai
- grid.5591.80000 0001 2294 6276Laboratory of Nonlinear Chemical Dynamics, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/A, H-1117 Budapest, Hungary
| |
Collapse
|
2
|
Wang J, Shen J. Turing instability mechanism of short-memory formation in multilayer FitzHugh-Nagumo network. Front Psychiatry 2023; 14:1083015. [PMID: 37051165 PMCID: PMC10083418 DOI: 10.3389/fpsyt.2023.1083015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Accepted: 02/14/2023] [Indexed: 04/14/2023] Open
Abstract
Introduction The study of brain function has been favored by scientists, but the mechanism of short-term memory formation has yet to be precise. Research problem Since the formation of short-term memories depends on neuronal activity, we try to explain the mechanism from the neuron level in this paper. Research contents and methods Due to the modular structures of the brain, we analyze the pattern properties of the FitzHugh-Nagumo model (FHN) on a multilayer network (coupled by a random network). The conditions of short-term memory formation in the multilayer FHN model are obtained. Then the time delay is introduced to more closely match patterns of brain activity. The properties of periodic solutions are obtained by the central manifold theorem. Conclusion When the diffusion coeffcient, noise intensity np, and network connection probability p reach a specific range, the brain forms a relatively vague memory. It is found that network and time delay can induce complex cluster dynamics. And the synchrony increases with the increase of p. That is, short-term memory becomes clearer.
Collapse
Affiliation(s)
- Junjie Wang
- School of Mathematics and Statistics, Zhengzhou University, Zhengzhou, China
| | - Jianwei Shen
- School of Mathematics and Statistics, North China University of Water Resources and Electric Power, Zhengzhou, China
- *Correspondence: Jianwei Shen
| |
Collapse
|
3
|
Van Gorder RA. Pattern formation from spatially heterogeneous reaction-diffusion systems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20210001. [PMID: 34743604 DOI: 10.1098/rsta.2021.0001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
First proposed by Turing in 1952, the eponymous Turing instability and Turing pattern remain key tools for the modern study of diffusion-driven pattern formation. In spatially homogeneous Turing systems, one or a few linear Turing modes dominate, resulting in organized patterns (peaks in one dimension; spots, stripes, labyrinths in two dimensions) which repeats in space. For a variety of reasons, there has been increasing interest in understanding irregular patterns, with spatial heterogeneity in the underlying reaction-diffusion system identified as one route to obtaining irregular patterns. We study pattern formation from reaction-diffusion systems which involve spatial heterogeneity, by way of both analytical and numerical techniques. We first extend the classical Turing instability analysis to track the evolution of linear Turing modes and the nascent pattern, resulting in a more general instability criterion which can be applied to spatially heterogeneous systems. We also calculate nonlinear mode coefficients, employing these to understand how each spatial mode influences the long-time evolution of a pattern. Unlike for the standard spatially homogeneous Turing systems, spatially heterogeneous systems may involve many Turing modes of different wavelengths interacting simultaneously, with resulting patterns exhibiting a high degree of variation over space. We provide a number of examples of spatial heterogeneity in reaction-diffusion systems, both mathematical (space-varying diffusion parameters and reaction kinetics, mixed boundary conditions, space-varying base states) and physical (curved anisotropic domains, apical growth of space domains, chemicalsimmersed within a flow or a thermal gradient), providing a qualitative understanding of how spatial heterogeneity can be used to modify classical Turing patterns. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.
Collapse
Affiliation(s)
- Robert A Van Gorder
- Department of Mathematics and Statistics, University of Otago, PO Box 56, Dunedin 9054, New Zealand
| |
Collapse
|
4
|
Krause AL, Gaffney EA, Maini PK, Klika V. Modern perspectives on near-equilibrium analysis of Turing systems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2021; 379:20200268. [PMID: 34743603 PMCID: PMC8580451 DOI: 10.1098/rsta.2020.0268] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 06/18/2021] [Indexed: 05/02/2023]
Abstract
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'.
Collapse
Affiliation(s)
- Andrew L. Krause
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
- Department of Mathematical Sciences, Durham University, Upper Mountjoy Campus, Stockton Rd, Durham DH1 3LE, UK
| | - Eamonn A. Gaffney
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Philip K. Maini
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford OX2 6GG, UK
| | - Václav Klika
- Department of Mathematics, FNSPE, Czech Technical University in Prague, Trojanova, 13, 12000 Praha, Czech Republic
| |
Collapse
|
5
|
Van Gorder RA. A theory of pattern formation for reaction–diffusion systems on temporal networks. Proc Math Phys Eng Sci 2021. [DOI: 10.1098/rspa.2020.0753] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Networks have become ubiquitous in the modern scientific literature, with recent work directed at understanding ‘temporal networks’—those networks having structure or topology which evolves over time. One area of active interest is pattern formation from reaction–diffusion systems, which themselves evolve over temporal networks. We derive analytical conditions for the onset of diffusive spatial and spatio-temporal pattern formation on undirected temporal networks through the Turing and Benjamin–Feir mechanisms, with the resulting pattern selection process depending strongly on the evolution of both global diffusion rates and the local structure of the underlying network. Both instability criteria are then extended to the case where the reaction–diffusion system is non-autonomous, which allows us to study pattern formation from time-varying base states. The theory we present is illustrated through a variety of numerical simulations which highlight the role of the time evolution of network topology, diffusion mechanisms and non-autonomous reaction kinetics on pattern formation or suppression. A fundamental finding is that Turing and Benjamin–Feir instabilities are generically transient rather than eternal, with dynamics on temporal networks able to transition between distinct patterns or spatio-temporal states. One may exploit this feature to generate new patterns, or even suppress undesirable patterns, over a given time interval.
Collapse
Affiliation(s)
- Robert A. Van Gorder
- Department of Mathematics and Statistics, University of Otago, PO Box 56, Dunedin 9054, New Zealand
| |
Collapse
|
6
|
Abstract
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.
Collapse
Affiliation(s)
- Thomas E Woolley
- Cardiff School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK.
| | - Andrew L Krause
- Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK
| | - Eamonn A Gaffney
- Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK
| |
Collapse
|
7
|
Van Gorder RA, Klika V, Krause AL. Turing conditions for pattern forming systems on evolving manifolds. J Math Biol 2021; 82:4. [PMID: 33475826 DOI: 10.1007/s00285-021-01552-y] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2019] [Revised: 07/29/2020] [Accepted: 10/16/2020] [Indexed: 12/18/2022]
Abstract
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.
Collapse
Affiliation(s)
- Robert A Van Gorder
- Department of Mathematics and Statistics, University of Otago, P.O. Box 56, Dunedin, 9054, New Zealand.
| | - Václav Klika
- Department of Mathematics FNSPE, Czech Technical University in Prague, Trojanova 13, 12000, Prague, Czech Republic
| | - Andrew L Krause
- Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG, UK
| |
Collapse
|
8
|
Abad E, Angstmann CN, Henry BI, McGann AV, Le Vot F, Yuste SB. Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains. Phys Rev E 2020; 102:032111. [PMID: 33075977 DOI: 10.1103/physreve.102.032111] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2020] [Accepted: 08/19/2020] [Indexed: 01/22/2023]
Abstract
Reaction-diffusion equations are widely used as the governing evolution equations for modeling many physical, chemical, and biological processes. Here we derive reaction-diffusion equations to model transport with reactions on a one-dimensional domain that is evolving. The model equations, which have been derived from generalized continuous time random walks, can incorporate complexities such as subdiffusive transport and inhomogeneous domain stretching and shrinking. Inhomogeneously growing domains are frequently encountered in biological phenomena involving stochastic transport, such as tumor growth and morphogen gradient formation. A method for constructing analytic expressions for short-time moments of the position of the particles is developed and moments calculated from this approach are shown to compare favorably with results from random walk simulations and numerical integration of the reaction transport equation. The results show the important role played by the initial condition. In particular, it strongly affects the time dependence of the moments in the short-time regime by introducing additional drift and diffusion terms. We also discuss how our reaction transport equation could be applied to study the spreading of a population on an evolving interface. From a more general perspective, our findings help to mitigate the scarcity of analytic results for reaction-diffusion problems in geometries displaying nonuniform growth. They are also expected to pave the way for further results, including the treatment of first-passage problems associated with encounter-controlled reactions in such domains.
Collapse
Affiliation(s)
- E Abad
- Departamento de Física Aplicada and Instituto de Computación Científica Avanzada, Centro Universitario de Mérida, Universidad de Extremadura, 06800 Mérida, Spain
| | - C N Angstmann
- School of Mathematics and Statistics, UNSW, Sydney New South Wales, 2052, Australia
| | | | | | - F Le Vot
- Departamento de Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain
| | - S B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain
| |
Collapse
|
9
|
Ledesma-Durán A, Ortiz-Durán EA, Aragón JL, Santamaría-Holek I. Eckhaus selection: The mechanism of pattern persistence in a reaction-diffusion system. Phys Rev E 2020; 102:032214. [PMID: 33076036 DOI: 10.1103/physreve.102.032214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2020] [Accepted: 09/01/2020] [Indexed: 11/07/2022]
Abstract
In this work, we show theoretically and numerically that a one-dimensional reaction-diffusion system, near the Turing bifurcation, produces different number of stripes when, in addition to random noise, the Fourier mode of a prepattern used to initialize the system changes. We also show that the Fourier modes that persist are inside the Eckhaus stability regions, while those outside this region follow a wave number selection process not predicted by the linear analysis. To test our results, we use the Brusselator reaction-diffusion system obtaining an excellent agreement between the weakly nonlinear predictions of the real Ginzburg-Landau equations and the numerical solutions near the bifurcation. Although the persistence of patterns is not relevant as a simple generating mechanism of self-organization, it is crucial to understand the formation of patterns that occurs in multiple stages. In this work, we discuss the relevance of our results on the robustness and diversity of solutions in multiple-steps mechanisms of biological pattern formation and auto-organization in growing domains.
Collapse
Affiliation(s)
- Aldo Ledesma-Durán
- Centro de Física Aplicada y Tecnología Avanzada, Universidad Nacional Autónoma de México, Boulevard Juriquilla 3001, Juriquilla, 76230 Querétaro, Mexico
| | - E A Ortiz-Durán
- Centro de Física Aplicada y Tecnología Avanzada, Universidad Nacional Autónoma de México, Boulevard Juriquilla 3001, Juriquilla, 76230 Querétaro, Mexico
| | - J L Aragón
- Centro de Física Aplicada y Tecnología Avanzada, Universidad Nacional Autónoma de México, Boulevard Juriquilla 3001, Juriquilla, 76230 Querétaro, Mexico
| | - Iván Santamaría-Holek
- Unidad Multidiciplinaria de Docencia e Investigación, Facutad de Ciencias, Universidad Nacional Autónoma de México, Boulevard Juriquilla 3001, Juriquilla, 76230 Querétaro, Mexico
| |
Collapse
|
10
|
Abstract
The Turing instability is one of the most commonly studied mechanisms leading to pattern formation in reaction–diffusion systems, yet there are still many open questions on the applicability of the Turing mechanism. Although experiments on pattern formation using chemical systems have shown that temperature differences play a role in pattern formation, there is far less theoretical work concerning the interplay between temperature and spatial instabilities. We consider a thermodynamically extended reaction–diffusion system, consisting of a pair of reaction–diffusion equations coupled to an energy equation for temperature, and use this to obtain a natural extension of the Turing instability accounting for temperature. We show that thermal contributions can restrict or enlarge the set of unstable modes possible under the instability, and in some cases may be used to completely shift the set of unstable modes, strongly modifying emergent Turing patterns. Spatial heterogeneity plays a role under several experimentally feasible configurations, and we give particular consideration to scenarios involving thermal gradients, thermodynamics of chemicals transported within a flow, and thermodiffusion. Control of Turing patterns is also an area of active interest, and we also demonstrate how patterns can be modified using time-dependent control of the boundary temperature.
Collapse
Affiliation(s)
- Robert A. Van Gorder
- Department of Mathematics and Statistics, University of Otago, PO Box 56, Dunedin 9054, New Zealand
| |
Collapse
|