title: Using control genes to correct for unwanted variation in microarray data.
abstract: Microarray expression studies suffer from the problem of batch effects and other unwanted variation. Many methods have been proposed to adjust microarray data to mitigate the problems of unwanted variation. Several of these methods rely on factor analysis to infer the unwanted variation from the data. A central problem with this approach is the difficulty in discerning the unwanted variation from the biological variation that is of interest to the researcher. We present a new method, intended for use in differential expression studies, that attempts to overcome this problem by restricting the factor analysis to negative control genes. Negative control genes are genes known a priori not to be differentially expressed with respect to the biological factor of interest. Variation in the expression levels of these genes can therefore be assumed to be unwanted variation. We name this method "Remove Unwanted Variation, 2-step" (RUV-2). We discuss various techniques for assessing the performance of an adjustment method and compare the performance of RUV-2 with that of other commonly used adjustment methods such as Combat and Surrogate Variable Analysis (SVA). We present several example studies, each concerning genes differentially expressed with respect to gender in the brain and find that RUV-2 performs as well or better than other methods. Finally, we discuss the possibility of adapting RUV-2 for use in studies not concerned with differential expression and conclude that there may be promise but substantial challenges remain.
authors.affiliation: Department of Statistics, University of California at Berkeley, Berkeley, CA 94720-3860, USA. johann@stat.berkeley.edu
authors.forename: Johann A
authors.initials: JA
authors.lastname: Gagnon-Bartsch
authors.forename: Terence P
authors.initials: TP
authors.lastname: Speed
mesh_terms: D003627:Data Interpretation, Statistical; D005260:Female; D020869:Gene Expression Profiling / Q000379:methods*; D006801:Humans; D008297:Male; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*
references.citation: Trends Genet. 2003 Jul;19(7):362-5
references.pmid: 12850439
references.citation: Biostatistics. 2007 Jan;8(1):118-27
references.pmid: 16632515
references.citation: Nucleic Acids Res. 2003 Feb 15;31(4):e15
references.pmid: 12582260
references.citation: Proc Natl Acad Sci U S A. 2010 Sep 21;107(38):16465-70
references.pmid: 20810919
references.citation: Nat Genet. 2010 Apr;42(4):348-54
references.pmid: 20208533
references.citation: Proc Natl Acad Sci U S A. 2000 Aug 29;97(18):10101-6
references.pmid: 10963673
references.citation: Genetics. 2008 Mar;178(3):1709-23
references.pmid: 18385116
references.citation: Genetics. 2008 Dec;180(4):1909-25
references.pmid: 18791227
references.citation: Nat Genet. 2006 Aug;38(8):904-9
references.pmid: 16862161
references.citation: PLoS Genet. 2007 Sep;3(9):1724-35
references.pmid: 17907809
references.citation: Stat Appl Genet Mol Biol. 2004;3:Article3
references.pmid: 16646809
references.citation: Bioinformatics. 2003 Jan 22;19(2):185-93
references.pmid: 12538238
references.citation: Nat Rev Genet. 2010 Oct;11(10):733-9
references.pmid: 20838408
references.citation: Neuropsychopharmacology. 2004 Feb;29(2):373-84
references.pmid: 14583743
references.citation: Bioinformatics. 2010 May 15;26(10):1308-15
references.pmid: 20363728
references.citation: BMC Res Notes. 2010 Dec 28;3:349
references.pmid: 21189145
references.citation: Biostatistics. 2003 Apr;4(2):249-64
references.pmid: 12925520
references.citation: Lancet. 2002 Apr 13;359(9314):1301-7
references.pmid: 11965276
references.citation: BMC Med Genomics. 2009 Aug 05;2:49
references.pmid: 19656400
references.citation: Nat Genet. 2006 Feb;38(2):203-8
references.pmid: 16380716
references.citation: Proc Natl Acad Sci U S A. 2008 Dec 2;105(48):18718-23
references.pmid: 19033188