Abstract
Aims/hypothesis
Two recent genome-wide association studies have identified several novel type 2 diabetes susceptibility variants in intron 15 of the KCNQ1 gene. We aimed to evaluate the effects of the variants in KCNQ1 on type 2 diabetes and metabolic traits in the population of mainland China.
Methods
Three candidate single nucleotide polymorphisms were genotyped in 1,912 individuals with type 2 diabetes and 2,041 normal controls using the ligase detection reaction method.
Results
We confirmed the association of KCNQ1 with type 2 diabetes in the population of mainland China. Allele frequency ORs of the three single nucleotide polymorphisms (SNPs) were: rs2237892 (OR 1.19, 95% CI 1.08–1.31, p = 3.0 × 10−4); rs2237895 (OR 1.20, 95% CI 1.09–1.32, p = 1.9 × 10−4); and rs2237897 (OR 1.24, 95% CI 1.13–1.36, p = 3.9 × 10−5). We also found a significant difference in the distribution of the global haplotypes between the type 2 diabetes group and the normal control group (p = 2.6 × 10−5). In addition, in the control group SNP rs2237892 was marginally associated with increasing fasting plasma glucose and SNPs rs2237892 and rs2237897 were associated with HbA1c. Furthermore, for all three variants, homozygous carriers of the diabetes-associated allele had significantly decreased BMI and waist circumferences.
Conclusions/interpretation
Our investigation confirmed the effects of KCNQ1 variants on type 2 diabetes risk in the Chinese population.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
Type 2 diabetes is a major health problem affecting more than 200 million individuals worldwide. Its prevalence is increasing in many countries, especially in the developing countries [1]. For example, the total number of people with diabetes in China is estimated to increase from 20.8 million in 2000 to 42.3 million by 2030. Type 2 diabetes is characterised by the presence of insulin resistance and pancreatic beta cell dysfunction, resulting from the interaction of genetic and environmental factors.
So far a few genes, identified through linkage scans or the candidate gene approach, have been confirmed as being associated with type 2 diabetes e.g. PPARG, KCNJ11 and TCF7L2 [2–4]. Recently, the use of genome-wide association scans as an investigative tool has led to a qualitative leap in identifying diabetes-related genes. Genome-wide association studies conducted in large-scale case–control samples by several independent European and American research groups have identified several novel genes and loci with modest effects on the risk of type 2 diabetes (OR 1.14–1.20), such as SLC30A8, HHEX, CDKAL1, CDKN2A, CDKN2B, IGF2BP2 and FTO [5–10]. These studies have been replicated and results confirmed in several other ethnic groups [11–17].
The two most recent genome-wide association studies, in the Japanese population, have identified the association of a novel gene (KCNQ1) with the risk of type 2 diabetes; a similar association was also detected in some Asian and European populations in these studies [18, 19]. Recently, Tan et al. also found variants in KCNQ1 to be associated with type 2 diabetes, fasting glucose and beta cell function in 3,734 individuals belonging to three ethnic groups living in Singapore [20]. It is well known that there are significant differences in the frequencies of some genetic variations among different ethnic groups and geographic regions. Therefore, we decided to investigate further the contribution of KCNQ1 to the aetiology of type 2 diabetes and to determine whether variants of the KCNQ1 gene were associated with the susceptibility to type 2 diabetes and diabetes-related metabolic traits in the population of mainland China. We also used haplotype analysis and the best-fitting model test in our investigation.
Methods
Participants
In this study, 3,953 Chinese Han participants were recruited from Shanghai, China; these comprised 1,912 unrelated type 2 diabetic individuals (785 men, 1,127 women; age 63.9 ± 9.5 years) and 2,041 control individuals (635 men, 1,406 women; age 58.1 ± 9.4 years). The study population overlapped completely with that used in our previous study [16]. Diabetic participants were defined in accordance with WHO criteria. Controls with a fasting plasma glucose concentration <6.1 mmol/l were enrolled from the same geographical region. A standard informed consent procedure was included in the protocol, and was reviewed and approved by the Ethics Committee of the Shanghai Institute for Biological Sciences. Participants gave their consent after the nature of study had been fully explained.
Blood samples were drawn for biochemical measurements (fasting plasma glucose, HbA1c, total cholesterol, triacylglycerol, HDL-cholesterol and LDL-cholesterol). Height, weight, waist and hip circumferences, and blood pressure were measured in all individuals. Data are showed as medians (25–75% range) or means ± SD (Table 1).
Single nucleotide polymorphism selection
For the KCNQ1 gene, five single nucleotide polymorphisms (SNPs) (rs2237892, rs2237895, rs2237897, rs2074196 and rs2283228) were identified as having the most significant association with type 2 diabetes based on previous studies [18–20]. After considering the linkage disequilibrium (LD) structure based on HapMap Han Chinese and the pairwise linkage disequilibrium D′ and r 2 values, we used rs2237892 to represent rs2074196 and rs2283228 (r 2 = 0.73 and 0.95, respectively) and thus selected three representative variants (rs2237892, rs2237895 and rs2237897) to be genotyped in our study.
Genotyping
High molecular weight genomic DNA was prepared from venous blood using the QuickGene 610L Automatic DNA/RNA Extraction System (Fujifilm, Tokyo, Japan). All the genotyping experiments were done using ligase detection reactions [21]. All genotyping success rates were above 98% and all mismatch rates were below 1% in 1641 duplicate samples.
Statistical analysis
We used SHEsis [22] to perform the Hardy–Weinberg equilibrium test, to compare the differences of allele, genotype and haplotype frequencies between cases and controls and to calculate LD coefficients (D′ and r 2) in the control group.
To establish the closest best-fit model for each SNP, we carried out the logistic regression analysis by comparing additive, dominant and recessive models with age, sex and loge BMI as covariants, and the model that had the lowest Akaike’s information criterion value was considered to be the best-fitting model for the respective SNP. In the additive model, homozygotes for the risk allele (1/1), heterozygotes (1/0) and homozygotes for the non-risk allele (0/0) were coded to an ordered categorical variable for the genotype (2, 1 and 0, respectively). The dominant model was defined as 1/1 + 1/0 vs 0/0 and the recessive model as 1/1 vs 1/0 + 0/0. The association of the SNPs with type 2 diabetes was assessed by logistic regression after adjusting for sex, age and log e BMI. The population attributable risk (PAR) [23] was calculated as:
where n c is the number if individuals in each case group and OR i is equal to β × Y (the estimated OR of the ith individual in each case group [β is the regression coefficient vector in terms of Y, and Y is the term for the respective SNP in the logistic regression model]). In addition, to analyse the independent effect of each SNP in the KCNQ1 gene, we carried out a conditional analysis by including the most significant SNP in the region and the SNP being tested, coded according to the most appropriate genetic model, combined with age, sex and log e BMI as covariates, in a logistic regression analysis.
For quantitative traits in relation to genotypes, a general linear statistical method was used, applying additive, dominant and recessive models while adjusting for the effect of age and sex (BMI, waist and waist to hip ratio), or age, sex and log e BMI (all other traits). The statistical analyses were performed using the SPSS (SPSS, Chicago, IL, USA) program. A p value <0.05 was considered significant.
Results
The allele and genotype distribution are summarised in Table 2. All three variants genotyped were in Hardy–Weinberg equilibrium in both the case and control groups. Allele frequency ORs of the three SNPs, rs2237892 (OR 1.19, 95% CI 1.08–1.31, p = 3.0 × 10−4), rs2237895 (OR 1.20, 95% CI 1.09–1.32, p = 1.9 × 10−4) and rs2237897 (OR 1.24, 95% CI 1.13–1.36, p = 3.9 × 10−5) were significantly associated with type 2 diabetes in our sample. The risk alleles were in the same direction as in previous studies [18, 19]. The best-fitting genetic model for all three was the additive model (see Electronic supplementary material [ESM] Table 1). Under the additive model with adjustment for age, sex and log e BMI, significant association was detected in all three polymorphisms (rs2237892, OR 1.23, 95% CI 1.11–1.36, p = 1.1 × 10−4, PAR = 23.9%; rs2237895, OR 1.23, 95% CI 1.11–1.37, p = 7.8 × 10−5, PAR = 12.6%; rs2237897, OR 1.28, 95% CI 1.16–1.41, p = 2.0 × 10−6, PAR = 27.3%; Table 2).
The LD structure (rs2237892/rs2237895/rs2237897) is shown in Fig. 1. LD coefficients (D′ and r 2) among the three variants (see ESM Table 2) suggested that the degree of linkage disequilibrium was relatively high. So haplotypes were constructed with all three variants and the frequencies of haplotypes were then analysed by excluding the rare haplotypes (those below 2% frequency in cases or controls); this revealed a significant difference in the distribution of the global haplotypes between the type 2 diabetes group and the normal controls (p = 2.6 × 10−5). The C-C-C haplotype frequency was significantly higher (p = 2.1 × 10−3) and the T-A-T haplotype frequency was significantly lower (p = 4.0 × 10−5) in the group with type 2 diabetes compared with the control group (see ESM Table 3).
As treatment for diabetes may have affected metabolic relationships, only non-diabetic individuals were included in the metabolic traits analysis. We found all three variants were associated with BMI and waist measurement, with the homozygous carriers of diabetes-associated allele showing the lowest value. Under the recessive model, rs2237892 and rs2237897 were associated with HbA1c levels and rs2237892 also showed association with fasting plasma glucose levels under the additive model (Table 3). Apart from these associations, the variants were not associated with any other measures in our study.
Discussion
We found the variants located at intron 15 of KCNQ1, rs2237892, rs2237895 and rs2237897, to be associated with type 2 diabetes. These results are consistent with previous studies [18–20]. The variant rs2237897 showed the strongest association with type 2 diabetes after testing for independent effects, which was consistent with the study of Tan et al. [20]. In addition, the conditional independent effects test of rs2237892, rs2237895 and rs2237897 confirmed that they were in a linkage disequilibrium block, as the effects of rs2237892 and rs2237895 on type 2 diabetes could be equally attributable to rs2237897 (see ESM Table 4).
KCNQ1, which encompasses 404 kb, is located at chromosome 11p15.5, not far from a candidate region at 11p13–p12 that two independent studies have linked to type 2 diabetes [24, 25]. KCNQ1 encodes a pore-forming subunit of the voltage-dependent K+ ion channel, which is expressed mainly in the heart and, to a lesser extent, in the brain, adipose tissue and pancreas. The two most important roles of potassium voltage-gated channel, KQT-like (KCNQ1) channels are: (1) repolarisation of the cardiac tissue following an action potential; and (2) water and salt transport in epithelial tissues. Moreover, mutations in KCNQ1 can lead to dysfunction of the channel and cause the cardiac long QT syndrome LQTS1 that, in turn, may lead to serious arrhythmias, ventricular fibrillation and cardiac arrest. Kcnq1-knockout mice have been reported to show cardiac dysfunctions such as prolonged QT interval, hypochlorhydria, hypergastrinaemia, gastric hyperplasia and vacuolation of the parietal cells [26, 27]. For type 2 diabetes, homozygous Kcnq1-null mice have been reported not to show hyperglycaemia or glucose intolerance, and the contribution of the KCNQ1-encoded protein to the molecular pathogenesis of type 2 diabetes remains unclear [19]. We found that homozygous carriers of diabetes-associated allele showed the lowest BMI in the control group, which contradicted the theory that weight gain may increase the risk of diabetes. This suggests that the KCNQ1 gene may impact on the pathogenesis of diabetes through pathways other than weight gain. The association of HbA1c levels and fasting plasma glucose levels in our study, and the association of higher fasting glucose levels and reduced corrected insulin response at 120 min in the study of Tan et al. [20], in conjunction with the lower homeostasis model assessment for beta cell function or corrected insulin response apparent for C-C homozygotes of rs2237892 in the study of Yasuda et al. [18], may reflect a functional role for this channel in beta cell function and the pathogenesis of type 2 diabetes. Thus, further studies will be useful to replicate these promising findings and to fully delineate the role of KCNQ1.
In the present study, the distribution of haplotypes constructed by rs2237892, rs2237895 and rs2237897 showed significant differences between cases and controls, particularly the C-C-C haplotype and the T-A-T haplotype. Thus, it can be seen that the three variants are in an LD block. However, no variant in this region has been detected to have a direct function, so we suggest that there may be one or more other functional variants in this region of KCNQ1. In a recent study, large gene deletions and duplications in KCNQ1 were identified in people with long QT syndrome [28]. This suggests that copy-number variants in KCNQ1 also play a role in associated complex diseases. To clarify the role of KCNQ1 in susceptibility to type 2 diabetes, further research is needed into the SNPs and copy-number variants.
In addition, compared with the controls in previous studies (see ESM Table 5), the minor allele frequencies of all three variants in our population from mainland China were higher or lower than in other populations; for example, the frequency of SNP rs2237897 in our Chinese sample was close to that in the Singapore Chinese population, higher than that in the Japanese population and completely different from the frequency in the Danish population (63.8%, 65%, 61% and 96%, respectively) [18–20]. These results confirmed the differences in frequencies of some genetic variations among different ethnic groups and different geographic regions. At the same time, these differences also highlight the need to extend genome-wide association scans and subsequent replication studies to different populations.
In summary, we detected the association of KCNQ1 with type 2 diabetes and some diabetes-related traits in the population of mainland China. In addition, the haplotype distribution was found to be significantly different between cases and controls. These findings may be useful in further research into the pathogenesis of type 2 diabetes, as well as in advancing clinical practice and public health genomics.
Abbreviations
- LD:
-
Linkage disequilibrium
- PAR:
-
Population attributable risk
- SNP:
-
Single nucleotide polymorphism
References
Wild S, Roglic G, Green A, Sicree R, King H (2004) Global prevalence of diabetes: estimates for the year 2000 and projections for 2030. Diabetes Care 27:1047–1053
Grant SF, Thorleifsson G, Reynisdottir I et al (2006) Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nat Genet 38:320–323
Gloyn AL, Weedon MN, Owen KR et al (2003) Large-scale association studies of variants in genes encoding the pancreatic beta-cell KATP channel subunits Kir6.2 (KCNJ11) and SUR1 (ABCC8) confirm that the KCNJ11 E23K variant is associated with type 2 diabetes. Diabetes 52:568–572
Altshuler D, Hirschhorn JN, Klannemark M et al (2000) The common PPARgamma Pro12Ala polymorphism is associated with decreased risk of type 2 diabetes. Nat Genet 26:76–80
Zeggini E, Weedon MN, Lindgren CM et al (2007) Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science 316:1336–1341
WTCCC (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447:661–678
Steinthorsdottir V, Thorleifsson G, Reynisdottir I et al (2007) A variant in CDKAL1 influences insulin response and risk of type 2 diabetes. Nat Genet 39:770–775
Sladek R, Rocheleau G, Rung J et al (2007) A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature 445:881–885
Scott LJ, Mohlke KL, Bonnycastle LL et al (2007) A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science 316:1341–1345
Saxena R, Voight BF, Lyssenko V et al (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316:1331–1336
Wu Y, Li H, Loos RJ et al (2008) Common variants in CDKAL1, CDKN2A/B, IGF2BP2, SLC30A8, and HHEX/IDE genes are associated with type 2 diabetes and impaired fasting glucose in a Chinese Han population. Diabetes 57:2834–2842
Ng MC, Park KS, Oh B et al (2008) Implication of genetic variants near TCF7L2, SLC30A8, HHEX, CDKAL1, CDKN2A/B, IGF2BP2, and FTO in type 2 diabetes and obesity in 6, 719 Asians. Diabetes 57:2226–2233
Lee YH, Kang ES, Kim SH et al (2008) Association between polymorphisms in SLC30A8, HHEX, CDKN2A/B, IGF2BP2, FTO, WFS1, CDKAL1, KCNQ1 and type 2 diabetes in the Korean population. J Hum Genet 53:11–12
Rong R, Hanson RL, Ortiz D et al (2008) 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 58:478–488
Hertel JK, Johansson S, Raeder H et al (2008) Genetic analysis of recently identified type 2 diabetes loci in 1,638 unselected patients with type 2 diabetes and 1,858 control participants from a Norwegian population-based cohort (the HUNT study). Diabetologia 51:971–977
Liu Y, Yu L, Zhang D et al (2008) Positive association between variations in CDKAL1 and type 2 diabetes in Han Chinese individuals. Diabetologia 51:2134–2137
Cauchi S, Meyre D, Durand E et al (2008) Post genome-wide association studies of novel genes associated with type 2 diabetes show gene–gene interaction and high predictive value. PLoS ONE 3:e2031
Yasuda K, Miyake K, Horikawa Y et al (2008) Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet 40:1092–1097
Unoki H, Takahashi A, Kawaguchi T et al (2008) SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations. Nat Genet 40:1098–1102
Tan JT, Nurbaya S, Gardner D, Sandra Y, Tai ES, Ng DP (2009) Genetic variation in KCNQ1 associates with fasting glucose and beta-cell function: a study of 3734 subjects comprising three ethnicities living in Singapore. Diabetes. doi:10.2337/db08-1138
Huang K, Tang W, Tang R et al (2008) Positive association between OLIG2 and schizophrenia in the Chinese Han population. Hum Genet 122:659–660
Shi YY, He L (2005) SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res 15:97–98
Bruzzi P, Green SB, Byar DP, Brinton LA, Schairer C (1985) Estimating the population attributable risk for multiple risk factors using case–control data. Am J Epidemiol 122:904–914
Nawata H, Shirasawa S, Nakashima N et al (2004) Genome-wide linkage analysis of type 2 diabetes mellitus reconfirms the susceptibility locus on 11p13–p12 in Japanese. J Hum Genet 49:629–634
Mori Y, Otabe S, Dina C et al (2002) Genome-wide search for type 2 diabetes in Japanese affected sib-pairs confirms susceptibility genes on 3q, 15q, and 20q and identifies two new candidate loci on 7p and 11p. Diabetes 51:1247–1255
Lee MP, Ravenel JD, Hu RJ et al (2000) Targeted disruption of the Kvlqt1 gene causes deafness and gastric hyperplasia in mice. J Clin Invest 106:1447–1455
Casimiro MC, Knollmann BC, Ebert SN et al (2001) Targeted disruption of the Kcnq1 gene produces a mouse model of Jervell and Lange–Nielsen Syndrome. Proc Natl Acad Sci USA 98:2526–2531
Eddy CA, MacCormick JM, Chung SK et al (2008) Identification of large gene deletions and duplications in KCNQ1 and KCNH2 in patients with long QT syndrome. Heart Rhythm 5:1275–1281
Acknowledgements
We thank the individuals who participated in the present study. This work was supported by grants from the National Key Technology R&D Program (2006BAI05A05), S973 Program (2007CB947300), 863 Program (2006AA02A407), S973 Program (2006CB910601), the Opening Project of Shanghai Key Laboratory of Complex Prescription (07DZ22917), the Shanghai–Unilever Research and Development Fund (06SU07007), the Shanghai Leading Academic Discipline Project (B205), the Chinese Nutrition Society (05015), the Shanghai Municipality Science & Technology Commission (05JC14090) and the Chinese Academy of Sciences (KSCX2-YW-R-01).
Duality of interest statement
The authors declare that there is no duality of interest associated with this manuscript.
Open Access
This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Author information
Authors and Affiliations
Corresponding authors
Additional information
Y. Liu and D. Z. Zhou contributed equally to this study.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Table S1
Genetic mode selection test for variants of KCNQ1 (PDF 31 kb)
Table S2
LD coefficients (D′ and r 2) among SNPs in KCNQ1 (PDF 26 kb)
Table S3
Estimated haplotype frequencies and association significance (PDF 51 kb)
Table S4
Tests of the independent effect of each SNP in KCNQ1 (PDF 26 kb)
Table S5
Genetic differences among different control populations (PDF 8 kb)
Rights and permissions
Open Access This is an open access article distributed under the terms of the Creative Commons Attribution Noncommercial License (https://creativecommons.org/licenses/by-nc/2.0), which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
About this article
Cite this article
Liu, Y., Zhou, D.Z., Zhang, D. et al. Variants in KCNQ1 are associated with susceptibility to type 2 diabetes in the population of mainland China. Diabetologia 52, 1315–1321 (2009). https://doi.org/10.1007/s00125-009-1375-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00125-009-1375-y