Statistical methods for analyzing gene expression data for cancer research

A microarray analysis of the temporal response of liver to methylprednisolone: a comparative analysis of two dosing regimens

Microarray analyses were performed on livers from adrenalectomized male Wistar rats chronically infused with methylprednisolone (MPL) (0.3 mg/kg.h) using Alzet mini-osmotic pumps for periods ranging from 6 h to 7 d. Four control and 40 drug-treated animals were killed at 10 different times during drug infusion. Total RNA preparations from the livers of these animals were hybridized to 44 individual Affymetrix REA230A gene chips, generating data for 15,967 different probe sets for each chip. A series of three filters were applied sequentially. These filters were designed to eliminate probe sets that were not expressed in the tissue, were not regulated by the drug, or did not meet defined quality control standards. These filters eliminated 13,978 probe sets (87.5%) leaving a remainder of 1989 probe sets for further consideration. We previously described a similar dataset obtained from animals after administration of a single dose of MPL (50 mg/kg given iv). That study involved 16 time points over a 72-h period. A similar filtering schema applied to the single-bolus-dose dataset identified 1519 probe sets as being regulated by MPL. A comparison of datasets from the two different dosing regimens identified 358 genes that were regulated by MPL in response to both dosing regimens. Regulated genes were grouped into 13 categories, mainly on gene product function. The temporal profiles of these common genes were subjected to detailed scrutiny. Examination of temporal profiles demonstrates that current perspectives on the mechanism of glucocorticoid action cannot entirely explain the temporal profiles of these regulated genes.

Pharmacogenomic responses of rat liver to methylprednisolone: an approach to mining a rich microarray time series

A data set was generated to examine global changes in gene expression in rat liver over time in response to a single bolus dose of methylprednisolone. Four control animals and 43 drug-treated animals were humanely killed at 16 different time points following drug administration. Total RNA preparations from the livers of these animals were hybridized to 47 individual Affymetrix RU34A gene chips, generating data for 8799 different probe sets for each chip. Data mining techniques that are applicable to gene array time series data sets in order to identify drug-regulated changes in gene expression were applied to this data set. A series of 4 sequentially applied filters were developed that were designed to eliminate probe sets that were not expressed in the tissue, were not regulated by the drug treatment, or did not meet defined quality control standards. These filters eliminated 7287 probe sets of the 8799 total (82%) from further consideration. Application of judiciously chosen filters is an effective tool for data mining of time series data sets. The remaining data can then be further analyzed by clustering and mathematical modeling techniques.

The genomic response of skeletal muscle to methylprednisolone using microarrays: tailoring data mining to the structure of the pharmacogenomic time series

High-throughput data collection using gene microarrays has great potential as a method for addressing the pharmacogenomics of complex biological systems. Similarly, mechanism-based pharmacokinetic/pharmacodynamic modeling provides a tool for formulating quantitative testable hypotheses concerning the responses of complex biological systems. As the response of such systems to drugs generally entails cascades of molecular events in time, a time series design provides the best approach to capturing the full scope of drug effects. A major problem in using microarrays for high-throughput data collection is sorting through the massive amount of data in order to identify probe sets and genes of interest. Due to its inherent redundancy, a rich time series containing many time points and multiple samples per time point allows for the use of less stringent criteria of expression, expression change and data quality for initial filtering of unwanted probe sets. The remaining probe sets can then become the focus of more intense scrutiny by other methods, including temporal clustering, functional clustering and pharmacokinetic/pharmacodynamic modeling, which provide additional ways of identifying the probes and genes of pharmacological interest.