Search

We review the process of reconstructing genome-scale metabolic models (GEMs) for single species and microbial communities. We present a detailed compendium of model-building platforms for experts and non-experts. We review and compare the approaches for testing GEMs towards engineering microbial systems. We discuss current challenges and future directions in model development, including the need for standards in microbial community models. Genome-scale metabolic models (GEMs) are consolidating as platforms for studying mixed microbial populations, by combining biological data and knowledge with mathematical rigor. However, deploying these models to answer research questions can be challenging due to the increasing number of available computational tools, the lack of universal standards, and their inherent limitations. Here, we present a comprehensive overview of foundational concepts for building and evaluating genome-scale models of microbial communities. We then compare tools in terms of requirements, capabilities, and applications. Next, we highlight the current pitfalls and open challenges to consider when adopting existing tools and developing new ones. Our compendium can be relevant for the expanding community of modelers, both at the entry and experienced levels. [ABSTRACT FROM AUTHOR]

Multi-strain probiotics are widely regarded as effective products for improving gut microbiota stability and host health, providing advantages over single-strain probiotics. However, in general, it is unclear to what extent different strains would cooperate or compete for resources, and how the establishment of a common biofilm microenvironment could influence their interactions. In this work, we develop an integrative experimental and computational approach to comprehensively assess the metabolic functionality and interactions of probiotics across growth conditions. Our approach combines co-culture assays with genome-scale modelling of metabolism and multivariate data analysis, thus exploiting complementary data- and knowledge-driven systems biology techniques. To show the advantages of the proposed approach, we apply it to the study of the interactions between two widely used probiotic strains of Lactobacillus reuteri and Saccharomyces boulardii , characterising their production potential for compounds that can be beneficial to human health. Our results show that these strains can establish a mixed cooperative-antagonistic interaction best explained by competition for shared resources, with an increased individual exchange but an often decreased net production of amino acids and short-chain fatty acids. Overall, our work provides a strategy that can be used to explore microbial metabolic fingerprints of biotechnological interest, capable of capturing multifaceted equilibria even in simple microbial consortia. • We propose an approach to predict interactions of probiotics across growth conditions. • Metabolism of L. reuteri and S. boulardii is studied in vitro and in silico. • Multi-strain probiotics can establish a mixed cooperative-antagonistic interaction. • Inter-species interaction are characterised in planktonic and biofilm growth. • In silico models well predict amino acid production across cultures. [ABSTRACT FROM AUTHOR]

Genetic heterogeneity is a common trait in microbial populations, caused by de novo mutations and changes in variant frequencies over time. Microbes can thus differ genetically within the same species and acquire different phenotypes. For instance, performance and stability of anaerobic reactors are linked to the composition of the microbiome involved in the digestion process and to the environmental parameters imposing selective pressure on the metagenome, shaping its evolution. Changes at the strain level have the potential to determine variations in microbial functions, and their characterization could provide new insight into ecological and evolutionary processes driving anaerobic digestion. In this work, single nucleotide variant dynamics were studied in two time-course biogas upgrading experiments, testing alternative carbon sources and the response to exogenous hydrogen addition. A cumulative total of 76,229 and 64,289 high-confidence single nucleotide variants were discerned in the experiments related to carbon substrate availability and hydrogen addition, respectively. By combining complementary bioinformatic approaches, the study reconstructed the precise strain count—two for both hydrogenotrophic archaea—and tracked their abundance over time, while also characterizing tens of genes under strong selection. Results in the dominant archaea revealed the presence of nearly 100 variants within genes encoding enzymes involved in hydrogenotrophic methanogenesis. In the bacterial counterparts, 119 mutations were identified across 23 genes associated with the Wood-Ljungdahl pathway, suggesting a possible impact on the syntrophic acetate-oxidation process. Strain replacement events took place in both experiments, confirming the trends suggested by the variants trajectories and providing a comprehensive understanding of the biogas upgrading microbiome at the strain level. Overall, this resolution level allowed us to reveal fine-scale evolutionary mechanisms, functional dynamics, and strain-level metabolic variation that could contribute to the selection of key species actively involved in the carbon dioxide fixation process. [Display omitted] • This is the first use of strain-resolved metagenomics in in-situ biogas upgrading. • A total of 97 high-quality metagenome assembled genomes were reconstructed. • The variants calling approach retrieved 140,518 variants in high abundance genomes. • Two strains were found in M. wolfeii and M. thermophilus at distinct time points. • Eight hydrogenotrophic methanogenesis genes were under selective pressure. [ABSTRACT FROM AUTHOR]

Motivation High-throughput biological data, thanks to technological advances, have become cheaper to collect, leading to the availability of vast amounts of omic data of different types. In parallel, the in silico reconstruction and modeling of metabolic systems is now acknowledged as a key tool to complement experimental data on a large scale. The integration of these model- and data-driven information is therefore emerging as a new challenge in systems biology, with no clear guidance on how to better take advantage of the inherent multisource and multiomic nature of these data types while preserving mechanistic interpretation. Results Here, we investigate different regularization techniques for high-dimensional data derived from the integration of gene expression profiles with metabolic flux data, extracted from strain-specific metabolic models, to improve cellular growth rate predictions. To this end, we propose ad-hoc extensions of previous regularization frameworks including group, view-specific and principal component regularization and experimentally compare them using data from 1143 Saccharomyces cerevisiae strains. We observe a divergence between methods in terms of regression accuracy and integration effectiveness based on the type of regularization employed. In multiomic regression tasks, when learning from experimental and model-generated omic data, our results demonstrate the competitiveness and ease of interpretation of multimodal regularized linear models compared to data-hungry methods based on neural networks. Availability and implementation All data, models and code produced in this work are available on GitHub at https://github.com/Angione-Lab/HybridGroupIPFLasso%5fpc2Lasso. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Cancer is considered a high‐risk condition for severe illness resulting from COVID‐19. The interaction between severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) and human metabolism is key to elucidating the risk posed by COVID‐19 for cancer patients and identifying effective treatments, yet it is largely uncharacterised on a mechanistic level. We present a genome‐scale map of short‐term metabolic alterations triggered by SARS‐CoV‐2 infection of cancer cells. Through transcriptomic‐ and proteomic‐informed genome‐scale metabolic modelling, we characterise the role of RNA and fatty acid biosynthesis in conjunction with a rewiring in energy production pathways and enhanced cytokine secretion. These findings link together complementary aspects of viral invasion of cancer cells, while providing mechanistic insights that can inform the development of treatment strategies. [ABSTRACT FROM AUTHOR]

The use of data mining and modeling methods in service industry is a promising avenue for optimizing current processes in a targeted manner, ultimately reducing costs and improving customer experience. However, the introduction of such tools in already established pipelines often must adapt to the way data is sampled and to its content. In this study, we tackle the challenge of characterizing and predicting customer experience having available only process log data with time-stamp information, without any ground truth feedback from the customers. As a case study, we consider the context of a contact center managed by TeleWare and analyze phone call logs relative to a two months span. We develop an approach to interpret the phone call process events registered in the logs and infer concrete points of improvement in the service management. Our approach is based on latent tree modeling and multi-class Naïve Bayes classification, which jointly allow us to infer a spectrum of customer experiences and test their predictability based on the current data sampling strategy. Moreover, such approach can overcome limitations in customer feedback collection and sharing across organizations, thus having wide applicability and being complementary to tools relying on more heavily constrained data. [ABSTRACT FROM AUTHOR]

Anaerobic and microaerophilic environments are pervasive in nature, providing essential contributions to the maintenance of human health, biogeochemical cycles and the Earth's climate. These ecological niches are characterised by low free oxygen and oxidants, or lack thereof. Under these conditions, interactions between species are essential for supporting the growth of syntrophic species and maintaining thermodynamic feasibility of anaerobic fermentation. Kinetic models provide a simplified view of complex metabolic networks, while genome-scale metabolic models and flux-balance analysis (FBA) aim to unravel these systems as a whole. The target of this review is to outline the main similarities, differences and challenges associated with kinetic and metabolic modelling, and describe state-of-the-art modelling practices for studying syntrophies in the anaerobic digestion (AD) case study. [Display omitted] • Anaerobic microbes play a significant role in the global biogeochemical cycle. • Mathematical models explain anaerobic processes, including bioconversion abilities, growth requirements and inhibitor susceptibility. • Kinetic modelling and FBA, are outlined, highlighting strengths and limitations. • The AD system is used as case of study to discuss modelling of syntrophies. [ABSTRACT FROM AUTHOR]

Omic data analysis is steadily growing as a driver of basic and applied molecular biology research. Core to the interpretation of complex and heterogeneous biological phenotypes are computational approaches in the fields of statistics and machine learning. In parallel, constraint-based metabolic modeling has established itself as the main tool to investigate large-scale relationships between genotype, phenotype, and environment. The development and application of these methodological frameworks have occurred independently for the most part, whereas the potential of their integration for biological, biomedical, and biotechnological research is less known. Here, we describe how machine learning and constraint-based modeling can be combined, reviewing recent works at the intersection of both domains and discussing the mathematical and practical aspects involved. We overlap systematic classifications from both frameworks, making them accessible to nonexperts. Finally, we delineate potential future scenarios, propose new joint theoretical frameworks, and suggest concrete points of investigation for this joint subfield. A multiview approach merging experimental and knowledge-driven omic data through machine learning methods can incorporate key mechanistic information in an otherwise biologically-agnostic learning process. [ABSTRACT FROM AUTHOR]

Background: The uncovering of genes linked to human diseases is a pressing challenge in molecular biology and precision medicine. This task is often hindered by the large number of candidate genes and by the heterogeneity of the available information. Computational methods for the prioritization of candidate genes can help to cope with these problems. In particular, kernel-based methods are a powerful resource for the integration of heterogeneous biological knowledge, however, their practical implementation is often precluded by their limited scalability. Results: We propose Scuba, a scalable kernel-based method for gene prioritization. It implements a novel multiple kernel learning approach, based on a semi-supervised perspective and on the optimization of the margin distribution. Scuba is optimized to cope with strongly unbalanced settings where known disease genes are few and large scale predictions are required. Importantly, it is able to efficiently deal both with a large amount of candidate genes and with an arbitrary number of data sources. As a direct consequence of scalability, Scuba integrates also a new efficient strategy to select optimal kernel parameters for each data source. We performed cross-validation experiments and simulated a realistic usage setting, showing that Scuba outperforms a wide range of state-of-the-art methods. Conclusions: Scuba achieves state-of-the-art performance and has enhanced scalability compared to existing kernel-based approaches for genomic data. This method can be useful to prioritize candidate genes, particularly when their number is large or when input data is highly heterogeneous. [ABSTRACT FROM AUTHOR]

Background: The anaerobic digestion process degrades organic matter into simpler compounds and occurs in strictly anaerobic and microaerophilic environments. The process is carried out by a diverse community of microorganisms where each species has a unique role and it has relevant biotechnological applications since it is used for biogas production. Some aspects of the microbiome, including its interaction with phages, remains still unclear: a better comprehension of the community composition and role of each species is crucial for a cured understanding of the carbon cycle in anaerobic systems and improving biogas production. Results: The primary objective of this study was to expand our understanding on the anaerobic digestion microbiome by jointly analyzing its prokaryotic and viral components. By integrating 192 additional datasets into a previous metagenomic database, the binning process generated 11,831 metagenome-assembled genomes from 314 metagenome samples published between 2014 and 2022, belonging to 4,568 non-redundant species based on ANI calculation and quality verification. CRISPR analysis on these genomes identified 76 archaeal genomes with active phage interactions. Moreover, single-nucleotide variants further pointed to archaea as the most critical members of the community. Among the MAGs, two methanogenic archaea, Methanothrix sp. 43zhSC_152 and Methanoculleus sp. 52maCN_3230, had the highest number of SNVs, with the latter having almost double the density of most other MAGs. Conclusions: This study offers a more comprehensive understanding of microbial community structures that thrive at different temperatures. The findings revealed that the fraction of archaeal species characterized at the genome level and reported in public databases is higher than that of bacteria, although still quite limited. The identification of shared spacers between phages and microbes implies a history of phage-bacterial interactions, and specifically lysogenic infections. A significant number of SNVs were identified, primarily comprising synonymous and nonsynonymous variants. Together, the findings indicate that methanogenic archaea are subject to intense selective pressure and suggest that genomic variants play a critical role in the anaerobic digestion process. Overall, this study provides a more balanced and diverse representation of the anaerobic digestion microbiota in terms of geographic location, temperature range and feedstock utilization. [ABSTRACT FROM AUTHOR]

Metabolic modeling and machine learning are key components in the emerging next generation of systems and synthetic biology tools, targeting the genotype-phenotype-environment relationship. Rather than being used in isolation, it is becoming clear that their value is maximized when they are combined. However, the potential of integrating these two frameworks for omic data augmentation and integration is largely unexplored. We propose, rigorously assess, and compare machine-learning-based data integration techniques, combining gene expression profiles with computationally generated metabolic flux data to predict yeast cell growth. To this end, we create strain-specific metabolic models for 1,143 Saccharomyces cerevisiae mutants and we test 27 machine-learning methods, incorporating state-of-the-art feature selection and multiview learning approaches. We propose a multiview neural network using fluxomic and transcriptomic data, showing that the former increases the predictive accuracy of the latter and reveals functional patterns that are not directly deducible from gene expression alone. We test the proposed neural network on a further 86 strains generated in a different experiment, therefore verifying its robustness to an additional independent dataset. Finally, we show that introducing mechanistic flux features improves the predictions also for knockout strains whose genes were not modeled in the metabolic reconstruction. Our results thus demonstrate that fusing experimental cues with in silico models, based on known biochemistry, can contribute with disjoint information toward biologically informed and interpretable machine learning. Overall, this study provides tools for understanding and manipulating complex phenotypes, increasing both the prediction accuracy and the extent of discernible mechanistic biological insights. [ABSTRACT FROM AUTHOR]

Motivation Gene regulation is responsible for controlling numerous physiological functions and dynamically responding to environmental fluctuations. Reconstructing the human network of gene regulatory interactions is thus paramount to understanding the cell functional organization across cell types, as well as to elucidating pathogenic processes and identifying molecular drug targets. Although significant effort has been devoted towards this direction, existing computational methods mainly rely on gene expression levels, possibly ignoring the information conveyed by mechanistic biochemical knowledge. Moreover, except for a few recent attempts, most of the existing approaches only consider the information of the organism under analysis, without exploiting the information of related model organisms. Results We propose a novel method for the reconstruction of the human gene regulatory network, based on a transfer learning strategy that synergically exploits information from human and mouse, conveyed by gene-related metabolic features generated in silico from gene expression data. Specifically, we learn a predictive model from metabolic activity inferred via tissue-specific metabolic modelling of artificial gene knockouts. Our experiments show that the combination of our transfer learning approach with the constructed metabolic features provides a significant advantage in terms of reconstruction accuracy, as well as additional clues on the contribution of each constructed metabolic feature. Availability and implementation The method, the datasets and all the results obtained in this study are available at: https://doi.org/10.6084/m9.figshare.c.5237687. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Cupriavidus necator is a bacterium with a high phenotypic diversity and versatile metabolic capabilities. It has been extensively studied as a model hydrogen oxidizer, as well as a producer of polyhydroxyalkanoates (PHA), plastic-like biopolymers with a high potential to substitute petroleum-based materials. Thanks to its adaptability to diverse metabolic lifestyles and to the ability to accumulate large amounts of PHA, C. necator is employed in many biotechnological processes, with particular focus on PHA production from waste carbon sources. The large availability of genomic information has enabled a characterization of C. necator 's metabolism, leading to the establishment of metabolic models which are used to devise and optimize culture conditions and genetic engineering approaches. In this work, the characteristics of available C. necator strains and genomes are reviewed, underlining how a thorough comprehension of the genetic variability of C. necator is lacking and it could be instrumental for wider application of this microorganism. The metabolic paradigms of C. necator and how they are connected to PHA production and accumulation are described, also recapitulating the variety of carbon substrates used for PHA accumulation, highlighting the most promising strategies to increase the yield. Finally, the review describes and critically analyzes currently available genome-scale metabolic models and reduced metabolic network applications commonly employed in the optimization of PHA production. Overall, it appears that the capacity of C. necator of performing CO 2 bioconversion to PHA is still underexplored, both in biotechnological applications and in metabolic modeling. However, the accurate characterization of this organism and the efforts in using it for gas fermentation can help tackle this challenging perspective in the future. • Taxonomy of Cupriavidus necator strains is examined by comparative genomics. • The metabolism of C. necator is suited to diverse aerobic and anaerobic lifestyles. • C. necator can convert many waste products to polyhydroxyalkanoates. • Metabolic modeling of C. necator is a powerful tool for bioprocess optimization. [ABSTRACT FROM AUTHOR]

Anaerobic digestion is a key biological process for renewable energy, yet the mechanistic knowledge on its hidden microbial dynamics is still limited. The present work charted the interaction network in the anaerobic digestion microbiome via the full characterization of pairwise interactions and the associated metabolite exchanges. To this goal, a novel collection of 836 genome-scale metabolic models was built to represent the functional capabilities of bacteria and archaea species derived from genome-centric metagenomics. Dominant microbes were shown to prefer mutualistic, parasitic and commensalistic interactions over neutralism, amensalism and competition, and are more likely to behave as metabolite importers and profiteers of the coexistence. Additionally, external hydrogen injection positively influences microbiome dynamics by promoting commensalism over amensalism. Finally, exchanges of glucogenic amino acids were shown to overcome auxotrophies caused by an incomplete tricarboxylic acid cycle. Our novel strategy predicted the most favourable growth conditions for the microbes, overall suggesting strategies to increasing the biogas production efficiency. In principle, this approach could also be applied to microbial populations of biomedical importance, such as the gut microbiome, to allow a broad inspection of the microbial interplays. Image 1 • Inspection of interactions among 836 microbial species with flux balance analysis. • Dominant microbes are often profiteers in the interaction. • Interactions based on amino acid exchanges are crucial for solving auxotrophies. • Positive interplays are more likely to occur interphylum than intraphylum. • External hydrogen supply decreases the rate of negative co-growth. [ABSTRACT FROM AUTHOR]