The influence of gene expression time delays on Gierer-Meinhardt pattern formation systems

There are numerous examples of morphogen gradients controlling long range signalling in developmental and cellular systems. The prospect of two such interacting morphogens instigating long range self-organisation in biological systems via a Turing bifurcation has been explored, postulated, or implicated in the context of numerous developmental processes. However, modelling investigations of cellular systems typically neglect the influence of gene expression on such dynamics, even though transcription and translation are observed to be important in morphogenetic systems. In particular, the influence of gene expression on a large class of Turing bifurcation models, namely those with pure kinetics such as the Gierer–Meinhardt system, is unexplored. Our investigations demonstrate that the behaviour of the Gierer–Meinhardt model profoundly changes on the inclusion of gene expression dynamics and is sensitive to the sub-cellular details of gene expression. Features such as concentration blow up, morphogen oscillations and radical sensitivities to the duration of gene expression are observed and, at best, severely restrict the possible parameter spaces for feasible biological behaviour. These results also indicate that the behaviour of Turing pattern formation systems on the inclusion of gene expression time delays may provide a means of distinguishing between possible forms of interaction kinetics. Finally, this study also emphasises that sub-cellular and gene expression dynamics should not be simply neglected in models of long range biological pattern formation via morphogens

Pattern formation and nonlocal logistic growth

Logistic growth process with nonlocal interactions is considered in one dimension. Spontaneous breakdown of translational invariance is shown to take place at some parameter region, and the bifurcation regime is identified for short and long range interactions. Domain walls between regions of different order parameter are expressed as soliton solutions of the reduced dynamics for nearest neighbor interactions. The analytic results are confirmed by numerical simulations

Patchiness and Demographic Noise in Three Ecological Examples

Understanding the causes and effects of spatial aggregation is one of the most fundamental problems in ecology. Aggregation is an emergent phenomenon arising from the interactions between the individuals of the population, able to sense only -at most- local densities of their cohorts. Thus, taking into account the individual-level interactions and fluctuations is essential to reach a correct description of the population. Classic deterministic equations are suitable to describe some aspects of the population, but leave out features related to the stochasticity inherent to the discreteness of the individuals. Stochastic equations for the population do account for these fluctuation-generated effects by means of demographic noise terms but, owing to their complexity, they can be difficult (or, at times, impossible) to deal with. Even when they can be written in a simple form, they are still difficult to numerically integrate due to the presence of the "square-root" intrinsic noise. In this paper, we discuss a simple way to add the effect of demographic stochasticity to three classic, deterministic ecological examples where aggregation plays an important role. We study the resulting equations using a recently-introduced integration scheme especially devised to integrate numerically stochastic equations with demographic noise. Aimed at scrutinizing the ability of these stochastic examples to show aggregation, we find that the three systems not only show patchy configurations, but also undergo a phase transition belonging to the directed percolation universality class.Comment: 20 pages, 5 figures. To appear in J. Stat. Phy

Turing instabilities in general systems

We present necessary and sufficient conditions on the stability matrix of a general n(S2)-dimensional reaction-diffusion system which guarantee that its uniform steady state can undergo a Turing bifurcation. The necessary (kinetic) condition, requiring that the system be composed of an unstable (or activator) and a stable (or inhibitor) subsystem, and the sufficient condition of sufficiently rapid inhibitor diffusion relative to the activator subsystem are established in three theorems which form the core of our results. Given the possibility that the unstable (activator) subsystem involves several species (dimensions), we present a classification of the analytically deduced Turing bifurcations into p (1 h p h (n m 1)) different classes. For n = 3 dimensions we illustrate numerically that two types of steady Turing pattern arise in one spatial dimension in a generic reaction-diffusion system. The results confirm the validity of an earlier conjecture [12] and they also characterise the class of so-called strongly stable matrices for which only necessary conditions have been known before [23, 24]. One of the main consequences of the present work is that biological morphogens, which have so far been expected to be single chemical species [1-9], may instead be composed of two or more interacting species forming an unstable subsystem

Multispecies vegetation patterns in arid to semi-arid lands

I vegetation patterns nei climi aridi e semi-aridi sono un fenomeno a tutt'oggi studiato. Benché una zonazione delle bande dal punto di vista della composizione botanica sia stata ampiamente documentata, il ruolo delle diverse specie vegetali non è stato ancora trattato formalmente nei modelli matematici che descrivono la formazione e la dinamica dei tiger bush. Lo scopo di questo lavoro è di indagare se meccanismi di competizione sono in grado di spiegare la zonazione delle band

Interacting Turing-Hopf Instabilities Drive Symmetry-Breaking Transitions in a Mean-Field Model of the Cortex: A Mechanism for the Slow Oscillation

Electrical recordings of brain activity during the transition from wake to anesthetic coma show temporal and spectral alterations that are correlated with gross changes in the underlying brain state. Entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (≲1  Hz) similar to the slow waves observed in natural sleep. Here we present a two-dimensional mean-field model of the cortex in which slow spatiotemporal oscillations arise spontaneously through a Turing (spatial) symmetry-breaking bifurcation that is modulated by a Hopf (temporal) instability. In our model, populations of neurons are densely interlinked by chemical synapses, and by interneuronal gap junctions represented as an inhibitory diffusive coupling. To demonstrate cortical behavior over a wide range of distinct brain states, we explore model dynamics in the vicinity of a general-anesthetic-induced transition from “wake” to “coma.” In this region, the system is poised at a codimension-2 point where competing Turing and Hopf instabilities coexist. We model anesthesia as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic response, producing a coma state that is characterized by emergent low-frequency oscillations whose dynamics is chaotic in time and space. The effect of long-range axonal white-matter connectivity is probed with the inclusion of a single idealized point-to-point connection. We find that the additional excitation from the long-range connection can provoke seizurelike bursts of cortical activity when inhibitory diffusion is weak, but has little impact on an active cortex. Our proposed dynamic mechanism for the origin of anesthetic slow waves complements—and contrasts with—conventional explanations that require cyclic modulation of ion-channel conductances. We postulate that a similar bifurcation mechanism might underpin the slow waves of natural sleep and comment on the possible consequences of chaotic dynamics for memory processing and learning

Excitable FitzHugh-Nagumo model with cross-diffusion: long-range activation instabilities

In this paper, we shall study a spatially extended version of the FitzHugh-Nagumo model, where one describes the motion of the species through cross-diffusion. The motivation comes from modeling biological species where reciprocal interaction influences spatial movement. We shall focus our analysis on the excitable regime of the system. In this case, we shall see how cross-diffusion terms can destabilize uniform equilibrium, allowing for the formation of close-to-equilibrium patterns; the species are out-of-phase spatially distributed, namely high concentration areas of one species correspond to a low density of the other (cross-Turing patterns). Moreover, depending on the magnitude of the inhibitor’s cross-diffusion, the pattern’s development can proceed in either case of the inhibitor/activator diffusivity ratio being higher or smaller than unity. This allows for spatial segregation of the species in both cases of short-range activation/long-range inhibition or long-range activation/short-range inhibition

Principles and theory of protein-based pattern formation

Biological systems perform functions by the orchestrated interplay of many small components without a "conductor." Such self-organization pervades life on many scales, from the subcellular level to populations of many organisms and whole ecosystems. On the intracellular level, protein-based pattern formation coordinates and instructs functions like cell division, differentiation and motility. A key feature of protein-based pattern formation is that the total numbers of the involved proteins remain constant on the timescale of pattern formation. The overarching theme of this thesis is the profound impact of this mass-conservation property on pattern formation and how one can harness mass conservation to understand the underlying physical principles. The central insight is that changes in local densities shift local reactive equilibria, and thus induce concentration gradients which, in turn, drive diffusive transport of mass. For two-component systems, this dynamic interplay can be captured by simple geometric objects in the (low-dimensional) phase space of chemical concentrations. On this phase-space level, physical insight can be gained from geometric criteria and graphical constructions. Moreover, we introduce the notion of regional (in)stabilities, which allows one to characterize the dynamics in the highly nonlinear regime reveals an inherent connection between Turing instability and stimulus-induced pattern formation. The insights gained for conceptual two-component systems can be generalized to systems with more components and several conserved masses. In the minimal setting of two diffusively coupled "reactors," the full dynamics can be embedded in the phase-space of redistributed masses where the phase space flow is organized by surfaces of local reactive equilibria. Building on the phase-space analysis for two component systems, we develop a new approach to the important open problem of wavelength selection in the highly nonlinear regime. We show that two-component reaction–diffusion systems always exhibit uninterrupted coarsening (the continual growth of the characteristic length scale) of patterns if they are strictly mass conserving. Selection of a finite wavelength emerges due to weakly broken mass-conservation, or coupling to additional components, which counteract and stop the competition instability that drives coarsening. For complex dynamical phenomena like wave patterns and the transition to spatiotemporal chaos, an analysis in terms of local equilibria and their stability properties provides a powerful tool to interpret data from numerical simulations and experiments, and to reveal the underlying physical mechanisms. In collaborations with different experimental labs, we studied the Min system of Escherichia coli. A central insight from these investigations is that bulk-surface coupling imparts a strong dependence of pattern formation on the geometry of the spatial confinement, which explains the qualitatively different dynamics observed inside cells compared to in vitro reconstitutions. By theoretically studying the polarization machinery in budding yeast and testing predictions in collaboration with experimentalists, we found that this functional module implements several redundant polarization mechanisms that depend on different subsets of proteins. Taken together, our work reveals unifying principles underlying biological self-organization and elucidates how microscopic interaction rules and physical constraints collectively lead to specific biological functions.Biologische Systeme führen Funktionen durch das orchestrierte Zusammenspiel vieler kleiner Komponenten ohne einen "Dirigenten" aus. Solche Selbstorganisation durchdringt das Leben auf vielen Skalen, von der subzellulären Ebene bis zu Populationen vieler Organismen und ganzen Ökosystemen. Auf der intrazellulären Ebene koordiniert und instruieren proteinbasierte Muster Funktionen wie Zellteilung, Differenzierung und Motilität. Ein wesentliches Merkmal der proteinbasierten Musterbildung ist, dass die Gesamtzahl der beteiligten Proteine auf der Zeitskala der Musterbildung konstant bleibt. Das übergreifende Thema dieser Arbeit ist es, den tiefgreifenden Einfluss dieser Massenerhaltung auf die Musterbildung zu untersuchen und Methoden zu entwickeln, die Massenerhaltung nutzen, um die zugrunde liegenden physikalischen Prinzipien von proteinbasierter Musterbildung zu verstehen. Die zentrale Erkenntnis ist, dass Änderungen der lokalen Dichten lokale reaktive Gleichgewichte verschieben und somit Konzentrationsgradienten induzieren, die wiederum den diffusiven Transport von Masse antreiben. Für Zweikomponentensysteme kann dieses dynamische Wechselspiel durch einfache geometrische Objekte im (niedrigdimensionalen) Phasenraum der chemischen Konzentrationen erfasst werden. Auf dieser Phasenraumebene können physikalische Erkenntnisse durch geometrische Kriterien und grafische Konstruktionen gewonnen werden. Darüber hinaus führen wir den Begriff der regionalen (In-)stabilität ein, der es erlaubt, die Dynamik im hochgradig nichtlinearen Regime zu charakterisieren und einen inhärenten Zusammenhang zwischen Turing-Instabilität und stimulusinduzierter Musterbildung aufzuzeigen. Die für konzeptionelle Zweikomponentensysteme gewonnenen Erkenntnisse können auf Systeme mit mehr Komponenten und mehreren erhaltenen Massen verallgemeinert werden. In der minimalen Fassung von zwei diffusiv gekoppelten "Reaktoren" kann die gesamte Dynamik in den Phasenraum umverteilter Massen eingebettet werden, wobei der Phasenraumfluss durch Flächen lokaler reaktiver Gleichgewichte organisiert wird. Aufbauend auf der Phasenraumanalyse für Zweikomponentensysteme entwickeln wir einen neuen Ansatz für die wichtige offene Fragestellung der Wellenängenselektion im hochgradig nichtlinearen Regime. Wir zeigen, dass "coarsening" (das stetige wachsen der charakteristischen Längenskala) von Mustern in Zweikomponentensystemen nie stoppt, wenn sie exakt massenerhaltend sind. Die Selektion einer endlichen Wellenlänge entsteht durch schwach gebrochene Massenerhaltung oder durch Kopplung an zusätzliche Komponenten. Diese Prozesse wirken der Masseumverteilung, die coarsening treibt, entgegen und stoppen so das coarsening. Bei komplexen dynamischen Phänomenen wie Wellenmustern und dem Übergang zu raumzeitlichen Chaos bietet eine Analyse in Bezug auf lokale Gleichgewichte und deren Stabilitätseigenschaften ein leistungsstarkes Werkzeug, um Daten aus numerischen Simulationen und Experimenten zu interpretieren und die zugrunde liegenden physikalischen Mechanismen aufzudecken. In Zusammenarbeit mit verschiedenen experimentellen Labors haben wir das Min-System von Escherichia coli untersucht. Eine zentrale Erkenntnis aus diesen Untersuchungen ist, dass die Kopplung zwischen Volumen und Oberfläche zu einer starken Abhängigkeit der Musterbildung von der räumlichen Geometrie führt. Das erklärt die qualitativ unterschiedliche Dynamik, die in Zellen im Vergleich zu in vitro Rekonstitutionen beobachtet wird. Durch die theoretische Untersuchung der Polarisationsmaschinerie in Hefezellen, kombiniert mit experimentellen Tests theoretischer Vorhersagen, haben wir herausgefunden, dass dieses Funktionsmodul mehrere redundante Polarisationsmechanismen implementiert, die von verschiedenen Untergruppen von Proteinen abhängen. Zusammengenommen beleuchtet unsere Arbeit die vereinheitlichenden Prinzipien, die der intrazellulären Selbstorganisation zugrunde liegen, und zeigt, wie mikroskopische Interaktionsregeln und physikalische Bedingungen gemeinsam zu spezifischen biologischen Funktionen führen