Nucleic Acids Research Advance Access published online on October 13, 2009
Nucleic Acids Research, doi:10.1093/nar/gkp813
Genomics |
Local false discovery rate facilitates comparison of different microarray experiments
1Department of Medicine, 2Department of Biochemistry, 3Department of Statistics and 4Health Research and Policy, Stanford University Medical Center, Stanford, CA 94305, USA
*To whom correspondence should be addressed. Tel: +1 650 725 6442; Fax: +1 650 736 2282; Email: chu{at}stanford.edu
Received August 15, 2009. Accepted September 14, 2009.
The local false discovery rate (LFDR) estimates the probability of falsely identifying specific genes with changes in expression. In computer simulations, LFDR <10% successfully identified genes with changes in expression, while LFDR >90% identified genes without changes. We used LFDR to compare different microarray experiments quantitatively: (i) Venn diagrams of genes with and without changes in expression, (ii) scatter plots of the genes, (iii) correlation coefficients in the scatter plots and (iv) distributions of gene function. To illustrate, we compared three methods for pre-processing microarray data. Correlations between methods were high (r = 0.84–0.92). However, responses were often different in magnitude, and sometimes discordant, even though the methods used the same raw data. LFDR complements functional assessments like gene set enrichment analysis. To illustrate, we compared responses to ultraviolet radiation (UV), ionizing radiation (IR) and tobacco smoke. Compared to unresponsive genes, genes responsive to both UV and IR were enriched for cell cycle, mitosis, and DNA repair functions. Genes responsive to UV but not IR were depleted for cell adhesion functions. Genes responsive to tobacco smoke were enriched for detoxification functions. Thus, LFDR reveals differences and similarities among experiments.