Search tips
Search criteria 


Logo of plosbiolPLoS BiologySubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
PLoS Biol. 2013 September; 11(9): e1001661.
Published online 2013 September 17. doi:  10.1371/journal.pbio.1001661
PMCID: PMC3775722

Generalization and Dilution of Association Results from European GWAS in Populations of Non-European Ancestry: The PAGE Study

Greg Gibson, Academic Editor


The vast majority of genome-wide association study (GWAS) findings reported to date are from populations with European Ancestry (EA), and it is not yet clear how broadly the genetic associations described will generalize to populations of diverse ancestry. The Population Architecture Using Genomics and Epidemiology (PAGE) study is a consortium of multi-ancestry, population-based studies formed with the objective of refining our understanding of the genetic architecture of common traits emerging from GWAS. In the present analysis of five common diseases and traits, including body mass index, type 2 diabetes, and lipid levels, we compare direction and magnitude of effects for GWAS-identified variants in multiple non-EA populations against EA findings. We demonstrate that, in all populations analyzed, a significant majority of GWAS-identified variants have allelic associations in the same direction as in EA, with none showing a statistically significant effect in the opposite direction, after adjustment for multiple testing. However, 25% of tagSNPs identified in EA GWAS have significantly different effect sizes in at least one non-EA population, and these differential effects were most frequent in African Americans where all differential effects were diluted toward the null. We demonstrate that differential LD between tagSNPs and functional variants within populations contributes significantly to dilute effect sizes in this population. Although most variants identified from GWAS in EA populations generalize to all non-EA populations assessed, genetic models derived from GWAS findings in EA may generate spurious results in non-EA populations due to differential effect sizes. Regardless of the origin of the differential effects, caution should be exercised in applying any genetic risk prediction model based on tagSNPs outside of the ancestry group in which it was derived. Models based directly on functional variation may generalize more robustly, but the identification of functional variants remains challenging.

Author Summary

The number of known associations between human diseases and common genetic variants has grown dramatically in the past decade, most being identified in large-scale genetic studies of people of Western European origin. But because the frequencies of genetic variants can differ substantially between continental populations, it's important to assess how well these associations can be extended to populations with different continental ancestry. Are the correlations between genetic variants, disease endpoints, and risk factors consistent enough for genetic risk models to be reliably applied across different ancestries? Here we describe a systematic analysis of disease outcome and risk-factor–associated variants (tagSNPs) identified in European populations, in which we test whether the effect size of a tagSNP is consistent across six populations with significant non-European ancestry. We demonstrate that although nearly all such tagSNPs have effects in the same direction across all ancestries (i.e., variants associated with higher risk in Europeans will also be associated with higher risk in other populations), roughly a quarter of the variants tested have significantly different magnitude of effect (usually lower) in at least one non-European population. We therefore advise caution in the use of tagSNP-based genetic disease risk models in populations that have a different genetic ancestry from the population in which original associations were first made. We then show that this differential strength of association can be attributed to population-dependent variations in the correlation between tagSNPs and the variant that actually determines risk—the so-called functional variant. Risk models based on functional variants are therefore likely to be more robust than tagSNP-based models.


In the past six years, genome-wide association studies (GWAS) have revealed thousands of common polymorphisms (tagSNPs) associated with a wide variety of traits and diseases, particularly as study sample sizes have increased from thousands to hundreds of thousands of subjects. Typically GWAS analyses stratify on genetic ancestry, because many polymorphism allele frequencies differ by ancestral group, easily producing false positive associations for traits that also correlate with genetic ancestry. The large majority of GWAS results reported to date derive from analyses in populations of European ancestry (EA) [1],[2]. Although GWAS in Asian populations in particular are becoming more common [3][6], it remains important to understand the degree to which the magnitude and direction of allelic effects generalize across diverse populations [7][10]. The multi-ethnic PAGE consortium [11] provides a unique opportunity to assess GWAS generalization across multiple non-EA populations and multiple traits.

Results and Discussion

Subject and genotyping panel selection for the PAGE consortium have been described elsewhere [11],[12]. In brief, a panel of 68 common polymorphisms previously reported to associate with body mass index (BMI) [13], type 2 diabetes (T2D) [14], or lipid levels [15] was genotyped in up to 14,492 self-reported African Americans (AA), 8,202 Hispanic Americans (HA), 5,425 Asian Americans (AS), 6,186 Native Americans (NA), 1,801 Pacific Islanders (PI), and 37,061 EA (for details, see Materials and Methods, Table S1 and Table S2). We also analyzed a subset of 5863 AA from PAGE who were genotyped on the Illumina Metabochip, which contains approximately 200,000 SNPs densely focused on 257 regions with reported GWAS associations to traits that include lipids, BMI, and T2D [16].

For a replication analysis it would be overly conservative to use the Bonferroni correction, so the Benjamini-Hochberg method [17] was applied to assess replication of previous EA reports in the PAGE EA population. Reported effects in EA were replicated for 51 out of the 68 index SNPs at a 5% FDR. Power to replicate at most of these 68 SNPs far exceeded 80%; 16 of the 17 SNPs that did not replicate exceeded 80% power to replicate the reported effect size, and the 17th exceeded 70% power, as described previously [13][15]. The originally reported effect sizes tend to be less extreme for these seventeen index SNPs, but in 63 out of 79 comparisons between non-EA and EA populations involving these 17 SNPs, the direction of effect was the same in EA and non-EA groups (p<10−5 for the null hypothesis of random effects in either direction, data in Table 1 column “Index SNPs Not Replicated in EA”). Only 79 of the 85 possible pairwise comparisons against EA were assessed, because some of the 17 SNPs were not genotyped in all five non-EA populations. Thus, it appears likely that most of the 17 failures to replicate represent weak effects that were underpowered in PAGE EA, rather than false-positive primary reports. Therefore, all 68 index SNPs were carried forward in the generalization analysis.

Table 1
Summary of direction and strength of β relative to EA.

In all non-EA groups, we observe significantly more effects in the same direction as in EA than expected under the null hypothesis, ranging from 68% in Asians to 88% in Hispanics (p<0.001 in all non-EA groups, Table 1 and Figure 1). Even in the relatively small Pacific Islander population (N = 1801), where only four index SNPs were significantly associated with reported traits, 48 out of 62 effects were in the same direction as EA (p<0.001), so in larger samples from this population we would expect additional loci to generalize. Although a higher proportion of effects in the opposite direction of EA was observed in Asians and Pacific Islanders, the opposite effects were neither significantly different from no effect, nor significantly different from the observed effect in the EA population. This suggests that the greater number of effects in the opposite direction observed in these smallest groups simply reflects greater uncertainty in estimating effect sizes for these populations, rather than any true trend toward opposite effects. The proportion of effects in the same direction as EA was similar across all non-EA populations, suggesting that for at least 70% of index SNPs, a significant effect in a consistent direction will ultimately be observed in non-EA populations of adequate size.

Figure 1
Generalization analysis in the PAGE populations.

Whereas the direction of effect was consistent between EA and non-EA populations, the magnitude of effect varied considerably, consistent with prior meta-analyses of generalization [18]. Because effect sizes were correlated among non-EA populations, we applied the Benjamini-Hochberg method within each population to identify index SNPs with significantly inconsistent effects between EA and non-EA populations. Inconsistent effects (βpopβEA at 5% FDR) were observed for 17 of 68 index SNPs in at least one non-EA population (Table 2 and Table S2, see Box 1 for definitions). Inconsistent effects were most frequent in the AA population (12 out of 68 loci), but examples were also observed in Pacific Islanders and Native Americans. Although most effects were consistent between EA and non-EA populations, the relatively high frequency with which differential effects were observed in non-EA populations suggests that genetic risk models derived from GWAS in EA will predict risk less reliably in non-EA populations, particularly AA. Consequently, caution should be exercised in applying risk models based upon risk variants genotyped outside of the ethnic background in which they were derived [19], regardless of the factors causing the observed variation between populations,.

Table 2
Summary of generalization results.

Box 1. Definitions

ßpop: The effect size of a given SNP in linear or logistic regression models for a specific PAGE population. Where available and when allowed by the informed consent protocols, effect sizes were estimated in models that included estimated genetic ancestry, as previously reported (see Text S1).

ßEA: The effect size of a given SNP in the PAGE EA population. We use the PAGE EA effect size for comparisons to PAGE non-EA populations rather than the original report for two reasons: to minimize the impact of winner's curse on these comparisons, and because several of the SNPs genotyped in PAGE were proxies strongly correlated with the original tagSNP, and might not match the reported effect size.

We define replicated SNPs as SNPs with direction of effect consistent with the original report in EA, and significant ßEA in PAGE (using α = 0.05 as the threshold for hypothesis rejection, unadjusted for multiple testing, as these are considered specific prior hypothesis being validated).

When comparing two populations, the direction of effect can be either the samepop and ßEA are either both positive, or both negative) or opposite (either ßpop or ßEA is positive, and the other is negative). Magnitude of effect was evaluated only for SNPs that replicated in EA and can be either stronger (|ßEA|<|ßpop|), the same (|ßEA| = pop|), or weaker (|ßEA|>|ßpop|).

In order to describe the generalization of EA findings to non-EA populations, SNPs are categorized in terms of (a) significance in the non-EA population and (b) consistency between non-EA and EA populations. Here we use the Benjamini-Hochberg procedure to adjust for testing up to 68 SNPs in each non-EA population.

Significant SNPs reject the null hypothesis of no effect in the non-EA population (ßpop≠0) at q = 0.05.

Inconsistent SNPs reject the hypothesis of equal effect size in EA and non-EA populations (ßpop≠ßEA) at q = 0.05.

Combining these parameters yields four categories of generalization:

  • Ambiguous SNPs are neither significant in non-EA, nor inconsistent between non-EA and EA.
  • Differential SNPs are not significant in non-EA, but are inconsistent between non-EA and EA.
  • Differentially Generalized SNPs are significant in non-EA, and inconsistent between non-EA and EA.
  • Strictly Generalized SNPs are significant in non-EA, and consistent between non-EA and EA.

Four index SNPs showed differentially generalized effects (ßpop≠ßEA and ßpop≠0). Two of these did not replicate in EA (rs7578597 and rs7961581 for T2D in NA) so consistency of direction cannot accurately be inferred. Direction of effect in EA and non-EA was the same for the remaining two index SNPs; rs3764261 was significantly weaker for HDL in AA, and rs28927680 was significantly stronger for TG in Pacific Islanders. There were no observations of opposite effects where both the EA effect and the non-EA effect were significant.

Considering only the 15 SNPs with a significantly inconsistent effect between EA and at least one non-EA population, 14 of 15 diluted toward the null (p<0.01, Table 2), a trend driven by the AA population, where all 12 out of 12 significant inconsistencies were diluted. Expanding analysis to all 51 loci replicated in EA, regardless of whether a significant difference was observed between EA and non-EA at a given SNP, we observed a significant excess of effects diluted toward the null (ßpopEA<1) in AA, HA, and NA populations (Table S5). Comparisons between non-EA populations revealed that diluted effect sizes were significantly more likely in AA than in any other non-EA population.

Given that differential effect sizes were observed for many tagSNPs, we sought to leverage the data in order to assess the relative contributions of several factors that might contribute to the significant trend toward diluted effects, including gene–environment interaction with an exposure that varies across populations (differential environment), differences in the correlation between the index SNP and the functional variant across populations (differential tagging), modulation of the index SNP effect by additional, population-specific polymorphism (differential genetic background), population-specific synthetic alleles (combinations of rare, functional alleles tagged by a single common tagSNP [20]), or some combination of these factors. It seems unlikely that differential environments would be much more frequent in AA than other non-EA populations, or that differential environment would consistently bias toward the null within AA. Differential tagging is consistent with differentially diluted effects in AA; because linkage disequilibrium extends over significantly shorter distances in African populations than in non-African populations [21],[22], common functional variants (or synthetic alleles) are likely to be less strongly tagged by the index tagSNPs in AA. Differential genetic background effects in AA would also be consistent with the high nucleotide diversity known to exist in this population. The rare functional variants contributing to synthetic alleles will tend to be younger than common variants, and therefore are more likely to be population-specific, so synthetic alleles are compatible with the trend toward dilution. Thus, although differential environmental effects cannot be excluded, the observed data are more consistent with differential tagging and/or differential genetic background effects, and synthetic alleles cannot be excluded.

Genetic background effects can be subdivided into modifying effects, where variants elsewhere in the genome directly alter the effect associated with a given index SNP, and interference effects, where secondary variants change the proportion of variance explained by the index SNP. Interfering functional variants with effects in the same direction as the index SNP would tend to dilute the apparent effect size at the index variant. The most likely source of such variants is the region surrounding an index SNP, as demonstrably functional variants already exist in that region. Although examples have been described of genes carrying both risk and protective mutations [23][25], others clearly exhibit trends toward risk alleles with similar effects (e.g., preferentially toward breast cancer risk alleles at BRCA1 [26]). If the direction of effect for functional variants in a given region is consistently biased, then an increase in the number of interfering variants within a given population would be consistent with a trend toward dilution of index effects. The higher nucleotide diversity observed in African populations relative to non-African populations [27],[28] would be consistent with a greater burden of secondary functional variants in AA than other populations.

In order to assess contribution of the factors outlined above to differential effect sizes between EA and AA in the index tagSNP associations, high density genotype data were collected from a subset of the PAGE African American sample. The number of AA individuals used for index tagSNP analyses varied by phenotype, with an average of 7501 (Table S3). Similar data on other populations are currently unavailable, so only loci showing differential effects between EA and AA could be analyzed. Genotype data were collected using the Metabochip, a high density genotyping array commercially available from Illumina. Detailed methods for the Metabochip genotype data collection, calling, and quality control are available elsewhere [12].

In order to measure the contribution of differential LD to dilution, we need a model of how changes in LD between tagSNP and a functional variant would be expected to alter the observed effect size at the tagSNP, assuming that the effect size at the functional variant is the same in both populations. Given a functional SNP (fSNP) and an associated tagSNP, linkage disequilibrium between the two SNPs can be described as the measurement error introduced by genotyping the tagSNP, rather than genotyping the fSNP directly. As such, by appealing to prior work on regression dilution bias, it can be shown that the effect size β′ at the tagSNP is related to the effect size β at the fSNP by the following equation: An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e016.jpg (see Text S1 for details). Thus, assuming that the effect size at the fSNP is constant between populations, when linkage disequilibrium between tagSNP and fSNP is weaker in a given population, we expect to see a greater degree of dilution bias for the estimated tagSNP effect size. Rearranging this equation, An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e017.jpg. Extrapolating to compare the degree of dilution bias between AA and EA populations, we expect changes in linkage disequilibrium across populations to be reflected by changes in relative effect size:

equation image

Assuming the effect size of the functional variant is the same in both populations, this reduces to:

equation image

The above equation allows us to directly compare the observed distribution of relative effect sizes at the tagSNPs in AA and EA (An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e020.jpg) against the relative strength of tagging in AA and EA (An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e021.jpg). Considering the subset of index tagSNPs in regions that were present on the Metabochip, we observed 51 index tagSNPs that fell into 47 independent loci on the Metabochip. We identified the set of SNPs tagged by each index tagSNP at r2>0 .8 in an EA population [29],[30], yielding a total of 1,093 tagged SNPs for the 51 index tagSNPs. For each of these 1,144 SNPs, we then calculated An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e022.jpg. Let this represent the expected distribution of differential LD between AA and EA. Next, we calculated An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e023.jpg for the subset of 40 of the 51 index tagSNPs that replicated at q = 0.05 in EA, truncating at 0 if the signs were opposite between populations. These two distributions (An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e024.jpg in all 1,144 SNPs versus An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e025.jpg for the 40 index tagSNPs) were not significantly different by two-tailed t test. Thus, we cannot reject the hypothesis that the observed dilution bias in AA effect sizes at the index tagSNPs is consistent with the observed distribution of differential LD between the two populations. A single-locus example of the potential for differential LD to contribute to diluted effect sizes is shown in Figure 2.

Figure 2
Dilution of effect size at PSRC1 for LDL.

Considering the 12 SNPs showing differential effect size in AA, regions spanning 11 were present on the Metabochip (Table S3). Before comparison with EA, we compared the observed effect sizes at the index tagSNPs in the full AA sample and the subsample of AAs genotyped on the Metabochip (AAmchip). Three of the index tagSNPs failed to genotype on the Metabochip, leaving eight index tagSNPs for this direct comparison (Table S4). No significant allele frequency differences were observed between the AAmchip subset and the full AA population, consistent with AAmchip being a representative subsample. A significantly inconsistent and diluted effect size in AAmchip compared to EA was still observed for five of these eight tagSNPs (p<0.05, Table S4). The index tagSNPs without a significant difference likely reflect reduced power to detect the differential effect size in the AAmchip subsample, as these three index tagSNPs also had the least significant differential effect when comparing the full PAGE AA subpopulation against EA.

The Metabochip genotype data allowed us to evaluate regions spanning each of the 11 variants for the underlying contributions of population-specific alleles, differential tagging, and secondary alleles to differential effect sizes. Detailed discussion of each locus is provided in Text S1. In summary, the 11 SNPs fell in 10 Metabochip regions, so all SNPs in each of the 10 regions were assessed for association with the reported trait in AAmchip. The threshold level for significance within each region was conservatively adjusted for multiple testing by Bonferroni adjustment for the number of SNPs successfully genotyped on the Metabochip within the region, with minor allele frequency greater than 1% in the AAmchip sample. For example, the Metabochip region spanning CETP contained 84 SNPs, so our significance threshold for that region was p<0.05/84 = 1.1*10−4. One locus (APOE) could not be dissected confidently as LD data for the index tagSNP were not available in EA, and two loci were underpowered to draw strong conclusions, as evidenced by the failure of any variant in the region to show a significantly inconsistent effect with the index tagSNP effect in EA. Among the remaining seven loci, we observed one clear example of a diluted signal consistent with EA-specific functional alleles, either common or synthetic (Figure 3a), and five loci showed patterns consistent with fine-mapping of the index tagSNP bin (Figure 3d–f, Figure 4a, 4d). One of these fine-mapped the EA association to a variant that was not strongly associated with the index tagSNP in EA (r2<0.5, Figure 3f), potentially consistent with a synthetic allele in EA. We also observed statistically significant secondary functional alleles at three loci (Figure 4).

Figure 3
Examples of loci without evidence of association in AAmchip or fine mapping EA signal.
Figure 4
Examples of secondary alleles in the AA population.

Thus, although the overall pattern of effect dilution in AA is consistent with expectations on the basis of differential LD patterns between AA and EA populations, putative examples of EA-specific alleles and secondary alleles in AA were also observed. A contribution from synthetic alleles cannot be excluded, and may well account for the EA-specific allele at CILP2 (Figure 3a). However, at half of the 10 loci we observed at least one of the tagged SNPs in EA that showed an effect size in AA consistent with the effect size at the tagSNP in EA. These examples of fine-mapping EA signal suggest that at least half of EA GWAS signals tag a common, functional variant. The observed excess of dilution effects in AA (as compared to other non-EA populations) suggests that African-descended populations will be the most useful single subpopulation for fine-mapping of EA GWAS associations, although the significant trend toward excess dilution in HA and NA populations (Table S5) suggests that trans-ethnic fine-mapping may prove more powerful than fine-mapping with any single non-EA population.

In conclusion, we have assessed the generalization of GWAS associations from EA populations across five clinically relevant traits, in five non-EA populations. Our results demonstrate that although most EA GWAS findings can be expected to show an effect in the same direction for non-EA populations, a significant fraction of GWAS-identified variants from EA will exhibit differential effect sizes in at least one non-EA population, and these differential results will be far more frequent in the AA population. These findings suggest that expanded GWAS and fine-mapping efforts focused on non-EA populations, especially AA, will substantially enhance our understanding of the genetic architecture of common traits within non-EA populations. It will be particularly important to extend GWAS discovery efforts to non-EA populations if genetic risk prediction models using tagSNP genotypes demonstrate clinical utility, because risk estimates derived from European GWAS clearly generalize imperfectly to non-EA populations. Our analyses suggest that variable LD in its many guises accounts for much of the heterogeneity of effect size at index tagSNPs, rather than any “true” differences in effect size between populations for the functional variants that were tagged. Thus, risk models derived directly from genotypes at functional variants (rather than tagSNPs) may generalize more effectively to non-EA populations.

Materials and Methods

Selecting Index SNPs from Prior Reports

Traits considered were those for which more than 10 GWAS-identified variants were genotyped in the first year of PAGE. Variants considered for this analysis included 13 previously reported to associate with body mass index, 20 for type 2 diabetes, 27 for HDL, 19 for LDL, and 14 for triglycerides. Eleven of these GWAS-identified variants were previously reported to associate with more than one trait in EA (Table S1), so we constrained the analysis of each such SNP to whichever trait had the most significant association (smallest p value) in the PAGE EA population, leaving a panel of 82 unique variants.

Because highly correlated SNPs might overweight specific results toward a specific trait or gene, we extracted a subset of minimally correlated index GWAS-identified variants from this panel of 82. At each step, we added the SNP with the most significant association in PAGE EA to a list of index SNPs, and then filtered the remaining SNPs not yet in the index list to exclude those exceeding r2 = 0.2 in the PAGE EA population with any index tagSNP. The panel of 82 SNPs was recursively filtered in this manner, leaving a final panel of 69 index SNPs, each of which was minimally associated with any other index SNP in the PAGE EA (r2<0.2). One additional SNP (rs11084753) was removed from the analysis due to concerns regarding power to replicate, leaving a final panel of 68 index SNPs for analysis, including seven index SNPs for BMI, 18 for HDL, 15 for LDL, nine for triglycerides, and 19 for T2D (Table S2).

Assessing Power to Replicate Previous Reports

Power estimates are taken directly from Fesinmeyer et al. [13] for BMI, Dumitrescu et al. [15] for lipids, and Haiman et al. [14] for T2D. For details, see the original publications.

Defining Generalization results for Each SNP

In order to assess the generalization of effects to each population, we used effect sizes (β) and standard errors derived from minimally adjusted (age, sex, and study), ancestry-specific meta-analyses described in the primary PAGE publications [15],[13],[14]. Using these data, we tested two hypotheses: first, that the GWAS-identified variant has no effect in the non-EA population (i.e., the coefficient ßpop = 0 in a linear or logistic regression model), and second, that the effect size in the non-EA population is the same as the effect size in EA (An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e026.jpg). The first hypothesis was tested by assuming the estimate An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e027.jpg is normally distributed (which is reasonable as sample sizes exceeded 500 for all populations) and calculating the probability that An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e028.jpg given An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e029.jpg and the standard error of An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e030.jpg. The second hypothesis was tested by defining An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e031.jpg and calculating the probability that An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e032.jpg, again assuming An external file that holds a picture, illustration, etc.
Object name is pbio.1001661.e033.jpg to be normally distributed. These tests are all carried out at a nominal significance level of 0.05, as we see them as the (single) test that an investigator may carry out to validate a result first observed in EA in another ethnic group, and then significance was assigned using the Benjamini-Hochberg method at a false discovery rate of 5%.

A reasonable concern in these analyses is that population stratification can distort effect size estimates in some circumstances. Some of the effect sizes from trait-specific PAGE manuscripts were not adjusted for genetic ancestry, due to either availability of data [15] or informed consent in specific populations [13],[14]; where available and allowed we have used the ancestry adjusted effect sizes. Both ancestry adjusted and unadjusted data were available for the PAGE obesity analysis [13], where ancestry adjustment did not significantly alter effect size estimates.

Supporting Information

Figure S1

Generalization analysis in the PAGE populations. We plot the ratio of ßpopEA on the y-axis as an indicator of both consistency of direction (positive values are consistent with effects in the same direction) and relative magnitude of effect (consistent but weaker effects in the non-EA will have ratios between 0 and 1). We plot the coded allele frequency (CAF) in the non-EA population on the x-axis, as a proxy for power to replicate in the available sample size. Data points are colored as follows: ambiguous SNPs are light blue (ßpop = 0 and ßpop = ßEA), strictly generalized SNPs are dark blue (ßpop≠0 and ßpop = ßEA), differentially generalized SNPs are dark red (ßpop≠0 and ßpop≠ßEA), and differential SNPs are pink (ßpop = 0 and ßpop≠ßEA). The y-axis has been constrained to (−4.4) for illustrative purposes; some loci yielded ßpopEA ratios outside this range, but pEA>0.05 for all of these. As expected, larger non-EA populations show less scatter in ßpopEA than the smaller non-EA populations (particularly Pacific Islanders), consistent with more precise estimates of ßpop in the larger non-EA populations. A clear trend is observed toward ßpopEA ratios greater than zero in all populations, reflecting consistency of direction between EA and non-EA populations. An additional trend toward ratios greater than zero but less than one is observed in African Americans, representing the trend toward dilution in this population, relative to EA. No such trend is apparent in the other non-EA populations. Differential and ambiguous SNPs are observed throughout the CAF range, consistent with the assertion that these categories do not reflect a systematic bias toward underpowered, low-frequency variants.


Figure S2

Multiple associations at CETP. Two index tagSNPs with differential effect size were observed at CETP: rs3764261 and rs9989419. (a) CETP regional LocusZoom plot with LD in EA samples color coded relative to rs3764261. Although this index tagSNP exhibited differential effect (purple point indicated with arrow), several of the tagged SNPs (red data points near −log10(p) = 34, rs247616, rs247617, rs183130) exhibited effect sizes consistent with fine-mapping of this association (r2 = 0.99 in EA, r2 = 0.74 in AA). Interestingly, the strongest effect observed in the region was at a SNP uncorrelated with rs3764261 (rs17231520, r2EA<0.001). (b) Same plot, but adjusting genotype at rs274616 (the best signal from a tagged SNP in the EA rs3764261 bin). The signal from tagged SNPs has clearly been reduced to background levels, and residual signal is clearly visible for untagged SNPs. (c) Same plot, but now adjusting for genotype at rs274616 and rs193695 (the best signal from a tagged SNP in the EA rs9989419 bin). Again, significant residual signal is observed. Figures S2d–f show the same data, but with LD in EA samples color coded relative to rs9989419. Although rs9989419 failed to genotype on the Metabochip, a strongly tagged SNP is visible in (d). Although this variant was weakly tagged by rs3764261 (compare panel a with d), the association signal does not appear to be independent of rs3764261, as residual association is not significant at this variant after adjustment for rs247616 (e). Thus, there is clearly residual association at this locus after adjusting for both of the strongest EAtaggedSNPs, consistent with either additional functional variation in the region, differential tagging, or differential synthetic alleles at this locus.


Table S1

Categorization of tagSNPs. TagSNPs are categorized on the basis of the primary phenotype in the original GWAS report, and by whether these SNPs were categorized as index SNPs or not in the present analysis.


Table S2

Raw meta-analysis data for each of the index SNPs is given in each PAGE subpopulation, including the number of individuals in the subpopulation, the observed effect size and s.e., the p value associated with testing the hypothesis that the observed effect was significant within the subpopulation (p_Beta_<subpopulation>_ = _0), and the p value associated with testing the hypothesis of equal effect size in European and non-European populations (p_Beta_<subpopulation>_equal_Beta_Eur.Am.), as well as the generalization category for that SNP in the subpopulation (<subpopulation>_snpcat).


Table S3

Details of the effect sizes and generalization results for the 12 tagSNPs with inconsistent effect size observed in the EA and AA PAGE populations, extracted from Table S2.


Table S4

Summary of Metabochip variants that passed QC with minor allele frequency greater than 1% in the PAGE AAmchip subpopulation is given, along with the generalization results comparing only the index tagSNP from the AAmchip against PAGE EA.


Table S5

Summary of dilution effects in subpopulations. Data are shown first for the subset of 51 tagSNPs replicated in EA, then for all 68 tagSNPs in Table S2. Within each subpopulation, we first assessed the probability of the observed frequency of dilution (ßpopEA<1) against the null hypothesis of no such trend, by chi-square test. Then we compared the frequency of dilution effects between other subpopulations and the AA subpopulation. As noted, we observed a significant excess of effects diluted toward the null (ßpopEA<1) in AA, HA, and NA populations, and this trend toward diluted effects was significantly stronger in the AA subpopulation than in any other subpopulation.


Text S1

Supplemental methods.



The data and materials included in this report result from a collaboration between the following studies: The Epidemiologic Architecture for Genes Linked to Environment (EAGLE), the Multiethnic Cohort study (MEC), the Epidemiology of putative genetic variants: The Women's Health Initiative (WHI), the Genetic Epidemiology of Causal Variants Across the Life Course (CALiCo), and the PAGE Coordinating Center. The PAGE consortium thanks the staff and participants of all PAGE studies for their important contributions, as well as Dr. Teri Manolio for her support and leadership within the consortium. Active PAGE investigators at the time of this analysis included: Coordinating Center: Rutgers University, Piscataway, NJ: Tara Matise, Steve Buyske, Julia Higashio, Andrew Nato; University of Southern California, Los Angeles, CA: Jose Luis Ambite, Ewa Deelman. NHGRI: Division of Genomic Medicine NHGRI, NIH, Bethesda, MD: Teri Manolio, Lucia Hindorff, Heather Junkins, and Erin Ramos. CALiCo: University of North Carolina, Chapel Hill, NC: Kari E. North, Gerardo Heiss, Kira Taylor, Nora Franceschini, Christy Avery, Misa Graff, Danyu Lin, Miguel Quibrera; Baylor College of Medicine, Houston, TX: Barbara Cochran; Johns Hopkins Bloomberg School of Public Health, Baltimore, MD: Linda Kao; Penn Medical Lab, Washington DC: Jason Umans; SW Foundation for BioMedical Research, San Antonio, TX: Shelley Cole, Jean MacCluer; University of Alabama at Birmingham, Birmingham, AB: Sharina Person; University of Minnesota, Minneapolis, MN: James Pankow, Myron Gross; University of Texas Health Science Center, Houston: Eric Boerwinkle, Myriam Fornage; University of Vermont, Burlington, VT: Peter Durda, Nancy Jenny; University of Washington, Seattle, WA: Bruce Patsy, Alice Arnold, Petra Buzkova. EAGLE: Vanderbilt University, Nashville, TN: Dana Crawford, Jonathan Haines, Deborah Murdock, Kim Glenn, Kristin Brown-Gentry, Tricia Thornton-Wells, Logan Dumitrescu, Janina Jeff, William S. Bush, Sabrina L. Mitchell, Robert Goodloe, Sarah Wilson, Jonathan Boston, Jennifer Malinowski, Nicole Restrepo, Matthew Oetjens, Jay Fowke, Wei Zheng; Heidelberg University, Tiffin, OH: Kylee Spencer; Pennsylvania State University, State College, PA: Marylyn Ritchie, Sarah Pendergrass. MEC: University of Hawaii, Honolulu, HI: Loïc Le Marchand, Lynne Wilkens, Lani Park, Maarit Tiirikainen, Laurence Kolonel, Unhee Lim, Iona Cheng, Hansong Wang, Ralph Shohet; Keck School of Medicine, University of Southern California, Los Angeles, CA: Christopher Haiman, Daniel Stram, Brian Henderson, Kristine Monroe, Fredrick Schumacher. WHI: Fred Hutchinson Cancer Research Institute (FHCRC), Seattle, WA: Charles Kooperberg, Ulrike Peters, Garnet Anderson, Chris Carlson, Ross Prentice, Andrea LaCroix, Chunyuan Wu, Cara Carty, Jian Gong, Stephanie Rosse, Alicia Young, Jeff Haessler, Jonathan Kocarnik, Megan Fesinmeyer, Yi Lin; Ohio State Medical Center, Columbus, OH: Rebecca Jackson; Translational Genomic Science Institute (TGen): David Duggan; University of Pittsburgh, Pittsburgh, PA: Lew Kuller. The authors thank the WHI investigators and staff for their dedication, and the study participants for making the program possible. EAGLE would like to thank Dr. Geraldine McQuillan and Jody McLean for their help in accessing the Genetic NHANES data. The Vanderbilt University Center for Human Genetics Research, Computational Genomics Core provided computational and/or analytical support for this work. The EAGLE/NHANES DNA samples are stored and plated by the Vanderbilt DNA Resources Core. Genotyping was performed by Ping Mayo, Melissa Allen, and Dr. Nathalie Schnetz-Boutaud in the laboratory of Dr. Jonathan Haines and Hailing Jin and Nila Gillani under the direction of Dr. Holli Dilks in the Vanderbilt DNA Resources Core.


African Americans
Apolipoprotein AI, CIII, AIV, AV gene cluster
Apolipoprotein E, CI, CIV gene cluster
Asian Americans
body mass index
causal variants across the life course
cadherin, EGF LAG seven-pass G-type receptor 2
cholesteryl ester transfer protein, plasma
cartilage intermediate layer protein 2
European ancestry
Epidemiologic Architecture for Genes Linked to Environment
false discovery rate
fat mass and obesity associated
glucokinase (hexokinase 4) regulator
genome-wide association study
Hispanic Americans
high density lipoprotein
potassium inwardly-rectifying channel, subfamily J, member 11
linkage disequilibrium
low density lipoprotein
lipase, hepatic
multi-ethnic cohort
MLX interacting protein-like
Native Americans
Population Architecture using Genomics and Epidemiology
pre-B-cell leukemia homeobox 4
Pacific Islanders
proline/serine-rich coiled-coil 1
simple nucleotide polymorphism
Type 2 Diabetes
tribbles homolog 1
Women's Health Initiative

Funding Statement

The Population Architecture Using Genomics and Epidemiology (PAGE) program is funded by the National Human Genome Research Institute (NHGRI). The data and materials included in this report result from a collaboration between the following studies: The “Epidemiologic Architecture for Genes Linked to Environment (EAGLE)” is funded through the NHGRI PAGE program (U01HG004798-01). Genotyping services for select NHANES III SNPs presented here were also provided by the Johns Hopkins University under federal contract number (N01-HV-48195) from NHLBI. The study participants derive from the National Health and Nutrition Examination Surveys (NHANES), and these studies are supported by the Centers for Disease Control and Prevention. The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the Centers for Disease Control and Prevention. The Multiethnic Cohort study (MEC) characterization of epidemiological architecture is funded through the NHGRI PAGE program (U01HG004802). The MEC study is funded through the National Cancer Institute (R37CA54281, R01 CA63, P01CA33619, U01CA136792, and U01CA98758). Funding support for the “Epidemiology of putative genetic variants: The Women's Health Initiative” study is provided through the NHGRI PAGE program (U01HG004790). The WHI program is funded by the National Heart, Lung, and Blood Institute; NIH; and US Department of Health and Human Services through contracts N01WH22110, 24152, 32100-2, 32105-6, 32108-9, 32111-13, 32115, 32118-32119, 32122, 42107-26, 42129-32, and 44221. Funding support for the Genetic Epidemiology of Causal Variants Across the Life Course (CALiCo) program was provided through the NHGRI PAGE program (U01HG004803). The following studies contributed to this manuscript and are funded by the following agencies: The Atherosclerosis Risk in Communities (ARIC) Study is carried out as a collaborative study supported by National Heart, Lung, and Blood Institute contracts N01-HC-55015, N01-HC-55016, N01-HC-55018, N01-HC-55019, N01-HC-55020, N01-HC-55021, N01-HC-55022. The Coronary Artery Risk Development in Young Adults (CARDIA) study is supported by the following National Institutes of Health, National Heart, Lung and Blood Institute contracts: N01-HC-95095; N01-HC-48047; N01-HC-48048; N01-HC-48049; N01-HC-48050; N01-HC-45134; N01-HC-05187; and N01-HC-45205. The Cardiovascular Health Study (CHS) is supported by NHLBI contracts HHSN268201200036C, N01-HC-85239, N01-HC-85079 through N01-HC-85086; N01-HC-35129, N01 HC-15103, N01 HC-55222, N01-HC-75150, N01-HC-45133 and NHLBI grant HL080295, with additional contribution from NINDS. Additional support was provided through AG-023629, AG-15928, AG-20098, and AG-027058 from the NIA. See also CHS GWAS DNA handling and genotyping was supported in part by National Center for Research Resources grant M01-RR00425 to the Cedars-Sinai General Clinical Research Center Genotyping core and National Institute of Diabetes and Digestive and Kidney Diseases grant DK063491 to the Southern California Diabetes Endocrinology Research Center. The Strong Heart Study (SHS) is supported by NHLBI grants U01 HL65520, U01 HL41642, U01 HL41652, U01 HL41654, and U01 HL65521. The opinions expressed in this paper are those of the author(s) and do not necessarily reflect the views of the Indian Health Service. Assistance with phenotype harmonization, SNP selection and annotation, data cleaning, data management, integration and dissemination, and general study coordination was provided by the PAGE Coordinating Center (U01HG004801-01). The National Institutes of Mental Health also contributes to the support for the Coordinating Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Need AC, Goldstein DB (2009) Next generation disparities in human genomics: concerns and remedies. Trends Genet 25 (11) 489–494. [PubMed]
2. Hindorff LA, et al. (2009) Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A 106 (23) 9362–9367. [PubMed]
3. Lu X, et al. (2012) Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease. Nat Genet 44 (8) 890–894. [PubMed]
4. Okada Y, et al. (2012) Common variants at CDKAL1 and KLF9 are associated with body mass index in east Asian populations. Nat Genet 44 (3) 302–306. [PubMed]
5. Prasad P, et al. (2012) Caucasian and Asian specific rheumatoid arthritis risk loci reveal limited replication and apparent allelic heterogeneity in north Indians. PLoS ONE 7 (2) e31584.doi:10.1371/journal.pone.0031584. [PMC free article] [PubMed]
6. Takeuchi F, et al. (2012) Genome-wide association study of coronary artery disease in the Japanese. Eur J Hum Genet 20 (3) 333–340. [PMC free article] [PubMed]
7. Chen R, et al. (2012) Type 2 diabetes risk alleles demonstrate extreme directional differentiation among human populations, compared to other diseases. PLoS Genet 8 (4) e1002621.doi:10.1371/journal.pgen.1002621. [PMC free article] [PubMed]
8. Dorajoo R, et al. (2012) Replication of 13 obesity loci among Singaporean Chinese, Malay and Asian-Indian populations. Int J Obes (Lond) 36 (1) 159–163. [PubMed]
9. Ntzani EE, et al. (2012) Consistency of genome-wide associations across major ancestral groups. Hum Genet 131 (7) 1057–1071. [PubMed]
10. Sharma M, et al. (2012) Large-scale replication and heterogeneity in Parkinson disease genetic loci. Neurology 79 (7) 659–667. [PMC free article] [PubMed]
11. Matise TC, et al. (2011) The Next PAGE in understanding complex traits: design for the analysis of Population Architecture Using Genetics and Epidemiology (PAGE) Study. Am J Epidemiol 174 (7) 849–859. [PMC free article] [PubMed]
12. Buyske S, et al. (2012) Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study. PLoS ONE 7 (4) e35651.doi:10.1371/journal.pone.0035651. [PMC free article] [PubMed]
13. Fesinmeyer MD, et al. (2012) Genetic risk factors for BMI and obesity in an ethnically diverse population: Results from the Population Architecture Using Genomics and Epidemiology (PAGE) Study. Obesity (Silver Spring) [PubMed]
14. Haiman CA, et al. (2012) Consistent directions of effect for established type 2 diabetes risk variants across populations: the Population Architecture using Genomics and Epidemiology (PAGE) Consortium. Diabetes 61 (6) 1642–1647. [PMC free article] [PubMed]
15. Dumitrescu L, et al. (2011) Genetic determinants of lipid traits in diverse populations from the population architecture using genomics and epidemiology (PAGE) study. PLoS Genet 7 (6) e1002138.doi:10.1371/journal.pgen.1002138. [PMC free article] [PubMed]
16. Wu Y, Waite LL, Jackson AU, Sheu WH-H, Buyske S, et al. . (2013) Trans-ethnic fine-mapping of lipid loci identifies population-specific signals and allelic heterogeneity that increases the trait variance explained. PLoS Genet 9 ((3) ) e1003379.10.1371/journal.pgen.1003379. [PMC free article] [PubMed]
17. Hochberg Y, Benjamini Y (1990) More powerful procedures for multiple significance testing. Stat Med 9 (7) 811–818. [PubMed]
18. Ioannidis JP, Ntzani EE, Trikalinos TA (2004) “Racial” differences in genetic effects for complex diseases. Nat Genet 36 (12) 1312–1318. [PubMed]
19. Kiryluk K, et al. (2012) Geographic differences in genetic susceptibility to IgA nephropathy: GWAS replication study and geospatial risk analysis. PLoS Genet 8 (6) e1002765.doi:10.1371/journal.pgen.1002765. [PMC free article] [PubMed]
20. Dickson SP, et al. (2010) Rare variants create synthetic genome-wide associations. PLoS Biol 8 (1) e1000294.doi:10.1371/journal.pbio.1000294. [PMC free article] [PubMed]
21. Reich DE, et al. (2001) Linkage disequilibrium in the human genome. Nature 411 (6834) 199–204. [PubMed]
22. Carlson CS, et al. (2004) Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am J Hum Genet 74 (1) 106–120. [PubMed]
23. Abifadel M, et al. (2003) Mutations in PCSK9 cause autosomal dominant hypercholesterolemia. Nat Genet 34 (2) 154–156. [PubMed]
24. Kotowski IK, et al. (2006) A spectrum of PCSK9 alleles contributes to plasma levels of low-density lipoprotein cholesterol. Am J Hum Genet 78 (3) 410–422. [PubMed]
25. Cohen J, et al. (2005) Low LDL cholesterol in individuals of African descent resulting from frequent nonsense mutations in PCSK9. Nat Genet 37 (2) 161–165. [PubMed]
26. Hogervorst FB, et al. (1995) Rapid detection of BRCA1 mutations by the protein truncation test. Nat Genet 10 (2) 208–212. [PubMed]
27. Crawford DC, et al. (2004) Haplotype diversity across 100 candidate genes for inflammation, lipid metabolism, and blood pressure regulation in two populations. Am J Hum Genet 74 (4) 610–622. [PubMed]
28. The 1000 Genomes Project Consortium (2010) A map of human genome variation from population-scale sequencing. Nature 467 (7319) 1061–1073. [PMC free article] [PubMed]
29. Voight BF, et al. (2012) The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits. PLoS Genet 8 (8) e1002793.doi:10.1371/journal.pgen.1002793. [PMC free article] [PubMed]
30. Berglund G, et al. (1993) The Malmo Diet and Cancer Study. Design and feasibility. J Intern Med 233 (1) 45–51. [PubMed]
31. Teslovich TM, et al. (2010) Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466 (7307) 707–713. [PMC free article] [PubMed]

Articles from PLoS Biology are provided here courtesy of Public Library of Science