Incorporating prior knowledge to facilitate discoveries in a genome-wide association study on age-related macular degeneration
© Lin and Lee; licensee BioMed Central Ltd. 2010
Received: 15 December 2009
Accepted: 28 January 2010
Published: 28 January 2010
Substantial genotyping data produced by current high-throughput technologies have brought opportunities and difficulties. With the number of single-nucleotide polymorphisms (SNPs) going into millions comes the harsh challenge of multiple-testing adjustment. However, even with the false discovery rate (FDR) control approach, a genome-wide association study (GWAS) may still fall short of discovering any true positive gene, particularly when it has a relatively small sample size.
To counteract such a harsh multiple-testing penalty, in this report, we incorporate findings from previous linkage and association studies to re-analyze a GWAS on age-related macular degeneration. While previous Bonferroni correction and the traditional FDR approach detected only one significant SNP (rs380390), here we have been able to detect seven significant SNPs with an easy-to-implement prioritized subset analysis (PSA) with the overall FDR controlled at 0.05. These include SNPs within three genes: CFH, CFHR4, and SGCD.
Based on the success of this example, we advocate using the simple method of PSA to facilitate discoveries in future GWASs.
Substantial genotyping data produced by current high-throughput technologies have brought opportunities and difficulties. High-density genotyping platforms have been developed in a hope that underlying disease-associated genes can be identified through denser and denser collections of single-nucleotide polymorphism (SNP) data. However with the number of SNPs going into millions comes the harsh challenge of multiple-testing adjustment. To counteract multiple-testing penalty incurred by testing such a large number of SNPs, some genome-wide association studies (GWASs) responded by taking a large sample size--with the number of study subjects soaring into thousands, tens of thousands, or even more .
There are two approaches for multiple-testing adjustments. One is controlling the family-wise error rate (FWER), the other is controlling the false discovery rate (FDR) [2, 3]. The FWER is defined as the probability of at least one type I error. Among methods for controlling FWER, the Bonferroni correction is the best known approach, although it is very conservative. Holm's step-down procedure  is less conservative than the classical Bonferroni correction. The FWER can also be controlled by the resampling-based P-value adjustment procedure. Compared with controlling the FWER, controlling the FDR is usually a more powerful approach. However, even with the FDR approach, a GWAS may still fall short of discovering any true positive gene, particularly when it has a relatively small sample size. When testing simultaneously for a huge number of SNPs, even true positive SNPs would have difficulty in standing out among all the noise, based on a straight (and brutal) comparison of their p values. GWAS on age-related macular degeneration (AMD) is a good example, and we will show this in this paper.
The above simple FDR approach has been further extended to dependent tests and to tests with prior information . The false discovery control with P-value weighting [5, 6] can improve power when the assignment of weights (based on previous linkage evidence ) is adequate, but there is some power loss when the weights are poorly assigned. Sun et al.'s  stratified false discovery control is another approach. They partitioned all SNPs into two subsets based on minor-allele frequencies (MAFs), and then the FDR control is applied to the two subsets respectively. However, as pointed out by Li et al. , MAFs have little relevance with biological information and thus partitioning SNPs based on MAFs probably may not improve much power. To address this issue, Li et al.  proposed a 'prioritized subset analysis' (PSA). The PSA makes clever use of available prior knowledge, either of the linkage information, the biological information or both. We will show that the PSA can greatly facilitate discoveries in GWASs, with a demonstration on an AMD data.
Materials: a GWAS on Age-related Macular Degeneration (AMD)
AMD is a genetically complex disorder. The heritability was estimated to range from 46% to 71%. Klein et al.  reported an AMD data set containing 96 AMD cases and 50 controls. Of all the 116,204 genotyped SNPs, 99,317 SNPs were informative (MAF ≥ 1%) and conformed to Hardy-Weinberg equilibrium (with Hardy-Weinberg exact p value ≥ 0.05 in the 50 controls). Following Klein et al. , we test for allelic association with disease status on each SNP.
Prioritized Subset Analysis
To facilitate discoveries in GWASs, we turned to a new method of 'prioritized subset analysis' (PSA) . To perform a PSA, a researcher based on his/her prior biological knowledge first picks from among all SNPs under study, a certain number of SNPs likely to be the true positives. He/she then places those selected SNPs in a 'prioritized subset' and those remaining in a 'non-prioritized' subset. The FDR control is then applied to these two subsets separately, and the significant results are harvested from both the two subsets.
We took findings from previous genome-wide linkage and association studies on AMD as our prior knowledge to prioritize SNPs. Our prioritization process is detailed below.
We first incorporated evidence of linkage (with LOD score >2.0) based on previous linkage studies [10–16]. We obtained the physical position of each D-number marker (listed in Table 1) from the Gene Location website http://genecards.weizmann.ac.il/geneloc/index.shtml. Then SNPs within 500 kb from each D-number marker were prioritized.
Genes or markers to be prioritized, in the prioritized subset analysis
No. of SNPs in the prioritized subset
D2S1356, D2S1394, D2S1384
D3S1768, D3S1304, D3S3045
D5S820, GATA12A08, D5S1506
HLA, C2-CFB, VEGF, ELOVL4, SOD2
In the end, a total of 639 SNPs were prioritized, and the remaining 98,678 SNPs, non-prioritized. We then applied the PSA with the FDR being controlled at 0.05, for both the prioritized subset and the non-prioritized subset. We used Storey and Tibshirani's  smoothing spline approach provided by the package 'fdrtool'  to estimate the proportions of true negative SNPs.
Bonferroni Correction and Traditional FDR Approach
Results of the AMD data set
P value *
5.12 × 10-4
3.01 × 10-4
5.40 × 10-8
1.59 × 10-4
7.20 × 10-5
3.69 × 10-4
Prioritized Subset Analysis
The PSA identified a total of seven significant SNPs (all from the prioritized subset) (Table 2). These include SNPs within three genes: CFH, CFHR4, and SGCD. By using the PSA method, we have been able to detect six additional significant SNPs (in two additional genes), compared to the Bonferroni approach (the method used by Klein et al. ) or the traditional FDR approach. Two of the three significant genes found in this study, CFH and CFHR4, are located in a chromosomal region (1q31-1q32) having been most replicated in previous AMD studies. The remaining one significant SGCD gene had not been previously reported to be AMD-related, though. However, we notice that previous animal studies showed the SGCD gene is related to vascular abnormalities in mice . This might suggest a link of SGCD to neovascular AMD in humans.
All the seven significant SNPs are from the prioritized subset. To evaluate how well the FDR is controlled in our prioritized subset, we further estimated the permutation-based FDR  in this subset. We randomly permuted the data and calculated the null P values - for the i th SNP in the b th permutation (i = 1,...,639). Through B permutations, the number of false positives (FP) is estimated as , where d = 5.12 × 10-4 is the largest P value of the seven significant SNPs (see Table 2). We took B = 100,000 and obtained = 0.225. The permutation-based FDR in the prioritized subset is thus estimated as 0.225/7 = 0.032, which is still less than our FDR control level of 0.05, suggesting a satisfactory FDR control in this subset.
Prior information can come from a researcher's biological knowledge, or findings of data other than that provided in the current study. But one should not 'snoop' his/her data at hand for the prior knowledge. If one naively prioritizes those SNPs with the smallest p values in the study data, the actual overall FDR would no longer be properly controlled. To avoid such bias, we searched findings of other data to build our 'prior knowledge', before seeing the analysis results of individual SNPs in the current AMD data set. At that time, we did know that rs380390 is a significant SNP in the AMD data set which can withstand a FWER control of 0.05 . But the chromosomal region around rs380390 had already been replicated by many previous linkage studies [10–15] (all published before Klein et al. ). And so, prioritizing chromosomal region around rs380390 won't constitute an act of data snooping.
Around a particular gene, how large a chromosomal region should be prioritized is also an issue. Because of the consistent findings in the CFH gene, both from genome-wide linkage analyses [10–15] and case-control studies [19–21], we prioritized SNPs within 1 Mb from the CFH gene. Other evidence of linkage and associations are relatively unconfirmed by prior studies, so we prioritized SNPs within 500 kb and 50 kb, respectively. Because linkage is a coarse mapping whereas association is a fine mapping, in general a wider region of SNPs should be prioritized for a linkage peak. Admittedly, there is no absolute criterion for choosing the sizes of prioritized regions. No matter how large a chromosomal region is prioritized, the FDR within subsets should be controlled at the desired level, and this can be verified by estimating the permutation-based FDR .
In recent GWASs, a commonly used approach to incorporate prior knowledge is to calculate the Bayes factors [1, 30]. However, to estimate the Bayes factors, the prior distributions and the effect sizes should be carefully specified . This may limit its applicability. By contrast, the PSA method used in this paper can feed on prior knowledge that is only rudimentary (we need only to decide beforehand whether a particular SNP is more likely a true positive or a true negative, but don't need to know exactly how likely). And there is almost no penalty for poor guessing . In this paper, we demonstrated that such a simple dichotomization followed by a simple PSA can greatly facilitate discoveries in a GWAS on AMD.
Note that we did not recruit any more subjects or type any more SNPs beyond what Klein et al.  had done. The only thing we did is to incorporate prior knowledge about AMD into the analysis. And we see this input of knowledge is rather powerful (six/two additional significant SNPs/genes were identified in the same AMD case-control data). One may question that our input of knowledge and the subsequent partition of SNPs into two subsets to be tested separately and harvested combinedly are making easier (and perhaps too easier) for the SNPs to come out. But we should emphasize that we did not loosen our FDR control in any way. The total seven significant SNPs found in this re-analysis have an overall 0.05 FDR attached to them, much the same way with the one SNP rs380390 originally found in Klein et al.  had a 0.05 FDR attached to it. And we believe that researchers will find no difficulties to choose seven SNPs or just one--that is, under the same FDR criteria.
The PSA approach is rather powerful and is easy to implement. Based on the success of our re-analysis of Klein et al's GWAS on AMD, we advocate using PSA to facilitate discoveries in future GWASs.
We would like to thank the anonymous reviewers for their constructive comments. We also thank Dr. Josephine Hoh for kindly providing the AMD data set. This study was supported by National Science Councils, Taiwan.
- WTCCC: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678. 10.1038/nature05911.View ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.Google Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-9445. 10.1073/pnas.1530509100.PubMed CentralPubMedView ArticleGoogle Scholar
- Holm S: A simple sequentially rejective multiple test procedure. Scand J Statist. 1979, 6: 65-70.Google Scholar
- Genovese C, Roeder K, Wasserman L: False discovery control with P-value weighting. Biometrika. 2006, 93: 509-524. 10.1093/biomet/93.3.509.View ArticleGoogle Scholar
- Roeder K, Bacanu SA, Wasserman L, Devlin B: Using linkage genome scans to improve power of association in genome scans. Am J Hum Genet. 2006, 78 (2): 243-252. 10.1086/500026.PubMed CentralPubMedView ArticleGoogle Scholar
- Sun L, Craiu RV, Paterson AD, Bull SB: Stratified false discovery control for large-scale hypothesis testing with application to genome-wide association studies. Genet Epidemiol. 2006, 30 (6): 519-530. 10.1002/gepi.20164.PubMedView ArticleGoogle Scholar
- Li C, Li M, Lange EM, Watanabe RM: Prioritized subset analysis: improving power in genome-wide association studies. Hum Hered. 2008, 65 (3): 129-141. 10.1159/000109730.PubMedView ArticleGoogle Scholar
- Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, Haynes C, Henning AK, SanGiovanni JP, Mane SM, Mayne ST: Complement factor H polymorphism in age-related macular degeneration. Science. 2005, 308 (5720): 385-389. 10.1126/science.1109557.PubMed CentralPubMedView ArticleGoogle Scholar
- Abecasis GR, Yashar BM, Zhao Y, Ghiasvand NM, Zareparsi S, Branham KE, Reddick AC, Trager EH, Yoshida S, Bahling J: Age-related macular degeneration: a high-resolution genome scan for susceptibility loci in a population enriched for late-stage disease. Am J Hum Genet. 2004, 74 (3): 482-494. 10.1086/382786.PubMed CentralPubMedView ArticleGoogle Scholar
- Iyengar SK, Song D, Klein BE, Klein R, Schick JH, Humphrey J, Millard C, Liptak R, Russo K, Jun G: Dissection of genomewide-scan data in extended families reveals a major locus and oligogenic susceptibility for age-related macular degeneration. Am J Hum Genet. 2004, 74 (1): 20-39. 10.1086/380912.PubMed CentralPubMedView ArticleGoogle Scholar
- Klein ML, Schultz DW, Edwards A, Matise TC, Rust K, Berselli CB, Trzupek K, Weleber RG, Ott J, Wirtz MK: Age-related macular degeneration. Clinical features in a large family and linkage to chromosome 1q. Arch Ophthalmol. 1998, 116 (8): 1082-1088.PubMedView ArticleGoogle Scholar
- Majewski J, Schultz DW, Weleber RG, Schain MB, Edwards AO, Matise TC, Acott TS, Ott J, Klein ML: Age-related macular degeneration--a genome scan in extended families. Am J Hum Genet. 2003, 73 (3): 540-550. 10.1086/377701.PubMed CentralPubMedView ArticleGoogle Scholar
- Seddon JM, Santangelo SL, Book K, Chong S, Cote J: A genomewide scan for age-related macular degeneration provides evidence for linkage to several chromosomal regions. Am J Hum Genet. 2003, 73 (4): 780-790. 10.1086/378505.PubMed CentralPubMedView ArticleGoogle Scholar
- Weeks DE, Conley YP, Tsai HJ, Mah TS, Schmidt S, Postel EA, Agarwal A, Haines JL, Pericak-Vance MA, Rosenfeld PJ: Age-related maculopathy: a genomewide scan with continued evidence of susceptibility loci within the 1q31, 10q26, and 17q25 regions. Am J Hum Genet. 2004, 75 (2): 174-189. 10.1086/422476.PubMed CentralPubMedView ArticleGoogle Scholar
- Jun G, Klein BE, Klein R, Fox K, Millard C, Capriotti J, Russo K, Lee KE, Elston RC, Iyengar SK: Genome-wide analyses demonstrate novel loci that predispose to drusen formation. Invest Ophthalmol Vis Sci. 2005, 46 (9): 3081-3088. 10.1167/iovs.04-1360.PubMedView ArticleGoogle Scholar
- Scholl HP, Fleckenstein M, Charbel Issa P, Keilhauer C, Holz FG, Weber BH: An update on the genetics of age-related macular degeneration. Mol Vis. 2007, 13: 196-205.PubMed CentralPubMedGoogle Scholar
- Haddad S, Chen CA, Santangelo SL, Seddon JM: The genetics of age-related macular degeneration: a review of progress to date. Surv Ophthalmol. 2006, 51 (4): 316-363. 10.1016/j.survophthal.2006.05.001.PubMedView ArticleGoogle Scholar
- Edwards AO, Ritter R, Abel KJ, Manning A, Panhuysen C, Farrer LA: Complement factor H polymorphism and age-related macular degeneration. Science. 2005, 308 (5720): 421-424. 10.1126/science.1110189.PubMedView ArticleGoogle Scholar
- Haines JL, Hauser MA, Schmidt S, Scott WK, Olson LM, Gallins P, Spencer KL, Kwan SY, Noureddine M, Gilbert JR: Complement factor H variant increases the risk of age-related macular degeneration. Science. 2005, 308 (5720): 419-421. 10.1126/science.1110359.PubMedView ArticleGoogle Scholar
- Hageman GS, Anderson DH, Johnson LV, Hancox LS, Taiber AJ, Hardisty LI, Hageman JL, Stockman HA, Borchardt JD, Gehrs KM: A common haplotype in the complement regulatory gene factor H (HF1/CFH) predisposes individuals to age-related macular degeneration. Proc Natl Acad Sci USA. 2005, 102 (20): 7227-7232. 10.1073/pnas.0501536102.PubMed CentralPubMedView ArticleGoogle Scholar
- Rivera A, Fisher SA, Fritsche LG, Keilhauer CN, Lichtner P, Meitinger T, Weber BH: Hypothetical LOC387715 is a second major susceptibility gene for age-related macular degeneration, contributing independently of complement factor H to disease risk. Hum Mol Genet. 2005, 14 (21): 3227-3236. 10.1093/hmg/ddi353.PubMedView ArticleGoogle Scholar
- Jakobsdottir J, Conley YP, Weeks DE, Mah TS, Ferrell RE, Gorin MB: Susceptibility genes for age-related maculopathy on chromosome 10q26. Am J Hum Genet. 2005, 77 (3): 389-407. 10.1086/444437.PubMed CentralPubMedView ArticleGoogle Scholar
- Gold B, Merriam JE, Zernant J, Hancox LS, Taiber AJ, Gehrs K, Cramer K, Neel J, Bergeron J, Barile GR: Variation in factor B (BF) and complement component 2 (C2) genes is associated with age-related macular degeneration. Nat Genet. 2006, 38 (4): 458-462. 10.1038/ng1750.PubMed CentralPubMedView ArticleGoogle Scholar
- Jakobsdottir J, Conley YP, Weeks DE, Ferrell RE, Gorin MB: C2 and CFB genes in age-related maculopathy and joint action with CFH and LOC387715 genes. PLoS ONE. 2008, 3 (5): e2199-10.1371/journal.pone.0002199.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen YH, Liu CK, Chang SC, Lin YJ, Tsai MF, Chen YT, Yao A: GenoWatch: a disease gene mining browser for association study. Nucleic Acids Res. 2008, W336-340. 10.1093/nar/gkn214. 36 Web Server
- Strimmer K: A unified approach to false discovery rate estimation. BMC Bioinformatics. 2008, 9: 303-10.1186/1471-2105-9-303.PubMed CentralPubMedView ArticleGoogle Scholar
- Dye WW, Gleason RL, Wilson E, Humphrey JD: Altered biomechanical properties of carotid arteries in two mouse models of muscular dystrophy. J Appl Physiol. 2007, 103 (2): 664-672. 10.1152/japplphysiol.00118.2007.PubMedView ArticleGoogle Scholar
- Xie Y, Pan W, Khodursky AB: A note on using permutation-based false discovery rate estimates to compare different analysis methods for microarray data. Bioinformatics. 2005, 21 (23): 4280-4288. 10.1093/bioinformatics/bti685.PubMedView ArticleGoogle Scholar
- Wakefield J: Reporting and interpretation in genome-wide association studies. Int J Epidemiol. 2008, 37 (3): 641-653. 10.1093/ije/dym257.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.