Correction for hidden confounders in the genetic analysis of gene expression.
Understanding the genetic underpinnings of disease is important for screening, treatment, drug development, and basic biological insight. One way of getting at such an understanding is to find out which parts of our DNA, such as single-nucleotide polymorphisms, affect particular intermediary processes such as gene expression. Naively, such associations can be identified using a simple statistical test on all paired combinations of genetic variants and gene transcripts. However, a wide variety of confounders lie hidden in the data, leading to both spurious associations and missed associations if not properly addressed. We present a statistical model that jointly corrects for two particular kinds of hidden structure--population structure (e.g., race, family-relatedness), and microarray expression artifacts (e.g., batch effects), when these confounders are unknown. Applying our method to both real and synthetic, human and mouse data, we demonstrate the need for such a joint correction of confounders, and also the disadvantages of other possible approaches based on those in the current literature. In particular, we show that our class of models has maximum power to detect eQTL on synthetic data, and has the best performance on a bronze standard applied to real data. Lastly, our software and the associations we found with it are available at http://www.microsoft.com/science.
Listgarten
Jennifer
J
Microsoft Research, 1100 Glendon Avenue, Suite PH1, Los Angeles, CA, USA. jennl@microsoft.com
Kadie
Carl
C
Schadt
Eric E
EE
Heckerman
David
D
D000818:Animals; D030541:Databases, Genetic; D015870:Gene Expression; D055106:Genome-Wide Association Study / Q000706:statistics & numerical data; D023281:Genomics / Q000706:statistics & numerical data*; D006801:Humans; D016014:Linear Models; D051379:Mice; D008957:Models, Genetic; D015233:Models, Statistical; D020411:Oligonucleotide Array Sequence Analysis / Q000706:statistics & numerical data; D020641:Polymorphism, Single Nucleotide; D040641:Quantitative Trait Loci; D012984:Software
Nat Rev Genet. 2009 Mar;10(3):184-94
19223927
Nature. 2008 Mar 27;452(7186):429-35
18344982
PLoS Comput Biol. 2007 Apr 13;3(4):e69
17432931
Genet Epidemiol. 2008 Apr;32(3):215-26
18161052
Methods. 2009 Jul;48(3):265-76
19303049
Hum Mol Genet. 2008 Oct 15;17(R2):R129-34
18852201
PLoS Comput Biol. 2008 Nov;4(11):e1000225
19023406
Genetics. 2008 Dec;180(4):1909-25
18791227
Nat Genet. 2006 Aug;38(8):904-9
16862161
Genome Res. 2009 Nov;19(11):1942-52
19654370
PLoS Biol. 2008 May 6;6(5):e107
18462017
Mol Syst Biol. 2008;4:162
18319721
Nat Rev Genet. 2006 Oct;7(10):781-91
16983374
Trends Genet. 2008 Aug;24(8):408-15
18597885
Am J Hum Genet. 2000 Jul;67(1):170-81
10827107
PLoS Genet. 2006 Dec;2(12):e190
17194218
Genetics. 2007 May;176(1):675-83
17409088
Biometrics. 1999 Dec;55(4):997-1004
11315092
Genetics. 2008 Mar;178(3):1709-23
18385116
Mamm Genome. 2010 Apr;21(3-4):143-52
20135320
Proc Natl Acad Sci U S A. 2003 Aug 5;100(16):9440-5
12883005
PLoS Genet. 2007 Sep;3(9):1724-35
17907809
Nature. 2003 Mar 20;422(6929):297-302
12646919
Nat Genet. 2005 Jul;37(7):710-7
15965475
PLoS Genet. 2009 Jan;5(1):e1000358
19180192
Genet Epidemiol. 2001 Jan;20(1):4-16
11119293
Genet Epidemiol. 2010 Jan;34(1):51-9
19455578
PLoS Genet. 2007 Jan 19;3(1):e4
17238287
Nat Genet. 2006 Feb;38(2):203-8
16380716
Nat Rev Genet. 2009 Sep;10(9):595-604
19636342