metaGEM: reconstruction of genome scale metabolic models directly from metagenomes

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Revealing microbiota characteristics and predicting flavor-producing sub-communities in Nongxiangxing baijiu pit mud through metagenomic analysis and metabolic modeling

The microorganisms of the pit mud (PM) of Nongxiangxing baijiu (NXXB) have an important role in the synthesis of flavor substances, and they determine attributes and quality of baijiu. Herein, we utilize metagenomics and genome-scale metabolic models (GSMMs) to investigate the microbial composition, metabolic functions in PM microbiota, as well as to identify microorganisms and communities linked to flavor compounds. Metagenomic data revealed that the most prevalent assembly of bacteria and archaea was Proteiniphilum, Caproicibacterium, Petrimonas, Lactobacillus, Clostridium, Aminobacterium, Syntrophomonas, Methanobacterium, Methanoculleus, and Methanosarcina. The important enzymes ofPMwere in bothGH and GT familymetabolism. A total of 38 high-quality metagenome-assembled genomes (MAGs) were obtained, including those at the family level (n = 13), genus level (n = 17), and species level (n = 8). GSMMs of the 38 MAGs were then constructed. From the GSMMs, individual and community capabilities respectively were predicted to be able to produce 111 metabolites and 598 metabolites. Twenty-three predicted metabolites were consistent with the metabonomics detected flavors and served as targets. Twelve sub-community of were screened by cross-feeding of 38 GSMMs. Of them, Methanobacterium, Sphaerochaeta, Muricomes intestini, Methanobacteriaceae, Synergistaceae, and Caloramator were core microorganisms for targets in each sub-community. Overall, this study of metagenomic and target-community screening could help our understanding of the metabolite-microbiome association and further bioregulation of baijiu.

Engineering agricultural soil microbiomes and predicting plant phenotypes

Plant growth-promoting rhizobacteria (PGPR) can improve crop yields, nutrient use efficiency, plant tolerance to stressors, and confer benefits to future generations of crops grown in the same soil. Unlocking the potential of microbial communities in the rhizosphere and endosphere is therefore of great interest for sustainable agriculture advancements. Before plant microbiomes can be engineered to confer desirable phenotypic effects on their plant hosts, a deeper understanding of the interacting factors influencing rhizosphere community structure and function is needed. Dealing with this complexity is becoming more feasible using computational approaches. In this review, we discuss recent advances at the intersection of experimental and computational strategies for the investigation of plant-microbiome interactions and the engineering of desirable soil microbiomes.

A benchmark of optimization solvers for genome-scale metabolic modeling of organisms and communities

Genome-scale metabolic modeling is a powerful framework for predicting metabolic phenotypes of any organism with an annotated genome. For two decades, this framework has been used for the rational design of microbial cell factories. In the last decade, the range of applications has exploded, and new frontiers have emerged, including the study of the gut microbiome and its health implications and the role of microbial communities in global ecosystems. However, all the critical steps in this framework, from model construction to simulation, require the use of powerful linear optimization solvers, with the choice often relying on commercial solvers for their well-known computational efficiency. In this work, I benchmark a total of six solvers (two commercial and four open source) and measure their performance to solve linear and mixed-integer linear problems of increasing complexity. Although commercial solvers are still the fastest, at least two open-source solvers show comparable performance. These results show that genome-scale metabolic modeling does not need to be hindered by commercial licensing schemes and can become a truly open science framework for solving urgent societal challenges.IMPORTANCEModeling the metabolism of organisms and communities allows for computational exploration of their metabolic capabilities and testing their response to genetic and environmental perturbations. This holds the potential to address multiple societal issues related to human health and the environment. One of the current limitations is the use of commercial optimization solvers with restrictive licenses for academic and non-academic use. This work compares the performance of several commercial and open-source solvers to solve some of the most complex problems in the field. Benchmarking results show that, although commercial solvers are indeed faster, some of the open-source options can also efficiently tackle the hardest problems, showing great promise for the development of open science applications.

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.

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Microbial interactions shape cheese flavour formation

Cheese fermentation and flavour formation are the result of complex biochemical reactions driven by the activity of multiple microorganisms. Here, we studied the roles of microbial interactions in flavour formation in a year-long Cheddar cheese making process, using a commercial starter culture containing Streptococcus thermophilus and Lactococcus strains. By using an experimental strategy whereby certain strains were left out from the starter culture, we show that S. thermophilus has a crucial role in boosting Lactococcus growth and shaping flavour compound profile. Controlled milk fermentations with systematic exclusion of single Lactococcus strains, combined with genomics, genome-scale metabolic modelling, and metatranscriptomics, indicated that S. thermophilus proteolytic activity relieves nitrogen limitation for Lactococcus and boosts de novo nucleotide biosynthesis. While S. thermophilus had large contribution to the flavour profile, Lactococcus cremoris also played a role by limiting diacetyl and acetoin formation, which otherwise results in an off-flavour when in excess. This off-flavour control could be attributed to the metabolic re-routing of citrate by L. cremoris from diacetyl and acetoin towards α-ketoglutarate. Further, closely related Lactococcus lactis strains exhibited different interaction patterns with S. thermophilus, highlighting the significance of strain specificity in cheese making. Our results highlight the crucial roles of competitive and cooperative microbial interactions in shaping cheese flavour profile.

Metabolic interactions affect the biomass of synthetic bacterial biofilm communities

IMPORTANCE

Co-occurrence network analysis is an effective tool for predicting complex networks of microbial interactions in the natural environment. Using isolates from a rhizosphere, we constructed multi-species biofilm communities and investigated co-occurrence patterns between microbial species in genome-scale metabolic models and in vitro experiments. According to our results, metabolic exchanges and resource competition may partially explain the co-occurrence network analysis results found in synthetic bacterial biofilm communities.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

ILIAD: a suite of automated Snakemake workflows for processing genomic data for downstream applications

BACKGROUND

Processing raw genomic data for downstream applications such as imputation, association studies, and modeling requires numerous third-party bioinformatics software tools. It is highly time-consuming and resource-intensive with computational demands and storage limitations that pose significant challenges that increase cost. The use of software tools independent of one another, in a disjointed stepwise fashion, increases the difficulty and sets forth higher error rates because of fragmented job executions in alignment, variant calling, and/or build conversion complications. As sequencing data availability grows, the ability for biologists to process it using stable, automated, and reproducible workflows is paramount as it significantly reduces the time to generate clean and reliable data.

RESULTS

The Iliad suite of genomic data workflows was developed to provide users with seamless file transitions from raw genomic data to a quality-controlled variant call format (VCF) file for downstream applications. Iliad benefits from the efficiency of the Snakemake best practices framework coupled with Singularity and Docker containers for repeatability, portability, and ease of installation. This feat is accomplished from the onset with download acquisitions of any raw data type (FASTQ, CRAM, IDAT) straight through to the generation of a clean merged data file that can combine any user-preferred datasets using robust programs such as BWA, Samtools, and BCFtools. Users can customize and direct their workflow with one straightforward configuration file. Iliad is compatible with Linux, MacOS, and Windows platforms and scalable from a local machine to a high-performance computing cluster.

CONCLUSION

Iliad offers automated workflows with optimized time and resource management that are comparable to other workflows available but generates analysis-ready VCF files from the most common datatypes using a single command. The storage footprint challenge of genomic data is overcome by utilizing temporary intermediate files before the final VCF is generated. This file is ready for use in imputation, genome-wide association study (GWAS) pipelines, high-throughput population genetics studies, select gene candidate studies, and more. Iliad was developed to be portable, compatible, scalable, robust, and repeatable with a simplistic setup, so biologists that are less familiar with programming can manage their own big data with this open-source suite of workflows.

Disease-specific loss of microbial cross-feeding interactions in the human gut

Many gut microorganisms critical to human health rely on nutrients produced by each other for survival; however, these cross-feeding interactions are still challenging to quantify and remain poorly characterized. Here, we introduce a Metabolite Exchange Score (MES) to quantify those interactions. Using metabolic models of prokaryotic metagenome-assembled genomes from over 1600 individuals, MES allows us to identify and rank metabolic interactions that are significantly affected by a loss of cross-feeding partners in 10 out of 11 diseases. When applied to a Crohn's disease case-control study, our approach identifies a lack of species with the ability to consume hydrogen sulfide as the main distinguishing microbiome feature of disease. We propose that our conceptual framework will help prioritize in-depth analyses, experiments and clinical targets, and that targeting the restoration of microbial cross-feeding interactions is a promising mechanism-informed strategy to reconstruct a healthy gut ecosystem.

Rethinking Aerobic Respiration in the Hyporheic Zone under Variation in Carbon and Nitrogen Stoichiometry

Hyporheic zones (HZs)─zones of groundwater-surface water mixing─are hotspots for dissolved organic matter (DOM) and nutrient cycling that can disproportionately impact aquatic ecosystem functions. However, the mechanisms affecting DOM metabolism through space and time in HZs remain poorly understood. To resolve this gap, we investigate a recently proposed theory describing trade-offs between carbon (C) and nitrogen (N) limitations as a key regulator of HZ metabolism. We propose that throughout the extent of the HZ, a single process like aerobic respiration (AR) can be limited by both DOM thermodynamics and N content due to highly variable C/N ratios over short distances (centimeter scale). To investigate this theory, we used a large flume, continuous optode measurements of dissolved oxygen (DO), and spatially and temporally resolved molecular analysis of DOM. Carbon and N limitations were inferred from changes in the elemental stoichiometric ratio. We show sequential, depth-stratified relationships of DO with DOM thermodynamics and organic N that change across centimeter scales. In the shallow HZ with low C/N, DO was associated with the thermodynamics of DOM, while deeper in the HZ with higher C/N, DO was associated with inferred biochemical reactions involving organic N. Collectively, our results suggest that there are multiple competing processes that limit AR in the HZ. Resolving this spatiotemporal variation could improve predictions from mechanistic models, either via more highly resolved grid cells or by representing AR colimitation by DOM thermodynamics and organic N.

Metabolic engineering of human gut microbiome: Recent developments and future perspectives

Many studies have demonstrated that the gut microbiota is associated with human health and disease. Manipulation of the gut microbiota, e.g. supplementation of probiotics, has been suggested to be feasible, but subject to limited therapeutic efficacy. To develop efficient microbiota-targeted diagnostic and therapeutic strategies, metabolic engineering has been applied to construct genetically modified probiotics and synthetic microbial consortia. This review mainly discusses commonly adopted strategies for metabolic engineering in the human gut microbiome, including the use of in silico, in vitro, or in vivo approaches for iterative design and construction of engineered probiotics or microbial consortia. Especially, we highlight how genome-scale metabolic models can be applied to advance our understanding of the gut microbiota. Also, we review the recent applications of metabolic engineering in gut microbiome studies as well as discuss important challenges and opportunities.

Practical sampling of constraint-based models: Optimized thinning boosts CHRR performance

Thinning is a sub-sampling technique to reduce the memory footprint of Markov chain Monte Carlo. Despite being commonly used, thinning is rarely considered efficient. For sampling constraint-based models, a highly relevant use-case in systems biology, we here demonstrate that thinning boosts computational and, thereby, sampling efficiencies of the widely used Coordinate Hit-and-Run with Rounding (CHRR) algorithm. By benchmarking CHRR with thinning with simplices and genome-scale metabolic networks of up to thousands of dimensions, we find a substantial increase in computational efficiency compared to unthinned CHRR, in our examples by orders of magnitude, as measured by the effective sample size per time (ESS/t), with performance gains growing with polytope (effective network) dimension. Using a set of benchmark models we derive a ready-to-apply guideline for tuning thinning to efficient and effective use of compute resources without requiring additional coding effort. Our guideline is validated using three (out-of-sample) large-scale networks and we show that it allows sampling convex polytopes uniformly to convergence in a fraction of time, thereby unlocking the rigorous investigation of hitherto intractable models. The derivation of our guideline is explained in detail, allowing future researchers to update it as needed as new model classes and more training data becomes available. CHRR with deliberate utilization of thinning thereby paves the way to keep pace with progressing model sizes derived with the constraint-based reconstruction and analysis (COBRA) tool set. Sampling and evaluation pipelines are available at https://jugit.fz-juelich.de/IBG-1/ModSim/fluxomics/chrrt.

Food and drug design for gut microbiota-directed regulation: Current experimental landscape and future innovation

Most diets and medications enhance host health via microbiota-dependent ways, but it is in the present situation of untargeted regulation. Non-targeted regulation may lead to the ineffectiveness of dietary supplements or drug treatment. Microbiota-directed food, aiming to improve diseases by targeting specific microbes without affecting other bacteria, have been proposed to deal with this problem. However, there is currently no universally applicable method to explore such foods or drugs. In this review, thirty studies on recent efforts in microbiota directed diets and medications are summarized from various databases. The methods used to find new foods and medications are primarily divided into four groups depending on the experimental models: in vivo and in vitro, as well as predictions based on bioinformatics. We also discuss their implementation, interpretation, and respective limitations, and describe the present situation. We further put forward a framework for microbiota-directed foods and medicine according to above methods and other microbiome manipulation, which will spur precision medicine.

Flux balance analysis-based metabolic modeling of microbial secondary metabolism: Current status and outlook

In microorganisms, different from primary metabolism for cellular growth, secondary metabolism is for ecological interactions and stress responses and an important source of natural products widely used in various areas such as pharmaceutics and food additives. With advancements of sequencing technologies and bioinformatics tools, a large number of biosynthetic gene clusters of secondary metabolites have been discovered from microbial genomes. However, due to challenges from the difficulty of genome-scale pathway reconstruction and the limitation of conventional flux balance analysis (FBA) on secondary metabolism, the quantitative modeling of secondary metabolism is poorly established, in contrast to that of primary metabolism. This review first discusses current efforts on the reconstruction of secondary metabolic pathways in genome-scale metabolic models (GSMMs), as well as related FBA-based modeling techniques. Additionally, potential extensions of FBA are suggested to improve the prediction accuracy of secondary metabolite production. As this review posits, biosynthetic pathway reconstruction for various secondary metabolites will become automated and a modeling framework capturing secondary metabolism onset will enhance the predictive power. Expectedly, an improved FBA-based modeling workflow will facilitate quantitative study of secondary metabolism and in silico design of engineering strategies for natural product production.

Recent advances in data- and knowledge-driven approaches to explore primary microbial metabolism

With the rapid progress in metabolomics and sequencing technologies, more data on the metabolome of single microbes and their communities become available, revealing the potential of microorganisms to metabolize a broad range of chemical compounds. The analysis of microbial metabolomics datasets remains challenging since it inherits the technical challenges of metabolomics analysis, such as compound identification and annotation, while harboring challenges in data interpretation, such as distinguishing metabolite sources in mixed samples. This review outlines the recent advances in computational methods to analyze primary microbial metabolism: knowledge-based approaches that take advantage of metabolic and molecular networks and data-driven approaches that employ machine/deep learning algorithms in combination with large-scale datasets. These methods aim at improving metabolite identification and disentangling reciprocal interactions between microbes and metabolites. We also discuss the perspective of combining these approaches and further developments required to advance the investigation of primary metabolism in mixed microbial samples.

Environmental insults and compensative responses: when microbiome meets cancer

Tumor microenvironment has recently been ascribed a new hallmark-the polymorphic microbiome. Accumulating evidence regarding the tissue specific territories of tumor-microbiome have opened new and interesting avenues. A pertinent question is regarding the functional consequence of the interface between host-microbiome and cancer. Given microbial communities have predominantly been explored through an ecological perspective, it is important that the foundational aspects of ecological stress and the fight to 'survive and thrive' are accounted for tumor-micro(b)environment as well. Building on existing evidence and classical microbial ecology, here we attempt to characterize the ecological stresses and the compensative responses of the microorganisms inside the tumor microenvironment. What insults would microbes experience inside the cancer jungle? How would they respond to these insults? How the interplay of stress and microbial quest for survival would influence the fate of tumor? This work asks these questions and tries to describe this underdiscussed ecological interface of the tumor and its microbiota. It is hoped that a larger scientific thought on the importance of microbial competition sensing vis-à-vis tumor-microenvironment would be stimulated.

Moving beyond DNA: towards functional analysis of the vaginal microbiome by non-sequencing-based methods

Over the last two decades, sequencing-based methods have revolutionised our understanding of niche-specific microbial complexity. In the lower female reproductive tract, these approaches have enabled identification of bacterial compositional structures associated with health and disease. Application of metagenomics and metatranscriptomics strategies have provided insight into the putative function of these communities but it is increasingly clear that direct measures of microbial and host cell function are required to understand the contribution of microbe-host interactions to pathophysiology. Here we explore and discuss current methods and approaches, many of which rely upon mass-spectrometry, being used to capture functional insight into the vaginal mucosal interface. In addition to improving mechanistic understanding, these methods offer innovative solutions for the development of diagnostic and therapeutic strategies designed to improve women's health.

Microbes with higher metabolic independence are enriched in human gut microbiomes under stress

A wide variety of human diseases are associated with loss of microbial diversity in the human gut, inspiring a great interest in the diagnostic or therapeutic potential of the microbiota. However, the ecological forces that drive diversity reduction in disease states remain unclear, rendering it difficult to ascertain the role of the microbiota in disease emergence or severity. One hypothesis to explain this phenomenon is that microbial diversity is diminished as disease states select for microbial populations that are more fit to survive environmental stress caused by inflammation or other host factors. Here, we tested this hypothesis on a large scale, by developing a software framework to quantify the enrichment of microbial metabolisms in complex metagenomes as a function of microbial diversity. We applied this framework to over 400 gut metagenomes from individuals who are healthy or diagnosed with inflammatory bowel disease (IBD). We found that high metabolic independence (HMI) is a distinguishing characteristic of microbial communities associated with individuals diagnosed with IBD. A classifier we trained using the normalized copy numbers of 33 HMI-associated metabolic modules not only distinguished states of health versus IBD, but also tracked the recovery of the gut microbiome following antibiotic treatment, suggesting that HMI is a hallmark of microbial communities in stressed gut environments.

Flux Sampling in Genome-scale Metabolic Modeling of Microbial Communities

Microbial communities play a crucial role in ecosystem function through metabolic interactions. Genome-scale modeling is a promising method to understand these interactions. 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 may capture additional heterogeneity across cells, especially when cells exhibit sub-maximal growth rates. In this study, we simulate the metabolism of microbial communities and compare the metabolic characteristics found with FBA and flux sampling. We find significant differences in the predicted metabolism with sampling, including increased cooperative interactions and pathway-specific changes in predicted flux. Our results suggest the importance of sampling-based and objective function-independent approaches to evaluate metabolic interactions and emphasize their utility in quantitatively studying interactions between cells and organisms.

The media composition as a crucial element in high-throughput metabolic network reconstruction

In recent years, metagenome-assembled genomes (MAGs) have provided glimpses into the intra- and interspecies genetic diversity and interactions that form the bases of complex microbial communities. High-throughput reconstruction of genome-scale metabolic networks (GEMs) from MAGs is a promising avenue to disentangle the myriad trophic interactions stabilizing these communities. However, high-throughput reconstruction of GEMs relies on accurate gap filling of metabolic pathways using automated algorithms. Here, we systematically explore how the composition of the media (specification of the available nutrients and metabolites) during gap filling influences the resulting GEMs concerning predicted auxotrophies for fully sequenced model organisms and environmental isolates. We expand this analysis by using 106 MAGs from the same species with differing quality. We find that although the completeness of MAGs influences the fraction of gap-filled reactions, the composition of the media plays the dominant role in the accurate prediction of auxotrophies that form the basis of myriad community interactions. We propose that constraining the media composition for gap filling through both experimental approaches and computational approaches will increase the reliability of high-throughput reconstruction of genome-scale metabolic models from MAGs and paves the way for culture independent prediction of trophic interactions in complex microbial communities.

Metatranscriptomics-guided genome-scale metabolic modeling of microbial communities

Multi-omics data integration via mechanistic models of metabolism is a scalable and flexible framework for exploring biological hypotheses in microbial systems. However, although most microorganisms are unculturable, such multi-omics modeling is limited to isolate microbes or simple synthetic communities. Here, we developed an approach for modeling microbial activity and interactions that leverages the reconstruction of metagenome-assembled genomes and associated genome-centric metatranscriptomes. At its core, we designed a method for condition-specific metabolic modeling of microbial communities through the integration of metatranscriptomic data. Using this approach, we explored the behavior of anaerobic digestion consortia driven by hydrogen availability and human gut microbiota dysbiosis associated with Crohn's disease, identifying condition-dependent amino acid requirements in archaeal species and a reduced short-chain fatty acid exchange network associated with disease, respectively. Our approach can be applied to complex microbial communities, allowing a mechanistic contextualization of multi-omics data on a metagenome scale.

In silico prediction of potential indigenous microbial biomarkers in Penaeus vannamei identified through meta-analysis and genome-scale metabolic modelling

BACKGROUND

Understanding the microbiome is crucial as it contributes to the metabolic health of the host and, upon dysbiosis, may influence disease development. With the recent surge in high-throughput sequencing technology, the availability of microbial genomic data has increased dramatically. Amplicon sequence-based analyses majorly profile microbial abundance and determine taxonomic markers. Furthermore, the availability of genome sequences for various microbial organisms has prompted the integration of genome-scale metabolic modelling that provides insights into the metabolic interactions influencing host health. However, the analysis from a single study may not be consistent, necessitating a meta-analysis.

RESULTS

We conducted a meta-analysis and integrated with constraint-based metabolic modelling approach, focusing on the microbiome of pacific white shrimp Penaeus vannamei, an extensively cultured marine candidate species. Meta-analysis revealed that Acinetobacter and Alteromonas are significant indicators of "health" and "disease" specific taxonomic biomarkers, respectively. Further, we enumerated metabolic interactions among the taxonomic biomarkers by applying a constraint-based approach to the community metabolic models (4416 pairs). Under different nutrient environments, a constraint-based flux simulation identified five beneficial species: Acinetobacter spWCHA55, Acinetobacter tandoii SE63, Bifidobacterium pseudolongum 49 D6, Brevundimonas pondensis LVF1, and Lutibacter profundi LP1 mediating parasitic interactions majorly under sucrose environment in the pairwise community. The study also reports the healthy biomarkers that can co-exist and have functionally dependent relationships to maintain a healthy state in the host.

CONCLUSIONS

Toward this, we collected and re-analysed the amplicon sequence data of P. vannamei (encompassing 117 healthy and 142 disease datasets). By capturing the taxonomic biomarkers and modelling the metabolic interaction between them, our study provides a valuable resource, a first-of-its-kind analysis in aquaculture scenario toward a sustainable shrimp farming.

Interspecies Metabolic Interactions in a Synergistic Consortium Drive Efficient Degradation of the Herbicide Bromoxynil Octanoate

Microbial communities play vital roles in biogeochemical cycles, allowing biodegradation of a wide range of pollutants. Although many studies have shown the importance of interspecies interactions on activities of communities, fully elucidating the complex interactions in microbial communities is still challenging. Here, we isolated a consortium containing two bacterial strains (Acinetobacter sp. AG3 and Bacillus sp. R45), which could mineralize bromoxynil octanoate (BO) with higher efficiency than either strain individually. The BO degradation pathway by the synergistic consortium was elucidated, and interspecies interactions in the consortium were explored using genome-scale metabolic models (GSMMs). Modeling showed that growth and degradation enhancements were driven by metabolic interactions, such as syntrophic exchanges of small metabolites in the consortium. Besides, nutritional enhancers were predicted to improve BO degradation, which were tested experimentally. Overall, our results will enhance our understanding of microbial mineralization of BO by consortia and promote the application of microbial communities for bioremediation.

Competition-cooperation in the chemoautotrophic ecosystem of Movile Cave: first metagenomic approach on sediments

BACKGROUND

Movile Cave (SE Romania) is a chemoautotrophically-based ecosystem fed by hydrogen sulfide-rich groundwater serving as a primary energy source analogous to the deep-sea hydrothermal ecosystems. Our current understanding of Movile Cave microbiology has been confined to the sulfidic water and its proximity, as most studies focused on the water-floating microbial mat and planktonic accumulations likely acting as the primary production powerhouse of this unique subterranean ecosystem. By employing comprehensive genomic-resolved metagenomics, we questioned the spatial variation, chemoautotrophic abilities, ecological interactions and trophic roles of Movile Cave's microbiome thriving beyond the sulfidic-rich water.

RESULTS

A customized bioinformatics pipeline led to the recovery of 106 high-quality metagenome-assembled genomes from 7 cave sediment metagenomes. Assemblies' taxonomy spanned 19 bacterial and three archaeal phyla with Acidobacteriota, Chloroflexota, Proteobacteria, Planctomycetota, Ca. Patescibacteria, Thermoproteota, Methylomirabilota, and Ca. Zixibacteria as prevalent phyla. Functional gene analyses predicted the presence of CO₂ fixation, methanotrophy, sulfur and ammonia oxidation in the explored sediments. Species Metabolic Coupling Analysis of metagenome-scale metabolic models revealed the highest competition-cooperation interactions in the sediments collected away from the water. Simulated metabolic interactions indicated autotrophs and methanotrophs as major donors of metabolites in the sediment communities. Cross-feeding dependencies were assumed only towards 'currency' molecules and inorganic compounds (O₂, PO₄^3-, H⁺, Fe²⁺, Cu²⁺) in the water proximity sediment, whereas hydrogen sulfide and methanol were assumedly traded exclusively among distant gallery communities.

CONCLUSIONS

These findings suggest that the primary production potential of Movile Cave expands way beyond its hydrothermal waters, enhancing our understanding of the functioning and ecological interactions within chemolithoautotrophically-based subterranean ecosystems.

Ecological modelling approaches for predicting emergent properties in microbial communities

Recent studies have brought forward the critical role of emergent properties in shaping microbial communities and the ecosystems of which they are a part. Emergent properties-patterns or functions that cannot be deduced linearly from the properties of the constituent parts-underlie important ecological characteristics such as resilience, niche expansion and spatial self-organization. While it is clear that emergent properties are a consequence of interactions within the community, their non-linear nature makes mathematical modelling imperative for establishing the quantitative link between community structure and function. As the need for conservation and rational modulation of microbial ecosystems is increasingly apparent, so is the consideration of the benefits and limitations of the approaches to model emergent properties. Here we review ecosystem modelling approaches from the viewpoint of emergent properties. We consider the scope, advantages and limitations of Lotka-Volterra, consumer-resource, trait-based, individual-based and genome-scale metabolic models. Future efforts in this research area would benefit from capitalizing on the complementarity between these approaches towards enabling rational modulation of complex microbial ecosystems.

Toward merging bottom-up and top-down model-based designing of synthetic microbial communities

The increasing interest of microbial communities as promising biocatalyst is leading an intense effort into the development of computational frameworks assisting the analysis and rational engineering of such complex ecosystems. Here, we critically review the recent computational and model-guided advances in the system-level engineering of microbiome, including both the rational bottom-up and the evolutionary top-down approaches. Furthermore, we highlight modeling and computational methods supporting both engineering paradigms. Finally, we discuss the advantages of combining both strategies into a hybrid top-down/bottom-up (middle-out) strategy to engineer synthetic microbial communities with improved performance and scope.

Cancer's second genome: Microbial cancer diagnostics and redefining clonal evolution as a multispecies process: Humans and their tumors are not aseptic, and the multispecies nature of cancer modulates clinical care and clonal evolution: Humans and their tumors are not aseptic, and the multispecies nature of cancer modulates clinical care and clonal evolution

The presence and role of microbes in human cancers has come full circle in the last century. Tumors are no longer considered aseptic, but implications for cancer biology and oncology remain underappreciated. Opportunities to identify and build translational diagnostics, prognostics, and therapeutics that exploit cancer's second genome-the metagenome-are manifold, but require careful consideration of microbial experimental idiosyncrasies that are distinct from host-centric methods. Furthermore, the discoveries of intracellular and intra-metastatic cancer bacteria necessitate fundamental changes in describing clonal evolution and selection, reflecting bidirectional interactions with non-human residents. Reconsidering cancer clonality as a multispecies process similarly holds key implications for understanding metastasis and prognosing therapeutic resistance while providing rational guidance for the next generation of bacterial cancer therapies. Guided by these new findings and challenges, this Review describes opportunities to exploit cancer's metagenome in oncology and proposes an evolutionary framework as a first step towards modeling multispecies cancer clonality. Also see the video abstract here: https://youtu.be/-WDtIRJYZSs.

The hidden simplicity of metabolic networks is revealed by multireaction dependencies

Understanding the complexity of metabolic networks has implications for manipulation of their functions. The complexity of metabolic networks can be characterized by identifying multireaction dependencies that are challenging to determine due to the sheer number of combinations to consider. Here, we propose the concept of concordant complexes that captures multireaction dependencies and can be efficiently determined from the algebraic structure and operational constraints of metabolic networks. The concordant complexes imply the existence of concordance modules based on which the apparent complexity of 12 metabolic networks of organisms from all kingdoms of life can be reduced by at least 78%. A comparative analysis against an ensemble of randomized metabolic networks shows that the metabolic network of Escherichia coli contains fewer concordance modules and is, therefore, more tightly coordinated than expected by chance. Together, our findings demonstrate that metabolic networks are considerably simpler than what can be perceived from their structure alone.

Systems Biology of Gut Microbiota-Human Receptor Interactions: Toward Anti-inflammatory Probiotics

The incidence and prevalence of inflammatory disorders have increased globally, and is projected to double in the next decade. Gut microbiome-based therapeutics have shown promise in ameliorating chronic inflammation. However, they are largely experimental, context- or strain-dependent and lack a clear mechanistic basis. This hinders precision probiotics and poses significant risk, especially to individuals with pre-existing conditions. Molecules secreted by gut microbiota act as ligands to several health-relevant receptors expressed in human gut, such as the G-protein coupled receptors (GPCRs), Toll-like receptor 4 (TLR4), pregnane X receptor (PXR), and aryl hydrocarbon receptor (AhR). Among these, the human AhR expressed in different tissues exhibits anti-inflammatory effects and shows activity against a wide range of ligands produced by gut bacteria. However, different AhR ligands induce varying host responses and signaling in a tissue/organ-specific manner, which remain mostly unknown. The emerging systems biology paradigm, with its powerful in silico tool repertoire, provides opportunities for comprehensive and high-throughput strain characterization. In particular, combining metabolic models with machine learning tools can be useful to delineate tissue and ligand-specific signaling and thus their causal mechanisms in disease and health. The knowledge of such a mechanistic basis is indispensable to account for strain heterogeneity and actualize precision probiotics.

KEMET - a python tool for KEGG Module evaluation and microbial genome annotation expansion

•

KEMET allows evaluating and expanding microbial genome annotations.

•

KEMET can identify unannotated genes having partial sequence truncations.

•

Annotation expansion can improve contextual draft genome-scale metabolic model reconstruction.

Background

The rapid accumulation of sequencing data from metagenomic studies is enabling the generation of huge collections of microbial genomes, with new challenges for mapping their functional potential. In particular, metagenome-assembled genomes are typically incomplete and harbor partial gene sequences that can limit their annotation from traditional tools. New scalable solutions are thus needed to facilitate the evaluation of functional potential in microbial genomes.

Methods

To resolve annotation gaps in microbial genomes, we developed KEMET, an open-source Python library devised for the analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) functional units. KEMET focuses on the in-depth analysis of metabolic reaction networks to identify missing orthologs through hidden Markov model profiles.

Results

We evaluate the potential of KEMET for expanding functional annotations by simulating the effect of assembly issues on real gene sequences and showing that our approach can identify missing KEGG orthologs. Additionally, we show that recovered gene annotations can sensibly increase the quality of draft genome-scale metabolic models obtained from metagenome-assembled genomes, in some cases reaching the accuracy of models generated from complete genomes.

Conclusions

KEMET therefore allows expanding genome annotations by targeted searches for orthologous sequences, enabling a better qualitative and quantitative assessment of metabolic capabilities in novel microbial organisms.

METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks

BACKGROUND

Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling.

RESULTS

We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric "MW-score" (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut.

CONCLUSION

METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC . Video abstract.

Plastic-Degrading Potential across the Global Microbiome Correlates with Recent Pollution Trends

Biodegradation is a plausible route toward sustainable management of the millions of tons of plastic waste that have accumulated in terrestrial and marine environments. However, the global diversity of plastic-degrading enzymes remains poorly understood. Taking advantage of global environmental DNA sampling projects, here we constructed hidden Markov models from experimentally verified enzymes and mined ocean and soil metagenomes to assess the global potential of microorganisms to degrade plastics. By controlling for false positives using gut microbiome data, we compiled a catalogue of over 30,000 nonredundant enzyme homologues with the potential to degrade 10 different plastic types. While differences between the ocean and soil microbiomes likely reflect the base compositions of these environments, we find that ocean enzyme abundance increases with depth as a response to plastic pollution and not merely taxonomic composition. By obtaining further pollution measurements, we observed that the abundance of the uncovered enzymes in both ocean and soil habitats significantly correlates with marine and country-specific plastic pollution trends. Our study thus uncovers the earth microbiome's potential to degrade plastics, providing evidence of a measurable effect of plastic pollution on the global microbial ecology as well as a useful resource for further applied research. IMPORTANCE Utilization of synthetic biology approaches to enhance current plastic degradation processes is of crucial importance, as natural plastic degradation processes are very slow. For instance, the predicted lifetime of a polyethylene terephthalate (PET) bottle under ambient conditions ranges from 16 to 48 years. Moreover, although there is still unexplored diversity in microbial communities, synergistic degradation of plastics by microorganisms holds great potential to revolutionize the management of global plastic waste. To this end, the methods and data on novel plastic-degrading enzymes presented here can help researchers by (i) providing further information about the taxonomic diversity of such enzymes as well as understanding of the mechanisms and steps involved in the biological breakdown of plastics, (ii) pointing toward the areas with increased availability of novel enzymes, and (iii) giving a basis for further application in industrial plastic waste biodegradation. Importantly, our findings provide evidence of a measurable effect of plastic pollution on the global microbial ecology.