A Bayesian framework to account for complex non-genetic factors in gene expression levels greatly increases power in eQTL studies.
Gene expression measurements are influenced by a wide range of factors, such as the state of the cell, experimental conditions and variants in the sequence of regulatory regions. To understand the effect of a variable of interest, such as the genotype of a locus, it is important to account for variation that is due to confounding causes. Here, we present VBQTL, a probabilistic approach for mapping expression quantitative trait loci (eQTLs) that jointly models contributions from genotype as well as known and hidden confounding factors. VBQTL is implemented within an efficient and flexible inference framework, making it fast and tractable on large-scale problems. We compare the performance of VBQTL with alternative methods for dealing with confounding variability on eQTL mapping datasets from simulations, yeast, mouse, and human. Employing Bayesian complexity control and joint modelling is shown to result in more precise estimates of the contribution of different confounding factors resulting in additional associations to measured transcript levels compared to alternative approaches. We present a threefold larger collection of cis eQTLs than previously found in a whole-genome eQTL scan of an outbred human population. Altogether, 27% of the tested probes show a significant genetic association in cis, and we validate that the additional eQTLs are likely to be real by replicating them in different sets of individuals. Our method is the next step in the analysis of high-dimensional phenotype data, and its application has revealed insights into genetic regulation of gene expression by demonstrating more abundant cis-acting eQTLs in human than previously shown. Our software is freely available online at http://www.sanger.ac.uk/resources/software/peer/.
Stegle
Oliver
O
Max Planck Institutes Tübingen, Tübingen, Germany. oliver.stegle@tuebingen.mpg.de
Parts
Leopold
L
Durbin
Richard
R
Winn
John
J
D000818:Animals; D001499:Bayes Theorem; D030541:Databases, Genetic; D015870:Gene Expression; D006801:Humans; D020407:Internet; D008390:Markov Chains; D051379:Mice; D008957:Models, Genetic; D015233:Models, Statistical; D010641:Phenotype; D040641:Quantitative Trait Loci; D015203:Reproducibility of Results; D012984:Software; D015003:Yeasts
Nature. 2008 Mar 27;452(7186):423-8
18344981
Proc Natl Acad Sci U S A. 2005 Feb 1;102(5):1572-7
15659551
PLoS Genet. 2008 Dec;4(12):e1000294
19057673
BMC Bioinformatics. 2008 May 20;9:244
18492285
Nat Genet. 2007 Feb;39(2):226-31
17206142
Genome Res. 2007 Dec;17(12):1707-16
18063559
Nature. 2005 Oct 27;437(7063):1299-320
16255080
Nat Genet. 2007 Oct;39(10):1217-24
17873874
Science. 2007 Feb 9;315(5813):848-53
17289997
Genetics. 2008 Dec;180(4):1909-25
18791227
PLoS One. 2008 Aug 13;3(8):e2966
18698422
PLoS Genet. 2007 Sep;3(9):1724-35
17907809
Nat Genet. 2005 Jul;37(7):710-7
15965475
Nat Rev Genet. 2008 Aug;9(8):575-81
18574472
Science. 2002 Apr 26;296(5568):752-5
11923494
Genetics. 1989 Jan;121(1):185-99
2563713
Nucleic Acids Res. 2007 Jul;35(Web Server issue):W193-200
17478515
Brief Bioinform. 2006 Mar;7(1):37-47
16761363
PLoS Biol. 2008 Apr 15;6(4):e83
18416601
Nature. 2008 Mar 27;452(7186):429-35
18344982
Nucleic Acids Res. 2007 Jan;35(Database issue):D760-5
17099226
Hum Mol Genet. 2006 Apr 15;15 Spec No 1:R9-16
16651375