Optimal two-stage genome-wide association designs based on false discovery rate
Genome-wide association studies are likely to be conducted in large scale in the near future. In such studies, searching over hundreds of thousands of markers for the few ones that are associated with disease brings out the multiple-hypothesis testing problem in its severe form. We explore, in a two-stage design, how the use of false discovery rate (FDR) can alleviate the burden of a prohibitively strict significance level for single marker tests and still control the number of false positive findings, when there is more than one causal variant. FDR is the expected proportion of false positives among all significant findings. It can be approximated by (1-p0 p 0 )? ? /[(1-p0 p 0 )?+p0 ? + p 0 (1-? ? )], where p0 p 0 is the proportion of true causal markers, ? ? is the type I error rate and 1-? ? the power of a two-stage study. When 500,000 SNPs are genotyped in the first stage with fixed SNP array and the most significant SNPs are genotyped in the second stage with standard but 20 times more expensive high-throughput techniques, up to 20% savings in the minimum genotyping cost is achieved for p0 p 0 in the range of 10-5 10 - 5 to 510-4 5  10 - 4 and FDR in the range of 0.05 to 0.7, compared to when Bonferroni-corrected significance level is used. In terms of sample size, the saving is up to 60%. However, these savings come at a cost of more false positive findings.
multiple comparison false positive report probability whole-genome scan linkage disequilibrium (ld) cost-efficient