The idiosyncrasy of spatial structure in bacterial competition
© Hol et al. 2015
Received: 22 April 2014
Accepted: 13 May 2015
Published: 17 June 2015
The spatial structure of a habitat can have a strong impact on community dynamics. Different experimental approaches exist to explore the effect of spatial structure on bacterial communities. To investigate the effect of ‘space’, a single implementation of spatial structure is often contrasted to bacterial community dynamics in well-mixed cultures. While such comparisons are useful, it is likely that the observed dynamics will be particular to the specific experimental implementation of spatial structure. In order to address this question, we track the community dynamics of a two-strain Escherichia coli community in various spatial habitats and relate the observed dynamics to the structure of a habitat.
By tracking the community dynamics of rpoS wild-type and mutant E. coli in radially expanding colonies on solid and semi-solid agar plates, we find that the mutant strain outcompetes the wild-type on semi-solid agar plates, whereas the two strains coexist on solid agar. We compare these results to previous studies in which the same two strains were shown to coexist in habitats spatially structured by microfabrication, while the mutant outcompeted the wild-type in well-mixed batch cultures. Together, these observations show that different implementations of space may result in qualitatively different community dynamics. Furthermore, we argue that the same competitive outcome (e.g. coexistence) may arise from distinct underlying dynamics in different experimental implementations of spatial structure.
Our observations demonstrate that different experimental implementations of spatial structure may not only lead to quantitatively different communities (changes in the relative abundance of types) but can also lead to qualitatively different outcomes of long-term community dynamics (coexistence versus extinction and loss of biodiversity).
Ecological as well as physiological interactions between cells are strongly determined by their adjacency relationships. The spatial structure of an environment sets the geometry of cell–cell interactions , and community dynamics therefore depend critically on spatial structure . Traditionally, well-mixed flasks have been the method of choice to culture microbial populations. While invaluable for microbiological experiments, this culture method is chemically homogeneous and spatially unstructured, imposing a network of interactions where all cells interact with one another in a similar (mean-field) fashion, thus lacking ecological realism. The importance of spatial structure in governing microbial interactions has long been recognized, and various experimental approaches to implement spatial structure have emerged. Experimental implementations of spatial structure are diverse and range from e.g. the periodic dispersal of unstructured subpopulations [3–5], to radially expanding colonies growing on an agar surface [6–10], and synthetic ecosystems created by microfabrication [11–18]. In contrast to cultures growing in a well-mixed flask, these spatial constructions have in common that cells interact with a limited number of partners over relatively short spatial scales. However, there are important differences between the various implementations of spatial structure, for instance concerning cellular motility and dispersal mode. The particular experimental implementation of spatial structure is thus likely to affect the competitive dynamics of its inhabitants. Yet a systematic study of the same biological community in various spatial systems is currently lacking. In order to investigate whether different implementations of spatial structure result in qualitatively different community dynamics, we study a defined biological community (wild-type and rpoS mutant Escherichia coli bacteria) in different spatial constructions while keeping parameters like growth medium and temperature constant.
The influence of spatial structure is of particular interest for bacteria that are engaged in a social dilemma, as a multitude of theoretical work suggests that spatial structure can have a profound influence on social interactions [19–27]. Although various theoretical studies suggest that spatial structure stabilizes coexistence in competitive social interactions, in practice, different implementations of spatial structure will induce distinct interaction networks and therefore may give rise to different community dynamics and long-term community structure. We investigate a community consisting of an E. coli rpoS wild-type (WT) and an rpoS819 mutant strain. The rpoS gene encodes for E. coli’s general stress regulator \(\sigma ^S\). Activity of \(\sigma ^S\) is triggered by nutrient deprivation and various other stress factors. The rpoS819 allele confers a growth advantage in stationary phase (GASP) phenotype which is characterized by a delayed entry into stationary phase [28–30]. The WT and GASP strain are engaged in a social dilemma concerning resource use and management: in a well-mixed batch culture, WT cells cooperate by collectively going into stationary phase and restrain growth before all resources are depleted, ensuring that enough resources remain available to sustain a viable population for long periods. GASP mutants on the other hand, defect on this social contract by ignoring the collective decision to restrict growth, using the scarce resources for growth instead of investing in maintenance [29, 31]. RpoS WT and GASP mutants have previously been shown to interact according to a prisoner’s dilemma game in unstructured cultures (well-mixed flasks) [12, 29]. In the absence of spatial structure, GASP cheaters were found to drive WT cooperators extinct as predicted by theoretical studies of the prisoner’s dilemma. Previous experimental studies of this E. coli community furthermore have confirmed that the WT cooperator strain can coexist with the GASP cheater strain in habitats that were spatially structured using microfabrication: Refs. [12, 32, 33] followed the community dynamics of a WT-GASP E. coli community in engineered microhabitats and demonstrated that the WT and GASP strain indeed coexist at intermediate frequencies in this type of environments. To investigate whether different implementations of spatial structure indeed lead to different dynamics of the same community, we investigated WT-GASP competition in colonies growing on agar plates.
We inoculated solid agar plates (LB medium + 1.5% agar) with droplets having various WT-GASP fractions, and tracked the abundance of both strains during the sectored growth phase (see Figure 1b). Figure 1c shows the time evolution of the GASP fraction of colonies inoculated with starting GASP fractions of 0.1, 0.5, and 0.9. These time traces reveal two points: (1) for the inoculated GASP fractions of 0.1 and 0.5, the GASP fraction has increased when sectored growth commences, (2) the GASP fraction stays fairly constant during sectored growth.
The increased GASP fraction at the onset of sectored growth relative to the inoculated GASP fraction, suggests that the GASP strain has a competitive advantage over the WT when the homeland is formed. In order to test whether the competitive advantage of the GASP strain during formation of the homeland was caused by growth rate differences, we measured the growth rate of pure WT and GASP cultures during exponential growth in liquid batch cultures. The average doubling time of the GASP strain was 45 ± 4 min (mean ± SEM, n = 8), whereas the WT strain doubled every 49 ± 6 min (mean ± SEM, n = 7). This (marginal) difference is not significant (student’s t test), and growth rate differences thus likely do not explain the increase in GASP fraction during formation of the homeland. Previous studies using well-mixed batch cultures [12, 29] have shown that the GASP mutant delays its entry into stationary phase and saturates at higher population densities relative to the WT. It is possible that a similar effect is at play during formation of the homeland, allowing the GASP strain to sustain growth while the WT limits its growth in response to crowding. In the absence of clear growth rate differences, prolonged growth by the GASP strain, and thus higher (saturation) densities, could explain the increased GASP fraction at the onset of sectored growth.
The competitive superiority of the GASP strain during homeland-growth initially translates into an increased GASP fraction, however, when growth proceeds in the sectored phase, the increase in the population share of the GASP strain levels off. Figure 1c shows that in this regime the interaction between the WT and GASP strain is neutral, i.e. the two strains coexist at constant fractions. It is informative to consider the growth state that characterizes the sectored phase. During sectored growth, the colony’s radius grows linearly  resulting in an approximately quadratic increase in population size over time. A quadratic accumulation of biomass holds the middle ground between exponential growth (governed by E. coli's housekeeping sigma factor rpoD) and stationary-phase growth (regulated by rpoS). Furthermore, as the colony expands its territory fresh nutrients are available beyond the colony’s edge, ruling out a scenario of severe nutrient deprivation at the expanding population front. The presence of fresh nutrients, together with a growth phase that is not exponential nor stationary, suggest that the rpoS regulon may not be (strongly) induced during sectored growth, rendering the rpoS WT and rpoS819 GASP mutant strains neutral. This hypothesis is in agreement with the observed neutrality of the WT and GASP strain in the sectored growth phase.
Also in the two spatial constructions where the GASP strain drives the WT strain extinct (well-mixed flasks and semi-solid plates) different mechanisms are responsible for the GASP strain’s superiority. In well-mixed flasks GASP cells outcompete WT cells by delaying their entry into stationary phase, whereas differences in swarming ability are the likely the cause of GASP superiority on semi-solid plates.
Radially expanding colonies and microhabitats both restrict cell–cell interactions to relatively small scales (when compared to a well-mixed flask), yet they differ in important respects. On solid agar, E. coli cells are immotile and dispersal in colonies results from growth and division only. The microhabitats used in Refs. [12, 32, 33] on the other hand, consist of micro-scale patches filled with liquid medium that support swimming motility. Bacteria inhabiting a microhabitat can thus be sessile (surface associated) or planktonic (free swimming) and switch between those modes, allowing cells to aggregate in biofilms and disperse from those dynamically (see Figure 3). Such dynamic self-structuring of the community in space and time is absent in colonies growing on solid agar. Other differences exist between the various implementations of spatial structure considered here. The availability of resources for instance, was varied in the three studies that used microhabitats: in Refs. [32, 33], the bacteria inhabited chambers that were connected to external nutrient reservoirs. In Ref.  and Figure 3 such external resource supply is absent (see Table 1). The situation of colony growth on an agar plate is reminiscent of the microhabitats with an external nutrient supply, as a source of fresh nutrients extends beyond the colony edge.
The observations presented here, and in previous work [12, 32, 33], emphasize the idiosyncrasy of spatial community dynamics: different implementations of space achieved by different culturing methods may not only lead to quantitatively different communities (changes in the relative abundance of types) but can also lead to qualitatively different outcomes of the long-term community dynamics (coexistence versus extinction and loss of biodiversity). Conversely, similar outcomes in different habitats could be the result of distinct competitive interactions. Although this work focused on one particular set of strains, we expect similar implications to hold in other microbial systems. It is therefore important to keep these subtleties in mind when evaluating the long-term community structure resulting from spatial competition across different laboratory constructions representing the spatial ecology of microbes.
Strains and growth conditions
All strains used in this study have been described before in Refs. [12, 32]. Wild-type strains JEK1036 (green) and JEK1037 (red) are Escherichia coli W3110 carrying the WT rpoS allele, and are lacYZ::GFPmut2 and lacYZ::mRFP, respectively. GASP mutants JEK1032 (green) and JEK1033 (red) have the GASP phenotype and carry the rpoS819 allele (described in Ref. ) and the same fluorescent markers and genetic background as strains JEK1036 and JEK1037, respectively. Prior to all experiments cells were taken from a –80°C glycerol stock and grown overnight (30°C, 200 rpm) in lysogeny broth (LB).
Growth rates of pure WT and GASP cultures were obtained by measuring the absorbance of 200 μL cultures growing at 30°C, shaken at 600 rpm in a 96 well plate using a BMG FluoStar Optima plate reader. Growth curves were log-transformed and fitted with a first-order polynomial in Matlab to obtain doubling times.
The four E. coli strains (JEK1032, JEK1033, JEK1036 and JEK1037) were grown overnight in liquid cultures (30°C, shaken at 200 rpm) and separately diluted 1/500 in LB medium containing 100 μM IPTG and grown to OD600 = 0.25. LB agar plates (1.5% agar for solid plates, 0.5% agar for semi-solid plates; supplemented with 100 μM IPTG) were inoculated with 1 μL of WT-GASP mixtures (GASP fractions of: 0.1, 0.5 and 0.9) and incubated at 30°C for 72 h. Colonies were imaged using an inverted Olympus IX81 microscope equipped with a 10× (NA = 0.3) objective, an ORCA-R2 camera (Hamamatsu) and a motorized stage controlled using μManager software . The sample was illuminated using an X-cite 120 Q light source (Lumen dynamics).
Images were background corrected and converted to polar coordinates using the center of the homeland as the origin through a custom Matlab script (see Figure 1b). The converted image was used to evaluate the GASP fraction over time during the sectored growth phase. A previous study  demonstrated that the radius of E. coli colonies expands linearly over time, we therefore mapped the obtained fraction onto a linear time axis.
FJHH, PG, RGW, CD and JEK conceived and designed the study, and drafted the manuscript. FJHH, PG and RGW performed experiments. FJHH analyzed the data. All authors read and approved the final manuscript.
We thank F. Wu, S. v. Vliet, and other members of the Bionanoscience Department for fruitful discussions and two anonymous reviewers for insightful comments. This research was supported by the European Research Council Advanced Grant #247072 to C.D.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Dieckmann U, Law R, Metz JA (2000) The The geometry of ecological interactions: simplifying spatial complexity. Cambridge University Press, CambridgeView ArticleGoogle Scholar
- Durrett R, Levin S (1994) The importance of being discrete (and spatial). Theor Popul Biol 46(3):363–394View ArticleGoogle Scholar
- Griffin A, West S, Buckling A (2004) Cooperation and competition in pathogenic bacteria. Nature 430(7003):1024–1027PubMedView ArticleGoogle Scholar
- Nahum JR, Harding BN, Kerr B (2011) Evolution of restraint in a structured rock-paper-scissors community. Proc Natl Acad Sci USA 108(Suppl 2):10831–8PubMed CentralPubMedView ArticleGoogle Scholar
- Datta MS, Korolev KS, Cvijovic I, Dudley C, Gore J (2013) Range expansion promotes cooperation in an experimental microbial metapopulation. Proc Natl Acad Sci 110(18):7354–7359PubMed CentralPubMedView ArticleGoogle Scholar
- Hallatschek O, Hersen P, Ramanathan S, Nelson DR (2007) Genetic drift at expanding frontiers promotes gene segregation. Proc Natl Acad Sci 104(50):19926–19930PubMed CentralPubMedView ArticleGoogle Scholar
- Van Dyken JD, Müller MJ, Mack KM, Desai MM (2013) Spatial population expansion promotes the evolution of cooperation in an experimental prisoner’s dilemma. Curr Biol 23:919–923PubMed CentralPubMedGoogle Scholar
- Müller MJ, Neugeboren BI, Nelson DR, Murray AW (2014) Genetic drift opposes mutualism during spatial population expansion. Proc Natl Acad Sci 111:1037–1042PubMed CentralPubMedView ArticleGoogle Scholar
- Korolev KS, Xavier JB, Nelson DR, Foster KR (2011) A quantitative test of population genetics using spatiogenetic patterns in bacterial colonies. Am Nat 178(4):538PubMed CentralPubMedView ArticleGoogle Scholar
- Rudge TJ, Federici F, Steiner PJ, Kan A, Haseloff J (2013) Cell polarity-driven instability generates self-organized, fractal patterning of cell layers. ACS Synth Biol 2:705–714PubMedView ArticleGoogle Scholar
- Keymer J, Galajda P, Muldoon C, Park S, Austin R (2006) Bacterial metapopulations in nanofabricated landscapes. Proc Natl Acad Sci USA 103(46):17290–17295PubMed CentralPubMedView ArticleGoogle Scholar
- Hol FJ, Galajda P, Nagy K, Woolthuis RG, Dekker C, Keymer JE (2013) Spatial structure facilitates cooperation in a social dilemma: Empirical evidence from a bacterial community. PloS One 8(10):77042View ArticleGoogle Scholar
- Connell JL, Ritschdorff ET, Whiteley M, Shear JB (2013) 3D printing of microscopic bacterial communities. Proc Natl Acad Sci USA 110(46):18380–5PubMed CentralPubMedView ArticleGoogle Scholar
- Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL (2014) Solutions to the public goods dilemma in bacterial biofilms. Curr Biol 24(1):50–5PubMed CentralPubMedView ArticleGoogle Scholar
- van Vliet S, Hol FJH, Weenink T, Galajda P, Keymer JE (2014) The effects of chemical interactions and culture history on the colonization of structured habitats by competing bacterial populations. BMC Microbiol 14:116PubMed CentralPubMedView ArticleGoogle Scholar
- Hol FJ, Voges MJ, Dekker C, Keymer JE (2014) Nutrient-responsive regulation determines biodiversity in a colicin-mediated bacterial community. BMC Biol 12(1):68PubMed CentralPubMedView ArticleGoogle Scholar
- Weber MF, Poxleitner G, Hebisch E, Frey E, Opitz M (2014) Chemical warfare and survival strategies in bacterial range expansions. J R Soc Interface 11(96):20140172PubMed CentralPubMedView ArticleGoogle Scholar
- Hol FJH, Dekker C (2014) Zooming in to see the bigger picture: microfluidic and nanofabrication tools to study bacteria. Science 346(6208):1251821PubMedView ArticleGoogle Scholar
- Nowak M, May R (1992) Evolutionary games and spatial chaos. Nature 359(6398):826–829View ArticleGoogle Scholar
- Lion S, Baalen M (2008) Self-structuring in spatial evolutionary ecology. Ecol Lett 11(3):277–295PubMedView ArticleGoogle Scholar
- Kreft J-U (2004) Biofilms promote altruism. Microbiology 150(8):2751–2760PubMedView ArticleGoogle Scholar
- Santos F, Pacheco J, Lenaerts T (2006) Evolutionary dynamics of social dilemmas in structured heterogeneous populations. Proc Natl Acad Sci USA 103(9):3490–3494PubMed CentralPubMedView ArticleGoogle Scholar
- Doebeli M, Hauert C (2005) Models of cooperation based on the prisoner’s dilemma and the snowdrift game. Ecol Lett 8(7):748–766View ArticleGoogle Scholar
- Nadell C, Foster K, Xavier J (2010) Emergence of spatial structure in cell groups and the evolution of cooperation. PLoS Comput Biol 6(3):1000716View ArticleGoogle Scholar
- Hauert C, Doebeli M (2004) Spatial structure often inhibits the evolution of cooperation in the snowdrift game. Nature 428(6983):643–646PubMedView ArticleGoogle Scholar
- Verbruggen E, El Mouden C, Jansa J, Akkermans G, Bücking H, West SA, Kiers ET (2012) Spatial structure and interspecific cooperation: theory and an empirical test using the mycorrhizal mutualism. Am Nat 179(5):133–146View ArticleGoogle Scholar
- Allen B, Gore J, Nowak MA (2013) Spatial dilemmas of diffusible public goods. eLife 2Google Scholar
- Zambrano M, Siegele D, Almiron M, Tormo A, Kolter R (1993) Microbial competition: Escherichia coli mutants that take over stationary phase cultures. Science 259(5102):1757–1760PubMedView ArticleGoogle Scholar
- Vulic M, Kolter R (2001) Evolutionary cheating in Escherichia coli stationary phase cultures. Genetics 158(2):519–26PubMed CentralPubMedGoogle Scholar
- Finkel SE, Kolter R (1999) Evolution of microbial diversity during prolonged starvation. Proc Natl Acad Sci 96(7):4023–4027PubMed CentralPubMedView ArticleGoogle Scholar
- Finkel S (2006) Long-term survival during stationary phase: evolution and the gasp phenotype. Nat Rev Microbiol 4(2):113–120PubMedView ArticleGoogle Scholar
- Keymer JE, Galajda P, Lambert G, Liao D, Austin RH (2008) Computation of mutual fitness by competing bacteria. Proc Natl Acad Sci USA 105(51):20269–73PubMed CentralPubMedView ArticleGoogle Scholar
- Lambert G, Liao D, Vyawahare S, Austin RH (2011) Anomalous spatial redistribution of competing bacteria under starvation conditions. J Bacteriol 193(8):1878–83PubMed CentralPubMedView ArticleGoogle Scholar
- Kussell E (2013) Evolution in microbes. Ann Rev Biophys 42:493–514View ArticleGoogle Scholar
- Patten C, Kirchhof M, Schertzberg M, Morton R, Schellhorn H (2004) Microarray analysis of RpoS-mediated gene expression in Escherichia coli K-12. Mol Genet Genomics 272(5):580–591PubMedView ArticleGoogle Scholar
- Dong T, Kirchhof M, Schellhorn H (2008) Rpos regulation of gene expression during exponential growth of Escherichia coli k12. Mol Genet Genomics 279(3):267–277PubMedView ArticleGoogle Scholar
- Edelstein A, Amodaj N, Hoover K, Vale R, Stuurman N (2010) Computer control of microscopes using μManager. Curr Protoc Mol Biol 14.20.1–14.20.17Google Scholar