Resource-allocation constraint governs structure and function of microbial communities in metabolic modeling

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.

Metabolic model-based ecological modeling for probiotic design

The microbial community composition in the human gut has a profound effect on human health. This observation has lead to extensive use of microbiome therapies, including over-the-counter 'probiotic' treatments intended to alter the composition of the microbiome. Despite so much promise and commercial interest, the factors that contribute to the success or failure of microbiome-targeted treatments remain unclear. We investigate the biotic interactions that lead to successful engraftment of a novel bacterial strain introduced to the microbiome as in probiotic treatments. We use pairwise genome-scale metabolic modeling with a generalized resource allocation constraint to build a network of interactions between taxa that appear in an experimental engraftment study. We create induced sub-graphs using the taxa present in individual samples and assess the likelihood of invader engraftment based on network structure. To do so, we use a generalized Lotka-Volterra model, which we show has strong ability to predict if a particular invader or probiotic will successfully engraft into an individual's microbiome. Furthermore, we show that the mechanistic nature of the model is useful for revealing which microbe-microbe interactions potentially drive engraftment.

From metagenomes to metabolism: Systematically assessing the metabolic flux feasibilities for "Candidatus Accumulibacter" species during anaerobic substrate uptake

With the rapid growing availability of metagenome assembled genomes (MAGs) and associated metabolic models, the identification of metabolic potential in individual community members has become possible. However, the field still lacks an unbiassed systematic evaluation of the generated metagenomic information to uncover not only metabolic potential, but also feasibilities of these models under specific environmental conditions. In this study, we present a systematic analysis of the metabolic potential in species of "Candidatus Accumulibacter", a group of polyphosphate-accumulating organisms (PAOs). We constructed a metabolic model of the central carbon metabolism and compared the metabolic potential among available MAGs for "Ca. Accumulibacter" species. By combining Elementary Flux Modes Analysis (EFMA) with max-min driving force (MDF) optimization, we obtained all possible flux distributions of the metabolic network and calculated their individual thermodynamic feasibility. Our findings reveal significant variations in the metabolic potential among "Ca. Accumulibacter" MAGs, particularly in the presence of anaplerotic reactions. EFMA revealed 700 unique flux distributions in the complete metabolic model that enable the anaerobic uptake of acetate and its conversion into polyhydroxyalkanoates (PHAs), a well-known phenotype of "Ca. Accumulibacter". However, thermodynamic constraints narrowed down this solution space to 146 models that were stoichiometrically and thermodynamically feasible (MDF > 0 kJ/mol), of which only 8 were strongly feasible (MDF > 7 kJ/mol). Notably, several novel flux distributions for the metabolic model were identified, suggesting putative, yet unreported, functions within the PAO communities. Overall, this work provides valuable insights into the metabolic variability among "Ca. Accumulibacter" species and redefines the anaerobic metabolic potential in the context of phosphate removal. More generally, the integrated workflow presented in this paper can be applied to any metabolic model obtained from a MAG generated from microbial communities to objectively narrow the expected phenotypes from community members.

Top-down and bottom-up microbiome engineering approaches to enable biomanufacturing from waste biomass

Growing environmental concerns and the need to adopt a circular economy have highlighted the importance of waste valorization for resource recovery. Microbial consortia-enabled biotechnologies have made significant developments in the biomanufacturing of valuable resources from waste biomass that serve as suitable alternatives to petrochemical-derived products. These microbial consortia-based processes are designed following a top-down or bottom-up engineering approach. The top-down approach is a classical method that uses environmental variables to selectively steer an existing microbial consortium to achieve a target function. While high-throughput sequencing has enabled microbial community characterization, the major challenge is to disentangle complex microbial interactions and manipulate the structure and function accordingly. The bottom-up approach uses prior knowledge of the metabolic pathway and possible interactions among consortium partners to design and engineer synthetic microbial consortia. This strategy offers some control over the composition and function of the consortium for targeted bioprocesses, but challenges remain in optimal assembly methods and long-term stability. In this review, we present the recent advancements, challenges, and opportunities for further improvement using top-down and bottom-up approaches for microbiome engineering. As the bottom-up approach is relatively a new concept for waste valorization, this review explores the assembly and design of synthetic microbial consortia, ecological engineering principles to optimize microbial consortia, and metabolic engineering approaches for efficient conversion. Integration of top-down and bottom-up approaches along with developments in metabolic modeling to predict and optimize consortia function are also highlighted.

ONE-SENTENCE SUMMARY

This review highlights the microbial consortia-driven waste valorization for biomanufacturing through top-down and bottom-up design approaches and describes strategies, tools, and unexplored opportunities to optimize the design and stability of such consortia.

Combining compositional data sets introduces error in covariance network reconstruction

Microbial communities are diverse biological systems that include taxa from across multiple kingdoms of life. Notably, interactions between bacteria and fungi play a significant role in determining community structure. However, these statistical associations across kingdoms are more difficult to infer than intra-kingdom associations due to the nature of the data involved using standard network inference techniques. We quantify the challenges of cross-kingdom network inference from both theoretical and practical points of view using synthetic and real-world microbiome data. We detail the theoretical issue presented by combining compositional data sets drawn from the same environment, e.g. 16S and ITS sequencing of a single set of samples, and we survey common network inference techniques for their ability to handle this error. We then test these techniques for the accuracy and usefulness of their intra- and inter-kingdom associations by inferring networks from a set of simulated samples for which a ground-truth set of associations is known. We show that while the two methods mitigate the error of cross-kingdom inference, there is little difference between techniques for key practical applications including identification of strong correlations and identification of possible keystone taxa (i.e. hub nodes in the network). Furthermore, we identify a signature of the error caused by transkingdom network inference and demonstrate that it appears in networks constructed using real-world environmental microbiome data.

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.

Evaluating the prebiotic effect of oligosaccharides on gut microbiome wellness using in vitro fecal fermentation

We previously proposed the Gut Microbiome Wellness Index (GMWI), a predictor of disease presence based on a gut microbiome taxonomic profile. As an application of this index for food science research, we applied GMWI as a quantitative tool for measuring the prebiotic effect of oligosaccharides. Mainly, in an in vitro anaerobic batch fermentation system, fructooligosaccharides (FOS), galactooligosaccharides (GOS), xylooligosaccharides (XOS), inulin (IN), and 2'-fucosyllactose (2FL), were mixed separately with fecal samples obtained from healthy adult volunteers. To find out how 24 h prebiotic fermentation influenced the GMWI values in their respective microbial communities, changes in species-level relative abundances were analyzed in the five prebiotics groups, as well as in two control groups (no substrate addition at 0 h and for 24 h). The GMWI of fecal microbiomes treated with any of the five prebiotics (IN (0.48 ± 0.06) > FOS (0.47 ± 0.03) > XOS (0.33 ± 0.02) > GOS (0.26 ± 0.02) > 2FL (0.16 ± 0.06)) were positive, which indicates an increase of relative abundances of microbial species previously found to be associated with a healthy, disease-free state. In contrast, the GMWI of samples without substrate addition for 24 h (-0.60 ± 0.05) reflected a non-healthy, disease-harboring microbiome state. Compared to the original prebiotic index (PI) and α-diversity metrics, GMWI provides a more data-driven, evidence-based indexing system for evaluating the prebiotic effect of food components. This study demonstrates how GMWI can be applied as a novel PI in dietary intervention studies, with wider implications for designing personalized diets based on their impact on gut microbiome wellness.

TaxiBGC: a Taxonomy-Guided Approach for Profiling Experimentally Characterized Microbial Biosynthetic Gene Clusters and Secondary Metabolite Production Potential in Metagenomes

Biosynthetic gene clusters (BGCs) in microbial genomes encode bioactive secondary metabolites (SMs), which can play important roles in microbe-microbe and host-microbe interactions. Given the biological significance of SMs and the current profound interest in the metabolic functions of microbiomes, the unbiased identification of BGCs from high-throughput metagenomic data could offer novel insights into the complex chemical ecology of microbial communities. Currently available tools for predicting BGCs from shotgun metagenomes have several limitations, including the need for computationally demanding read assembly, predicting a narrow breadth of BGC classes, and not providing the SM product. To overcome these limitations, we developed taxonomy-guided identification of biosynthetic gene clusters (TaxiBGC), a command-line tool for predicting experimentally characterized BGCs (and inferring their known SMs) in metagenomes by first pinpointing the microbial species likely to harbor them. We benchmarked TaxiBGC on various simulated metagenomes, showing that our taxonomy-guided approach could predict BGCs with much-improved performance (mean F₁ score, 0.56; mean PPV score, 0.80) compared with directly identifying BGCs by mapping sequencing reads onto the BGC genes (mean F₁ score, 0.49; mean PPV score, 0.41). Next, by applying TaxiBGC on 2,650 metagenomes from the Human Microbiome Project and various case-control gut microbiome studies, we were able to associate BGCs (and their SMs) with different human body sites and with multiple diseases, including Crohn's disease and liver cirrhosis. In all, TaxiBGC provides an in silico platform to predict experimentally characterized BGCs and their SM production potential in metagenomic data while demonstrating important advantages over existing techniques. IMPORTANCE Currently available bioinformatics tools to identify BGCs from metagenomic sequencing data are limited in their predictive capability or ease of use to even computationally oriented researchers. We present an automated computational pipeline called TaxiBGC, which predicts experimentally characterized BGCs (and infers their known SMs) in shotgun metagenomes by first considering the microbial species source. Through rigorous benchmarking techniques on simulated metagenomes, we show that TaxiBGC provides a significant advantage over existing methods. When demonstrating TaxiBGC on thousands of human microbiome samples, we associate BGCs encoding bacteriocins with different human body sites and diseases, thereby elucidating a possible novel role of this antibiotic class in maintaining the stability of microbial ecosystems throughout the human body. Furthermore, we report for the first time gut microbial BGC associations shared among multiple pathologies. Ultimately, we expect our tool to facilitate future investigations into the chemical ecology of microbial communities across diverse niches and pathologies.