Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease. Nat Commun. 2014;5:3083. [PMID: 24419221 DOI: 10.1038/ncomms4083]

Downregulation of PHGDH expression and hepatic serine level contribute to the development of fatty liver disease

OBJECTIVE

Supplementation with serine attenuates alcoholic fatty liver by regulating homocysteine metabolism and lipogenesis. However, little is known about serine metabolism in fatty liver disease (FLD). We aimed to investigate the changes in serine biosynthetic pathways in humans and animal models of fatty liver and their contribution to the development of FLD.

METHODS

High-fat diet (HFD)-induced steatosis and methionine-choline-deficient diet-induced steatohepatitis animal models were employed. Human serum samples were obtained from patients with FLD whose proton density fat fraction was estimated by magnetic resonance imaging. 3-Phosphoglycerate dehydrogenase (Phgdh)-knockout mouse embryonic fibroblasts (MEF) and transgenic mice overexpressing Phgdh (Tg-phgdh) were used to evaluate the role of serine metabolism in the development of FLD.

RESULTS

Expression of Phgdh was markedly reduced in the animal models. There were significant negative correlations of the serum serine with the liver fat fraction, serum alanine transaminase, and triglyceride levels among patients with FLD. Increased lipid accumulation and reduced NAD⁺ and SIRT1 activity were observed in Phgdh-knockout MEF and primary hepatocytes incubated with free fatty acids; these effects were reversed by overexpression of Phgdh. Tg-Phgdh mice showed significantly reduced hepatic triglyceride accumulation compared with wild-type littermates fed a HFD, which was accompanied by increased SIRT1 activity and reduced expression of lipogenic genes and proteins.

CONCLUSIONS

Human and experimental data suggest that reduced Phgdh expression and serine levels are closely associated with the development of FLD.

Metabolic Network Reconstructions to Predict Drug Targets and Off-Target Effects

The drug development pipeline has stalled because of the difficulty in identifying new drug targets while minimizing off-target effects. Computational methods, such as the use of metabolic network reconstructions, may provide a cost-effective platform to test new hypotheses for drug targets and prevent off-target effects. Here, we summarize available methods to identify drug targets and off-target effects using either reaction-centric, gene-centric, or metabolite-centric approaches with genome-scale metabolic network reconstructions.

Compartment and hub definitions tune metabolic networks for metabolomic interpretations

BACKGROUND

Metabolic networks represent all chemical reactions that occur between molecular metabolites in an organism's cells. They offer biological context in which to integrate, analyze, and interpret omic measurements, but their large scale and extensive connectivity present unique challenges. While it is practical to simplify these networks by placing constraints on compartments and hubs, it is unclear how these simplifications alter the structure of metabolic networks and the interpretation of metabolomic experiments.

RESULTS

We curated and adapted the latest systemic model of human metabolism and developed customizable tools to define metabolic networks with and without compartmentalization in subcellular organelles and with or without inclusion of prolific metabolite hubs. Compartmentalization made networks larger, less dense, and more modular, whereas hubs made networks larger, more dense, and less modular. When present, these hubs also dominated shortest paths in the network, yet their exclusion exposed the subtler prominence of other metabolites that are typically more relevant to metabolomic experiments. We applied the non-compartmental network without metabolite hubs in a retrospective, exploratory analysis of metabolomic measurements from 5 studies on human tissues. Network clusters identified individual reactions that might experience differential regulation between experimental conditions, several of which were not apparent in the original publications.

CONCLUSIONS

Exclusion of specific metabolite hubs exposes modularity in both compartmental and non-compartmental metabolic networks, improving detection of relevant clusters in omic measurements. Better computational detection of metabolic network clusters in large data sets has potential to identify differential regulation of individual genes, transcripts, and proteins.

Metabolic rearrangements in primary liver cancers: cause and consequences

Primary liver cancer (PLC) is the fourth most frequent cause of cancer-related death. The high mortality rates arise from late diagnosis and the limited accuracy of diagnostic and prognostic biomarkers. The liver is a major regulator, orchestrating the clearance of toxins, balancing glucose, lipid and amino acid uptake, managing whole-body metabolism and maintaining metabolic homeostasis. Tumour onset and progression is frequently accompanied by rearrangements of metabolic pathways, leading to dysregulation of metabolism. The limitation of current therapies targeting PLCs, such as hepatocellular carcinoma and cholangiocarcinoma, points towards the importance of deciphering this metabolic complexity. In this Review, we discuss the role of metabolic liver disruptions and the implications of these processes in PLCs, emphasizing their clinical relevance and value in early diagnosis and prognosis and as putative therapeutic targets. We also describe system biology approaches able to reconstruct the metabolic complexity of liver diseases. We also discuss whether metabolic rearrangements are a cause or consequence of PLCs, emphasizing the opportunity to clinically exploit the rewired metabolism. In line with this idea, we discuss circulating metabolites as promising biomarkers for PLCs.

Network Modeling Approaches and Applications to Unravelling Non-Alcoholic Fatty Liver Disease

Non-alcoholic fatty liver disease (NAFLD) is a progressive condition of the liver encompassing a range of pathologies including steatosis, non-alcoholic steatohepatitis (NASH), cirrhosis, and hepatocellular carcinoma. Research into this disease is imperative due to its rapid growth in prevalence, economic burden, and current lack of FDA approved therapies. NAFLD involves a highly complex etiology that calls for multi-tissue multi-omics network approaches to uncover the pathogenic genes and processes, diagnostic biomarkers, and potential therapeutic strategies. In this review, we first present a basic overview of disease pathogenesis, risk factors, and remaining knowledge gaps, followed by discussions of the need and concepts of multi-tissue multi-omics approaches, various network methodologies and application examples in NAFLD research. We highlight the findings that have been uncovered thus far including novel biomarkers, genes, and biological pathways involved in different stages of NAFLD, molecular connections between NAFLD and its comorbidities, mechanisms underpinning sex differences, and druggable targets. Lastly, we outline the future directions of implementing network approaches to further improve our understanding of NAFLD in order to guide diagnosis and therapeutics.

Complex I is bypassed during high intensity exercise

Human muscles are tailored towards ATP synthesis. When exercising at high work rates muscles convert glucose to lactate, which is less nutrient efficient than respiration. There is hence a trade-off between endurance and power. Metabolic models have been developed to study how limited catalytic capacity of enzymes affects ATP synthesis. Here we integrate an enzyme-constrained metabolic model with proteomics data from muscle fibers. We find that ATP synthesis is constrained by several enzymes. A metabolic bypass of mitochondrial complex I is found to increase the ATP synthesis rate per gram of protein compared to full respiration. To test if this metabolic mode occurs in vivo, we conduct a high resolved incremental exercise tests for five subjects. Their gas exchange at different work rates is accurately reproduced by a whole-body metabolic model incorporating complex I bypass. The study therefore shows how proteome allocation influences metabolism during high intensity exercise.

Reconstruction and analysis of a Kluyveromyces marxianus genome-scale metabolic model

BACKGROUND

Kluyveromyces marxianus is a thermotolerant yeast with multiple biotechnological potentials for industrial applications, which can metabolize a broad range of carbon sources, including less conventional sugars like lactose, xylose, arabinose and inulin. These phenotypic traits are sustained even up to 45 °C, what makes it a relevant candidate for industrial biotechnology applications, such as ethanol production. It is therefore of much interest to get more insight into the metabolism of this yeast. Recent studies suggested, that thermotolerance is achieved by reducing the number of growth-determining proteins or suppressing oxidative phosphorylation. Here we aimed to find related factors contributing to the thermotolerance of K. marxianus.

RESULTS

Here, we reported the first genome-scale metabolic model of Kluyveromyces marxianus, iSM996, using a publicly available Kluyveromyces lactis model as template. The model was manually curated and refined to include the missing species-specific metabolic capabilities. The iSM996 model includes 1913 reactions, associated with 996 genes and 1531 metabolites. It performed well to predict the carbon source utilization and growth rates under different growth conditions. Moreover, the model was coupled with transcriptomics data and used to perform simulations at various growth temperatures.

CONCLUSIONS

K. marxianus iSM996 represents a well-annotated metabolic model of thermotolerant yeast, which provides a new insight into theoretical metabolic profiles at different temperatures of K. marxianus. This could accelerate the integrative analysis of multi-omics data, leading to model-driven strain design and improvement.

Mitochondrial Dysfunction in the Transition from NASH to HCC

The liver constantly adapts to meet energy requirements of the whole body. Despite its remarkable adaptative capacity, prolonged exposure of liver cells to harmful environmental cues (such as diets rich in fat, sugar, and cholesterol) results in the development of chronic liver diseases (including non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH)) that can progress to hepatocellular carcinoma (HCC). The pathogenesis of these diseases is extremely complex, multifactorial, and poorly understood. Emerging evidence suggests that mitochondrial dysfunction or maladaptation contributes to detrimental effects on hepatocyte bioenergetics, reactive oxygen species (ROS) homeostasis, endoplasmic reticulum (ER) stress, inflammation, and cell death leading to NASH and HCC. The present review highlights the potential contribution of altered mitochondria function to NASH-related HCC and discusses how agents targeting this organelle could provide interesting treatment strategies for these diseases.

Serum metabolites detect the presence of advanced fibrosis in derivation and validation cohorts of patients with non-alcoholic fatty liver disease

OBJECTIVE

Non-invasive and accurate diagnostic tests for the screening of disease severity in non-alcoholic fatty liver disease (NAFLD) remain a major unmet need. Therefore, we aimed to examine if a combination of serum metabolites can accurately predict the presence of advanced fibrosis.

DESIGN

This is a cross-sectional analysis of a prospective derivation cohort including 156 well-characterised patients with biopsy-proven NAFLD and two validation cohorts, including (1) 142 patients assessed using MRI elastography (MRE) and(2) 59 patients with biopsy-proven NAFLD with untargeted serum metabolome profiling.

RESULTS

In the derivation cohort, 23 participants (15%) had advanced fibrosis and 32 of 652 analysed metabolites were significantly associated with advanced fibrosis after false-discovery rate adjustment. Among the top 10 metabolites, 8 lipids (5alpha-androstan-3beta monosulfate, pregnanediol-3-glucuronide, androsterone sulfate, epiandrosterone sulfate, palmitoleate, dehydroisoandrosterone sulfate, 5alpha-androstan-3beta disulfate, glycocholate), one amino acid (taurine) and one carbohydrate (fucose) were identified. The combined area under the receiver operating characteristic curve (AUROC) of the top 10 metabolite panel was higher than FIB--4 and NAFLD Fibrosis Score (NFS) for the detection of advanced fibrosis: 0.94 (95% CI 0.897 to 0.982) versus 0.78 (95% CI0.674 to 0.891), p=0.002 and versus 0.84 (95% CI 0.724 to 0.929), p=0.017, respectively. The AUROC of the top 10 metabolite panel remained excellent in the independent validation cohorts assessed by MRE or liver biopsy: c-statistic of 0.94 and 0.84, respectively.

CONCLUSION

A combination of 10 serum metabolites demonstrated excellent discriminatory ability for the detection of advanced fibrosis in an derivation and two independent validation cohorts with greater diagnostic accuracy than the FIB-4-index and NFS. This proof-of-concept study demonstrates that a non-invasive blood-based diagnostic test can provide excellent performance characteristics for the detection of advanced fibrosis.

Differences in Serum Amino Acid Phenotypes Among Patients with Diabetic Nephropathy, Hypertensive Nephropathy, and Chronic Nephritis

Background

We assessed levels of circulating amino acids in different etiologies of chronic kidney disease (CKD) and the association of amino acids with risk factors of CKD progression.

Material/Methods

High-performance liquid chromatography-based analysis was used to determine amino acid profiles in patients with diabetic nephropathy (DN, n=20), hypertensive nephropathy (HN, n=26), and chronic nephritis (CN, n=33), and in healthy controls (HC, n=25).

Results

All 3 types of CKD patients displayed decreased serum levels of serine, glycine, GABA, and tryptophan compared with healthy controls. Moreover, serine and tryptophan were positively correlated with glucose in DN cohorts. Total cholesterol was positively correlated with tryptophan levels in the DN cohort and negatively correlated with serine levels in the CN cohort. In the HN cohort, glycine was negatively correlated with triglyceride levels, and systolic blood pressure (SBP) was negatively correlated with GABA levels.

Conclusions

Patients with different etiologies of CKD have significantly different amino acids profiles, and this indicates specific supplementary nutritional needs in CKD patients.

Pan-cancer analysis of the metabolic reaction network

Metabolic reprogramming is considered a hallmark of malignant transformation. However, it is not clear whether the network of metabolic reactions expressed by cancers of different origin differ from each other or from normal human tissues. In this study, we reconstructed functional and connected genome-scale metabolic models for 917 primary tumor samples across 13 types based on the probability of expression for 3765 reference metabolic genes in the sample. This network-centric approach revealed that tumor metabolic networks are largely similar in terms of accounted reactions, despite diversity in the expression of the associated genes. On average, each network contained 4721 reactions, of which 74% were core reactions (present in >95% of all models). Whilst 99.3% of the core reactions were classified as housekeeping also in normal tissues, we identified reactions catalyzed by ARG2, RHAG, SLC6 and SLC16 family gene members, and PTGS1 and PTGS2 as core exclusively in cancer. These findings were subsequently replicated in an independent validation set of 3388 genome-scale metabolic models. The remaining 26% of the reactions were contextual reactions. Their inclusion was dependent in one case (GLS2) on the absence of TP53 mutations and in 94.6% of cases on differences in cancer types. This dependency largely resembled differences in expression patterns in the corresponding normal tissues, with some exceptions like the presence of the NANP-encoded reaction in tumors not from the female reproductive system or of the SLC5A9-encoded reaction in kidney-pancreatic-colorectal tumors. In conclusion, tumors expressed a metabolic network virtually overlapping the matched normal tissues, raising the possibility that metabolic reprogramming simply reflects cancer cell plasticity to adapt to varying conditions thanks to redundancy and complexity of the underlying metabolic networks. At the same time, the here uncovered exceptions represent a resource to identify selective liabilities of tumor metabolism.

Stoichiometric gene-to-reaction associations enhance model-driven analysis performance: Metabolic response to chronic exposure to Aldrin in prostate cancer

BACKGROUND

Genome-scale metabolic models (GSMM) integrating transcriptomics have been widely used to study cancer metabolism. This integration is achieved through logical rules that describe the association between genes, proteins, and reactions (GPRs). However, current gene-to-reaction formulation lacks the stoichiometry describing the transcript copies necessary to generate an active catalytic unit, which limits our understanding of how genes modulate metabolism. The present work introduces a new state-of-the-art GPR formulation that considers the stoichiometry of the transcripts (S-GPR). As case of concept, this novel gene-to-reaction formulation was applied to investigate the metabolic effects of the chronic exposure to Aldrin, an endocrine disruptor, on DU145 prostate cancer cells. To this aim we integrated the transcriptomic data from Aldrin-exposed and non-exposed DU145 cells through S-GPR or GPR into a human GSMM by applying different constraint-based-methods.

RESULTS

Our study revealed a significant improvement of metabolite consumption/production predictions when S-GPRs are implemented. Furthermore, our computational analysis unveiled important alterations in carnitine shuttle and prostaglandine biosynthesis in Aldrin-exposed DU145 cells that is supported by bibliographic evidences of enhanced malignant phenotype.

CONCLUSIONS

The method developed in this work enables a more accurate integration of gene expression data into model-driven methods. Thus, the presented approach is conceptually new and paves the way for more in-depth studies of aberrant cancer metabolism and other diseases with strong metabolic component with important environmental and clinical implications.

Development of a Plasma Screening Panel for Pediatric Nonalcoholic Fatty Liver Disease Using Metabolomics

Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in children, but diagnosis is challenging due to limited availability of noninvasive biomarkers. Machine learning applied to high-resolution metabolomics and clinical phenotype data offers a novel framework for developing a NAFLD screening panel in youth. Here, untargeted metabolomics by liquid chromatography-mass spectrometry was performed on plasma samples from a combined cross-sectional sample of children and adolescents ages 2-25 years old with NAFLD (n = 222) and without NAFLD (n = 337), confirmed by liver biopsy or magnetic resonance imaging. Anthropometrics, blood lipids, liver enzymes, and glucose and insulin metabolism were also assessed. A machine learning approach was applied to the metabolomics and clinical phenotype data sets, which were split into training and test sets, and included dimension reduction, feature selection, and classification model development. The selected metabolite features were the amino acids serine, leucine/isoleucine, and tryptophan; three putatively annotated compounds (dihydrothymine and two phospholipids); and two unknowns. The selected clinical phenotype variables were waist circumference, whole-body insulin sensitivity index (WBISI) based on the oral glucose tolerance test, and blood triglycerides. The highest performing classification model was random forest, which had an area under the receiver operating characteristic curve (AUROC) of 0.94, sensitivity of 73%, and specificity of 97% for detecting NAFLD cases. A second classification model was developed using the homeostasis model assessment of insulin resistance substituted for the WBISI. Similarly, the highest performing classification model was random forest, which had an AUROC of 0.92, sensitivity of 73%, and specificity of 94%. Conclusion: The identified screening panel consisting of both metabolomics and clinical features has promising potential for screening for NAFLD in youth. Further development of this panel and independent validation testing in other cohorts are warranted.

TCSBN: a database of tissue and cancer specific biological networks

Biological networks provide new opportunities for understanding the cellular biology in both health and disease states. We generated tissue specific integrated networks (INs) for liver, muscle and adipose tissues by integrating metabolic, regulatory and protein-protein interaction networks. We also generated human co-expression networks (CNs) for 46 normal tissues and 17 cancers to explore the functional relationships between genes as well as their relationships with biological functions, and investigate the overlap between functional and physical interactions provided by CNs and INs, respectively. These networks can be employed in the analysis of omics data, provide detailed insight into disease mechanisms by identifying the key biological components and eventually can be used in the development of efficient treatment strategies. Moreover, comparative analysis of the networks may allow for the identification of tissue-specific targets that can be used in the development of drugs with the minimum toxic effect to other human tissues. These context-specific INs and CNs are presented in an interactive website http://inetmodels.com without any limitation.

Quantifying the impact of public omics data

The amount of omics data in the public domain is increasing every year. Modern science has become a data-intensive discipline. Innovative solutions for data management, data sharing, and for discovering novel datasets are therefore increasingly required. In 2016, we released the first version of the Omics Discovery Index (OmicsDI) as a light-weight system to aggregate datasets across multiple public omics data resources. OmicsDI aggregates genomics, transcriptomics, proteomics, metabolomics and multiomics datasets, as well as computational models of biological processes. Here, we propose a set of novel metrics to quantify the attention and impact of biomedical datasets. A complete framework (now integrated into OmicsDI) has been implemented in order to provide and evaluate those metrics. Finally, we propose a set of recommendations for authors, journals and data resources to promote an optimal quantification of the impact of datasets.

Increasing amount of public omics data are important and valuable resources for the research community. Here, the authors develop a set of metrics to quantify the attention and impact of biomedical datasets and integrate them into the framework of Omics Discovery Index (OmicsDI).

Elucidating the Reprograming of Colorectal Cancer Metabolism Using Genome-Scale Metabolic Modeling

Colorectal cancer is the third most incidental cancer worldwide, and the response rate of current treatment for colorectal cancer is very low. Genome-scale metabolic models (GEMs) are systems biology platforms, and they had been used to assist researchers in understanding the metabolic alterations in different types of cancer. Here, we reconstructed a generic colorectal cancer GEM by merging 374 personalized GEMs from the Human Pathology Atlas and used it as a platform for systematic investigation of the difference between tumor and normal samples. The reconstructed model revealed the metabolic reprogramming in glutathione as well as the arginine and proline metabolism in response to tumor occurrence. In addition, six genes including ODC1, SMS, SRM, RRM2, SMOX, and SAT1 associated with arginine and proline metabolism were found to be key players in this metabolic alteration. We also investigated these genes in independent colorectal cancer patients and cell lines and found that many of these genes showed elevated level in colorectal cancer and exhibited adverse effect in patients. Therefore, these genes could be promising therapeutic targets for treatment of a specific colon cancer patient group.

The Potential Use of Metabolic Cofactors in Treatment of NAFLD

Non-alcoholic fatty liver disease (NAFLD) is caused by the imbalance between lipid deposition and lipid removal from the liver, and its global prevalence continues to increase dramatically. NAFLD encompasses a spectrum of pathological conditions including simple steatosis and non-alcoholic steatohepatitis (NASH), which can progress to cirrhosis and liver cancer. Even though there is a multi-disciplinary effort for development of a treatment strategy for NAFLD, there is not an approved effective medication available. Single or combined metabolic cofactors can be supplemented to boost the metabolic processes altered in NAFLD. Here, we review the dosage and usage of metabolic cofactors including l-carnitine, Nicotinamide riboside (NR), l-serine, and N-acetyl-l-cysteine (NAC) in human clinical studies to improve the altered biological functions associated with different human diseases. We also discuss the potential use of these substances in treatment of NAFLD and other metabolic diseases including neurodegenerative and cardiovascular diseases of which pathogenesis is linked to mitochondrial dysfunction.

Genome-scale metabolic model of the rat liver predicts effects of diet restriction

Mapping network analysis in cells and tissues can provide insights into metabolic adaptations to changes in external environment, pathological conditions, and nutrient deprivation. Here, we reconstructed a genome-scale metabolic network of the rat liver that will allow for exploration of systems-level physiology. The resulting in silico model (iRatLiver) contains 1,882 reactions, 1,448 metabolites, and 994 metabolic genes. We then used this model to characterize the response of the liver’s energy metabolism to a controlled perturbation in diet. Transcriptomics data were collected from the livers of Sprague Dawley rats at 4 or 14 days of being subjected to 15%, 30%, or 60% diet restriction. These data were integrated with the iRatLiver model to generate condition-specific metabolic models, allowing us to explore network differences under each condition. We observed different pathway usage between early and late time points. Network analysis identified several highly connected “hub” genes (Pklr, Hadha, Tkt, Pgm1, Tpi1, and Eno3) that showed differing trends between early and late time points. Taken together, our results suggest that the liver’s response varied with short- and long-term diet restriction. More broadly, we anticipate that the iRatLiver model can be exploited further to study metabolic changes in the liver under other conditions such as drug treatment, infection, and disease.

Assessing key decisions for transcriptomic data integration in biochemical networks

To gain insights into complex biological processes, genome-scale data (e.g., RNA-Seq) are often overlaid on biochemical networks. However, many networks do not have a one-to-one relationship between genes and network edges, due to the existence of isozymes and protein complexes. Therefore, decisions must be made on how to overlay data onto networks. For example, for metabolic networks, these decisions include (1) how to integrate gene expression levels using gene-protein-reaction rules, (2) the approach used for selection of thresholds on expression data to consider the associated gene as "active", and (3) the order in which these steps are imposed. However, the influence of these decisions has not been systematically tested. We compared 20 decision combinations using a transcriptomic dataset across 32 tissues and showed that definition of which reaction may be considered as active (i.e., reactions of the genome-scale metabolic network with a non-zero expression level after overlaying the data) is mainly influenced by thresholding approach used. To determine the most appropriate decisions, we evaluated how these decisions impact the acquisition of tissue-specific active reaction lists that recapitulate organ-system tissue groups. These results will provide guidelines to improve data analyses with biochemical networks and facilitate the construction of context-specific metabolic models.

Therapeutic Effects of Amino Acids in Liver Diseases: Current Studies and Future Perspectives

Hepatocellular carcinoma (HCC) is the most common primary malignant tumor of the liver and the third most common cause of cancer-related death worldwide. HCC is caused by infection of hepatitis B/C virus and liver dysfunctions, such as alcoholic liver disease, nonalcoholic fatty liver disease, and cirrhosis. Amino acids are organic substances containing amine and carboxylic acid functional groups. There are over 700 kinds of amino acids in nature, but only about 20 of them are used to synthesize proteins in cells. Liver is an important organ for protein synthesis, degradation and detoxification as well as amino acid metabolism. In the liver, there are abundant non-essential amino acids, such as alanine, aspartate, glutamate, glycine, and serine and essential amino acids, such as histidine and threonine. These amino acids are involved in various cellular metabolisms, the synthesis of lipids and nucleotides as well as detoxification reactions. Understanding the role of amino acids in the pathogenesis of liver and the effects of amino acid intake on liver disease can be a promising strategy for the prevention and treatment of liver disease. In this review, we describe the biochemical properties and functions of amino acids and to review how they have been applied to treatment of liver diseases.

Integration of Metabolomic and Other Omics Data in Population-Based Study Designs: An Epidemiological Perspective

It is not controversial that study design considerations and challenges must be addressed when investigating the linkage between single omic measurements and human phenotypes. It follows that such considerations are just as critical, if not more so, in the context of multi-omic studies. In this review, we discuss (1) epidemiologic principles of study design, including selection of biospecimen source(s) and the implications of the timing of sample collection, in the context of a multi-omic investigation, and (2) the strengths and limitations of various techniques of data integration across multi-omic data types that may arise in population-based studies utilizing metabolomic data.

Glycine Metabolism and Its Alterations in Obesity and Metabolic Diseases

Glycine is the proteinogenic amino-acid of lowest molecular weight, harboring a hydrogen atom as a side-chain. In addition to being a building-block for proteins, glycine is also required for multiple metabolic pathways, such as glutathione synthesis and regulation of one-carbon metabolism. Although generally viewed as a non-essential amino-acid, because it can be endogenously synthesized to a certain extent, glycine has also been suggested as a conditionally essential amino acid. In metabolic disorders associated with obesity, type 2 diabetes (T2DM), and non-alcoholic fatty liver disease (NAFLDs), lower circulating glycine levels have been consistently observed, and clinical studies suggest the existence of beneficial effects induced by glycine supplementation. The present review aims at synthesizing the recent advances in glycine metabolism, pinpointing its main metabolic pathways, identifying the causes leading to glycine deficiency-especially in obesity and associated metabolic disorders-and evaluating the potential benefits of increasing glycine availability to curb the progression of obesity and obesity-related metabolic disturbances. This study focuses on the importance of diet, gut microbiota, and liver metabolism in determining glycine availability in obesity and associated metabolic disorders.

Current status and applications of genome-scale metabolic models

Genome-scale metabolic models (GEMs) computationally describe gene-protein-reaction associations for entire metabolic genes in an organism, and can be simulated to predict metabolic fluxes for various systems-level metabolic studies. Since the first GEM for Haemophilus influenzae was reported in 1999, advances have been made to develop and simulate GEMs for an increasing number of organisms across bacteria, archaea, and eukarya. Here, we review current reconstructed GEMs and discuss their applications, including strain development for chemicals and materials production, drug targeting in pathogens, prediction of enzyme functions, pan-reactome analysis, modeling interactions among multiple cells or organisms, and understanding human diseases.

Broad Views of Non-alcoholic Fatty Liver Disease

Multi-omics multi-tissue data are used to interpret genome-wide association study results from mice to identify key driver genes of non-alcoholic fatty liver disease. Non-alcoholic fatty liver disease (NAFLD) is the accumulation of fat (steatosis) in the liver due to causes other than excessive alcohol consumption. The disease may progress to more severe forms of liver diseases, including non-alcoholic steatohepatitis, cirrhosis, and hepatocellular carcinoma. In this issue of Cell Systems, Krishnan et al. (2018) reveal mechanisms underlying NAFLD by generating multi-omics data using liver and adipose tissues obtained from the Hybrid Mouse Diversity Panel, consisting of 113 mouse strains with various degrees of NAFLD. The study identified key driver genes of NAFLD that can be used in the development of efficient treatment strategies and illustrates the potential utility of systematic analysis of multi-layer biological networks.

Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.

Identifying and targeting cancer-specific metabolism with network-based drug target prediction

Background

Metabolic rewiring allows cancer cells to sustain high proliferation rates. Thus, targeting only the cancer-specific cellular metabolism will safeguard healthy tissues.

Methods

We developed the very efficient FASTCORMICS RNA-seq workflow (rFASTCORMICS) to build 10,005 high-resolution metabolic models from the TCGA dataset to capture metabolic rewiring strategies in cancer cells. Colorectal cancer (CRC) was used as a test case for a repurposing workflow based on rFASTCORMICS.

Findings

Alternative pathways that are not required for proliferation or survival tend to be shut down and, therefore, tumours display cancer-specific essential genes that are significantly enriched for known drug targets. We identified naftifine, ketoconazole, and mimosine as new potential CRC drugs, which were experimentally validated.

Interpretation

The here presented rFASTCORMICS workflow successfully reconstructs a metabolic model based on RNA-seq data and successfully predicted drug targets and drugs not yet indicted for colorectal cancer.

Fund

This study was supported by the University of Luxembourg (IRP grant scheme; R-AGR-0755-12), the Luxembourg National Research Fund (FNR PRIDE PRIDE15/10675146/CANBIO), the Fondation Cancer (Luxembourg), the European Union‘s Horizon2020 research and innovation programme under the Marie Sklodowska- Curie grant agreement No 642295 (MEL-PLEX), and the German Federal Ministry of Education and Research (BMBF) within the project MelanomSensitivity (BMBF/BM/7643621).

The Fate of Glutamine in Human Metabolism. The Interplay with Glucose in Proliferating Cells

Genome-scale models of metabolism (GEM) are used to study how metabolism varies in different physiological conditions. However, the great number of reactions involved in GEM makes it difficult to understand these variations. In order to have a more understandable tool, we developed a reduced metabolic model of central carbon and nitrogen metabolism, C2M2N with 77 reactions, 54 internal metabolites, and 3 compartments, taking into account the actual stoichiometry of the reactions, including the stoichiometric role of the cofactors and the irreversibility of some reactions. In order to model oxidative phosphorylation (OXPHOS) functioning, the proton gradient through the inner mitochondrial membrane is represented by two pseudometabolites DPH (∆pH) and DPSI (∆ψ). To illustrate the interest of such a reduced and quantitative model of metabolism in mammalian cells, we used flux balance analysis (FBA) to study all the possible fates of glutamine in metabolism. Our analysis shows that glutamine can supply carbon sources for cell energy production and can be used as carbon and nitrogen sources to synthesize essential metabolites. Finally, we studied the interplay between glucose and glutamine for the formation of cell biomass according to ammonia microenvironment. We then propose a quantitative analysis of the Warburg effect.

The Fate of Glutamine in Human Metabolism. The Interplay with Glucose in Proliferating Cells

Genome-scale models of metabolism (GEM) are used to study how metabolism varies in different physiological conditions. However, the great number of reactions involved in GEM makes it difficult to understand these variations. In order to have a more understandable tool, we developed a reduced metabolic model of central carbon and nitrogen metabolism, C2M2N with 77 reactions, 54 internal metabolites, and 3 compartments, taking into account the actual stoichiometry of the reactions, including the stoichiometric role of the cofactors and the irreversibility of some reactions. In order to model oxidative phosphorylation (OXPHOS) functioning, the proton gradient through the inner mitochondrial membrane is represented by two pseudometabolites DPH (∆pH) and DPSI (∆ψ). To illustrate the interest of such a reduced and quantitative model of metabolism in mammalian cells, we used flux balance analysis (FBA) to study all the possible fates of glutamine in metabolism. Our analysis shows that glutamine can supply carbon sources for cell energy production and can be used as carbon and nitrogen sources to synthesize essential metabolites. Finally, we studied the interplay between glucose and glutamine for the formation of cell biomass according to ammonia microenvironment. We then propose a quantitative analysis of the Warburg effect.

A benchmark-driven approach to reconstruct metabolic networks for studying cancer metabolism

Genome-scale metabolic modeling has emerged as a promising way to study the metabolic alterations underlying cancer by identifying novel drug targets and biomarkers. To date, several computational methods have been developed to integrate high-throughput data with existing human metabolic reconstructions to generate context-specific cancer metabolic models. Despite a number of studies focusing on benchmarking the context-specific algorithms, no quantitative assessment has been made to compare the predictive performance of these methods. Here, we integrated various and different datasets used in previous works to design a quantitative platform to examine functional and consistency performance of several existing genome-scale cancer modeling approaches. Next, we used the results obtained here to develop a method for the reconstruction of context-specific metabolic models. We then compared the predictive power and consistency of networks generated by our method to other computational approaches investigated here. Our results showed a satisfactory performance of the developed method in most of the benchmarks. This benchmarking platform is of particular use in algorithm selection and assessing the performance of newly developed algorithms. More importantly, it can serve as guidelines for designing and developing new methods focusing on weaknesses and strengths of existing algorithms.

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.

Mechanistic identification of biofluid metabolite changes as markers of acetaminophen-induced liver toxicity in rats

Acetaminophen (APAP) is the most commonly used analgesic and antipyretic drug in the world. Yet, it poses a major risk of liver injury when taken in excess of the therapeutic dose. Current clinical markers do not detect the early onset of liver injury associated with excess APAP-information that is vital to reverse injury progression through available therapeutic interventions. Hence, several studies have used transcriptomics, proteomics, and metabolomics technologies, both independently and in combination, in an attempt to discover potential early markers of liver injury. However, the casual relationship between these observations and their relation to the APAP mechanism of liver toxicity are not clearly understood. Here, we used Sprague-Dawley rats orally gavaged with a single dose of 2 g/kg of APAP to collect tissue samples from the liver and kidney for transcriptomic analysis and plasma and urine samples for metabolomic analysis. We developed and used a multi-tissue, metabolism-based modeling approach to integrate these data, characterize the effect of excess APAP levels on liver metabolism, and identify a panel of plasma and urine metabolites that are associated with APAP-induced liver toxicity. Our analyses, which indicated that pathways involved in nucleotide-, lipid-, and amino acid-related metabolism in the liver were most strongly affected within 10 h following APAP treatment, identified a list of potential metabolites in these pathways that could serve as plausible markers of APAP-induced liver injury. Our approach identifies toxicant-induced changes in endogenous metabolism, is applicable to other toxicants based on transcriptomic data, and provides a mechanistic framework for interpreting metabolite alterations.

From sugar to liver fat and public health: systems biology driven studies in understanding non-alcoholic fatty liver disease pathogenesis

Non-alcoholic fatty liver disease (NAFLD) is now a major public health concern with an estimated prevalence of 25-30% of adults in many countries. Strongly associated with obesity and the metabolic syndrome, the pathogenesis of NAFLD is dependent on complex interactions between genetic and environmental factors that are not completely understood. Weight loss through diet and lifestyle modification underpins clinical management; however, the roles of individual dietary nutrients (e.g. saturated and n-3 fatty acids; fructose, vitamin D, vitamin E) in the pathogenesis or treatment of NAFLD are only partially understood. Systems biology offers valuable interdisciplinary methods that are arguably ideal for application to the studying of chronic diseases such as NAFLD, and the roles of nutrition and diet in their molecular pathogenesis. Although present in silico models are incomplete, computational tools are rapidly evolving and human metabolism can now be simulated at the genome scale. This paper will review NAFLD and its pathogenesis, including the roles of genetics and nutrition in the development and progression of disease. In addition, the paper introduces the concept of systems biology and reviews recent work utilising genome-scale metabolic networks and developing multi-scale models of liver metabolism relevant to NAFLD. A future is envisioned where individual genetic, proteomic and metabolomic information can be integrated computationally with clinical data, yielding mechanistic insight into the pathogenesis of chronic diseases such as NAFLD, and informing personalised nutrition and stratified medicine approaches for improving prognosis.

Non-Alcoholic Fatty Liver Disease and Risk of Incident Type 2 Diabetes: Role of Circulating Branched-Chain Amino Acids

Non-alcoholic fatty liver disease (NAFLD) is likely to be associated with elevated plasma branched-chain amino acids (BCAAs) and may precede the development of type 2 diabetes (T2D). We hypothesized that BCAAs may be involved in the pathogenesis of T2D attributable to NAFLD and determined the extent to which plasma BCAAs influence T2D development in NAFLD. We evaluated cross-sectional associations of NAFLD with fasting plasma BCAAs (nuclear magnetic resonance spectroscopy), and prospectively determined the extent to which the influence of NAFLD on incident T2D is attributable to BCAA elevations. In the current study, 5791 Prevention of REnal and Vascular ENd-stage Disease (PREVEND) cohort participants without T2D at baseline were included. Elevated fatty liver index (FLI) ≥60, an algorithm based on triglycerides, gamma-glutamyltransferase, body mass index (BMI) and waist circumference, was used as proxy of NAFLD. Elevated FLI ≥ 60 was present in 1671 (28.9%) participants. Cross-sectionally, BCAAs were positively associated with FLI ≥ 60 (β = 0.208, p < 0.001). During a median follow-up of 7.3 years, 276 participants developed T2D, of which 194 (70.2%) had an FLI ≥ 60 (log-rank test, p < 0.001). Cox regression analyses revealed that both FLI ≥60 (hazard ratio (HR) 3.46, 95% CI 2.45⁻4.87, p < 0.001) and higher BCAA levels (HR 1.19, 95% CI 1.03⁻1.37, p = 0.01) were positively associated with incident T2D. Mediation analysis showed that the association of FLI with incident T2D was in part attributable to elevated BCAAs (proportion mediated 19.6%). In conclusion, both elevated FLI and elevated plasma BCAA levels are associated with risk of incident T2D. The association of NAFLD with T2D development seems partly mediated by elevated BCAAs.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

The OMICs Window into Nonalcoholic Fatty Liver Disease (NAFLD)

Nonalcoholic fatty liver disease (NAFLD) is a common cause of hepatic abnormalities worldwide. Nonalcoholic steatohepatitis (NASH) is part of the spectrum of NAFLD and leads to progressive liver disease, such as cirrhosis and hepatocellular carcinoma. In NASH patient, fibrosis represents the major predictor of liver-related mortality; therefore, it is important to have an early and accurate diagnosis of NASH. The current gold standard for the diagnosis of NASH is still liver biopsy. The development of biomarkers able to predict disease severity, prognosis, as well as response to therapy without the need for a biopsy is the focus of most up-to-date genomic, transcriptomic, proteomic, and metabolomic research. In the future, patients might be diagnosed and treated according to their molecular signatures. In this short review, we discuss how information from genomics, proteomics, and metabolomics contribute to the understanding of NAFLD pathogenesis.

Metabolic Modeling of Human Gut Microbiota on a Genome Scale: An Overview

There is growing interest in the metabolic interplay between the gut microbiome and host metabolism. Taxonomic and functional profiling of the gut microbiome by next-generation sequencing (NGS) has unveiled substantial richness and diversity. However, the mechanisms underlying interactions between diet, gut microbiome and host metabolism are still poorly understood. Genome-scale metabolic modeling (GSMM) is an emerging approach that has been increasingly applied to infer diet⁻microbiome, microbe⁻microbe and host⁻microbe interactions under physiological conditions. GSMM can, for example, be applied to estimate the metabolic capabilities of microbes in the gut. Here, we discuss how meta-omics datasets such as shotgun metagenomics, can be processed and integrated to develop large-scale, condition-specific, personalized microbiota models in healthy and disease states. Furthermore, we summarize various tools and resources available for metagenomic data processing and GSMM, highlighting the experimental approaches needed to validate the model predictions.

Genetic basis for plasma amino acid concentrations based on absolute quantification: a genome-wide association study in the Japanese population

To assess the use of plasma free amino acids (PFAAs) as biomarkers for metabolic disorders, it is essential to identify genetic factors that influence PFAA concentrations. PFAA concentrations were absolutely quantified by liquid chromatography–mass spectrometry using plasma samples from 1338 Japanese individuals, and genome-wide quantitative trait locus (QTL) analysis was performed for the concentrations of 21 PFAAs. We next conducted a conditional QTL analysis using the concentration of each PFAA adjusted by the other 20 PFAAs as covariates to elucidate genetic determinants that influence PFAA concentrations. We identified eight genes that showed a significant association with PFAA concentrations, of which two, SLC7A2 and PKD1L2, were identified. SLC7A2 was associated with the plasma levels of arginine and ornithine, and PKD1L2 with the level of glycine. The significant associations of these two genes were revealed in the conditional QTL analysis, but a significant association between serine and the CPS1 gene disappeared when glycine was used as a covariate. We demonstrated that conditional QTL analysis is useful for determining the metabolic pathways predominantly used for PFAA metabolism. Our findings will help elucidate the physiological roles of genetic components that control the metabolism of amino acids.

SLIMEr: probing flexibility of lipid metabolism in yeast with an improved constraint-based modeling framework

BACKGROUND

A recurrent problem in genome-scale metabolic models (GEMs) is to correctly represent lipids as biomass requirements, due to the numerous of possible combinations of individual lipid species and the corresponding lack of fully detailed data. In this study we present SLIMEr, a formalism for correctly representing lipid requirements in GEMs using commonly available experimental data.

RESULTS

SLIMEr enhances a GEM with mathematical constructs where we Split Lipids Into Measurable Entities (SLIME reactions), in addition to constraints on both the lipid classes and the acyl chain distribution. By implementing SLIMEr on the consensus GEM of Saccharomyces cerevisiae, we can represent accurate amounts of lipid species, analyze the flexibility of the resulting distribution, and compute the energy costs of moving from one metabolic state to another.

CONCLUSIONS

The approach shows potential for better understanding lipid metabolism in yeast under different conditions. SLIMEr is freely available at https://github.com/SysBioChalmers/SLIMEr .

1-Deoxysphingolipids

Sphingolipids (SLs) are fundamental components of eukaryotic cells. 1-Deoxysphingolipids differ structurally from canonical SLs as they lack the essential C1-OH group. Consequently, 1-deoxysphingolipids cannot be converted to complex sphingolipids and are not degraded over the canonical catabolic pathways. Pathologically elevated 1-deoxySLs are involved in several disease conditions. Within this review, we will provide an up-to-date overview on the metabolic, physiological and pathophysiological aspects of this enigmatic class of "headless" sphingolipids.

Pyruvate kinase L/R is a regulator of lipid metabolism and mitochondrial function

The pathogenesis of non-alcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC) has been associated with altered expression of liver-specific genes including pyruvate kinase liver and red blood cell (PKLR), patatin-like phospholipase domain containing 3 (PNPLA3) and proprotein convertase subtilisin/kexin type 9 (PCSK9). Here, we inhibited and overexpressed the expression of these three genes in HepG2 cells, generated RNA-seq data before and after perturbation and revealed the altered global biological functions with the modulation of these genes using integrated network (IN) analysis. We found that modulation of these genes effects the total triglycerides levels within the cells and viability of the cells. Next, we generated IN for HepG2 cells, identified reporter transcription factors based on IN and found that the modulation of these genes affects key metabolic pathways associated with lipid metabolism (steroid biosynthesis, PPAR signalling pathway, fatty acid synthesis and oxidation) and cancer development (DNA replication, cell cycle and p53 signalling) involved in the progression of NAFLD and HCC. Finally, we observed that inhibition of PKLR lead to decreased glucose uptake and decreased mitochondrial activity in HepG2 cells. Hence, our systems level analysis indicated that PKLR can be targeted for development efficient treatment strategy for NAFLD and HCC.

Genome-Scale Metabolic Modeling from Yeast to Human Cell Models of Complex Diseases: Latest Advances and Challenges

Genome-scale metabolic models (GEMs) are mathematical models that enable systematic analysis of metabolism. This modeling concept has been applied to study the metabolism of many organisms including the eukaryal model organism, the yeast Saccharomyces cerevisiae, that also serves as an important cell factory for production of fuels and chemicals. With the application of yeast GEMs, our knowledge of metabolism is increasing. Therefore, GEMs have also been used for modeling human cells to study metabolic diseases. Here we introduce the concept of GEMs and provide a protocol for reconstructing GEMs. Besides, we show the historic development of yeast GEMs and their applications. Also, we review human GEMs as well as their uses in the studies of complex diseases.

A Model Integration Pipeline for the Improvement of Human Genome-Scale Metabolic Reconstructions

Metabolism has been a major field of study in the last years, mainly due to its importance in understanding cell physiology and certain disease phenotypes due to its deregulation. Genome-scale metabolic models (GSMMs) have been established as important tools to help achieve a better understanding of human metabolism. Towards this aim, advances in systems biology and bioinformatics have allowed the reconstruction of several human GSMMs, although some limitations and challenges remain, such as the lack of external identifiers for both metabolites and reactions. A pipeline was developed to integrate multiple GSMMs, starting by retrieving information from the main human GSMMs and evaluating the presence of external database identifiers and annotations for both metabolites and reactions. Information from metabolites was included into a graph database with omics data repositories, allowing clustering of metabolites through their similarity regarding database cross-referencing. Metabolite annotation of several older GSMMs was enriched, allowing the identification and integration of common entities. Using this information, as well as other metrics, we successfully integrated reactions from these models. These methods can be leveraged towards the creation of a unified consensus model of human metabolism.

A simplified metabolic network reconstruction to promote understanding and development of flux balance analysis tools

GEnome-scale Network REconstructions (GENREs) mathematically describe metabolic reactions of an organism or a specific cell type. GENREs can be used with a number of constraint-based reconstruction and analysis (COBRA) methods to make computational predictions on how a system changes in different environments. We created a simplified GENRE (referred to as iSIM) that captures central energy metabolism with nine metabolic reactions to illustrate the use of and promote the understanding of GENREs and constraint-based methods. We demonstrate the simulation of single and double gene deletions, flux variability analysis (FVA), and test a number of metabolic tasks with the GENRE. Code to perform these analyses is provided in Python, R, and MATLAB. Finally, with iSIM as a guide, we demonstrate how inaccuracies in GENREs can limit their use in the interrogation of energy metabolism.

Metabolic network-based stratification of hepatocellular carcinoma reveals three distinct tumor subtypes

Hepatocellular carcinoma (HCC) is one of the most frequent forms of liver cancer, and effective treatment methods are limited due to tumor heterogeneity. There is a great need for comprehensive approaches to stratify HCC patients, gain biological insights into subtypes, and ultimately identify effective therapeutic targets. We stratified HCC patients and characterized each subtype using transcriptomics data, genome-scale metabolic networks and network topology/controllability analysis. This comprehensive systems-level analysis identified three distinct subtypes with substantial differences in metabolic and signaling pathways reflecting at genomic, transcriptomic, and proteomic levels. These subtypes showed large differences in clinical survival associated with altered kynurenine metabolism, WNT/β-catenin-associated lipid metabolism, and PI3K/AKT/mTOR signaling. Integrative analyses indicated that the three subtypes rely on alternative enzymes (e.g., ACSS1/ACSS2/ACSS3, PKM/PKLR, ALDOB/ALDOA, MTHFD1L/MTHFD2/MTHFD1) to catalyze the same reactions. Based on systems-level analysis, we identified 8 to 28 subtype-specific genes with pivotal roles in controlling the metabolic network and predicted that these genes may be targeted for development of treatment strategies for HCC subtypes by performing in silico analysis. To validate our predictions, we performed experiments using HepG2 cells under normoxic and hypoxic conditions and observed opposite expression patterns between genes expressed in high/moderate/low-survival tumor groups in response to hypoxia, reflecting activated hypoxic behavior in patients with poor survival. In conclusion, our analyses showed that the heterogeneous HCC tumors can be stratified using a metabolic network-driven approach, which may also be applied to other cancer types, and this stratification may have clinical implications to drive the development of precision medicine.

Large-Scale Modeling Approach Reveals Functional Metabolic Shifts during Hepatic Differentiation

Being able to explore the metabolism of broad metabolizing cells is of critical importance in many research fields. This article presents an original modeling solution combining metabolic network and omics data to identify modulated metabolic pathways and changes in metabolic functions occurring during differentiation of a human hepatic cell line (HepaRG). Our results confirm the activation of hepato-specific functionalities and newly evidence modulation of other metabolic pathways, which could not be evidenced from transcriptomic data alone. Our method takes advantage of the network structure to detect changes in metabolic pathways that do not have gene annotations and exploits flux analyses techniques to identify activated metabolic functions. Compared to the usual cell-specific metabolic network reconstruction approaches, it limits false predictions by considering several possible network configurations to represent one phenotype rather than one arbitrarily selected network. Our approach significantly enhances the comprehensive and functional assessment of cell metabolism, opening further perspectives to investigate metabolic shifts occurring within various biological contexts.

Aberrant Metabolism in Hepatocellular Carcinoma Provides Diagnostic and Therapeutic Opportunities

Hepatocellular carcinoma (HCC) accounts for over 80% of liver cancer cases and is highly malignant, recurrent, drug-resistant, and often diagnosed in the advanced stage. It is clear that early diagnosis and a better understanding of molecular mechanisms contributing to HCC progression is clinically urgent. Metabolic alterations clearly characterize HCC tumors. Numerous clinical parameters currently used to assess liver functions reflect changes in both enzyme activity and metabolites. Indeed, differences in glucose and acetate utilization are used as a valid clinical tool for stratifying patients with HCC. Moreover, increased serum lactate can distinguish HCC from normal subjects, and serum lactate dehydrogenase is used as a prognostic indicator for HCC patients under therapy. Currently, the emerging field of metabolomics that allows metabolite analysis in biological fluids is a powerful method for discovering new biomarkers. Several metabolic targets have been identified by metabolomics approaches, and these could be used as biomarkers in HCC. Moreover, the integration of different omics approaches could provide useful information on the metabolic pathways at the systems level. In this review, we provided an overview of the metabolic characteristics of HCC considering also the reciprocal influences between the metabolism of cancer cells and their microenvironment. Moreover, we also highlighted the interaction between hepatic metabolite production and their serum revelations through metabolomics researches.