Mining gold dust under the genome wide significance level: a two-stage approach to analysis of GWAS.
We propose a two-stage approach to analyze genome-wide association data in order to identify a set of promising single-nucleotide polymorphisms (SNPs). In stage one, we select a list of top signals from single SNP analyses by controlling false discovery rate. In stage two, we use the least absolute shrinkage and selection operator (LASSO) regression to reduce false positives. The proposed approach was evaluated using simulated quantitative traits based on genome-wide SNP data on 8,861 Caucasian individuals from the Atherosclerosis Risk in Communities (ARIC) Study. Our first stage, targeted at controlling false negatives, yields better power than using Bonferroni-corrected significance level. The LASSO regression reduces the number of significant SNPs in stage two: it reduces false-positive SNPs and it reduces true-positive SNPs also at simulated causal loci due to linkage disequilibrium. Interestingly, the LASSO regression preserves the power from stage one, i.e., the number of causal loci detected from the LASSO regression in stage two is almost the same as in stage one, while reducing false positives further. Real data on systolic blood pressure in the ARIC study was analyzed using our two-stage approach which identified two significant SNPs, one of which was reported to be genome-significant in a meta-analysis containing a much larger sample size. On the other hand, a single SNP association scan did not yield any significant results.
Division of Biostatistics, Washington University School of Medicine, Saint Louis, Missouri 63110-1093, USA. gang@wubios.wustl.edu
Gang
G
Shi
Eric
E
Boerwinkle
Alanna C
AC
Morrison
C Charles
CC
Gu
Aravinda
A
Chakravarti
D C
DC
Rao
D003198:Computer Simulation; D005188:False Negative Reactions; D005189:False Positive Reactions; D055106:Genome-Wide Association Study; D006801:Humans; D015810:Linkage Disequilibrium; D015233:Models, Statistical; D017720:Molecular Epidemiology / Q000379:methods; D020641:Polymorphism, Single Nucleotide; D012044:Regression Analysis; D015203:Reproducibility of Results
Genetics. 2005 Oct;171(2):783-90
15956674
Genet Epidemiol. 2008 May;32(4):381-5
18348202
Nat Genet. 2009 Jun;41(6):677-87
19430479
Nat Genet. 2009 Dec;41(12):1325-9
19915573
BMC Proc. 2009 Dec 15;3 Suppl 7:S25
20018015
PLoS Genet. 2009 Jul;5(7):e1000564
19609347
Nat Genet. 2009 Apr;41(4):407-14
19305409
Hypertension. 2003 Feb;41(2):207-10
12574083
BMC Proc. 2009 Dec 15;3 Suppl 7:S2
20018009
Am J Epidemiol. 1989 Apr;129(4):687-702
2646917
BMC Proc. 2009 Dec 15;3 Suppl 7:S67
20018061
J Comput Chem. 2002 Jan 15;23(1):77-83
11913391
Nat Genet. 2009 Jan;41(1):56-65
19060906
Proc Natl Acad Sci U S A. 2009 Jun 9;106(23):9362-7
19474294
Bioinformatics. 2003 Feb 12;19(3):368-75
12584122
Hum Mol Genet. 2010 Feb 1;19(3):553-62
19897590
BMC Proc. 2009 Dec 15;3 Suppl 7:S4
20018031
Blood. 2010 Jan 7;115(1):94-6
19880490
BMC Proc. 2009 Dec 15;3 Suppl 7:S89
20018085
BMC Proc. 2007;1 Suppl 1:S60
18466561
Genet Epidemiol. 2009 Jan;33(1):45-53
18618761
BMC Proc. 2009 Dec 15;3 Suppl 7:S27
20018017
Pac Symp Biocomput. 2008;:190-200
18229686
BMC Proc. 2009 Dec 15;3 Suppl 7:S62
20018056
Am J Hum Genet. 2007 Sep;81(3):559-75
17701901
Proc Natl Acad Sci U S A. 2001 Apr 24;98(9):5116-21
11309499
Proc Natl Acad Sci U S A. 2003 Aug 5;100(16):9440-5
12883005
Nat Genet. 2008 Feb;40(2):189-97
18193044
Bioinformatics. 2009 Mar 15;25(6):714-21
19176549
Am J Hum Genet. 2007 Sep;81(3):540-51
17701899
Nat Genet. 2009 Jun;41(6):666-76
19430483
BMC Bioinformatics. 2007 Jun 29;8:230
17603887
Am J Hum Genet. 2008 Feb;82(2):375-85
18252218
Nature. 2001 May 10;411(6834):199-204
11346797