An improved empirical bayes approach to estimating differential gene expression in microarray time-course data: BETR (Bayesian Estimation of Temporal Regulation)
1 Department of Biostatistics, Harvard School of Public Health, 655 Huntington Avenue, Boston, Massachusetts 02115, USA
2 Department of Immunology and Infectious Diseases, Harvard School of Public Health, 655 Huntington Avenue, Boston, Massachusetts 02115, USA
3 Department of Statistics and Probability, Michigan State University, East Lansing, Massachusetts 48824, USA
4 Department of Biostatistics and Computational Biology and Department of Cancer Biology, Dana-Farber Cancer Institute, 44 Binney St, Boston, Massachusetts 02115, USA
BMC Bioinformatics 2009, 10:409 doi:10.1186/1471-2105-10-409Published: 10 December 2009
Microarray gene expression time-course experiments provide the opportunity to observe the evolution of transcriptional programs that cells use to respond to internal and external stimuli. Most commonly used methods for identifying differentially expressed genes treat each time point as independent and ignore important correlations, including those within samples and between sampling times. Therefore they do not make full use of the information intrinsic to the data, leading to a loss of power.
We present a flexible random-effects model that takes such correlations into account, improving our ability to detect genes that have sustained differential expression over more than one time point. By modeling the joint distribution of the samples that have been profiled across all time points, we gain sensitivity compared to a marginal analysis that examines each time point in isolation. We assign each gene a probability of differential expression using an empirical Bayes approach that reduces the effective number of parameters to be estimated.
Based on results from theory, simulated data, and application to the genomic data presented here, we show that BETR has increased power to detect subtle differential expression in time-series data. The open-source R package betr is available through Bioconductor. BETR has also been incorporated in the freely-available, open-source MeV software tool available from http://www.tm4.org/mev.html webcite.