Integrating time-course microarray gene expression profiles with cytotoxicity for identification of biomarkers in primary rat hepatocytes exposed to cadmium

MOTIVATION DNA microarrays can provide information about the expression levels of thousands of genes simultaneously at the transcriptomic level, while conventional cell viability and cytotoxicity measurement methods provide information about the biological functions at the cellular level. Integrating these data at different levels provides a promising approach for evaluating or predicting how cells respond to chemical exposure. It is important to investigate the multi-scale biological system in a systematic way to better understand the gene regulation networks and signal transduction pathways involved in the cellular responses to environmental factors. RESULTS Primary rat hepatocytes were exposed to cadmium acetate at 0, 1.25 and 2 microM. mRNA expression profiles at 0, 3, 6, 12 and 24 h were measured using the Affymetrix RatTox U34 GeneChip arrays. Simultaneously, cytotoxicity was assessed by lactase dehydrogenase leakage assay. Gene expression profiles at different time points were used to evaluate cytotoxicity at subsequent time points using partial least squares, and it was found that gene expression profiles at 0 h had the best prediction accuracy for the cytotoxicity observed at 12 h. Some biomarkers whose expression profiles showed strong relationship with cytotoxicity were identified and the underlying pathways were reconstructed to illustrate how hepatocytes respond to cadmium exposure. Permutation studies were also applied to assess the reliability of the predictive models. AVAILABILITY Matlab source code is available upon request and DNA microarray data are available at GEO (http://www.ncbi.nlm.nih.gov/geo).

[1]  M. Waalkes,et al.  Chronic toxic and carcinogenic effects of cadmium chloride in male DBA/2NCr and NFS/NCr mice: strain-dependent association with tumors of the hematopoietic system, injection site, liver, and lung. , 1994, Fundamental and applied toxicology : official journal of the Society of Toxicology.

[2]  D. Longo,et al.  Bcl-2 and Bcl-XL Block Thapsigargin-Induced Nitric Oxide Generation, c-Jun NH2-Terminal Kinase Activity, and Apoptosis , 1999, Molecular and Cellular Biology.

[3]  S. Sheweita,et al.  Drug-metabolizing enzymes: mechanisms and functions. , 2000, Current drug metabolism.

[4]  Charles Wang,et al.  Multi-class tumor classification by discriminant partial least squares using microarray gene expression data and assessment of classification models , 2004, Comput. Biol. Chem..

[5]  Timothy C. Elston,et al.  Complementary Roles for Toxicologic Pathology and Mathematics in Toxicogenomics, With Special Reference to Data Interpretation and Oscillatory Dynamics , 2004, Toxicologic pathology.

[6]  安部 哲哉.  Induction of heat shock 70 mRNA by cadmium is mediated by glutathione suppressive and non-suppressive triggers , 1995 .

[7]  L. Schiaffonati,et al.  Gene expression in liver after toxic injury: analysis of heat shock response and oxidative stress-inducible genes. , 2008, Liver.

[8]  John M Frazier,et al.  Involvement of apoptosis in hydrazine induced toxicity in rat primary hepatocytes. , 2003, Toxicology in vitro : an international journal published in association with BIBRA.

[9]  A. Kong,et al.  Pharmacogenomics, regulation and signaling pathways of phase I and II drug metabolizing enzymes. , 2002, Current drug metabolism.

[10]  M. Su,et al.  Cleavage of Poly(ADP-ribose) Polymerase by Interleukin-1β Converting Enzyme and Its Homologs TX and Nedd-2 (*) , 1995, The Journal of Biological Chemistry.

[11]  S. Orrenius,et al.  Isolation and use of liver cells. , 1978, Methods in enzymology.

[12]  V. Dixit,et al.  Apoptosis control by death and decoy receptors. , 1999, Current opinion in cell biology.

[13]  Danh V. Nguyen,et al.  Partial least squares proportional hazard regression for application to DNA microarray survival data , 2002, Bioinform..

[14]  B. Trottier,et al.  Studies on lipid peroxidation in rat tissues following administration of low and moderate doses of cadmium chloride. , 1991, Toxicology.

[15]  Xin Zhou,et al.  LS Bound based gene selection for DNA microarray data , 2005, Bioinform..

[16]  G. Salvesen,et al.  FLICE Induced Apoptosis in a Cell-free System , 1997, The Journal of Biological Chemistry.

[17]  S. Wold,et al.  The Collinearity Problem in Linear Regression. The Partial Least Squares (PLS) Approach to Generalized Inverses , 1984 .

[18]  V. H. Liao,et al.  Cadmium-regulated genes from the nematode Caenorhabditis elegans. Identification and cloning of new cadmium-responsive genes by differential display. , 1998, The Journal of biological chemistry.

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

[20]  Danh V. Nguyen,et al.  Tumor classification by partial least squares using microarray gene expression data , 2002, Bioinform..

[21]  John Lehmann.  Toxicogenomics takes center stage. , 2003, Drug news & perspectives.

[22]  R. Manne.  Analysis of two partial-least-squares algorithms for multivariate calibration , 1987 .

[23]  Danh V. Nguyen,et al.  Multi-class cancer classification via partial least squares with gene expression profiles , 2002, Bioinform..

[24]  J. Kyriakis,et al.  Direct Activation of Mitogen-Activated Protein Kinase Kinase Kinase MEKK1 by the Ste20p Homologue GCK and the Adapter Protein TRAF2 , 2002, Molecular and Cellular Biology.

[25]  H. Lodish,et al.  Fetal Anemia and Apoptosis of Red Cell Progenitors in Stat5a−/−5b−/− Mice A Direct Role for Stat5 in Bcl-XL Induction , 1999, Cell.

[26]  J. Prehn,et al.  Activation of Nuclear Factor κb and bcl-x Survival Gene Expression by Nerve Growth Factor Requires Tyrosine Phosphorylation of IκBα , 2001, The Journal of cell biology.

[27]  Leroy Hood,et al.  The impact of systems approaches on biological problems in drug discovery , 2004, Nature Biotechnology.

[28]  James T. Mac Grego.  Toxicology and Genetic Toxicology in the New Era of "Toxicogenomics": Impact of "-omics" Technologies , 2003 .

[29]  Philip Lazarus,et al.  Glucuronidation: an important mechanism for detoxification of benzo[a]pyrene metabolites in aerodigestive tract tissues. , 2002, Drug metabolism and disposition: the biological fate of chemicals.

[30]  Sten Orrenius,et al.  [4] Isolation and use of liver cells , 1978 .

[31]  Michael E Phelps,et al.  Systems Biology and New Technologies Enable Predictive and Preventative Medicine , 2004, Science.

[32]  R. Goyer,et al.  Toxicology of Metals , 1995, Handbook of Experimental Pharmacology.

[33]  G L Johnson,et al.  Differential activation of ERK and JNK mitogen-activated protein kinases by Raf-1 and MEKK. , 1994, Science.

[34]  R. Gentleman,et al.  Classification Using Generalized Partial Least Squares , 2005 .

[35]  Inyoung Kim,et al.  Statistical methods of translating microarray data into clinically relevant diagnostic information in colorectal cancer , 2005, Bioinform..

[36]  D. Townsend,et al.  The role of glutathione-S-transferase in anti-cancer drug resistance , 2003, Oncogene.

[37]  Guido Kroemer,et al.  Hsp27 negatively regulates cell death by interacting with cytochrome c , 2000, Nature Cell Biology.

[38]  R D Walter,et al.  Identification of stress-responsive genes in Caenorhabditis elegans using RT-PCR differential display. , 1998, Nucleic acids research.

[39]  J. Goldstein,et al.  Advances in mechanisms of activation and deactivation of environmental chemicals. , 1993, Environmental health perspectives.

[40]  Richard S Paules,et al.  Application of Toxicogenomics to Toxicology: Basic Concepts in the Analysis of Microarray Data , 2004, Toxicologic pathology.

[41]  P. Goering,et al.  Mercury, cadmium, and arsenite enhance heat shock protein synthesis in chick embryos prior to embryotoxicity. , 2003, Birth defects research. Part B, Developmental and reproductive toxicology.

[42]  Yumiko Sakai,et al.  Applications of microarrays with toxicologically relevant genes (tox genes) for the evaluation of chemical toxicants in Sprague Dawley rats in vivo and human hepatocytes in vitro. , 2004, Mutation research.

[43]  J. Hinson,et al.  Phase II enzymes and bioactivation. , 1995, Canadian journal of physiology and pharmacology.

[44]  C J Omiecinski,et al.  Epoxide hydrolases: biochemistry and molecular biology. , 2000, Chemico-biological interactions.

[45]  Beryllium, cadmium, mercury, and exposures in the glass manufacturing industry. Working Group views and expert opinions, Lyon, 9-16 February 1993. , 1993, IARC monographs on the evaluation of carcinogenic risks to humans.

[46]  G. Nemerow,et al.  MEK kinase 1 is critically required for c-Jun N-terminal kinase activation by proinflammatory stimuli and growth factor-induced cell migration. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[47]  N. DelRaso,et al.  EFFECT OF CULTURE CONDITIONS PRIOR TO EXPOSURE ON CADMIUM CYTOTOXICITY IN PRIMARY RAT HEPATOCYTES , 1999 .

[48]  Daniel L Farkas,et al.  Transcriptomic fingerprinting of bone marrow-derived hepatic beta2m-/Thy-1+ stem cells. , 2005, Biochemical and biophysical research communications.

[49]  D. Rubinsztein,et al.  Heat shock protein 27 prevents cellular polyglutamine toxicity and suppresses the increase of reactive oxygen species caused by huntingtin. , 2002, Human molecular genetics.

[50]  Richard S Paules,et al.  An overview of toxicogenomics. , 2002, Current issues in molecular biology.

[51]  N. Holbrook,et al.  Cellular response to oxidative stress: Signaling for suicide and survival * , 2002, Journal of cellular physiology.

[52]  J. Woodgett,et al.  The stress-activated protein kinase pathways , 1999, Cellular and Molecular Life Sciences CMLS.

[53]  D. Bleyl,et al.  IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. Overall Evaluations of Carcinogenicity: An Updating of IARC Monographs vol. 1 to 42. Supplement 7. 440 Seiten. International Agency for Research on Cancer, Lyon 1987. Preis: 65, – s.Fr , 1989 .

[54]  W. Kolch.  Meaningful relationships: the regulation of the Ras/Raf/MEK/ERK pathway by protein interactions. , 2000, The Biochemical journal.

[55]  R. Tibshirani,et al.  Diagnosis of multiple cancer types by shrunken centroids of gene expression , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[56]  L. Babiss,et al.  Toxicogenomics in predictive toxicology in drug development. , 2004, Chemistry & biology.

[57]  Gersende Fort,et al.  Classification using partial least squares with penalized logistic regression , 2005, Bioinform..

[58]  M. Tenenhaus,et al.  Prediction of clinical outcome with microarray data: a partial least squares discriminant analysis (PLS-DA) approach , 2003, Human Genetics.

[59]  A. Höskuldsson.  PLS regression methods , 1988 .

[60]  J. Mesirov,et al.  Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. , 1999, Science.

[61]  Weida Tong,et al.  Multi-class cancer classification by total principal component regression (TPCR) using microarray gene expression data , 2005, Nucleic acids research.

[62]  J Remacle,et al.  DNA microarrays as a tool in toxicogenomics. , 2004, Combinatorial chemistry & high throughput screening.

[63]  R Simon,et al.  Diagnostic and prognostic prediction using gene expression profiles in high-dimensional microarray data , 2003, British Journal of Cancer.

[64]  M. Waalkes.  Cadmium carcinogenesis in review. , 2000, Journal of inorganic biochemistry.

[65]  Geoffrey J McLachlan,et al.  Selection bias in gene extraction on the basis of microarray gene-expression data , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[66]  Wei Pan,et al.  Linear regression and two-class classification with gene expression data , 2003, Bioinform..

[67]  M. Radmacher,et al.  Pitfalls in the use of DNA microarray data for diagnostic and prognostic classification. , 2003, Journal of the National Cancer Institute.

[68]  Jacques Robert,et al.  Molecular Determinants of the Cytotoxicity of Platinum Compounds , 2004, Cancer Research.

[69]  Carlo Di Bello,et al.  PCA disjoint models for multiclass cancer analysis using gene expression data , 2003, Bioinform..

[70]  Brent D. Foy,et al.  Tailored gene array databases: applications in mechanistic toxicology , 2004, Bioinform..

[71]  J. Mesirov,et al.  Chemosensitivity prediction by transcriptional profiling , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[72]  A. Cederbaum,et al.  Requirement for iron for the production of hydroxyl radicals by rat liver quinone reductase. , 1993, The Journal of pharmacology and experimental therapeutics.

[73]  V. Dixit,et al.  Death receptors: signaling and modulation. , 1998, Science.

[74]  Albert J Fornace,et al.  Gadd45a protects against UV irradiation-induced skin tumors, and promotes apoptosis and stress signaling via MAPK and p53. , 2002, Cancer research.

[75]  Hiroshi Yamada,et al.  Investigation of a hepatotoxicity screening system in primary cell cultures --"what biomarkers would need to be addressed to estimate toxicity in conventional and new approaches?". , 2005, The Journal of toxicological sciences.

[76]  Emad S. Alnemri,et al.  Activation of a CrmA-insensitive, p35-sensitive Pathway in Ionizing Radiation-induced Apoptosis* , 1997, The Journal of Biological Chemistry.

[77]  S. Wold,et al.  PLS-regression: a basic tool of chemometrics , 2001 .

[78]  Lu Tian,et al.  Linking gene expression data with patient survival times using partial least squares , 2002, ISMB.

[79]  R. Barouki,et al.  Repression of gene expression by oxidative stress. , 1999, The Biochemical journal.

[80]  George Stephanopoulos,et al.  Mapping physiological states from microarray expression measurements , 2002, Bioinform..

[81]  V. H. Liao,et al.  Molecular Characterization of a Novel, Cadmium-inducible Gene from the Nematode Caenorhabditis elegans , 2002, The Journal of Biological Chemistry.

[82]  Tsuguyoshi Suzuki,et al.  Toxicology of metals , 1996 .

[83]  H. Ichijo,et al.  From receptors to stress-activated MAP kinases , 1999, Oncogene.

[84]  N. Ahn,et al.  Signal transduction through MAP kinase cascades. , 1998, Advances in cancer research.

[85]  D. Templeton,et al.  Protective elevations of glutathione and metallothionein in cadmium-exposed mesangial cells. , 1993, Toxicology.

[86]  R. Young,et al.  Biomedical Discovery with DNA Arrays , 2000, Cell.

[87]  R. Tibshirani,et al.  Significance analysis of microarrays applied to the ionizing radiation response , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[88]  Albert-László Barabási,et al.  Life's Complexity Pyramid , 2002, Science.

[89]  Yuri Lazebnik,et al.  Multiple species of CPP32 and Mch2 are the major active caspases present in apoptotic cells , 1997, The EMBO journal.

[90]  A. Boone,et al.  Constitutive heat shock protein 70 (HSC70) expression in rainbow trout hepatocytes: effect of heat shock and heavy metal exposure. , 2002, Comparative biochemistry and physiology. Toxicology & pharmacology : CBP.

[91]  F. Guéraud,et al.  Glucuronidation: a dual control. , 1998, General pharmacology.

[92]  M. Mumtaz,et al.  Development of a priority list of chemical mixtures occurring at 1188 hazardous waste sites, using the HazDat database. , 1996, Food and chemical toxicology : an international journal published for the British Industrial Biological Research Association.

[93]  Richard Baumgartner,et al.  Class prediction and discovery using gene microarray and proteomics mass spectroscopy data: curses, caveats, cautions , 2003, Bioinform..

[94]  G. Natoli,et al.  Activation of SAPK/JNK by TNF Receptor 1 Through a Noncytotoxic TRAF2-Dependent Pathway , 1997, Science.

[95]  Xuefeng Bruce Ling,et al.  Multiclass cancer classification and biomarker discovery using GA-based algorithms , 2005, Bioinform..

[96]  Shigeo Ohno,et al.  Tumor Necrosis Factor Signaling to Stress-activated Protein Kinase ( SAPK ) / Jun , 1998 .

[97]  R. Weichselbaum,et al.  Proteolytic activation of protein kinase C delta by an ICE‐like protease in apoptotic cells. , 1995, The EMBO journal.