Genetic variation near the hepatocyte nuclear factor-4 alpha gene predicts susceptibility to type 2 diabetes

The Finland-United States Investigation Of NIDDM Genetics (FUSION) study aims to identify genetic variants that predispose to type 2 diabetes by studying affected sibling pair families from Finland. Chromosome 20 showed our strongest initial evidence for linkage. It currently has a maximum logarithm of odds (LOD) score of 2.48 at 70 cM in a set of 495 families. In this study, we searched for diabetes susceptibility variant(s) at 20q13 by genotyping single nucleotide polymorphism (SNP) markers in case and control DNA pools. Of 291 SNPs successfully typed in a 7.5-Mb interval, the strongest association confirmed by individual genotyping was with SNP rs2144908, located 1.3 kb downstream of the primary beta-cell promoter P2 of hepatocyte nuclear factor-4 alpha (HNF4A). This SNP showed association with diabetes disease status (odds ratio [OR] 1.33, 95% CI 1.06-1.65, P = 0.011) and with several diabetes-related traits. Most of the evidence for linkage at 20q13 could be attributed to the families carrying the risk allele. We subsequently found nine additional associated SNPs spanning a 64-kb region, including the P2 and P1 promoters and exons 1-3. Our results and the independent observation of association of SNPs near the P2 promoter with diabetes in a separate study population of Ashkenazi Jewish origin suggests that variant(s) located near or within HNF4A increases susceptibility to type 2 diabetes.     

Variants in hepatocyte nuclear factor 4alpha are modestly associated with type 2 diabetes in Pima Indians

Single nucleotide polymorphisms (SNPs) within the hepatocyte nuclear factor 4alpha (HNF4alpha) gene are associated with type 2 diabetes in Finns and Ashkenazi Jews. Previous studies in both populations have reported linkage to type 2 diabetes near the HNF4alpha locus on chromosome 20q12-13. To investigate whether HNF4alpha is a diabetes susceptibility gene in Pima Indians, a population with the highest reported prevalence of type 2 diabetes but with no evidence for linkage of the disease on chromosome 20q, 19 SNPs across the promoter and coding region of HNF4alpha were genotyped for association analysis. In a group of 1,037 Pima Indians (573 diabetic and 464 nondiabetic subjects), three SNPs in HNF4alpha (rs3212183 and rs2071197 located in introns 3 and 1, respectively, and rs6031558, an extremely rare SNP located in the P2 promoter region) were modestly associated with type 2 diabetes (rs3212183 odds ratio [OR] 1.34 [95% CI 1.07-1.67], P = 0.009; rs2071197 1.34 [1.07-1.66], P = 0.008; and rs6031558 3.18 [1.03-9.84], P = 0.04, adjusted for age, sex, birth year, heritage, and family membership). We conclude that variants in HNF4alpha do not appear to be major determinants for type 2 diabetes in Pima Indians; however, HNF4alpha may have a minor role in type 2 diabetes susceptibility within this Native American population.     

Identification of PVT1 as a candidate gene for end-stage renal disease in type 2 diabetes using a pooling-based genome-wide single nucleotide polymorphism association study

To identify genetic variants contributing to end-stage renal disease (ESRD) in type 2 diabetes, we performed a genome-wide analysis of 115,352 single nucleotide polymorphisms (SNPs) in pools of 105 unrelated case subjects with ESRD and 102 unrelated control subjects who have had type 2 diabetes for > or =10 years without macroalbuminuria. Using a sliding window statistic of ranked SNPs, we identified a 200-kb region on 8q24 harboring three SNPs showing substantial differences in allelic frequency between case and control pools. These SNPs were genotyped in individuals comprising each pool, and strong evidence for association was found with rs2720709 (P = 0.000021; odds ratio 2.57 [95% CI 1.66-3.96]), which is located in the plasmacytoma variant translocation gene PVT1. We sequenced all exons, exon-intron boundaries, and the promoter of PVT1 and identified 47 variants, 11 of which represented nonredundant markers with minor allele frequency > or =0.05. We subsequently genotyped these 11 variants and an additional 87 SNPs identified through public databases in 319-kb flanking rs2720709 ( approximately 1 SNP/3.5 kb); 23 markers were associated with ESRD at P < 0.01. The strongest evidence for association was found for rs2648875 (P = 0.0000018; 2.97 [1.90-4.65]), which maps to intron 8 of PVT1. Together, these results suggest that PVT1 may contribute to ESRD susceptibility in diabetes.     

Calsquestrin 1 (CASQ1) gene polymorphisms under chromosome 1q21 linkage peak are associated with type 2 diabetes in Northern European Caucasians

Genome-wide scans in multiple populations have identified chromosome 1q21-q24 as one susceptibility region for type 2 diabetes. To map the susceptibility genes, we first placed a dense single nucleotide polymorphism (SNP) map across the linked region. We identified two SNPs that showed strong associations, and both mapped to within intron 2 of the calsequestrin 1 (CASQ1) gene. We tested the hypothesis that sequence variation in or near CASQ1 contributed to type 2 diabetes susceptibility in Northern European Caucasians by identifying additional SNPs from the public database and by screening the CASQ1 gene for additional variation. In addition to 15 known SNPs in this region, we found 8 new SNPs, 3 of which were in exons. A single rare nonsynonymous SNP in exon 11 (A348V) was not associated with type 2 diabetes. The associated SNPs were localized to the region between -1,404 in the 5' flanking region and 2,949 in intron 2 (P = 0.002 to P = 0.034). No SNP 3' to intron 2, including the adjacent gene PEA15, showed an association. The strongest associations were restricted to individuals of Northern European ancestry ascertained in Utah. A six-marker haplotype was also associated with type 2 diabetes (P = 0.008), but neither transmission disequilibrium test nor family-based association studies were significant for the most strongly associated SNP in intron 2 (SNP CASQ2312). An independent association of SNPs in introns 2 and 4 with type 2 diabetes is reported in Amish families with linkage to chromosome 1q21-q24. Our findings suggest that noncoding SNPs in CASQ1 alter diabetes susceptibility, either by a direct effect on CASQ1 gene expression or perhaps by regulating a nearby gene such as PEA15.     

Evidence that Rho guanine nucleotide exchange factor 11 (ARHGEF11) on 1q21 is a type 2 diabetes susceptibility gene in the Old Order Amish

Rho guanine nucleotide exchange factor 11 (ARHGEF11), located on chromosome 1q21, is involved in G protein signaling and is a pathway known to play a role in both insulin secretion and action. We genotyped 52 single nucleotide polymorphims (SNPs) in ARHGEF11 and compared the genotype frequencies of subjects with type 2 diabetes (n = 145) or type 2 diabetes/impaired glucose tolerance (IGT) (n = 293) with those of control subjects with normal glucose tolerance (NGT) (n = 358). Thirty SNPs, spanning the entire gene, were significantly associated with type 2 diabetes or type 2 diabetes/IGT. The most significantly associated SNP was rs6427340 (intron 2), in which the less common allele was the risk allele (odds ratio [OR] 1.82 [95% CI 1.20-2.70], P = 0.005 for type 2 diabetes vs. NGT and 1.79 [1.27-2.50], P = 0.0008 for type 2 diabetes/IGT vs. NGT). In an expanded set of nondiabetic subjects (n = 754), most of the type 2 diabetes-and IGT-associated SNPs were significantly associated with glucose levels during an oral glucose tolerance test, with the same SNP (rs6427340) showing the most significant associations (P = 0.007). All type 2 diabetes-and IGT-associated SNPs were in high linkage disequilibrium and constitute a single 133-kb haplotype block. These results, coupled with similar findings in Pima Indians, suggest that sequence variation in ARHGEF11 may influence risk of type 2 diabetes.     

A polymorphism in the AMPKalpha2 subunit gene is associated with insulin resistance and type 2 diabetes in the Japanese population

AMP-activated protein kinase (AMPK) acts as a fuel gauge for glucose and lipid metabolism. The gene encoding the alpha2 isoform of the catalytic subunit of AMPK (PRKAA2) is located at one of the Japanese type 2 diabetes loci mapped by our previous genome scan (1p36-32). PRKAA2 is, therefore, a good candidate gene for insulin resistance and type 2 diabetes. We screened all nine exons, their exon-intron boundaries, and the 5' and 3' flanking regions of PRKAA2 to identify single nucleotide polymorphisms (SNPs), and we genotyped 192 type 2 diabetic patients and 272 nondiabetic subjects to assess possible associations between genotypes or haplotypes and type 2 diabetes. None of the 10 SNPs genotyped was associated with type 2 diabetes, but the haplotype analysis, consisting of six representative SNPs, revealed one haplotype, with the A (minor) allele for rs2051040 and a major allele for the other five SNPs, to be associated with type 2 diabetes (P = 0.009). This finding was confirmed in two larger replication samples (657 case and 360 control subjects, P = 0.021; and 356 case and 192 control subjects from the same area in Japan, P = 0.007) and a significant P value was obtained in the joint haplotype analysis of all samples (1,205 case and 824 control subjects, P = 0.0001). Furthermore, insulin resistance was associated with rs2051040 in nondiabetic subjects, and those with the A (minor) allele had a higher homeostasis model assessment of insulin resistance index than those who did not (initial control subjects [n = 272], P = 0.002; and joint replication control subjects [n = 552], P = 0.037). We speculate that the PRKAA2 gene influences insulin resistance and susceptibility to type 2 diabetes in the Japanese population.     

Searching for type 2 diabetes genes on chromosome 20

Genome scans in families with type 2 diabetes identified a putative locus on chromosome 20q. For this study, linkage disequilibrium mapping was used in an effort to narrow a 7.3-Mb region in an Ashkenazi type 2 diabetic population. The region encompassed a 1-logarithm of odds (LOD) interval around the microsatellite marker D20S107, which gave a LOD score of >3 in linkage analysis of a combined Caucasian population. This 7.3-Mb region contained 25 known and 99 predicted genes. Predicted single nucleotide polymorphisms (SNPs) were chosen from public databases and validated. Two SNPs were unique to the Ashkenazi. Here, 91 SNPs with a minor allele frequency of >or=10% were genotyped in pooled DNA from 150 case subjects and 150 control subjects of Ashkenazi Jewish descent. The SNP association study showed that SNP rs2664537 in the TIX1 gene had a significant P value of 0.035, but the finding did not replicate in an additional case pool. In addition, HNF4a and Mybl2 were screened for mutations and new polymorphisms. No mutations were identified, and a new nonsynonymous SNP (R687C in exon 14 of Mybl2) was found. The limits to this type of association study are discussed.     

EIF4A2 is a positional candidate gene at the 3q27 locus linked to type 2 diabetes in French families

One of the most replicated loci influencing type 2 diabetes-related quantitative traits (quantitative trait loci [QTL]) is on chromosome 3q27 and modulates both type 2 diabetes-and metabolic syndrome-associated phenotypes. A QTL for type 2 diabetes age of onset (logarithm of odds [LOD] score = 3.01 at D3S3686, P = 0.0001) was identified in a set of French families. To assess genetic variation underlying both age-of-onset QTL and our previous type 2 diabetes linkage in a 3.87-Mb interval, we explored 36 single nucleotide polymorphisms (SNPs) in two biologically relevant candidate genes for glucose homeostasis, kininogen (KNG1), and eukaryotic translation initiation factor 4alpha2 (EIF4A2). Analysis of 148 families showed significant association of a frequent SNP, rs266714, located 2.47 kb upstream of EIF4A2, with familial type 2 diabetes (family-based association test, P = 0.0008) and early age of onset (P = 0.0008). This SNP also contributes to both age-of-onset QTL (1.13 LOD score decrease P = 0.02) and type 2 diabetes linkage (genotype identical-by-descent sharing test, P = 0.02). However, no association was observed in three independent European diabetic cohorts. EIF4A2 controls specific mRNA translation and protein synthesis rate in pancreatic beta-cells, and our data indicates that EIF4A2 is downregulated by high glucose in rat beta-INS832/13 cells. The potential role of EIF4A2 in glucose homeostasis and its putative contribution to type 2 diabetes in the presence of metabolic stress will require further investigation.     

Liver pyruvate kinase polymorphisms are associated with type 2 diabetes in northern European Caucasians

Pyruvate kinase is a key glycolytic enzyme. Isoforms that are expressed in the red cell, liver, pancreatic beta-cells, small intestine, and proximal renal tubule are encoded by the 12 exons of the PKLR gene, which maps to chromosome 1q23. We hypothesized that common variants of the PKLR gene could account for the linkage of diabetes to this region. We screened the promoter regions, exons and surrounding introns, and the 3' untranslated region for mutations. We identified five single-nucleotide polymorphisms (SNPs), and only one (V506I, exon 11) altered the coding sequence. We tested the five SNPs, a poly-T insertion-deletion polymorphism, and an ATT triplet repeat in 131 unrelated diabetic patients and 118 nondiabetic control subjects. The V506I variant was rare and not associated with type 2 diabetes. The four SNPs and the insertion-deletion polymorphism were associated with diabetes, with a 10% difference between individuals with diabetes and nondiabetic individuals (P = 0.001-0.011, relative risk for minor allele 1.85). The same trend was found for the ATT repeat (P = 0.029). Common variants in the PKLR are associated with increased risk of type 2 diabetes, but because of strong linkage disequilibrium between variants, the actual susceptibility allele may be in a different gene.     

Molecular analysis of KCNJ10 on 1q as a candidate gene for Type 2 diabetes in Pima Indians

The KCNJ10 gene is located within a region on chromosome 1q linked to type 2 diabetes in the Pima Indians and six other populations. We therefore investigated this gene as a potential type 2 diabetes candidate gene in Pima Indians. KCNJ10 consists of two exons, spans approximately 33 kb, and we identified eight single-nucleotide polymorphisms (SNPs), including one (SNP2) in the coding region leading to a Glu359Lys substitution. Association studies were carried out in a case-control group composed of 149 affected and 150 unaffected Pimas, and the linkage analysis was performed in a linkage set of 1,338 Pimas. SNP1 in the promoter and SNP2 in the intron, which were in a complete linkage disequilibrium, and SNP5 in the 3' untranslated region showed association with diabetes in the case-control group (P = 0.02 and P = 0.01, respectively). When genotyped in the linkage set, only the KCNJ10-SNP1 variant showed a modest association with type 2 diabetes (P = 0.01). KCNJ10-SNP1 is in a strong linkage disquilibrium with SNP14 of the adjacent KCNJ9 locus, which we previously found to be associated with type 2 diabetes. After adjustment for KCNJ10-SNP1, the original linkage score at this locus was marginally reduced from 3.1 to 2.9. We conclude that these variants in KCNJ10 are unlikely to be the cause of linkage of type 2 diabetes with 1q in Pima Indians.     

Common variation in the LMNA gene (encoding lamin A/C) and type 2 diabetes: association analyses in 9,518 subjects

Mutations in the LMNA gene (encoding lamin A/C) underlie familial partial lipodystrophy, a syndrome of monogenic insulin resistance and diabetes. LMNA maps to the well-replicated diabetes-linkage region on chromosome 1q, and there are reported associations between LMNA single nucleotide polymorphisms (SNPs) (particularly rs4641; H566H) and metabolic syndrome components. We examined the relationship between LMNA variation and type 2 diabetes (using six tag SNPs capturing >90% of common variation) in several large datasets. Analysis of 2,490 U.K. diabetic case and 2,556 control subjects revealed no significant associations at either genotype or haplotype level: the minor allele at rs4641 was no more frequent in case subjects (allelic odds ratio [OR] 1.07 [95% CI 0.98-1.17], P = 0.15). In 390 U.K. trios, family-based association analyses revealed nominally significant overtransmission of the major allele at rs12063564 (P = 0.01), which was not corroborated in other samples. Finally, genotypes for 2,817 additional subjects from the International 1q Consortium revealed no consistent case-control or family-based associations with LMNA variants. Across all our data, the OR for the rs4641 minor allele approached but did not attain significance (1.07 [0.99-1.15], P = 0.08). Our data do not therefore support a major effect of LMNA variation on diabetes risk. However, in a meta-analysis including other available data, there is evidence that rs4641 has a modest effect on diabetes susceptibility (1.10 [1.04-1.16], P = 0.001).     

Genetic analysis of ADIPOR1 and ADIPOR2 candidate polymorphisms for type 2 diabetes in the Caucasian population

Adiponectin is a metabolic link between adipose tissue and insulin action, mediating part of obesity-associated insulin resistance and type 2 diabetes. Two adiponectin receptors have been identified, and we investigated whether sequence variations in adiponectin receptor 1 (ADIPOR1) and adiponectin receptor 2 (ADIPOR2) genes could contribute to the genetic risk for type 2 diabetes in a case-control study of 1,498 Caucasian subjects. We sequenced the putative functional regions of the two genes in 48 subjects and selected single nucleotide polymorphisms (SNPs) from the public database. Five SNPs in ADIPOR1 and 12 in ADIPOR2 were tested for association with type 2 diabetes. No SNP of ADIPOR1 showed association in any of the samples from the French population. In contrast, three SNPs of ADIPOR2 showed nominal evidence for association with type 2 diabetes before correction for multiple testing (odds ratio [OR] 1.29-1.37, P = 0.034-0.014); only rs767870, located in intron 6, was replicated in an additional diabetes dataset (n = 636, OR 1.29, P = 0.020) with significant allelic association from the overall meta-analysis of 2,876 subjects (adjusted OR 1.25 [95% CI 1.07-1.45], P = 0.0051). In conclusion, our data suggest a modest contribution of ADIPOR2 variants in diabetes risk in the French population.     

Linkage of plasma adiponectin levels to 3q27 explained by association with variation in the APM1 gene

We performed a genome-wide linkage scan of plasma adiponectin levels in 569 nondiabetic participants in the Amish Family Diabetes Study. The highest logarithm of odds (LOD) score (2.13; P = 0.0009) occurred on chromosome 3q27 between markers D3S1602 and D3S1580, which flank APM1/ACDC, the adiponectin gene. The APM1 +2019 A/- insertion/deletion polymorphism in the 3' untranslated region (single nucleotide polymorphism [SNP] +2019; deletion allele frequency 0.30 in Amish) showed strong association with adiponectin levels in a dosage-dependent manner in a direction consistent with that reported in previous studies, with deletion heterozygosity increasing adiponectin levels by 1.3 +/- 0.5 microg/ml and deletion homozygosity increasing levels by 3.0 +/- 0.8 microg/ml (P < 0.0001). Two other SNPs, rs2241766 and rs1501299, showed moderate association. In a subset of 523 subjects genotyped for both SNP +2019 and rs2241766, including the APM1 SNP +2019 genotype as a covariate reduced the linkage signal at 3q27 by 1.26 LOD units (from 2.22 to 0.96) and including both SNPs reduced the signal by 1.51 LOD units (to 0.71). These findings, combined with a two-point LOD score of 2.35 for SNP +2019, provide evidence that variation in APM1 is responsible for linkage of adiponectin levels to 3q27 in the Old Order Amish.