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
. 2010 Jun 29;107(26):11715-20.
doi: 10.1073/pnas.1001994107. Epub 2010 May 24.

Arrested phase separation in reproducing bacteria creates a generic route to pattern formation

Affiliations

Arrested phase separation in reproducing bacteria creates a generic route to pattern formation

M E Cates et al. Proc Natl Acad Sci U S A. .

Abstract

We present a generic mechanism by which reproducing microorganisms, with a diffusivity that depends on the local population density, can form stable patterns. For instance, it is known that a decrease of bacterial motility with density can promote separation into bulk phases of two coexisting densities; this is opposed by the logistic law for birth and death that allows only a single uniform density to be stable. The result of this contest is an arrested nonequilibrium phase separation in which dense droplets or rings become separated by less dense regions, with a characteristic steady-state length scale. Cell division predominates in the dilute regions and cell death in the dense ones, with a continuous flux between these sustained by the diffusivity gradient. We formulate a mathematical model of this in a case involving run-and-tumble bacteria and make connections with a wider class of mechanisms for density-dependent motility. No chemotaxis is assumed in the model, yet it predicts the formation of patterns strikingly similar to some of those believed to result from chemotactic behavior.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Fig. 1.
Fig. 1.
Growth of the instability in the supercritical (Left) and subcritical cases (Right). The three lines correspond to three successive times. A small perturbation around ρ0 (red line) growth toward harmonic or anharmonic patterns in the supercritical or subcritical case, respectively. (Left) Supercritical case (α = κ = 0.01, λ = 0.02, ρ0 = 15, formula image; times: 102, 103, 104). (Right) Subcritical case (α = κ = 0.005, λ = 0.02, ρ0 = 11, formula image; times: 3.102, 3.103, 105).
Fig. 2.
Fig. 2.
Three plots of Λq(q) for formula image, 2, 3 (from bottom to top). At the transition, only one critical mode q = qc is unstable.
Fig. 3.
Fig. 3.
(Top) Phase diagram in the (R,Φ) plane. The outer region corresponds to stable behavior, whereas within the curve, patterning occurs. The solid line is the theoretical phase boundary—Eq. 7—which accurately fits the numerics (black squares). The blue and red sections correspond to continuous and discontinuous transitions, respectively. The two magenta dots correspond to two 2D simulations that show ordered harmonic patterns close to supercriticality and amorphous patterns otherwise. (Bottom Left) Transition in the supercritical regime. The blues lines correspond to the theory—Eq. 8—whereas the squares come from simulations (Φ = 1.5, 1.35, 1.2 from top to bottom). (Bottom Right) Transition in the subcritical regime for Φ = 1.06 and Φ = 1.7 (bottom to top).
Fig. 4.
Fig. 4.
Numerical results for a 2D simulation with size equal to (in dimensionless units, see Eq. 5) 28 × 28, R ≃ 316, and Φ = 1.5. Times corresponding to the snapshots are (from left to right) formula image, 0.12, 0.20, and 0.51. The color bar shows values of the dimensionless density u; see Eq. 5.
Fig. 5.
Fig. 5.
Plot of the characteristic domain size, formula image, as a function of time, formula image (both in dimensionless units), for a system in the inhomogeneous phase, with initially random density fluctuations around ρ0. Parameters were R ∼ 316 and Φ = 1.35, whereas the system size was 35.5 × 35.5 (in dimensionless units). The solid line corresponds to a single run, whereas the dashed line is an average over six runs. The steps in the single run curve correspond to evaporation-condensation events, highlighted by black squares in the snapshots shown in the figure (before and after one of the steps, respectively, arrows indicate positions on the plot corresponding to the two snapshots). The color bar shows values of the dimensionless density u.
Fig. 6.
Fig. 6.
Dynamics of formation of patterns in 2D, starting from a single small bacterial droplet in the middle of the simulation sample. (Top) Formation of rings in a system with R = 100 and Φ = 1.65. The simulation box has size 125 × 125 (in dimensionless units). The snapshots correspond to times equal to (from left to right) formula image, 5, 10, and 27. (Bottom) Breakage of rings into dots. The four snapshots correspond to the time evolution of a system with R = 100 and Φ = 1.3. We show a 125 × 125 fraction of the simulation box, with the boundaries far away and not affecting the pattern. The snapshots correspond to times equal to (from left to right) formula image, 7, 29, and 122. For both rows, the color bar shows values of the dimensionless density u.

Comment in

  • Chemotactic patterns without chemotaxis.
    Brenner MP. Brenner MP. Proc Natl Acad Sci U S A. 2010 Jun 29;107(26):11653-4. doi: 10.1073/pnas.1006363107. Epub 2010 Jun 18. Proc Natl Acad Sci U S A. 2010. PMID: 20562344 Free PMC article. No abstract available.

References

    1. Shapiro JA. The significance of bacterial colony patterns. BioEssays. 1995;17:597–607. - PubMed
    1. Harshey RM. Bacterial motility on a surface: Many ways to a common goal. Annu Rev Microbiol. 2003;57:249–273. - PubMed
    1. Budrene EO, Berg HC. Complex patterns formed by motile cells of Escherichia coli. Nature. 1991;349:630–633. - PubMed
    1. Budrene EO, Berg HC. Dynamics of formation of symmetrical patterns by chemotactic bacteria. Nature. 1995;376:49–53. - PubMed
    1. Woodward DE, et al. Spatio-temporal patterns generated by Salmonella typhimurium. Biophys J. 1995;68:2181–2189. - PMC - PubMed

Publication types

MeSH terms

LinkOut - more resources