Genome-scale Metabolic Network Research Articles

Insight into the biotechnology potential of Alicyclobacillus tolerans from whole genome sequence analysis and genome-scale metabolic network modeling.

Extremophilic bacteria have numerous uncovered biotechnological potentials. Acidophilic bacteria are important iron oxidizers that are valuable in bioleaching and in studying extreme environments on earth and in space. Despite their obvious potential, little is known about the genetic traits that underpin their metabolic functions, which are equally poorly understood from a mechanistic perspective. Novel bioinformatics and computational biology pipelines can be used to analyze whole genomes to obtain insights into the phenotypic potential of organisms as well as develop a mathematical model representation of metabolism. Whole-genome sequence analysis and a genome-scale metabolic network model was curated for an iron-oxidizing bacterium initially isolated from an acid mine drainage in Turkey, previously identified as Alicyclobacillus tolerans. The genome contained a high proportion of genes for energy generation from carbohydrates, amino acids synthesis and conversion, nucleic acid metabolism and repair which contribute to robust adaption to their extreme environments. Several candidate genes for pyrite metabolism, iron uptake, regulation and storage, as well as genes for resistance to important heavy metals were annotated. A curated genome-scale metabolic network analysis accurately predicted facultative anaerobic growth, heterotrophic characteristics, and growth on a wide variety of carbon sources. This is the first in-depth in silico analysis of A. tolerans to the best of our knowledge which is expected to lay the groundwork for future research and drive innovations in environmental microbiology and biotechnological applications. The genomic data and mechanistic framework will have applications in biomining, synthetic geomicrobiology on earth, as well as for space exploration and settlement.

Pinpointing Cancer Sub-Type Specific Metabolic Tasks Facilitates Identification of Anti-cancer Targets.

Metabolic reprogramming is one of the hallmarks of tumorigenesis. Understanding the metabolic changes in cancer cells may provide attractive therapeutic targets and new strategies for cancer therapy. The metabolic states are not the same in different cancer types or subtypes, even within the same sample of solid tumors. In order to understand the heterogeneity of cancer cells, we used the Pareto tasks inference method to analyze the metabolic tasks of different cancers, including breast cancer, lung cancer, digestive organ cancer, digestive tract cancer, and reproductive cancer. We found that cancer subtypes haves different propensities toward metabolic tasks, and the biological significance of these metabolic tasks also varies greatly. Normal cells treat metabolic tasks uniformly, while different cancer cells focus on different pathways. We then integrated the metabolic tasks into the multi-objective genome-scale metabolic network model, which shows higher accuracy in the in silico prediction of cell states after gene knockout than the conventional biomass maximization model. The predicted potential single drug targets could potentially turn into biomarkers or drug design targets. We further implemented the multi-objective genome-scale metabolic network model to predict synthetic lethal target pairs of the Basal and Luminal B subtypes of breast cancer. By analyzing the predicted synthetic lethal targets, we found that mitochondrial enzymes are potential targets for drug combinations. Our study quantitatively analyzes the metabolic tasks of cancer and establishes cancer type-specific metabolic models, which opens a new window for the development of specific anti-cancer drugs and provides promising treatment plans for specific cancer subtypes.

Reconstruction of a generic genome-scale metabolic network for chicken: Investigating network connectivity and finding potential biomarkers.

Chicken is the first sequenced avian that has a crucial role in human life for its meat and egg production. Because of various metabolic disorders, study the metabolism of chicken cell is important. Herein, the first genome-scale metabolic model of a chicken cell named iES1300, consists of 2427 reactions, 2569 metabolites, and 1300 genes, was reconstructed manually based on KEGG, BiGG, CHEBI, UNIPROT, REACTOME, and MetaNetX databases. Interactions of metabolic genes for growth were examined for E. coli, S. cerevisiae, human, and chicken metabolic models. The results indicated robustness to genetic manipulation for iES1300 similar to the results for human. iES1300 was integrated with transcriptomics data using algorithms and Principal Component Analysis was applied to compare context-specific models of the normal, tumor, lean and fat cell lines. It was found that the normal model has notable metabolic flexibility in the utilization of various metabolic pathways, especially in metabolic pathways of the carbohydrate metabolism, compared to the others. It was also concluded that the fat and tumor models have similar growth metabolisms and the lean chicken model has a more active lipid and carbohydrate metabolism.

Development of metabolic models with multiple constraints: a review

Constraint-based genome-scale metabolic network models (genome-scale metabolic models, GEMs) have been widely used to predict metabolic phenotypes. In addition to stoichiometric constraints, other constraints such as enzyme availability and thermodynamic feasibility may also limit the cellular phenotype solution space. Recently, extended GEM models considering either enzymatic or thermodynamic constraints have been developed to improve model prediction accuracy. This review summarizes the recent progresses on metabolic models with multiple constraints (MCGEMs). We presented the construction methods and various applications of MCGEMs including the simulation of gene knockout, prediction of biologically feasible pathways and identification of bottleneck steps. By integrating multiple constraints in a consistent modeling framework, MCGEMs can predict the metabolic bottlenecks and key controlling and modification targets for pathway optimization more precisely, and thus may provide more reliable design results to guide metabolic engineering of industrially important microorganisms.

CobraMod: a pathway-centric curation tool for constraint-based metabolic models.

SummaryCOnstraint-Based Reconstruction and Analysis of genome-scale metabolic models has become a widely used tool to understand metabolic network behavior at a large scale. However, existing reconstruction tools lack functionalities to address modellers' common objective to study metabolic networks on the pathway level. Thus, we developed CobraMod—a Python package for pathway-centric modification and extension of genome-scale metabolic networks. CobraMod can integrate data from various metabolic pathway databases as well as user-curated information. Our tool tests newly added metabolites, reactions and pathways against multiple curation criteria, suggests manual curation steps and provides the user with records of changes to ensure high quality metabolic reconstructions. CobraMod uses the visualization tool Escher for pathway representation and offers simple customization options for comparison of pathways and flux distributions. Our package enables coherent and reproducible workflows as it can be seamlessly integrated with COBRApy and Escher.Availability and implementationThe source code can be found at https://github.com/Toepfer-Lab/cobramod/ and can be installed with pip. The documentation including tutorials is available at https://cobramod.readthedocs.io/.

Quantifying cumulative phenotypic and genomic evidence for procedural generation of metabolic network reconstructions.

Genome-scale metabolic network reconstructions (GENREs) are valuable tools for understanding microbial metabolism. The process of automatically generating GENREs includes identifying metabolic reactions supported by sufficient genomic evidence to generate a draft metabolic network. The draft GENRE is then gapfilled with additional reactions in order to recapitulate specific growth phenotypes as indicated with associated experimental data. Previous methods have implemented absolute mapping thresholds for the reactions automatically included in draft GENREs; however, there is growing evidence that integrating annotation evidence in a continuous form can improve model accuracy. There is a need for flexibility in the structure of GENREs to better account for uncertainty in biological data, unknown regulatory mechanisms, and context-specificity associated with data inputs. To address this issue, we present a novel method that provides a framework for quantifying combined genomic, biochemical, and phenotypic evidence for each biochemical reaction during automated GENRE construction. Our method, Constraint-based Analysis Yielding reaction Usage across metabolic Networks (CANYUNs), generates accurate GENREs with a quantitative metric for the cumulative evidence for each reaction included in the network. The structuring of CANYUNs allows for the simultaneous integration of three data inputs while maintaining all supporting evidence for biochemical reactions that may be active in an organism. CANYUNs is designed to maximize the utility of experimental and annotation datasets and to ultimately assist in the curation of the reference datasets used for the automatic construction of metabolic networks. We validated CANYUNs by generating an E. coli K-12 model and compared it to the manually curated reconstruction iML1515. Finally, we demonstrated the use of CANYUNs to build a model by generating an E. coli Nissle CANYUNs model using novel phenotypic data that we collected. This method may address key challenges for the procedural construction of metabolic networks by leveraging uncertainty and redundancy in biological data.

Genome-Scale Modeling Specifies the Metabolic Capabilities of Rhizophagus irregularis.

ABSTRACTRhizophagus irregularis is one of the most extensively studied arbuscular mycorrhizal fungi (AMF) that forms symbioses with and improves the performance of many crops. Lack of transformation protocol for R. irregularis renders it challenging to investigate molecular mechanisms that shape the physiology and interactions of this AMF with plants. Here, we used all published genomics, transcriptomics, and metabolomics resources to gain insights into the metabolic functionalities of R. irregularis by reconstructing its high-quality genome-scale metabolic network that considers enzyme constraints. Extensive validation tests with the enzyme-constrained metabolic model demonstrated that it can be used to (i) accurately predict increased growth of R. irregularis on myristate with minimal medium; (ii) integrate enzyme abundances and carbon source concentrations that yield growth predictions with high and significant Spearman correlation ( = 0.74) to measured hyphal dry weight; and (iii) simulate growth rate increases with tighter association of this AMF with the host plant across three fungal structures. Based on the validated model and system-level analyses that integrate data from transcriptomics studies, we predicted that differences in flux distributions between intraradical mycelium and arbuscles are linked to changes in amino acid and cofactor biosynthesis. Therefore, our results demonstrated that the enzyme-constrained metabolic model can be employed to pinpoint mechanisms driving developmental and physiological responses of R. irregularis to different environmental cues. In conclusion, this model can serve as a template for other AMF and paves the way to identify metabolic engineering strategies to modulate fungal metabolic traits that directly affect plant performance.IMPORTANCE Mounting evidence points to the benefits of the symbiotic interactions between the arbuscular mycorrhiza fungus Rhizophagus irregularis and crops; however, the molecular mechanisms underlying the physiological responses of this fungus to different host plants and environments remain largely unknown. We present a manually curated, enzyme-constrained, genome-scale metabolic model of R. irregularis that can accurately predict experimentally observed phenotypes. We show that this high-quality model provides an entry point into better understanding the metabolic and physiological responses of this fungus to changing environments due to the availability of different nutrients. The model can be used to design metabolic engineering strategies to tailor R. irregularis metabolism toward improving the performance of host plants.

Nonlinear Model Predictive Control based on multi-scale models: is it worth the complexity?

Nonlinear model predictive control (NMPC) has been used to control a variety of bioprocesses. As NMPC requires full state information, a nonlinear moving horizon estimator (NMHE) is used to reconstruct the states from the available measurements. In context of bioprocesses, both NMPC and NMHE typically use macroscopic mass balance type models. These models ignore all intracellular information that is available via, e.g., Genome-Scale Metabolic networks. However, the use of these models adds significantly to the computational complexity of the controller. In this paper, we aim to assess whether incorporating metabolic networks in an NMPC framework leads to improvement in the tracking performance of a continuous bioreactor. To avoid complex bilevel optimisation problems, our NMPC scheme utilises the dynamic Metabolic Flux Analysis framework. In this framework, the intracellular fluxes are estimated from the measurements of extracellular metabolites. This is achieved by an NMHE scheme. The tracking performance of this metabolic network-based NMPC-NMHE is compared to a simple macroscopic model based NMPC in which all the states, i.e., the extracellular metabolites are assumed to be measured. In our case study, the tracking performance of the metabolic network-based NMPC is marginally better than the performance of the macroscopic model based NMPC. The added advantage of metabolic network based NMPC comes from the insight they offer from the estimated fluxes. Moreover, the performance and the insights obtained depend on the complexity of the metabolic network used.

Metabolic signatures of regulation by phosphorylation and acetylation

SummaryAcetylation and phosphorylation are highly conserved posttranslational modifications (PTMs) that regulate cellular metabolism, yet how metabolic control is shared between these PTMs is unknown. Here we analyze transcriptome, proteome, acetylome, and phosphoproteome datasets in E. coli, S. cerevisiae, and mammalian cells across diverse conditions using CAROM, a new approach that uses genome-scale metabolic networks and machine learning to classify targets of PTMs. We built a single machine learning model that predicted targets of each PTM in a condition across all three organisms based on reaction attributes (AUC>0.8). Our model predicted phosphorylated enzymes during a mammalian cell-cycle, which we validate using phosphoproteomics. Interpreting the machine learning model using game theory uncovered enzyme properties including network connectivity, essentiality, and condition-specific factors such as maximum flux that differentiate targets of phosphorylation from acetylation. The conserved and predictable partitioning of metabolic regulation identified here between these PTMs may enable rational rewiring of regulatory circuits.

Model-based dynamic engineering of Escherichia coli for N-acetylglucosamine overproduction

N-acetylglucosamine (GlcNAc), a glucosamine derivative, has a wide range of applications in pharmaceutical fields, and there is an increasing interest in the efficient production of GlcNAc genetic engineered bacteria. In this work, Escherichia coli ATCC 25947 (DE3) strain was engineered by a model-based dynamic regulation strategy achieving GlcNAc overproduction. First, the GlcNAc synthetic pathway was introduced into E. coli, and through flux balance analysis of the genome-scale metabolic network model, metabolic engineering strategies were generated to further increase GlcNAc yield. Knock-out of genes poxB and ldhA, encoding pyruvate oxidase and lactate dehydrogenase, increased GlcNAc titer by 5.1%. Furthermore, knocking out N-acetylmuramic acid 6-phosphate etherase encoded by murQ and enhancing glutamine synthetase encoded by glnA gene further increased GlcNAc titer to 130.8 g/L. Analysis of metabolic flux balance showed that GlcNAc production maximization requires the strict dynamic restriction of the reactions catalyzed by pfkA and zwf to balance cell growth and product synthesis. Hence, a dynamic regulatory system was constructed by combining the CRISPRi (clustered regularly interspaced short palindromic repeats interference) system with the lactose operon lacI and the transcription factor pdhR, allowing the cell to respond to the concentration of pyruvate and IPTG to dynamically repress pfkA and zwf transcription. Finally, the engineered bacteria with the dynamic regulatory system produced 143.8 g/L GlcNAc in a 30-L bioreactor in 55 h with a yield reaching 0.539 g/g glucose. Taken together, this work significantly enhanced the GlcNAc production of E. coli. Moreover, it provides a systematic, effective, and universal way to improve the synthetic ability of other engineered strains.

Relative flux trade-offs and optimization of metabolic network functionalities

Trade-offs between traits are present across different levels of biological systems and ultimately reflect constraints imposed by physicochemical laws and the structure of underlying biochemical networks. Yet, mechanistic explanation of how trade-offs between molecular traits arise and how they relate to optimization of fitness-related traits remains elusive. Here, we introduce the concept of relative flux trade-offs and propose a constraint-based approach, termed FluTOr, to identify metabolic reactions whose fluxes are in relative trade-off with respect to an optimized fitness-related cellular task, like growth. We then employed FluTOr to identify relative flux trade-offs in the genome-scale metabolic networks of Escherichia coli, Saccharomyces cerevisiae, and Arabidopsis thaliana. For the metabolic models of E. coli and S. cerevisiae we showed that: (i) the identified relative flux trade-offs depend on the carbon source used and that (ii) reactions that participated in relative trade-offs in both species were implicated in cofactor biosynthesis. In contrast to the two microorganisms, the relative flux trade-offs for the metabolic model of A. thaliana did not depend on the available nitrogen sources, reflecting the differences in the underlying metabolic network as well as the considered environments. Lastly, the established connection between relative flux trade-offs allowed us to identify overexpression targets that can be used to optimize fitness-related traits. Altogether, our computational approach and findings demonstrate how relative flux trade-offs can shape optimization of metabolic tasks, important in biotechnological applications.

Identification of flux trade-offs in metabolic networks

Trade-offs are inherent to biochemical networks governing diverse cellular functions, from gene expression to metabolism. Yet, trade-offs between fluxes of biochemical reactions in a metabolic network have not been formally studied. Here, we introduce the concept of absolute flux trade-offs and devise a constraint-based approach, termed FluTO, to identify and enumerate flux trade-offs in a given genome-scale metabolic network. By employing the metabolic networks of Escherichia coli and Saccharomyces cerevisiae, we demonstrate that the flux trade-offs are specific to carbon sources provided but that reactions involved in the cofactor and prosthetic group biosynthesis are present in trade-offs across all carbon sources supporting growth. We also show that absolute flux trade-offs depend on the biomass reaction used to model the growth of Arabidopsis thaliana under different carbon and nitrogen conditions. The identified flux trade-offs reflect the tight coupling between nitrogen, carbon, and sulphur metabolisms in leaves of C3 plants. Altogether, FluTO provides the means to explore the space of alternative metabolic routes reflecting the constraints imposed by inherent flux trade-offs in large-scale metabolic networks.

Global connectivity in genome-scale metabolic networks revealed by comprehensive FBA-based pathway analysis

BackgroundGraph-based analysis (GBA) of genome-scale metabolic networks has revealed system-level structures such as the bow-tie connectivity that describes the overall mass flow in a network. However, many pathways obtained by GBA are biologically impossible, making it difficult to study how the global structures affect the biological functions of a network. New method that can calculate the biologically relevant pathways is desirable for structural analysis of metabolic networks.ResultsHere, we present a new method to determine the bow-tie connectivity structure by calculating possible pathways between any pairs of metabolites in the metabolic network using a flux balance analysis (FBA) approach to ensure that the obtained pathways are biologically relevant. We tested this method with 15 selected high-quality genome-scale metabolic models from BiGG database. The results confirmed the key roles of central metabolites in network connectivity, locating in the core part of the bow-tie structure, the giant strongly connected component (GSC). However, the sizes of GSCs revealed by GBA are significantly larger than those by FBA approach. A great number of metabolites in the GSC from GBA actually cannot be produced from or converted to other metabolites through a mass balanced pathway and thus should not be in GSC but in other subsets of the bow-tie structure. In contrast, the bow-tie structural classification of metabolites obtained by FBA is more biologically relevant and suitable for the study of the structure-function relationships of genome scale metabolic networks.ConclusionsThe FBA based pathway calculation improve the biologically relevant classification of metabolites in the bow-tie connectivity structure of the metabolic network, taking us one step further toward understanding how such system-level structures impact the biological functions of an organism.

An updated genome-scale metabolic network reconstruction of Pseudomonas aeruginosa PA14 to characterize mucin-driven shifts in bacterial metabolism

Mucins are present in mucosal membranes throughout the body and play a key role in the microbe clearance and infection prevention. Understanding the metabolic responses of pathogens to mucins will further enable the development of protective approaches against infections. We update the genome-scale metabolic network reconstruction (GENRE) of one such pathogen, Pseudomonas aeruginosa PA14, through metabolic coverage expansion, format update, extensive annotation addition, and literature-based curation to produce iPau21. We then validate iPau21 through MEMOTE, growth rate, carbon source utilization, and gene essentiality testing to demonstrate its improved quality and predictive capabilities. We then integrate the GENRE with transcriptomic data in order to generate context-specific models of P. aeruginosa metabolism. The contextualized models recapitulated known phenotypes of unaltered growth and a differential utilization of fumarate metabolism, while also revealing an increased utilization of propionate metabolism upon MUC5B exposure. This work serves to validate iPau21 and demonstrate its utility for providing biological insights.

TRIMER: Transcription Regulation Integrated with Metabolic Regulation

SummaryThere has been extensive research in predictive modeling of genome-scale metabolic reaction networks. Living systems involve complex stochastic processes arising from interactions among different biomolecules. For more accurate and robust prediction of target metabolic behavior under different conditions, not only metabolic reactions but also the genetic regulatory relationships involving transcription factors (TFs) affecting these metabolic reactions should be modeled. We have developed a modeling and simulation pipeline enabling the analysis of Transcription Regulation Integrated with Metabolic Regulation: TRIMER. TRIMER utilizes a Bayesian network (BN) inferred from transcriptomes to model the transcription factor regulatory network. TRIMER then infers the probabilities of the gene states relevant to the metabolism of interest, and predicts the metabolic fluxes and their changes that result from the deletion of one or more transcription factors at the genome scale. We demonstrate TRIMER’s applicability to both simulated and experimental data and provide performance comparison with other existing approaches.

A Systematic Strategy to Find Potential Therapeutic Targets for Pseudomonas aeruginosa Using Integrated Computational Models.

Pseudomonas aeruginosa is an opportunistic human pathogen that has been a constant global health problem due to its ability to cause infection at different body sites and its resistance to a broad spectrum of clinically available antibiotics. The World Health Organization classified multidrug-resistant Pseudomonas aeruginosa among the top-ranked organisms that require urgent research and development of effective therapeutic options. Several approaches have been taken to achieve these goals, but they all depend on discovering potential drug targets. The large amount of data obtained from sequencing technologies has been used to create computational models of organisms, which provide a powerful tool for better understanding their biological behavior. In the present work, we applied a method to integrate transcriptome data with genome-scale metabolic networks of Pseudomonas aeruginosa. We submitted both metabolic and integrated models to dynamic simulations and compared their performance with published in vitro growth curves. In addition, we used these models to identify potential therapeutic targets and compared the results to analyze the assumption that computational models enriched with biological measurements can provide more selective and (or) specific predictions. Our results demonstrate that dynamic simulations from integrated models result in more accurate growth curves and flux distribution more coherent with biological observations. Moreover, identifying drug targets from integrated models is more selective as the predicted genes were a subset of those found in the metabolic models. Our analysis resulted in the identification of 26 non-host homologous targets. Among them, we highlighted five top-ranked genes based on lesser conservation with the human microbiome. Overall, some of the genes identified in this work have already been proposed by different approaches and (or) are already investigated as targets to antimicrobial compounds, reinforcing the benefit of using integrated models as a starting point to selecting biologically relevant therapeutic targets.

An extended reconstruction of human gut microbiota metabolism of dietary compounds

Understanding how diet and gut microbiota interact in the context of human health is a key question in personalized nutrition. Genome-scale metabolic networks and constraint-based modeling approaches are promising to systematically address this complex problem. However, when applied to nutritional questions, a major issue in existing reconstructions is the limited information about compounds in the diet that are metabolized by the gut microbiota. Here, we present AGREDA, an extended reconstruction of diet metabolism in the human gut microbiota. AGREDA adds the degradation pathways of 209 compounds present in the human diet, mainly phenolic compounds, a family of metabolites highly relevant for human health and nutrition. We show that AGREDA outperforms existing reconstructions in predicting diet-specific output metabolites from the gut microbiota. Using 16S rRNA gene sequencing data of faecal samples from Spanish children representing different clinical conditions, we illustrate the potential of AGREDA to establish relevant metabolic interactions between diet and gut microbiota.

Transcriptomic and Metabolic Network Analysis of Metabolic Reprogramming and IGF-1 Modulation in SCA3 Transgenic Mice.

Spinocerebellar ataxia type 3 (SCA3) is a genetic neurodegenerative disease for which a cure is still needed. Growth hormone (GH) therapy has shown positive effects on the exercise behavior of mice with cerebellar atrophy, retains more Purkinje cells, and exhibits less DNA damage after GH intervention. Insulin-like growth factor 1 (IGF-1) is the downstream mediator of GH that participates in signaling and metabolic regulation for cell growth and modulation pathways, including SCA3-affected pathways. However, the underlying therapeutic mechanisms of GH or IGF-1 in SCA3 are not fully understood. In the present study, tissue-specific genome-scale metabolic network models for SCA3 transgenic mice were proposed based on RNA-seq. An integrative transcriptomic and metabolic network analysis of a SCA3 transgenic mouse model revealed that metabolic signaling pathways were activated to compensate for the metabolic remodeling caused by SCA3 genetic modifications. The effect of IGF-1 intervention on the pathology and balance of SCA3 disease was also explored. IGF-1 has been shown to invoke signaling pathways and improve mitochondrial function and glycolysis pathways to restore cellular functions. As one of the downregulated factors in SCA3 transgenic mice, IGF-1 could be a potential biomarker and therapeutic target.

Optimization of a modeling platform to predict oncogenes from genome-scale metabolic networks of non-small-cell lung cancers.

Cancer cell dysregulations result in the abnormal regulation of cellular metabolic pathways. By simulating this metabolic reprogramming using constraint‐based modeling approaches, oncogenes can be predicted, and this knowledge can be used in prognosis and treatment. We introduced a trilevel optimization problem describing metabolic reprogramming for inferring oncogenes. First, this study used RNA‐Seq expression data of lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) samples and their healthy counterparts to reconstruct tissue‐specific genome‐scale metabolic models and subsequently build the flux distribution pattern that provided a measure for the oncogene inference optimization problem for determining tumorigenesis. The platform detected 45 genes for LUAD and 84 genes for LUSC that lead to tumorigenesis. A high level of differentially expressed genes was not an essential factor for determining tumorigenesis. The platform indicated that pyruvate kinase (PKM), a well‐known oncogene with a low level of differential gene expression in LUAD and LUSC, had the highest fitness among the predicted oncogenes based on computation. By contrast, pyruvate kinase L/R (PKLR), an isozyme of PKM, had a high level of differential gene expression in both cancers. Phosphatidylserine synthase 1 (PTDSS1), an oncogene in LUAD, was inferred to have a low level of differential gene expression, and overexpression could significantly reduce survival probability. According to the factor analysis, PTDSS1 characteristics were close to those of the template, but they were unobvious in LUSC. Angiotensin‐converting enzyme 2 (ACE2) has recently garnered widespread interest as the SARS‐CoV‐2 virus receptor. Moreover, we determined that ACE2 is an oncogene of LUSC but not of LUAD. The platform developed in this study can identify oncogenes with low levels of differential expression and be used to identify potential therapeutic targets for cancer treatment.

Genome-scale metabolic network reconstruction of model animals as a platform for translational research

Genome-scale metabolic models (GEMs) are used extensively for analysis of mechanisms underlying human diseases and metabolic malfunctions. However, the lack of comprehensive and high-quality GEMs for model organisms restricts translational utilization of omics data accumulating from the use of various disease models. Here we present a unified platform of GEMs that covers five major model animals, including Mouse1 (Mus musculus), Rat1 (Rattus norvegicus), Zebrafish1 (Danio rerio), Fruitfly1 (Drosophila melanogaster), and Worm1 (Caenorhabditis elegans). These GEMs represent the most comprehensive coverage of the metabolic network by considering both orthology-based pathways and species-specific reactions. All GEMs can be interactively queried via the accompanying web portal Metabolic Atlas. Specifically, through integrative analysis of Mouse1 with RNA-sequencing data from brain tissues of transgenic mice we identified a coordinated up-regulation of lysosomal GM2 ganglioside and peptide degradation pathways which appears to be a signature metabolic alteration in Alzheimer's disease (AD) mouse models with a phenotype of amyloid precursor protein overexpression. This metabolic shift was further validated with proteomics data from transgenic mice and cerebrospinal fluid samples from human patients. The elevated lysosomal enzymes thus hold potential to be used as a biomarker for early diagnosis of AD. Taken together, we foresee that this evolving open-source platform will serve as an important resource to facilitate the development of systems medicines and translational biomedical applications.