Proteome- and transcriptome-driven reconstruction of the human myocyte metabolic network and its use for identification of markers for diabetes.

Skeletal myocytes are metabolically active and susceptible to insulin resistance and are thus implicated in type 2 diabetes (T2D). This complex disease involves systemic metabolic changes, and their elucidation at the systems level requires genome-wide data and biological networks. Genome-scale metabolic models (GEMs) provide a network context for the integration of high-throughput data. We generated myocyte-specific RNA-sequencing data and investigated their correlation with proteome data. These data were then used to reconstruct a comprehensive myocyte GEM. Next, we performed a meta-analysis of six studies comparing muscle transcription in T2D versus healthy subjects. Transcriptional changes were mapped on the myocyte GEM, revealing extensive transcriptional regulation in T2D, particularly around pyruvate oxidation, branched-chain amino acid catabolism, and tetrahydrofolate metabolism, connected through the downregulated dihydrolipoamide dehydrogenase. Strikingly, the gene signature underlying this metabolic regulation successfully classifies the disease state of individual samples, suggesting that regulation of these pathways is a ubiquitous feature of myocytes in response to T2D.

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

[2]  Aarash Bordbar,et al.  A multi-tissue type genome-scale metabolic network for analysis of whole-body systems physiology , 2011, BMC Systems Biology.

[3]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[4]  P. Zimmet,et al.  The worldwide epidemiology of type 2 diabetes mellitus—present and future perspectives , 2012, Nature Reviews Endocrinology.

[5]  B. Pedersen,et al.  Deficient leukemia inhibitory factor signaling in muscle precursor cells from patients with type 2 diabetes. , 2012, American journal of physiology. Endocrinology and metabolism.

[6]  Dorothy D. Sears,et al.  Mechanisms of human insulin resistance and thiazolidinedione-mediated insulin sensitization , 2009, Proceedings of the National Academy of Sciences.

[7]  Jussi Paananen,et al.  Hyperglycemia and a Common Variant of GCKR Are Associated With the Levels of Eight Amino Acids in 9,369 Finnish Men , 2012, Diabetes.

[8]  M. Patti,et al.  Increased SRF transcriptional activity in human and mouse skeletal muscle is a signature of insulin resistance. , 2011, The Journal of clinical investigation.

[9]  S. Teichmann,et al.  RNA sequencing reveals two major classes of gene expression levels in metazoan cells , 2011, Molecular systems biology.

[10]  Cheng-chin Hsu,et al.  Histidine and carnosine delay diabetic deterioration in mice and protect human low density lipoprotein against oxidation and glycation. , 2005, European journal of pharmacology.

[11]  E. Lundberg,et al.  Towards a knowledge-based Human Protein Atlas , 2010, Nature Biotechnology.

[12]  E. Oakeley,et al.  Transcriptional profiling of myotubes from patients with type 2 diabetes: no evidence for a primary defect in oxidative phosphorylation genes , 2008, Diabetologia.

[13]  Sangsoo Kim,et al.  Combining multiple microarray studies and modeling interstudy variation , 2003, ISMB.

[14]  J. Bergström.  Percutaneous Needle Biopsy of Skeletal Muscle in Physiological and Clinical Research , 1975 .

[15]  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.

[16]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[17]  Douglas G Altman,et al.  Key Issues in Conducting a Meta-Analysis of Gene Expression Microarray Datasets , 2008, PLoS medicine.

[18]  N. M. van den Broek,et al.  Early or advanced stage type 2 diabetes is not accompanied by in vivo skeletal muscle mitochondrial dysfunction. , 2008, European journal of endocrinology.

[19]  P. Poulsen,et al.  Body Composition Is the Main Determinant for the Difference in Type 2 Diabetes Pathophysiology Between Japanese and Caucasians , 2014, Diabetes Care.

[20]  Jens Nielsen,et al.  Kiwi: a tool for integration and visualization of network topology and gene-set analysis , 2014, BMC Bioinformatics.

[21]  Zachary A. King,et al.  Constraint-based models predict metabolic and associated cellular functions , 2014, Nature Reviews Genetics.

[22]  J. Nielsen,et al.  Uncovering transcriptional regulation of metabolism by using metabolic network topology. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[23]  J. Sowers,et al.  The metabolic syndrome: role of skeletal muscle metabolism. , 2006, Annals of medicine.

[24]  M. Mann,et al.  Defining the transcriptome and proteome in three functionally different human cell lines , 2010, Molecular systems biology.

[25]  M. Roden,et al.  The role of mitochondria in insulin resistance and type 2 diabetes mellitus , 2012, Nature Reviews Endocrinology.

[26]  J. Nielsen,et al.  Identification of anticancer drugs for hepatocellular carcinoma through personalized genome‐scale metabolic modeling , 2014, Molecular systems biology.

[27]  A. Butte,et al.  Coordinated reduction of genes of oxidative metabolism in humans with insulin resistance and diabetes: Potential role of PGC1 and NRF1 , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[28]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[29]  M. Szyf,et al.  Gestational Diabetes Alters Offspring DNA Methylation Profiles in Human and Rat: Identification of Key Pathways Involved in Endocrine System Disorders, Insulin Signaling, Diabetes Signaling, and ILK Signaling. , 2015, Endocrinology.

[30]  M. Laakso,et al.  Downregulation of Diacylglycerol Kinase Delta Contributes to Hyperglycemia-Induced Insulin Resistance , 2008, Cell.

[31]  Richard D. Smith,et al.  Brain insulin lowers circulating BCAA levels by inducing hepatic BCAA catabolism. , 2014, Cell metabolism.

[32]  R. DeFronzo,et al.  Pathogenesis of Insulin Resistance in Skeletal Muscle , 2010, Journal of biomedicine & biotechnology.

[33]  E. Marcotte,et al.  Insights into the regulation of protein abundance from proteomic and transcriptomic analyses , 2012, Nature Reviews Genetics.

[34]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[35]  Maria Keays,et al.  ArrayExpress update—trends in database growth and links to data analysis tools , 2012, Nucleic Acids Res..

[36]  L. Mandarino,et al.  Increased Reactive Oxygen Species Production and Lower Abundance of Complex I Subunits and Carnitine Palmitoyltransferase 1B Protein Despite Normal Mitochondrial Respiration in Insulin-Resistant Human Skeletal Muscle , 2010, Diabetes.

[37]  Jens Nielsen,et al.  Genome‐scale modeling of human metabolism – a systems biology approach , 2013, Biotechnology journal.

[38]  B. Pedersen,et al.  Satellite Cells Derived from Obese Humans with Type 2 Diabetes and Differentiated into Myocytes In Vitro Exhibit Abnormal Response to IL-6 , 2012, PloS one.

[39]  Wolfgang Wurst,et al.  Expression of the splicing factor gene SFRS10 is reduced in human obesity and contributes to enhanced lipogenesis. , 2011, Cell metabolism.

[40]  M. Febbraio,et al.  Muscles, exercise and obesity: skeletal muscle as a secretory organ , 2012, Nature Reviews Endocrinology.

[41]  Nathan D. Price,et al.  Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE , 2012, BMC Systems Biology.

[42]  Claes Wahlestedt,et al.  Integration of microRNA changes in vivo identifies novel molecular features of muscle insulin resistance in type 2 diabetes , 2010, Genome Medicine.

[43]  I. Nookaew,et al.  Enriching the gene set analysis of genome-wide data by incorporating directionality of gene expression and combining statistical hypotheses and methods , 2013, Nucleic acids research.

[44]  N. M. van den Broek,et al.  Physical activity is the key determinant of skeletal muscle mitochondrial function in type 2 diabetes. , 2012, The Journal of clinical endocrinology and metabolism.

[45]  B. Pedersen,et al.  Elevated NF-κB Activation Is Conserved in Human Myocytes Cultured From Obese Type 2 Diabetic Patients and Attenuated by AMP-Activated Protein Kinase , 2011, Diabetes.

[46]  G. Getz,et al.  Inferring tumour purity and stromal and immune cell admixture from expression data , 2013, Nature Communications.

[47]  Intawat Nookaew,et al.  Novel insights into obesity and diabetes through genome-scale metabolic modeling , 2013, Front. Physiol..

[48]  J. Nielsen,et al.  Analysis of the Human Tissue-specific Expression by Genome-wide Integration of Transcriptomics and Antibody-based Proteomics* , 2013, Molecular & Cellular Proteomics.

[49]  M. Selbach,et al.  Global quantification of mammalian gene expression control , 2011, Nature.

[50]  V. M. Zat︠s︡iorskiĭ,et al.  Biomechanics of Skeletal Muscles , 2012 .

[51]  Y. Kido,et al.  Histidine Augments the Suppression of Hepatic Glucose Production by Central Insulin Action , 2013, Diabetes.

[52]  S. Shoelson,et al.  Type 2 diabetes as an inflammatory disease , 2011, Nature Reviews Immunology.

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

[54]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[55]  Martin Kircher,et al.  Deep proteome and transcriptome mapping of a human cancer cell line , 2011, Molecular systems biology.

[56]  M. Daly,et al.  PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes , 2003, Nature Genetics.

[57]  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..

[58]  Rong Wang,et al.  Integrating shotgun proteomics and mRNA expression data to improve protein identification , 2009, Bioinform..

[59]  C. Lynch,et al.  Branched-chain amino acids in metabolic signalling and insulin resistance , 2014, Nature Reviews Endocrinology.