Empirical bayes gene screening tool for time-course or dose-response microarray data

Abstract

An efficient method to reduce the dimensionality of microarray gene expression data from thousands or tens of thousands of cDNA clones down to a subset of the most differentially expressed cDNA clones is essential in order to simplify the massive amount of data generated from microarray experiments. An extension to the methods of Efron et al. [Efron, B., Tibshirani, R., Storey, J., Tusher, V. (2001). Empirical Bayes analysis of a microarray experiment. J. Am. Statist. Assoc. 96:1151-1160] is applied to a differential time-course experiment to determine a subset of cDNAs that have the largest probability of being differentially expressed with respect to treatment conditions across a set of unequally spaced time points. The proposed extension, which is advocated to be a screening tool, allows for inference across a continuous variable in addition to incorporating a more complex experimental design and allowing for multiple design replications. With the current data the focus is on a time-course experiment; however, the proposed methods can easily be implemented on a dose-response experiment, or any other microarray experiment that contains a continuous variable of interest. The proposed empirical Bayes gene-screening tool is compared with the Efron et al. (2001) method in addition to an adjusted model-based t-value using a time-course data set where the toxicological effect of a specific mixture of chemicals is being studied.