Bioinformatics Advance Access originally published online on October 25, 2005
Bioinformatics 2006 22(1):77-87; doi:10.1093/bioinformatics/bti737
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Integrating time-course microarray gene expression profiles with cytotoxicity for identification of biomarkers in primary rat hepatocytes exposed to cadmium

,*
1Cedars-Sinai Research Institute, Cedars-Sinai Medical Center Los Angeles, CA 90048, USA
2Center for Toxicoinformatics, National Center for Toxicological Research FDA, Jefferson, AR 72079, USA
3Applied Biotechnology Branch, US Air Force Research Laboratory WPAFB, OH 45434, USA
4David Geffen School of Medicine, UCLA Los Angeles, CA 90048, USA
*To whom correspondence should be addressed.
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 µM. 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).
Contact: cwang61{at}ucla.edu
Supplementary information: Supplementary data are available at Bioinformatics online.
Received on June 10, 2005; revised on September 28, 2005; accepted on October 20, 2005