Genome-scale metabolic models: reconstruction and analysis

Arduino Soft Sensor for Monitoring Schizochytrium sp. Fermentation, a Proof of Concept for the Industrial Application of Genome-Scale Metabolic Models in the Context of Pharma 4.0

Schizochytrium sp. is a microorganism cultured for producing docosahexaenoic acid (DHA). Genome-scale metabolic modeling (GEM) is a promising technique for describing gen-protein-reactions in cells, but with still limited industrial application due to its complexity and high computation requirements. In this work, we simplified GEM results regarding the relationship between the specific oxygen uptake rate (−rO2), the specific growth rate (µ), and the rate of lipid synthesis (rL) using an evolutionary algorithm for developing a model that can be used by a soft sensor for fermentation monitoring. The soft sensor estimated the concentration of active biomass (X), glutamate (N), lipids (L), and DHA in a Schizochytrium sp. fermentation using the dissolved oxygen tension (DO) and the oxygen mass transfer coefficient (kLa) as online input variables. The soft sensor model described the biomass concentration response of four reported experiments characterized by different kLa values. The average range normalized root-mean-square error for X, N, L, and DHA were equal to 1.1, 1.3, 1.1, and 3.2%, respectively, suggesting an acceptable generalization capacity. The feasibility of implementing the soft sensor over a low-cost electronic board was successfully tested using an Arduino UNO, showing a novel path for applying GEM-based soft sensors in the context of Pharma 4.0.

Genome-scale modeling of Chinese hamster ovary cells by hybrid semi-parametric flux balance analysis

Flux balance analysis (FBA) is currently the standard method to compute metabolic fluxes in genome-scale networks. Several FBA extensions employing diverse objective functions and/or constraints have been published. Here we propose a hybrid semi-parametric FBA extension that combines mechanistic-level constraints (parametric) with empirical constraints (non-parametric) in the same linear program. A CHO dataset with 27 measured exchange fluxes obtained from 21 reactor experiments served to evaluate the method. The mechanistic constraints were deduced from a reduced CHO-K1 genome-scale network with 686 metabolites, 788 reactions and 210 degrees of freedom. The non-parametric constraints were obtained by principal component analysis of the flux dataset. The two types of constraints were integrated in the same linear program showing comparable computational cost to standard FBA. The hybrid FBA is shown to significantly improve the specific growth rate prediction under different constraints scenarios. A metabolically efficient cell growth feed targeting minimal byproducts accumulation was designed by hybrid FBA. It is concluded that integrating parametric and nonparametric constraints in the same linear program may be an efficient approach to reduce the solution space and to improve the predictive power of FBA methods when critical mechanistic information is missing.

Personalized Immunotherapy in Colorectal Cancers: Where Do We Stand?

Colorectal cancer (CRC) is the second leading cause of cancer death in the world. Immunotherapy using monoclonal antibodies, immune-checkpoint inhibitors, adoptive cell therapy, and cancer vaccines has raised great hopes for treating poor prognosis metastatic CRCs that are resistant to the conventional therapies. However, high inter-tumor and intra-tumor heterogeneity hinder the success of immunotherapy in CRC. Patients with a similar tumor phenotype respond differently to the same immunotherapy regimen. Mutation-based classification, molecular subtyping, and immunoscoring of CRCs facilitated the multi-aspect grouping of CRC patients and improved immunotherapy. Personalized immunotherapy using tumor-specific neoantigens provides the opportunity to consider each patient as an independent group deserving of individualized immunotherapy. In the recent decade, the development of sequencing and multi-omics techniques has helped us classify patients more precisely. The expansion of such advanced techniques along with the neoantigen-based immunotherapy could herald a new era in treating heterogeneous tumors such as CRC. In this review article, we provided the latest findings in immunotherapy of CRC. We elaborated on the heterogeneity of CRC patients as a bottleneck of CRC immunotherapy and reviewed the latest advances in personalized immunotherapy to overcome CRC heterogeneity.

Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.

Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism-targeting as a promising antiviral strategy.

Unraveling the Genome of a High Yielding Colombian Sugarcane Hybrid

Recent developments in High Throughput Sequencing (HTS) technologies and bioinformatics, including improved read lengths and genome assemblers allow the reconstruction of complex genomes with unprecedented quality and contiguity. Sugarcane has one of the most complicated genomes among grassess with a haploid length of 1Gbp and a ploidies between 8 and 12. In this work, we present a genome assembly of the Colombian sugarcane hybrid CC 01-1940. Three types of sequencing technologies were combined for this assembly: PacBio long reads, Illumina paired short reads, and Hi-C reads. We achieved a median contig length of 34.94 Mbp and a total genome assembly of 903.2 Mbp. We annotated a total of 63,724 protein coding genes and performed a reconstruction and comparative analysis of the sucrose metabolism pathway. Nucleotide evolution measurements between orthologs with close species suggest that divergence between Saccharum officinarum and Saccharum spontaneum occurred <2 million years ago. Synteny analysis between CC 01-1940 and the S. spontaneum genome confirms the presence of translocation events between the species and a random contribution throughout the entire genome in current sugarcane hybrids. Analysis of RNA-Seq data from leaf and root tissue of contrasting sugarcane genotypes subjected to water stress treatments revealed 17,490 differentially expressed genes, from which 3,633 correspond to genes expressed exclusively in tolerant genotypes. We expect the resources presented here to serve as a source of information to improve the selection processes of new varieties of the breeding programs of sugarcane.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Biochemical Characteristics and a Genome-Scale Metabolic Model of an Indian Euryhaline Cyanobacterium with High Polyglucan Content

Marine cyanobacteria are promising microbes to capture and convert atmospheric CO₂ and light into biomass and valuable industrial bio-products. Yet, reports on metabolic characteristics of non-model cyanobacteria are scarce. In this report, we show that an Indian euryhaline Synechococcus sp. BDU 130192 has biomass accumulation comparable to a model marine cyanobacterium and contains approximately double the amount of total carbohydrates, but significantly lower protein levels compared to Synechococcus sp. PCC 7002 cells. Based on its annotated chromosomal genome sequence, we present a genome scale metabolic model (GSMM) of this cyanobacterium, which we have named as iSyn706. The model includes 706 genes, 908 reactions, and 900 metabolites. The difference in the flux balance analysis (FBA) predicted flux distributions between Synechococcus sp. PCC 7002 and Synechococcus sp. BDU130192 strains mimicked the differences in their biomass compositions. Model-predicted oxygen evolution rate for Synechococcus sp. BDU130192 was found to be close to the experimentally-measured value. The model was analyzed to determine the potential of the strain for the production of various industrially-useful products without affecting growth significantly. This model will be helpful to researchers interested in understanding the metabolism as well as to design metabolic engineering strategies for the production of industrially-relevant compounds.

Robust Analytical Methods for the Accurate Quantification of the Total Biomass Composition of Mammalian Cells

Biomass composition is an important input for genome-scale metabolic models and has a big impact on their predictive capabilities. However, researchers often rely on generic data for biomass composition, e.g. collected from similar organisms. This leads to inaccurate predictions, because biomass composition varies between different cell lines, conditions, and growth phases. In this chapter we present protocols for the determination of the biomass composition of Chinese Hamster Ovary (CHO) cells. These methods can easily be adapted to other types of mammalian cells. The protocols include the quantification of cell dry mass and of the main biomass components, namely protein, lipid, DNA, RNA, and carbohydrates. Cell dry mass is determined gravimetrically by weighing a defined number of cells. Amino acid composition and protein content are measured by gas chromatography mass spectrometry. Lipids are quantified by shotgun mass spectrometry, which provides quantities for the different lipid classes and also the distribution of fatty acids. RNA is purified and then quantified spectrophotometrically. The methods for DNA and carbohydrates are simple fluorometric and colorimetric assays adapted to a 96-well plate format. To ensure quantitative results, internal standards or spike-in controls are used in all methods, e.g. to account for possible matrix effects or loss of material. Finally, the last section provides a guide on how to convert the measured data into biomass equations, which can then be integrated into a metabolic model.

Bioinoculant capability enhancement through metabolomics and systems biology approaches

Bioinoculants are eco-friendly microorganisms, and their products are utilized for improving the potential of soil and fulfill the nutrients requirement for the host plant. The agricultural yield has increased due to the use of bioinoculants over chemical-based fertilizers, and thus it generates interest in understanding the innovation process by various methods. By gene-editing tool, the desired gene product can be changed for engineered microbial inoculants. We have also described various modern biotechnological tools like constraint-based modeling, OptKnock, flux balance analysis and modeling of the biological network for enhancing the bioinoculant capability. These fluxes give the fascinating perception of the metabolic network in the absence of comprehensive kinetic information. These tools also help in the stimulation of the metabolic networks by incorporation of enzyme-encoding genes. The present review explains the use of systems biology and gene-editing tools for improving the capability of bioinoculants. Moreover, this review also emphasizes on the challenges and future perspective of systems biology and its multidisciplinary facets.

Integration of transcriptomic data in a genome-scale metabolic model to investigate the link between obesity and breast cancer

BACKGROUND

Obesity is a complex disorder associated with an increased risk of developing several comorbid chronic diseases, including postmenopausal breast cancer. Although many studies have investigated this issue, the link between body weight and either risk or poor outcome of breast cancer is still to characterize. Systems biology approaches, based on the integration of multiscale models and data from a wide variety of sources, are particularly suitable for investigating the underlying molecular mechanisms of complex diseases. In this scenario, GEnome-scale metabolic Models (GEMs) are a valuable tool, since they represent the metabolic structure of cells and provide a functional scaffold for simulating and quantifying metabolic fluxes in living organisms through constraint-based mathematical methods. The integration of omics data into the structural information described by GEMs allows to build more accurate descriptions of metabolic states.

RESULTS

In this work, we exploited gene expression data of postmenopausal breast cancer obese and lean patients to simulate a curated GEM of the human adipocyte, available in the Human Metabolic Atlas database. To this aim, we used a published algorithm which exploits a data-driven approach to overcome the limitation of defining a single objective function to simulate the model. The flux solutions were used to build condition-specific graphs to visualise and investigate the reaction networks and their properties. In particular, we performed a network topology differential analysis to search for pattern differences and identify the principal reactions associated with significant changes across the two conditions under study.

CONCLUSIONS

Metabolic network models represent an important source to study the metabolic phenotype of an organism in different conditions. Here we demonstrate the importance of exploiting Next Generation Sequencing data to perform condition-specific GEM analyses. In particular, we show that the qualitative and quantitative assessment of metabolic fluxes modulated by gene expression data provides a valuable method for investigating the mechanisms associated with the phenotype under study, and can foster our interpretation of biological phenomena.

The intertwined metabolism during symbiotic nitrogen fixation elucidated by metabolic modelling

Genome-scale metabolic network models can be used for various analyses including the prediction of metabolic responses to changes in the environment. Legumes are well known for their rhizobial symbiosis that introduces nitrogen into the global nutrient cycle. Here, we describe a fully compartmentalised, mass and charge-balanced, genome-scale model of the clover Medicago truncatula, which has been adopted as a model organism for legumes. We employed flux balance analysis to demonstrate that the network is capable of producing biomass components in experimentally observed proportions, during day and night. By connecting the plant model to a model of its rhizobial symbiont, Sinorhizobium meliloti, we were able to investigate the effects of the symbiosis on metabolic fluxes and plant growth and could demonstrate how oxygen availability influences metabolic exchanges between plant and symbiont, thus elucidating potential benefits of inter organism amino acid cycling. We thus provide a modelling framework, in which the interlinked metabolism of plants and nodules can be studied from a theoretical perspective.

Large-scale computational models of liver metabolism: How far from the clinics?

Understanding the dynamics of human liver metabolism is fundamental for effective diagnosis and treatment of liver diseases. This knowledge can be obtained with systems biology/medicine approaches that account for the complexity of hepatic responses and their systemic consequences in other organs. Computational modeling can reveal hidden principles of the system by classification of individual components, analyzing their interactions and simulating the effects that are difficult to investigate experimentally. Herein, we review the state-of-the-art computational models that describe liver dynamics from metabolic, gene regulatory, and signal transduction perspectives. We focus especially on large-scale liver models described either by genome scale metabolic networks or an object-oriented approach. We also discuss the benefits and limitations of each modeling approach and their value for clinical applications in diagnosis, therapy, and prevention of liver diseases as well as precision medicine in hepatology. (Hepatology 2017;66:1323-1334).

Towards personalized medicine of colorectal cancer

Efforts in colorectal cancer (CRC) research aim to improve early detection and treatment for metastatic stages which could translate into better prognosis of this disease. One of the major challenges that hinder these efforts is the heterogeneous nature of CRC and involvement of diverse molecular pathways. New large-scale 'omics' technologies are making it possible to generate, analyze and interpret biological data from molecular determinants of CRC. The developments of sophisticated computational analyses would allow information from different omics platforms to be integrated, thus providing new insights into the biology of CRC. Together, these technological advances and an improved mechanistic understanding might allow CRC to be clinically managed at the level of the individual patient. This review provides an account of the current challenges in CRC management and an insight into how new technologies could allow the development of personalized medicine for CRC.

Computational modelling of genome-scale metabolic networks and its application to CHO cell cultures

Genome-scale metabolic models (GEMs) have become increasingly important in recent years. Currently, GEMs are the most accurate in silico representation of the genotype-phenotype link. They allow us to study complex networks from the systems perspective. Their application may drastically reduce the amount of experimental and clinical work, improve diagnostic tools and increase our understanding of complex biological phenomena. GEMs have also demonstrated high potential for the optimisation of bio-based production of recombinant proteins. Herein, we review the basic concepts, methods, resources and software tools used for the reconstruction and application of GEMs. We overview the evolution of the modelling efforts devoted to the metabolism of Chinese Hamster Ovary (CHO) cells. We present a case study on CHO cell metabolism under different amino acid depletions. This leads us to the identification of the most influential as well as essential amino acids in selected CHO cell lines.

BioSIMP: Using Software Testing Techniques for Sampling and Inference in Biological Organisms

Years of research in software engineering has given us novel ways to reason about, test, and predict the behavior of complex software systems that contain hundreds of thousands of lines of code. Many of these techniques have been inspired by nature such as genetic algorithms, swarm intelligence, and ant colony optimization. In this paper we reverse the direction and present BioSIMP, a process that models and predicts the behavior of biological organisms to aid in the emerging field of systems biology. It utilizes techniques from testing and modeling of highly-configurable software systems. Using both experimental and simulation data we show that BioSIMP can find important environmental factors in two microbial organisms. However, we learn that in order to fully reason about the complexity of biological systems, we will need to extend existing or create new software engineering techniques.

Genome-Scale Reconstruction of the Human Astrocyte Metabolic Network

Astrocytes are the most abundant cells of the central nervous system; they have a predominant role in maintaining brain metabolism. In this sense, abnormal metabolic states have been found in different neuropathological diseases. Determination of metabolic states of astrocytes is difficult to model using current experimental approaches given the high number of reactions and metabolites present. Thus, genome-scale metabolic networks derived from transcriptomic data can be used as a framework to elucidate how astrocytes modulate human brain metabolic states during normal conditions and in neurodegenerative diseases. We performed a Genome-Scale Reconstruction of the Human Astrocyte Metabolic Network with the purpose of elucidating a significant portion of the metabolic map of the astrocyte. This is the first global high-quality, manually curated metabolic reconstruction network of a human astrocyte. It includes 5,007 metabolites and 5,659 reactions distributed among 8 cell compartments, (extracellular, cytoplasm, mitochondria, endoplasmic reticle, Golgi apparatus, lysosome, peroxisome and nucleus). Using the reconstructed network, the metabolic capabilities of human astrocytes were calculated and compared both in normal and ischemic conditions. We identified reactions activated in these two states, which can be useful for understanding the astrocytic pathways that are affected during brain disease. Additionally, we also showed that the obtained flux distributions in the model, are in accordance with literature-based findings. Up to date, this is the most complete representation of the human astrocyte in terms of inclusion of genes, proteins, reactions and metabolic pathways, being a useful guide for in-silico analysis of several metabolic behaviors of the astrocyte during normal and pathologic states.

Reconstruction of the Fatty Acid Biosynthetic Pathway of Exiguobacterium antarcticum B7 Based on Genomic and Bibliomic Data

Exiguobacterium antarcticum B7 is extremophile Gram-positive bacteria able to survive in cold environments. A key factor to understanding cold adaptation processes is related to the modification of fatty acids composing the cell membranes of psychrotrophic bacteria. In our study we show the in silico reconstruction of the fatty acid biosynthesis pathway of E. antarcticum B7. To build the stoichiometric model, a semiautomatic procedure was applied, which integrates genome information using KEGG and RAST/SEED. Constraint-based methods, namely, Flux Balance Analysis (FBA) and elementary modes (EM), were applied. FBA was implemented in the sense of hexadecenoic acid production maximization. To evaluate the influence of the gene expression in the fluxome analysis, FBA was also calculated using the log₂⁡FC values obtained in the transcriptome analysis at 0°C and 37°C. The fatty acid biosynthesis pathway showed a total of 13 elementary flux modes, four of which showed routes for the production of hexadecenoic acid. The reconstructed pathway demonstrated the capacity of E. antarcticum B7 to de novo produce fatty acid molecules. Under the influence of the transcriptome, the fluxome was altered, promoting the production of short-chain fatty acids. The calculated models contribute to better understanding of the bacterial adaptation at cold environments.

Cyanobacterial Biofuels: Strategies and Developments on Network and Modeling

Cyanobacteria, the phototrophic microorganisms, have attracted much attention recently as a promising source for environmentally sustainable biofuels production. However, barriers for commercial markets of cyanobacteria-based biofuels concern the economic feasibility. Miscellaneous strategies for improving the production performance of cyanobacteria have thus been developed. Among these, the simple ad hoc strategies resulting in failure to optimize fully cell growth coupled with desired product yield are explored. With the advancement of genomics and systems biology, a new paradigm toward systems metabolic engineering has been recognized. In particular, a genome-scale metabolic network reconstruction and modeling is a crucial systems-based tool for whole-cell-wide investigation and prediction. In this review, the cyanobacterial genome-scale metabolic models, which offer a system-level understanding of cyanobacterial metabolism, are described. The main process of metabolic network reconstruction and modeling of cyanobacteria are summarized. Strategies and developments on genome-scale network and modeling through the systems metabolic engineering approach are advanced and employed for efficient cyanobacterial-based biofuels production.

Systems Biology of Microbial Exopolysaccharides Production

Exopolysaccharides (EPSs) produced by diverse group of microbial systems are rapidly emerging as new and industrially important biomaterials. Due to their unique and complex chemical structures and many interesting physicochemical and rheological properties with novel functionality, the microbial EPSs find wide range of commercial applications in various fields of the economy such as food, feed, packaging, chemical, textile, cosmetics and pharmaceutical industry, agriculture, and medicine. EPSs are mainly associated with high-value applications, and they have received considerable research attention over recent decades with their biocompatibility, biodegradability, and both environmental and human compatibility. However, only a few microbial EPSs have achieved to be used commercially due to their high production costs. The emerging need to overcome economic hurdles and the increasing significance of microbial EPSs in industrial and medical biotechnology call for the elucidation of the interrelations between metabolic pathways and EPS biosynthesis mechanism in order to control and hence enhance its microbial productivity. Moreover, a better understanding of biosynthesis mechanism is a significant issue for improvement of product quality and properties and also for the design of novel strains. Therefore, a systems-based approach constitutes an important step toward understanding the interplay between metabolism and EPS biosynthesis and further enhances its metabolic performance for industrial application. In this review, primarily the microbial EPSs, their biosynthesis mechanism, and important factors for their production will be discussed. After this brief introduction, recent literature on the application of omics technologies and systems biology tools for the improvement of production yields will be critically evaluated. Special focus will be given to EPSs with high market value such as xanthan, levan, pullulan, and dextran.

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

BACKGROUND

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

RESULTS

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG's annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii's metabolism in comparison to the one of S. cerevisiae (post-WGD - whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

CONCLUSIONS

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii's metabolism.

Merging multiple omics datasets in silico: statistical analyses and data interpretation

By the combinations of high-throughput analytical technologies in the fields of transcriptomics, proteomics, and metabolomics, we are now able to gain comprehensive and quantitative snapshots of the intracellular processes. Dynamic intracellular activities and their regulations can be elucidated by systematic observation of these multi-omics data. On the other hand, careful statistical analysis is necessary for such integration, since each of the omics layers as well as the specific analytical methodologies harbor different levels of noise and variations. Moreover, interpretation of such multitude of data requires an intuitive pathway context. Here we describe such statistical methods for the integration and comparison of multi-omics data, as well as the computational methods for pathway reconstruction, ID conversion, mapping, and visualization that play key roles for the efficient study of multi-omics information.

Reconciliation of metabolites and biochemical reactions for metabolic networks

Genome-scale metabolic network reconstructions are now routinely used in the study of metabolic pathways, their evolution and design. The development of such reconstructions involves the integration of information on reactions and metabolites from the scientific literature as well as public databases and existing genome-scale metabolic models. The reconciliation of discrepancies between data from these sources generally requires significant manual curation, which constitutes a major obstacle in efforts to develop and apply genome-scale metabolic network reconstructions. In this work, we discuss some of the major difficulties encountered in the mapping and reconciliation of metabolic resources and review three recent initiatives that aim to accelerate this process, namely BKM-react, MetRxn and MNXref (presented in this article). Each of these resources provides a pre-compiled reconciliation of many of the most commonly used metabolic resources. By reducing the time required for manual curation of metabolite and reaction discrepancies, these resources aim to accelerate the development and application of high-quality genome-scale metabolic network reconstructions and models.