Reconciled rat and human metabolic networks for comparative toxicogenomics and biomarker predictions

The laboratory rat has been used as a surrogate to study human biology for more than a century. Here we present the first genome-scale network reconstruction of Rattus norvegicus metabolism, iRno, and a significantly improved reconstruction of human metabolism, iHsa. These curated models comprehensively capture metabolic features known to distinguish rats from humans including vitamin C and bile acid synthesis pathways. After reconciling network differences between iRno and iHsa, we integrate toxicogenomics data from rat and human hepatocytes, to generate biomarker predictions in response to 76 drugs. We validate comparative predictions for xanthine derivatives with new experimental data and literature-based evidence delineating metabolite biomarkers unique to humans. Our results provide mechanistic insights into species-specific metabolism and facilitate the selection of biomarkers consistent with rat and human biology. These models can serve as powerful computational platforms for contextualizing experimental data and making functional predictions for clinical and basic science applications.

[1]  Monica L. Mo,et al.  Global reconstruction of the human metabolic network based on genomic and bibliomic data , 2007, Proceedings of the National Academy of Sciences.

[2]  T. Ozawa,et al.  A missense mutation of L-gulono-gamma-lactone oxidase causes the inability of scurvy-prone osteogenic disorder rats to synthesize L-ascorbic acid. , 1992, The Journal of biological chemistry.

[3]  R. Matthiesen,et al.  Differential expression of alpha-2,3-sialyltransferases and alpha-1,3/4-fucosyltransferases regulates the levels of sialyl Lewis a and sialyl Lewis x in gastrointestinal carcinoma cells. , 2010, The international journal of biochemistry & cell biology.

[4]  K. Maurer,et al.  Third national health and nutrition examination survey , 1985 .

[5]  Bernhard O. Palsson,et al.  Optimizing genome-scale network reconstructions , 2014, Nature Biotechnology.

[6]  Yun‐Sil Lee,et al.  Acetaminophen selectively suppresses peripheral prostaglandin E2 release and increases COX-2 gene expression in a clinical model of acute inflammation , 2007, PAIN.

[7]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[8]  Hiroshi Yamada,et al.  Open TG-GATEs: a large-scale toxicogenomics database , 2014, Nucleic Acids Res..

[9]  Gregory D. Schuler,et al.  Database resources of the National Center for Biotechnology Information , 2008, Nucleic Acids Res..

[10]  G. Benga,et al.  Amino acid composition of rat and human liver microsomes in normal and pathological conditions , 1995, Bioscience reports.

[11]  Aarash Bordbar,et al.  Pharmacogenomic and clinical data link non-pharmacokinetic metabolic dysregulation to drug side effect pathogenesis , 2015, Nature Communications.

[12]  Jason A. Papin,et al.  Comparative Metabolic Systems Analysis of Pathogenic Burkholderia , 2013, Journal of bacteriology.

[13]  Jason A Papin,et al.  Linking genome-scale metabolic modeling and genome annotation. , 2013, Methods in molecular biology.

[14]  E. Fukusaki,et al.  GC‐MS‐based metabolomics reveals mechanism of action for hydrazine induced hepatotoxicity in rats , 2011, Journal of applied toxicology : JAT.

[15]  Susumu Goto,et al.  Data, information, knowledge and principle: back to metabolism in KEGG , 2013, Nucleic Acids Res..

[16]  Michael D. Waters,et al.  Toxicogenomics and systems toxicology: aims and prospects , 2004, Nature Reviews Genetics.

[17]  M. Lewis,et al.  Structure–function relationships in 3α-hydroxysteroid dehydrogenases: a comparison of the rat and human isoforms , 2003, The Journal of Steroid Biochemistry and Molecular Biology.

[18]  Intawat Nookaew,et al.  The RAVEN Toolbox and Its Use for Generating a Genome-scale Metabolic Model for Penicillium chrysogenum , 2013, PLoS Comput. Biol..

[19]  N. Christensen,et al.  Theophylline enhances glucose recovery after hypoglycemia in healthy man and in type I diabetic patients. , 1994, Metabolism: clinical and experimental.

[20]  S. Dong,et al.  Dedifferentiation of conditionally immortalized hepatocytes with long-term in vitro passage , 2000, Experimental & Molecular Medicine.

[21]  E. Ruppin,et al.  Predicting metabolic biomarkers of human inborn errors of metabolism , 2009, Molecular systems biology.

[22]  Yizeng Liang,et al.  Joint MS-based platforms for comprehensive comparison of rat plasma and serum metabolic profiling. , 2014, Biomedical chromatography : BMC.

[23]  Evan Bolton,et al.  Database resources of the National Center for Biotechnology Information , 2017, Nucleic Acids Res..

[24]  L. Simchowitz,et al.  Mechanism of action of colchicine in acute urate crystal-induced arthritis. , 1979, The Journal of clinical investigation.

[25]  Jason A. Papin,et al.  Functional integration of a metabolic network model and expression data without arbitrary thresholding , 2011, Bioinform..

[26]  E. Want,et al.  Systemic gut microbial modulation of bile acid metabolism in host tissue compartments , 2010, Proceedings of the National Academy of Sciences.

[27]  M. Uhlén,et al.  Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease , 2014, Nature Communications.

[28]  Lisa M. D'Souza,et al.  Genome sequence of the Brown Norway rat yields insights into mammalian evolution , 2004, Nature.

[29]  Tom M. Conrad,et al.  Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models , 2010, Molecular systems biology.

[30]  C. Gille,et al.  HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology , 2010, Molecular systems biology.

[31]  Riccardo Percudani,et al.  Completing the uric acid degradation pathway through phylogenetic comparison of whole genomes , 2006, Nature chemical biology.

[32]  Erhan Bilal,et al.  Understanding the limits of animal models as predictors of human biology: lessons learned from the sbv IMPROVER Species Translation Challenge , 2014, Bioinform..

[33]  Hyon K. Choi,et al.  Coffee, tea, and caffeine consumption and serum uric acid level: the third national health and nutrition examination survey. , 2007, Arthritis and rheumatism.

[34]  J. Reifman,et al.  Prediction of Metabolic Flux Distribution from Gene Expression Data Based on the Flux Minimization Principle , 2014, PloS one.

[35]  H. Yamada,et al.  The Japanese toxicogenomics project: application of toxicogenomics. , 2010, Molecular nutrition & food research.

[36]  D. Muzny,et al.  Two independent mutational events in the loss of urate oxidase during hominoid evolution , 2004, Journal of Molecular Evolution.

[37]  K. Brouwer,et al.  Endogenous bile acid disposition in rat and human sandwich-cultured hepatocytes. , 2012, Toxicology and applied pharmacology.

[38]  Y. Wang,et al.  Transplantation of immortalized human fetal hepatocytes prevents acute liver failure in 90% hepatectomized mice. , 2010, Transplantation proceedings.

[39]  Christopher J. Gregg,et al.  A red meat-derived glycan promotes inflammation and cancer progression , 2014, Proceedings of the National Academy of Sciences.

[40]  Svetoslav H. Slavov,et al.  Identification of a metabolic biomarker panel in rats for prediction of acute and idiosyncratic hepatotoxicity , 2014, Computational and structural biotechnology journal.

[41]  Bernhard O. Palsson,et al.  A detailed genome-wide reconstruction of mouse metabolism based on human Recon 1 , 2010, BMC Systems Biology.

[42]  Kyongbum Lee,et al.  Metabolic flux analysis of cultured hepatocytes exposed to plasma , 2003, Biotechnology and bioengineering.

[43]  Jason A. Papin,et al.  MetDraw: automated visualization of genome-scale metabolic network reconstructions and high-throughput data , 2014, Bioinform..

[44]  Jens Nielsen,et al.  The gut microbiota modulates host amino acid and glutathione metabolism in mice , 2015 .

[45]  W. Mattes,et al.  Detection of hepatotoxicity potential with metabolite profiling (metabolomics) of rat plasma. , 2014, Toxicology letters.

[46]  E. Nickerson,et al.  A mutation in human CMP-sialic acid hydroxylase occurred after the Homo-Pan divergence. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[47]  Ronan M. T. Fleming,et al.  A community-driven global reconstruction of human metabolism , 2013, Nature Biotechnology.

[48]  P. Roos,et al.  A novel CYP3 gene from female rats. , 1995, Biochimica et biophysica acta.

[49]  J. Lindon,et al.  Pharmaco-metabonomic phenotyping and personalized drug treatment , 2006, Nature.

[50]  D. Waxman,et al.  Detoxification of lithocholic acid. Elucidation of the pathways of oxidative metabolism in rat liver microsomes. , 1989, Journal of lipid research.

[51]  H. T. Parsons THE ANTISCORBUTIC CONTENT OF CERTAIN BODY TISSUES OF THE RAT THE PERSISTENCE OF THE ANTISCORBUTIC SUBSTANCE IN THE LIVER OF THE RAT AFTER LONG INTERVALS ON A SCORBUTIC DIET , 1920 .

[52]  M. Costache,et al.  Evolution of Fucosyltransferase Genes in Vertebrates* , 1997, The Journal of Biological Chemistry.

[53]  D. Searls Pharmacophylogenomics: genes, evolution and drug targets , 2003, Nature Reviews Drug Discovery.

[54]  M. Portolés,et al.  Effect of caffeine on urea biosynthesis and some related processes, ketone bodies, ATP and liver amino acids. , 1989, Biochemical pharmacology.

[55]  Natapol Pornputtapong,et al.  Reconstruction of Genome-Scale Active Metabolic Networks for 69 Human Cell Types and 16 Cancer Types Using INIT , 2012, PLoS Comput. Biol..

[56]  M. Whitehouse,et al.  Pharmacology of uricosuric drugs. , 1974, Annals of the rheumatic diseases.

[57]  Karthik Raman,et al.  Critical assessment of genome-scale metabolic networks: the need for a unified standard , 2015, Briefings Bioinform..

[58]  N. Kamatani,et al.  Theophylline increases serum uric acid levels. , 1984, The Journal of allergy and clinical immunology.

[59]  Biologie Rat Genome Database , 2011 .

[60]  J. Aubrecht,et al.  Evaluation of serum bile acid profiles as biomarkers of liver injury in rodents. , 2014, Toxicological sciences : an official journal of the Society of Toxicology.

[61]  B. Palsson,et al.  Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110 , 1994, Applied and environmental microbiology.

[62]  G. Stephanopoulos,et al.  Profiling of dynamic changes in hypermetabolic livers , 2003, Biotechnology and bioengineering.

[63]  Weisong Liu,et al.  The Rat Genome Database 2015: genomic, phenotypic and environmental variations and disease , 2014, Nucleic Acids Res..

[64]  J. Klinenberg,et al.  Effect of drugs on urate binding to plasma proteins , 1969, British medical journal.

[65]  I. Nookaew,et al.  Integration of clinical data with a genome-scale metabolic model of the human adipocyte , 2013, Molecular systems biology.

[66]  María Martín,et al.  UniProt: A hub for protein information , 2015 .

[67]  The Uniprot Consortium,et al.  UniProt: a hub for protein information , 2014, Nucleic Acids Res..

[68]  A. Ledbetter,et al.  Pulmonary Oxidative Stress, Inflammation, and Dysregulated Iron Homeostasis in Rat Models of Cardiovascular Disease , 2010, Journal of toxicology and environmental health. Part A.

[69]  Jason A. Papin,et al.  Reconciliation of Genome-Scale Metabolic Reconstructions for Comparative Systems Analysis , 2011, PLoS Comput. Biol..

[70]  M. T. Donato,et al.  Targeted profiling of circulating and hepatic bile acids in human, mouse, and rat using a UPLC-MRM-MS-validated method[S] , 2012, Journal of Lipid Research.