Traceability, reproducibility and wiki-exploration for "à-la-carte" reconstructions of genome-scale metabolic models

Modification and analysis of context-specific genome-scale metabolic models: methane-utilizing microbial chassis as a case study

Context-specific genome-scale model (CS-GSM) reconstruction is becoming an efficient strategy for integrating and cross-comparing experimental multi-scale data to explore the relationship between cellular genotypes, facilitating fundamental or applied research discoveries. However, the application of CS modeling for non-conventional microbes is still challenging. Here, we present a graphical user interface that integrates COBRApy, EscherPy, and RIPTiDe, Python-based tools within the BioUML platform, and streamlines the reconstruction and interrogation of the CS genome-scale metabolic frameworks via Jupyter Notebook. The approach was tested using -omics data collected for Methylotuvimicrobium alcaliphilum 20ZR, a prominent microbial chassis for methane capturing and valorization. We optimized the previously reconstructed whole genome-scale metabolic network by adjusting the flux distribution using gene expression data. The outputs of the automatically reconstructed CS metabolic network were comparable to manually optimized iIA409 models for Ca-growth conditions. However, the CS model questions the reversibility of the phosphoketolase pathway and suggests higher flux via primary oxidation pathways. The model also highlighted unresolved carbon partitioning between assimilatory and catabolic pathways at the formaldehyde-formate node. Only a very few genes and only one enzyme with a predicted function in C1 metabolism, a homolog of the formaldehyde oxidation enzyme (fae1-2), showed a significant change in expression in La-growth conditions. The CS-GSM predictions agreed with the experimental measurements under the assumption that the Fae1-2 is a part of the tetrahydrofolate-linked pathway. The cellular roles of the tungsten (W)-dependent formate dehydrogenase (fdhAB) and fae homologs (fae1-2 and fae3) were investigated via mutagenesis. The phenotype of the fdhAB mutant followed the model prediction. Furthermore, a more significant reduction of the biomass yield was observed during growth in La-supplemented media, confirming a higher flux through formate. M. alcaliphilum 20ZR mutants lacking fae1-2 did not display any significant defects in methane or methanol-dependent growth. However, contrary to fae1, the fae1-2 homolog failed to restore the formaldehyde-activating enzyme function in complementation tests. Overall, the presented data suggest that the developed computational workflow supports the reconstruction and validation of CS-GSM networks of non-model microbes.

IMPORTANCE

The interrogation of various types of data is a routine strategy to explore the relationship between genotype and phenotype. An efficient approach for integrating and cross-comparing experimental multi-scale data in the context of whole-genome-based metabolic network reconstruction becomes a powerful tool that facilitates fundamental and applied research discoveries. The present study describes the reconstruction of a context-specific (CS) model for the methane-utilizing bacterium, Methylotuvimicrobium alcaliphilum 20ZR. M. alcaliphilum 20ZR is becoming an attractive microbial platform for the production of biofuels, chemicals, pharmaceuticals, and bio-sorbents for capturing atmospheric methane. We demonstrate that this pipeline can help reconstruct metabolic models that are similar to manually curated networks. Furthermore, the model is able to highlight previously overlooked pathways, thus advancing fundamental knowledge of non-model microbial systems or promoting their development toward biotechnological or environmental implementations.

Evolutionary genomics of the emergence of brown algae as key components of coastal ecosystems

Brown seaweeds are keystone species of coastal ecosystems, often forming extensive underwater forests, and are under considerable threat from climate change. In this study, analysis of multiple genomes has provided insights across the entire evolutionary history of this lineage, from initial emergence, through later diversification of the brown algal orders, down to microevolutionary events at the genus level. Emergence of the brown algal lineage was associated with a marked gain of new orthologous gene families, enhanced protein domain rearrangement, increased horizontal gene transfer events, and the acquisition of novel signaling molecules and key metabolic pathways, the latter notably related to biosynthesis of the alginate-based extracellular matrix, and halogen and phlorotannin biosynthesis. We show that brown algal genome diversification is tightly linked to phenotypic divergence, including changes in life cycle strategy and zoid flagellar structure. The study also showed that integration of large viral genomes has had a significant impact on brown algal genome content throughout the emergence of the lineage.

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.

pan-Draft: automated reconstruction of species-representative metabolic models from multiple genomes

The accurate reconstruction of genome-scale metabolic models (GEMs) for unculturable species poses challenges due to the incomplete and fragmented genetic information typical of metagenome-assembled genomes (MAGs). While existing tools leverage sequence homology from single genomes, this study introduces pan-Draft, a pan-reactome-based approach exploiting recurrent genetic evidence to determine the solid core structure of species-level GEMs. By comparing MAGs clustered at the species-level, pan-Draft addresses the issues due to the incompleteness and contamination of individual genomes, providing high-quality draft models and an accessory reactions catalog supporting the gapfilling step. This approach will improve our comprehension of metabolic functions of uncultured species.

From nets to networks: tools for deciphering phytoplankton metabolic interactions within communities and their global significance

Our oceans are populated with a wide diversity of planktonic organisms that form complex dynamic communities at the base of marine trophic networks. Within such communities are phytoplankton, unicellular photosynthetic taxa that provide an estimated half of global primary production and support biogeochemical cycles, along with other essential ecosystem services. One of the major challenges for microbial ecologists has been to try to make sense of this complexity. While phytoplankton distributions can be well explained by abiotic factors such as temperature and nutrient availability, there is increasing evidence that their ecological roles are tightly linked to their metabolic interactions with other plankton members through complex mechanisms (e.g. competition and symbiosis). Therefore, unravelling phytoplankton metabolic interactions is the key for inferring their dependency on, or antagonism with, other taxa and better integrating them into the context of carbon and nutrient fluxes in marine trophic networks. In this review, we attempt to summarize the current knowledge brought by ecophysiology, organismal imaging, in silico predictions and co-occurrence networks using 'omics data, highlighting successful combinations of approaches that may be helpful for future investigations of phytoplankton metabolic interactions within their complex communities.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

Machine learning for the advancement of genome-scale metabolic modeling

Constraint-based modeling (CBM) has evolved as the core systems biology tool to map the interrelations between genotype, phenotype, and external environment. The recent advancement of high-throughput experimental approaches and multi-omics strategies has generated a plethora of new and precise information from wide-ranging biological domains. On the other hand, the continuously growing field of machine learning (ML) and its specialized branch of deep learning (DL) provide essential computational architectures for decoding complex and heterogeneous biological data. In recent years, both multi-omics and ML have assisted in the escalation of CBM. Condition-specific omics data, such as transcriptomics and proteomics, helped contextualize the model prediction while analyzing a particular phenotypic signature. At the same time, the advanced ML tools have eased the model reconstruction and analysis to increase the accuracy and prediction power. However, the development of these multi-disciplinary methodological frameworks mainly occurs independently, which limits the concatenation of biological knowledge from different domains. Hence, we have reviewed the potential of integrating multi-disciplinary tools and strategies from various fields, such as synthetic biology, CBM, omics, and ML, to explore the biochemical phenomenon beyond the conventional biological dogma. How the integrative knowledge of these intersected domains has improved bioengineering and biomedical applications has also been highlighted. We categorically explained the conventional genome-scale metabolic model (GEM) reconstruction tools and their improvement strategies through ML paradigms. Further, the crucial role of ML and DL in omics data restructuring for GEM development has also been briefly discussed. Finally, the case-study-based assessment of the state-of-the-art method for improving biomedical and metabolic engineering strategies has been elaborated. Therefore, this review demonstrates how integrating experimental and in silico strategies can help map the ever-expanding knowledge of biological systems driven by condition-specific cellular information. This multiview approach will elevate the application of ML-based CBM in the biomedical and bioengineering fields for the betterment of society and the environment.

Genome-scale community modelling elucidates the metabolic interaction in Indian type-2 diabetic gut microbiota

Type-2 diabetes (T2D) is a rapidly growing multifactorial metabolic disorder that induces the onset of various diseases in the human body. The compositional and metabolic shift of the gut microbiota is a crucial factor behind T2D. Hence, gaining insight into the metabolic profile of the gut microbiota is essential for revealing their role in regulating the metabolism of T2D patients. Here, we have focused on the genome-scale community metabolic model reconstruction of crucial T2D-associated gut microbes. The model-based analysis of biochemical flux in T2D and healthy gut conditions showed distinct biochemical signatures and diverse metabolic interactions in the microbial community. The metabolic interactions encompass cross-feeding of short-chain fatty acids, amino acids, and vitamins among individual microbes within the community. In T2D conditions, a reduction in the metabolic flux of acetate, butyrate, vitamin B5, and bicarbonate was observed in the microbial community model, which can impact carbohydrate metabolism. The decline in butyrate levels is correlated with both insulin resistance and diminished glucose metabolism in T2D patients. Compared to the healthy gut, an overall reduction in glucose consumption and SCFA production flux was estimated in the T2D gut environment. Moreover, the decreased consumption profiles of branch chain amino acids (BCAAs) and aromatic amino acids (AAAs) in the T2D gut microbiota can be a distinct biomarker for T2D. Hence, the flux-level analysis of the microbial community model can provide insights into the metabolic reprogramming in diabetic gut microbiomes, which may be helpful in personalized therapeutics and diet design against T2D.

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.

Emerging methods for genome-scale metabolic modeling of microbial communities

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.

Using metabolic networks to predict cross-feeding and competition interactions between microorganisms

Understanding the interactions between microorganisms and their impact on bacterial behavior at the community level is a key research topic in microbiology. Different methods, relying on experimental or mathematical approaches based on the diverse properties of bacteria, are currently employed to study these interactions. Recently, the use of metabolic networks to understand the interactions between bacterial pairs has increased, highlighting the relevance of this approach in characterizing bacteria. In this study, we leverage the representation of bacteria through their metabolic networks to build a predictive model aimed at reducing the number of experimental assays required for designing bacterial consortia with specific behaviors. Our novel method for predicting cross-feeding or competition interactions between pairs of microorganisms utilizes metabolic network features. Machine learning classifiers are employed to determine the type of interaction from automatically reconstructed metabolic networks. Several algorithms were assessed and selected based on comprehensive testing and careful separation of manually compiled data sets obtained from literature sources. We used different classification algorithms, including K Nearest Neighbors, XGBoost, Support Vector Machine, and Random Forest, tested different parameter values, and implemented several data curation approaches to reduce the biological bias associated with our data set, ultimately achieving an accuracy of over 0.9. Our method holds substantial potential to advance the understanding of community behavior and contribute to the development of more effective approaches for consortia design.IMPORTANCEUnderstanding bacterial interactions at the community level is critical for microbiology, and leveraging metabolic networks presents an efficient and effective approach. The introduction of this novel method for predicting interactions through machine learning classifiers has the potential to advance the field by reducing the number of experimental assays required and contributing to the development of more effective bacterial consortia.

Genome-scale metabolic network models for industrial microorganisms metabolic engineering: Current advances and future prospects

The construction of high-performance microbial cell factories (MCFs) is the centerpiece of biomanufacturing. However, the complex metabolic regulatory network of microorganisms poses great challenges for the efficient design and construction of MCFs. The genome-scale metabolic network models (GSMs) can systematically simulate the metabolic regulation process of microorganisms in silico, providing effective guidance for the rapid design and construction of MCFs. In this review, we summarized the development status of 16 important industrial microbial GSMs, and further outline the technologies or methods that continuously promote high-quality GSMs construction from five aspects: I) Databases and modeling tools facilitate GSMs reconstruction; II) evolving gap-filling technologies; III) constraint-based model reconstruction; IV) advances in algorithms; and V) developed visualization tools. In addition, we also summarized the applications of GSMs in guiding metabolic engineering from four aspects: I) exploring and explaining metabolic features; II) predicting the effects of genetic perturbations on metabolism; III) predicting the optimal phenotype; IV) guiding cell factories construction in practical experiment. Finally, we discussed the development of GSMs, aiming to provide a reference for efficiently reconstructing GSMs and guiding metabolic engineering.

mergem: merging, comparing, and translating genome-scale metabolic models using universal identifiers

Numerous methods exist to produce and refine genome-scale metabolic models. However, due to the use of incompatible identifier systems for metabolites and reactions, computing and visualizing the metabolic differences and similarities of such models is a current challenge. Furthermore, there is a lack of automated tools that can combine the strengths of multiple reconstruction pipelines into a curated single comprehensive model by merging different drafts, which possibly use incompatible namespaces. Here we present mergem, a novel method to compare, merge, and translate two or more metabolic models. Using a universal metabolic identifier mapping system constructed from multiple metabolic databases, mergem robustly can compare models from different pipelines, merge their common elements, and translate their identifiers to other database systems. mergem is implemented as a command line tool, a Python package, and on the web-application Fluxer, which allows simulating and visually comparing multiple models with different interactive flux graphs. The ability to merge, compare, and translate diverse genome scale metabolic models can facilitate the curation of comprehensive reconstructions and the discovery of unique and common metabolic features among different organisms.

Development and applications of genome-scale metabolic network models

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

Strategies for tailoring functional microbial synthetic communities

Natural ecosystems harbor a huge reservoir of taxonomically diverse microbes that are important for plant growth and health. The vast diversity of soil microorganisms and their complex interactions make it challenging to pinpoint the main players important for the life support functions microbes can provide to plants, including enhanced tolerance to (a)biotic stress factors. Designing simplified microbial synthetic communities (SynComs) helps reduce this complexity to unravel the molecular and chemical basis and interplay of specific microbiome functions. While SynComs have been successfully employed to dissect microbial interactions or reproduce microbiome-associated phenotypes, the assembly and reconstitution of these communities have often been based on generic abundance patterns or taxonomic identities and co-occurrences but have only rarely been informed by functional traits. Here, we review recent studies on designing functional SynComs to reveal common principles and discuss multidimensional approaches for community design. We propose a strategy for tailoring the design of functional SynComs based on integration of high-throughput experimental assays with microbial strains and computational genomic analyses of their functional capabilities.

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.

Genome-Scale Model of Rhizopus microsporus: Metabolic integration of a fungal holobiont with its bacterial and viral endosymbionts

Rhizopus microsporus often lives in association with bacterial and viral symbionts that alter its biology. This fungal model represents an example of the complex interactions established among diverse organisms in functional holobionts. We constructed a Genome-Scale Model (GSM) of the fungal-bacterial-viral holobiont (iHol). We employed a constraint-based method to calculate the metabolic fluxes to decipher the metabolic interactions of the symbionts with their host. Our computational analyses of iHol simulate the holobiont's growth and the production of the toxin rhizoxin. Analyses of the calculated fluxes between R. microsporus in symbiotic (iHol) versus asymbiotic conditions suggest that changes in the lipid and nucleotide metabolism of the host are necessary for the functionality of the holobiont. Glycerol plays a pivotal role in the fungal-bacterial metabolic interaction, as its production does not compromise fungal growth, and Mycetohabitans bacteria can efficiently consume it. Narnavirus RmNV-20S and RmNV-23S affected the nucleotide metabolism without impacting the fungal-bacterial symbiosis. Our analyses highlighted the metabolic stability of Mycetohabitans throughout its co-evolution with the fungal host. We also predicted changes in reactions of the bacterial metabolism required for the active production of rhizoxin. This iHol is the first GSM of a fungal holobiont.

State-of the-Art Constraint-Based Modeling of Microbial Metabolism: From Basics to Context-Specific Models with a Focus on Methanotrophs

Methanotrophy is the ability of an organism to capture and utilize the greenhouse gas, methane, as a source of energy-rich carbon. Over the years, significant progress has been made in understanding of mechanisms for methane utilization, mostly in bacterial systems, including the key metabolic pathways, regulation and the impact of various factors (iron, copper, calcium, lanthanum, and tungsten) on cell growth and methane bioconversion. The implementation of -omics approaches provided vast amount of heterogeneous data that require the adaptation or development of computational tools for a system-wide interrogative analysis of methanotrophy. The genome-scale mathematical modeling of its metabolism has been envisioned as one of the most productive strategies for the integration of muti-scale data to better understand methane metabolism and enable its biotechnological implementation. Herein, we provide an overview of various computational strategies implemented for methanotrophic systems. We highlight functional capabilities as well as limitations of the most popular web resources for the reconstruction, modification and optimization of the genome-scale metabolic models for methane-utilizing bacteria.

Review: Towards the next-generation models of the rumen microbiome for enhancing predictive power and guiding sustainable production strategies

The rumen ecosystem harbours a galaxy of microbes working in syntrophy to carry out a metabolic cascade of hydrolytic and fermentative reactions. This fermentation process allows ruminants to harvest nutrients from a wide range of feedstuff otherwise inaccessible to the host. The interconnection between the ruminant and its rumen microbiota shapes key animal phenotypes such as feed efficiency and methane emissions and suggests the potential of reducing methane emissions and enhancing feed conversion into animal products by manipulating the rumen microbiota. Whilst significant technological progress in omics techniques has increased our knowledge of the rumen microbiota and its genome (microbiome), translating omics knowledge into effective microbial manipulation strategies remains a great challenge. This challenge can be addressed by modelling approaches integrating causality principles and thus going beyond current correlation-based approaches applied to analyse rumen microbial genomic data. However, existing rumen models are not yet adapted to capitalise on microbial genomic information. This gap between the rumen microbiota available omics data and the way microbial metabolism is represented in the existing rumen models needs to be filled to enhance rumen understanding and produce better predictive models with capabilities for guiding nutritional strategies. To fill this gap, the integration of computational biology tools and mathematical modelling frameworks is needed to translate the information of the metabolic potential of the rumen microbes (inferred from their genomes) into a mathematical object. In this paper, we aim to discuss the potential use of two modelling approaches for the integration of microbial genomic information into dynamic models. The first modelling approach explores the theory of state observers to integrate microbial time series data into rumen fermentation models. The second approach is based on the genome-scale network reconstructions of rumen microbes. For a given microorganism, the network reconstruction produces a stoichiometry matrix of the metabolism. This matrix is the core of the so-called genome-scale metabolic models which can be exploited by a plethora of methods comprised within the constraint-based reconstruction and analysis approaches. We will discuss how these methods can be used to produce the next-generation models of the rumen microbiome.

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.

Synthetic biology tools for environmental protection

Synthetic biology transforms the way we perceive biological systems. Emerging technologies in this field affect many disciplines of science and engineering. Traditionally, synthetic biology approaches were commonly aimed at developing cost-effective microbial cell factories to produce chemicals from renewable sources. Based on this, the immediate beneficial impact of synthetic biology on the environment came from reducing our oil dependency. However, synthetic biology is starting to play a more direct role in environmental protection. Toxic chemicals released by industries and agriculture endanger the environment, disrupting ecosystem balance and biodiversity loss. This review highlights synthetic biology approaches that can help environmental protection by providing remediation systems capable of sensing and responding to specific pollutants. Remediation strategies based on genetically engineered microbes and plants are discussed. Further, an overview of computational approaches that facilitate the design and application of synthetic biology tools in environmental protection is presented.

Unveiling abundance-dependent metabolic phenotypes of microbial communities

IMPORTANCE

In nature, organisms live in communities and not as isolated species, and their interactions provide a source of resilience to environmental disturbances. Despite their importance in ecology, human health, and industry, understanding how organisms interact in different environments remains an open question. In this work, we provide a novel approach that, only using genomic information, studies the metabolic phenotype exhibited by communities, where the exploration of suboptimal growth flux distributions and the composition of a community allows to unveil its capacity to respond to environmental changes, shedding light of the degrees of metabolic plasticity inherent to the community.

Reconciliation and evolution of Penicillium rubens genome-scale metabolic networks-What about specialised metabolism?

In recent years, genome sequencing of filamentous fungi has revealed a high proportion of specialised metabolites with growing pharmaceutical interest. However, detecting such metabolites through in silico genome analysis does not necessarily guarantee their expression under laboratory conditions. However, one plausible strategy for enabling their production lies in modifying the growth conditions. Devising a comprehensive experimental design testing in different culture environments is time-consuming and expensive. Therefore, using in silico modelling as a preliminary step, such as Genome-Scale Metabolic Network (GSMN), represents a promising approach to predicting and understanding the observed specialised metabolite production in a given organism. To address these questions, we reconstructed a new high-quality GSMN for the Penicillium rubens Wisconsin 54-1255 strain, a commonly used model organism. Our reconstruction, iPrub22, adheres to current convention standards and quality criteria, incorporating updated functional annotations, orthology searches with different GSMN templates, data from previous reconstructions, and manual curation steps targeting primary and specialised metabolites. With a MEMOTE score of 74% and a metabolic coverage of 45%, iPrub22 includes 5,192 unique metabolites interconnected by 5,919 reactions, of which 5,033 are supported by at least one genomic sequence. Of the metabolites present in iPrub22, 13% are categorised as belonging to specialised metabolism. While our high-quality GSMN provides a valuable resource for investigating known phenotypes expressed in P. rubens, our analysis identifies bottlenecks related, in particular, to the definition of what is a specialised metabolite, which requires consensus within the scientific community. It also points out the necessity of accessible, standardised and exhaustive databases of specialised metabolites. These questions must be addressed to fully unlock the potential of natural product production in P. rubens and other filamentous fungi. Our work represents a foundational step towards the objective of rationalising the production of natural products through GSMN modelling.

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.

Inferring and comparing metabolism across heterogeneous sets of annotated genomes using AuCoMe

Comparative analysis of genome-scale metabolic networks (GSMNs) may yield important information on the biology, evolution, and adaptation of species. However, it is impeded by the high heterogeneity of the quality and completeness of structural and functional genome annotations, which may bias the results of such comparisons. To address this issue, we developed AuCoMe, a pipeline to automatically reconstruct homogeneous GSMNs from a heterogeneous set of annotated genomes without discarding available manual annotations. We tested AuCoMe with three data sets, one bacterial, one fungal, and one algal, and showed that it successfully reduces technical biases while capturing the metabolic specificities of each organism. Our results also point out shared and divergent metabolic traits among evolutionarily distant algae, underlining the potential of AuCoMe to accelerate the broad exploration of metabolic evolution across the tree of life.

Reconstructor: a COBRApy compatible tool for automated genome-scale metabolic network reconstruction with parsimonious flux-based gap-filling

MOTIVATION

Genome-scale metabolic network reconstructions (GENREs) are valuable for understanding cellular metabolism in silico. Several tools exist for automatic GENRE generation. However, these tools frequently (i) do not readily integrate with some of the widely-used suites of packaged methods available for network analysis, (ii) lack effective network curation tools, (iii) are not sufficiently user-friendly, and (iv) often produce low-quality draft reconstructions.

RESULTS

Here, we present Reconstructor, a user-friendly, COBRApy-compatible tool that produces high-quality draft reconstructions with reaction and metabolite naming conventions that are consistent with the ModelSEED biochemistry database and includes a gap-filling technique based on the principles of parsimony. Reconstructor can generate SBML GENREs from three input types: annotated protein .fasta sequences (Type 1 input), a BLASTp output (Type 2), or an existing SBML GENRE that can be further gap-filled (Type 3). While Reconstructor can be used to create GENREs of any species, we demonstrate the utility of Reconstructor with bacterial reconstructions. We demonstrate how Reconstructor readily generates high-quality GENRES that capture strain, species, and higher taxonomic differences in functional metabolism of bacteria and are useful for further biological discovery.

AVAILABILITY AND IMPLEMENTATION

The Reconstructor Python package is freely available for download. Complete installation and usage instructions and benchmarking data are available at http://github.com/emmamglass/reconstructor.

Insights into the potential for mutualistic and harmful host-microbe interactions affecting brown alga freshwater acclimation

Microbes can modify their hosts' stress tolerance, thus potentially enhancing their ecological range. An example of such interactions is Ectocarpus subulatus, one of the few freshwater-tolerant brown algae. This tolerance is partially due to its (un)cultivated microbiome. We investigated this phenomenon by modifying the microbiome of laboratory-grown E. subulatus using mild antibiotic treatments, which affected its ability to grow in low salinity. Low salinity acclimation of these algal-bacterial associations was then compared. Salinity significantly impacted bacterial and viral gene expression, albeit in different ways across algal-bacterial communities. In contrast, gene expression of the host and metabolite profiles were affected almost exclusively in the freshwater-intolerant algal-bacterial communities. We found no evidence of bacterial protein production that would directly improve algal stress tolerance. However, vitamin K synthesis is one possible bacterial service missing specifically in freshwater-intolerant cultures in low salinity. In this condition, we also observed a relative increase in bacterial transcriptomic activity and the induction of microbial genes involved in the biosynthesis of the autoinducer AI-1, a quorum-sensing regulator. This could have resulted in dysbiosis by causing a shift in bacterial behaviour in the intolerant algal-bacterial community. Together, these results provide two promising hypotheses to be examined by future targeted experiments. Although they apply only to the specific study system, they offer an example of how bacteria may impact their host's stress response.

Genome-Scale Modeling and Systems Metabolic Engineering of Vibrio natriegens for the Production of 1,3-Propanediol

The fastest-growing bacterium Vibrio natriegens is a highly promising next-generation workhorse for molecular biology and industrial biotechnology. In this work, we described the workflows for developing genome-scale metabolic models and genome-editing protocols for engineering Vibrio natriegens. A case study for metabolic engineering of Vibrio natriegens for the production of 1,3-propanediol was also presented.

Multi-genome metabolic modeling predicts functional inter-dependencies in the Arabidopsis root microbiome

BACKGROUND

From a theoretical ecology point of view, microbiomes are far more complex than expected. Besides competition and competitive exclusion, cooperative microbe-microbe interactions have to be carefully considered. Metabolic dependencies among microbes likely explain co-existence in microbiota.

METHODOLOGY

In this in silico study, we explored genome-scale metabolic models (GEMs) of 193 bacteria isolated from Arabidopsis thaliana roots. We analyzed their predicted producible metabolites under simulated nutritional constraints including "root exudate-mimicking growth media" and assessed the potential of putative metabolic exchanges of by- and end-products to avoid those constraints.

RESULTS

We found that the genome-encoded metabolic potential is quantitatively and qualitatively clustered by phylogeny, highlighting metabolic differentiation between taxonomic groups. Random, synthetic combinations of increasing numbers of strains (SynComs) indicated that the number of producible compounds by GEMs increased with average phylogenetic distance, but that most SynComs were centered around an optimal phylogenetic distance. Moreover, relatively small SynComs could reflect the capacity of the whole community due to metabolic redundancy. Inspection of 30 specific end-product metabolites (i.e., target metabolites: amino acids, vitamins, phytohormones) indicated that the majority of the strains had the genetic potential to produce almost all the targeted compounds. Their production was predicted (1) to depend on external nutritional constraints and (2) to be facilitated by nutritional constraints mimicking root exudates, suggesting nutrient availability and root exudates play a key role in determining the number of producible metabolites. An answer set programming solver enabled the identification of numerous combinations of strains predicted to depend on each other to produce these targeted compounds under severe nutritional constraints thus indicating a putative sub-community level of functional redundancy.

CONCLUSIONS

This study predicts metabolic restrictions caused by available nutrients in the environment. By extension, it highlights the importance of the environment for niche potential, realization, partitioning, and overlap. Our results also suggest that metabolic dependencies and cooperation among root microbiota members compensate for environmental constraints and help maintain co-existence in complex microbial communities. Video Abstract.

ChiMera: an easy to use pipeline for bacterial genome based metabolic network reconstruction, evaluation and visualization

BACKGROUND

Genome-scale metabolic reconstruction tools have been developed in the last decades. They have helped to reconstruct eukaryotic and prokaryotic metabolic models, which have contributed to fields, e.g., genetic engineering, drug discovery, prediction of phenotypes, and other model-driven discoveries. However, the use of these programs requires a high level of bioinformatic skills. Moreover, the functionalities required to build models are scattered throughout multiple tools, requiring knowledge and experience for utilizing several tools.

RESULTS

Here we present ChiMera, which combines tools used for model reconstruction, prediction, and visualization. ChiMera uses CarveMe in the reconstruction module, generating a gap-filled draft reconstruction able to produce growth predictions using flux balance analysis for gram-positive and gram-negative bacteria. ChiMera also contains two modules for metabolic network visualization. The first module generates maps for the most important pathways, e.g., glycolysis, nucleotides and amino acids biosynthesis, fatty acid oxidation and biosynthesis and core-metabolism. The second module produces a genome-wide metabolic map, which can be used to retrieve KEGG pathway information for each compound in the model. A module to investigate gene essentiality and knockout is also present.

CONCLUSIONS

Overall, ChiMera uses automation algorithms to combine a variety of tools to automatically perform model creation, gap-filling, flux balance analysis (FBA), and metabolic network visualization. ChiMera models readily provide metabolic insights that can aid genetic engineering projects, prediction of phenotypes, and model-driven discoveries.

Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology

Computational models are often employed in systems biology to study the dynamic behaviours of complex systems. With the rise in the number of computational models, finding ways to improve the reusability of these models and their ability to reproduce virtual experiments becomes critical. Correct and effective model annotation in community-supported and standardised formats is necessary for this improvement. Here, we present recent efforts toward a common framework for annotated, accessible, reproducible and interoperable computational models in biology, and discuss key challenges of the field.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

merlin, an improved framework for the reconstruction of high-quality genome-scale metabolic models

Genome-scale metabolic models have been recognised as useful tools for better understanding living organisms' metabolism. merlin (https://www.merlin-sysbio.org/) is an open-source and user-friendly resource that hastens the models' reconstruction process, conjugating manual and automatic procedures, while leveraging the user's expertise with a curation-oriented graphical interface. An updated and redesigned version of merlin is herein presented. Since 2015, several features have been implemented in merlin, along with deep changes in the software architecture, operational flow, and graphical interface. The current version (4.0) includes the implementation of novel algorithms and third-party tools for genome functional annotation, draft assembly, model refinement, and curation. Such updates increased the user base, resulting in multiple published works, including genome metabolic (re-)annotations and model reconstructions of multiple (lower and higher) eukaryotes and prokaryotes. merlin version 4.0 is the only tool able to perform template based and de novo draft reconstructions, while achieving competitive performance compared to state-of-the art tools both for well and less-studied organisms.

COMMIT: Consideration of metabolite leakage and community composition improves microbial community reconstructions

Composition and functions of microbial communities affect important traits in diverse hosts, from crops to humans. Yet, mechanistic understanding of how metabolism of individual microbes is affected by the community composition and metabolite leakage is lacking. Here, we first show that the consensus of automatically generated metabolic reconstructions improves the quality of the draft reconstructions, measured by comparison to reference models. We then devise an approach for gap filling, termed COMMIT, that considers metabolites for secretion based on their permeability and the composition of the community. By applying COMMIT with two soil communities from the Arabidopsis thaliana culture collection, we could significantly reduce the gap-filling solution in comparison to filling gaps in individual reconstructions without affecting the genomic support. Inspection of the metabolic interactions in the soil communities allows us to identify microbes with community roles of helpers and beneficiaries. Therefore, COMMIT offers a versatile fully automated solution for large-scale modelling of microbial communities for diverse biotechnological applications.

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

Modeling approaches for probing cross-feeding interactions in the human gut microbiome

Microbial communities perform emergent activities that are essentially different from those carried by their individual members. The gut microbiome and its metabolites have a significant impact on the host, contributing to homeostasis or disease. Food molecules shape this community, being fermented through cross-feeding interactions of metabolites such as lactate, acetate, and amino acids, or products derived from macromolecule degradation. Mathematical and experimental approaches have been applied to understand and predict the interactions between microorganisms in complex communities such as the gut microbiota. Rational and mechanistic understanding of microbial interactions is essential to exploit their metabolic activities and identify keystone taxa and metabolites. The latter could be used in turn to modulate or replicate the metabolic behavior of the community in different contexts. This review aims to highlight recent experimental and modeling approaches for studying cross-feeding interactions within the gut microbiome. We focus on short-chain fatty acid production and fiber fermentation, which are fundamental processes in human health and disease. Special attention is paid to modeling approaches, particularly kinetic and genome-scale stoichiometric models of metabolism, to integrate experimental data under different diet and health conditions. Finally, we discuss limitations and challenges for the broad application of these modeling approaches and their experimental verification for improving our understanding of the mechanisms of microbial interactions.

Advances in Genome-Scale Metabolic Modeling toward Microbial Community Analysis of the Human Microbiome

A genome-scale metabolic model (GEM) represents metabolic pathways of an organism in a mathematical form and can be built using biochemistry and genome annotation data. GEMs are invaluable for understanding organisms since they analyze the metabolic capabilities and behaviors quantitatively and can predict phenotypes. The development of high-throughput data collection techniques led to an immense increase in omics data such as metagenomics, which expand our knowledge on the human microbiome, but this also created a need for systematic analysis of these data. In recent years, GEMs have also been reconstructed for microbial species, including human gut microbiota, and methods for the analysis of microbial communities have been developed to examine the interaction between the organisms or the host. The purpose of this review is to provide a comprehensive guide for the applications of GEMs in microbial community analysis. Starting with GEM repositories, automatic GEM reconstruction tools, and quality control of models, this review will give insights into microbe-microbe and microbe-host interaction predictions and optimization of microbial community models. Recent studies that utilize microbial GEMs and personalized models to infer the influence of microbiota on human diseases such as inflammatory bowel diseases (IBD) or Parkinson's disease are exemplified. Being powerful system biology tools for both species-level and community-level analysis of microbes, GEMs integrated with omics data and machine learning techniques will be indispensable for studying the microbiome and their effects on human physiology as well as for deciphering the mechanisms behind human diseases.

Genome Scale Modeling to Study the Metabolic Competition between Cells in the Tumor Microenvironment

The tumor's physiology emerges from the dynamic interplay of numerous cell types, such as cancer cells, immune cells and stromal cells, within the tumor microenvironment. Immune and cancer cells compete for nutrients within the tumor microenvironment, leading to a metabolic battle between these cell populations. Tumor cells can reprogram their metabolism to meet the high demand of building blocks and ATP for proliferation, and to gain an advantage over the action of immune cells. The study of the metabolic reprogramming mechanisms underlying cancer requires the quantification of metabolic fluxes which can be estimated at the genome-scale with constraint-based or kinetic modeling. Constraint-based models use a set of linear constraints to simulate steady-state metabolic fluxes, whereas kinetic models can simulate both the transient behavior and steady-state values of cellular fluxes and concentrations. The integration of cell- or tissue-specific data enables the construction of context-specific models that reflect cell-type- or tissue-specific metabolic properties. While the available modeling frameworks enable limited modeling of the metabolic crosstalk between tumor and immune cells in the tumor stroma, future developments will likely involve new hybrid kinetic/stoichiometric formulations.

Comparative Transcriptome Analysis Provides Insights into Response of Ulva compressa to Fluctuating Salinity Conditions

Ulva compressa, a green tide-forming species, can adapt to hypo-salinity conditions, such as estuaries and brackish lakes. To understand the underlying molecular mechanisms of hypo-salinity stress tolerance, transcriptome-wide gene expression profiles in U. compressa were created using digital gene expression profiles. The RNA-seq data were analyzed based on the comparison of differently expressed genes involved in specific pathways under hypo-salinity and recovery conditions. The up-regulation of genes in photosynthesis and glycolysis pathways may contribute to the recovery of photosynthesis and energy metabolism, which could provide sufficient energy for the tolerance under long-term hyposaline stress. Multiple strategies, such as ion transportation and osmolytes metabolism, were performed to maintain the osmotic homeostasis. Additionally, several long noncoding RNA were differently expressed during the stress, which could play important roles in the osmotolerance. Our work will serve as an essential foundation for the understanding of the tolerance mechanism of U. compressa under the fluctuating salinity conditions.

Modelling microbial communities: Harnessing consortia for biotechnological applications

Microbes propagate and thrive in complex communities, and there are many benefits to studying and engineering microbial communities instead of single strains. Microbial communities are being increasingly leveraged in biotechnological applications, as they present significant advantages such as the division of labour and improved substrate utilisation. Nevertheless, they also present some interesting challenges to surmount for the design of efficient biotechnological processes. In this review, we discuss key principles of microbial interactions, followed by a deep dive into genome-scale metabolic models, focussing on a vast repertoire of constraint-based modelling methods that enable us to characterise and understand the metabolic capabilities of microbial communities. Complementary approaches to model microbial communities, such as those based on graph theory, are also briefly discussed. Taken together, these methods provide rich insights into the interactions between microbes and how they influence microbial community productivity. We finally overview approaches that allow us to generate and test numerous synthetic community compositions, followed by tools and methodologies that can predict effective genetic interventions to further improve the productivity of communities. With impending advancements in high-throughput omics of microbial communities, the stage is set for the rapid expansion of microbial community engineering, with a significant impact on biotechnological processes.

Metabolic control of nitrogen fixation in rhizobium-legume symbioses

Rhizobia induce nodule formation on legume roots and differentiate into bacteroids, which catabolize plant-derived dicarboxylates to reduce atmospheric N₂ into ammonia. Despite the agricultural importance of this symbiosis, the mechanisms that govern carbon and nitrogen allocation in bacteroids and promote ammonia secretion to the plant are largely unknown. Using a metabolic model derived from genome-scale datasets, we show that carbon polymer synthesis and alanine secretion by bacteroids facilitate redox balance in microaerobic nodules. Catabolism of dicarboxylates induces not only a higher oxygen demand but also a higher NADH/NAD⁺ ratio than sugars. Modeling and ¹³C metabolic flux analysis indicate that oxygen limitation restricts the decarboxylating arm of the tricarboxylic acid cycle, which limits ammonia assimilation into glutamate. By tightly controlling oxygen supply and providing dicarboxylates as the energy and electron source donors for N₂ fixation, legumes promote ammonia secretion by bacteroids. This is a defining feature of rhizobium-legume symbioses.

Elucidation of trophic interactions in an unusual single-cell nitrogen-fixing symbiosis using metabolic modeling

Marine nitrogen-fixing microorganisms are an important source of fixed nitrogen in oceanic ecosystems. The colonial cyanobacterium Trichodesmium and diatom symbionts were thought to be the primary contributors to oceanic N₂ fixation until the discovery of the unusual uncultivated symbiotic cyanobacterium UCYN-A (Candidatus Atelocyanobacterium thalassa). UCYN-A has atypical metabolic characteristics lacking the oxygen-evolving photosystem II, the tricarboxylic acid cycle, the carbon-fixation enzyme RuBisCo and de novo biosynthetic pathways for a number of amino acids and nucleotides. Therefore, it is obligately symbiotic with its single-celled haptophyte algal host. UCYN-A receives fixed carbon from its host and returns fixed nitrogen, but further insights into this symbiosis are precluded by both UCYN-A and its host being uncultured. In order to investigate how this syntrophy is coordinated, we reconstructed bottom-up genome-scale metabolic models of UCYN-A and its algal partner to explore possible trophic scenarios, focusing on nitrogen fixation and biomass synthesis. Since both partners are uncultivated and only the genome sequence of UCYN-A is available, we used the phylogenetically related Chrysochromulina tobin as a proxy for the host. Through the use of flux balance analysis (FBA), we determined the minimal set of metabolites and biochemical functions that must be shared between the two organisms to ensure viability and growth. We quantitatively investigated the metabolic characteristics that facilitate daytime N₂ fixation in UCYN-A and possible oxygen-scavenging mechanisms needed to create an anaerobic environment to allow nitrogenase to function. This is the first application of an FBA framework to examine the tight metabolic coupling between uncultivated microbes in marine symbiotic communities and provides a roadmap for future efforts focusing on such specialized systems.

Reduction of dinitrogen gas to biologically useful forms via nitrogen fixation is a key component of the biogeochemical cycle. In the marine environment, the cyanobacteria UCYN-A (Candidatus Atelocyanobacterium thalassa) has been found to be a primary contributor to biological nitrogen fixation at a global scale. UCYN-A exhibits a highly streamlined genome which lacks genes coding for essential cyanobacterial processes such as the energy-generating TCA cycle, oxygen-producing photosystem II, the carbon-fixing RuBisCo and de novo production pathways for numerous amino acids and nucleotides. Thus, it exists in a symbiosis with unicellular planktonic algae where it exchanges fixed nitrogen for fixed carbon with its host. However, both UCYN-A and its symbiotic partner remain uncultured under laboratory conditions. This necessitates implementing a computational approach to glean insights into UCYN-A’s unique physiology and metabolic processes governing the symbiotic association. To this end, we develop a constraints-based framework that infers all possible trophic scenarios consistent with the observed data. Possible mechanisms employed by UCYN-A to accommodate diazotrophy with daytime carbon fixation by the host (i.e., two mutually incompatible processes) are also elucidated. We envision that the developed framework using UCYN-A and its algal host will be used as a roadmap and motivate the study of similarly unique microbial systems in the future.

Robustness analysis of metabolic predictions in algal microbial communities based on different annotation pipelines

Animals, plants, and algae rely on symbiotic microorganisms for their development and functioning. Genome sequencing and genomic analyses of these microorganisms provide opportunities to construct metabolic networks and to analyze the metabolism of the symbiotic communities they constitute. Genome-scale metabolic network reconstructions rest on information gained from genome annotation. As there are multiple annotation pipelines available, the question arises to what extent differences in annotation pipelines impact outcomes of these analyses. Here, we compare five commonly used pipelines (Prokka, MaGe, IMG, DFAST, RAST) from predicted annotation features (coding sequences, Enzyme Commission numbers, hypothetical proteins) to the metabolic network-based analysis of symbiotic communities (biochemical reactions, producible compounds, and selection of minimal complementary bacterial communities). While Prokka and IMG produced the most extensive networks, RAST and DFAST networks produced the fewest false positives and the most connected networks with the fewest dead-end metabolites. Our results underline differences between the outputs of the tested pipelines at all examined levels, with small differences in the draft metabolic networks resulting in the selection of different microbial consortia to expand the metabolic capabilities of the algal host. However, the consortia generated yielded similar predicted producible compounds and could therefore be considered functionally interchangeable. This contrast between selected communities and community functions depending on the annotation pipeline needs to be taken into consideration when interpreting the results of metabolic complementarity analyses. In the future, experimental validation of bioinformatic predictions will likely be crucial to both evaluate and refine the pipelines and needs to be coupled with increased efforts to expand and improve annotations in reference databases.

gapseq: informed prediction of bacterial metabolic pathways and reconstruction of accurate metabolic models

Genome-scale metabolic models of microorganisms are powerful frameworks to predict phenotypes from an organism's genotype. While manual reconstructions are laborious, automated reconstructions often fail to recapitulate known metabolic processes. Here we present gapseq ( https://github.com/jotech/gapseq ), a new tool to predict metabolic pathways and automatically reconstruct microbial metabolic models using a curated reaction database and a novel gap-filling algorithm. On the basis of scientific literature and experimental data for 14,931 bacterial phenotypes, we demonstrate that gapseq outperforms state-of-the-art tools in predicting enzyme activity, carbon source utilisation, fermentation products, and metabolic interactions within microbial communities.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

Semi-Quantitative Targeted Gas Chromatography-Mass Spectrometry Profiling Supports a Late Side-Chain Reductase Cycloartenol-to-Cholesterol Biosynthesis Pathway in Brown Algae

Sterols are biologically important molecules that serve as membrane fluidity regulators and precursors of signaling molecules, either endogenous or involved in biotic interactions. There is currently no model of their biosynthesis pathways in brown algae. Here, we benefit from the availability of genome data and gas chromatography-mass spectrometry (GC-MS) sterol profiling using a database of internal standards to build such a model. We expand the set of identified sterols in 11 species of red, brown, and green macroalgae and integrate these new data with genomic data. Our analyses suggest that some metabolic reactions may be conserved despite the loss of canonical eukaryotic enzymes, like the sterol side-chain reductase (SSR). Our findings are consistent with the principle of metabolic pathway drift through enzymatic replacement and show that cholesterol synthesis from cycloartenol may be a widespread but variable pathway among chlorophyllian eukaryotes. Among the factors contributing to this variability, one could be the recruitment of cholesterol biosynthetic intermediates to make signaling molecules, such as the mozukulins. These compounds were found in some brown algae belonging to Ectocarpales, and we here provide a first mozukulin biosynthetic model. Our results demonstrate that integrative approaches can already be used to infer experimentally testable models, which will be useful to further investigate the biological roles of those newly identified algal pathways.

Genome-Scale Metabolic Modeling of Escherichia coli and Its Chassis Design for Synthetic Biology Applications

Genome-scale metabolic modeling is and will continue to play a central role in computational systems metabolic engineering and synthetic biology applications for the productions of chemicals and antibiotics. To that end, a survey and workflows of methods used for the development of high-quality genome-scale metabolic models (GEMs) and chassis design for synthetic biology are described here. The chapter consists of two parts (a) the methods of development of GEMs (Escherichia coli as a case study) and (b) E. coli chassis design for synthetic production of 1,4-butanediol (BDO). The methods described here can guide existing and future development of GEMs coupled with host chassis design for synthetic productions of novel antibiotics.

Metage2Metabo, microbiota-scale metabolic complementarity for the identification of key species

To capture the functional diversity of microbiota, one must identify metabolic functions and species of interest within hundreds or thousands of microorganisms. We present Metage2Metabo (M2M) a resource that meets the need for de novo functional screening of genome-scale metabolic networks (GSMNs) at the scale of a metagenome, and the identification of critical species with respect to metabolic cooperation. M2M comprises a flexible pipeline for the characterisation of individual metabolisms and collective metabolic complementarity. In addition, M2M identifies key species, that are meaningful members of the community for functions of interest. We demonstrate that M2M is applicable to collections of genomes as well as metagenome-assembled genomes, permits an efficient GSMN reconstruction with Pathway Tools, and assesses the cooperation potential between species. M2M identifies key organisms by reducing the complexity of a large-scale microbiota into minimal communities with equivalent properties, suitable for further analyses.