A Three Stage Integrative Pathway Search (TIPS©) framework to identify toxicity relevant genes and pathways

BackgroundThe ability to obtain profiles of gene expressions, proteins and metabolites with the advent of high throughput technologies has advanced the study of pathway and network reconstruction. Genome-wide network reconstruction requires either interaction measurements or large amount of perturbation data, often not available for mammalian cell systems. To overcome these shortcomings, we developed a Three Stage Integrative Pathway Search (TIPS©) approach to reconstruct context-specific active pathways involved in conferring a specific phenotype, from limited amount of perturbation data. The approach was tested on human liver cells to identify pathways that confer cytotoxicity.ResultsThis paper presents a systems approach that integrates gene expression and cytotoxicity profiles to identify a network of pathways involved in free fatty acid (FFA) and tumor necrosis factor-α (TNF-α) induced cytotoxicity in human hepatoblastoma cells (HepG2/C3A). Cytotoxicity relevant genes were first identified and then used to reconstruct a network using Bayesian network (BN) analysis. BN inference was used subsequently to predict the effects of perturbing a gene on the other genes in the network and on the cytotoxicity. These predictions were subsequently confirmed through the published literature and further experiments.ConclusionThe TIPS© approach is able to reconstruct active pathways that confer a particular phenotype by integrating gene expression and phenotypic profiles. A web-based version of TIPS© that performs the analysis described herein can be accessed at http://www.egr.msu.edu/tips.

[1]  Wolfram Liebermeister,et al.  Linear modes of gene expression determined by independent component analysis , 2002, Bioinform..

[2]  F. Nargang,et al.  The Tim8-Tim13 Complex of Neurospora crassa Functions in the Assembly of Proteins into Both Mitochondrial Membranes* , 2004, Journal of Biological Chemistry.

[3]  Christina Chan,et al.  Hydrogen peroxide and hydroxyl radicals mediate palmitate-induced cytotoxicity to hepatoma cells: Relation to mitochondrial permeability transition , 2007, Free radical research.

[4]  Takashi Kadowaki,et al.  [Molecular mechanism of insulin resistance]. , 2005, Nihon rinsho. Japanese journal of clinical medicine.

[5]  Tom Burr,et al.  Causation, Prediction, and Search , 2003, Technometrics.

[6]  J. Ntambi,et al.  Peroxisome proliferators induce mouse liver stearoyl-CoA desaturase 1 gene expression. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Nir Friedman,et al.  Inferring Cellular Networks Using Probabilistic Graphical Models , 2004, Science.

[8]  Michael A. Savageau,et al.  Application of Biochemical Systems Theory to Metabolism in Human Red Blood Cells , 1996, The Journal of Biological Chemistry.

[9]  Raymond T Chung,et al.  Protein kinase R is increased and is functional in hepatitis C virus–related hepatocellular carcinoma , 2003, American Journal of Gastroenterology.

[10]  Christina Chan,et al.  Integrating Gene Expression and Metabolic Profiles* , 2004, Journal of Biological Chemistry.

[11]  M. Savageau Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[12]  Max Henrion,et al.  Propagating uncertainty in bayesian networks by probabilistic logic sampling , 1986, UAI.

[13]  Zheng Li,et al.  Inferring pathways and networks with a Bayesian framework , 2004, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[14]  Gregory F. Cooper,et al.  The Computational Complexity of Probabilistic Inference Using Bayesian Belief Networks , 1990, Artif. Intell..

[15]  Akihiko Satoh,et al.  PKC-delta and -epsilon regulate NF-kappaB activation induced by cholecystokinin and TNF-alpha in pancreatic acinar cells. , 2004, American journal of physiology. Gastrointestinal and liver physiology.

[16]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[17]  M. Savageau Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[18]  Kathleen E. Sullivan,et al.  A role for PKC-δ and PI 3-kinase in TNF-α-mediated antiapoptotic signaling in the human neutrophil , 2002 .

[19]  Takahiro Doi,et al.  Tumor necrosis factor-alpha-induced IKK phosphorylation of NF-kappaB p65 on serine 536 is mediated through the TRAF2, TRAF5, and TAK1 signaling pathway. , 2003, The Journal of biological chemistry.

[20]  Joachim Rassow,et al.  Bcl-2 and porin follow different pathways of TOM-dependent insertion into the mitochondrial outer membrane. , 2002, Journal of molecular biology.

[21]  Dirk Husmeier,et al.  Sensitivity and specificity of inferring genetic regulatory interactions from microarray experiments with dynamic Bayesian networks , 2003, Bioinform..

[22]  J. Senior,et al.  Drug-related hepatotoxicity. , 2006, The New England journal of medicine.

[23]  T. Jaakkola,et al.  Validation and refinement of gene-regulatory pathways on a network of physical interactions , 2005, Genome Biology.

[24]  Chris Wiggins,et al.  ARACNE: An Algorithm for the Reconstruction of Gene Regulatory Networks in a Mammalian Cellular Context , 2004, BMC Bioinformatics.

[25]  David J. C. MacKay,et al.  A decomposition model to track gene expression signatures: preview on observer-independent classification of ovarian cancer , 2002, Bioinform..

[26]  J. Ntambi,et al.  The regulation of stearoyl-CoA desaturase (SCD). , 1995, Progress in lipid research.

[27]  Takahiro Doi,et al.  Tumor Necrosis Factor-α-induced IKK Phosphorylation of NF-κB p65 on Serine 536 Is Mediated through the TRAF2, TRAF5, and TAK1 Signaling Pathway* , 2003, Journal of Biological Chemistry.

[28]  M. Mathews,et al.  Human breast cancer cells contain elevated levels and activity of the protein kinase, PKR , 2000, Oncogene.

[29]  J. Haddad,et al.  On the antioxidant mechanisms of Bcl-2: a retrospective of NF-kappaB signaling and oxidative stress. , 2004, Biochemical and biophysical research communications.

[30]  M. Mumby,et al.  Reversible Phosphorylation of Bcl2 following Interleukin 3 or Bryostatin 1 Is Mediated by Direct Interaction with Protein Phosphatase 2A* , 1998, The Journal of Biological Chemistry.

[31]  Adam A. Margolin,et al.  Reverse engineering of regulatory networks in human B cells , 2005, Nature Genetics.

[32]  Akihiko Satoh,et al.  PKC-δ and -ε regulate NF-κB activation induced by cholecystokinin and TNF-α in pancreatic acinar cells , 2004 .

[33]  Peter Vandenabeele,et al.  Translation Inhibition in Apoptosis , 2001, The Journal of Biological Chemistry.

[34]  Dale Schuurmans,et al.  Learning Bayesian Nets that Perform Well , 1997, UAI.

[35]  M. Leroy,et al.  A role for PKCζ in the LPS-induced translocation NF-κB p65 subunit in cultured myometrial cells , 2005 .

[36]  J. Collins,et al.  Chemogenomic profiling on a genome-wide scale using reverse-engineered gene networks , 2005, Nature Biotechnology.

[37]  K. Sullivan,et al.  A role for PKC-delta and PI 3-kinase in TNF-alpha-mediated antiapoptotic signaling in the human neutrophil. , 2002, American journal of physiology. Cell physiology.

[38]  Haruo Okado,et al.  Tumor Necrosis Factor Induces Bcl-2 and Bcl-x Expression through NFκB Activation in Primary Hippocampal Neurons* , 1999, The Journal of Biological Chemistry.

[39]  Sib Sankar Roy,et al.  Molecular mechanism of insulin resistance , 2007, Journal of Biosciences.

[40]  Satoru Miyano,et al.  Dynamic Bayesian Network and Nonparametric Regression for Nonlinear Modeling of Gene Networks from Time Series Gene Expression Data , 2003, CMSB.

[41]  Nir Friedman,et al.  Bayesian Network Classifiers , 1997, Machine Learning.

[42]  C. Hoppel,et al.  Fatty acid import into mitochondria. , 2000, Biochimica et biophysica acta.

[43]  Byung Chul Kim,et al.  Tumor Necrosis Factor Induces Apoptosis in Hepatoma Cells by Increasing Ca2+ Release from the Endoplasmic Reticulum and Suppressing Bcl-2 Expression* , 2002, The Journal of Biological Chemistry.

[44]  David A. Bell,et al.  Learning Bayesian networks from data: An information-theory based approach , 2002, Artif. Intell..

[45]  Harini Sampath,et al.  Stearoyl-CoA Desaturase-1 Mediates the Pro-lipogenic Effects of Dietary Saturated Fat* , 2007, Journal of Biological Chemistry.

[46]  Daniel S. Ory,et al.  Palmitate-induced Apoptosis Can Occur through a Ceramide-independent Pathway* , 2001, The Journal of Biological Chemistry.

[47]  Jeremy K Nicholson,et al.  An integrated reverse functional genomic and metabolic approach to understanding orotic acid-induced fatty liver. , 2004, Physiological genomics.

[48]  Lawrence Carin,et al.  Using Functional Genomic Units to Corroborate User Experiments with the Rosetta Compendium , 2002 .

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

[50]  Jie-Jun Tseng,et al.  Reconstructing gene regulatory networks from time-series microarray data , 2005 .

[51]  P. Vandenabeele,et al.  Translation inhibition in apoptosis: caspase-dependent PKR activation and eIF2-alpha phosphorylation. , 2001, The Journal of biological chemistry.

[52]  J. Haddad,et al.  Oxygen-sensing mechanisms and the regulation of redox-responsive transcription factors in development and pathophysiology , 2002, Respiratory research.

[53]  Weiru Liu,et al.  Learning belief networks from data: an information theory based approach , 1997, CIKM '97.

[54]  A. Golay,et al.  Pathways from obesity to diabetes , 2002, International Journal of Obesity.

[55]  M. Holland,et al.  Transcript Abundance in Yeast Varies over Six Orders of Magnitude* , 2002, The Journal of Biological Chemistry.

[56]  Paul P. Wang,et al.  Advances to Bayesian network inference for generating causal networks from observational biological data , 2004, Bioinform..

[57]  Michal Linial,et al.  Using Bayesian Networks to Analyze Expression Data , 2000, J. Comput. Biol..

[58]  M. Leroy,et al.  A role for PKCzeta in the LPS-induced translocation NF-kappaB p65 subunit in cultured myometrial cells. , 2005, Biochimie.

[59]  C. Ciudad,et al.  Differential induction of stearoyl-CoA desaturase and acyl-CoA oxidase genes by fibrates in HepG2 cells. , 2001, Biochemical pharmacology.

[60]  P. Spirtes,et al.  Causation, prediction, and search , 1993 .

[61]  K. Sachs,et al.  Causal Protein-Signaling Networks Derived from Multiparameter Single-Cell Data , 2005, Science.

[62]  G. Kroemer,et al.  Palmitate induces apoptosis via a direct effect on mitochondria , 1999, Apoptosis.

[63]  M A Savageau,et al.  Model assessment and refinement using strategies from biochemical systems theory: application to metabolism in human red blood cells. , 1996, Journal of theoretical biology.

[64]  Zan Xu,et al.  The B56α Regulatory Subunit of Protein Phosphatase 2A Is a Target for Regulation by Double-Stranded RNA-Dependent Protein Kinase PKR , 2000, Molecular and Cellular Biology.

[65]  Byung-Chul Kim,et al.  Tumor necrosis factor induces apoptosis in hepatoma cells by increasing Ca(2+) release from the endoplasmic reticulum and suppressing Bcl-2 expression. , 2002, The Journal of biological chemistry.

[66]  Korbinian Strimmer,et al.  An empirical Bayes approach to inferring large-scale gene association networks , 2005, Bioinform..

[67]  Benno Schwikowski,et al.  Discovering regulatory and signalling circuits in molecular interaction networks , 2002, ISMB.

[68]  H. Tilg,et al.  Cytokines and liver diseases. , 2001, Canadian journal of gastroenterology = Journal canadien de gastroenterologie.

[69]  D. Heckerman,et al.  A Bayesian Approach to Causal Discovery , 2006 .

[70]  J. Friedman,et al.  Stearoyl-CoA desaturase 1 deficiency increases fatty acid oxidation by activating AMP-activated protein kinase in liver. , 2004, Proceedings of the National Academy of Sciences of the United States of America.