DOI: 10.2337/diabetes.55.04.06.db05-0727 © 2006 by the American Diabetes Association
A Polymorphism in the AMPK
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
2 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.
Insulin resistance has been well demonstrated to be a fundamental element in the etiology of type 2 diabetes. The molecular mechanisms and genetic factors involved in insulin resistance have been extensively investigated. AMP-activated protein kinase (AMPK) is an enzyme that is activated by physiological and pharmacological stimuli, including muscle contraction (1), the antidiabetic agent metformin (2), and adipokines, such as leptin (3) or adiponectin (4). It is a heterotrimeric protein composed of a catalytic subunit (
) and two regulatory subunits (ß and
). The catalytic subunit
has two isoforms (
1 and
2), and while the
1 subunit isoform is ubiquitously expressed, the
2 subunit isoform, PRKAA2, is predominantly found in skeletal muscle and liver. Activation of AMPK in skeletal muscle leads to increased glucose uptake and enhanced insulin sensitivity (5), whereas in the liver AMPK activation inhibits glucose production (6). These characteristics of AMPK make this enzyme one of the key regulators of insulin sensitivity and glucose homeostasis. While AMPK
1–/– mice have no apparent metabolic defects, AMPK
2–/– mice exhibit insulin resistance (7). Moreover, PRKAA2 is located on 1p36-32 (8), which is reportedly linked to type 2 diabetes in the Japanese population (9). PRKAA2 is therefore a good candidate for the susceptibility gene to insulin resistance and type 2 diabetes.
| RESEARCH DESIGN AND METHODS |
|---|
|
|
|---|
Screening and genotyping of PRKAA2.
We screened all nine exons, including their exon-intron boundaries, 5' untranslated region, 3' untranslated region, and 2,000-bp upstream region of PRKAA2 in 32 diabetic patients to detect new single nucleotide polymorphisms (SNPs). The primer sequences are available from the authors. SNPs were genotyped by direct sequencing. Sequencing reactions were performed with the BigDye terminator (Applied Biosystems, Foster City, CA) and resolved using an ABI 3700 automated DNA sequencer (Applied Biosystems). The results were integrated using a Sequencher (Gene Codes, Ann Arbor, MI), and individual SNPs were manually genotyped. Ambiguous base identifications were genotyped twice. SNPs were identified based on the sequence reported in GenBank (PRKAA2: NT_032977).
Statistical analysis.
Genotype or allele associations for each SNP with type 2 diabetes were tested by Fishers exact test. Homeostasis model assessment (HOMA) of insulin resistance was calculated as described previously for the nondiabetic subjects (11) in the initial samples because insulin therapy or treatment with oral hypoglycemic agents for type 2 diabetes would presumably alter their insulin levels. A multiple regression analysis was used to test for associations between SNPs and insulin resistance after adjustment for age, sex, and BMI. The statistical analyses were performed using JMP for Windows version 4.00 software (SAS Institute, Cary, NC). P values were corrected by Bonferroni adjustment, and a P value <0.005 (i.e., 0.05 divided by the total number of SNPs) was considered significant in the initial study. The statistical power was calculated based on a test for differences in proportions of alleles between case and control subjects (described in detail by Ohashi et al. [12]).
Haplotype analysis.
To examine the linkage disequilibrium (LD) structure, pairwise LD, D', and r2 between SNPs and haplotype frequencies were estimated via the method of maximum likelihood from two-locus genotype data using the E-M algorithm under the assumption of Hardy-Weinberg equilibrium (13). For the estimation of haplotype frequencies, we selected one of the SNPs as a tagging SNP from every set of SNPs with r2 > 0.80. All haplotypes were jointly tested for association with disease status by performing a 2 x n
2 test of independence in a permutation procedure, where n indicates the number of haplotypes with a frequency >0. Individual haplotypes were also tested for association with disease status with a 2 x 2
2 test of independence in a permutation procedure. In the permutation procedure, to account for the variability introduced by the haplotype frequency estimation, significance was assessed by permuting case and control status and recalculating the test statistic 1,000 times for each of the sample sets and 10,000 times for the combined sample set. The above calculations were performed with SNPAlyze V3.2 Pro software (Dynacom, Yokohama, Japan).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
|
We next investigated whether the polymorphisms in PRKAA2 were associated with insulin resistance, as assessed by HOMA of insulin resistance. The association was compared between subjects with and without the minor allele. Among the 10 SNPs investigated, only rs2051040 was associated with insulin resistance. Subjects with the A (minor) allele for rs2051040 had a higher HOMA of insulin resistance than those without it (AA/AG vs. GG 1.95 ± 0.08 vs. 1.49 ± 0.10, regression coefficient = 0.18, P = 0.002). This finding was also confirmed in the joint replication samples (AA/AG vs. GG 1.83 ± 0.07 vs. 1.59 ± 0.09, regression coefficient = 0.099, P = 0.037). No other SNPs were associated with insulin resistance or any other diabetes-related quantitative traits, such as age, sex, BMI, fasting glucose, fasting insulin, HbA1c, and HOMA-ß.
PRKAA2 is a good candidate for the susceptibility gene to insulin resistance and type 2 diabetes. We therefore focused on PRKAA2 and genotyped 10 SNPs spanning from the promoter region to the 3' untranslated region. Even though alleles with an odds ratio of 1.2–1.3 could be detected with 80% power for the total sample for allele frequencies of 0.2 and 0.4, which correspond to the commonly observed allele frequencies in this study, we were unable to find associations between any single polymorphisms and type 2 diabetes. However, one common haplotype with the A (minor) allele for rs2051040 and a major allele for all the other SNPs was associated with type 2 diabetes. To jointly test all common haplotypes for associations with disease status, we also used PHASE version 2.1, performing a case-control permutation test with the default setting. In agreement with the haplotype analysis using SNPAlyze, haplotype analysis using PHASE software revealed a significant difference (P = 0.01) in haplotype frequencies between case and control subjects (data not shown). Likewise, SNPAlyze also revealed a significant difference (global P = 0.004) in a permutation test for differences in haplotype frequencies between case and control subjects. We therefore conclude that there is a significant difference in haplotype frequencies between case and control subjects and that this difference may be attributable to the haplotype AGTAAT, i.e., the haplotype with a minor (A) allele for rs2051040 and a major allele for all the other SNPs.
This intronic SNP rs2051040 alone was not associated with type 2 diabetes, but it was significantly associated with insulin resistance (P = 0.002), consistent with the known function of AMPK. The association was still significant even when the conservative Bonferroni adjustment was taken into account. Subjects with the A allele for rs2051040 had more marked insulin resistance, and a haplotype containing the A allele for rs2051040 was observed more frequently in case compared with control subjects. The initial results were further confirmed in the joint replication samples.
The risk haplotype frequency is quite close to the frequency of rs2051040 in the case subjects, for example 0.38 and 0.41 in the initial sample, but quite different, 0.29 and 0.37, respectively, in the control subjects. The A allele of rs2051040 is present on a rare haplotype, AGCAAT, whose frequency is higher in the control than in the case subjects, and the presence of this haplotype appears to be responsible for the discrepancy between the A allele of rs2051040 and risk haplotype frequencies in the case and control subjects (online appendix Table 4).
We note two possible reasons why the association between rs2051040 and type 2 diabetes is evident in haplotypic, but not individual SNP, association analyses. One possibility is that an unidentified SNP, which is in LD with this risk haplotype but is in weaker LD with rs2051040, is the SNP actually causing type 2 diabetes. The other possibility is that rs2051040 and another SNP, contained in this risk haplotype, function in a coordinate manner to increase the risk of type 2 diabetes. The HapMap shows additional SNPs in the Chinese/Japanese sample between exon 1 and 2 of PRKAA2, where SNP finding was scarce in our study. An intronic SNP that is associated with type 2 diabetes in Japanese subjects may lie in this region. Further research to identify either the true causative SNP or as-yet-unidentified SNPs functioning in a coordinate manner with rs2051040 is needed.
| ACKNOWLEDGMENTS |
|---|
We thank Ms. Y. Okada for technical assistance.
| FOOTNOTES |
|---|
Additional information for this article can be found in an online appendix available at http://diabetes.diabetesjournals.org.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received for publication June 7, 2005 and accepted in revised form November 14, 2005
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Arad, C. E. Seidman, and J.G. Seidman AMP-Activated Protein Kinase in the Heart: Role During Health and Disease Circ. Res., March 2, 2007; 100(4): 474 - 488. [Abstract] [Full Text] [PDF] |
||||
![]() |
N J Spencer-Jones, D Ge, H Snieder, U Perks, R Swaminathan, T D Spector, N D Carter, and S D O'Dell AMP-kinase {alpha}2 subunit gene PRKAA2 variants are associated with total cholesterol, low-density lipoprotein-cholesterol and high-density lipoprotein-cholesterol in normal women J. Med. Genet., December 1, 2006; 43(12): 936 - 942. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum |