Evolution of species interactions in a biofilm community

Low-abundant species facilitates specific spatial organization that promotes multispecies biofilm formation

Microorganisms frequently co-exist in matrix-embedded multispecies biofilms. Within biofilms, interspecies interactions influence the spatial organization of member species, which likely play an important role in shaping the development, structure and function of these communities. Here, a reproducible four-species biofilm, composed of Stenotrophomonas rhizophila, Xanthomonas retroflexus, Microbacterium oxydans and Paenibacillus amylolyticus, was established to study the importance of individual species spatial organization during multispecies biofilm development. We found that the growth of species that are poor biofilm formers, M. oxydans and P. amylolyticus, were highly enhanced when residing in the four-species biofilm. Interestingly, the presence of the low-abundant M. oxydans (0.5% of biomass volume) was observed to trigger changes in the composition of the four-species community. The other three species were crucially needed for the successful inclusion of M. oxydans in the four-species biofilm, where X. retroflexus was consistently positioned in the top layer of the mature four-species biofilm. These findings suggest that low abundance key species can significantly impact the spatial organization and hereby stabilize the function and composition of complex microbiomes.

Effect of Antimicrobial and Physical Treatments on Growth of Multispecies Staphylococcal Biofilms

The prevalence and structure of Staphylococcus aureus and Staphylococcus epidermidis within multispecies biofilms were found to depend sensitively on physical environment and antibiotic dosage. Although these species commonly infect similar sites, such as orthopedic implants, little is known about their behavior in multispecies communities, particularly in response to treatment. This research establishes that S. aureus is much more prevalent than S. epidermidis when simultaneously seeded and grown under unstressed conditions (pH 7, 37°C) in both laboratory and clinical strains. In multispecies communities, S. epidermidis is capable of growing a more confluent biofilm when the addition of S. aureus is delayed 4 to 6 h during 18 h of growth. Different vancomycin dosages generate various behaviors: S. epidermidis is more prevalent at a dose of 1.0 μg/ml vancomycin, but reduced growth of both species occurs at 1.9 μg/ml vancomycin. This variability is consistent with the different MICs of S. aureus and S. epidermidis Growth at higher temperature (45°C) results in an environment where S. aureus forms porous biofilms. This porosity allows S. epidermidis to colonize more of the surface, resulting in detectable S. epidermidis biomass. Variations in pH result in increased prevalence of S. epidermidis at low pH (pH 5 and 6), while S. aureus remains dominant at high pH (pH 8 and 9). This work establishes the structural variability of multispecies staphylococcal biofilms as they undergo physical and antimicrobial treatments. It provides a basis for understanding the structure of these communities at infection sites and how treatments disrupt their multispecies behaviors.IMPORTANCEStaphylococcus aureus and Staphylococcus epidermidis are two species of bacteria that are commonly responsible for biofilm infections on medical devices. Biofilms are structured communities of bacteria surrounded by polysaccharides, proteins, and DNA; bacteria are more resistant to antimicrobials as part of a biofilm than as individual cells. This work investigates the structure and prevalence of these two organisms when grown together in multispecies biofilms and shows shifts in the behavior of the polymicrobial community when grown in various concentrations of vancomycin (an antibiotic commonly used to treat staphylococcal infections), in a high-temperature environment (a condition previously shown to lead to cell disruption and death), and at low and high pH (a change that has been previously shown to soften the mechanical properties of staphylococcal biofilms). These shifts in community structure demonstrate the effect such treatments may have on multispecies staphylococcal infections.

Biofilms: an emergent form of bacterial life

Bacterial biofilms are formed by communities that are embedded in a self-produced matrix of extracellular polymeric substances (EPS). Importantly, bacteria in biofilms exhibit a set of 'emergent properties' that differ substantially from free-living bacterial cells. In this Review, we consider the fundamental role of the biofilm matrix in establishing the emergent properties of biofilms, describing how the characteristic features of biofilms - such as social cooperation, resource capture and enhanced survival of exposure to antimicrobials - all rely on the structural and functional properties of the matrix. Finally, we highlight the value of an ecological perspective in the study of the emergent properties of biofilms, which enables an appreciation of the ecological success of biofilms as habitat formers and, more generally, as a bacterial lifestyle.

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

In phenotypically heterogeneous microbial populations, the decision to adopt one or another phenotype is often stochastically regulated. However, how this stochasticity affects interactions between competing microbes in mixed communities is difficult to assess. One example of such an interaction system is the competition of an Escherichia coli strain C, which performs division of labor between reproducers and self-sacrificing toxin producers, with a toxin-sensitive strain S. The decision between reproduction or toxin production within a single C cell is inherently stochastic. Here, combining experimental and theoretical approaches, we demonstrate that this stochasticity in the initial phase of colony formation is the crucial determinant for the competition outcome. In the initial phase (t < 12h), stochasticity influences the formation of viable C clusters at the colony edge. In the subsequent phase, the effective fitness differences (set primarily by the degree of division of labor in the C strain population) dictate the deterministic population dynamics and consequently competition outcome. In particular, we observe that competitive success of the C strain is only found if (i) a C edge cluster has formed at the end of the initial competition phase and (ii) the beneficial and detrimental effects of toxin production are balanced, which is the case at intermediate toxin producer fractions. Our findings highlight the importance of stochastic processes during the initial phase of colony formation, which might be highly relevant for other microbial community interactions in which the random choice between phenotypes can have long-lasting consequences for community fate.

Comparative genomics and metabolomics analyses of the adaptation mechanism in Ketogulonicigenium vulgare-Bacillus thuringiensis consortium

Adaptive evolution by serial subcultivation of co-cultured Bacillus thuringiensis and Ketogulonicigenium vulgare significantly enhanced the productivity of 2-keto-L-gulonic acid in two-step vitamin C production. The adaptation mechanism in K. vulgare-B. thuringiensis consortium was investigated in this study based on comparative genomics and metabolomics studies. It was found that the growth, anti-oxidation, transcription and regulation were significantly enhanced in the adapted consortium. The mutation of the genes, which encode amidohydrolase in adapted K. vulgare (K150) and amino acid permease in adapted B. thuringiensis (B150), resulted in the increase of some amino acids levels in each species, and further enhanced the metabolic exchange and growth ability of the two species. Besides, the mutation of the gene encoding spore germination protein enhanced the metabolic levels of tricarboxylic acid cycle, and decreased the sporulation in B150, which induced its growth. The mutation of the genes, which encode NADPH nitroreductase in K150 and NADPH-dependent FMN reductase in B150, may enhance the ability of anti-oxidation. Overall, the long-term adaptation of K. vulgare and B. thuringiensis influenced the global regulation and made them more inseparable in metabolite exchange. Our work will provide ideas for the molecular design and optimization in microbial consortium.

Draft Genome Sequence of Acinetobacter johnsonii C6, an Environmental Isolate Engaging in Interspecific Metabolic Interactions

Acinetobacter johnsonii C6 originates from creosote-polluted groundwater and performs ecological and evolutionary interactions with Pseudomonas putida in biofilms. The draft genome of A. johnsonii C6 is 3.7 Mbp and was shaped by mobile genetic elements. It reveals genes facilitating the biodegradation of aromatic hydrocarbons and resistance to antimicrobials and metals.

Coevolution with bacteria drives the evolution of aerobic fermentation in Lachancea kluyveri

The Crabtree positive yeasts, such as Saccharomyces cerevisiae, prefer fermentation to respiration, even under fully aerobic conditions. The selective pressures that drove the evolution of this trait remain controversial because of the low ATP yield of fermentation compared to respiration. Here we propagate experimental populations of the weak-Crabtree yeast Lachancea kluyveri, in competitive co-culture with bacteria. We find that L. kluyveri adapts by producing quantities of ethanol lethal to bacteria and evolves several of the defining characteristics of Crabtree positive yeasts. We use precise quantitative analysis to show that the rate advantage of fermentation over aerobic respiration is insufficient to provide an overall growth advantage. Thus, the rapid consumption of glucose and the utilization of ethanol are essential for the success of the aerobic fermentation strategy. These results corroborate that selection derived from competition with bacteria could have provided the impetus for the evolution of the Crabtree positive trait.

Cooperation in carbon source degradation shapes spatial self-organization of microbial consortia on hydrated surfaces

Mounting evidence suggests that natural microbial communities exhibit a high level of spatial organization at the micrometric scale that facilitate ecological interactions and support biogeochemical cycles. Microbial patterns are difficult to study definitively in natural environments due to complex biodiversity, observability and variable physicochemical factors. Here, we examine how trophic dependencies give rise to self-organized spatial patterns of a well-defined bacterial consortium grown on hydrated surfaces. The model consortium consisted of two Pseudomonas putida mutant strains that can fully degrade the aromatic hydrocarbon toluene. We demonstrated that obligate cooperation in toluene degradation (cooperative mutualism) favored convergence of 1:1 partner ratio and strong intermixing at the microscale (10–100 μm). In contrast, competition for benzoate, a compound degraded independently by both strains, led to distinct segregation patterns. Emergence of a persistent spatial pattern has been predicted for surface attached microbial activity in liquid films that mediate diffusive exchanges while permitting limited cell movement (colony expansion). This study of a simple microbial consortium offers mechanistic glimpses into the rules governing the assembly and functioning of complex sessile communities, and points to general principles of spatial organization with potential applications for natural and engineered microbial systems.

Cell morphology drives spatial patterning in microbial communities

The clearest phenotypic characteristic of microbial cells is their shape, but we do not understand how cell shape affects the dense communities, known as biofilms, where many microbes live. Here, we use individual-based modeling to systematically vary cell shape and study its impact in simulated communities. We compete cells with different cell morphologies under a range of conditions and ask how shape affects the patterning and evolutionary fitness of cells within a community. Our models predict that cell shape will strongly influence the fate of a cell lineage: we describe a mechanism through which coccal (round) cells rise to the upper surface of a community, leading to a strong spatial structuring that can be critical for fitness. We test our predictions experimentally using strains of Escherichia coli that grow at a similar rate but differ in cell shape due to single amino acid changes in the actin homolog MreB. As predicted by our model, cell types strongly sort by shape, with round cells at the top of the colony and rod cells dominating the basal surface and edges. Our work suggests that cell morphology has a strong impact within microbial communities and may offer new ways to engineer the structure of synthetic communities.

Eco-evolutionary Model of Rapid Phenotypic Diversification in Species-Rich Communities

Evolutionary and ecosystem dynamics are often treated as different processes –operating at separate timescales– even if evidence reveals that rapid evolutionary changes can feed back into ecological interactions. A recent long-term field experiment has explicitly shown that communities of competing plant species can experience very fast phenotypic diversification, and that this gives rise to enhanced complementarity in resource exploitation and to enlarged ecosystem-level productivity. Here, we build on progress made in recent years in the integration of eco-evolutionary dynamics, and present a computational approach aimed at describing these empirical findings in detail. In particular we model a community of organisms of different but similar species evolving in time through mechanisms of birth, competition, sexual reproduction, descent with modification, and death. Based on simple rules, this model provides a rationalization for the emergence of rapid phenotypic diversification in species-rich communities. Furthermore, it also leads to non-trivial predictions about long-term phenotypic change and ecological interactions. Our results illustrate that the presence of highly specialized, non-competing species leads to very stable communities and reveals that phenotypically equivalent species occupying the same niche may emerge and coexist for very long times. Thus, the framework presented here provides a simple approach –complementing existing theories, but specifically devised to account for the specificities of the recent empirical findings for plant communities– to explain the collective emergence of diversification at a community level, and paves the way to further scrutinize the intimate entanglement of ecological and evolutionary processes, especially in species-rich communities.

Population ecology and evolutionary biology have been traditionally studied as separate disciplines, even if feedbacks between community and evolutionary processes are known to exist, having been empirically characterized in recent years in different types of communities (from microbes to plants and vertebrates), and theoretically analyzed with novel and powerful mathematical tools. Recent long-term field experiments with plants have proven that rapid co-evolution and diversification of species traits results in an overall enhancement of the ecosystem productivity, with important consequences for agriculture and conservation. Here, we propose a relatively simple computational eco-evolutionary model specifically devised to describe rapid phenotypic diversification in this type of species-rich communities. Our model captures the main phenomenology observed experimentally, and it also makes non-trivial predictions for long term phenotypic change and ecological interactions, such as the stable coexistence of highly specialized species or the possible emergence of phenotypically equivalent species occupying the same niche. Finally, the model is easily generalizable to analyze different eco-evolutionary problems within a relatively simple and unified computational framework.

Metabolic Complementation in Bacterial Communities: Necessary Conditions and Optimality

Bacterial communities may display metabolic complementation, in which different members of the association partially contribute to the same biosynthetic pathway. In this way, the end product of the pathway is synthesized by the community as a whole. However, the emergence and the benefits of such complementation are poorly understood. Herein, we present a simple model to analyze the metabolic interactions among bacteria, including the host in the case of endosymbiotic bacteria. The model considers two cell populations, with both cell types encoding for the same linear biosynthetic pathway. We have found that, for metabolic complementation to emerge as an optimal strategy, both product inhibition and large permeabilities are needed. In the light of these results, we then consider the patterns found in the case of tryptophan biosynthesis in the endosymbiont consortium hosted by the aphid Cinara cedri. Using in-silico computed physicochemical properties of metabolites of this and other biosynthetic pathways, we verified that the splitting point of the pathway corresponds to the most permeable intermediate.

Diversification of Gene Expression during Formation of Static Submerged Biofilms by Escherichia coli

Many bacteria primarily exist in nature as structured multicellular communities, so called biofilms. Biofilm formation is a highly regulated process that includes the transition from the motile planktonic to sessile biofilm lifestyle. Cellular differentiation within a biofilm is a commonly accepted concept but it remains largely unclear when, where and how exactly such differentiation arises. Here we used fluorescent transcriptional reporters to quantitatively analyze spatio-temporal expression patterns of several groups of genes during the formation of submerged Escherichia coli biofilms in an open static system. We first confirm that formation of such submerged biofilms as well as pellicles at the liquid-air interface requires the major matrix component, curli, and flagella-mediated motility. We further demonstrate that in this system, diversification of gene expression leads to emergence of at least three distinct subpopulations of E. coli, which differ in their levels of curli and flagella expression, and in the activity of the stationary phase sigma factor σS. Our study reveals mutually exclusive expression of curli fibers and flagella at the single cell level, with high curli levels being confined to dense cell aggregates/microcolonies and flagella expression showing an opposite expression pattern. Interestingly, despite the known σS-dependence of curli induction, there was only a partial correlation between the σS activity and curli expression, with subpopulations of cells having high σS activity but low curli expression and vice versa. Finally, consistent with different physiology of the observed subpopulations, we show striking differences between the growth rates of cells within and outside of aggregates.

A crucial role for spatial distribution in bacterial quorum sensing

Quorum sensing (QS) is a process that enables bacteria to communicate using secreted signaling molecules, and then makes a population of bacteria to regulate gene expression collectively and control behavior on a community-wide scale. Theoretical studies of efficiency sensing have suggested that both mass-transfer performance in the local environment and the spatial distribution of cells are key factors affecting QS. Here, an experimental model based on hydrogel microcapsules with a three-dimensional structure was established to investigate the influence of the spatial distribution of cells on bacterial QS. Vibrio harveyi cells formed different spatial distributions in the microcapsules, i.e., they formed cell aggregates with different structures and sizes. The cell aggregates displayed stronger QS than did unaggregated cells even when equal numbers of cells were present. Large aggregates (LA) of cells, with a size of approximately 25 μm, restricted many more autoinducers (AIs) than did small aggregates (SA), with a size of approximately 10 μm, thus demonstrating that aggregate size significantly affects QS. These findings provide a powerful demonstration of the fact that the spatial distribution of cells plays a crucial role in bacterial QS.

Long-term performance of rapid oxidation of arsenite in simulated groundwater using a population of arsenite-oxidizing microorganisms in a bioreactor

A population of arsenite-oxidizing microorganisms enriched from the tailing of the Shimen realgar mine was used to generate biofilms on the surfaces of perlites. This bioreactor is able to completely oxidize 1100 μg/L As(III) dissolved in simulated groundwater into As(V) within 10 min; after 140 days of operation, approximately 20 min were required to completely oxidize the same concentration of As(III). Analysis for the 16S rRNA genes of the microbial community showed that Bacteroidetes and Proteobacteria are dominant in the reactor. Six different bacterial strains were randomly isolated from the reactor. Function and gene analysis indicated that all the isolates possess arsenite-oxidizing activity, and five of them are chemoautotrophic. Further analysis showed that a large diversity of AioAs and two types of RuBisCOs are present in the microbial community. This suggests that many chemoautotrophic arsenite-oxidizing microorganisms were responsible for quick oxidation of arsenite in the reactor. We also found that the reactor is easily regenerated and its number is readily expanded. To the best of our knowledge, the arsenite-oxidizing efficiency, which was expressed as the minimum time for complete oxidization of a certain concentration of As(III) under a single operation, of this bioreactor is the highest among the described bioreactors; it is also the most stable, economic and environment-friendly.

Laboratory Evolution of Microbial Interactions in Bacterial Biofilms

Microbial adaptation is conspicuous in essentially every environment, but the mechanisms of adaptive evolution are poorly understood. Studying evolution in the laboratory under controlled conditions can be a tractable approach, particularly when new, discernible phenotypes evolve rapidly. This is especially the case in the spatially structured environments of biofilms, which promote the occurrence and stability of new, heritable phenotypes. Further, diversity in biofilms can give rise to nascent social interactions among coexisting mutants and enable the study of the emerging field of sociomicrobiology. Here, we review findings from laboratory evolution experiments with either Pseudomonas fluorescens or Burkholderia cenocepacia in spatially structured environments that promote biofilm formation. In both systems, ecotypes with overlapping niches evolve and produce competitive or facilitative interactions that lead to novel community attributes, demonstrating the parallelism of adaptive processes captured in the lab.

Interspecific Bacterial Interactions are Reflected in Multispecies Biofilm Spatial Organization

Interspecies interactions are essential for the persistence and development of any kind of complex community, and microbial biofilms are no exception. Multispecies biofilms are structured and spatially defined communities that have received much attention due to their omnipresence in natural environments. Species residing in these complex bacterial communities usually interact both intra- and interspecifically. Such interactions are considered to not only be fundamental in shaping overall biomass and the spatial distribution of cells residing in multispecies biofilms, but also to result in coordinated regulation of gene expression in the different species present. These communal interactions often lead to emergent properties in biofilms, such as enhanced tolerance against antibiotics, host immune responses, and other stresses, which have been shown to provide benefits to all biofilm members not only the enabling sub-populations. However, the specific molecular mechanisms of cellular processes affecting spatial organization, and vice versa, are poorly understood and very complex to unravel. Therefore, detailed description of the spatial organization of individual bacterial cells in multispecies communities can be an alternative strategy to reveal the nature of interspecies interactions of constituent species. Closing the gap between visual observation and biological processes may become crucial for resolving biofilm related problems, which is of utmost importance to environmental, industrial, and clinical implications. This review briefly presents the state of the art of studying interspecies interactions and spatial organization of multispecies communities, aiming to support theoretical and practical arguments for further advancement of this field.

New Is Old, and Old Is New: Recent Advances in Antibiotic-Based, Antibiotic-Free and Ethnomedical Treatments against Methicillin-Resistant Staphylococcus aureus Wound Infections

Staphylococcus aureus is the most common pathogen of wound infections. Thus far, methicillin-resistant S. aureus (MRSA) has become the major causative agent in wound infections, especially for nosocomial infections. MRSA infections are seldom eradicated by routine antimicrobial therapies. More concerning, some strains have become resistant to the newest antibiotics of last resort. Furthermore, horizontal transfer of a polymyxin resistance gene, mcr-1, has been identified in Enterobacteriaceae, by which resistance to the last group of antibiotics will likely spread rapidly. The worst-case scenario, "a return to the pre-antibiotic era", is likely in sight. A perpetual goal for antibiotic research is the discovery of an antibiotic that lacks resistance potential, such as the recent discovery of teixobactin. However, when considering the issue from an ecological and evolutionary standpoint, it is evident that it is insufficient to solve the antibiotic dilemma through the use of antibiotics themselves. In this review, we summarized recent advances in antibiotic-based, antibiotic-free and ethnomedical treatments against MRSA wound infections to identify new clues to solve the antibiotic dilemma. One potential solution is to use ethnomedical drugs topically. Some ethnomedical drugs have been demonstrated to be effective antimicrobials against MRSA. A decline in antibiotic resistance can therefore be expected, as has been demonstrated when antibiotic-free treatments were used to limit the use of antibiotics. It is also anticipated that these drugs will have low resistance potential, although there is only minimal evidence to support this claim to date. More clinical trials and animal tests should be conducted on this topic.

Building the microbiome in health and disease: niche construction and social conflict in bacteria

Microbes collectively shape their environment in remarkable ways via the products of their metabolism. The diverse environmental impacts of macro-organisms have been collated and reviewed under the banner of ‘niche construction’. Here, we identify and review a series of broad and overlapping classes of bacterial niche construction, ranging from biofilm production to detoxification or release of toxins, enzymes, metabolites and viruses, and review their role in shaping microbiome composition, human health and disease. Some bacterial niche-constructing traits can be seen as extended phenotypes, where individuals actively tailor their environment to their benefit (and potentially to the benefit of others, generating social dilemmas). Other modifications can be viewed as non-adaptive by-products from a producer perspective, yet they may lead to remarkable within-host environmental changes. We illustrate how social evolution and niche construction perspectives offer complementary insights into the dynamics and consequences of these traits across distinct timescales. This review highlights that by understanding the coupled bacterial and biochemical dynamics in human health and disease we can better manage host health.

Mechanisms of Competition in Biofilm Communities

Bacterial biofilms are dense and often mixed-species surface-attached communities in which bacteria coexist and compete for limited space and nutrients. Here we present the different antagonistic interactions described in biofilm environments and their underlying molecular mechanisms, along with ecological and evolutionary insights as to how competitive interactions arise and are maintained within biofilms.

Lineage Tracking for Probing Heritable Phenotypes at Single-Cell Resolution

Determining the phenotype and genotype of single cells is central to understand microbial evolution. DNA sequencing technologies allow the detection of mutants at high resolution, but similar approaches for phenotypic analyses are still lacking. We show that a drop-based millifluidic system enables the detection of heritable phenotypic changes in evolving bacterial populations. At time intervals, cells were sampled and individually compartmentalized in 100 nL drops. Growth through 15 generations was monitored using a fluorescent protein reporter. Amplification of heritable changes–via growth–over multiple generations yields phenotypically distinct clusters reflecting variation relevant for evolution. To demonstrate the utility of this approach, we follow the evolution of Escherichia coli populations during 30 days of starvation. Phenotypic diversity was observed to rapidly increase upon starvation with the emergence of heritable phenotypes. Mutations corresponding to each phenotypic class were identified by DNA sequencing. This scalable lineage-tracking technology opens the door to large-scale phenotyping methods with special utility for microbiology and microbial population biology.

New Technologies for Studying Biofilms

Bacteria have traditionally been studied as single-cell organisms. In laboratory settings, aerobic bacteria are usually cultured in aerated flasks, where the cells are considered essentially homogenous. However, in many natural environments, bacteria and other microorganisms grow in mixed communities, often associated with surfaces. Biofilms are comprised of surface-associated microorganisms, their extracellular matrix material, and environmental chemicals that have adsorbed to the bacteria or their matrix material. While this definition of a biofilm is fairly simple, biofilms are complex and dynamic. Our understanding of the activities of individual biofilm cells and whole biofilm systems has developed rapidly, due in part to advances in molecular, analytical, and imaging tools and the miniaturization of tools designed to characterize biofilms at the enzyme level, cellular level, and systems level.

An algorithm for designing minimal microbial communities with desired metabolic capacities

MOTIVATION

Recent efforts to manipulate various microbial communities, such as fecal microbiota transplant and bioreactor systems' optimization, suggest a promising route for microbial community engineering with numerous medical, environmental and industrial applications. However, such applications are currently restricted in scale and often rely on mimicking or enhancing natural communities, calling for the development of tools for designing synthetic communities with specific, tailored, desired metabolic capacities.

RESULTS

Here, we present a first step toward this goal, introducing a novel algorithm for identifying minimal sets of microbial species that collectively provide the enzymatic capacity required to synthesize a set of desired target product metabolites from a predefined set of available substrates. Our method integrates a graph theoretic representation of network flow with the set cover problem in an integer linear programming (ILP) framework to simultaneously identify possible metabolic paths from substrates to products while minimizing the number of species required to catalyze these metabolic reactions. We apply our algorithm to successfully identify minimal communities both in a set of simple toy problems and in more complex, realistic settings, and to investigate metabolic capacities in the gut microbiome. Our framework adds to the growing toolset for supporting informed microbial community engineering and for ultimately realizing the full potential of such engineering efforts.

AVAILABILITY AND IMPLEMENTATION

The algorithm source code, compilation, usage instructions and examples are available under a non-commercial research use only license at https://github.com/borenstein-lab/CoMiDA CONTACT: elbo@uw.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Experimental evolution in biofilm populations

Biofilms are a major form of microbial life in which cells form dense surface associated communities that can persist for many generations. The long-life of biofilm communities means that they can be strongly shaped by evolutionary processes. Here, we review the experimental study of evolution in biofilm communities. We first provide an overview of the different experimental models used to study biofilm evolution and their associated advantages and disadvantages. We then illustrate the vast amount of diversification observed during biofilm evolution, and we discuss (i) potential ecological and evolutionary processes behind the observed diversification, (ii) recent insights into the genetics of adaptive diversification, (iii) the striking degree of parallelism between evolution experiments and real-life biofilms and (iv) potential consequences of diversification. In the second part, we discuss the insights provided by evolution experiments in how biofilm growth and structure can promote cooperative phenotypes. Overall, our analysis points to an important role of biofilm diversification and cooperation in bacterial survival and productivity. Deeper understanding of both processes is of key importance to design improved antimicrobial strategies and diagnostic techniques.

This review paper provides an overview of (i) the different experimental models used to study biofilm evolution, (ii) the vast amount of diversification observed during biofilm evolution (including potential causes and consequences) and (iii) recent insights in how growth in biofilms can lead to the evolution of cooperative phenotypes.

Large-scale biofilm cultivation of Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 for physiologic studies and drug discovery

Microbial biofilms are mainly studied due to detrimental effects on human health but they are also well established in industrial biotechnology for the production of chemicals. Moreover, biofilm can be considered as a source of novel drugs since the conditions prevailing within biofilm can allow the production of specific metabolites. Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 when grown in biofilm condition produces an anti-biofilm molecule able to inhibit the biofilm of the opportunistic pathogen Staphylococcus epidermidis. In this paper we set up a P. haloplanktis TAC125 biofilm cultivation methodology in automatic bioreactor. The biofilm cultivation was designated to obtain two goals: (1) the scale up of cell-free supernatant production in an amount necessary for the anti-biofilm molecule/s purification; (2) the recovery of P. haloplanktis TAC125 cells grown in biofilm for physiological studies. We set up a fluidized-bed reactor fermentation in which floating polystyrene supports were homogeneously mixed, exposing an optimal air-liquid interface to let bacterium biofilm formation. The proposed methodology allowed a large-scale production of anti-biofilm molecule and paved the way to study differences between P. haloplanktis TAC125 cells grown in biofilm and in planktonic conditions. In particular, the modifications occurring in the lipopolysaccharide of cells grown in biofilm were investigated.

Tools for the Microbiome: Nano and Beyond

The microbiome presents great opportunities for understanding and improving the world around us and elucidating the interactions that compose it. The microbiome also poses tremendous challenges for mapping and manipulating the entangled networks of interactions among myriad diverse organisms. Here, we describe the opportunities, technical needs, and potential approaches to address these challenges, based on recent and upcoming advances in measurement and control at the nanoscale and beyond. These technical needs will provide the basis for advancing the largely descriptive studies of the microbiome to the theoretical and mechanistic understandings that will underpin the discipline of microbiome engineering. We anticipate that the new tools and methods developed will also be more broadly useful in environmental monitoring, medicine, forensics, and other areas.

Population-Dynamic Modeling of Bacterial Horizontal Gene Transfer by Natural Transformation

Natural transformation is a major mechanism of horizontal gene transfer (HGT) and plays an essential role in bacterial adaptation, evolution, and speciation. Although its molecular underpinnings have been increasingly revealed, natural transformation is not well characterized in terms of its quantitative ecological roles. Here, by using Neisseria gonorrhoeae as an example, we developed a population-dynamic model for natural transformation and analyzed its dynamic characteristics with nonlinear tools and simulations. Our study showed that bacteria capable of natural transformation can display distinct population behaviors ranging from extinction to coexistence and to bistability, depending on their HGT rate and selection coefficient. With the model, we also illustrated the roles of environmental DNA sources-active secretion and passive release-in impacting population dynamics. Additionally, by constructing and utilizing a stochastic version of the model, we examined how noise shapes the steady and dynamic behaviors of the system. Notably, we found that distinct waiting time statistics for HGT events, namely a power-law distribution, an exponential distribution, and a mix of the both, are associated with the dynamics in the regimes of extinction, coexistence, and bistability accordingly. This work offers a quantitative illustration of natural transformation by revealing its complex population dynamics and associated characteristics, therefore advancing our ecological understanding of natural transformation as well as HGT in general.

Privatization of cooperative benefits stabilizes mutualistic cross-feeding interactions in spatially structured environments

Metabolic cross-feeding interactions are ubiquitous in natural microbial communities. However, it remains generally unclear whether the production and exchange of metabolites incurs fitness costs to the producing cells and if so, which ecological mechanisms can facilitate a cooperative exchange of metabolites among unrelated individuals. We hypothesized that positive assortment within structured environments can maintain mutualistic cross-feeding. To test this, we engineered Acinetobacter baylyi and Escherichia coli to reciprocally exchange essential amino acids. Interspecific coculture experiments confirmed that non-cooperating types were selectively favoured in spatially unstructured (liquid culture), yet disfavoured in spatially structured environments (agar plates). Both an individual-based model and experiments with engineered genotypes indicated that a segregation of cross-feeders and non-cooperating auxotrophs stabilized cooperative cross-feeding in spatially structured environments. Chemical imaging confirmed that auxotrophs were spatially excluded from cooperative benefits. Together, these results demonstrate that cooperative cross-feeding between different bacterial species is favoured in structured environments such as bacterial biofilms, suggesting this type of interactions might be common in natural bacterial communities.

Microstencils to generate defined, multi-species patterns of bacteria

Microbial communities are complex heterogeneous systems that are influenced by physical and chemical interactions with their environment, host, and community members. Techniques that facilitate the quantitative evaluation of how microscale organization influences the morphogenesis of multispecies communities could provide valuable insights into the dynamic behavior and organization of natural communities, the design of synthetic environments for multispecies culture, and the engineering of artificial consortia. In this work, we demonstrate a method for patterning microbes into simple arrangements that allow the quantitative measurement of growth dynamics as a function of their proximity to one another. The method combines parylene-based liftoff techniques with microfluidic delivery to simultaneously pattern multiple bacterial species with high viability using low-cost, customizable methods. Quantitative measurements of bacterial growth for two competing isolates demonstrate that spatial coordination can play a critical role in multispecies growth and structure.

Single gene locus changes perturb complex microbial communities as much as apex predator loss

Many bacterial species are highly social, adaptively shaping their local environment through the production of secreted molecules. This can, in turn, alter interaction strengths among species and modify community composition. However, the relative importance of such behaviours in determining the structure of complex communities is unknown. Here we show that single-locus changes affecting biofilm formation phenotypes in Bacillus subtilis modify community structure to the same extent as loss of an apex predator and even to a greater extent than loss of B. subtilis itself. These results, from experimentally manipulated multitrophic microcosm assemblages, demonstrate that bacterial social traits are key modulators of the structure of their communities. Moreover, they show that intraspecific genetic variability can be as important as strong trophic interactions in determining community dynamics. Microevolution may therefore be as important as species extinctions in shaping the response of microbial communities to environmental change.

Some species of social bacteria can chemically modify their nutrient environments, which may influence community interactions. Here, McClean et al. show that changes at a single gene locus in a biofilm-forming bacteria can perturb community structure to the same extent as the loss of an apex predator.

Development of Spatial Distribution Patterns by Biofilm Cells

Confined spatial patterns of microbial distribution are prevalent in nature, such as in microbial mats, soil communities, and water stream biofilms. The symbiotic two-species consortium of Pseudomonas putida and Acinetobacter sp. strain C6, originally isolated from a creosote-polluted aquifer, has evolved a distinct spatial organization in the laboratory that is characterized by an increased fitness and productivity. In this consortium, P. putida is reliant on microcolonies formed by Acinetobacter sp. C6, to which it attaches. Here we describe the processes that lead to the microcolony pattern by Acinetobacter sp. C6. Ecological spatial pattern analyses revealed that the microcolonies were not entirely randomly distributed and instead were arranged in a uniform pattern. Detailed time-lapse confocal microscopy at the single-cell level demonstrated that the spatial pattern was the result of an intriguing self-organization: small multicellular clusters moved along the surface to fuse with one another to form microcolonies. This active distribution capability was dependent on environmental factors (carbon source and oxygen) and historical contingency (formation of phenotypic variants). The findings of this study are discussed in the context of species distribution patterns observed in macroecology, and we summarize observations about the processes involved in coadaptation between P. putida and Acinetobacter sp. C6. Our results contribute to an understanding of spatial species distribution patterns as they are observed in nature, as well as the ecology of engineered communities that have the potential for enhanced and sustainable bioprocessing capacity.

Fast-moving bacteria self-organize into active two-dimensional crystals of rotating cells

We investigate a new form of collective dynamics displayed by Thiovulum majus, one of the fastest-swimming bacteria known. Cells spontaneously organize on a surface into a visually striking two-dimensional hexagonal lattice of rotating cells. As each constituent cell rotates its flagella, it creates a tornadolike flow that pulls neighboring cells towards and around it. As cells rotate against their neighbors, they exert forces on one another, causing the crystal to rotate and cells to reorganize. We show how these dynamics arise from hydrodynamic and steric interactions between cells. We derive the equations of motion for a crystal, show that this model explains several aspects of the observed dynamics, and discuss the stability of these active crystals.

Slow and steady wins the race: a bacterial exploitative competition strategy in fluctuating environments

One promising frontier for synthetic biology is the development of synthetic ecologies, whereby interacting species form an additional layer of connectivity for engineered gene circuits. Toward this goal, an important step is to understand different types of bacterial interactions in natural settings, among which competition is the most prevalent. By constructing a two-species population dynamics model, here, we mimicked bacterial growth in nature with resource-limited fluctuating environments and searched for optimal strategies for bacterial exploitative competition. In a simple game with two strategy options (constant or susceptible growth), we found that the species playing the constant growth strategy always outplays or is evenly matched with its competitor, suggesting that constant growth is a "no-loss" good bet. We also showed that adoption of sophisticated strategies enables a species to maximize its fitness when its competitor grows susceptibly. The pursuit of fitness maximization is, however, associated with potential loss if both species are capable of strategy adjustment, indicating an intrinsic risk-return trade-off. These findings offer new insights into bacterial competition and may also facilitate the engineering of microbial consortia for synthetic biology applications.

Staphylococcus aureus dry stress survivors have a heritable fitness advantage in subsequent dry exposure

Staphylococcus aureus is a major cause of hospital-acquired infections. The ability to survive on abiotic surfaces is an important characteristic that facilitates transmission between human hosts. We found that S. aureus survivors of dry surface incubation are resistant to subsequent dry stress exposure. Survivors also had reduced sensitivity to the disinfectant chlorhexidine gluconate, but not to ethanol. By using a set of mutants in cardiolipin synthase genes, we further demonstrated that the housekeeping cardiolipin synthase, Cls2, was significant for survival on dry surface. Taken together, this study provides insights into S. aureus survival outside of a host.

Computational Studies of the Intestinal Host-Microbiota Interactome

A large and growing body of research implicates aberrant immune response and compositional shifts of the intestinal microbiota in the pathogenesis of many intestinal disorders. The molecular and physical interaction between the host and the microbiota, known as the host-microbiota interactome, is one of the key drivers in the pathophysiology of many of these disorders. This host-microbiota interactome is a set of dynamic and complex processes, and needs to be treated as a distinct entity and subject for study. Disentangling this complex web of interactions will require novel approaches, using a combination of data-driven bioinformatics with knowledge-driven computational modeling. This review describes the computational approaches for investigating the host-microbiota interactome, with emphasis on the human intestinal tract and innate immunity, and highlights open challenges and existing gaps in the computation methodology for advancing our knowledge about this important facet of human health.

Microfluidic static droplet array for analyzing microbial communication on a population gradient

Quorum sensing (QS) is a type of cell-cell communication using signal molecules that are released and detected by cells, which respond to changes in their population density. A few studies explain that QS may operate in a density-dependent manner; however, due to experimental challenges, this fundamental hypothesis has never been investigated. Here, we present a microfluidic static droplet array (SDA) that combines a droplet generator with hydrodynamic traps to independently generate a bacterial population gradient into a parallel series of droplets under complete chemical and physical isolation. The SDA independently manipulates both a chemical concentration gradient and a bacterial population density. In addition, the bacterial population gradient in the SDA can be tuned by a simple change in the number of sample plug loading. Finally, the method allows the direct analysis of complicated biological events in an addressable droplet to enable the characterization of bacterial communication in response to the ratio of two microbial populations, including two genetically engineered QS circuits, such as the signal sender for acyl-homoserine lactone (AHL) production and the signal receiver bacteria for green fluorescent protein (GFP) expression induced by AHL. For the first time, we found that the population ratio of the signal sender and receiver indicates a significant and potentially interesting partnership between microbial communities. Therefore, we envision that this simple SDA could be a useful platform in various research fields, including analytical chemistry, combinatorial chemistry, synthetic biology, microbiology, and molecular biology.

Fap2 of Fusobacterium nucleatum is a galactose-inhibitable adhesin involved in coaggregation, cell adhesion, and preterm birth

Fusobacterium nucleatum is a common oral anaerobe involved in periodontitis that is known to translocate and cause intrauterine infections. In the oral environment, F. nucleatum adheres to a large diversity of species, facilitating their colonization and creating biological bridges that stabilize the multispecies dental biofilm. Many of these interactions (called coadherences or coaggregations) are galactose sensitive. Galactose-sensitive interactions are also involved in the binding of F. nucleatum to host cells. Hemagglutination of some F. nucleatum strains is also galactose sensitive, suggesting that a single galactose-sensitive adhesin might mediate the interaction of fusobacteria with many partners and targets. In order to identify the fusobacterial galactose-sensitive adhesin, a system for transposon mutagenesis in fusobacteria was created. The mutant library was screened for hemagglutination deficiency, and three clones were isolated. All three clones were found to harbor the transposon in the gene coding for the Fap2 outer membrane autotransporter. The three fap2 mutants failed to show galactose-inhibitable coaggregation with Porphyromonas gingivalis and were defective in cell binding. A fap2 mutant also showed a 2-log reduction in murine placental colonization compared to that of the wild type. Our results suggest that Fap2 is a galactose-sensitive hemagglutinin and adhesin that is likely to play a role in the virulence of fusobacteria.

Biofilm growth mode promotes maximum carrying capacity and community stability during product inhibition syntrophy

Sulfate-reducing bacteria (SRB) can interact syntrophically with other community members in the absence of sulfate, and interactions with hydrogen-consuming methanogens are beneficial when these archaea consume potentially inhibitory H₂ produced by the SRB. A dual continuous culture approach was used to characterize population structure within a syntrophic biofilm formed by the SRB Desulfovibrio vulgaris Hildenborough and the methanogenic archaeum Methanococcus maripaludis. Under the tested conditions, monocultures of D. vulgaris formed thin, stable biofilms, but monoculture M. maripaludis did not. Microscopy of intact syntrophic biofilm confirmed that D. vulgaris formed a scaffold for the biofilm, while intermediate and steady-state images revealed that M. maripaludis joined the biofilm later, likely in response to H₂ produced by the SRB. Close interactions in structured biofilm allowed efficient transfer of H₂ to M. maripaludis, and H₂ was only detected in cocultures with a mutant SRB that was deficient in biofilm formation (ΔpilA). M. maripaludis produced more carbohydrate (uronic acid, hexose, and pentose) as a monoculture compared to total coculture biofilm, and this suggested an altered carbon flux during syntrophy. The syntrophic biofilm was structured into ridges (∼300 × 50 μm) and models predicted lactate limitation at ∼50 μm biofilm depth. The biofilm had structure that likely facilitated mass transfer of H₂ and lactate, yet maximized biomass with a more even population composition (number of each organism) when compared to the bulk-phase community. Total biomass protein was equivalent in lactate-limited and lactate-excess conditions when a biofilm was present, but in the absence of biofilm, total biomass protein was significantly reduced. The results suggest that multispecies biofilms create an environment conducive to resource sharing, resulting in increased biomass retention, or carrying capacity, for cooperative populations.

The demographic determinants of human microbiome health

The human microbiome is a vast reservoir of microbial diversity and increasingly recognized to have a fundamental role in human health. In polymicrobial communities, the presence of one species can modulate the demography (i.e., growth and distribution) of other species. These demographic impacts generate feedbacks in multispecies interactions, which can be magnified in spatially structured populations (e.g., host-associated communities). Here, we argue that demographic feedbacks between species are central to microbiome development, shaping whether and how potential metabolic interactions come to be realized between expanding lineages of bacteria. Understanding how demographic feedbacks tune metabolic interactions and in turn shape microbiome structure and function is now a key challenge to our abilities to better manage microbiome health.

Bacteria dialog with Santa Rosalia: Are aggregations of cosmopolitan bacteria mainly explained by habitat filtering or by ecological interactions?

BACKGROUND

Since the landmark Santa Rosalia paper by Hutchinson, niche theory addresses the determinants of biodiversity in terms of both environmental and biological aspects. Disentangling the role of habitat filtering and interactions with other species is critical for understanding microbial ecology. Macroscopic biogeography explores hypothetical ecological interactions through the analysis of species associations. These methods have started to be incorporated into microbial ecology relatively recently, due to the inherent experimental difficulties and the coarse grained nature of the data.

RESULTS

Here we investigate the influence of environmental preferences and ecological interactions in the tendency of bacterial taxa to either aggregate or segregate, using a comprehensive dataset of bacterial taxa observed in a wide variety of environments. We assess significance of taxa associations through a null model that takes into account habitat preferences and the global distribution of taxa across samples. The analysis of these associations reveals a surprisingly large number of significant aggregations between taxa, with a marked community structure and a strong propensity to aggregate for cosmopolitan taxa. Due to the coarse grained nature of our data we cannot conclusively reject the hypothesis that many of these aggregations are due to environmental preferences that the null model fails to reproduce. Nevertheless, some observations are better explained by ecological interactions than by habitat filtering. In particular, most pairs of aggregating taxa co-occur in very different environments, which makes it unlikely that these associations are due to habitat preferences, and many are formed by cosmopolitan taxa without well defined habitat preferences. Moreover, known cooperative interactions are retrieved as aggregating pairs of taxa. As observed in similar studies, we also found that phylogenetically related taxa are much more prone to aggregate than to segregate, an observation that may play a role in bacterial speciation.

CONCLUSIONS

We hope that these results stimulate experimental verification of the putative cooperative interactions between cosmopolitan bacteria, and we suggest several groups of aggregated cosmopolitan bacteria that are interesting candidates for such an investigation.

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

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.