PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2014; 4: 6757.
Published online 2014 October 24. doi:  10.1038/srep06757
PMCID: PMC5381366

Trophic interactions induce spatial self-organization of microbial consortia on rough surfaces

Gang Wanga,1,* and Dani Orb,1

Abstract

The spatial context of microbial interactions common in natural systems is largely absent in traditional pure culture-based microbiology. The understanding of how interdependent microbial communities assemble and coexist in limited spatial domains remains sketchy. A mechanistic model of cell-level interactions among multispecies microbial populations grown on hydrated rough surfaces facilitated systematic evaluation of how trophic dependencies shape spatial self-organization of microbial consortia in complex diffusion fields. The emerging patterns were persistent irrespective of initial conditions and resilient to spatial and temporal perturbations. Surprisingly, the hydration conditions conducive for self-assembly are extremely narrow and last only while microbial cells remain motile within thin aqueous films. The resulting self-organized microbial consortia patterns could represent optimal ecological templates for the architecture that underlie sessile microbial colonies on natural surfaces. Understanding microbial spatial self-organization offers new insights into mechanisms that sustain small-scale soil microbial diversity; and may guide the engineering of functional artificial microbial consortia.

Soil microorganisms are constrained to life in a fragmented and highly dynamic aqueous network in complex pore spaces that limit nutrient diffusion pathways, control cell dispersion rates, and shape interactions among microbial populations1,2,3,4,5,6. The wide range of trophic interactions and the inherent variability of nutrient fluxes in soil give rise to the formation of microbial consortia that are considered important for maintaining stable ecological interactions within a complex and dynamic soil environment7,8,9. Unlike interactions in well-mixed cultures, the large heterogeneity of diffusion pathways coupled with trophic interdependencies may dictate spatial self-organization for efficient exploitation of physico-chemical niches and for carrying out complex functions such as the biodegradation of compounds that require cooperation of multiple microbial species8,9,10,11. Not surprisingly, the spatial context of microbial interactions and the resulting self-assembly that commonly occur in natural systems remain poorly understood due to limitations of traditional pure culture-based microbiology11,12. In addition, eecological principles conceived for macro-ecosystems (e.g., involving plants and animals) where trophic interactions shape coexistence, biodiversity and ecological function13,14 may not capture the high dimensional nature of the nuanced and highly localized trophic interdependencies with rapidly changing microbial community composition that give rise to emergence of stable consortia11,12,15. The establishment of a quantitative framework for systematically linking environmental biophysical processes with multispecies microbial interactions in their natural spatial and temporal context is considered essential to advance our understanding of environmental microbiology4,12,16 and serves as the main objective of this study. We report a mechanistic modeling study of cell-level tropic interactions among microbial populations inhabiting hydrated soil surfaces, focusing on how trophic interdependencies are reflected in cell-level interactions and give rise to spatial self-organization of motile microbial communities. The model resolves local nutrient diffusion in space and time that shape cell-level nutrient interception, growth, and interactions among neighboring cells and their local environment. The fragmented and highly dynamic aqueous phase that support microbial community function presents radically different constraints to cell motion and interactions than the more familiar water replete biofilm studies17,18,19,20. The study expands on the recent observations and models for spatial pattern formation21 by considerations of life and interactions on physically constraining hydrated surfaces, and by considering various trophic interactions that drive formation of stable community patterns20. These expand the range of interactions previously studied such as the emergence and dispersion of successful mutants21,22, or protective segregation of cooperative organisms studied in Hallatschek et al.21.

A systematic evaluation of processes giving rise to spatially self-organized microbial consortia was based on a mechanistic modeling framework built on previous studies6,23. The present study combines concepts from microbial biology with physical representation of hydrated soil surfaces enabling mechanistic quantification of microbial nutrient uptake, growth, movement, and interactions with neighboring cells at their local (and dynamic) environments. These are essential ingredients for understanding microbial life in the natural concourse11. We first studied the simple case of a consortium consisting of two species (sp1 and sp2) that utilize two obligatory nutrients (N1 and N2) in complementary fashion. We then added complexity with a third species (sp3) that utilizes a by-product (nutrient 3 − N3) excreted by the two species as a food source in commensal or mutualistic relations (see Table 1).

Table 1
Schematic presentation of the trophic interactions and consortia

Results and Discussion

We first considered trophic interactions between two microbial species (sp1 and sp2) with different apparent yields to two obligatory nutrients (N1 and N2) present in the aqueous phase of a homogenous rough surface. The mechanistic modeling of growth, multiplication and chemotactic motility of microbial cells of the two populations spontaneously gave rise to formation of spatially segregated sectors populated by single species, as shown in Figs. 1 and and2.2. The resulting spatial patterns were remarkably stable irrespective of population growth within a sector as evidenced by the persistent angular distributions of microbial populations (seeing Fig. S1, and Supplementary Movie 1), and by calculated segregation index27 values of sp1 and sp2 (for consortium I as an example), which decreased rapidly with elapsed time, and remained stable since 10 hours after inoculation on homogeneous and heterogeneous surfaces (Fig. S2). The segregation and spatial patterns were linked to differences in physiological yields for the two nutrients and the spatial patterns of local concentrations that resulted from preferred nutrient uptake. Interestingly, differences in nutrient utilization resulted in leaky trails of nutrients (see Supplementary Moves 2 and 3) indicative of consumption efficiency commensurate with the prescribed stoichiometry and complementary nutrient utilization by the dominant species in each sector. These nutrient trails, and the presence of both nutrients at sector boundaries, induced the formation of internal microbial population bands (Fig. 3, and Supplementary Movie 2). The populated bands of the two microbial species maintained a relatively stable separation distance that reflected optimal concentration peaks based on the stoichiometric relations of each species and their respective nutrient yields. In other words, each species flourished at the optimal combination of concentrations of the obligatory N1 and N2 along opposing gradients of these two nutrients in adjacent sectors. The results depicted in Fig. 3 show agreement between the ratios of local nutrient concentrations of N1 and N2 at the sector boundaries due to spontaneous self-assembly, and the theoretical values based on trophic stoichiometry illustrating the strong coupling between physiological traits, trophic dependencies and spatial organization.

Figure 1
Simulated microbial spatial patterns on rough surfaces of different consortium types.
Figure 2
Simulated microbial spatial patterns on homogeneous and heterogeneous rough surfaces.
Figure 3
Angular distributions of microbial population fraction and associating nutrient concentrations forming during self-assembly of simulated consortia patterns.

Next, we introduced a third species (sp3) into the virtual inoculum where sp3 relies exclusively on a by-product (N3) produced by both sp1 and sp2, forming either commensal25 or mutualistic26,28 trophic interactions among the three species (see Table 1). For the commensal scenario, the addition of sp3 did not alter the original spatial self-organization of sp1 and sp2 (in the absence of sp3), and the sp3 population simply followed sp1 and sp2 (Fig. 1). The simulated patterns are reminiscent of spatial self-assembly of microbial consortia observed in other microbial systems such as on moist leaf surfaces29 or agar surfaces30. A scenario where the increased concentration of the by-product plays an inhibitory role24 for the simulated growth of sp1 and sp2 (i.e., a transition into mutualistic interactions) was clearly manifested in the resulting spatial pattern of the stable consortium. A narrow and loosely populated front consisting of sp1 and sp2 only formed ahead of the expanding main population front, and the spatial segregation between trailing internal population bands of sp1 and sp2 has been reduced (Fig. 1). The results illustrate the consequences of mutualistic trophic dynamics in which inhibition at high by-product concentrations depresses growth of sp1 and sp2, “dilutes” population density at the front, and forces a tighter and denser association with sp3 at the perimeter of the consortium. The role of mutualism in promoting soil microbial diversity has been considered previously28, and the results here illustrate a spatial manifestation of mutualistic trophic interactions.

The additional constraint generally retarded consortium expansion and resulted in proximal internal population bands of sp1 and sp2 that persisted at a close distance from each other due to the depressed “leaky” concentrations (Supplementary Movie 3). Additionally, to quantify the role of chemotaxis on microbial spatial self-organization, we conducted simulations with random cell motions. The results depicted in Fig. S3 showed complete mixing and limited spatial segregation among trophically interdependent populations for the same conditions listed in Table 1. The results highlight the importance of chemotaxis for the emergence of spatially self-organized microbial consortia.

Surface heterogeneity and patchy nutrient concentrations due to hydration dynamics and complex diffusion pathways are the rules in natural terrestrial systems1,2,3,6,12. We evaluated effects of surface heterogeneity and hydration status on consortia spatial patterns for similar trophic interactions considered above. The results illustrated emergence of spatially self-organized microbial spatial patterns similar to those forming on homogeneous surfaces (Fig. 1, and Supplementary Movie 1) albeit with suppressed population sizes and constrained dispersion even for mildly wet (−2.0 kPa) and heterogeneous surfaces (Fig. S4). No organized microbial spatial patterns emerged under drier conditions (−5.0 kPa, Fig. S4c) presumably due to suppression of microbial motility, an essential ingredient for growth and self-organization on hydrated surfaces3,6,23. The results suggest a narrow range of hydration conditions (a few kPa of matric potential) that support microbial self-motion in aqueous films31 and enable self-assembly of spatially ordered consortia. Equally important, the results indicate that the drivers for self-assembly of microbial populations on heterogeneous surfaces remain dominant despite limitations to motility and fragmentation of nutrient diffusion fields. In other words, trophic interactions and nutrient supply boundary conditions are likely to support microbial consortia, whose spatial pattern optimizes nutrient interceptions even within heterogeneous and patchy nutrient fields as seen, for example, consortium II on a heterogeneous surface (Fig. S5a). We conducted additional tests for spatial self-organization of trophically-interacting microbial populations based on experimentally derived parameters for commensal25 (Consortium IV) or mutualistic26 (Consortium V) microbial consortia. The modelling results showed well-defined spatial self-arrangement for both commensal and mutualistic microbial consortia (Figs. S5b and S5c) lending support to the generality of the principle of emergence of spatial self-assembly of interacting microbial populations. The results reflect the importance of specific metabolic stoichiometry within multi-nutrients environments that result in trophic-induced niche partitioning32,33 for promoting stable microbial spatial patterns.

The study contributes to the understanding of mechanisms for spatial self-organization of microbial communities34, by elucidating the joint roles of a threshold nutrient levels34 (in excess of cell maintenance requirements), and of microbial chemotactic motion35. Even without resolving the ubiquity of these conditions for microbial pattern formation in the broader ecological context, we can identify a theoretical biophysical characteristic length (maximal separation distance) for the onset of microbial community self-assembly. We derived an analytic prediction of such characteristic length termed the trophic interactions distance (TID) that defines the maximal initial separation distance between consortium members for the activation of trophic interactions and self-organization (Fig. 4a). The TID links hydration-mediated nutrient diffusivity with microbial motility ranges, and thus provides a predictive metric for the onset of consortium self-assembly on rough surfaces. Conditions resulting in large TID values (namely rapid cell migration or low threshold concentrations) suggest potential for self-assembly even for low densities or large separation distances among community partners; whereas low TID values indicate that consortium members must reside within this limited range to trigger self-assembly (Fig. 4b). This biophysical metric could be useful for design of bioremediation activities including selection of optimal water content, and volumetric densities and strengths of bio-stimulation. Similar estimates may be used for estimating quorum sensing ranges and thresholds36,37 in soil and other domains.

Figure 4
Analytical predictions of microbial Trophic Interactions Distance (TID) for the onset of consortium self-assembly.

The growing interest in microbial spatial association as an important ingredient for better understanding of natural microbial ecology has not yet been fully integrated into practical models due to limited experimental information8,9,11. The proposed modeling framework offers an exploratory platform for guiding information gathering and systematic evaluation of trophic dependencies and potential spatial patterns in the context of natural hydrated surfaces that could contribute to anchoring environmental microbiology back in the natural concourse11. The results show that cell-level interactions among multispecies with different trophic dependencies induce dynamic and heterogeneous localized nutrient patterns that have been observed in natural systems2,3,4,11. Considering the relatively limited temporal window of favorable hydration conditions that support microbial self-organization (for many geographic and climatic regions lasting only several hours a few times per year), the resulting stable spatial patterns could be viewed as rudimentary ecological templates for the more permanent microbial colonies forming on newly inhabited soil surfaces2,3,5,6. The results may also offer guidance for acquiring experimental observations to enhance understanding of functioning microbial communities11,12. The proposed biophysical TID for triggering spatial patterns11 may have practical applications for the design of artificial microbial consortia in the context of synthetic biology22, and for improving efficiency of bioremediation activities in the natural environment12,38.

Methods

Representing hydrated soil surfaces

We developed a spatially-explicit and individual-based model17,23 to systematically study the spatial and temporal dynamics of multispecies cell-level trophic interactions in the context of the self-assembly of microbial consortia on hydrated surfaces. The simulation domain abstracts natural soil surfaces into a two dimensional network of roughness/capillary elements arranged on a regular lattice6,23 comprised of 100 × 87 sites spanning a domain with physical size of 17.2 × 17.2 mm. We considered homogeneous roughness networks (HM) consisting of identical roughness elements/channels, and equivalent heterogeneous networks (HT) consisting of heterogeneous roughness elements drawn from a statistical distribution of sizes.

Modeling microbial growth kinetics and self-motion of individual cells

Microbial growth rate and metabolic reactions for two nutrients limiting growth were described by Monod-type kinetics as39,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m2.jpg

and

An external file that holds a picture, illustration, etc.
Object name is srep06757-m1.jpg

where μi (i = 1, 2) is effective microbial specific growth rate, μ1i and μ2i are the actual specific growth rates, and K1i and K2i are half-saturation constants for the first and second nutrients of N1i and N2i for species i, respectively (see details in Supplementary Information). Diffusion of nutrient within the hydrated roughness network was solved based on Fick's law for the domain updated by microbial nutrient uptake within each roughness element according to the reaction-diffusion model6,23,40.

The self-motion of microbial cells is an important trait that confers advantages for survival in patchy and heterogeneous environments1. Self-motion also promotes other biophysical interactions such as self-organization in response to chemotactic gradients23,41,42. Flagellated and other forms of cell motions41 on soil surfaces become rapidly restricted with reduction in soil aqueous phase content. These physical restrictions are attributed to enhanced cell-wall viscous drag in thin films, followed by capillary pinning as air-water interfaces interact with microbial cells in unsaturated soil3,23. The effects of surface hydration state on individual cell motion and on population dispersion rates were expressed by relationship between cell size and effective water film thickness23, d(ψ). For a given matric potential value (ψ), the resulting cell velocity (V) considering capillary and hydrodynamic limitations is obtained as6: V(ψ) = V0(FMFC(d(ψ)) − Fλ(d(ψ)))/FM (details see Supplementary Information). For cell motions in response to chemotactic gradients1,42, we first evaluate the hydration-constrained mean cell velocity, V(ψ), as a function of local aqueous film thickness. Next, we assign a displacement vector that depends on nutrient (chemo-attractant) gradient by weighing chemotactic and random motility components using complementary weight factors, ζ and 1-ζ, where ζ is the normalized dimensionless nutrient gradient defined as the ratio of local to maximal nutrient gradients, with ζ = 0 for entirely random cell motility6 (details see Supplementary Information).

Simulated hydration and heterogeneity scenarios

Microbial population interactions on homogeneous or on heterogeneous hydrated surface roughness networks6,31 were simulated for different hydration conditions expressed by water matric potential values of −0.5, −2.0 and −5.0 kPa, each with three replicates (we reported −0.5 kPa only for homogeneous scenarios). These hydration conditions mimic soil hydration ranging from very wet to mildly dry surfaces that limit microbial motility in aqueous films. Mixed microbial populations were randomly introduced onto a 2 × 2 mm region at the center of the domain with each species consisting of 100 individual cells (or 30 cells for each inoculation site for the scenario depicted in Fig. S5a). The initial and boundary conditions included uniform initial distributions of nutrients N1 and N2 concentrations (2 mg l−1) throughout the simulation domain (of physical size 17.2 × 17.2 mm), while maintaining zero nutrient fluxes across the domain boundaries. The exceptions to these conditions were: (i) the scenario presented in Fig. S5a with fixed point nutrient concentrations of N2 at the region's interior marked by solid circles; and (ii) for the scenarios in Figs. S5b and S5c – consortia IV and V, where only N1 was initialized throughout the simulation domain, with zero fluxes at the boundaries of the domain.

Analysis of microbial spatial segregation

Motile microbial cells may relocate within the aqueous network towards positions that improve their nutrient acquisition18. We focus on modelling results for consortium II considering the resulting spatial distributions of sp1 and sp2 (Fig. 1b, with population bands of sp1 and sp2 marked by red and green arrows). Interestingly, complementary population bands of the two microbial species form along the boundaries of the occupied sectors. The population bands were separated by relatively constant distance that reflects the nutrient utilization efficiencies of each species and their specific stoichiometric relations to each of the nutrients. Specifically, each species identified an optimal combination of the two obligatory nutrients concentrations (expressed as [N1i,res] and [N2i,res]) to satisfy the following condition,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m3.jpg

Rearranging equation (2) yields,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m4.jpg

Substituting parameter values for consortium II into equation (3) one obtains,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m6.jpg

and

An external file that holds a picture, illustration, etc.
Object name is srep06757-m5.jpg

Because the ratio An external file that holds a picture, illustration, etc.
Object name is srep06757-m12.jpg (i = 1, 2), the resulting value of [N11,res] < 0.001[mg l−1], and that of [N22,res] < 0.001[mg l−1]. The concentration of residual nutrients within each segregated band (occupied by a single microbial species) can thus be estimated, according to the respective apparent yield for a specific nutrient, as,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m8.jpg

and

An external file that holds a picture, illustration, etc.
Object name is srep06757-m7.jpg

where [N10] and [N20] are initial nutrient concentrations of N1 and N2, respectively. Substituting the parameter values into equation (5), one obtains [N2res] = [N1res] = 0.67 mg l−1. Taking into account that [N11,res] < 0.001[mg l−1] and [N22,res] < 0.001[mg l−1], the consumption amount of N2 for persisting population of sp1 (or N1 for sp2) is negligible as compared to the values of [N2res] = [N1res] = 0.67 mg l−1 according to microbial nutrient consumption stoichiometry (see Table 2). Therefore, one may set [N21,res] = [N2res] = 0.67 for sp1, and [N12,res] = [N1res] = 0.67 for sp2; and obtain [N11,res] = 0.000997 for sp1, and [N22,res] = 0.000997 for sp2 by solving equation (4). Substituting these values into equation (4), one obtains An external file that holds a picture, illustration, etc.
Object name is srep06757-m13.jpg for sp1, and An external file that holds a picture, illustration, etc.
Object name is srep06757-m14.jpg for sp2. These theoretical predictions (based entirely on microbial growth kinetics) were in very good agreement with simulated values associated with spontaneous spatial self-organization that resulted in ratios of 0.0062 for sp1, and 412 for sp2 (see Fig. 3). Considering more general scenarios for environments with low nutrient concentrations, e.g., [N1i] [double less-than sign] K1i and [N2i] [double less-than sign] K2i, common in natural environments3,7,44, equation (3) can be simplified as,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m9.jpg

The degree of microbial spatial segregation can be quantified according to a segregation index by Belmonte et al.27,

An external file that holds a picture, illustration, etc.
Object name is srep06757-m10.jpg

where n and n = are numbers of neighboring channels dominated by different and same population of species i, respectively (with species i dominates the target channel), and An external file that holds a picture, illustration, etc.
Object name is srep06757-m16.jpg donates an average over the channels that are dominated by species i.

Table 2
Physiological parameters for microbial growth, metabolism and nutrient concentrations

Analytical prediction of critical trophic interactions distance (TID)

The spatial self-organization of microbial consortia emerges through collective interactions among individual cells of consortium members and their local aqueous and nutrient environments. These interactions are shaped by acquisition of essential nutrients and by other environmental stimuli45. The spatial separation between the initially unorganized yet trophically interlinked consortium members is critical for the triggering of subsequent spatial self-organization. The key biophysical quantify of interest here, is the critical distance for activation of trophic interactions. We refer this characteristic distance as trophic interactions distance (TID), and it defines the maximal initial separation distance between consortium members for activation of trophic interactions calculated as (see details in Supplementary Information),

An external file that holds a picture, illustration, etc.
Object name is srep06757-m11.jpg

The TID reflects the interplay of hydration-mediated diffusion and motility, and threshold concentrations for setting the conditions for self-assembly and formation of consortium on heterogeneous rough surfaces.

Supplementary Material

Supplementary Information:

Movie 1

Supplementary Information:

Movie 2

Supplementary Information:

Movie 3

Supplementary Information:

Supplementary Information

Acknowledgments

We thank Martin Ackermann, Robin Tecon and Stan Schymanski (ETH, Zurich) for many useful comments and discussions. The financial support of the Swiss National Science Foundation (200021-113442), the European Research Council (ERC) Advanced Grant (320499 - SoilLife), SystemsX.ch (MicroScapes project), and the National Science Foundation of China (41401265) are gratefully acknowledged.

Footnotes

The authors declare no competing financial interests.

Author Contributions D.O. and G.W. conceived the research, G.W. performed the research, and G.W. and D.O. wrote the manuscript.

References

  • Fenchel T. Microbial behavior in a heterogeneous world. Science 296, 1068–1071 (2002). [PubMed]
  • Ramette A. & Tiedje J. M. Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc Natl Acad Sci U S A 104, 2761–2766 (2007). [PubMed]
  • Or D., Smets B. F., Wraith J. M., Dechesne A. & Friedman S. P. Physical constraints affecting bacterial habitats and activity in unsaturated porous media - a review. Adv Water Resour 30, 1505–1527 (2007).
  • Vos M., Wolf A. B., Jennings S. J. & Kowalchuk G. A. Micro-scale determinants of bacterial diversity in soil. FEMS Microbiol Rev 37, 936–954 (2013). [PubMed]
  • Pion M. et al. Gains of bacterial flagellar motility in a fungal world. Appl Environ Microbiol 79, 6862–6867 (2013). [PMC free article] [PubMed]
  • Wang G. & Or D. Hydration dynamics promote bacterial coexistence on rough surfaces. ISME J 7, 395–404 (2013). [PMC free article] [PubMed]
  • Young M. & Crawford J. W. Interactions and self-organization in the soil-microbe complex. Science 304, 1634–1637 (2004). [PubMed]
  • Woyke T. et al. Symbiosis insights through metagenomic analysis of a microbial consortium. Nature 443, 950–955 (2006). [PubMed]
  • VerBerkmoes N. C., Denef V. J., Hettich R. L. & Banfield J. F. Systems biology: functional analysis of natural microbial consortia using community proteomics. Nat Rev Microbiol 7, 196–205 (2009). [PubMed]
  • Breugelmans P. et al. Architecture and spatial organization in a triple-species bacterial biofilm synergistically degrading the phenylurea herbicide linuron. FEMS Microbiol Ecol 64, 271–282 (2008). [PubMed]
  • Wintermute E. H. & Silver P. A. Dynamics in the mixed microbial concourse. Genes Dev 24, 2603–2614 (2010). [PubMed]
  • Prosser J. I. et al. The role of ecological theory in microbial ecology. Nature 5, 384–392 (2007). [PubMed]
  • Knight T. M., McCoy M. W., Chase J. M., McCoy K. A. & Holt R. D. Trophic cascades across ecosystems. Nature 437, 880–883 (2005). [PubMed]
  • Alexandrou M. A. et al. Competition and phylogeny determine community structure in Müllerian co-mimics. Nature 469, 84–88 (2011). [PubMed]
  • Clark J. S. et al. Resolving the biodiversity paradox. Ecol Lett 10, 647–662 (2007). [PubMed]
  • Gonzalez A. et al. Our microbial selves: what ecology can teach us. EMBO Rep 12, 775–784 (2011). [PubMed]
  • Kreft J. U., Booth G. & Wimpenny J. W. T. BacSim, a simulator for individual-based modelling of bacterial colony growth. Microbiology 144, 3275–3287 (1998). [PubMed]
  • Skolimowski M. et al. Microfluidic dissolved oxygen gradient generator biochip as a useful tool in bacterial biofilm studies. Lab Chip 10, 2162–2169 (2010). [PubMed]
  • Nadell C. D., Foster K. R. & Xavier J. B. Emergence of spatial structure in cell groups and the evolution of cooperation. PLoS Comput Biol 6, e1000716 (2010). [PMC free article] [PubMed]
  • Hansen S. K., Rainey P. B., Haagensen J. A. & Molin S. Evolution of species interactions in a biofilm community. Nature 445, 533–536 (2007). [PubMed]
  • Hallatschek O., Hersen P., Ramanathan S. & Nelson D. R. Genetic drift at expanding frontiers promotes gene segregation. Proc Natl Acad Sci U S A 104, 19926–19930 (2007). [PubMed]
  • Kim H. J., Boedicker J. Q., Choi J. W. & Ismagilov R. F. Defined spatial structure stabilizes a synthetic multispecies bacterial community. Proc Natl Acad Sci U S A 105, 18188–18193 (2008). [PubMed]
  • Dechesne A., Wang G., Gülez G., Or D. & Smets B. F. Hydration-controlled bacterial motility and dispersal on surfaces. Proc Natl Acad Sci U S A 107, 14369–14372 (2010). [PubMed]
  • Leskovac V. Comprehensive Enzyme Kinetics. Kulwer Academic/Plenum Publisher: New York (2003).
  • Nielsen A. T., Tolker-Nielsen T., Barken K. B. & Molin S. Role of commensal relationships on the spatial structure of a surface-attached microbial consortium. Environ Microbiol 2, 59–68 (2000). [PubMed]
  • Kawai F. Microbial degradation of polyethers. Appl Microbiol Biotechnol 58, 30–38 (2002). [PubMed]
  • Belmonte J. M., Thomas G. L., Brunnet L. G., de Almeida R. M. C. & Chaté H. Self-propelled particle model for cell-sorting phenomena. Phys Rev Lett 100, 248702 (2008). [PubMed]
  • Pachepsky E., Crawford J. W., Bown J. L. & Squire G. Towards a general theory of biodiversity. Nature 410, 923–926 (2001). [PubMed]
  • Monier J. M. & Lindow S. E. Spatial organization of dual-species bacterial aggregates on leaf surfaces. Appl Environ Microbiol 71, 5484–5493 (2005). [PMC free article] [PubMed]
  • Momeni B., Brileya K. A., Fields M. W. & Shou W. Strong inter-population cooperation leads to partner intermixing in microbial communities. eLife 2, e00230 (2013). [PMC free article] [PubMed]
  • Wang G. & Or D. A hydration-based biophysical index for the onset of soil microbial coexistence. Sci Rep 2, 881 (2012). [PMC free article] [PubMed]
  • Krone S. M. & Guan Y. Spatial self-organization in a cyclic resource-species model. J Theor Biol 241, 14–25 (2006). [PubMed]
  • Oakley B. B., Carbonero F., Dowd S. E., Hawkins R. J. & Purdy K. J. Contrasting patterns of niche partitioning between two anaerobic terminal oxidizers of organic matter. ISME J 6, 905–914 (2012). [PMC free article] [PubMed]
  • Hibbing M. E., Fuqua C., Parsek M. R. & Peterson S. B. Bacterial competition: surviving and thriving in the microbial jungle. Nature Rev Microbiol 8, 15–25 (2010). [PMC free article] [PubMed]
  • Turnbull G. A., Morgan J. A. W., Whipps J. M. & Saunders J. R. The role of bacterial motility in the survival and spread of Pseudomonas fluorescens in soil and in the attachment and colonization of wheat roots. FEMS Microbiol Ecol 36, 21–31 (2001). [PubMed]
  • Park S. et al. Motion to form a quorum. Science 301, 188 (2003). [PubMed]
  • DeAngelis K. M., Lindow S. E. & Firestone M. K. Bacterial quorum sensing and nitrogen cycling in rhizosphere soil. FEMS Microbiol Ecol 66, 197–207 (2008). [PubMed]
  • Crawford J. W. et al. Microbial diversity affects self-organization of the soil-microbe system with consequences for function. J R Soc Interface 9, 1302–1310 (2012). [PMC free article] [PubMed]
  • Zinn M., Witholt B. & Egli T. Dual nutrient limited growth: models, experimental observations, and applications. J Biotechnol 113, 263–279 (2004). [PubMed]
  • Golding I., Kozlovsky Y., Cohen I. & Ben-Jacob E. Studies of bacterial branching growth using reaction-diffusion models for colonial development. Physica A 260, 510–554 (1998).
  • Darnton N. C. & Berg H. C. Bacterial flagella are firmly anchored. J Bacteriol 190, 8223–8224 (2008). [PMC free article] [PubMed]
  • Paul D., Singh R. & Jain R. K. Chemotaxis of Ralstonia sp. SJ98 towards p-nitrophenol in soil. Environ Microbiol 8, 1797–1804 (2006). [PubMed]
  • Lipson D. A., Monson R. K., Schmidt S. K. & Weintraub M. N. The trade-off between growth rate and yield in microbial communities and the consequences for under-snow soil respiration in a high elevation coniferous forest. Biogeochemistry 95, 23–35 (2009).
  • Hoehler T. M. & Jørgensen B. B. Microbial life under extreme energy limitation. Nat Rev Microbiol 11, 83–94 (2013). [PubMed]
  • Phelan V. V., Liu W. T., Pogliano K. & Dorrestein P. C. Microbial metabolic exchange – the chemotype-to-phenotype link. Nat Chem Biol 8, 26–35 (2012). [PMC free article] [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group