Assessing Concordance of Drug-Induced Transcriptional Response in Rodent Liver and Cultured Hepatocytes

The effect of drugs, disease and other perturbations on mRNA levels are studied using gene expression microarrays or RNA-seq, with the goal of understanding molecular effects arising from the perturbation. Previous comparisons of reproducibility across laboratories have been limited in scale and focused on a single model. The use of model systems, such as cultured primary cells or cancer cell lines, assumes that mechanistic insights derived from the models would have been observed via in vivo studies. We examined the concordance of compound-induced transcriptional changes using data from several sources: rat liver and rat primary hepatocytes (RPH) from Drug Matrix (DM) and open TG-GATEs (TG), human primary hepatocytes (HPH) from TG, and mouse liver / HepG2 results from the Gene Expression Omnibus (GEO) repository. Gene expression changes for treatments were normalized to controls and analyzed with three methods: 1) gene level for 9071 high expression genes in rat liver, 2) gene set analysis (GSA) using canonical pathways and gene ontology sets, 3) weighted gene co-expression network analysis (WGCNA). Co-expression networks performed better than genes or GSA when comparing treatment effects within rat liver and rat vs. mouse liver. Genes and modules performed similarly at Connectivity Map-style analyses, where success at identifying similar treatments among a collection of reference profiles is the goal. Comparisons between rat liver and RPH, and those between RPH, HPH and HepG2 cells reveal lower concordance for all methods. We observe that the baseline state of untreated cultured cells relative to untreated rat liver shows striking similarity with toxicant-exposed cells in vivo, indicating that gross systems level perturbation in the underlying networks in culture may contribute to the low concordance.

[1]  E R Weibel,et al.  Distribution of Organelles and Membranes between Hepatocytes and Nonhepatocytes in a Stereological Study , 2022 .

[2]  R. Myers,et al.  Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data , 2005, Nucleic acids research.

[3]  Stephen C. Harris,et al.  Rat toxicogenomic study reveals analytical consistency across microarray platforms , 2006, Nature Biotechnology.

[4]  Alexander A. Morgan,et al.  Discovery and Preclinical Validation of Drug Indications Using Compendia of Public Gene Expression Data , 2011, Science Translational Medicine.

[5]  J. Kleinjans,et al.  A transcriptomics-based in vitro assay for predicting chemical genotoxicity in vivo. , 2012, Carcinogenesis.

[6]  Jan G. Hengstler,et al.  Pharmacokinetics explain in vivo/in vitro discrepancies of carcinogen-induced gene expression alterations in rat liver and cultivated hepatocytes , 2012, Archives of Toxicology.

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

[8]  S. Horvath,et al.  A General Framework for Weighted Gene Co-Expression Network Analysis , 2005, Statistical applications in genetics and molecular biology.

[9]  J. Zhang,et al.  Data mining reveals a network of early-response genes as a consensus signature of drug-induced in vitro and in vivo toxicity , 2013, The Pharmacogenomics Journal.

[10]  A. Butte,et al.  Leveraging models of cell regulation and GWAS data in integrative network-based association studies , 2012, Nature Genetics.

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

[12]  T. Zacharewski,et al.  In vivo – in vitro toxicogenomic comparison of TCDD-elicited gene expression in Hepa1c1c7 mouse hepatoma cells and C57BL/6 hepatic tissue , 2006, BMC Genomics.

[13]  Yudong D. He,et al.  Systems analysis of eleven rodent disease models reveals an inflammatome signature and key drivers , 2012, Molecular systems biology.

[14]  D. Friend,et al.  HIGH-YIELD PREPARATION OF ISOLATED RAT LIVER PARENCHYMAL CELLS , 1969, The Journal of cell biology.

[15]  J. Fraczek,et al.  Primary hepatocyte cultures for pharmaco-toxicological studies: at the busy crossroad of various anti-dedifferentiation strategies , 2012, Archives of Toxicology.

[16]  J. Kleinjans,et al.  Evaluating microRNA profiles reveals discriminative responses following genotoxic or non-genotoxic carcinogen exposure in primary mouse hepatocytes. , 2015, Mutagenesis.

[17]  David P. Kreil,et al.  The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance , 2014, Nature Biotechnology.

[18]  Eric T. Wang,et al.  An Abundance of Ubiquitously Expressed Genes Revealed by Tissue Transcriptome Sequence Data , 2009, PLoS Comput. Biol..

[19]  Andrej Bugrim,et al.  Identification of interspecies concordance of mechanisms of arsenic-induced bladder cancer. , 2007, Toxicology in vitro : an international journal published in association with BIBRA.

[20]  J. Lo-Guidice,et al.  Gene Expression Profiling of Systems Involved in the Metabolism and the Disposition of Xenobiotics: Comparison between Human Intestinal Biopsy Samples and Colon Cell Lines , 2012, Drug Metabolism and Disposition.

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

[22]  Maqc Consortium.  The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements , 2006, Nature Biotechnology.

[23]  Trey Ideker,et al.  Boosting Signal-to-Noise in Complex Biology: Prior Knowledge Is Power , 2011, Cell.

[24]  Qianqian Zhu,et al.  Preferred analysis methods for Affymetrix GeneChips. II. An expanded, balanced, wholly-defined spike-in dataset , 2010, BMC Bioinformatics.

[25]  N. Shanks,et al.  Are animal models predictive for humans? , 2009, Philosophy, ethics, and humanities in medicine : PEHM.

[26]  David Z. Chen,et al.  Architecture of the human regulatory network derived from ENCODE data , 2012, Nature.

[27]  Liam O'Connor,et al.  Colorectal cancer cell lines are representative models of the main molecular subtypes of primary cancer. , 2014, Cancer research.

[28]  Apurva R. Patel,et al.  AlgiMatrix™ Based 3D Cell Culture System as an In-Vitro Tumor Model for Anticancer Studies , 2013, PloS one.

[29]  Paul A Clemons,et al.  The Connectivity Map: Using Gene-Expression Signatures to Connect Small Molecules, Genes, and Disease , 2006, Science.

[30]  K. Hillgren,et al.  Chimeric TK-NOG Mice: A Predictive Model for Cholestatic Human Liver Toxicity , 2015, The Journal of Pharmacology and Experimental Therapeutics.

[31]  Ravi Iyengar,et al.  Network analyses in systems pharmacology , 2009, Bioinform..

[32]  Rui Luo,et al.  Is My Network Module Preserved and Reproducible? , 2011, PLoS Comput. Biol..

[33]  Daniel R. Salomon,et al.  Strategies for aggregating gene expression data: The collapseRows R function , 2011, BMC Bioinformatics.

[34]  Jeroen L A Pennings,et al.  A transcriptomics-based hepatotoxicity comparison between the zebrafish embryo and established human and rodent in vitro and in vivo models using cyclosporine A, amiodarone and acetaminophen. , 2015, Toxicology letters.

[35]  V Vécsei,et al.  Dedifferentiation-associated changes in morphology and gene expression in primary human articular chondrocytes in cell culture. , 2002, Osteoarthritis and cartilage.

[36]  C. Austin,et al.  Improving the Human Hazard Characterization of Chemicals: A Tox21 Update , 2013, Environmental health perspectives.

[37]  Adam S. Hayward,et al.  Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME , 2013, Archives of Toxicology.

[38]  HighWire Press,et al.  Drug metabolism and disposition : the biological fate of chemicals. , 1973 .

[39]  David P. Kreil,et al.  A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control consortium , 2014, Nature Biotechnology.

[40]  J. Kleinjans,et al.  Comparison of genotoxicant-modified transcriptomic responses in conventional and epigenetically stabilized primary rat hepatocytes with in vivo rat liver data , 2012, Archives of Toxicology.

[41]  Silvio Albertini,et al.  Gene expression in two hepatic cell lines, cultured primary hepatocytes, and liver slices compared to the in vivo liver gene expression in rats: possible implications for toxicogenomics use of in vitro systems. , 2003, Toxicological sciences : an official journal of the Society of Toxicology.

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

[43]  Seon-Young Kim,et al.  PAGE: Parametric Analysis of Gene Set Enrichment , 2005, BMC Bioinform..

[44]  Benjamin Haibe-Kains,et al.  Characterization of Conserved Toxicogenomic Responses in Chemically Exposed Hepatocytes across Species and Platforms , 2015, Environmental health perspectives.

[45]  D. Bergstrom,et al.  Evaluating patient-derived colorectal cancer xenografts as preclinical models by comparison with patient clinical data. , 2015, Cancer research.

[46]  R. Schwabe,et al.  Gene expression profiles during hepatic stellate cell activation in culture and in vivo. , 2007, Gastroenterology.

[47]  Patrik D'haeseleer,et al.  How does gene expression clustering work? , 2005, Nature Biotechnology.

[48]  Robert Petryszak,et al.  ArrayExpress update—simplifying data submissions , 2014, Nucleic Acids Res..

[49]  J. Kleinjans,et al.  Gene expression profiling in primary mouse hepatocytes discriminates true from false-positive genotoxic compounds. , 2010, Mutagenesis.

[50]  M. Fielden,et al.  Development of a large-scale chemogenomics database to improve drug candidate selection and to understand mechanisms of chemical toxicity and action. , 2005, Journal of biotechnology.

[51]  Weida Tong,et al.  Of text and gene – using text mining methods to uncover hidden knowledge in toxicogenomics , 2014, BMC Systems Biology.

[52]  P. Seglen.  Preparation of rat liver cells. I. Effect of Ca 2+ on enzymatic dispersion of isolated, perfused liver. , 1972, Experimental cell research.

[53]  Atsushi Ono,et al.  Profiling of gene expression in rat liver and rat primary cultured hepatocytes treated with peroxisome proliferators. , 2006, The Journal of toxicological sciences.

[54]  Maarten L J Coonen,et al.  Classification of hepatotoxicants using HepG2 cells: A proof of principle study. , 2014, Chemical research in toxicology.

[55]  T. Hartung.  Toxicology for the twenty-first century , 2009, Nature.

[56]  Yudong D. He,et al.  Common Handling Procedures Conducted in Preclinical Safety Studies Result in Minimal Hepatic Gene Expression Changes in Sprague-Dawley Rats , 2014, PloS one.

[57]  T. Zacharewski,et al.  Differences in TCDD-elicited gene expression profiles in human HepG2, mouse Hepa1c1c7 and rat H4IIE hepatoma cells , 2011, BMC Genomics.

[58]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[59]  S. Horvath,et al.  Divergence of human and mouse brain transcriptome highlights Alzheimer disease pathways , 2010, Proceedings of the National Academy of Sciences.

[60]  J. Kleinjans,et al.  Benzo[a]pyrene-induced transcriptomic responses in primary hepatocytes and in vivo liver: Toxicokinetics is essential for in vivo–in vitro comparisons , 2012, Archives of Toxicology.

[61]  L. Tran,et al.  Integrated Systems Approach Identifies Genetic Nodes and Networks in Late-Onset Alzheimer’s Disease , 2013, Cell.

[62]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[63]  Julio Saez-Rodriguez,et al.  DvD: An R/Cytoscape pipeline for drug repurposing using public repositories of gene expression data , 2012, Bioinform..

[64]  Deepak K Rajpal,et al.  Applications of Connectivity Map in drug discovery and development. , 2012, Drug discovery today.

[65]  Yong Qian,et al.  Multi-walled carbon nanotube-induced gene expression in vitro: concordance with in vivo studies. , 2015, Toxicology.

[66]  J. Kelm,et al.  Multi-cell type human liver microtissues for hepatotoxicity testing , 2012, Archives of Toxicology.

[67]  Per Artursson,et al.  Rapid measurement of intracellular unbound drug concentrations. , 2013, Molecular pharmaceutics.