Skip to main content
  • Research article
  • Open access
  • Published:

Disease-associated alleles in genome-wide association studies are enriched for derived low frequency alleles relative to HapMap and neutral expectations

Abstract

Background

Genome-wide association studies give insight into the genetic basis of common diseases. An open question is whether the allele frequency distributions and ancestral vs. derived states of disease-associated alleles differ from the rest of the genome. Characteristics of disease-associated alleles can be used to increase the yield of future studies.

Methods

The set of all common disease-associated alleles found in genome-wide association studies prior to January 2010 was analyzed and compared with HapMap and theoretical null expectations. In addition, allele frequency distributions of different disease classes were assessed. Ages of HapMap and disease-associated alleles were also estimated.

Results

The allele frequency distribution of HapMap alleles was qualitatively similar to neutral expectations. However, disease-associated alleles were more likely to be low frequency derived alleles relative to null expectations. 43.7% of disease-associated alleles were ancestral alleles. The mean frequency of disease-associated alleles was less than randomly chosen CEU HapMap alleles (0.394 vs. 0.610, after accounting for probability of detection). Similar patterns were observed for the subset of disease-associated alleles that have been verified in multiple studies. SNPs implicated in genome-wide association studies were enriched for young SNPs compared to randomly selected HapMap loci. Odds ratios of disease-associated alleles tended to be less than 1.5 and varied by frequency, confirming previous studies.

Conclusions

Alleles associated with genetic disease differ from randomly selected HapMap alleles and neutral expectations. The evolutionary history of alleles (frequency and ancestral vs. derived state) influences whether they are implicated in genome-wide assocation studies.

Peer Review reports

Background

The onset of affordable high-throughput genotyping technology has enabled association studies to be conducted on a genome-wide scale, and multiple successes have occurred using this approach [1, 2]. Notable examples include genes associated with LDL cholesterol levels [3], colorectal cancer [4], and type 1 diabetes [5]. Despite these successes, genome-wide association studies (GWAS) have been unable to account for the majority of heritable variation for most diseases [6, 7]. One reason for the mixed success of GWAS is that the efficacy of such studies depends upon the underlying genetic architecture of traits [8–10]. Also relevant is the accuracy of the common disease-common variant hypothesis. Under this formulation, complex diseases are caused by high frequency alleles. By contrast, the genetic heterogeneity (rare allele-major effect) hypothesis proposes that distinct low-frequency alleles are responsible for the same trait in different individuals. Regardless of the relative validity of each of these hypotheses [11, 12], there are now enough genome-wide association studies to be able to say something about the frequencies and ancestral or derived state of disease-associated alleles relative to HapMap and neutral expectations. A previous study gave an initial glimpse into the nature of disease-associated alleles, finding a median allele frequency of 0.40 [13]. However, less than 50 replicated SNPs were analyzed in the review article by Iles. Multiple studies have found that the allele frequencies of disease-associated SNPs do not significantly differ across populations compared to random SNPs [14, 15]. Additionally, alleles associated with genetic disease are underrepresented in intergenic regions and overrepresented in nonsynonymous sites and 5 kb promoter regions [16]. An open question is whether the allele frequency distribution and ancestral vs. derived states of disease-associated alleles differ from the rest of the genome. What types of disease-associated alleles have been found in GWAS to date?

Alleles can be classified by frequency, ancestral vs. derived state, and relative disease risk. Disease-associated marker alleles need not be causal; they may be linked to alleles that are actually responsible for increased risk. Importantly, the alleles characterized in this paper are marker alleles. Because statistical power is maximized at intermediate allele frequencies [17], alleles detected in a GWAS are unlikely to be rare. Alleles can also be classified as ancestral or derived, where ancestral alleles are shared with closely related species and derived alleles are not [18]. There is currently a lack of published studies indicating whether disease-associated alleles are enriched for ancestral or derived states. Finally, alleles can be classified by their relative risk (measured as an odds ratio). In many GWAS, disease-associated alleles increase risk by only a modest amount, i.e. odds ratios are less than 1.5 [1]. However, it is unknown whether these odds ratios vary by ancestral vs. derived state. In addition, the first study to detect a particular association can overestimate the odds ratio, a phenomenon dubbed the "winner's curse" [19].

The aim of this study was to determine the characteristics of disease-associated alleles and compare this data to null expectations from HapMap data and the neutral theory. Allele frequency distributions and ancestral vs. derived states were obtained for every disease-associated allele found in genome-wide association studies prior to January 2010. The hypothesis that disease-associated alleles do not differ from the rest of the genome [20] was also tested.

Methods

Null expectations: HapMap data

Alleles from the HapMap dataset were obtained to serve as a baseline of genomic diversity. These alleles were subsequently weighted by the probability of detection in a GWAS. The set of all HapMap SNPs from data release #27 were downloaded via the HapMart tool [21]. This build included merged data from phases II and III of the International HapMap Project [22]. A Perl script was then used to randomly select 1000 unique SNPs from this file. The positions of HapMap SNPs were tested to ensure that linkage disequilibrium was minimal (distances between SNPs were at least 200 kb). Because allele frequencies can vary from population to population, and the majority of GWAS use European and European-American cases and controls, CEU allele frequencies were used in this study. This allowed HapMap dataset to act as a control for demographic processes. While disease prevalence and allele frequencies vary among populations [23], population-level differences in allele frequencies are similar for disease-associated SNPs and random genomic SNPs [14]. An additional consideration is that SNP discovery protocols may bias the HapMap dataset towards high frequency alleles [24, 25]. For each SNP, alleles were chosen at random after weighting by allele frequency. Thus, a SNP with allele frequencies of 0.80 and 0.20 would have an 80% chance of yielding the major allele. Because high frequency alleles are more likely to be ancestral [26], randomly chosen alleles from the HapMap dataset are more likely to be ancestral than derived.

Outgroups (such as chimpanzees) enable SNPs to be polarized and ancestral states to be inferred. Ancestral alleles were determined via BLAST searches of disease-associated SNP regions against the chimpanzee genome. In addition, the single nucleotide polymorphism database, dbSNP, contains information on the ancestral state of SNPs at http://www.ncbi.nlm.nih.gov/projects/SNP[27]. Ancestral alleles in dbSNP were inferred via parsimony [28]. SNPs were only selected if ancestral allele states could be inferred. HapMap alleles were then binned by allele frequency and ancestral vs. derived state. Unlike disease-associated alleles, randomly selected HapMap alleles are not weighted by their probability of detection. Because of this, comparisons between the allele frequency distribution of disease-associated alleles and null expectations incorporated the probability that a particular allele is detectable in a GWAS. Phase three of the HapMap project used the Illumina Human1 M and Affymetrix SNP 6.0 platforms (the same platforms used in many GWAS). This indicates that the HapMap dataset served as an appropriate control for GWAS data.

Null expectations: neutral theory

Population genetic theory was used to test whether disease-associated SNPs differ from neutral expectations. These expectations were subsequently weighted by the probability of detection in a GWAS. An infinite sites model was used to obtain the theoretical allele frequency distribution of neutral loci. Under this model, novel mutations occur at different nucleotide sites [29]. Marker loci were assumed to be diallelic and lack recurring mutations. In constant sized populations the probability of observing a derived neutral allele at a particular frequency is inversely proportional to allele frequency [30]. The parameter C in the equations below is a normalizing constant, and "unweighted" refers to the fact that these probability density functions do not incorporate the probability of detection in a GWAS. x is the frequency of the disease-associated allele at a marker locus.

P ( x = X , unweighted | derived ) = C x
(1)
P ( x = X , unweighted | ancestral ) = C 1 − x
(2)

The probability density of Equation 1 goes to infinity as allele frequencies go to zero. Because of this, polymorphisms were only considered if the minor allele frequency was above some arbitrary threshhold frequency, d. Allele frequencies were allowed to range from d to 1-d, with the parameter d arbitrarily set equal to 0.025. MATLAB [31] simulations verified the accuracy of Equation 1 (see Figure 1A). In these simulations, alleles were binomially sampled each generation and frequencies of derived alleles were recorded. Upon fixation or loss, a single derived allele was allowed to enter the population. Simulations were run for 107 time steps with a population size of 104 individuals.

Figure 1
figure 1

Neutral expectations, statistical power, and probability of detection. A) Unweighted probability density of derived allele frequencies under the neutral theory. Calculations use Equation 1 and are represented by a black line. MATLAB simulation data is represented by grey circles. B) Statistical power of GWAS plotted as a function of odds ratio and allele frequency. Darker shading indicates greater statistical power. Parameter values: 2500 cases and controls, a p-value of 5 × 10-8, complete linkage, multiplicative dominance. C) Probability of detecting a disease-associated allele. Three distributions of odds ratios are considered: uniform (evenly distributed odds ratios between 1.0 and 1.5, dashed grey line), Dirac delta (every allele has an odds ratio of 1.25, solid grey line), and normal (mean 1.25, standard deviation = 0.1, dashed black line). Equation 6 is plotted as a solid black line.

The probability that an allele is ancestral is equal to its probability of fixation [32, 33]. For neutral loci, the probability that a randomly chosen allele is ancestral is simply its allele frequency [26].

P ( ancestral | x = X ) = x
(3)

The right-hand sides of Equations 2 and 3 were multiplied and integrated from d to 1-d. This expression was then normalized by dividing by the total probability density, giving the overall probability that a neutral allele is ancestral (unweighted by the chance of detection in a GWAS).

P (ancestral, unweighted) = ∫ d 1 − d x 1-x d x ∫ d 1 − d 1 1-x d x
(4)

After integration and extensive algebra:

P (ancestral, unweighted) = 1 + 2 d − 1 ln ( 1 − d ) − ln ( d )
(5)

Probability of detection (statistical power calculations)

Because statistical power varies by allele frequency [17], computer simulations were used to calculate the probability of detecting an association between disease and a marker allele at a particular frequency. Statistical power calculations were obtained via the Windows program QUANTO [34, 35]. This program numerically calculates power for a variety of experimental designs, but it does not explicitly take into account ages of SNPs and ancestral vs. derived states. Statistical power is a function of linkage disequilibrium between causal and marker alleles (see Appendix). The following parameter values were used: 2500 cases and controls, a p-value of 5 × 10-8, complete linkage, multiplicative dominance. Statistical power was calculated for odds ratios ranging from 1 to 2 (at increments of 0.01) and allele frequencies ranging from 0 to 1 (at increments of 0.01). As indicated by Figure 1B, alleles with low odds ratios can only be detected at intermediate allele frequencies.

Although the true distribution of odds ratios for the set of all disease-associated alleles is unknown, most disease-associated alleles increase relative risk by small amount (i.e. odds ratios are less than 1.5) [1]. Three different distributions of odds ratios were considered: uniform (evenly distributed odds ratios between 1.0 and 1.5), Dirac delta (every allele has an odds ratio of 1.25), and normal (mean 1.25, standard deviation = 0.1). Plots of statistical power vs. allele frequency were similar for each of these distributions (Figure 1C). A simple expression yields a reasonable estimate of statistical power given the parameter values listed above. For mathematical simplicity, subsequent calculations assume that the probability of detecting a disease-associated allele is:

P ( detection| x = X , parameter values listed above ) ≈ 2 x ( 1 − x )
(6)

It is coincidental that the expression in Equation 6 also gives the expected heterozygosity. Different sample sizes and/or odds ratios would yield different expressions for the probability of detection.

Null expectations: weighted frequency distributions

HapMap and neutral expectations were weighted by the probability of detection to enable fair comparisons with disease-associated alleles. Probability densities of each allele frequency bin were multiplied by the expression in Equation 6 and normalized. This resulted in elevated probabilities of observing intermediate frequency alleles. For neutral expectations.

P ( x = X , weighted | derived ) = 2 ( 1 − x )
(7)
P ( x = X , weighted | ancestral ) = 2 x
(8)

The right-hand sides of Equations 3 and 8 were multiplied and integrated from d to 1-d. This expression was then normalized by dividing by the total probability density, giving the overall probability that a neutral allele is ancestral (weighted by the chance of detection in a GWAS).

P (ancestral, weighted) = ∫ d 1 − d 2x 2 d x ∫ d 1 − d 2x  d x
(9)

After integration and extensive algebra:

P (ancestral, weighted) = 2 3 ( d 2 − d + 1 )
(10)

Empirical data: disease-associated alleles

The set of all disease-associated alleles found prior to January 1, 2010 was obtained to investigate whether these alleles differed from the rest of the genome. An excellent database of GWAS and disease-associated SNPs exists online at http://www.genome.gov/gwastudies and it was used in this study [16, 36]. This Catalog of Published Genome-Wide Association Studies includes every disease-associated SNP to date. Criteria for inclusion in this database included p-values < 10-5 and a minimum of 100,000 SNPs tested in the initial stage of a study [16]. The set of all genome-wide association studies prior to January 1, 2010 spans 486 papers and a total of 2186 disease-associated SNPs. Some of these SNPs were present in multiple studies, and in many cases the disease-associated allele was not listed in the database. Allele frequencies in control populations were obtained from NHGRI's Catalog of Published Genome-Wide Association Studies. When allele frequency data were absent from this database, CEU HapMap frequencies were used. If a particular SNP was associated with diseases in multiple studies, mean allele frequencies were calculated. The ancestral vs. derived state of each allele was determined via dbSNP. When ancestral allele state could not be inferred, SNPs were omitted from the dataset. After taking into account ancestral vs. derived states of alleles, a total of 1143 disease-associated SNPs remained. Of these SNPs, 530 had odds ratio data.

For comparisons with null expectations, disease-associated SNPs were binned into 10% allele frequency intervals (see Figure 3). Differences between the allele-frequency distributions of disease-associated alleles and null expectations were assessed via χ2 goodness-of-fit tests. Relative magnitudes of disease-associated and control allele frequencies were compared via Mann-Whitney U tests. Proportions of ancestral alleles were compared via binomial tests, with a null hypothesis that equal proportions of disease-associated and control alleles were ancestral. In addition, 95% confidence intervals of mean allele frequency were found by sampling with replacement. This bootstrap analysis was performed in MATLAB (100000 replicates) [31]. It is possible that alleles found in multiple studies have different characteristics than alleles found in a single study. Because of this, the mean frequency and evolutionary history of replicated alleles (alleles implicated in multiple studies) were compared with the overall patterns of disease-associated alleles. A total of 142 replicated alleles had frequency and evolutionary history data.

Alleles implicated in GWAS were also sorted into seven different phenotypic classes (cancer, cardiovascular, metabolism, miscellaneous disease, morphological, and neurological). Some alleles were associated with multiple phenotypic classes. The morphological class included alleles that were not technically associated with any disease. Instead, they were associated with traits such as height and hair color. 92 of 1143 GWAS alleles were implicated in studies of non-European populations, and the remaining 1051 alleles were re-analyzed to determine whether this had any effect. Because of the small number of alleles implicated in studies of non-European populations, additional analysis was not conducted on these 92 GWAS alleles. See Additional file 1 for a list of SNPs, allele frequencies, ancestral vs. derived states, odds ratios, and phenotypic class.

To test whether the properties of disease-associated alleles changed over time, disease-associated were binned into six-month intervals by date of first publication. Mean frequency, probability ancestral, and odds ratio data were calculated for each time interval. To determine whether genotyping platforms biased the properties of alleles, data from Affymetrix, Illumina, and Perlegen arrays were compared. Many studies used multiple genotyping platforms, and characteristics of disease-associated alleles in these studies were also analyzed. The number of genotyped SNPs passing quality control varied in each study. To assess whether this had any effect, disease-associated alleles were binned into three intervals corresponding to the number of genotyped SNPs ( < 500,000, between 500,000 and 1,000,000, and >1,000,000). These numbers include imputed SNPs.

Estimated ages of SNPs

Ages of SNPs were calculated to determine if detectable associations occur more often in young or old SNPs. Estimates of SNP ages were calculated from the allele frequency distributions of HapMap and GWAS alleles. Under the neutral theory [37], the expected age of a SNP (Ï„) is

E ( τ ) = − 2 x 1 − x ln ( x )
(11)

x in Equation 11 refers to the frequency of the derived allele and time is measured in units of 2Ne generations (where Ne is the effective population size). SNPs with low frequency derived alleles tend to be younger SNPs. However, variance in Ï„ tends to be quite large and allele frequencies only give a rough estimate of the SNP age. Because of this, the cumulative probability density function [38, 39] was used.

P ( τ ≤ t) ≅ ( 1 − x ) − 1 + n / ( 1 + n t / 2 )
(12)

Sample size (n) in Equation 12 was arbitrarily set equal to 2500. The derivative of Equation 12 with respect to t was taken for a range of allele frequencies (0.05 to 0.95 at intervals of 0.10) and SNP ages (0 to 8 Ne generations at intervals of 0.04 Ne generations). Allele frequency distributions in Figures 2B and 3A were then used to generate the expected distributions of SNP ages for weighted HapMap loci and GWAS loci.

Figure 2
figure 2

Allele frequency distributions for null expectations. Allele frequencies are binned into 0.10 intervals. Derived probabilities are labelled in black and ancestral probabilities are labelled in grey. Neutral expectations use a polymorphism threshold (d) of 0.025. A) Allele frequency distributions for HapMap alleles prior to detection (n = 1000). B) Allele frequency distributions for HapMap alleles after weighting by the probability of detection (n = 1000). C) Theoretical allele frequency distributions for neutral alleles prior to detection. D) Theoretical allele frequency distributions for neutral alleles after weighting by the probability of detection.

Results

Null expectations

HapMap alleles were qualitatively similar to neutral alleles. In both cases, randomly selected alleles were likely to be ancestral and high frequency (see Table 1). High frequency (0.80-1.00) ancestral alleles and low frequency (0.00-0.10) derived alleles were the most common types of alleles for neutral and HapMap datasets (Figure 2). Unweighted proportions of ancestral alleles were comparable to estimates from whole-genome data (0.707 for a French individual) [40]. Weighting by the probability of detection in a GWAS increased the probability of observing intermediate frequency alleles and decreased the probability that an allele was ancestral.

Table 1 Disease-associated alleles vs. null expectations

Despite these similarities, there were important quantitative differences between HapMap alleles and neutral expectations. HapMap alleles were less likely to be ancestral than alleles under the neutral theory (0.623 vs. 0.741 for the unweighted scenario, and 0.568 vs. 0.650 for the weighted scenario). HapMap alleles were also more likely to be found at intermediate frequencies, and goodness-of-fit tests indicate that differences existed between the allele frequency distributions of HapMap and neutral alleles (p-value < 10-10 for both the unweighted and weighted scenarios, χ2 test with 19 d.f.). The relative lack of HapMap loci with a minor allele frequency < 0.1 may be due to ascertainment bias, as SNPs identified from a small panel of individuals are known to have an excess of intermediate frequency alleles [25].

Sensitivity analysis of the HapMap dataset suggests the absence of selection bias in the 1000 randomly selected SNPs. The mean frequency of minor alleles was 0.219 for the HapMap dataset used in this paper, compared to 0.214 and 0.199 for additional sets of 1000 SNPs from the HapMap (CEU) and Perlegen (EUR) databases, respectively.

Disease-associated alleles

Empirical data from genome-wide association studies indicated that a majority of disease-associated alleles were derived. Out of 1143 unique SNPs, disease-associated alleles were ancestral in 499 cases and derived in 644 cases. As shown in Table 1 the proportion of ancestral alleles was less than HapMap and neutral theory expectations (p-value < 0.0001, binomial test for each comparison). Alleles shared with chimpanzees were less likely to be associated with genetic disease than alleles that are not shared with chimpanzees.

The majority of disease-associated alleles had frequencies below 0.50 (Figure 3). Goodness of fit tests indicated that the empirical allele frequency distribution differed from both the HapMap and neutral expectations (p-value < 0.0001 in both cases, χ2 test with 19 d.f.). Allele frequency bins containing the highest proportion of disease-associated alleles were derived alleles with frequencies between 0.00 and 0.50 and ancestral alleles with frequencies between 0.30 and 0.50. The mean allele frequency of disease-associated alleles was 0.394, and comparisons disease-associated alleles had lower frequencies than HapMap and neutral expectations (p-value < 0.0001, Mann-Whitney U test). This is consistent with population genetics theory that predicts neutral variants linked to deleterious alleles should be found at lower frequencies [41]. 95% bootstrap confidence intervals of mean allele frequency ranged from 0.3802 to 0.4076. When data from studies of non-European populations were excluded, patterns were largely unchanged (mean frequency was 0.397 and the proportion of ancestral alleles was 0.428). Observed patterns were largely insensitive to p-value thresholds of GWAS: the correlation between allele frequency and negative log p-value was -0.0733, and disease-associated alleles found at p-values < 10-8 had a mean frequency of 0.404. Mean allele frequencies for the complete set of GWAS alleles were 0.463 for ancestral and 0.340 for derived alleles, indicating that disease-associated alleles tended to be minor alleles. This is in contrast to null expectations, where ancestral alleles tended to be major alleles. The shape of the disease-associated allele frequency distribution did not resemble either null expectation, suggesting that additional factors were involved. Misidentification of ancestral states can result in an excess of high frequency alleles [42]. Because the dataset of disease-associated alleles had a deficiency of high-frequency alleles, this suggests that ancestral states were correctly inferred.

Figure 3
figure 3

Allele frequency distributions for GWAS data. Allele frequencies are binned into 0.10 intervals. Derived probability densities are labelled in black and ancestral probability densities are labelled in grey. A) All disease-associated alleles (n = 1143). B) Disease-associated alleles that have been implicated in multiple studies (n = 142).

Alleles implicated in multiple studies showed similar patterns to the overall set of disease-associated alleles. Replicated alleles had a mean frequency of 0.396 and a 0.437 probability of being ancestral. The allele frequency distribution of replicated alleles also exhibited an excess of rare alleles relative to null expectations (Figure 3B).

Differences between disease-associated alleles and the rest of the genome can be due either to properties of loci or properties of alleles. Characteristics of loci were revealed in derived frequency distributions (Figure 4A). By contrast, the frequency distributions in (Figures 2 and 3) revealed characteristics of both alleles and loci. Goodness of fit tests indicated that derived frequency distributions differed for neutral expectations, HapMap SNPs, and disease-associated SNPs (p-value < 0.0001 for each pairwise comparison, χ2 tests with 9 d.f.). However, derived frequency distributions were more similar than allele frequency distributions (compare Figures 2, 3, and 4). In addition, disease-associated SNPs and weighted HapMap SNPs had similar mean derived frequencies (0.426 vs. 0.432). This suggests that much of the difference between disease-associated alleles and the rest of the genome was due to properties of alleles rather than loci.

Figure 4
figure 4

Characteristics of different types of loci. A) Derived frequency distributions. Derived allele frequencies are binned into 0.10 intervals and probability densities for different types of loci are indicated by shading (neutral expectations in light grey, HapMap SNPs in black, and disease-associated GWAS SNPs in dark grey). Similar numbers of loci were analyzed for each data type (1000 for neutral expectations and HapMap SNPs, and 1143 for GWAS SNPs). B) Estimated ages of SNPs. The probability density for HapMap SNPs (weighted by probability of detection) is labelled black and the probability density for disease-associated GWAS SNPs is labelled dark grey. Probability densities were obtained via Equation 12, and calculated at intervals of 0.04 Ne generations.

Different phenotypic classes

Similar patterns were observed for each of the phenotypic classes (Table 1). In each case, disease-associated alleles had lower frequencies than HapMap and neutral expectations (p-value < 0.0001 for each comparison, Mann-Whitney U test). Regardless of phenotypic class, disease-associated alleles were more likely to be derived alleles than randomly selected HapMap alleles and neutral expectations (p-value < 0.05 for each comparison, binomial test). Alleles associated with cardiovascular disease were most likely to be low frequency derived alleles. Alleles associated with neurological disease had the highest mean allele frequency, and alleles associated with morphological traits were more likely to be ancestral. However, differences among phenotypic classes were smaller than the differences between each phenotypic class and null expectations (p-value < 0.05 for allele frequency and ancestral vs. derived data, One-way ANOVA).

Additional controls

Characteristics of disease-associated alleles were independent of publication date and genotyping platform. Mean allele frequencies for each six-month interval were between 0.363 and 0.421. Similarly, the proportion of ancestral alleles ranged from 0.414 to 0.481. Median odds ratios for each six-month interval were between 1.24 to 1.365. Temporal trends were not observed for any of these characteristics. Allele frequencies of disease-associated alleles and proportion of ancestral alleles were similar for different genotyping platforms (Table 2). Median odds ratios were also similar for each manufacturer (1.28 for Affymetrix, 1.26 for Illumina, and 1.425 for Perlegen). Although the number of genotyped SNPs did not appear to affect mean allele frequency, the proportion ancestral alleles was slightly less for studies with >1,000,000 genotyped SNPs (Table 2). Overall, differences between platforms were smaller than differences between disease-associated alleles and null expectations.

Table 2 Disease-associated alleles from different genotyping platforms

Estimated ages of SNPs

Genome-wide association studies were enriched for young SNPs compared to randomly selected HapMap loci (Figure 4B). Mean ages of SNPs were estimated to be 2.78 Ne generations for HapMap loci and 2.74 Ne generations for GWAS loci. However, the spread around the mean was quite large for each locus type. The probability densities in Figure 4B reveal that SNPs arising in the last 1 Ne generations were over-represented in the GWAS dataset. This occurred because disease-associated alleles had an excess of low frequency derived alleles (the same sorts of alleles that tend to occur in young SNPs). Recall that these calculations assumed that SNPs were neutral. Directional selection would reduce the expected ages of SNPs for both types of loci [39]. Similarly, ascertainment bias due to the small sample size of the SNP discovery panel [43] can affect both types of SNPs.

Odds ratios

The findings of previous studies [1, 13, 16] were confirmed: most disease-associated alleles only increase disease risk by only a moderate amount (Figure 5). This is consistent with expectations from population genetics theory as alleles with high odds ratios are expected to have a higher fitness burden [44]. In addition, odds ratios of disease-associated alleles varied by frequency. The median odds ratio of ancestral alleles was 1.28 and the median odds ratio of derived alleles was 1.32. This indicates that the average effect of derived alleles was slightly larger than ancestral alleles (p-value < 0.05, Mann-Whitney U test). Low frequency alleles with high odds ratios tended to be derived alleles, but this was simply a byproduct of GWAS being enriched for derived alleles. Few intermediate frequency disease-associated alleles had high odds ratios. 26.8% of disease-associated SNPs with a minor allele frequency ≤ 0.2 had an odds ratio > 2, while only 8.0% of SNPs with a minor allele frequency > 0.2 had an odds ratio > 2 (p-value < 0.0001, Fisher's exact test). It is difficult to detect statistical associations between diseases and low frequency marker alleles, suggesting that the high odds ratios observed for SNPs with a minor allele frequency ≤ 0.2 were indicative of a "winner's curse." As only those alleles that are statistically significant were reported, published odds ratios may be overestimated [45].

Figure 5
figure 5

Odds ratios for ancestral and derived alleles as a function of allele frequency. Derived alleles are represented by black circles, and ancestral alleles are represented by grey circles. A total of 530 alleles have odds ratio data.

Discussion

Disease-associated alleles were more likely to be low frequency derived alleles than neutral and CEU HapMap expectations. Patterns were similar for alleles associated with different phenotypic classes. These findings were independent of publication date and genotyping platform. SNPs implicated in genome-wide association studies were enriched for young SNPs compared to randomly selected HapMap loci. One caveat is that that the majority of published studies to date have used european populations, and it is unclear whether these patterns will apply to other populations.

Statistical power, sample sizes, and allele frequencies

Because statistical power is minimal at extreme allele frequencies, it is not surprising that most disease-associated alleles have minor allele frequencies greater than 0.1 (Figure 3). The relative inability of GWAS to explain the high heritabilities of many diseases [7] suggests that many genes responsible for common diseases might actually be at undetectably low frequencies. Simulations reveal that much of the fitness variance associated with genetic diseases can be due to very low frequency, large-effect alleles [46]. Theoretical work also indicates that rare causal alleles can create associations that are credited to common marker alleles (a phenomenon that has been called "synthetic association") [47]. Alternatively, common genetic diseases may be due to multiple allele of small effect, synergistic epistasis, and/or genotype-by-environment interactions [6]. To avoid only detecting associations with intermediate frequency alleles, larger sample sizes are needed. Increasing the number of genotyped SNPs also results in a greater likelihood of detecting an association [48], but all SNPs are not equally informative. In addition, the results of this study suggest that future GWAS may benefit from the inclusion of many young SNPs with low frequency derived alleles.

Genetic background and linkage phase of causal and marker alleles

Genetic background and linkage phase may explain why disease-associated alleles were enriched for derived alleles. Consider the following thought experiment: Two alleles already segregate at a marker locus when a causal mutation occurs at a nearby locus. If the causal mutation occurs in an ancestral genetic background, only a small proportion of disease-associated marker alleles will be in phase with the causal mutation. This is because ancestral alleles tend to have higher frequencies than derived alleles [26]. As a result, the (P(B|A) - P(B|a))2 term in Equation A1 (see Appendix) tends to be smaller when causal alleles are in phase with an ancestral allele at the marker locus. All other things being equal, linkage disequilibrium (r 2) and statistical power are greater when causal mutations occur in a derived genetic background.

Recombination breaks down linkage disequilibrium between causal and marker alleles over time, reducing the likelihood of statistical associations. This is consistent with the finding that SNPs showing detectable associations with genetic disease are younger than randomly selected HapMap SNPs. Genome-wide association studies are also less likely to be successful if mutations occur multiple times at a causal locus. This is because the causal alleles can be found in multiple genetic backgrounds, reducing statistical associations between causal alleles and marker alleles. Population heterogeneity can also be an issue, as causal mutations can occur in different genetic backgrounds in different populations [49, 50].

Natural selection against disease alleles

Natural selection against deleterious alleles may also cause disease-associated alleles to differ from the rest of the genome. Population genetic theory indicates that marker alleles linked to low fitness causal alleles are expected to be uncommon [30, 51]. This is in agreement with the finding that alleles associated with genetic disease tend to be found at lower frequencies than randomly selected HapMap alleles. Natural selection may also be able to explain differences in SNP age between disease-associated alleles and HapMap alleles.

The efficacy of selection varies for different genetic diseases. Because fitness refers to the expected contribution to the next generation's gene pool, diseases with late onset are likely to be found at higher frequencies. In addition, allele frequency distributions are shaped by the evolutionary history of a disease. Selection pressures that change over time can allow disease alleles to segregate at intermediate frequencies [52]. The genetic architecture of a disease also affects the strength of selection: Tajima's D and Dn/Ds ratios reveal that the signature of selection is stronger for Mendelian diseases than complex genetic diseases [53]. This may explain why different phenotypic classes have slightly different profiles (Table 1). However, alleles associated with morphological traits (as opposed to disease) differ from null expectations. Additional factors than natural selection may be required to explain why GWAS alleles differ from the rest of the genome.

Can population size changes explain the observed patterns?

Disease-associated alleles are enriched for derived low frequency alleles, a pattern that can occur when populations increase in size [54, 55]. However, population expansions affect the frequency distributions of all alleles, not just disease-associated alleles. Because HapMap data do not show an excess of derived low frequency alleles relative to neutral expectations (Figure 2), this indicates that population size changes alone cannot fully explain the characteristic patterns of disease-associated alleles.

Conclusion

At these initial stages, alleles found via genome-wide association studies tend to be low hanging fruit. However, there is strong evidence that disease-associated alleles differ from the rest of the genome. Costs of microarray-based genotyping platforms are decreasing, and as the number of SNPs analyzed per individual increases, so too does the chance of detecting direct associations between disease and causal SNPs (rather than merely linked marker SNPs). In addition, direct sequencing is becoming increasingly affordable and it allows previously unknown SNPs to be identified. Because of this, direct sequencing of large genomic regions adjacent to disease-associated marker alleles is advisable. Direct sequencing of GWAS-informed regions can also be combined with familial inheritance patterns to improve genetic linkage analyses. This brings up an intriguing question: are causal alleles for a particular trait more likely to be ancestral or derived? Also, how can GWAS be planned to maximize the likelihood of detecting candidate genes associated with a particular disease? Combining the perspectives of genetic epidemiology and evolutionary genetics allows these questions to be answered.

Appendix

Statistical power in GWAS is a function of odds ratios and the amount of linkage disequilibrium between causal alleles and marker alleles. The relevant measure of linkage disequilibrium in this case is r 2 [56]. Consider a causal locus with two segregating alleles (A and a), and a linked marker locus with two segregating alleles (B and b). r 2 is defined as:

r 2 = ( P ( B | A ) − P ( B | a ) ) 2 y ( 1 − y ) x ( 1 − x ) ,
(A1)

where P(B|A) and P(B|a) are the probabilities that a haplotype has the marker allele B given a linked causal allele of A or a, y is the frequency of the disease allele at the causal locus, and x is the frequency of the disease-associated allele at the marker locus.

References

  1. Altshuler D, Daly MJ, Lander ES: Genetic mapping in human disease. Science. 2008, 322 (5903): 881-888. 10.1126/science.1156409.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Wellcome Trust Case Control Consortium: 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.

    Article  Google Scholar 

  3. Kotowski IK, Pertsemlidis A, Luke A, Cooper RS, Vega GL, Cohen JC, Hobbs HH: A spectrum of PCSK9 Alleles contributes to plasma levels of low-density lipoprotein cholesterol. American Journal of Human Genetics. 2006, 78 (3): 410-422. 10.1086/500615.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Tomlinson I, Webb E, Carvajal-Carmona L, Broderick P, Kemp Z, Spain S, Penegar S, Chandler I, Gorman M, Wood W, et al: A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21. Nature Genetics. 2007, 39 (8): 984-988. 10.1038/ng2085.

    Article  CAS  PubMed  Google Scholar 

  5. Todd JA, Walker NM, Cooper JD, Smyth DJ, Downes K, Plagnol V, Bailey R, Nejentsev S, Field SF, Payne F, et al: Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nature Genetics. 2007, 39 (7): 857-864. 10.1038/ng2068.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 2009, 461 (7265): 747-753. 10.1038/nature08494.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Maher B: Personal genomes: The case of the missing heritability. Nature. 2008, 456 (7218): 18-21. 10.1038/456018a.

    Article  CAS  PubMed  Google Scholar 

  8. Clark AG, Boerwinkle E, Hixson J, Sing CF: Determinants of the success of whole-genome association testing. Genome Res. 2005, 15 (11): 1463-1467. 10.1101/gr.4244005.

    Article  CAS  PubMed  Google Scholar 

  9. Clarke AJ, Cooper DN: GWAS: heritability missing in action?. European Journal of Human Genetics. 2010, 18: 859-861. 10.1038/ejhg.2010.35.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Moore JH, Williams SM: Epistasis and its implications for personal genetics. Am J Hum Genet. 2009, 85 (3): 309-320. 10.1016/j.ajhg.2009.08.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Goldstein DB: Common genetic variation and human traits. N Engl J Med. 2009, 360 (17): 1696-1698. 10.1056/NEJMp0806284.

    Article  CAS  PubMed  Google Scholar 

  12. Hirschhorn JN: Genomewide association studies--illuminating biologic pathways. N Engl J Med. 2009, 360 (17): 1699-1701. 10.1056/NEJMp0808934.

    Article  CAS  PubMed  Google Scholar 

  13. Iles MM: What can genome-wide association studies tell us about the genetics of common disease?. PLoS Genet. 2008, 4 (2): e33-10.1371/journal.pgen.0040033.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Myles S, Davison D, Barrett J, Stoneking M, Timpson N: Worldwide population differentiation at disease-associated SNPs. BMC Med Genomics. 2008, 1: 22-10.1186/1755-8794-1-22.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Lohmueller KE, Mauney MM, Reich D, Braverman JM: Variants associated with common disease are not unusually differentiated in frequency across populations. Am J Hum Genet. 2006, 78 (1): 130-136. 10.1086/499287.

    Article  CAS  PubMed  Google Scholar 

  16. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA. 2009, 106 (23): 9362-9367. 10.1073/pnas.0903103106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang WYS, Barratt BJ, Clayton DG, Todd JA: Genome-wide association studies: Theoretical and practical concerns. Nature Reviews Genetics. 2005, 6 (2): 109-118. 10.1038/nrg1522.

    Article  CAS  PubMed  Google Scholar 

  18. Hacia JG, Fan JB, Ryder O, Jin L, Edgemon K, Ghandour G, Mayer RA, Sun B, Hsie L, Robbins C, et al: Determination of ancestral alleles for human single-nucleotide polymorphisms using high-density oligonucleotide arrays. Nat Genet. 1999, 22 (2): 164-167. 10.1038/9674.

    Article  CAS  PubMed  Google Scholar 

  19. Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN: Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003, 33 (2): 177-182. 10.1038/ng1071.

    Article  CAS  PubMed  Google Scholar 

  20. Wang WY, Pike N: The allelic spectra of common diseases may resemble the allelic spectrum of the full genome. Med Hypotheses. 2004, 63 (4): 748-751. 10.1016/j.mehy.2003.12.057.

    Article  CAS  PubMed  Google Scholar 

  21. HapMart. [http://hapmart.hapmap.org/BioMart/martview/]

  22. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal S, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449 (7164): 851-861. 10.1038/nature06258.

    Article  CAS  PubMed  Google Scholar 

  23. Rotimi CN, Jorde LB: Ancestry and disease in the age of genomic medicine. N Engl J Med. 2010, 363 (16): 1551-1558. 10.1056/NEJMra0911564.

    Article  PubMed  Google Scholar 

  24. Ganapathy G, Uyenoyama MK: Site frequency spectra from genomic SNP surveys. Theor Popul Biol. 2009, 75 (4): 346-354. 10.1016/j.tpb.2009.04.003.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Nielsen R, Hubisz MJ, Clark AG: Reconstituting the frequency spectrum of ascertained single-nucleotide polymorphism data. Genetics. 2004, 168 (4): 2373-2382. 10.1534/genetics.104.031039.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Watterson GA, Guess HA: Is the most frequent allele the oldest?. Theor Popul Biol. 1977, 11 (2): 141-160. 10.1016/0040-5809(77)90023-5.

    Article  CAS  PubMed  Google Scholar 

  27. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29 (1): 308-311. 10.1093/nar/29.1.308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Spencer CC, Deloukas P, Hunt S, Mullikin J, Myers S, Silverman B, Donnelly P, Bentley D, McVean G: The influence of recombination on human genetic diversity. PLoS Genet. 2006, 2 (9): e148-10.1371/journal.pgen.0020148.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Kimura M: The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations. Genetics. 1969, 61 (4): 893-903.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Sethupathy P, Hannenhalli S: A Tutorial of the Poisson Random Field Model in Population Genetics. Advances in Bioinformatics. 2008, 2008: Article ID 257864-10.1155/2008/257864.

    Article  Google Scholar 

  31. Mathworks: MATLAB 7. 2005, Natick, MA: The Mathworks

    Google Scholar 

  32. Choi SC, Stone EA, Kishino H, Thorne JL: Estimates of natural selection due to protein tertiary structure inform the ancestry of biallelic loci. Gene. 2008

    Google Scholar 

  33. Taylor JE: The common ancestor process for a Wright-Fisher diffusion. Electron J Probab. 2007, 12: 808-847.

    Article  Google Scholar 

  34. Gauderman WJ: Sample size requirements for matched case-control studies of gene-environment interaction. Stat Med. 2002, 21 (1): 35-50. 10.1002/sim.973.

    Article  PubMed  Google Scholar 

  35. Gauderman WJ, Morrison JM: QUANTO 1.2: A computer program for power and sample size calculations for genetic epidemiology studies. 2007

    Google Scholar 

  36. Johnson AD, O'Donnell CJ: An open access database of genome-wide association results. BMC Med Genet. 2009, 10: 6-10.1186/1471-2350-10-6.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Kimura M, Ohta T: The age of a neutral mutant persisting in a finite population. Genetics. 1973, 75 (1): 199-212.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Griffiths RC, Tavaré S: The age of a mutation in a general coalescent tree. Stoch Models. 1998, 14: 273-295. 10.1080/15326349808807471.

    Article  Google Scholar 

  39. Slatkin M, Rannala B: Estimating the age of alleles by use of intraallelic variability. Am J Hum Genet. 1997, 60 (2): 447-458.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Green R, Krause J, Briggs A, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz M, et al: A Draft Sequence of the Neanderthal Genome. Science. 2010, 328: 710-722. 10.1126/science.1188021.

    Article  CAS  PubMed  Google Scholar 

  41. Bamshad M, Wooding SP: Signatures of natural selection in the human genome. Nat Rev Genet. 2003, 4 (2): 99-111. 10.1038/nrg999.

    Article  CAS  PubMed  Google Scholar 

  42. Hernandez RD, Williamson SH, Bustamante CD: Context dependence, ancestral misidentification, and spurious signatures of natural selection. Mol Biol Evol. 2007, 24 (8): 1792-1800. 10.1093/molbev/msm108.

    Article  CAS  PubMed  Google Scholar 

  43. Clark AG, Hubisz MJ, Bustamante CD, Williamson SH, Nielsen R: Ascertainment bias in studies of human genome-wide polymorphism. Genome Res. 2005, 15 (11): 1496-1502. 10.1101/gr.4107905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Brookfield JF: Q&A: promise and pitfalls of genome-wide association studies. BMC Biol. 2010, 8: 41.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Zhong H, Prentice RL: Correcting "winner's curse" in odds ratios from genomewide association findings for major complex human diseases. Genet Epidemiol. 2009

    Google Scholar 

  46. Eyre-Walker A: Genetic architecture of a complex trait and its implications for fitness and genome-wide association studies. Proc Natl Acad Sci USA. 2010, 107 (Suppl 1): 1752-1756. 10.1073/pnas.0906182107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Dickson SP, Wang K, Krantz I, Hakonarson H, Goldstein DB: Rare variants create synthetic genome-wide associations. PLoS Biol. 2010, 8 (1): e1000294-10.1371/journal.pbio.1000294.

    Article  PubMed  PubMed Central  Google Scholar 

  48. McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008, 9 (5): 356-369. 10.1038/nrg2344.

    Article  CAS  PubMed  Google Scholar 

  49. McClellan J, King MC: Genetic heterogeneity in human disease. Cell. 2010, 141 (2): 210-217. 10.1016/j.cell.2010.03.032.

    Article  CAS  PubMed  Google Scholar 

  50. Campbell MC, Tishkoff SA: African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping. Annu Rev Genomics Hum Genet. 2008, 9: 403-433. 10.1146/annurev.genom.9.081307.164258.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Pritchard JK, Cox NJ: The allelic architecture of human disease genes: common disease-common variant...or not?. Hum Mol Genet. 2002, 11 (20): 2417-2423. 10.1093/hmg/11.20.2417.

    Article  CAS  PubMed  Google Scholar 

  52. Neel JV: Diabetes mellitus: a 'thrifty' genotype rendered detrimental by 'progress'?. Am J Hum Genet. 1962, 14: 353-362.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Blekhman R, Man O, Herrmann L, Boyko AR, Indap A, Kosiol C, Bustamante CD, Teshima KM, Przeworski M: Natural selection on genes that underlie human disease susceptibility. Curr Biol. 2008, 18 (12): 883-889. 10.1016/j.cub.2008.04.074.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129 (2): 555-562.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Griffiths RC, Tavare S: Sampling theory for neutral alleles in a varying environment. Philos Trans R Soc Lond B Biol Sci. 1994, 344 (1310): 403-410. 10.1098/rstb.1994.0079.

    Article  CAS  PubMed  Google Scholar 

  56. Pritchard JK, Przeworski M: Linkage disequilibrium in humans: models and data. Am J Hum Genet. 2001, 69 (1): 1-14. 10.1086/321275.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references

Acknowledgements

I thank D. Dykhuizen, W. Eanes, E. Hatchwell, L. Jung, E. Sezgin, M. Talbert, J. True, and S. Yeh for helpful comments and suggestions. The constructive comments of two reviewers (J. Zhao and N. Timpson) greatly improved the quality of this paper. This work was supported by a National Institutes of Health Predoctoral Training Grant (5 T32 GM007964-24).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Joseph Lachance.

Additional information

Competing interests

The author declares that they have no competing interests.

Authors' contributions

JL conceived the study, participated in the analysis of the data, and wrote the manuscript.

Electronic supplementary material

12920_2010_192_MOESM1_ESM.XLS

Additional file 1: GWAS data. This file is a Microsoft Excel spreadsheet that contains allele frequencies, ancestral vs. derived state, and phenotypic class for each disease-associated allele analyzed in this study. (XLS 150 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lachance, J. Disease-associated alleles in genome-wide association studies are enriched for derived low frequency alleles relative to HapMap and neutral expectations. BMC Med Genomics 3, 57 (2010). https://doi.org/10.1186/1755-8794-3-57

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1755-8794-3-57

Keywords