Global epistasis and the emergence of function in microbial consortia.

The many functions of microbial communities emerge from a complex web of interactions between organisms and their environment. This poses a significant obstacle to engineering microbial consortia, hindering our ability to harness the potential of microorganisms for biotechnological applications. In this study, we demonstrate that the collective effect of ecological interactions between microbes in a community can be captured by simple statistical models that predict how adding a new species to a community will affect its function. These predictive models mirror the patterns of global epistasis reported in genetics, and they can be quantitatively interpreted in terms of pairwise interactions between community members. Our results illuminate an unexplored path to quantitatively predicting the function of microbial consortia from their composition, paving the way to optimizing desirable community properties and bringing the tasks of predicting biological function at the genetic, organismal, and ecological scales under the same quantitative formalism.

Genomic structure predicts metabolite dynamics in microbial communities.

The metabolic activities of microbial communities play a defining role in the evolution and persistence of life on Earth, driving redox reactions that give rise to global biogeochemical cycles. Community metabolism emerges from a hierarchy of processes, including gene expression, ecological interactions, and environmental factors. In wild communities, gene content is correlated with environmental context, but predicting metabolite dynamics from genomes remains elusive. Here, we show, for the process of denitrification, that metabolite dynamics of a community are predictable from the genes each member of the community possesses. A simple linear regression reveals a sparse and generalizable mapping from gene content to metabolite dynamics for genomically diverse bacteria. A consumer-resource model correctly predicts community metabolite dynamics from single-strain phenotypes. Our results demonstrate that the conserved impacts of metabolic genes can predict community metabolite dynamics, enabling the prediction of metabolite dynamics from metagenomes, designing denitrifying communities, and discovering how genome evolution impacts metabolism.

Atlas of clinically distinct cell states and ecosystems across human solid tumors.

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

Bifunctional Immunity Proteins Protect Bacteria against FtsZ-Targeting ADP-Ribosylating Toxins.

ADP-ribosylation of proteins can profoundly impact their function and serves as an effective mechanism by which bacterial toxins impair eukaryotic cell processes. Here, we report the discovery that bacteria also employ ADP-ribosylating toxins against each other during interspecies competition. We demonstrate that one such toxin from Serratia proteamaculans interrupts the division of competing cells by modifying the essential bacterial tubulin-like protein, FtsZ, adjacent to its protomer interface, blocking its capacity to polymerize. The structure of the toxin in complex with its immunity determinant revealed two distinct modes of inhibition: active site occlusion and enzymatic removal of ADP-ribose modifications. We show that each is sufficient to support toxin immunity; however, the latter additionally provides unprecedented broad protection against non-cognate ADP-ribosylating effectors. Our findings reveal how an interbacterial arms race has produced a unique solution for safeguarding the integrity of bacterial cell division machinery against inactivating post-translational modifications.

Microbial Interkingdom Interactions in Roots Promote Arabidopsis Survival.

Roots of healthy plants are inhabited by soil-derived bacteria, fungi, and oomycetes that have evolved independently in distinct kingdoms of life. How these microorganisms interact and to what extent those interactions affect plant health are poorly understood. We examined root-associated microbial communities from three Arabidopsis thaliana populations and detected mostly negative correlations between bacteria and filamentous microbial eukaryotes. We established microbial culture collections for reconstitution experiments using germ-free A. thaliana. In plants inoculated with mono- or multi-kingdom synthetic microbial consortia, we observed a profound impact of the bacterial root microbiota on fungal and oomycetal community structure and diversity. We demonstrate that the bacterial microbiota is essential for plant survival and protection against root-derived filamentous eukaryotes. Deconvolution of 2,862 binary bacterial-fungal interactions ex situ, combined with community perturbation experiments in planta, indicate that biocontrol activity of bacterial root commensals is a redundant trait that maintains microbial interkingdom balance for plant health.

Species-Independent Attraction to Biofilms through Electrical Signaling.

Bacteria residing within biofilm communities can coordinate their behavior through cell-to-cell signaling. However, it remains unclear if these signals can also influence the behavior of distant cells that are not part of the community. Using a microfluidic approach, we find that potassium ion channel-mediated electrical signaling generated by a Bacillus subtilis biofilm can attract distant cells. Integration of experiments and mathematical modeling indicates that extracellular potassium emitted from the biofilm alters the membrane potential of distant cells, thereby directing their motility. This electrically mediated attraction appears to be a generic mechanism that enables cross-species interactions, as Pseudomonas aeruginosa cells also become attracted to the electrical signal released by the B. subtilis biofilm. Cells within a biofilm community can thus not only coordinate their own behavior but also influence the behavior of diverse bacteria at a distance through long-range electrical signaling. PAPERCLIP.

Culturing Human Gut Microbiomes in the Laboratory.

The human gut microbiota is a complex community of prokaryotic and eukaryotic microbes and viral particles that is increasingly associated with many aspects of host physiology and health. However, the classical microbiology approach of axenic culture cannot provide a complete picture of the complex interactions between microbes and their hosts in vivo. As such, recently there has been much interest in the culture of gut microbial ecosystems in the laboratory as a strategy to better understand their compositions and functions. In this review, we discuss the model platforms and methods available in the contemporary microbiology laboratory to study human gut microbiomes, as well as current knowledge surrounding the isolation of human gut microbes for the potential construction of defined communities for use in model systems.

Community interactions drive the evolution of antibiotic tolerance in bacteria.

The emergence of antibiotic tolerance (prolonged survival against exposure) in natural bacterial populations is a major concern. Since it has been studied primarily in isogenic populations, we do not yet understand how ecological interactions in a diverse community impact the evolution of tolerance. To address this, we studied the evolutionary dynamics of a synthetic bacterial community composed of two interacting strains. In this community, an antibiotic-resistant strain protected the other, susceptible strain by degrading the antibiotic ampicillin in the medium. Surprisingly, we found that in the presence of antibiotics, the susceptible strain evolved tolerance. Tolerance was typified by an increase in survival as well as an accompanying decrease in the growth rate, highlighting a trade-off between the two. A simple mathematical model explained that the observed decrease in the death rate, even when coupled with a decreased growth rate, is beneficial in a community with weak protective interactions. In the presence of strong interactions, the model predicted that the trade-off would instead be detrimental, and tolerance would not emerge, which we experimentally verified. By whole genome sequencing the evolved tolerant isolates, we identified two genetic hot spots which accumulated mutations in parallel lines, suggesting their association with tolerance. Our work highlights that ecological interactions can promote antibiotic tolerance in bacterial communities, which has remained understudied.

Environmental fluctuations explain the universal decay of species-abundance correlations with phylogenetic distance.

Multiple ecological forces act together to shape the composition of microbial communities. Phyloecology approaches-which combine phylogenetic relationships between species with community ecology-have the potential to disentangle such forces but are often hard to connect with quantitative predictions from theoretical models. On the other hand, macroecology, which focuses on statistical patterns of abundance and diversity, provides natural connections with theoretical models but often neglects interspecific correlations and interactions. Here, we propose a unified framework combining both such approaches to analyze microbial communities. In particular, by using both cross-sectional and longitudinal metagenomic data for species abundances, we reveal the existence of an empirical macroecological law establishing that correlations in species-abundance fluctuations across communities decay from positive to null values as a function of phylogenetic dissimilarity in a consistent manner across ecologically distinct microbiomes. We formulate three variants of a mechanistic model-each relying on alternative ecological forces-that lead to radically different predictions. From these analyses, we conclude that the empirically observed macroecological pattern can be quantitatively explained as a result of shared population-independent fluctuating resources, i.e., environmental filtering and not as a consequence of, e.g., species competition. Finally, we show that the macroecological law is also valid for temporal data of a single community and that the properties of delayed temporal correlations can be reproduced as well by the model with environmental filtering.

Intransitivity as a dynamic assembly engine of competitive communities.

Historically, those ecological communities thought to be dominated by competitive interactions among their component species have been assumed to exhibit transitive competition, that is, a hierarchy of competitive strength from most dominant to most submissive. A surge of recent literature takes issue with this assumption and notes that some species in some communities are intransitive, where a rock/scissors/paper arrangement characterizes some components of some communities. We here propose a merging of these two ideas, wherein an intransitive subgroup of species connects with a distinct subcomponent that is organized hierarchically, such that the expected eventual takeover by the dominant competitor in the hierarchy is thwarted, and the entire community can be sustained. This means that the combination of transitive and intransitive structures can maintain many species even when competition is strong. Here, we develop this theoretical framework using a simple variant on the Lotka-Volterra competition equations to illustrate the process. We also present data for the ant community in a coffee agroecosystem in Puerto Rico, that appears to be organized in this way. A detailed study on one typical coffee farm illustrates an intransitive loop of three species that seems to maintain a distinct competitive community of at least 13 additional species.

Minorities drive growth resumption in cross-feeding microbial communities.

Microbial communities are fundamental to life on Earth. Different strains within these communities are often connected by a highly connected metabolic network, where the growth of one strain depends on the metabolic activities of other community members. While distributed metabolic functions allow microbes to reduce costs and optimize metabolic pathways, they make them metabolically dependent. Here, we hypothesize that such dependencies can be detrimental in situations where the external conditions change rapidly, as they often do in natural environments. After a shift in external conditions, microbes need to remodel their metabolism, but they can only resume growth once partners on which they depend have also adapted to the new conditions. It is currently not well understood how microbial communities resolve this dilemma and how metabolic interactions are reestablished after an environmental shift. To address this question, we investigated the dynamical responses to environmental perturbation by microbial consortia with distributed anabolic functions. By measuring the regrowth times at the single-cell level in spatially structured communities, we found that metabolic dependencies lead to a growth delay after an environmental shift. However, a minority of cells-those in the immediate neighborhood of their metabolic partners-can regrow quickly and come to numerically dominate the community after the shift. The spatial arrangement of a microbial community is thus a key factor in determining the communities' ability to maintain metabolic interactions and growth in fluctuating conditions. Our results suggest that environmental fluctuations can limit the emergence of metabolic dependencies between microorganisms.

Tracing human life trajectory using gut microbial communities by context-aware deep learning.

The gut microbial communities are highly plastic throughout life, and the human gut microbial communities show spatial-temporal dynamic patterns at different life stages. However, the underlying association between gut microbial communities and time-related factors remains unclear. The lack of context-awareness, insufficient data, and the existence of batch effect are the three major issues, making the life trajection of the host based on gut microbial communities problematic. Here, we used a novel computational approach (microDELTA, microbial-based deep life trajectory) to track longitudinal human gut microbial communities' alterations, which employs transfer learning for context-aware mining of gut microbial community dynamics at different life stages. Using an infant cohort, we demonstrated that microDELTA outperformed Neural Network for accurately predicting the age of infant with different delivery mode, especially for newborn infants of vaginal delivery with the area under the receiver operating characteristic curve of microDELTA and Neural Network at 0.811 and 0.436, respectively. In this context, we have discovered the influence of delivery mode on infant gut microbial communities. Along the human lifespan, we also applied microDELTA to a Chinese traveler cohort, a Hadza hunter-gatherer cohort and an elderly cohort. Results revealed the association between long-term dietary shifts during travel and adult gut microbial communities, the seasonal cycling of gut microbial communities for the Hadza hunter-gatherers, and the distinctive microbial pattern of elderly gut microbial communities. In summary, microDELTA can largely solve the issues in tracing the life trajectory of the human microbial communities and generate accurate and flexible models for a broad spectrum of microbial-based longitudinal researches.

Evolution and stability of complex microbial communities driven by trade-offs.

Microbial communities are ubiquitous in nature and play an important role in ecology and human health. Cross-feeding is thought to be core to microbial communities, though it remains unclear precisely why it emerges. Why have multi-species microbial communities evolved in many contexts and what protects microbial consortia from invasion? Here, we review recent insights into the emergence and stability of coexistence in microbial communities. A particular focus is the long-term evolutionary stability of coexistence, as observed for microbial communities that spontaneously evolved in the E. coli long-term evolution experiment (LTEE). We analyze these findings in the context of recent work on trade-offs between competing microbial objectives, which can constitute a mechanistic basis for the emergence of coexistence. Coexisting communities, rather than monocultures of the 'fittest' single strain, can form stable endpoints of evolutionary trajectories. Hence, the emergence of coexistence might be an obligatory outcome in the evolution of microbial communities. This implies that rather than embodying fragile metastable configurations, some microbial communities can constitute formidable ecosystems that are difficult to disrupt.

Fast growth can counteract antibiotic susceptibility in shaping microbial community resilience to antibiotics.

Microbial communities often face external perturbations that can induce lasting changes in their composition and functions. Our understanding of how multispecies communities respond to perturbations such as antibiotics is limited, with susceptibility assays performed on individual, isolated species our primary guide in predicting community transitions. Here, we studied how bacterial growth dynamics can overcome differences in antibiotic susceptibility in determining community resilience: the recovery of the original community state following antibiotic exposure. We used an experimental community containing Corynebacterium ammoniagenes and Lactobacillus plantarum that displays two alternative stable states as a result of mutual inhibition. Although C. ammoniagenes was more susceptible to chloramphenicol in monocultures, we found that chloramphenicol exposure nonetheless led to a transition from the L. plantarum-dominated to the C. ammoniagenes-dominated community state. Combining theory and experiments, we demonstrated that growth rate differences between the two species made the L. plantarum-dominated community less resilient to several antibiotics with different mechanisms of action. Taking advantage of an observed cooperativity­a dependence on population abundance­in the growth of C. ammoniagenes, we next analyzed in silico scenarios that could compromise the high resilience of the C. ammoniagenes-dominated state. The model predicted that lowering the dispersal rate, through interacting with the growth at low population densities, could make the C. ammoniagenes state fragile against virtually any kind of antibiotic, a prediction that we confirmed experimentally. Our results highlight that species susceptibility to antibiotics is often uninformative of community resilience, as growth dynamics in the wake of antibiotic exposure can play a dominant role.

Factors influencing terrestriality in primates of the Americas and Madagascar.

Among mammals, the order Primates is exceptional in having a high taxonomic richness in which the taxa are arboreal, semiterrestrial, or terrestrial. Although habitual terrestriality is pervasive among the apes and African and Asian monkeys (catarrhines), it is largely absent among monkeys of the Americas (platyrrhines), as well as galagos, lemurs, and lorises (strepsirrhines), which are mostly arboreal. Numerous ecological drivers and species-specific factors are suggested to set the conditions for an evolutionary shift from arboreality to terrestriality, and current environmental conditions may provide analogous scenarios to those transitional periods. Therefore, we investigated predominantly arboreal, diurnal primate genera from the Americas and Madagascar that lack fully terrestrial taxa, to determine whether ecological drivers (habitat canopy cover, predation risk, maximum temperature, precipitation, primate species richness, human population density, and distance to roads) or species-specific traits (body mass, group size, and degree of frugivory) associate with increased terrestriality. We collated 150,961 observation hours across 2,227 months from 47 species at 20 sites in Madagascar and 48 sites in the Americas. Multiple factors were associated with ground use in these otherwise arboreal species, including increased temperature, a decrease in canopy cover, a dietary shift away from frugivory, and larger group size. These factors mostly explain intraspecific differences in terrestriality. As humanity modifies habitats and causes climate change, our results suggest that species already inhabiting hot, sparsely canopied sites, and exhibiting more generalized diets, are more likely to shift toward greater ground use.

Dynamics and drivers of fungal communities in a multipartite ant-plant association.

BACKGROUND: Fungi and ants belong to the most important organisms in terrestrial ecosystems on Earth. In nutrient-poor niches of tropical rainforests, they have developed steady ecological relationships as a successful survival strategy. In tropical ant-plant mutualisms worldwide, where resident ants provide the host plants with defense and nutrients in exchange for shelter and food, fungi are regularly found in the ant nesting space, inhabiting ant-made dark-colored piles ("patches"). Unlike the extensively investigated fungus-growing insects, where the fungi serve as the primary food source, the purpose of this ant-fungi association is less clear. To decipher the roles of fungi in these structures within ant nests, it is crucial to first understand the dynamics and drivers that influence fungal patch communities during ant colony development. RESULTS: In this study, we investigated how the ant colony age and the ant-plant species affect the fungal community in the patches. As model we selected one of the most common mutualisms in the Tropics of America, the Azteca-Cecropia complex. By amplicon sequencing of the internal transcribed spacer 2 (ITS2) region, we analyzed the patch fungal communities of 93 Azteca spp. colonies inhabiting Cecropia spp. trees. Our study demonstrates that the fungal diversity in patches increases as the ant colony grows and that a change in the prevalent fungal taxa occurs between initial and established patches. In addition, the ant species significantly influences the composition of the fungal community in established ant colonies, rather than the host plant species. CONCLUSIONS: The fungal patch communities become more complex as the ant colony develops, due to an acquisition of fungi from the environment and a substrate diversification. Our results suggest a successional progression of the fungal communities in the patches during ant colony growth and place the ant colony as the main driver shaping such communities. The findings of this study demonstrate the unexpectedly complex nature of ant-plant mutualisms in tropical regions at a micro scale.

A new approach to describe the taxonomic structure of microbiome and its application to assess the relationship between microbial niches.

BACKGROUND: Data from microbiomes from multiple niches is often collected, but methods to analyse these often ignore associations between niches. One interesting case is that of the oral microbiome. Its composition is receiving increasing attention due to reports on its associations with general health. While the oral cavity includes different niches, multi-niche microbiome data analysis is conducted using a single niche at a time and, therefore, ignores other niches that could act as confounding variables. Understanding the interaction between niches would assist interpretation of the results, and help improve our understanding of multi-niche microbiomes. METHODS: In this study, we used a machine learning technique called latent Dirichlet allocation (LDA) on two microbiome datasets consisting of several niches. LDA was used on both individual niches and all niches simultaneously. On individual niches, LDA was used to decompose each niche into bacterial sub-communities unveiling their taxonomic structure. These sub-communities were then used to assess the relationship between microbial niches using the global test. On all niches simultaneously, LDA allowed us to extract meaningful microbial patterns. Sets of co-occurring operational taxonomic units (OTUs) comprising those patterns were then used to predict the original location of each sample. RESULTS: Our approach showed that the per-niche sub-communities displayed a strong association between supragingival plaque and saliva, as well as between the anterior and posterior tongue. In addition, the LDA-derived microbial signatures were able to predict the original sample niche illustrating the meaningfulness of our sub-communities. For the multi-niche oral microbiome dataset we had an overall accuracy of 76%, and per-niche sensitivity of up to 83%. Finally, for a second multi-niche microbiome dataset from the entire body, microbial niches from the oral cavity displayed stronger associations to each other than with those from other parts of the body, such as niches within the vagina and the skin. CONCLUSION: Our LDA-based approach produces sets of co-occurring taxa that can describe niche composition. LDA-derived microbial signatures can also be instrumental in summarizing microbiome data, for both descriptions as well as prediction.

Flux sampling in genome-scale metabolic modeling of microbial communities.

BACKGROUND: Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions and identify strategies to optimize the community. Flux balance analysis (FBA) is most often used to predict the flux through all reactions in a genome-scale model; however, the fluxes predicted by FBA depend on a user-defined cellular objective. Flux sampling is an alternative to FBA, as it provides the range of fluxes possible within a microbial community. Furthermore, flux sampling can capture additional heterogeneity across a population, especially when cells exhibit sub-maximal growth rates. RESULTS: In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. With sampling, we find significant differences in the predicted metabolism, including an increase in cooperative interactions and pathway-specific changes in predicted flux. CONCLUSIONS: Our results suggest the importance of sampling-based approaches to evaluate metabolic interactions. Furthermore, we emphasize the utility of flux sampling in quantitatively studying interactions between cells and organisms.

When can higher-order interactions produce stable coexistence?

Most ecological models are based on the assumption that species interact in pairs. Diverse communities, however, can have higher-order interactions, in which two or more species jointly impact the growth of a third species. A pitfall of the common pairwise approach is that it misses the higher-order interactions potentially responsible for maintaining natural diversity. Here, we explore the stability properties of systems where higher-order interactions guarantee that a specified set of abundances is a feasible equilibrium of the dynamics. Even these higher-order interactions which lead to equilibria do not necessarily produce stable coexistence. Instead, these systems are more likely to be stable when the pairwise interactions are weak or facilitative. Correlations between the pairwise and higher-order interactions, however, do permit robust coexistence even in diverse systems. Our work not only reveals the challenges in generating stable coexistence through higher-order interactions but also uncovers interaction patterns that can enable diversity.

The three-species problem: Incorporating competitive asymmetry and intransitivity in modern coexistence theory.

While natural communities can contain hundreds of species, modern coexistence theory focuses primarily on species pairs. Alternatively, the structural stability approach considers the feasibility of equilibria, gaining scalability to larger communities but sacrificing information about dynamic stability. Three-species competitive communities are a bridge to more-diverse communities. They display novel phenomena while remaining amenable to mathematical analysis, but remain incompletely understood. Here, we combine these approaches to identify the key quantities that determine three-species competition outcomes. We show that pairwise niche overlap and fitness differences are insufficient to completely characterize competitive outcomes, which requires a strictly triplet-wise quantity: cyclic asymmetry, which underlies intransitivity. Low pairwise niche overlap stabilizes the triplet, while high fitness differences promote competitive exclusion. The effect of cyclic asymmetry on stability is complex and depends on pairwise niche overlap. In summary, we elucidate how pairwise niche overlap, fitness differences and cyclic asymmetry determine three-species competition outcomes.