Drivers of Spatial Structure in Social Microbial Communities

Higher-order interactions and emergent properties of microbial communities: The power of synthetic ecology

Humans have long relied on microbial communities to create products, produce energy, and treat waste. The microbiota residing within our bodies directly impacts our health, while the soil and rhizosphere microbiomes influence the productivity of our crops. However, the complexity and diversity of microbial communities make them challenging to study and difficult to develop into applications, as they often exhibit the emergence of unpredictable higher-order phenomena. Synthetic ecology aims at simplifying complexity by constituting synthetic or semi-natural microbial communities with reduced diversity that become easier to study and analyze. This strategy combines methodologies that simplify existing complex systems (top-down approach) or build the system from its constituent components (bottom-up approach). Simplified communities are studied to understand how interactions among populations shape the behavior of the community and to model and predict their response to external stimuli. By harnessing the potential of synthetic microbial communities through a multidisciplinary approach, we can advance knowledge of ecological concepts and address critical public health, agricultural, and environmental issues more effectively.

Phage predation accelerates the spread of plasmid-encoded antibiotic resistance

Phage predation is generally assumed to reduce microbial proliferation while not contributing to the spread of antibiotic resistance. However, this assumption does not consider the effect of phage predation on the spatial organization of different microbial populations. Here, we show that phage predation can increase the spread of plasmid-encoded antibiotic resistance during surface-associated microbial growth by reshaping spatial organization. Using two strains of the bacterium Escherichia coli, we demonstrate that phage predation slows the spatial segregation of the strains during growth. This increases the number of cell-cell contacts and the extent of conjugation-mediated plasmid transfer between them. The underlying mechanism is that phage predation shifts the location of fastest growth from the biomass periphery to the interior where cells are densely packed and aligned closer to parallel with each other. This creates straighter interfaces between the strains that are less likely to merge together during growth, consequently slowing the spatial segregation of the strains and enhancing plasmid transfer between them. Our results have implications for the design and application of phage therapy and reveal a mechanism for how microbial functions that are deleterious to human and environmental health can proliferate in the absence of positive selection.

Adhesive interactions within microbial consortia can be differentiated at the single-cell level through expansion microscopy

Investigating microbe-microbe interactions at the single-cell level is critical to unraveling the ecology and dynamics of microbial communities. In many situations, microbes assemble themselves into densely packed multi-species biofilms. The density and complexity pose acute difficulties for visualizing individual cells and analyzing their interactions. Here, we address this problem through an unconventional application of expansion microscopy, which allows for the 'decrowding' of individual bacterial cells within a multispecies community. Expansion microscopy generally has been carried out under isotropic expansion conditions and used as a resolution-enhancing method. In our variation of expansion microscopy, we carry out expansion under heterotropic conditions; that is, we expand the space between bacterial cells but not the space within individual cells. The separation of individual bacterial cells from each other reflects the competition between the expansion force pulling them apart and the adhesion force holding them together. We employed heterotropic expansion microscopy to study the relative strength of adhesion in model biofilm communities. These included mono and dual-species Streptococcus biofilms, and a three-species synthetic community (Fusobacterium nucleatum, Streptococcus mutans, and Streptococcus sanguinis) under conditions that facilitated interspecies coaggregation. Using adhesion mutants, we investigated the interplay between F. nucleatum outer membrane protein RadD and different Streptococcus species. We also examined the Schaalia-TM7 epibiont association. Quantitative proximity analysis was used to evaluate the separation of individual microbial members. Our study demonstrates that heterotropic expansion microscopy can 'decrowd' dense biofilm communities, improve visualization of individual bacterial members, and enable analysis of microbe-microbe adhesive interactions at the single-cell level.

Multiscale selection in spatially structured populations

The spatial structure of populations is key to many (eco-)evolutionary processes. In such cases, the strength and sign of selection on a trait may depend on the spatial scale considered. An example is the evolution of altruism: selection in local environments often favours cheaters over altruists, but this can be outweighed by selection at larger scales, favouring clusters of altruists over clusters of cheaters. For populations subdivided into distinct groups, this effect is described formally by multilevel selection theory. However, many populations do not consist of non-overlapping groups but rather (self-)organize into other ecological patterns. We therefore present a mathematical framework for multiscale selection. This framework decomposes natural selection into two parts: local selection, acting within environments of a certain size, and interlocal selection, acting among them. Varying the size of the local environments subsequently allows one to measure the contribution to selection of each spatial scale. To illustrate the use of this framework, we apply it to models of the evolution of altruism and pathogen transmissibility. The analysis identifies how and to what extent ecological processes at different spatial scales contribute to selection and compete, thus providing a rigorous underpinning to eco-evolutionary intuitions.

Ecological drift during colonization drives within-host and between-host heterogeneity in an animal-associated symbiont

Specialized host-microbe symbioses canonically show greater diversity than expected from simple models, both at the population level and within individual hosts. To understand how this heterogeneity arises, we utilize the squash bug, Anasa tristis, and its bacterial symbionts in the genus Caballeronia. We modulate symbiont bottleneck size and inoculum composition during colonization to demonstrate the significance of ecological drift, the noisy fluctuations in community composition due to demographic stochasticity. Consistent with predictions from the neutral theory of biodiversity, we found that ecological drift alone can account for heterogeneity in symbiont community composition between hosts, even when 2 strains are nearly genetically identical. When acting on competing strains, ecological drift can maintain symbiont genetic diversity among different hosts by stochastically determining the dominant strain within each host. Finally, ecological drift mediates heterogeneity in isogenic symbiont populations even within a single host, along a consistent gradient running the anterior-posterior axis of the symbiotic organ. Our results demonstrate that symbiont population structure across scales does not necessarily require host-mediated selection, as it can emerge as a result of ecological drift acting on both isogenic and unrelated competitors. Our findings illuminate the processes that might affect symbiont transmission, coinfection, and population structure in nature, which can drive the evolution of host-microbe symbioses and microbe-microbe interactions within host-associated microbiomes.

The effect of cooperator recognition on competition among clones in spatially structured microbial communities

In spatially structured microbial communities, clonal growth of stationary cells passively generates clusters of related individuals. This can lead to stable cooperation without the need for recognition mechanisms. However, recent research suggests that some biofilm-forming microbes may have mechanisms of kin recognition. To explore this unexpected observation, we studied the effects of different types of cooperation in a microbial colony using spatially explicit, agent-based simulations of two interacting strains. We found scenarios that favor a form of kin recognition in spatially structured microbial communities. In the presence of a "cheater" strain, a strain with greenbeard cooperation was able to increase in frequency more than a strain with obligate cooperation. This effect was most noticeable in high density colonies and when the cooperators were not as abundant as the cheaters. We also studied whether a polychromatic greenbeard, in which cells only cooperate with their own type, could provide a numerical benefit beyond a simple, binary greenbeard. We found the greatest benefit to a polychromatic greenbeard when cooperation is highly effective. These results suggest that in some ecological scenarios, recognition mechanisms may be beneficial even in spatially structured communities.

Interspecific interactions facilitate keystone species in a multispecies biofilm that promotes plant growth

Microorganisms colonizing plant roots co-exist in complex, spatially structured multispecies biofilm communities. However, little is known about microbial interactions and the underlying spatial organization within biofilm communities established on plant roots. Here, a well-established four-species biofilm model (Stenotrophomonas rhizophila, Paenibacillus amylolyticus, Microbacterium oxydans, and Xanthomonas retroflexus, termed as SPMX) was applied to Arabidopsis roots to study the impact of multispecies biofilm on plant growth and the community spatial dynamics on the roots. SPMX co-culture notably promoted root development and plant biomass. Co-cultured SPMX increased root colonization and formed multispecies biofilms, structurally different from those formed by monocultures. By combining 16S rRNA gene amplicon sequencing and fluorescence in situ hybridization with confocal laser scanning microscopy, we found that the composition and spatial organization of the four-species biofilm significantly changed over time. Monoculture P. amylolyticus colonized plant roots poorly, but its population and root colonization were highly enhanced when residing in the four-species biofilm. Exclusion of P. amylolyticus from the community reduced overall biofilm production and root colonization of the three species, resulting in the loss of the plant growth-promoting effects. Combined with spatial analysis, this led to identification of P. amylolyticus as a keystone species. Our findings highlight that weak root colonizers may benefit from mutualistic interactions in complex communities and hereby become important keystone species impacting community spatial organization and function. This work expands the knowledge on spatial organization uncovering interspecific interactions in multispecies biofilm communities on plant roots, beneficial for harnessing microbial mutualism promoting plant growth.

Optogenetic spatial patterning of cooperation in yeast populations

Microbial communities are shaped by complex metabolic interactions such as cooperation and competition for resources. Methods to control such interactions could lead to major advances in our ability to better engineer microbial consortia for synthetic biology applications. Here, we use optogenetics to control SUC2 invertase production in yeast, thereby shaping spatial assortment of cooperator and cheater cells. Yeast cells behave as cooperators (i.e., transform sucrose into hexose, a public good) upon blue light illumination or cheaters (i.e., consume hexose produced by cooperators to grow) in the dark. We show that cooperators benefit best from the hexoses they produce when their domain size is constrained between two cut-off length-scales. From an engineering point of view, the system behaves as a bandpass filter. The lower limit is the trace of cheaters' competition for hexoses, while the upper limit is defined by cooperators' competition for sucrose. Cooperation mostly occurs at the frontiers with cheater cells, which not only compete for hexoses but also cooperate passively by letting sucrose reach cooperators. We anticipate that this optogenetic method could be applied to shape metabolic interactions in a variety of microbial ecosystems.

Spatial constraints and stochastic seeding subvert microbial arms race

Surface attached communities of microbes grow in a wide variety of environments. Often, the size of these microbial community is constrained by their physical surroundings. However, little is known about how size constraints of a colony impact the outcome of microbial competitions. Here, we use individual-based models to simulate contact killing between two bacterial strains with different killing rates in a wide range of community sizes. We found that community size has a substantial impact on outcomes; in fact, in some competitions the identity of the most fit strain differs in large and small environments. Specifically, when at a numerical disadvantage, the strain with the slow killing rate is more successful in smaller environments than in large environments. The improved performance in small spaces comes from finite size effects; stochastic fluctuations in the initial relative abundance of each strain in small environments lead to dramatically different outcomes. However, when the slow killing strain has a numerical advantage, it performs better in large spaces than in small spaces, where stochastic fluctuations now aid the fast killing strain in small communities. Finally, we experimentally validate these results by confining contact killing strains of Vibrio cholerae in transmission electron microscopy grids. The outcomes of these experiments are consistent with our simulations. When rare, the slow killing strain does better in small environments; when common, the slow killing strain does better in large environments. Together, this work demonstrates that finite size effects can substantially modify antagonistic competitions, suggesting that colony size may, at least in part, subvert the microbial arms race.

Evolutionary consequences of nascent multicellular life cycles

A key step in the evolutionary transition to multicellularity is the origin of multicellular groups as biological individuals capable of adaptation. Comparative work, supported by theory, suggests clonal development should facilitate this transition, although this hypothesis has never been tested in a single model system. We evolved 20 replicate populations of otherwise isogenic clonally reproducing 'snowflake' yeast (Δace2/∆ace2) and aggregative 'floc' yeast (GAL1p::FLO1 /GAL1p::FLO1) with daily selection for rapid growth in liquid media, which favors faster cell division, followed by selection for rapid sedimentation, which favors larger multicellular groups. While both genotypes adapted to this regime, growing faster and having higher survival during the group-selection phase, there was a stark difference in evolutionary dynamics. Aggregative floc yeast obtained nearly all their increased fitness from faster growth, not improved group survival; indicating that selection acted primarily at the level of cells. In contrast, clonal snowflake yeast mainly benefited from higher group-dependent fitness, indicating a shift in the level of Darwinian individuality from cells to groups. Through genome sequencing and mathematical modeling, we show that the genetic bottlenecks in a clonal life cycle also drive much higher rates of genetic drift-a result with complex implications for this evolutionary transition. Our results highlight the central role that early multicellular life cycles play in the process of multicellular adaptation.

Effectiveness of Pseudomonas aeruginosa type VI secretion system relies on toxin potency and type IV pili-dependent interaction

The type VI secretion system (T6SS) is an antibacterial weapon that is used by numerous Gram-negative bacteria to gain competitive advantage by injecting toxins into adjacent prey cells. Predicting the outcome of a T6SS-dependent competition is not only reliant on presence-absence of the system but instead involves a multiplicity of factors. Pseudomonas aeruginosa possesses 3 distinct T6SSs and a set of more than 20 toxic effectors with diverse functions including disruption of cell wall integrity, degradation of nucleic acids or metabolic impairment. We generated a comprehensive collection of mutants with various degrees of T6SS activity and/or sensitivity to each individual T6SS toxin. By imaging whole mixed bacterial macrocolonies, we then investigated how these P. aeruginosa strains gain a competitive edge in multiple attacker/prey combinations. We observed that the potency of single T6SS toxin varies significantly from one another as measured by monitoring the community structure, with some toxins acting better in synergy or requiring a higher payload. Remarkably the degree of intermixing between preys and attackers is also key to the competition outcome and is driven by the frequency of contact as well as the ability of the prey to move away from the attacker using type IV pili-dependent twitching motility. Finally, we implemented a computational model to better understand how changes in T6SS firing behaviours or cell-cell contacts lead to population level competitive advantages, thus providing conceptual insight applicable to all types of contact-based competition.

Indirect Enrichment of Desirable, but Less Fit Phenotypes, from a Synthetic Microbial Community Using Microdroplet Confinement

Spatial structure within microbial communities can provide nearly limitless opportunities for social interactions and are an important driver for evolution. As metabolites are often molecular signals, metabolite diffusion within microbial communities can affect the composition and dynamics of the community in a manner that can be challenging to deconstruct. We used encapsulation of a synthetic microbial community within microdroplets to investigate the effects of spatial structure and metabolite diffusion on population dynamics and to examine the effects of cheating by one member of the community. The synthetic community was composed of three strains: a "Producer" that makes the diffusible quorum sensing molecule (N-(3-oxododecanoyl)-l-homoserine lactone, C12-oxo-HSL) or AHL; a "Receiver" that is killed by AHL; and a Non-Producer or "cheater" that benefits from the extinction of the Receivers, but without the costs associated with the AHL synthesis. We demonstrate that despite rapid diffusion of AHL between microdroplets, the spatial structure imposed by the microdroplets allows a more efficient but transient enrichment of more rare and slower-growing Producer subpopulations. Eventually, the Non-Producer population drove the Producers to extinction. By including fluorescence-activated microdroplet sorting and providing sustained competition by the Receiver strain, we demonstrate a strategy for indirect enrichment of a rare and unlabeled Producer. The ability to screen and enrich metabolite Producers from a much larger population under conditions of rapid diffusion provides an important framework for the development of applications in synthetic ecology and biotechnology.

Formation and emergent dynamics of spatially organized microbial systems

Spatial organization is the norm rather than the exception in the microbial world. While the study of microbial physiology has been dominated by studies in well-mixed cultures, there is now increasing interest in understanding the role of spatial organization in microbial physiology, coexistence and evolution. Where studied, spatial organization has been shown to influence all three of these aspects. In this mini review and perspective article, we emphasize that the dynamics within spatially organized microbial systems (SOMS) are governed by feedbacks between local physico-chemical conditions, cell physiology and movement, and evolution. These feedbacks can give rise to emergent dynamics, which need to be studied through a combination of spatio-temporal measurements and mathematical models. We highlight the initial formation of SOMS and their emergent dynamics as two open areas of investigation for future studies. These studies will benefit from the development of model systems that can mimic natural ones in terms of species composition and spatial structure.

Scaling up and down: movement ecology for microorganisms

Movement is critical for the fitness of organisms, both large and small. It dictates how individuals acquire resources, evade predators, exchange genetic material, and respond to stressful environments. Movement also influences ecological and evolutionary dynamics at higher organizational levels, such as populations and communities. However, the links between individual motility and the processes that generate and maintain microbial diversity are poorly understood. Movement ecology is a framework linking the physiological and behavioral properties of individuals to movement patterns across scales of space, time, and biological organization. By synthesizing insights from cell biology, ecology, and evolution, we expand theory from movement ecology to predict the causes and consequences of microbial movements.

Assembled mixed co-cultures for emerging pollutant removal using native microorganisms from sewage sludge

The global pharmaceutical pollution caused by drug consumption (>100,000 tonnes) and its disposal into the environment is an issue which is currently being addressed by bioremediation techniques, using single or multiple microorganisms. Nevertheless, the low efficiency and the selection of non-compatible species interfere with the success of this methodology. This paper proposes a novel way of obtaining an effective multi-domain co-culture, with the capacity to degrade multi-pharmaceutical compounds simultaneously. To this end, seven microorganisms (fungi and bacteria) previously isolated from sewage sludge were investigated to enhance their degradation performance. All seven strains were factorially mixed and used to assemble different artificial co-cultures. Consequently, 127 artificial co-cultures were established and ranked, based on their fitness performance, by using the BSocial analysis web tool. The individual strains were categorized according to their social behaviour, whose net effect over the remaining strains was defined as 'Positive', 'Negative' or 'Neutral'. To evaluate the emerging-pollutant degradation rate, the best 10 co-cultures, and those which contained the social strains were then challenged with three different Pharmaceutical Active compounds (PhACs): diclofenac, carbamazepine and ketoprofen. The co-cultures with the fungi Penicillium oxalicum XD-3.1 and Penicillium rastrickii were able to degrade PhACs. However, the highest performance (>80% degradation) was obtained by the minimal active microbial consortia consisting of both Penicillium spp., Cladosporium cladosporoides and co-existing bacteria. These consortia transformed the PhACs to derivate molecules through hydroxylation and were released to the media, resulting in a low ecotoxicity effect. High-throughput screening of co-cultures provides a quick, reliable and efficient method to narrow down suitable degradation co-cultures for emerging PhAC contaminants while avoiding toxic metabolic derivatives.

Eco-evolutionary modelling of microbial syntrophy indicates the robustness of cross-feeding over cross-facilitation

Syntrophic cooperation among prokaryotes is ubiquitous and diverse. It relies on unilateral or mutual aid that may be both catalytic and metabolic in nature. Hypotheses of eukaryotic origins claim that mitochondrial endosymbiosis emerged from mutually beneficial syntrophy of archaeal and bacterial partners. However, there are no other examples of prokaryotic syntrophy leading to endosymbiosis. One potential reason is that when externalized products become public goods, they incite social conflict due to selfish mutants that may undermine any mutualistic interactions. To rigorously evaluate these arguments, here we construct a general mathematical framework of the ecology and evolution of different types of syntrophic partnerships. We do so both in a general microbial and in a eukaryogenetic context. Studying the case where partners cross-feed on each other's self-inhibiting waste, we show that cooperative partnerships will eventually dominate over selfish mutants. By contrast, systems where producers actively secrete enzymes that cross-facilitate their partners' resource consumption are not robust against cheaters over evolutionary time. We conclude that cross-facilitation is unlikely to provide an adequate syntrophic origin for endosymbiosis, but that cross-feeding mutualisms may indeed have played that role.

Indirect enrichment of desirable, but less fit phenotypes, from a synthetic microbial community using microdroplet confinement

Spatial structure within microbial communities can provide nearly limitless opportunities for social interactions and are an important driver for evolution. As metabolites are often molecular signals, metabolite diffusion within microbial communities can affect the composition and dynamics of the community in a manner that can be challenging to deconstruct. We used encapsulation of a synthetic microbial community within microdroplets to investigate the effects of spatial structure and metabolite diffusion on population dynamics and to examine the effects of cheating by one member of the community. The synthetic community was comprised of three strains: a 'Producer' that makes the diffusible quorum sensing molecule ( N -(3-Oxododecanoyl)-L-homoserine lactone, C12-oxo-HSL) or AHL; a 'Receiver' that is killed by AHL and a Non-Producer or 'cheater' that benefits from the extinction of the Receivers, but without the costs associated with the AHL synthesis. We demonstrate that despite rapid diffusion of AHL between microdroplets, the spatial structure imposed by the microdroplets allow a more efficient but transient enrichment of more rare and slower growing 'Producer' subpopulations. Eventually, the Non-Producer population drove the Producers to extinction. By including fluorescence-activated microdroplet sorting and providing sustained competition by the Receiver strain, we demonstrate a strategy for indirect enrichment of a rare and unlabeled Producer. The ability to screen and enrich metabolite Producers from a much larger population under conditions of rapid diffusion provides an important framework for the development of applications in synthetic ecology and biotechnology.

Abstract Figure

Multispecies biofilm architecture determines bacterial exposure to phages

Numerous ecological interactions among microbes-for example, competition for space and resources, or interaction among phages and their bacterial hosts-are likely to occur simultaneously in multispecies biofilm communities. While biofilms formed by just a single species occur, multispecies biofilms are thought to be more typical of microbial communities in the natural environment. Previous work has shown that multispecies biofilms can increase, decrease, or have no measurable impact on phage exposure of a host bacterium living alongside another species that the phages cannot target. The reasons underlying this variability are not well understood, and how phage-host encounters change within multispecies biofilms remains mostly unexplored at the cellular spatial scale. Here, we study how the cellular scale architecture of model 2-species biofilms impacts cell-cell and cell-phage interactions controlling larger scale population and community dynamics. Our system consists of dual culture biofilms of Escherichia coli and Vibrio cholerae under exposure to T7 phages, which we study using microfluidic culture, high-resolution confocal microscopy imaging, and detailed image analysis. As shown previously, sufficiently mature biofilms of E. coli can protect themselves from phage exposure via their curli matrix. Before this stage of biofilm structural maturity, E. coli is highly susceptible to phages; however, we show that these bacteria can gain lasting protection against phage exposure if they have become embedded in the bottom layers of highly packed groups of V. cholerae in co-culture. This protection, in turn, is dependent on the cell packing architecture controlled by V. cholerae biofilm matrix secretion. In this manner, E. coli cells that are otherwise susceptible to phage-mediated killing can survive phage exposure in the absence of de novo resistance evolution. While co-culture biofilm formation with V. cholerae can confer phage protection to E. coli, it comes at the cost of competing with V. cholerae and a disruption of normal curli-mediated protection for E. coli even in dual species biofilms grown over long time scales. This work highlights the critical importance of studying multispecies biofilm architecture and its influence on the community dynamics of bacteria and phages.

Mechanisms driving spatial distribution of residents in colony biofilms: an interdisciplinary perspective

Biofilms are consortia of microorganisms that form collectives through the excretion of extracellular matrix compounds. The importance of biofilms in biological, industrial and medical settings has long been recognized due to their emergent properties and impact on surrounding environments. In laboratory situations, one commonly used approach to study biofilm formation mechanisms is the colony biofilm assay, in which cell communities grow on solid-gas interfaces on agar plates after the deposition of a population of founder cells. The residents of a colony biofilm can self-organize to form intricate spatial distributions. The assay is ideally suited to coupling with mathematical modelling due to the ability to extract a wide range of metrics. In this review, we highlight how interdisciplinary approaches have provided deep insights into mechanisms causing the emergence of these spatial distributions from well-mixed inocula.

Calibrating spatiotemporal models of microbial communities to microscopy data: A review

Spatiotemporal models that account for heterogeneity within microbial communities rely on single-cell data for calibration and validation. Such data, commonly collected via microscopy and flow cytometry, have been made more accessible by recent advances in microfluidics platforms and data processing pipelines. However, validating models against such data poses significant challenges. Validation practices vary widely between modelling studies; systematic and rigorous methods have not been widely adopted. Similar challenges are faced by the (macrobial) ecology community, in which systematic calibration approaches are often employed to improve quantitative predictions from computational models. Here, we review single-cell observation techniques that are being applied to study microbial communities and the calibration strategies that are being employed for accompanying spatiotemporal models. To facilitate future calibration efforts, we have compiled a list of summary statistics relevant for quantifying spatiotemporal patterns in microbial communities. Finally, we highlight some recently developed techniques that hold promise for improved model calibration, including algorithmic guidance of summary statistic selection and machine learning approaches for efficient model simulation.

Spatial population genetics with fluid flow

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2D), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1D and 2D flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how thesegeneticdrops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

Evolution-proof inhibitors of public good cooperation: a screening strategy inspired by social evolution theory

Interference with public good cooperation provides a promising novel antimicrobial strategy since social evolution theory predicts that resistant mutants will be counter-selected if they share the public benefits of their resistance with sensitive cells in the population. Although this hypothesis is supported by a limited number of pioneering studies, an extensive body of more fundamental work on social evolution describes a multitude of mechanisms and conditions that can stabilize public behaviour, thus potentially allowing resistant mutants to thrive. In this paper we theorize on how these different mechanisms can influence the evolution of resistance against public good inhibitors. Based hereon, we propose an innovative 5-step screening strategy to identify novel evolution-proof public good inhibitors, which involves a systematic evaluation of the exploitability of public goods under the most relevant experimental conditions, as well as a careful assessment of the most optimal way to interfere with their action. Overall, this opinion paper is aimed to contribute to long-term solutions to fight bacterial infections.

Mechanisms Underlying Vibrio cholerae Biofilm Formation and Dispersion

Biofilms are a widely observed growth mode in which microbial communities are spatially structured and embedded in a polymeric extracellular matrix. Here, we focus on the model bacterium Vibrio cholerae and summarize the current understanding of biofilm formation, including initial attachment, matrix components, community dynamics, social interactions, molecular regulation, and dispersal. The regulatory network that orchestrates the decision to form and disperse from biofilms coordinates various environmental inputs. These cues are integrated by several transcription factors, regulatory RNAs, and second-messenger molecules, including bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP). Through complex mechanisms, V. cholerae weighs the energetic cost of forming biofilms against the benefits of protection and social interaction that biofilms provide. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Spatial-temporal dynamics of a microbial cooperative behavior resistant to cheating

Much of our understanding of bacterial behavior stems from studies in liquid culture. In nature, however, bacteria frequently live in densely packed spatially-structured communities. How does spatial structure affect bacterial cooperative behaviors? In this work, we examine rhamnolipid production—a cooperative and virulent behavior of Pseudomonas aeruginosa. Here we show that, in striking contrast to well-mixed liquid culture, rhamnolipid gene expression in spatially-structured colonies is strongly associated with colony specific growth rate, and is impacted by perturbation with diffusible quorum signals. To interpret these findings, we construct a data-driven statistical inference model which captures a length-scale of bacterial interaction that develops over time. Finally, we find that perturbation of P. aeruginosa swarms with quorum signals preserves the cooperating genotype in competition, rather than creating opportunities for cheaters. Overall, our data demonstrate that the complex response to spatial localization is key to preserving bacterial cooperative behaviors.

Bacteria often live in densely packed, spatially-structured communities; however, much of our understanding of their behavior stems from studies in liquid culture. Here, Monaco et al. show how spatial structure and quorum sensing modulate a cooperative behavior in colonies of Pseudomonas aeruginosa.

Environmental structure impacts microbial composition and secondary metabolism

Determining the drivers of microbial community assembly is a central theme of microbial ecology, and chemical ecologists seek to characterize how secondary metabolites mediate these assembly patterns. Environmental structure affects how communities assemble and what metabolic pathways aid in that assembly. Here, we bridged these two perspectives by addressing the chemical drivers of community assembly within a spatially structured landscape with varying oxygen availability. We hypothesized that structured environments would favor higher microbial diversity and metabolite diversity. We anticipated that the production of a compound would be more advantageous in a structured environment (less mixing) compared to an unstructured environment (more mixing), where the molecule would have a diminished local effect. We observed this to be partially true in our experiments: structured environments had similar microbial diversity compared to unstructured environments but differed significantly in the metabolites produced. We also found that structured environments selected for communities with higher evenness, rather than communities with higher richness. This supports the idea that when characterizing the drivers of community assembly, it matters less about who is there and more about what they are doing. Overall, these data contribute to a growing effort to approach microbial community assembly with interdisciplinary tools and perspectives.

Antagonism between killer yeast strains as an experimental model for biological nucleation dynamics

Antagonistic interactions are widespread in the microbial world and affect microbial evolutionary dynamics. Natural microbial communities often display spatial structure, which affects biological interactions, but much of what we know about microbial antagonism comes from laboratory studies of well-mixed communities. To overcome this limitation, we manipulated two killer strains of the budding yeast Saccharomyces cerevisiae, expressing different toxins, to independently control the rate at which they released their toxins. We developed mathematical models that predict the experimental dynamics of competition between toxin-producing strains in both well-mixed and spatially structured populations. In both situations, we experimentally verified theory's prediction that a stronger antagonist can invade a weaker one only if the initial invading population exceeds a critical frequency or size. Finally, we found that toxin-resistant cells and weaker killers arose in spatially structured competitions between toxin-producing strains, suggesting that adaptive evolution can affect the outcome of microbial antagonism in spatial settings.

Snowdrift game induces pattern formation in systems of self-propelled particles

Evolutionary games between species are known to lead to intriguing spatiotemporal patterns in systems of diffusing agents. However, the role of interspecies interactions is hardly studied when agents are (self-)propelled, as is the case in many biological systems. Here, we combine aspects from active matter and evolutionary game theory and study a system of two species whose individuals are (self-)propelled and interact through a snowdrift game. We derive hydrodynamic equations for the density and velocity fields of both species from which we identify parameter regimes in which one or both species form macroscopic orientational order as well as regimes of propagating wave patterns. Interestingly, we find simultaneous wave patterns in both species that result from the interplay between alignment and snowdrift interactions-a feedback mechanism that we call game-induced pattern formation. We test these results in agent-based simulations and confirm the different regimes of order and spatiotemporal patterns as well as game-induced pattern formation.

Evolutionary ecology theory - microbial population structure

Microbial populations typically show a large degree of intra-population diversity. This diversity is intertwined with the structure of the population. Here, we discuss endogenous and exogenous drivers of population structure in microbes and how the population structure can affect evolutionary dynamics and vice versa. Endogenous structure, which can be genetic or demographic, is driven by the ecology and evolutionary dynamics within the population. Exogenous structure is typically driven by the spatial and temporal properties of the environment. A particular interesting case arises when also this exogenous structure experiences feedbacks from the microbial population.

Geometrically regulating evolutionary dynamics in biofilms

The theoretical understanding of evolutionary dynamics in spatially structured populations often relies on nonspatial models. Biofilms are among such populations where a more accurate understanding is of theoretical interest and can reveal new solutions to existing challenges. Here, we studied how the geometry of the environment affects the evolutionary dynamics of expanding populations, using the Eden model. Our results show that fluctuations of subpopulations during range expansion in two- and three-dimensional environments are not Brownian. Furthermore, we found that the substrate's geometry interferes with the evolutionary dynamics of populations that grow upon it. Inspired by these findings, we propose a periodically wedged pattern on surfaces prone to develop biofilms. On such patterned surfaces, natural selection becomes less effective and beneficial mutants would have a harder time establishing. Additionally, this modification accelerates genetic drift and leads to less diverse biofilms. Both interventions are highly desired for biofilms.

Metabolic Profiling of Rhizobacteria Serratia plymuthica and Bacillus subtilis Revealed Intra- and Interspecific Differences and Elicitation of Plipastatins and Short Peptides Due to Co-cultivation

Rhizobacteria live in diverse and dynamic communities having a high impact on plant growth and development. Due to the complexity of the microbial communities and the difficult accessibility of the rhizosphere, investigations of interactive processes within this bacterial network are challenging. In order to better understand causal relationships between individual members of the microbial community of plants, we started to investigate the inter- and intraspecific interaction potential of three rhizobacteria, the S. plymuthica isolates 4Rx13 and AS9 and B. subtilis B2g, using high resolution mass spectrometry based metabolic profiling of structured, low-diversity model communities. We found that by metabolic profiling we are able to detect metabolite changes during cultivation of all three isolates. The metabolic profile of S. plymuthica 4Rx13 differs interspecifically to B. subtilis B2g and surprisingly intraspecifically to S. plymuthica AS9. Thereby, the release of different secondary metabolites represents one contributing factor of inter- and intraspecific variations in metabolite profiles. Interspecific co-cultivation of S. plymuthica 4Rx13 and B. subtilis B2g showed consistently distinct metabolic profiles compared to mono-cultivated species. Thereby, putative known and new variants of the plipastatin family are increased in the co-cultivation of S. plymuthica 4Rx13 and B. subtilis B2g. Interestingly, intraspecific co-cultivation of S. plymuthica 4Rx13 and S. plymuthica AS9 revealed a distinct interaction zone and showed distinct metabolic profiles compared to mono-cultures. Thereby, several putative short proline-containing peptides are increased in co-cultivation of S. plymuthica 4Rx13 with S. plymuthica AS9 compared to mono-cultivated strains. Our results demonstrate that the release of metabolites by rhizobacteria alters due to growth and induced by social interactions between single members of the microbial community. These results form a basis to elucidate the functional role of such interaction-triggered compounds in establishment and maintenance of microbial communities and can be applied under natural and more realistic conditions, since rhizobacteria also interact with the plant itself and many other members of plant and soil microbiota.

Modular Molecular Weaponry Plays a Key Role in Competition Within an Environmental Vibrio cholerae Population

The type VI secretion system (T6SS) operons of Vibrio cholerae contain extraordinarily diverse arrays of toxic effector and cognate immunity genes, which are thought to play an important role in the environmental lifestyle and adaptation of this human pathogen. Through the T6SS, proteinaceous "spears" tipped with antibacterial effectors are injected into adjacent cells, killing those not possessing immunity proteins to these effectors. Here, we investigate the T6SS-mediated dynamics of bacterial competition within a single environmental population of V. cholerae. We show that numerous members of a North American V. cholerae population possess strain-specific repertoires of cytotoxic T6SS effector and immunity genes. Using pairwise competition assays, we demonstrate that the vast majority of T6SS-mediated duels end in stalemates between strains with different T6SS repertoires. However, horizontally acquired effector and immunity genes can significantly alter the outcome of these competitions. Frequently observed horizontal gene transfer events can both increase or reduce competition between distantly related strains by homogenizing or diversifying the T6SS repertoire. Our results also suggest temperature-dependent outcomes in T6SS competition, with environmental isolates faring better against a pathogenic strain under native conditions than under those resembling a host-associated environment. Taken altogether, these interactions produce density-dependent fitness effects and a constant T6SS-mediated arms race in individual V. cholerae populations, which could ultimately preserve intraspecies diversity. Since T6SSs are widespread, we expect within-population diversity in T6SS repertoires and the resulting competitive dynamics to be a common theme in bacterial species harboring this machinery.

Nutrient complexity triggers transitions between solitary and colonial growth in bacterial populations

Microbial populations often experience fluctuations in nutrient complexity in their natural environment such as between high molecular weight polysaccharides and simple monosaccharides. However, it is unclear if cells can adopt growth behaviors that allow individuals to optimally respond to differences in nutrient complexity. Here, we directly control nutrient complexity and use quantitative single-cell analysis to study the growth dynamics of individuals within populations of the aquatic bacterium Caulobacter crescentus. We show that cells form clonal microcolonies when growing on the polysaccharide xylan, which is abundant in nature and degraded using extracellular cell-linked enzymes; and disperse to solitary growth modes when the corresponding monosaccharide xylose becomes available or nutrients are exhausted. We find that the cellular density required to achieve maximal growth rates is four-fold higher on xylan than on xylose, indicating that aggregating is advantageous on polysaccharides. When collectives on xylan are transitioned to xylose, cells start dispersing, indicating that colony formation is no longer beneficial and solitary behaviors might serve to reduce intercellular competition. Our study demonstrates that cells can dynamically tune their behaviors when nutrient complexity fluctuates, elucidates the quantitative advantages of distinct growth behaviors for individual cells and indicates why collective growth modes are prevalent in microbial populations.

Glucose confers protection to Escherichia coli against contact killing by Vibrio cholerae

Evolutionary arms races are broadly prevalent among organisms including bacteria, which have evolved defensive strategies against various attackers. A common microbial aggression mechanism is the type VI secretion system (T6SS), a contact-dependent bacterial weapon used to deliver toxic effector proteins into adjacent target cells. Sibling cells constitutively express immunity proteins that neutralize effectors. However, less is known about factors that protect non-sibling bacteria from T6SS attacks independently of cognate immunity proteins. In this study, we observe that human Escherichia coli commensal strains sensitive to T6SS attacks from Vibrio cholerae are protected when co-cultured with glucose. We confirm that glucose does not impair V. cholerae T6SS activity. Instead, we find that cells lacking the cAMP receptor protein (CRP), which regulates expression of hundreds of genes in response to glucose, survive significantly better against V. cholerae T6SS attacks even in the absence of glucose. Finally, we show that the glucose-mediated T6SS protection varies with different targets and killers. Our findings highlight the first example of an extracellular small molecule modulating a genetically controlled response for protection against T6SS attacks. This discovery may have major implications for microbial interactions during pathogen-host colonization and survival of bacteria in environmental communities.

Community lifespan, niche expansion and the evolution of interspecific cooperation

Microbes live in dense and diverse communities where they deploy many traits that promote the growth and survival of neighbouring species, all the while also competing for shared resources. Because microbial communities are highly dynamic, the costs and benefits of species interactions change over the growth cycle of a community. How mutualistic interactions evolve under such demographic and ecological conditions is still poorly understood. Here, we develop an eco-evolutionary model to explore how different forms of helping with distinct fitness effects (rate-enhancing and yield-enhancing) affect the multiple phases of community growth, and its consequences for the evolution of mutualisms. We specifically focus on a form of yield-enhancing trait in which cooperation augments the common pool of resources, termed niche expansion. We show that although mutualisms in which cooperation increases partners growth rate are generally favoured at early stages of community growth, niche expansion can evolve at later stages where densities are high. Further, we find that niche expansion can promote the evolution of reproductive restraint, in which a focal species adaptively reduces its own growth rate to increase the density of partner species. Our findings suggest that yield-enhancing mutualisms are more prevalent in stable habitats with a constant supply of resources, and where populations typically live at high densities. In general, our findings highlight the need to integrate different components of population growth in the analysis of mutualisms to understand the composition and function of microbial communities.

Accumulation of dead cells from contact killing facilitates coexistence in bacterial biofilms

Bacterial communities are governed by a wide variety of social interactions, some of which are antagonistic with potential significance for bacterial warfare. Several antagonistic mechanisms, such as killing via the type VI secretion system (T6SS), require killer cells to directly contact target cells. The T6SS is hypothesized to be a highly potent weapon, capable of facilitating the invasion and defence of bacterial populations. However, we find that the efficacy of contact killing is severely limited by the material consequences of cell death. Through experiments with Vibrio cholerae strains that kill via the T6SS, we show that dead cell debris quickly accumulates at the interface that forms between competing strains, preventing physical contact and thus preventing killing. While previous experiments have shown that T6SS killing can reduce a population of target cells by as much as 106-fold, we find that, as a result of the formation of dead cell debris barriers, the impact of contact killing depends sensitively on the initial concentration of killer cells. Killer cells are incapable of invading or eliminating competitors on a community level. Instead, bacterial warfare itself can facilitate coexistence between nominally antagonistic strains. While a variety of defensive strategies against microbial warfare exist, the material consequences of cell death provide target cells with their first line of defence.

The Evolution of Mutualistic Dependence

While the importance of mutualisms across the tree of life is recognized, it is not understood why some organisms evolve high levels of dependence on mutualistic partnerships, while other species remain autonomous or retain or regain minimal dependence on partners. We identify four main pathways leading to the evolution of mutualistic dependence. Then, we evaluate current evidence for three predictions: ( a) Mutualisms with different levels of dependence have distinct stabilizing mechanisms against exploitation and cheating, ( b) less dependent mutualists will return to autonomy more often than those that are highly dependent, and ( c) obligate mutualisms should be less context dependent than facultative ones. Although we find evidence supporting all three predictions, we stress that mutualistic partners follow diverse paths toward—and away from—dependence. We also highlight the need to better examine asymmetry in partner dependence. Recognizing how variation in dependence influences the stability, breakdown, and context dependence of mutualisms generates new hypotheses regarding how and why the benefits of mutualistic partnerships differ over time and space.

Nonreciprocity as a generic route to traveling states

We examine a nonreciprocally coupled dynamical model of a mixture of two diffusing species. We demonstrate that nonreciprocity, which is encoded in the model via antagonistic cross-diffusivities, provides a generic mechanism for the emergence of traveling patterns in purely diffusive systems with conservative dynamics. In the absence of nonreciprocity, the binary fluid mixture undergoes a phase transition from a homogeneous mixed state to a demixed state with spatially separated regions rich in one of the two components. Above a critical value of the parameter tuning nonreciprocity, the static demixed pattern acquires a finite velocity, resulting in a state that breaks both spatial and time-reversal symmetry, as well as the reflection parity of the static pattern. We elucidate the generic nature of the transition to traveling patterns using a minimal model that can be studied analytically. Our work has direct relevance to nonequilibrium assembly in mixtures of chemically interacting colloids that are known to exhibit nonreciprocal effective interactions, as well as to mixtures of active and passive agents where traveling states of the type predicted here have been observed in simulations. It also provides insight on transitions to traveling and oscillatory states seen in a broad range of nonreciprocal systems with nonconservative dynamics, from reaction-diffusion and prey-predators models to multispecies mixtures of microorganisms with antagonistic interactions.

The role of spatial structure in the evolution of viral innate immunity evasion: A diffusion-reaction cellular automaton model

Most viruses have evolved strategies for preventing interferon (IFN) secretion and evading innate immunity. Recent work has shown that viral shutdown of IFN secretion can be viewed as a social trait, since the ability of a given virus to evade IFN-mediated immunity depends on the phenotype of neighbor viruses. Following this idea, we investigate the role of spatial structure in the evolution of innate immunity evasion. For this, we model IFN signaling and viral spread using a spatially explicit approximation that combines a diffusion-reaction model and cellular automaton. Our results indicate that the benefits of preventing IFN secretion for a virus are strongly determined by spatial structure through paracrine IFN signaling. Therefore, innate immunity evasion can evolve as a cooperative or even altruistic trait based on indirect fitness effects that IFN shutdown exerts on other members of the viral population. We identify key factors determining whether evasion from IFN-mediated immunity should evolve, such as population bottlenecks occurring during viral transmission, the relative speed of cellular infection and IFN secretion, and the diffusion properties of the medium.

Cooperation in Microbial Populations: Theory and Experimental Model Systems

Cooperative behavior, the costly provision of benefits to others, is common across all domains of life. This review article discusses cooperative behavior in the microbial world, mediated by the exchange of extracellular products called public goods. We focus on model species for which the production of a public good and the related growth disadvantage for the producing cells are well described. To unveil the biological and ecological factors promoting the emergence and stability of cooperative traits we take an interdisciplinary perspective and review insights gained from both mathematical models and well-controlled experimental model systems. Ecologically, we include crucial aspects of the microbial life cycle into our analysis and particularly consider population structures where ensembles of local communities (subpopulations) continuously emerge, grow, and disappear again. Biologically, we explicitly consider the synthesis and regulation of public good production. The discussion of the theoretical approaches includes general evolutionary concepts, population dynamics, and evolutionary game theory. As a specific but generic biological example, we consider populations of Pseudomonas putida and its regulation and use of pyoverdines, iron scavenging molecules, as public goods. The review closes with an overview on cooperation in spatially extended systems and also provides a critical assessment of the insights gained from the experimental and theoretical studies discussed. Current challenges and important new research opportunities are discussed, including the biochemical regulation of public goods, more realistic ecological scenarios resembling native environments, cell-to-cell signaling, and multispecies communities.

Nucleation of antagonistic organisms and cellular competitions on curved, inflating substrates

We consider the dynamics of spatially distributed, diffusing populations of organisms with antagonistic interactions. These interactions are found on many length scales, ranging from kilometer-scale animal range dynamics with selection against hybrids to micron-scale interactions between poison-secreting microbial populations. We find that the dynamical line tension at the interface between antagonistic organisms suppresses survival probabilities of small clonal clusters: the line tension introduces a critical cluster size that an organism with a selective advantage must achieve before deterministically spreading through the population. We calculate the survival probability as a function of selective advantage δ and antagonistic interaction strength σ. Unlike a simple Darwinian selective advantage, the survival probability depends strongly on the spatial diffusion constant D_{s} of the strains when σ>0, with suppressed survival when both species are more motile. Finally, we study the survival probability of a single mutant cell at the frontier of a growing spherical cluster of cells, such as the surface of an avascular spherical tumor. Both the inflation and curvature of the frontier significantly enhance the survival probability by changing the critical size of the nucleating cell cluster.