Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Jun 26;115(26):6572-6577.
doi: 10.1073/pnas.1720770115. Epub 2018 Jun 11.

Stochastic Turing patterns in a synthetic bacterial population

Affiliations

Stochastic Turing patterns in a synthetic bacterial population

David Karig et al. Proc Natl Acad Sci U S A. .

Abstract

The origin of biological morphology and form is one of the deepest problems in science, underlying our understanding of development and the functioning of living systems. In 1952, Alan Turing showed that chemical morphogenesis could arise from a linear instability of a spatially uniform state, giving rise to periodic pattern formation in reaction-diffusion systems but only those with a rapidly diffusing inhibitor and a slowly diffusing activator. These conditions are disappointingly hard to achieve in nature, and the role of Turing instabilities in biological pattern formation has been called into question. Recently, the theory was extended to include noisy activator-inhibitor birth and death processes. Surprisingly, this stochastic Turing theory predicts the existence of patterns over a wide range of parameters, in particular with no severe requirement on the ratio of activator-inhibitor diffusion coefficients. To explore whether this mechanism is viable in practice, we have genetically engineered a synthetic bacterial population in which the signaling molecules form a stochastic activator-inhibitor system. The synthetic pattern-forming gene circuit destabilizes an initially homogenous lawn of genetically engineered bacteria, producing disordered patterns with tunable features on a spatial scale much larger than that of a single cell. Spatial correlations of the experimental patterns agree quantitatively with the signature predicted by theory. These results show that Turing-type pattern-forming mechanisms, if driven by stochasticity, can potentially underlie a broad range of biological patterns. These findings provide the groundwork for a unified picture of biological morphogenesis, arising from a combination of stochastic gene expression and dynamical instabilities.

Keywords: Turing patterns; biofilm; signaling molecules; stochastic gene expression; synthetic biology.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Fig. 1.
Fig. 1.
Design of a synthetic multicellular system for emergent pattern formation. (A) Abstractly, the system consists of two signaling species A3OC12HSL and IC4HSL. A3OC12HSL is an activator catalyzing synthesis of both species, while IC4HSL is an inhibitor repressing their synthesis, with additional repression by A3OC12HSL via competitive binding. (B) Genetic circuit implementation. Promoter regions are indicated by white boxes, while protein coding sequences are indicated by colored boxes. IPTG is an external inducer modulating system dynamics. (C, Top) Illustration of signaling species concentrations in 1D space. The dashed orange and blue lines correspond to A3OC12HSL and IC4HSL, respectively. (C, Middle) Spatial profiles of reporter proteins. RFP expression (red line) correlates with A3OC12HSL concentrations, while GFP expression (green line) roughly mirrors RFP expression. (C, Bottom) A vertical slice of cell lawn. Cells express fluorescence proteins according to the profiles above and produce a global multicellular pattern.
Fig. 2.
Fig. 2.
Experimental observations of emergent pattern formation measured by green fluorescence intensity (GFI) and red fluorescence intensity (RFI). (A) Representative microscope images (based on six technical replicates) of a typical field of view showing a fluorescent pattern formed by an initially homogeneous isogenic lawn of cells harboring the Turing circuit with no IPTG. Spots and voids appear in the red fluorescence (RF) and green fluorescence (GF) channels, respectively. (Scale bar: 100 μm.) (B) Microscope images of cell lawns with constitutive expression of fluorescent proteins. (Left) Cells expressing RFP, (Center) cells expressing GFP, and (Right) mixed population of red and green cells. (C) Fluorescence density plots computed from the images above (from left to right: red, green, red/green, and Turing). Color intensity is in log scale [arbitrary units (a.u.)].
Fig. 3.
Fig. 3.
Mathematical modeling and correlation between pattern modulation experiments and simulations. (A) Experimental results for IPTG modulation of pattern formation with microscopy images corresponding to specific IPTG concentrations in B. The same display mappings were used for all images in A. (B) Collectivity, metric parameter Θ is influenced by IPTG modulation. (C) Pattern statistics over IPTG modulation for experimental results. (D) Pattern obtained from simulating a deterministic reaction–diffusion model with Dv/Du=100. (E) Pattern statistics over IPTG modulation for deterministic modeling. (F) Patterns obtained from simulating our deterministic model (Upper) and stochastic spatiotemporal model (Lower) at the measured diffusion ratio of Dv/Du=21.6. (G) Pattern statistics over IPTG modulation for stochastic modeling. a.u., arbitrary unit.
Fig. 4.
Fig. 4.
Spectral analysis and parameter analysis. (A) Pattern-forming regimes in parameter space and estimated parameters for our system. Parameters above the green surface of neutral stochastic stability can form stochastic patterns, and parameters above the blue surface of deterministic neutral stability can form deterministic Turing patterns. The ratio of the diffusion coefficients ν/μ, the ratio of degradation rate to production rate d/p, and the ratio of production rates are estimated for our system by the yellow ellipsoid. The parameters for our system are mostly in the regime where stochastic patterns form and outside the region where deterministic Turing patterns form. Example stochastic simulations are shown for parameters drawn from a deterministic parameter region with Dν/Dμ=100 (Upper Right) and a stochastic region with Dν/Dμ=21.6 (Lower Right). (B) Radial power spectrum of green fluorescence and best fit power law tail with an exponent of 2.3±0.2. (C) Radial power spectrum for eight trials of our stochastic simulation, their mean, and the best fit power law tail.

References

    1. Gierer A, Meinhardt H. A theory of biological pattern formation. Biol Cybern. 1972;12:30–39. - PubMed
    1. Turing A. The chemical basis of morphogenesis. Philos Trans R Soc Lond Ser B Biol Sci. 1952;237:37–72. - PMC - PubMed
    1. Castets V, Dulos E, Boissonade J, De Kepper P. Experimental evidence of a sustained standing Turing-type nonequilibrium chemical pattern. Phys Rev Lett. 1990;64:2953–2956. - PubMed
    1. Raspopovic J, Marcon L, Russo L, Sharpe J. Digit patterning is controlled by a Bmp-Sox9-Wnt Turing network modulated by morphogen gradients. Science. 2014;345:566–570. - PubMed
    1. Butler T, Goldenfeld N. Robust ecological pattern formation induced by demographic noise. Phys Rev E. 2009;80:030902. - PubMed

Publication types

MeSH terms

LinkOut - more resources