skip to main content


Title: Model-based and phylogenetically adjusted quantification of metabolic interaction between microbial species
Microbial community members exhibit various forms of interactions. Taking advantage of the increasing availability of microbiome data, many computational approaches have been developed to infer bacterial interactions from the co-occurrence of microbes across diverse microbial communities. Additionally, the introduction of genome-scale metabolic models have also enabled the inference of cooperative and competitive metabolic interactions between bacterial species. By nature, phylogenetically similar microbial species are more likely to share common functional profiles or biological pathways due to their genomic similarity. Without properly factoring out the phylogenetic relationship, any estimation of the competition and cooperation between species based on functional/pathway profiles may bias downstream applications. To address these challenges, we developed a novel approach for estimating the competition and complementarity indices for a pair of microbial species, adjusted by their phylogenetic distance. An automated pipeline, PhyloMint, was implemented to construct competition and complementarity indices from genome scale metabolic models derived from microbial genomes. Application of our pipeline to 2,815 human-gut associated bacteria showed high correlation between phylogenetic distance and metabolic competition/cooperation indices among bacteria. Using a discretization approach, we were able to detect pairs of bacterial species with cooperation scores significantly higher than the average pairs of bacterial species with similar phylogenetic distances. A network community analysis of high metabolic cooperation but low competition reveals distinct modules of bacterial interactions. Our results suggest that niche differentiation plays a dominant role in microbial interactions, while habitat filtering also plays a role among certain clades of bacterial species.  more » « less
Award ID(s):
2025451
NSF-PAR ID:
10297357
Author(s) / Creator(s):
; ; ;
Editor(s):
Ouzounis, Christos A.
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
16
Issue:
10
ISSN:
1553-7358
Page Range / eLocation ID:
e1007951
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Campbell, Barbara J. (Ed.)
    ABSTRACT In nutrient-limited conditions, plants rely on rhizosphere microbial members to facilitate nutrient acquisition, and in return, plants provide carbon resources to these root-associated microorganisms. However, atmospheric nutrient deposition can affect plant-microbe relationships by changing soil bacterial composition and by reducing cooperation between microbial taxa and plants. To examine how long-term nutrient addition shapes rhizosphere community composition, we compared traits associated with bacterial (fast-growing copiotrophs, slow-growing oligotrophs) and plant (C 3 forb, C 4 grass) communities residing in a nutrient-poor wetland ecosystem. Results revealed that oligotrophic taxa dominated soil bacterial communities and that fertilization increased the presence of oligotrophs in bulk and rhizosphere communities. Additionally, bacterial species diversity was greatest in fertilized soils, particularly in bulk soils. Nutrient enrichment (fertilized versus unfertilized) and plant association (bulk versus rhizosphere) determined bacterial community composition; bacterial community structure associated with plant functional group (grass versus forb) was similar within treatments but differed between fertilization treatments. The core forb microbiome consisted of 602 unique taxa, and the core grass microbiome consisted of 372 unique taxa. Forb rhizospheres were enriched in potentially disease-suppressive bacterial taxa, and grass rhizospheres were enriched in bacterial taxa associated with complex carbon decomposition. Results from this study demonstrate that fertilization serves as a strong environmental filter on the soil microbiome, which leads to distinct rhizosphere communities and can shift plant effects on the rhizosphere microbiome. These taxonomic shifts within plant rhizospheres could have implications for plant health and ecosystem functions associated with carbon and nitrogen cycling. IMPORTANCE Over the last century, humans have substantially altered nitrogen and phosphorus cycling. Use of synthetic fertilizer and burning of fossil fuels and biomass have increased nitrogen and phosphorus deposition, which results in unintended fertilization of historically low-nutrient ecosystems. With increased nutrient availability, plant biodiversity is expected to decline, and the abundance of copiotrophic taxa is anticipated to increase in bacterial communities. Here, we address how bacterial communities associated with different plant functional types (forb, grass) shift due to long-term nutrient enrichment. Unlike other studies, results revealed an increase in bacterial diversity, particularly of oligotrophic bacteria in fertilized plots. We observed that nutrient addition strongly determines forb and grass rhizosphere composition, which could indicate different metabolic preferences in the bacterial communities. This study highlights how long-term fertilization of oligotroph-dominated wetlands could alter diversity and metabolism of rhizosphere bacterial communities in unexpected ways. 
    more » « less
  2. Abstract

    Most insects harbour influential, yet non‐essential heritable microbes in their hemocoel. Communities of these symbionts exhibit low diversity. But their frequent multi‐species nature raises intriguing questions on roles for symbiont–symbiont synergies in host adaptation, and on the stability of the symbiont communities, themselves. In this study, we build on knowledge of species‐defined symbiont community structure across US populations of the pea aphid,Acyrthosiphon pisum. Through extensive symbiont genotyping, we show that pea aphids' microbiomes can be more precisely defined at the symbiont strain level, with strain variability shaping five out of nine previously reported co‐infection trends. Field data provide a mixture of evidence for synergistic fitness effects and symbiont hitchhiking, revealing causes and consequences of these co‐infection trends. To test whether within‐host metabolic interactions predict common versus rare strain‐defined communities, we leveraged the high relatedness of our dominant, community‐defined symbiont strains vs. 12 pea aphid‐derived Gammaproteobacteria with sequenced genomes. Genomic inference, using metabolic complementarity indices, revealed high potential for cooperation among one pair of symbionts—Serratia symbioticaandRickettsiella viridis. Applying the expansion network algorithm, through additional use of pea aphid and obligateBuchnerasymbiont genomes,SerratiaandRickettsiellaemerged as the only symbiont community requiring both parties to expand holobiont metabolism. Through their joint expansion of the biotin biosynthesis pathway, these symbionts may span missing gaps, creating a multi‐party mutualism within their nutrient‐limited, phloem‐feeding hosts. Recent, complementary gene inactivation, within the biotin pathways ofSerratiaandRickettsiella, raises further questions on the origins of mutualisms and host–symbiont interdependencies.

     
    more » « less
  3. ABSTRACT Microbes face a trade-off between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A “gedanken” (thought) experiment to investigate this trade-off would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment using in silico genome-scale models. In particular, we used mixed-integer linear programming to identify trade-off solutions in communities of Escherichia coli strains. The strategies that we found revealed a large space of opportunities in nuanced and nonintuitive metabolic division of labor, including, for example, splitting the tricarboxylic acid (TCA) cycle into two separate halves. The systematic computation of possible solutions in division of labor for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displayed a nonlinear boundary, indicating that the loss of an intracellular reaction was not necessarily compensated for by a single imported metabolite. Different regions in this landscape were associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable but are outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor. IMPORTANCE Understanding how microbes assemble into communities is a fundamental open issue in biology, relevant to human health, metabolic engineering, and environmental sustainability. A possible mechanism for interactions of microbes is through cross-feeding, i.e., the exchange of small molecules. These metabolic exchanges may allow different microbes to specialize in distinct tasks and evolve division of labor. To systematically explore the space of possible strategies for division of labor, we applied advanced optimization algorithms to computational models of cellular metabolism. Specifically, we searched for communities able to survive under constraints (such as a limited number of reactions) that would not be sustainable by individual species. We found that predicted consortia partition metabolic pathways in ways that would be difficult to identify manually, possibly providing a competitive advantage over individual organisms. In addition to helping understand diversity in natural microbial communities, our approach could assist in the design of synthetic consortia. 
    more » « less
  4. Davies, Julian E. (Ed.)
    ABSTRACT Bacteria isolated from soils are major sources of specialized metabolites, including antibiotics and other compounds with clinical value that likely shape interactions among microbial community members and impact biogeochemical cycles. Yet, isolated lineages represent a small fraction of all soil bacterial diversity. It remains unclear how the production of specialized metabolites varies across the phylogenetic diversity of bacterial species in soils and whether the genetic potential for production of these metabolites differs with soil depth and vegetation type within a geographic region. We sampled soils and saprolite from three sites in a northern California Critical Zone Observatory with various vegetation and bedrock characteristics and reconstructed 1,334 metagenome-assembled genomes containing diverse biosynthetic gene clusters (BGCs) for secondary metabolite production. We obtained genomes for prolific producers of secondary metabolites, including novel groups within the Actinobacteria , Chloroflexi , and candidate phylum “ Candidatus Dormibacteraeota.” Surprisingly, one genome of a candidate phyla radiation (CPR) bacterium coded for a ribosomally synthesized linear azole/azoline-containing peptide, a capacity we found in other publicly available CPR bacterial genomes. Overall, bacteria with higher biosynthetic potential were enriched in shallow soils and grassland soils, with patterns of abundance of BGC type varying by taxonomy. IMPORTANCE Microbes produce specialized compounds to compete or communicate with one another and their environment. Some of these compounds, such as antibiotics, are also useful in medicine and biotechnology. Historically, most antibiotics have come from soil bacteria which can be isolated and grown in the lab. Though the vast majority of soil bacteria cannot be isolated, we can extract their genetic information and search it for genes which produce these specialized compounds. These understudied soil bacteria offer a wealth of potential for the discovery of new and important microbial products. Here, we identified the ability to produce these specialized compounds in diverse and novel bacteria in a range of soil environments. This information will be useful to other researchers who wish to isolate certain products. Beyond their use to humans, understanding the distribution and function of microbial products is key to understanding microbial communities and their effects on biogeochemical cycles. 
    more » « less
  5. Abstract Background Microbial colonization of subsurface shales following hydraulic fracturing offers the opportunity to study coupled biotic and abiotic factors that impact microbial persistence in engineered deep subsurface ecosystems. Shale formations underly much of the continental USA and display geographically distinct gradients in temperature and salinity. Complementing studies performed in eastern USA shales that contain brine-like fluids, here we coupled metagenomic and metabolomic approaches to develop the first genome-level insights into ecosystem colonization and microbial community interactions in a lower-salinity, but high-temperature western USA shale formation. Results We collected materials used during the hydraulic fracturing process (i.e., chemicals, drill muds) paired with temporal sampling of water produced from three different hydraulically fractured wells in the STACK ( S ooner T rend A nadarko Basin, C anadian and K ingfisher) shale play in OK, USA. Relative to other shale formations, our metagenomic and metabolomic analyses revealed an expanded taxonomic and metabolic diversity of microorganisms that colonize and persist in fractured shales. Importantly, temporal sampling across all three hydraulic fracturing wells traced the degradation of complex polymers from the hydraulic fracturing process to the production and consumption of organic acids that support sulfate- and thiosulfate-reducing bacteria. Furthermore, we identified 5587 viral genomes and linked many of these to the dominant, colonizing microorganisms, demonstrating the key role that viral predation plays in community dynamics within this closed, engineered system. Lastly, top-side audit sampling of different source materials enabled genome-resolved source tracking, revealing the likely sources of many key colonizing and persisting taxa in these ecosystems. Conclusions These findings highlight the importance of resource utilization and resistance to viral predation as key traits that enable specific microbial taxa to persist across fractured shale ecosystems. We also demonstrate the importance of materials used in the hydraulic fracturing process as both a source of persisting shale microorganisms and organic substrates that likely aid in sustaining the microbial community. Moreover, we showed that different physicochemical conditions (i.e., salinity, temperature) can influence the composition and functional potential of persisting microbial communities in shale ecosystems. Together, these results expand our knowledge of microbial life in deep subsurface shales and have important ramifications for management and treatment of microbial biomass in hydraulically fractured wells. 
    more » « less