PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of diabetesSubscribeSearchDiabetes JournalAmerican Diabetes Association
 
Diabetes. 2015 July; 64(7): 2646–2657.
Published online 2015 February 9. doi:  10.2337/db14-1715
PMCID: PMC4477349

Role of Established Type 2 Diabetes–Susceptibility Genetic Variants in a High Prevalence American Indian Population

Abstract

Several single nucleotide polymorphisms (SNPs) associated with type 2 diabetes mellitus (T2DM) have been identified, but there is little information on their role in populations at high risk for T2DM. We genotyped SNPs at 63 T2DM loci in 3,421 individuals from a high-risk American Indian population. Nominally significant (P < 0.05) associations were observed at nine SNPs in a direction consistent with the established association. A genetic risk score derived from all loci was strongly associated with T2DM (odds ratio 1.05 per risk allele, P = 6.2 × 10−6) and, in 292 nondiabetic individuals, with lower insulin secretion (by 4% per copy, P = 4.1 × 10−6). Genetic distances between American Indians and HapMap populations at T2DM markers did not differ significantly from genomic expectations. Analysis of U.S. national survey data suggested that 66% of the difference in T2DM prevalence between African Americans and European Americans, but none of the difference between American Indians and European Americans, was attributable to allele frequency differences at these loci. These analyses suggest that, in general, established T2DM loci influence T2DM in American Indians and that risk is mediated in part through an effect on insulin secretion. However, differences in allele frequencies do not account for the high population prevalence of T2DM.

Introduction

In recent years, more than 70 distinct genomic regions have been identified in which single nucleotide polymorphism (SNP) markers show reproducible association with type 2 diabetes mellitus (T2DM) at genome-wide statistical significance (P < 5 × 10−8) (117). Most of these variants were discovered by genome-wide association studies (GWAS) in European populations, and their effects are best characterized in populations of European ancestry. Studies in other ethnic groups suggest that effects on T2DM are similar to those seen in Europeans for most variants (18,19), but clear examples of heterogeneity in effects have been observed (8,10,20). There is limited information on the role of these established variants in populations at high risk for T2DM or on the extent to which differences in allele frequencies at these variants account for differences in population risk. In the present study, we analyze 63 established T2DM-susceptibility variants in Pima Indians, an American Indian population in whom the prevalence of T2DM is extraordinarily high (21).

Research Design and Methods

Participants

Subjects were participants in a longitudinal study conducted in the Gila River Indian Community in central Arizona, where most residents are Pima Indians (21). The present study consisted of 3,421 individuals whose self-reported heritage was full Pima, Tohono O’odham, or a mixture of these closely related tribes and who had DNA available. These individuals constituted 1,951 sibships. There were 1,964 women and 1,457 men; mean ± SD age at last examination was 40.6 ± 16.5 years. Height and weight were measured, and a 75-g oral glucose tolerance test was administered; diabetes was diagnosed in 1,615 individuals (47.2%) according to 1997 American Diabetes Association criteria (22), i.e., 2-h postload plasma glucose ≥11.1 mmol/L, fasting plasma glucose ≥7.0 mmol/L, or a diagnosis during routine clinical care.

A subset of individuals participated in detailed physiologic studies to assess metabolic predictors of T2DM. Body composition was measured by hydrodensitometry or by DEXA, as previously described (23), in 405 nondiabetic full-heritage Pimas (172 women and 233 men; mean ± SD age 26.7 ± 6.1 years). Insulin sensitivity was measured in these 405 individuals by the hyperinsulinemic-euglycemic clamp (23). Insulin was infused at physiologic levels (~130 µmol/L), and glucose was infused to maintain euglycemia. Rate of glucose uptake, normalized to estimated metabolic body size (EMBS), was taken as a measure of insulin sensitivity (M) (milligrams per kilogram EMBS per minute). Insulin secretion was measured as the acute insulin response (microunits per milliliter) 3–5 min after a 25-g intravenous glucose challenge (23) in 292 individuals (105 women and 187 men; mean ± SD age 26.7 ± 6.1 years) with normal glucose tolerance (2-hour postload glucose <7.8 mmol/L).

Genotyping

A sentinel SNP for each region was selected for genotyping from previously reported GWAS (117). Two SNPs were selected for KCNQ1 and CDC123, where two distinct sets of variants have been described. In addition, 45 ancestry-informative markers (24) were genotyped for estimation of the individual proportion of European heritage (25). Genotyping was conducted by the SNPplex method (Life Technologies, Carlsbad, CA) or the BeadXpress system (Illumina, San Diego, CA) according to the manufacturer’s instructions. Results for 18 SNPs were reported previously (20,2629). They are included here for a more complete characterization of the effects of established T2DM loci.

Association with T2DM and Related Traits

Association between genotype and T2DM at the last research examination was analyzed by a logistic regression model, which was fit by the generalized estimating equation procedure to account for sibship. Genotype was coded as a numeric variable representing number of risk alleles as defined in previous GWAS. Thus, an odds ratio (OR) >1 indicates association in the same direction as the established association and an OR <1 indicates association in the opposite direction. Continuous variables were analyzed with a linear mixed model in which genotype and other covariates were fixed effects and sibship was a random effect. The logarithm of each variable was analyzed, and the regression coefficient was exponentiated to obtain the effect per copy of the T2DM risk allele, expressed as a multiplier.

For assessment of whether associations in Pimas were consistent with those in Europeans, ORs were compared by the Cochran Q test of homogeneity, and heterogeneity was quantified by the I2 measure (30). ORs for Europeans were taken from previous publications (1,2,8,9,1417,3140). For assessment of whether GWAS-defined risk alleles contribute in aggregate to T2DM in Pimas, a multiallelic genetic risk score (GRS) was created by summing the number of risk alleles over all loci. To avoid reduction in sample size resulting from missing data at a few loci, we calculated the probability that an individual was of each possible genotype for each missing value from the genotypes in the individual’s relatives using MLINK (41); these probabilities were used in calculating the GRS.

To test for heterogeneity across all loci, we combined P values derived from the heterogeneity test for individual SNPs by constructing a signed Z score. The Z score was computed for each SNP as Zi = sign[ln(OREUi) − ln(ORPIi)]Φ−1(Pheti/2), where OREUi represents the OR for the ith SNP in Europeans, ORPIi represents the corresponding OR in Pimas, Pheti is the P value for heterogeneity, and Φ−1 represents the inverse of the cumulative normal probability function. The sum of the Z scores across all SNPs divided by the square root of the number of SNPs (Z*) was used to calculate a P value for the null hypothesis of homogeneity across all markers (42). If Z* is negative, it indicates that ORs on average are weaker in Pimas than in Europeans, whereas if Z* is positive it indicates that ORs are stronger in Pimas.

Differences in Allele Frequencies

Frequency of the risk allele was estimated by maximum likelihood methods using the ILINK program to account for family membership (41). Data for these 63 SNPs were obtained from the International HapMap Project (http://hapmap.ncbi.nlm.nih.gov/) or, if not available from HapMap, from the 1000 Genomes Project (http://www.1000genomes.org/). For comparison of allele frequencies with other major continental ethnic groups, data were obtained for individuals of European ancestry from Centre d’Etude du Polymorphisme Humain families in Utah (CEU), for East Asians from Han Chinese in Beijing (CHB), and for Africans from Yoruba in Ibadan, Nigeria (YRI), HapMap populations. The likelihood ratio test was used to test significance of the difference in risk allele frequency in Pimas (fR-Pima) with that in each HapMap population (fR-CEU, fR-CHB, and fR-YRI). For assessment of whether T2DM risk alleles were systematically higher in one population than in another, the mean of the GRS (μGRS) was compared between populations.

For a more general comparison of genetic distance between Pimas and other populations, the coancestry coefficient (FST) was calculated across all T2DM-susceptibility variants by the method of moments (43). Since interpretation of FST is most straightforward when sample sizes are equal, Pima allele frequencies used in these calculations were derived from a random sample of equal effective size to the corresponding HapMap population; effective sample size was estimated by the method of Yang et al. (44). For comparison of FST calculated across the T2DM markers with its genomic expectation, random markers were selected from a GWAS in Pimas (45). Since SNP characteristics may have influenced detection of the T2DM markers, each T2DM-associated SNP was matched to potential random SNPs by minor allele frequency in CEU, base pair type, and chromosome type (autosomal vs. X chromosome); to avoid selecting markers highly concordant with those for susceptibility to T2DM, we excluded a 2-Mb region on either side of the sentinel SNP from this selection. A total of 294,467 potentially matching random SNPs were thus identified. Significance of the difference between FST at T2DM variants and FST at random markers was calculated by a bootstrap procedure in which one random marker was selected for each T2DM variant in each iteration.

Excess Population Prevalence Attributable to Allele Frequency Differences

For quantification of the extent to which differences in T2DM risk allele frequencies can explain the difference in T2DM prevalence between Pimas and Europeans, standard multivariable epidemiologic methods for calculation of attributable fraction (46) were modified to calculate the genetic attributable fraction (GAF) for the population difference in prevalence. We define this as the proportion of the excess T2DM prevalence in a high-risk “target” population compared with a low-risk “reference” population attributable to differences in risk allele frequency. If P0 represents prevalence in the reference population (e.g., Europeans) and P1 is prevalence in the target population (e.g., Pimas), then

equation image

where P1adj is prevalence in the target population adjusted for the allele frequency differences (i.e., prevalence if the target population had the same risk allele frequencies as the reference population) (Eq. 1). Data from non-Hispanic white participants in the oral glucose tolerance subset of the U.S. National Health and Nutrition Examination Survey (NHANES) 2005–2010 were used for the reference population (http://www.cdc.gov/nchs/nhanes/nhanes_questionnaires.htm). These data were from 3,282 individuals age 12–84 years (1,585 women, 1,697 men; mean ± SD age 46.8 ± 20.1 years); 523 individuals (15.9%) had diabetes. These data were combined with Pimas of the same age range for calculation of GAF.

If genotypic data for all markers were available for all individuals, the quantities needed to calculate GAF could be derived from multivariable logistic regression. However, such data are not readily available for NHANES participants, so we developed an approximation that uses allele frequency and OR estimates from other sources. Adjusted prevalence in each population was obtained from the following logistic regression equation: logit (prevalence) = α0 + α1I + γ1(cov1) + . . . γm(covm), where I is an indicator variable that takes the value of 0 for the reference population and 1 for the target population and γ1–γm represent the coefficients corresponding to m covariates (centered at the mean values in the target population). Under an additive model with assumptions of Hardy-Weinberg equilibrium in both populations, independence among SNPs, and that the population OR changes as a function of the OR associated with each SNP and the difference in risk allele frequency, the expected value of α1, given that allele frequencies are the same as in the reference population, is as follows: α1adj = α1 − Σ[2βi(fR1i2 + fR1i[1 − fR1i] − fR0i2fRoi[1 − fRoi])], where βi is the logarithm of the OR for the ith SNP, fR1i is the risk allele frequency in the target population, and fR0i is the frequency in the reference population. For the present analyses, allele frequencies in the HapMap CEU population were taken as representative of the reference population. The values required for calculation of GAF (see Eq. 1) are as follows:

equation image

Simulation studies suggest that estimates of GAF derived by this method provide a good approximation of those derived from a multivariable regression in which all data are available for all individuals (Fig. 1). CIs and hypothesis tests for GAF were derived from a bootstrap procedure in which Pima, NHANES, and CEU populations were resampled and ORs were sampled from published values and standard errors.

Figure 1
Simulations comparing the approximate method for estimation of GAF with the “true” GAF, estimated from a multivariable logistic regression model containing all SNPs. Simulations were conducted generating 100,000 individuals from a reference ...

RESULTS

Association with T2DM

Eight SNPs (rs17106184 in FAF1, rs7578597 in THADA, rs3923113 in GRB14, rs831571 in PSMD6, rs6808574 in LPP, rs1531343 in HMGA2, rs7957197 in HNF1A, and rs17782313 in MC4R) were nearly monomorphic (minor allele frequency <0.01) in Pimas and were not analyzed for association. Table 1 shows the association for each of the remaining 55 SNPs with T2DM in Pima Indians, along with the test for heterogeneity in ORs between Pimas and Europeans. Nine SNPs, those in GCKR, ZBED3, CDKAL1, ZFAND3, KCNQ1, SPRY2, HMG20A, PRC1, and FTO, had nominally significant associations (P < 0.05) in Pimas in the same direction as the established association. The previously reported result with the KCNQ1 SNP rs2237892 was the strongest association. Ten SNPs, in IRS1, ADAMTS9, ARL15, ZFAND3, PTPRD, TCF7L2, MPHOSPH9, C2CD4A, SLC16A11, and DUSP9, showed nominally significant heterogeneity between Pimas and Europeans. Nonetheless, ORs were in the same direction as the established association for 39 of the 55 SNPs.

Table 1
Association of variants with T2DM in Pima Indians and comparison with association in Europeans

Associations with Metabolic Predictors of T2DM

Results for SNPs with nominally significant and directionally consistent associations with metabolic traits are shown in Table 2. Results for all SNPs are shown in Supplementary Table 1. The T2DM risk allele was associated with lower insulin secretion for SNPs in PROX1, IGF2BP2, ZBED3, DGKB-TMEM195, GLIS3, CDC123, HHEX, KCNQ1, and MNTR1B. The T2DM risk allele for SNPs in IRS1, PPARG, MNTR1B, PRC1, and SRR was associated with lower values of insulin sensitivity. Since the MNTR1B SNP was associated with both insulin secretion and total body insulin sensitivity, we further investigated its relationship with hepatic insulin sensitivity, measured in the clamp using radiolabeled glucose, and found that the risk allele was associated with lower sensitivity (r = −0.15, P = 0.002). The T2DM risk allele was significantly associated with higher percentage body fat for SNPs in PRC1 and ZFAND3. When BMI was analyzed in the larger population, the T2DM risk alleles for variants in GCK and FTO were associated with significantly higher BMI (Supplemental Table 2).

Table 2
Variants with significant (P < 0.05) and directionally consistent associations with metabolic traits

Multiallelic Association

Associations with the multiallelic GRS are shown in Fig. 2. The sum of the number of risk alleles over all 55 SNPs was significantly associated with T2DM (OR 1.05 per copy of a risk allele, P = 6.2 × 10−6). There was also a strong association between a greater number of T2DM risk alleles and lower values of insulin secretion such that each copy of a risk allele was associated with a 4% decrease in insulin secretion (P = 4.1 × 10−6). There was little association with insulin sensitivity or percentage body fat. When alleles were weighted by the logarithms of the published ORs in constructing the GRS, similar results were obtained (data not shown). When BMI was analyzed, results were similar to those seen with percentage body fat, but the inverse association was statistically significant (lower by 0.4% per risk allele, P = 1.2 × 10−5) (Supplementary Fig. 1). When the GRS was constructed using only the nine insulin secretion–associated SNPs, each risk allele was associated with a 13% decrease in insulin secretion; similarly, in a score constructed from the five insulin sensitivity SNPs, each risk allele was associated with a 7% decrease in insulin sensitivity (Supplementary Fig. 2). The insulin secretion score was associated with T2DM (OR 1.09, P = 2.7 × 10−4); when these nine SNPs were excluded from the global GRS, the T2DM association was modestly attenuated (OR 1.04, P = 1.1 × 10−3).

Figure 2
Relationship between the GRS across 55 loci and prevalence of T2DM, percentage body fat, insulin sensitivity, and insulin secretion in Pima Indians. The GRS is calculated as the number of risk alleles and shown in categories, plotted at the midpoint. ...

Heterogeneity

The test for heterogeneity in the effect on T2DM across all SNPs was statistically significant (P = 3.9 × 10−5) and negative in sign (Z* = −4.12); this indicates that the effects of these SNPs are on average weaker in Pimas than in Europeans. When the 18 SNPs with nominally significant association with T2DM or significant heterogeneity were excluded, the effects of the GRS on T2DM (OR 1.04, P = 4.9 × 10−4) and insulin secretion (effect −4%, P = 8.0 × 10−6) remained significant, as did evidence for heterogeneity (Z* = −2.82, P = 4.8 × 10−3).

Allele Frequency Differences between Pimas and Major Continental Populations

The difference in frequency of the T2DM risk allele between Pimas and HapMap populations is shown for each locus in Supplementary Fig. 3. The distribution of the GRS in each population is shown in Fig. 3. Mean GRS in Pimas (68.4) was slightly but significantly lower than in CEU (69.2, P = 0.049); mean GRS in Pimas was also significantly lower than in YRI (73.7, P = 4.4 × 10−38) but higher than in CHB (66.6, P = 9.6 × 10−5). When loci were weighted by the logarithms of the ORs in constructing the GRS, results were similar, except that the contrast in mean GRS between Pimas and CEU was more pronounced (P = 1.2 × 10−10).

Figure 3
Cumulative distribution of the GRS for T2DM in Pimas and in each of the HapMap populations (CHB, CEU, and YRI). In the left panel, the GRS was calculated as the sum of the number of risk alleles across all 63 loci, while in the right panel it is the sum ...

Genetic distances among populations across all 63 T2DM markers and across random markers are summarized in Fig. 4. FST across these T2DM loci was 0.163 (95% CI 0.154, 0.173) between Pimas and CEU, 0.138 (0.125, 0.152) between Pimas and CHB, and 0.232 (0.221, 0.244) between Pimas and YRI. These values were not significantly different from those derived from matched sets of SNPs randomly selected across the genome: FST 0.158 (0.106, 0.209) between Pimas and CEU (P = 0.83 for difference in FST), 0.129 (0.078, 0.180) between Pimas and CHB (P = 0.74), and 0.241 (0.173, 0.309) between Pimas and YRI (P = 0.80). Thus, differences in allele frequency are generally similar to those expected given genetic distances between populations.

Figure 4
Dendrogram summarizing genetic distances between Pima Indians and HapMap populations at 63 T2DM susceptibility loci and at randomly selected genomic markers. Genetic distance was calculated as FST, and the dendrogram was generated using PHYLIP. FST is ...

Population Risk Attributable to T2DM Loci

Results for the calculation of GAF for Pimas compared with Europeans are shown in Fig. 5A. The age-sex adjusted prevalence of T2DM was 48.2% in Pima Indians and 8.2% in non-Hispanic whites from NHANES (OR 10.5). The prevalence in Pimas adjusted to the frequency of risk alleles in CEU was slightly higher at 55.9%, resulting in a GAF of −0.19 (95% CI −0.34, −0.03); the low value of GAF reflects the lower value of the GRS in Pimas. When the 10 SNPs with statistically significant heterogeneity in the ORs were excluded from the calculation, GAF was −0.03 (−0.19, 0.08). Calculations were also conducted comparing non-Hispanic blacks in NHANES (n = 1,610, mean ± SD age 39.5 ± 20.1 years, 812 with diabetes) with non-Hispanic whites using allele frequencies derived from the African ancestry in the southwest U.S. (ASW) HapMap population. These analyses suggest that 66% of the excess prevalence in the black population is potentially attributable to allele frequency differences at these loci (GAF 0.66 [95% CI 0.32, 1.07]) (Figure 5B).

Figure 5
Calculation of the GAF in Pima Indians compared with non-Hispanic whites (NHW) from NHANES (left panel) and NHANES non-Hispanic blacks (NHB) compared with non-Hispanic whites (right panel). The age- and sex-adjusted prevalence of T2DM is shown on the ...

DISCUSSION

In recent years, many genetic variants reproducibly associated with T2DM have been identified. These have mostly been identified by GWAS in European populations. Many of these variants are also associated with T2DM in non-European populations, but there are instances of heterogeneity (8,10,20). The extent of association in high-risk populations, such as American Indians, is not well characterized. Our previous analyses in Pima Indians, with a much smaller number of SNPs, identified associations with SNPs in FTO and KCNQ1 (27,28); the KCNQ1 associations are subject to parent-of-origin effects and are particularly strong in Pimas (28). KLF14 variants also show parent-of-origin effects (28). Statistically significant heterogeneity between Pimas and Europeans at TCF7L2 was also observed, and a multiallelic score from eight SNPs was modestly associated with T2DM in Pimas and with diminished insulin secretion (20,27). In the present study, we have conducted a more complete survey of T2DM susceptibility variants in Pimas, including a total of 63 SNPs reproducibly associated with T2DM at genome-wide significance. These analyses identify additional SNPs that are nominally significantly associated with T2DM in Pimas in the same direction as in Europeans, including those in GCKR, ZBED3, CDKAL1, ZFAND3, SPRY2, HMG20A, and PRC1. Many of the T2DM susceptibility SNPs have effects in Pimas that are directionally consistent with those in Europeans, even if they were not individually statistically significant. Indeed, a multiallelic GRS that assesses effects of these variants in aggregate was statistically significant, even when SNPs with nominally significant effects or heterogeneity were excluded. The GRS was also strongly associated with diminished insulin secretion. Thus, the present findings suggest that the majority of T2DM-susceptibility variants do have modest effects on T2DM in this high-risk population but that some do not achieve statistical significance in the current sample size. Analyses in European populations suggest that the majority of T2DM-susceptibility variants influence T2DM risk through an effect on insulin secretion (2,47), and the current analyses suggest that this is also the case in Pimas.

Despite general consistency for most SNPs between the direction of association with T2DM in Pimas and that observed in the original GWAS, there were several SNPs that showed evidence for heterogeneity in effect between Pimas and Europeans. In addition to TCF7L2, nominally significant heterogeneity was observed at IRS1, ADAMTS9, ARL15, ZFAND3, PTPRD, C2CD4A, MPHOSPH9, SLC16A11, and DUSP9. With the exception of ZFAND3, which has previously been described as associated in East Asians but not Europeans (8), the effect in Pimas was weaker than that in Europeans. Furthermore, the combined test of heterogeneity across all loci indicated that effects were generally weaker in Pimas than in Europeans (even when SNPs with nominally significant association or heterogeneity were excluded). Thus, while most T2DM-susceptibility variants do have an effect on T2DM risk in Pimas, this effect is generally not as strong as it is in Europeans. It is possible that, despite the large sample sizes, this heterogeneity reflects overestimation of effects in Europeans. Given that functional variants at most of these loci have not been identified, however, some heterogeneity between Europeans and other populations might be expected on account of differing linkage disequilibrium patterns. Indeed, fine-mapping studies have suggested that population heterogeneity at GWAS signals derived from Europeans is at least partly due to differences in linkage disequilibrium patterns (18).

Recent studies have described divergence in allele frequency at T2DM-susceptibility variants between major continental populations that is greater than expected given genetic distances between these populations and a gradient in genetic risk for T2DM with risk alleles at highest frequency in Africans and at lowest frequency in East Asians (48,49). Such divergence in allele frequencies may reflect effects of natural selection in the different evolutionary histories of these populations. Prevalence of T2DM among Pima Indians is among the highest reported in the world, and if such evolutionary factors are responsible for this high prevalence, one might expect to see established T2DM risk alleles at high frequency in Pimas. The present analyses are consistent with previous studies, conducted with fewer SNPs (48), in that we observed the highest genetic risk scores in Africans (YRI) and the lowest in East Asians (CHB). However, genetic risk scores in Pimas were not particularly high and were comparable with, or lower than, those from low-risk populations, such as Europeans. The population differences in GRS observed here could reflect effects of genetic drift or natural selection. One study found that the Africa–East Asia gradient was greater than expected with random markers (49), which suggests natural selection, but a recent study that analyzed several global populations at 65 established T2DM-susceptibility loci suggested that T2DM-susceptibility alleles are generally evolutionarily neutral (50). Further work is needed to determine whether the high genetic risk scores for T2DM in African versus Asian populations is reflective of genetic drift or natural selection. However, in the present general analysis of genetic distances, we did not observe significant excess in the divergence between Pimas and other continental populations across established T2DM-susceptibility variants. This suggests that any overall effects of natural selection at these variants do not appear to have contributed to the high risk of T2DM in Pimas.

Regardless of the mechanisms by which population differences in risk allele frequency have arisen, such differences could explain population differences in prevalence of T2DM. The present analyses of GAF, however, suggest that differences in allele frequency at these established T2DM variants account for little of the increased population risk for T2DM in Pimas compared with European Americans. GWAS within Amerindian-derived populations may identify variants that are more likely to explain these population differences. Our recent GWAS comparing Pimas with young-onset T2DM to older nondiabetic individuals found association with a variant in DNER in Pimas but not in Europeans (45); this variant (rs1861612) shows little difference in allele frequency, however. A recent study suggested that the risk allele of rs75493593 in SLC16A11, which is more common in American Indians than Europeans, could explain ~20% of the excess risk in Mexican Americans compared with European Americans, ignoring the effects of all other loci (16). In the present study, we found that the risk allele at SLC16A11 is much more common in Pima Indians than in Europeans; however, its effect is outweighed by other variants at which the risk allele is less common Pimas, such that the overall extent to which established T2DM risk alleles can account for the excess prevalence in Pimas is negligible. In contrast, the present analyses suggest that 66% of the difference in T2DM prevalence between African Americans and European Americans is potentially attributable to allele frequency differences at these loci. Since transferability of European-derived T2DM variants to African Americans is somewhat uncertain given the highly divergent linkage disequilibrium patterns, the validity of the assumption that European-derived ORs represent causal effects may be questionable. Nonetheless, in light of the high proportion of excess prevalence between African Americans and Europeans that is attributable to differences in allele frequency at established T2DM variants, the fact that they account for none of the excess T2DM prevalence in Pimas seems remarkable.

In summary, the present analyses suggest that established T2DM variants are largely transferrable to high-risk populations, such as Pima Indians, albeit with weaker effects than in Europeans. However, differences in allele frequency across these established T2DM alleles account for little, if any, of the high T2DM prevalence in Pimas compared with populations of European ancestry. Thus, the high prevalence of T2DM in Pimas is likely the result of environmental factors or of genetic factors that remain largely unidentified.

Article Information

Acknowledgments. The authors thank the participants who volunteered for the study and the staff of the Phoenix Epidemiology and Clinical Research Branch who provided assistance.

Funding. This work was supported by the intramural research program of the National Institute of Diabetes and Digestive and Kidney Diseases.

Duality of Interest. No potential conflicts of interest relevant to this article were reported.

Author Contributions. R.L.H. wrote the manuscript, researched data, and contributed to discussion. R.R., S.K., Y.L.M., E.J.W., J.M.C., R.G.N., and L.J.B. researched data, contributed to discussion, and reviewed and edited the manuscript. R.L.H. is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Prior Presentation. Parts of this study were presented in abstract form at the Annual Meeting of the American Society of Human Genetics, San Francisco, CA, 6–10 November 2012.

Footnotes

This article contains Supplementary Data online at http://diabetes.diabetesjournals.org/lookup/suppl/doi:10.2337/db14-1715/-/DC1.

References

1. Zeggini E, Scott LJ, Saxena R, et al. .; Wellcome Trust Case Control Consortium . Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes. Nat Genet 2008;40:638–645 [PMC free article] [PubMed]
2. Voight BF, Scott LJ, Steinthorsdottir V, et al. .; MAGIC investigators; GIANT Consortium . Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet 2010;42:579–589 [PMC free article] [PubMed]
3. Morris AP, Voight BF, Teslovich TM, et al. .; Wellcome Trust Case Control Consortium; Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) Investigators; Genetic Investigation of ANthropometric Traits (GIANT) Consortium; Asian Genetic Epidemiology Network–Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium . Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet 2012;44:981–990 [PMC free article] [PubMed]
4. Unoki H, Takahashi A, Kawaguchi T, et al. . SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations. Nat Genet 2008;40:1098–1102 [PubMed]
5. Yasuda K, Miyake K, Horikawa Y, et al. . Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet 2008;40:1092–1097 [PubMed]
6. Tsai FJ, Yang CF, Chen CC, et al. . A genome-wide association study identifies susceptibility variants for type 2 diabetes in Han Chinese. PLoS Genet 2010;6:e1000847. [PMC free article] [PubMed]
7. Shu XO, Long J, Cai Q, et al. . Identification of new genetic risk variants for type 2 diabetes. PLoS Genet 2010;6:e1001127. [PMC free article] [PubMed]
8. Cho YS, Chen CH, Hu C, et al. .; DIAGRAM Consortium; MuTHER Consortium . Meta-analysis of genome-wide association studies identifies eight new loci for type 2 diabetes in east Asians. Nat Genet 2012;44:67–72 [PMC free article] [PubMed]
9. Kooner JS, Saleheen D, Sim X, et al. .; DIAGRAM; MuTHER . Genome-wide association study in individuals of South Asian ancestry identifies six new type 2 diabetes susceptibility loci. Nat Genet 2011;43:984–989 [PMC free article] [PubMed]
10. Saxena R, Saleheen D, Been LF, et al. .; DIAGRAM; MuTHER; AGEN . Genome-wide association study identifies a novel locus contributing to type 2 diabetes susceptibility in Sikhs of Punjabi origin from India. Diabetes 2013;62:1746–1755 [PMC free article] [PubMed]
11. Saxena R, Elbers CC, Guo Y, et al. .; Look AHEAD Research Group; DIAGRAM consortium . Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci. Am J Hum Genet 2012;90:410–425 [PubMed]
12. Tabassum R, Chauhan G, Dwivedi OP, et al. .; DIAGRAM; INDICO . Genome-wide association study for type 2 diabetes in Indians identifies a new susceptibility locus at 2q21. Diabetes 2013;62:977–986 [PMC free article] [PubMed]
13. Palmer ND, McDonough CW, Hicks PJ, et al. .; DIAGRAM Consortium; MAGIC Investigators . A genome-wide association search for type 2 diabetes genes in African Americans. PLoS ONE 2012;7:e29202. [PMC free article] [PubMed]
14. Dupuis J, Langenberg C, Prokopenko I, et al. .; DIAGRAM Consortium; GIANT Consortium; Global BPgen Consortium; Anders Hamsten on behalf of Procardis Consortium; MAGIC investigators . New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 2010;42:105–116 [PMC free article] [PubMed]
15. Gudmundsson J, Sulem P, Steinthorsdottir V, et al. . Two variants on chromosome 17 confer prostate cancer risk, and the one in TCF2 protects against type 2 diabetes. Nat Genet 2007;39:977–983 [PubMed]
16. Williams AL, Jacobs SB, Moreno-Macías H, et al. .; SIGMA Type 2 Diabetes Consortium . Sequence variants in SLC16A11 are a common risk factor for type 2 diabetes in Mexico. Nature 2014;506:97–101 [PMC free article] [PubMed]
17. Mahajan A, Go MJ, Zhang W, et al. .; DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium; Asian Genetic Epidemiology Network Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; Mexican American Type 2 Diabetes (MAT2D) Consortium; Type 2 Diabetes Genetic Exploration by Nex-generation sequencing in muylti-Ethnic Samples (T2D-GENES) Consortium . Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility. Nat Genet 2014;46:234–244 [PMC free article] [PubMed]
18. Carlson CS, Matise TC, North KE, et al. .; PAGE Consortium . Generalization and dilution of association results from European GWAS in populations of non-European ancestry: the PAGE study. PLoS Biol 2013;11:e1001661. [PMC free article] [PubMed]
19. Ng MC, Saxena R, Li J, et al. . Transferability and fine mapping of type 2 diabetes loci in African Americans: the Candidate Gene Association Resource Plus Study. Diabetes 2013;62:965–976 [PMC free article] [PubMed]
20. Guo T, Hanson RL, Traurig M, et al. . TCF7L2 is not a major susceptibility gene for type 2 diabetes in Pima Indians: analysis of 3,501 individuals. Diabetes 2007;56:3082–3088 [PubMed]
21. Knowler WC, Pettitt DJ, Saad MF, Bennett PH. Diabetes mellitus in the Pima Indians: incidence, risk factors and pathogenesis. Diabetes Metab Rev 1990;6:1–27 [PubMed]
22. American Diabetes Association Expert Committee on the Diagnosis and Classification of Diabetes Mellitus . Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Diabetes Care 1997;20:1183–1197 [PubMed]
23. Bunt JC, Krakoff J, Ortega E, Knowler WC, Bogardus C. Acute insulin response is an independent predictor of type 2 diabetes mellitus in individuals with both normal fasting and 2-h plasma glucose concentrations. Diabetes Metab Res Rev 2007;23:304–310 [PMC free article] [PubMed]
24. Tian C, Hinds DA, Shigeta R, et al. . A genomewide single-nucleotide-polymorphism panel for Mexican American admixture mapping. Am J Hum Genet 2007;80:1014–1023 [PubMed]
25. Hanis CL, Chakraborty R, Ferrell RE, Schull WJ. Individual admixture estimates: disease associations and individual risk of diabetes and gallbladder disease among Mexican-Americans in Starr County, Texas. Am J Phys Anthropol 1986;70:433–441 [PubMed]
26. Muller YL, Bogardus C, Beamer BA, Shuldiner AR, Baier LJ. A functional variant in the peroxisome proliferator-activated receptor gamma2 promoter is associated with predictors of obesity and type 2 diabetes in Pima Indians. Diabetes 2003;52:1864–1871 [PubMed]
27. Rong R, Hanson RL, Ortiz D, et al. . Association analysis of variation in/near FTO, CDKAL1, SLC30A8, HHEX, EXT2, IGF2BP2, LOC387761, and CDKN2B with type 2 diabetes and related quantitative traits in Pima Indians. Diabetes 2009;58:478–488 [PMC free article] [PubMed]
28. Hanson RL, Guo T, Muller YL, et al. . Strong parent-of-origin effects in the association of KCNQ1 variants with type 2 diabetes in American Indians. Diabetes 2013;62:2984–2991 [PMC free article] [PubMed]
29. Nair AK, Muller YL, McLean NA, et al. . Variants associated with type 2 diabetes identified by the transethnic meta-analysis study: assessment in American Indians and evidence for a new signal in LPP. Diabetologia 2014;57:2334–2338 [PMC free article] [PubMed]
30. Higgins JPT, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003;327:557–560 [PMC free article] [PubMed]
31. Gouda HN, Sagoo GS, Harding AH, Yates J, Sandhu MS, Higgins JP. The association between the peroxisome proliferator-activated receptor-gamma2 (PPARG2) Pro12Ala gene variant and type 2 diabetes mellitus: a HuGE review and meta-analysis. Am J Epidemiol 2010;171:645–655 [PMC free article] [PubMed]
32. Zhao Y, Ma YS, Fang Y, et al. . IGF2BP2 genetic variation and type 2 diabetes: a global meta-analysis. DNA Cell Biol 2012;31:713–720 [PubMed]
33. Franks PW, Rolandsson O, Debenham SL, et al. . Replication of the association between variants in WFS1 and risk of type 2 diabetes in European populations. Diabetologia 2008;51:458–463 [PMC free article] [PubMed]
34. Dehwah MA, Wang M, Huang QY. CDKAL1 and type 2 diabetes: a global meta-analysis. Genet Mol Res 2010;9:1109–1120 [PubMed]
35. Cauchi S, Del Guerra S, Choquet H, et al. . Meta-analysis and functional effects of the SLC30A8 rs13266634 polymorphism on isolated human pancreatic islets. Mol Genet Metab 2010;100:77–82 [PubMed]
36. Cugino D, Gianfagna F, Santimone I, et al. . Type 2 diabetes and polymorphisms on chromosome 9p21: a meta-analysis. Nutr Metab Cardiovasc Dis 2012;22:619–625 [PubMed]
37. Cai Y, Yi J, Ma Y, Fu D. Meta-analysis of the effect of HHEX gene polymorphism on the risk of type 2 diabetes. Mutagenesis 2011;26:309–314 [PubMed]
38. Tong Y, Lin Y, Zhang Y, et al. . Association between TCF7L2 gene polymorphisms and susceptibility to type 2 diabetes mellitus: a large Human Genome Epidemiology (HuGE) review and meta-analysis. BMC Med Genet 2009;10:15. [PMC free article] [PubMed]
39. Liu J, Wang F, Wu Y, et al. . Meta-analysis of the effect of KCNQ1 gene polymorphism on the risk of type 2 diabetes. Mol Biol Rep 2013;40:3557–3567 [PubMed]
40. Gong B, Yu J, Li H, Li W, Tong X. The effect of KCNJ11 polymorphism on the risk of type 2 diabetes: a global meta-analysis based on 49 case-control studies. DNA Cell Biol 2012;31:801–810 [PubMed]
41. Lathrop GM, Lalouel JM. Easy calculations of lod scores and genetic risks on small computers. Am J Hum Genet 1984;36:460–465 [PubMed]
42. Stouffer SA, Suchman EA, DeVinney LC, Star SA, Williams RM. How the volumes were produced. In The American Soldier, Volume 1: Adjustment During Army Life. Princeton, NJ, Princeton University Press, 1949, p 45
43. Weir BS, Hill WG. Estimating F-statistics. Annu Rev Genet 2002;36:721–750 [PubMed]
44. Yang Y, Remmers EF, Ogunwole CB, Kastner DL, Gregersen PK, Li W. Effective sample size: Quick estimation of the effect of related samples in genetic case-control association analyses. Comput Biol Chem 2011;35:40–49 [PMC free article] [PubMed]
45. Hanson RL, Muller YL, Kobes S, et al. . A genome-wide association study in American Indians implicates DNER as a susceptibility locus for type 2 diabetes. Diabetes 2014;63:369–376 [PMC free article] [PubMed]
46. Kooperberg C, Petitti DB. Using logistic regression to estimate the adjusted attributable risk of low birthweight in an unmatched case-control study. Epidemiology 1991;2:363–366 [PubMed]
47. Jonsson A, Ladenvall C, Ahluwalia TS, et al. . Effects of common genetic variants associated with type 2 diabetes and glycemic traits on α- and β-cell function and insulin action in humans. Diabetes 2013;62:2978–2983 [PMC free article] [PubMed]
48. Chen R, Corona E, Sikora M, et al. . Type 2 diabetes risk alleles demonstrate extreme directional differentiation among human populations, compared to other diseases. PLoS Genet 2012;8:e1002621. [PMC free article] [PubMed]
49. Corona E, Chen R, Sikora M, et al. . Analysis of the genetic basis of disease in the context of worldwide human relationships and migration. PLoS Genet 2013;9:e1003447. [PMC free article] [PubMed]
50. Ayub Q, Moutsianas L, Chen Y, et al. . Revisiting the thrifty gene hypothesis via 65 loci associated with susceptibility to type 2 diabetes. Am J Hum Genet 2014;94:176–185 [PubMed]

Articles from Diabetes are provided here courtesy of American Diabetes Association