Effects of stochasticity and division of labor in toxin production on two-strain bacterial competition in Escherichia coli

Resource limitation and population fluctuation drive spatiotemporal order in microbial communities

Microbial communities display complex spatiotemporal behaviors leading to spatially-structured and ordered organization driven by species interactions and environmental factors. Resource availability plays a pivotal role in shaping the dynamics of bacterial colonies. In this study, we delve into the intricate interplay between resource limitation and the emergent properties of a growing colony of two visually distinct bacterial strains having similar growth and mechanical properties. Employing an agent-based modeling and computer simulations, we analyze the resource-driven effect on segregation and sectoring, cell length regulation and nematic ordering within a growing colony. We introduce a dimensionless parameter referred to as the active layer thickness, derived from nutrient diffusion equations, indicating effective population participation due to local resource availability. Our results reveal that lower values of active layer thickness arising from decreased resource abundance lead to rougher colony fronts, fostering heightened population fluctuations within the colony and faster spatial genetic diversity loss. Our temporal analyses unveil the dynamics of mean cell length and fluctuations, showcasing how initial disturbances evolve as colonies are exposed to nutrients and subsequently settle. Furthermore, examining microscopic details, we find that lower resource levels yield diverse cell lengths and enhanced nematic ordering, driven by the increased prevalence of longer rod-shaped cells. Our investigation sheds light on the multifaceted relationship between resource constraints and bacterial colony dynamics, revealing insights into their spatiotemporal organization.

Spatial exclusion leads to "tug-of-war" ecological dynamics between competing species within microchannels

Competition is ubiquitous in microbial communities, shaping both their spatial and temporal structure and composition. Classical minimal models of competition, such as the Moran model, have been employed in ecology and evolutionary biology to understand the role of fixation and invasion in the maintenance of population diversity. Informed by recent experimental studies of cellular competition in confined spaces, we extend the Moran model to incorporate mechanical interactions between cells that divide within the limited space of a one-dimensional open microchannel. The model characterizes the skewed collective growth of the cells dividing within the channel, causing cells to be expelled at the channel ends. The results of this spatial exclusion model differ significantly from those of its classical well-mixed counterpart. The mean time to fixation of a species is greatly accelerated, scaling logarithmically, rather than algebraically, with the system size, and fixation/extinction probability sharply depends on the species' initial fractional abundance. By contrast, successful takeovers by invasive species, whether through mutation or immigration, are substantially less likely than in the Moran model. We also find that the spatial exclusion tends to attenuate the effects of fitness differences on the fixation times and probabilities. We find that these effects arise from the combination of the quasi-neutral "tug-of-war" diffusion dynamics of the inter-species boundary around an unstable equipoise point and the quasi-deterministic avalanche dynamics away from the fixed point. These results, which can be tested in microfluidic monolayer devices, have implications for the maintenance of species diversity in dense bacterial and cellular ecosystems where spatial exclusion is central to the competition, such as in organized biofilms or intestinal crypts.

Mutational meltdown of putative microbial altruists in Streptomyces coelicolor colonies

In colonies of the filamentous multicellular bacterium Streptomyces coelicolor, a subpopulation of cells arises that hyperproduces metabolically costly antibiotics, resulting in a division of labor that increases colony fitness. Because these cells contain large genomic deletions that cause massive reductions to individual fitness, their behavior is similar to altruistic worker castes in social insects or somatic cells in multicellular organisms. To understand these mutant cells’ reproductive and genomic fate after their emergence, we use experimental evolution by serially transferring populations via spore-to-spore transfer for 25 cycles, reflective of the natural mode of bottlenecked transmission for these spore-forming bacteria. We show that in contrast to wild-type cells, putatively altruistic mutant cells continue to decline in fitness during transfer while they lose more fragments from their chromosome ends. In addition, the base-substitution rate in mutants increases roughly 10-fold, possibly due to mutations in genes for DNA replication and repair. Ecological damage, caused by reduced sporulation, coupled with DNA damage due to point mutations and deletions, leads to an inevitable and irreversible type of mutational meltdown in these cells. Taken together, these results suggest the cells arising in the S. coelicolor division of labor are analogous to altruistic reproductively sterile castes of social insects.

In Streptomyces coelicolor, a subpopulation of cells can arise that produce metabolically costly antibiotics and a division of labor that maximizes colony fitness. This study uses experimental evolution to understand the reproductive and genomic fate of these mutant cells, showing that the arising altruistic cells are analogous to the reproductively sterile castes of social insects.

Modeling endonuclease colicin-like bacteriocin operons as 'genetic arms' in plasmid-genome conflicts

Plasmids are acellular propagating entities that depend on bacteria, as molecular parasites, for propagation. A 'tussle' between bacteria and plasmid ensues; bacteria for riddance of the plasmid and plasmid for persistence within a live host. Plasmid-maintenance systems such as endonuclease Colicin-Like Bacteriocins (CLBs) ensure plasmid propagation within the population; (i) the plasmid-cured cells are killed by the CLBs; (ii) damaged cells lyse and release the CLBs that eliminate the competitors, and (iii) the released plasmids invade new bacteria. Surprisingly, endonuclease CLB operons occur on bacterial genomes whose significance is unknown. Here, we study genetics, eco-evolutionary drive, and physiological relevance of genomic endonuclease CLB operons. We investigated plasmidic and genomic endonuclease CLB operons using sequence analyses from an eco-evolutionary perspective. We found 1266 genomic and plasmidic endonuclease CLB operons across 30 bacterial genera. Although 51% of the genomes harbor endonuclease CLB operons, the majority of the genomic endonuclease CLB operons lacked a functional lysis gene, suggesting the negative selection of lethal genes. The immunity gene of the endonuclease CLB operon protects the plasmid-cured host, eliminating the metabolic burden. We show mutual exclusivity of endonuclease CLB operons on genomes and plasmids. We propose an anti-addiction hypothesis for genomic endonuclease CLB operons. Using a stochastic hybrid agent-based model, we show that the endonuclease CLB operons on genomes confer an advantage to the host genome in terms of immunity to the toxin and elimination of plasmid burden. The conflict between bacterial genome and plasmids allows the emergence of 'genetic arms' such as CLB operons that regulate the ecological interplay of bacterial genomes and plasmids.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Coexistence of fluorescent Escherichia coli strains in millifluidic droplet reactors

Understanding competition and cooperation within microbiota is of high fundamental and clinical importance, helping to comprehend species' evolution and biodiversity. We co-encapsulated and cultured two isogenic Escherichia coli strains expressing blue (BFP) and yellow (YFP) fluorescent proteins into numerous emulsion droplets and quantified their growth by employing fluorescence measurements. To characterize and compare the bacterial growth kinetics and behavior in mono and co-culture, we compared the experimental observations with predictions from a simple growth model. Varying the initial ratio (R0) of both cell types injected, we observed a broad landscape from competition to cooperation between both strains in their confined microenvironments depending on start frequency: from a nearly symmetric situation at R0 = 1, up to the domination of one subpopulation when R0 ≫ 1 (or R0 ≪ 1). Due to competition between the strains, their doubling times and final biomass ratios (R1) continuously deviate from the monoculture behavior. The correlation map of the two strains' doubling times reveals that the R0 is one of the critical parameters affecting the competitive interaction between isogenic bacterial strains. Thanks to this strategy, different species of bacteria can be monitored simultaneously in real-time. Further advantages include high statistical output, unaffected bacteria growth, and long-time measurements in a well-mixed environment. We expect that the millifluidic droplet-based reactor can be utilized for practical clinical applications, such as bacterial antibiotic resistance and enzyme reaction kinetics studies.

Evolutionary Stabilization of Cooperative Toxin Production through a Bacterium-Plasmid-Phage Interplay

Colicins are toxins produced and released by Enterobacteriaceae to kill competitors in the gut. While group A colicins employ a division of labor strategy to liberate the toxin into the environment via colicin-specific lysis, group B colicin systems lack cognate lysis genes. In Salmonella enterica serovar Typhimurium (S. Tm), the group B colicin Ib (ColIb) is released by temperate phage-mediated bacteriolysis. Phage-mediated ColIb release promotes S. Tm fitness against competing Escherichia coli It remained unclear how prophage-mediated lysis is realized in a clonal population of ColIb producers and if prophages contribute to evolutionary stability of toxin release in S. Tm. Here, we show that prophage-mediated lysis occurs in an S. Tm subpopulation only, thereby introducing phenotypic heterogeneity to the system. We established a mathematical model to study the dynamic interplay of S. Tm, ColIb, and a temperate phage in the presence of a competing species. Using this model, we studied long-term evolution of phage lysis rates in a fluctuating infection scenario. This revealed that phage lysis evolves as bet-hedging strategy that maximizes phage spread, regardless of whether colicin is present or not. We conclude that the ColIb system, lacking its own lysis gene, is making use of the evolutionary stable phage strategy to be released. Prophage lysis genes are highly prevalent in nontyphoidal Salmonella genomes. This suggests that the release of ColIb by temperate phages is widespread. In conclusion, our findings shed new light on the evolution and ecology of group B colicin systems.IMPORTANCE Bacteria are excellent model organisms to study mechanisms of social evolution. The production of public goods, e.g., toxin release by cell lysis in clonal bacterial populations, is a frequently studied example of cooperative behavior. Here, we analyze evolutionary stabilization of toxin release by the enteric pathogen Salmonella The release of colicin Ib (ColIb), which is used by Salmonella to gain an edge against competing microbiota following infection, is coupled to bacterial lysis mediated by temperate phages. Here, we show that phage-dependent lysis and subsequent release of colicin and phage particles occurs only in part of the ColIb-expressing Salmonella population. This phenotypic heterogeneity in lysis, which represents an essential step in the temperate phage life cycle, has evolved as a bet-hedging strategy under fluctuating environments such as the gastrointestinal tract. Our findings suggest that prophages can thereby evolutionarily stabilize costly toxin release in bacterial populations.

Expansion, Exploitation and Extinction: Niche Construction in Ephemeral Landscapes

We aim to understand general consequences of niche construction on metapopulation dynamics in ephemeral landscapes. To this effect, a contact process-like stochastic spatial model is introduced where local populations colonize and go extinct on a dynamic landscape of habitable and destroyed patches. In contrast to previous models, where the extinction threshold is a consequence of available niche rendered by global rates of patch destruction/renewal, here we investigate how the metapopulation persists when they are the sole generators of their own niche. Niche construction is full-filled by localized populations through the transformation of destroyed patches in their neighborhood to viable habitat for future colonization. With this theoretical framework we are able to address the dual nature of niche construction by investigating the ephemerality of the landscape (destruction rate) and the continuum of population level strategies, where construction comes at a cost to colonization. Using mean field theory and Monte Carlo simulations of the model, we are able to quantify optimal population level strategies in a wide range of ephemeral landscapes. Interestingly, we observe qualitative differences at the extinction threshold between analytic and numeric results. Investigating this discrepancy further, we find that increasing niche construction neighborhood in the spatial model leads to two interrelated effects i) an increased rate in range expansion ii) a loss in resiliency and return of the discontinuous transition at the extinction threshold. Furthermore, in the discontinuous regime of the model, spatial clustering prior to a critical transition disappears. This is a significant finding as spatial clustering has been considered to be an early warning signal before ecosystems reach their 'tipping point'. In addition to maintaining stability, we find local niche construction strategies have an advantage when in scramble competition with an exploiter strategy because of their ability to monopolize the constructed niche due to spatial adjacency. As the niche construction neighborhood expands this advantage disappears and the exploiter strategy out-competes the niche constructor. In some cases the exploiter pushes the niche constructor to extinction, thus a tragedy of the commons ensues leading to 'ecological suicide' and a collapse of the niche.

A multidimensional perspective on microbial interactions

Beyond being simply positive or negative, beneficial or inhibitory, microbial interactions can involve a diverse set of mechanisms, dependencies and dynamical properties. These more nuanced features have been described in great detail for some specific types of interactions, (e.g. pairwise metabolic cross-feeding, quorum sensing or antibiotic killing), often with the use of quantitative measurements and insight derived from modeling. With a growing understanding of the composition and dynamics of complex microbial communities for human health and other applications, we face the challenge of integrating information about these different interactions into comprehensive quantitative frameworks. Here, we review the literature on a wide set of microbial interactions, and explore the potential value of a formal categorization based on multidimensional vectors of attributes. We propose that such an encoding can facilitate systematic, direct comparisons of interaction mechanisms and dependencies, and we discuss the relevance of an atlas of interactions for future modeling and rational design efforts.

Boosting functional response models for location, scale and shape with an application to bacterial competition

We extend generalized additive models for location, scale and shape (GAMLSS) to regression with functional response. This allows us to simultaneously model point-wise mean curves, variances and other distributional parameters of the response in dependence of various scalar and functional covariate effects. In addition, the scope of distributions is extended beyond exponential families. The model is fitted via gradient boosting, which offers inherent model selection and is shown to be suitable for both complex model structures and highly auto-correlated response curves. This enables us to analyse bacterial growth in Escherichia coli in a complex interaction scenario, fruitfully extending usual growth models.

Dynamics of ColicinE2 production and release determine the competitive success of a toxin-producing bacterial population

The release of toxins is one mechanism used by bacterial species to establish dominance over competitors, but how the dynamics of toxin expression determine the competitive success of a toxin-producing population is largely unknown. Here, we investigate how the expression dynamics of ColicinE2 - a toxic bacteriocin - affect competition between toxin-producing and toxin-sensitive strains of Escherichia coli. We demonstrate that, in addition to genetic modifications in the toxin expression system, alterations of the growth medium can be used to modulate the timing of toxin production and the amount of toxin released. Thus cells that release the toxin at later times can accumulate more colicin. In experiments, we found that delaying toxin release does not significantly alter competition outcome. However, our theoretical analysis allowed us to assess the relative contributions of release time and toxin level to the competitive success of the producer strain, that might counteract each other in experiments. The results reveal that the importance of delaying toxin release lies in increasing the toxin amount. This is a more effective strategy for the toxin-producing strain than prompt discharge of the colicin. In summary, our study shows how the toxin release dynamics influence the competitive success of the toxin-producing bacterial population.

High-resolution sweep metagenomics using fast probabilistic inference

Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.

Gene expression noise in a complex artificial toxin expression system

Gene expression is an intrinsically stochastic process. Fluctuations in transcription and translation lead to cell-to-cell variations in mRNA and protein levels affecting cellular function and cell fate. Here, using fluorescence time-lapse microscopy, we quantify noise dynamics in an artificial operon in Escherichia coli, which is based on the native operon of ColicinE2, a toxin. In the natural system, toxin expression is controlled by a complex regulatory network; upon induction of the bacterial SOS response, ColicinE2 is produced (cea gene) and released (cel gene) by cell lysis. Using this ColicinE2-based operon, we demonstrate that upon induction of the SOS response noise of cells expressing the operon is significantly lower for the (mainly) transcriptionally regulated gene cea compared to the additionally post-transcriptionally regulated gene cel. Likewise, we find that mutations affecting the transcriptional regulation by the repressor LexA do not significantly alter the population noise, whereas specific mutations to post-transcriptionally regulating units, strongly influence noise levels of both genes. Furthermore, our data indicate that global factors, such as the plasmid copy number of the operon encoding plasmid, affect gene expression noise of the entire operon. Taken together, our results provide insights on how noise in a native toxin-producing operon is controlled and underline the importance of post-transcriptional regulation for noise control in this system.

Microbial Interaction Network Inference in Microfluidic Droplets

Microbial interactions are major drivers of microbial community dynamics and functions but remain challenging to identify because of limitations in parallel culturing and absolute abundance quantification of community members across environments and replicates. To this end, we developed Microbial Interaction Network Inference in microdroplets (MINI-Drop). Fluorescence microscopy coupled to computer vision techniques were used to rapidly determine the absolute abundance of each strain in hundreds to thousands of droplets per condition. We showed that MINI-Drop could accurately infer pairwise and higher-order interactions in synthetic consortia. We developed a stochastic model of community assembly to provide insight into the heterogeneity in community states across droplets. Finally, we elucidated the complex web of interactions linking antibiotics and different species in a synthetic consortium. In sum, we demonstrated a robust and generalizable method to infer microbial interaction networks by random encapsulation of sub-communities into microfluidic droplets.

Evolutionary Stability of Salmonella Competition with the Gut Microbiota: How the Environment Fosters Heterogeneity in Exploitative and Interference Competition

Following ingestion, gastrointestinal pathogens compete against the gastrointestinal microbiota and overcome host immune defenses in order to cause infections. Besides employing direct killing mechanisms, the commensal microbiota occupies metabolic niches to outcompete invading pathogens. Salmonella enterica serovar Typhimurium (S. Typhimurium) uses several strategies to successfully colonize the gut and establish infection, of which an increasing number is based on phenotypic heterogeneity within the S. Typhimurium population. The utilization of myo-inositol (MI) and the production of colicin confer a selective advantage over the microbiota in terms of exploitative and interference competition, respectively. In this review, we summarize the genetic basis underlying bistability of MI catabolism and colicin production. As demonstrated by single-cell analyses, a stochastic switch in the expression of the genes responsible for colicin production and MI degradation constitutes the heterogeneity of the two phenotypes. Both genetic systems are tightly regulated to avoid their expression under non-appropriate conditions and possible detrimental effects on bacterial fitness. Moreover, evolutionary mechanisms underlying formation and stability of these phenotypes in S. Typhimurium are discussed. We propose that both MI catabolism and colicin production create a bet-hedging strategy, which provides an adaptive benefit for S. Typhimurium in the fluctuating environment of the mammalian gut.

Locality of interactions in three-strain bacterial competition in E. coli

The population dynamics that determine the composition and stability of ecosystems ultimately emerge from interactions between individual organisms. One well-studied system is the three-strain E. coli interaction of a heterogeneously toxin-producing C strain that interacts with a toxin-sensitive S and a toxin-resistant R strain. Here, we employ a multi-scale fluorescence microscopy approach, that has been proven useful in identifying previously unknown or underestimated stochastic effects in C-S competition. This approach allows us to investigate the microscopic interaction of the R strain and to quantify the role of stochastic effects in the spatially structured C-R-S interaction. We show that the early colony patterning at 12 h and at small length scales (near single cell level) is characterized by a number of microscopic variables (the number of C and R cell clusters and the area occupied by S) and is subject to random processes in positioning and toxin production. Then, in a second competition phase, mainly deterministic processes such as bacterial growth and global toxin action determine the following population dynamics. Consequently, together with environmental factors, the microscopic variables were predictive of the competition outcome. However, interactions of neighboring R and C clusters could amplify local variations. If R clusters originated near a C cell cluster, R could profit from the toxin produced by C without bearing the cost of production-a mechanism called cheating. By combining information from the micro- and macro-scale dynamics, we can estimate the distance at which the cheating interaction significantly changes to be in the order of 250 µm. In summary, after an initial phase influenced by stochastic patterning, largely deterministic growth dynamics follow, which are additionally affected by local interactions of neighboring clusters. As such, the results underline the importance of stochasticity and local effects in the context of ecological interactions.

Complex microbial systems across different levels of description

Complex biological systems offer a variety of interesting phenomena at the different physical scales. With increasing abstraction, details of the microscopic scales can often be extrapolated to average or typical macroscopic properties. However, emergent properties and cross-scale interactions can impede naïve abstractions and necessitate comprehensive investigations of these complex systems. In this review paper, we focus on microbial communities, and first, summarize a general hierarchy of relevant scales and description levels to understand these complex systems: (1) genetic networks, (2) single cells, (3) populations, and (4) emergent multi-cellular properties. Second, we employ two illustrating examples, microbial competition and biofilm formation, to elucidate how cross-scale interactions and emergent properties enrich the observed multi-cellular behavior in these systems. Finally, we conclude with pointing out the necessity of multi-scale investigations to understand complex biological systems and discuss recent investigations.

What Does Not Kill You Makes You Stronger

Colicin production is an extreme form of labor division; cells lyse after making the toxin! Stochastic phenotype switching allows producers to outcompete sensitive strains since colicin release frees up vacancy. If patch dynamics does not kill you, it stimulates adaptation to a dynamic habitat landscape which selects for rapid dispersal.