Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

An Exploratory Study of the Association between KCNB1 rs1051295 and Type 2 Diabetes and Its Related Traits in Chinese Han Population

  • Yu-Xiang Zhang ,

    Contributed equally to this work with: Yu-Xiang Zhang, Yan Liu, Jing Dong

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Yan Liu ,

    Contributed equally to this work with: Yu-Xiang Zhang, Yan Liu, Jing Dong

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Infection Control Office, Peking University People’s Hospital, Beijing, China

  • Jing Dong ,

    Contributed equally to this work with: Yu-Xiang Zhang, Yan Liu, Jing Dong

    Affiliation Health Medical Center, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • You-Xin Wang,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Jing Wang,

    Affiliation Department of Endocrinology, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • Guo-Qing Zhuang,

    Affiliation College of Life Science, Graduate University of Chinese Academy of Sciences, Beijing, China

  • Shu-Jing Han,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Qing-Qing Guo,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Yan-Xia Luo,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Jie Zhang,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Xiao-Xia Peng,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Ling Zhang,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Yu-Xiang Yan,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Xing-hua Yang,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Hong Wang,

    Affiliation Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China

  • Xu Han,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • Guang-Xu Liu,

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • You-Hou Kang,

    Affiliation Departments of Medicine and Physiology, University of Toronto, Toronto, Ontario, Canada

  • You-Qin Liu,

    Affiliation Health Medical Center, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • Sheng-Feng Weng,

    Affiliation Department of Clinical Laboratory, Beijing Geriatric Hospital, Beijing, China

  • Hong Zhang,

    Affiliation Department of Clinical Laboratory, Beijing Geriatric Hospital, Beijing, China

  • Xiao-Qiang Zhang,

    Affiliation Department of Clinical Laboratory, Beijing Geriatric Hospital, Beijing, China

  • Ke-Bao Jia,

    Affiliation Department of Endocrinology, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • Li Wang,

    Affiliation Department of Endocrinology, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • Lei Zhao,

    Affiliation Department of Endocrinology, Beijing Xuanwu Hospital, Capital Medical University, Beijing, China

  • Zhong-Xin Xiao,

    Affiliation Experimental Teaching Center, Capital Medical University, Beijing, China

  • Shu-Hua Zhang,

    Affiliation Experimental Teaching Center, Capital Medical University, Beijing, China

  • Hui-Hui Wu,

    Affiliation Experimental Teaching Center, Capital Medical University, Beijing, China

  • Qing-Xuan Lai,

    Affiliation College of Life Science, Graduate University of Chinese Academy of Sciences, Beijing, China

  • Na Qi,

    Affiliation Center for Clinical Laboratory, Capital Medical University, Beijing, China

  • Wei Wang,

    Affiliation College of Life Science, Graduate University of Chinese Academy of Sciences, Beijing, China

  • Herbert Gaisano ,

    herbert.gaisano@utoronto.ca (HG); liufen05@ccmu.edu.cn (FL); yanhe118@sina.com (YH)

    Affiliation Departments of Medicine and Physiology, University of Toronto, Toronto, Ontario, Canada

  • Fen Liu ,

    herbert.gaisano@utoronto.ca (HG); liufen05@ccmu.edu.cn (FL); yanhe118@sina.com (YH)

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  •  [ ... ],
  • Yan He

    herbert.gaisano@utoronto.ca (HG); liufen05@ccmu.edu.cn (FL); yanhe118@sina.com (YH)

    Affiliations Department of Epidemiology and Health Statistics, School of Public Health and Family Medicine, Capital Medical University, Beijing, China, Beijing Municipal Key Laboratory of Clinical Epidemiology, Beijing, China

  • [ view all ]
  • [ view less ]

Abstract

Since the KCNB1 encoding Kv2.1 channel accounts for the majority of Kv currents modulating insulin secretion by pancreatic islet beta-cells, we postulated that KCNB1 is a plausible candidate gene for genetic variation contributing to the variable compensatory secretory function of beta-cells in type-2 diabetes (T2D). We conducted two studies, a case-control study and a cross-section study, to investigate the association of common single-nucleotide polymorphisms (SNPs) in KCNB1 with T2D and its linking traits. In the case-control study, we first examined the association of 20 tag SNPs of KCNB1 with T2D in a population with 226 T2D patients and non-diabetic subjects (screening study). We then identified the association in an enlarged population of 412 T2D patients and non-diabetic subjects (replication study). In the cross-sectional study, we investigated the linkage between the candidate SNP rs1051295 and T2D by comparing beta-cell function and insulin sensitivity among rs1051295 genotypes in a general population of 1051 subjects at fasting and after glucose loading (oral glucose tolerance tests, OGTT) in 84 fasting glucose impaired subjects, and several T2D-related traits. We found that among the 19 available tag SNPs, only the KCNB1 rs1051295 was associated with T2D (P = 0.027), with the rs1051295 TT genotype associated with an increased risk of T2D compared with genotypes CC (P = 0.009). At fasting, rs1051295 genotype TT was associated with a 9.8% reduction in insulin sensitivity compared to CC (P = 0.008); along with increased plasma triglycerides (TG) levels (TT/CC: P = 0.046) and increased waist/hip (W/H) ratio (TT/CC: P = 0.013; TT/TC: P = 0.002). OGTT confirmed that genotype TT exhibited reduced insulin sensitivity by 16.3% (P = 0.030) compared with genotype TC+CC in a fasting glucose impaired population. The KCNB1 rs1051295 genotype TT in the Chinese Han population is associated with decreased insulin sensitivity and increased plasma TG and W/H ratio, which together contribute to an increased risk for T2D.

Introduction

Although the precise mechanisms underlying the development and progression of T2D remain unclear, a consensus is that pathophysiologic defects underlying T2D include insulin resistance of peripheral tissues and defects in pancreatic islet insulin secretory capacity [1], each influenced by environmental and genetic factors [2]. With respect to the latter, numerous genome-wide association studies (GWAS) have been conducted leading to identification of susceptibility loci for T2D. SNPs of ion channel genes contributing to beta-cell secretory defects in T2D include ATP-sensitive K+ (KATP) channels (KCNJ11) [3] and Ca2+ channels (CACNA1E) [4]. KATP channel closure causes cell depolarization which opens Ca2+ channels to enable Ca2+ influx that evokes insulin exocytosis, and thus defects in these ion channels would result in insulin secretory insufficiency [5]. Most recently, studies revealed SNPs in KCNQ1, encoding Kv7.1, associated with T2D in the Japanese population [6], and confirmed to be present in Chinese [7], [8], Koreans [6], and also Swedes [6] and Danes [9]. Furthermore, common variants of KCNQ1 have been shown to be associated with reduced insulin granule docking and depolarization-evoked insulin exocytosis [10], and impairment in insulin secretion during glucose loading [11].

Voltage-gated K+ (Kv) channels regulate cell membrane repolarisation that controls duration of Ca2+ channel opening [12], which in beta-cells influences duration of insulin secretion [13], [14]. Kv7.1 is however not the major Kv channel in pancreatic islets, but rather in the heart, where genetic defects account for several cardiac arrhythmic disorders [15], [16]. The dominant Kv channel accounting for ∼70% of membrane repolarisation in beta-cells is Kv2.1 [13], [14], encoded by KCNB1 located in chromosome 20q13.2 [17]. The Kv2.1 channel is also expressed in brain, atria, ventricle, skeletal muscle and other tissues [17]. Based on our previous work on the beta-cell Kv2.1 channel identifying its important role in modulating of insulin secretion [18], [19], we pursued the possibility and hypothesis that there may be a genetic variation in KCNB1 associated with T2D that could influence disease progression and/or compensatory capacity. This prompted us to search for candidate KCNB1 SNPs associated with T2D, employing a case-control study, followed by a cross-section study to examine underlying type 2 diabetic-related traits linking this association in the Chinese Han population.

Materials and Methods

This study was approved by the Ethics Committees of Capital Medical University (Beijing, China) and also the Beijing Geriatric Hospital and Beijing Xuanwu Hospital, and was conducted in accordance with the principles of the Helsinki Declaration II. Written consent was obtained from all participants.

Study Participants and Study Design

In the case–control study, the participants were composed of 412 Chinese Han participants. This included 176 unrelated individuals with T2D designated as Cases, identified by 1999 WHO criteria [20]. These subjects were recruited consecutively from the Department of Endocrinology of two Beijing hospitals, Xuanwu Hospital Capital Medical University and Beijing Geriatric Hospital. 236 unrelated non-diabetic individuals were recruited consecutively, designated as Controls, from the Departments of Otolaryngology and Ophthalmology of the above two hospitals during the same period as cases. These controls were identified by having a fasting plasma glucose level of <6.1 mmol/L and having not been diagnosed as T1D or T2D.

For the first stage of the case-control study, a KCNB1 gene-wide tag SNPs-based association study (screening study) was performed, aimed at screening for association between the presence of KCNB1 SNPs and T2D. 226 participants were randomly selected from the above 412 participants, which included 112 individuals with T2D as cases and 114 non-diabetic individuals as controls. In the initial screening study, the KCNB1 3′-UTR rs1051295 (T>C) was found to have a trend towards association with T2D, which prompted us to perform a replication study aimed at identifying a genuine association of the KCNB1 rs1051295 and T2D. The participants included all 412 Chinese Han participants.

Then in the cross-section study, we examined the underlying mechanism linking the association of KCNB1 rs1051295 and T2D. 1051 individuals participating in annual health examinations were consecutively recruited at the Physical Examination Center, Beijing Xuanwu Hospital. Here we tested for the association of some T2D-related traits with KCNB1 rs1051295 at basal condition. Individuals with previously diagnosed T1D or T2D or with fasting plasma glucose of ≥6.1 mmol/L were excluded.

We also investigated whether KCNB1 rs1051295 is associated with two major pathogenic contributing factors to T2D: beta-cell compensatory secretory function and peripheral insulin sensitivity at glucose load. 84 volunteer outpatients with fasting glucose between 6.1 and 7.0 mmol/L who willing to undertake an OGTT, were recruited consecutively from the Department of Endocrinology of the Beijing Xuanwu Hospital. For the OGTT study, participants drank a 75 g glucose solution, and blood samples were drawn at 0, 30, 60 and 120 minutes and sent for glucose and insulin determination, and genotyping.

From all participants, blood samples were drawn for DNA preparation and determination of fasting plasma glucose, fasting insulin, triglycerides (TG), total cholesterol (TCH), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C). A physical examination was performed on all participants to obtain their height, weight, waist circumference, hip circumference, and blood pressure.

Genotyping, PCR and Quality Control

We used genotyping data from the CHB panel (Han Chinese in Beijing) of the phase II HapMap Project to select 20 tag SNPs within the range of 5000 bp upstream of the initiation codon and 5000 bp downstream of the termination codon of KCNB1 on the basis of the following principal criteria: r2>0.8 and minor allele frequency (MAF)>0.1. Of those tag SNPs, 19 SNPs were located in introns, and 1 was in the 3′-UTR. Detailed information for all the tag SNPs (including genomic position, genic position, minor allele frequency, Hardy-Weinberg equilibrium) are shown in Table1.

Genomic DNA was extracted from peripheral white blood cells of participants using a DNeasy tissue kit (Qiagen) according to the manufacturer’s instructions. All SNPs were genotyped using MassArray (Sequenom, San Diego, CA) based on allele-specific MALDI-TOF mass spectrometry. DNA from cases and controls was randomly assigned to 96 well plates, and blinded genotyping was performed. The call rates for the genotyping of the SNPs were all >90% except for rs802952, for which the rate was 39.58%. Thus rs802952 was excluded from further analysis.

For the replication study and the cross section study, rs1051295 was genotyped using PCR-based pyro-sequencing technology and DNA sequencing with an ABI 3730 automated sequencer (Applied Biosystems, Foster city, CA). The primers are F: biotin-GCGCAAAACCCTTACTCA AAT and R: GCCAGGGGGCAATTAGAAT for PCR amplification and 5′- TGGTATCTCAA AATAAAATC-3′ for sequencing. The primers were designed according to the published sequence of KCNB1 (http://www.ncbi.nlm.nih.gov/pubmed/), the GenBank accession number is NT_002370.

PCR was conducted in a 50-µl reaction mixture containing 100 ng genomic DNA, 0.5 pmol primer, 2× master mix (mixture of reaction buffer, MgCl2 and DNA polymerase) (ToYoBo, Japan). PCR products were denatured for 2 min at 95°C and then thermal-cycled for 30 s at 95°C, 30 s at 59°C, and 60 s at 72°C, repeating the cycle 40 times. A final extension step at 72°C for 10 min completed the program. Pyrosequencing analysis was performed on Streptavidin Sepharose™ HP (Amersham, Sweden). ssDNA prepared from 50 ul of biotinylated PCR product was then subjected to pyrosequencing. The call rates was 96.02% (169/176) and 94.92% (224/236) for cases and control, respectively, for the replication study, and 97.91% (1029/1051) and 100% (84/84), respectively, for the first and second stages of the cross-section study.

For both the case-control and the cross-section studies, genotyping was repeated in 5% of random samples for verification and quality control, which all revealed the genotype data had an error rate of <1%.

Assessment of Islet Beta-cell Secretory Function and Insulin Sensitivity

We employed the homeostatic model to calculate beta-cell function (HOMA-B%) and insulin sensitivity (HOMA-S%) at basal condition (the software HOMA calculator downloaded from www.dtu.ox.ac.uk/homa). In the glucose loading condition (OGTT), insulin sensitivity was assessed by the Matsuda index = 10,000/. Ins0 is fasting insulin, Glu0 is fasting glucose. The islet beta-cell secretory function was assessed by the Insulin Secretion-Sensitivity Index-2 (ISSI-2) = AUC (insulin curve)/AUC (glucose curve)* Matsuda index. The Matsuda index is a validated OGTT-based measure of insulin sensitivity that is analogous to the disposition index obtained from iv glucose tolerance tests [21]. Area under curve (AUC) for glucose and insulin were calculated by the trapezoidal rule.

Statistical Analyses

In the case-control study, the χ2 test was used to examine the differences in gender between case and control groups. The differences between T2D-related traits in T2D case and control groups were compared by an independent–sample t test. The distribution of SNPs in Cases and Controls were compared using the χ2 test and an unconditional Logistic regression analysis, in which the association of a SNP with T2D was adjusted for age, gender and BMI. T2D-related traits among the three genotypes of rs1051295 were compared by one-way ANOVA test.

The program Haploview (http://www.broad.mit.edu/mpg/haploview/) was used to calculate pair-wise linkage disequilibrium statistics and to test allelic and haplotype associations with T2D.

In the cross-section study, we used independent-samples t test to compare T2D-related traits between genotypes (TT vs TC+CC) in OGTT, or one-way ANOVA test among 3 genotypes at fasting. Association of rs1051295 genotypes with HOMA-B%, HOMA-S%, ISSI-2, Matsuda index and other T2D-related traits were examined by multiple linear regression analysis, in which association was adjusted for age, gender, and BMI.

Hardy–Weinberg equilibrium was determined for each SNP distribution. All analysis was done on SPSS software, version 18.0 (purchased by Capital Medical University, China). All tests were two-tailed, with a significance level of 0.05. Data were expressed as means ± SD, unless the data did not conform to a normal distribution, in which case the data were expressed as median and quartiles and natural Log-transformed for analysis.

Results

Candidate-gene Association Study Identified KCNB1 3′-UTR rs1051295 is Likely to be Associated with T2D

The demographic and clinical characteristics of the participants for the case-control study are summarized in Table 2. Among the 19 available tag SNPs, the distribution of rs742759 (P = 0.16) and rs1051295 (P = 0.18) showed a possible association with T2D. The other 17 tag SNPs did not show a trend of association with T2D (P>0.20, data not shown). However, only rs1051295 had a trend towards association with T2D (P = 0.08) when adjusted for age, gender and BMI (Table 3).

thumbnail
Table 2. Characteristics of the participants in the case-control study.

https://doi.org/10.1371/journal.pone.0056365.t002

thumbnail
Table 3. Distribution of two tag SNPs of KCNB1 in the screening study.

https://doi.org/10.1371/journal.pone.0056365.t003

To determine whether these SNPs demonstrated any additional evidence of association with T2D when examined together, haplotypes for 19 SNPs were constructed using the Haploview program. We identified four haplotypes for the KCNB1 gene, and found no association between any of them and T2D (data not shown). It is noteworthy that rs1051295 was not in a haplotype block generated in this study due to weak linkage disequilibrium with other SNPs.

In order to identify this association, a replication study was performed in the source population of the above screening study. This population consisted of 176 cases and 236 controls. Hardy–Weinberg equilibrium testing showed P = 0.06 and 0.14 for rs1051295 genotypes in cases and controls, respectively. TT genotype of rs1051295 demonstrated a significantly increased risk for T2D compared with the CC (P = 0.009, OR = 2.58, 95% CI = (1.27, 5.23)), but not with TC (P = 0.226, OR = 1.36, 95% CI = (0.83, 2.25)) genotypes. In the dominant model, allele C (genotype CC and TC) did not show decreased risk for T2D compared with genotype TT (P = 0.071, OR = 0.64, 95% CI = (0.40, 1.04)); but in the recessive model, allele T(genotype TT+TC) showed increased risk for T2D compared with genotype CC(P = 0.02, OR = 2.10, 95% CI = 1.13,3.90) (Table 4).

thumbnail
Table 4. Distribution of KCNB1 rs1051295 SNPs in the replication study.

https://doi.org/10.1371/journal.pone.0056365.t004

We further compared a number of T2D-related traits among the three variants of rs1051295 including BMI, W/H ratio, TG, fasting plasma glucose, SBP, DBP as well as gender and age, but found no association of any genotype with these traits (Table S1 and Table S2).

KCNB1 rs1051295 Associated with T2D-related Traits Under Basal Condition and OGTT

Since the pathophysiology underlying T2D is characterized by defective insulin secretion and reduced insulin sensitivity, a cross-section study was designed to determine beta-cell secretory function and insulin sensitivity in the three KCNB1 rs1051295 variants under basal (fasting) condition to assess which of these traits could account for their differential risk for T2D. Table 5 shows that the insulin sensitivity of genotype TT (90.1%) was reduced by 9.8% compared to genotype CC (99.9%), which was a significant difference (P = 0.008, b = −0.09, 95% CI (−0.16, −0.02)) after adjusting for age, gender and BMI. Comparing TT to TC (95.2%), the difference was 5.1%, which did not reach statistical significance (P = 0.09, b = −0.04, 95% CI (−0.09, 0.01)). When we combined genotype CC and TC as one group of allele C and compared its insulin sensitivity (HOMA-S%) with that of genotype TT in the dominant model, the result showed genotype TT decreased insulin sensitivity (P = 0.05, b = −5.22, 95% CI = (−10.49, −0.00)), whereas allele T (genotype TC and TT) did not show decreased insulin sensitivity compared with genotype CC in the recessive model (P = 0.08, b = −5.87, 95% CI = (−12.53,0.79)).

thumbnail
Table 5. Comparison of T2D-related traits among rs1051295 genotypes in a general population at fasting.

https://doi.org/10.1371/journal.pone.0056365.t005

Importantly, the rs1051295 genotype TT was associated with an unfavorable W/H ratio and higher plasma TG levels compared to genotype CC (W/H ratio: P = 0.01, b = 0.01, 95% CI (0.00,0.02); TG: P = 0.04, b = 0.08, 95% CI (0.00, 0.16)) and TC (W/H ratio: P = 0.002, b = 0.01, 95% CI (0.00, 0.02); TG: P = 0.09, b = 0.05, 95% CI (−0.01, 0.11)) after adjustment for age, gender and BMI (Table 5). These results indicated that the reduced insulin sensitivity in rs1051295 genotype TT might at least be associated with the higher TG and more unfavorable W/H ratio.

Table 5 demonstrates that rs1051295 genotype TT exhibited increased fasting insulin level (P = 0.004, b = 0.86, 95% CI (0.27, 1.44)) and beta-cell secretory function (P = 0.004, b = 0.08, 95% CI (0.03, 0.13)) compared with genotype CC at basal condition. To investigate whether the insulin secretory function of beta-cell in the carrier of genotype TT is superior to genotype CC, we performed OGTT to assess the compensatory secretory function and insulin sensitivity of beta-cells after a glucose challenge. Because of the small size of the groups, and genotype CC and TC demonstrated the similar phenotype in insulin sensitivity in the dominant model at basal glucose, we combined CC and TC as one to compare its insulin sensitivity with genotype TT.

OGTT (Table 6) showed there are no significant difference between beta-cell function (ISSI-2) in genotype TT (32.23) and TC+CC (37.54, P = 0.13, b = −0.29, 95% CI (−0.67, 0.09)). OGTT also demonstrated genotype TT displayed decreased insulin sensitivity (Matsuda index: 5.05) by 16.3% compared to CC+TC (Matsuda index: 6.03) (P = 0.03, b = −1.34, 95% CI (−2.56, −0.12)). This is most dramatically seen at 120min where insulin secretion from genotype TT (76.45 mmol/L) was 28.2% higher than genotype CC+TC (54.84 mmol/L), whereas blood glucose levels at this time point were still 14.4% higher in genotype TT (9.9 mmol/L) than CC+TC (8.5 mmol/L) [Fig. 1]. These results indicate that subjects with rs1051295 genotype CC/TC exhibited greater insulin sensitivity than genotype TT during a glucose load and confirmed the finding at basal condition.

thumbnail
Figure 1. Relationship between plasma insulin and glucose levels during OGTT (0, 30th, 60th and 120th min).

n = 23 for rs1051295 genotype TT and 61 for rs1051295 genotypes CC+TC.

https://doi.org/10.1371/journal.pone.0056365.g001

thumbnail
Table 6. Comparison of T2D-related traits between rs1051295 genotype CC+TC and genotype TT in a type 2 diabetic-suspected population undergoing an OGTT.

https://doi.org/10.1371/journal.pone.0056365.t006

Discussion

In summary, we unexpectedly found that KCNB1 rs1051295 genotype TT was associated with decreased insulin sensitivity at basal conditions in a general population. Of the three variants in KCNB1 rs1051295, CC (99.9%) and to a less extent of TC (95.2%) exhibited normal insulin sensitivity at basal condition, whereas TT genotype (90.1%), exhibited reduced insulin sensitivity at the fasting state compared to CC (P = 0.008, b = −0.09, 95% CI (−0.16, −0.02)) and TC (P = 0.096, b = −0.04, 95% CI (−0.09, 0.01)). genotype TT also decreased insulin sensitivity compared with genotype(CC+TC) (allele C)(P = 0.05, b = −5.22, 95% CI = (−10.49, −0.00)). This was also confirmed in a fasting glucose impaired population at glucose loading condition. In the OGTT, the TT genotype also exhibited lower insulin sensitivity compared with genotype (TC+CC) (P = 0.03, b = −1.34, 95% CI (−2.56, −0.12)). Collectively, these results indicate that KCNB1 rs1051295 TT confers on its carriers the phenotype of decreased insulin sensitivity that is likely to increase the risk of T2D.

Kv2.1, encoded by KCNB1, is the major beta-cell Kv channel in humans [13], [22] and rodents [13], [14] accounting for >70% of outward K+ currents. Kv2.1 deletion in mice caused severe perturbation in insulin release and blood glucose [14]. It is reasonable to expect the rs1051295 genotype to be associated with an increased the risk of T2D by impairing islet beta-cell secretory function. We compared the beta-cell function among rs1051295 genotypes at basal condition, whereby HOMA-B% (in vivo assessment of beta-cell function) paradoxically showed a higher value in genotype TT than CC (P = 0.004, b = 0.08, 95% CI (0.03, 0.13)).This should not be interpreted to represent a superior beta cell secretory function of the TT genotype, but rather, the higher HOMA-B% in the TT genotype indicates a compensatory action of the pancreatic islet in response to the reduced insulin sensitivity. To confirm this thinking, we compared the beta-cell function under glucose loading condition, and employed the ISSI-2 assessment as a more accurate measurement of beta-cell compensatory insulin secretory function. ISS1-2 assessment showed the KCNB1 rs1051295 CC and TC genotypes was 16.31% higher than TT (Table 6) (P = 0.13, b = −0.29, 95% CI (−0.67, 0.09)), suggesting that the TT genotype might have actually reduced insulin secretory capacity, at least with increased glycemic demand during OGTT, and presumably after a meal.

In this study we unexpectedly found that insulin sensitivity was significantly reduced in the KCNB1 rs1051295 genotype TT compared with genotypes CC and TC. Major insulin sensitive tissues are skeletal muscle, adipose tissue and liver. While Kv2.1 is known to be present in skeletal muscle to regulate membrane excitability [23], it is not known to regulate GLUT4 vesicle transport and exocytosis. Recent work demonstrating that Kv2.1 directly interacts with exocytotic SNARE proteins to modulate exocytosis [19], [24], [25] raises the possibility of non-channel exocytotic function(s) of Kv2.1 in skeletal muscle regulating GLUT4 transport or function. Kv2.1 is abundant in neuronal tissue [12] which might influence insulin-sensitive tissues.

We also showed that rs1051295 genotype TT is associated with higher TG levels and affected fat distribution towards an unfavorable waist/hip ratio and abdominal obesity. It is not clear how Kv2.1 could affect adipose tissue metabolism and distribution. We demonstrated the mean HDL-C (TT: 1.65 versus CC: 1.67, P = 0.47) was lower and LDL-C (TT: 2.72 versus CC: 2.63, P = 0.46) was higher in the population with a TT genotype than those with a CC genotype, which may provide some clues for future investigation on how rs1051295 genotype TT might influence fat metabolism. In support of a role of Kv channels in adipocyte metabolism, Kv channel activity has been recorded in human adipocytes [26], which responded to insulin by increasing channel density [27]. However, the specific Kv isoform or its downstream actions in adipocytes have not been critically assessed, thus requiring much further study. Rs1051295 is located in the 3′-UTR region of KCNB1, and has been found implicated in rheumatoid arthritis [28]. Since this SNP is not located in the coding regions, we speculate that it might not influence channel pore kinetics per se, but rather it may influence the post translational processing of the Kv2.1 channel. Further study should be aimed at elucidating at identifying the target tissues that rs1051295 exert its effect on insulin sensitivity.

This work has thus identified KCNB1 rs1051295 genotype TT to be associated with reduced insulin sensitivity and increased plasma TG and W/H ratio, these together likely leading to increase the risk of T2D. Further confirmation will be needed for this study to be performed in a larger population and also to determine if these results could be replicated in other (non-Han) populations. Nonetheless, this initial study opens a new avenue to explore a possible role of Kv2.1 and other Kv channel in influencing fat metabolism and insulin sensitivity.

Supporting Information

Table S1.

T2D-related quantitative traits in different genotypes of KCNB1 rs1051295 in 226 case-control study.

https://doi.org/10.1371/journal.pone.0056365.s001

(DOC)

Table S2.

T2D-related quantitative traits in different genotypes of KCNB1 rs1051295 in 412 validation study.

https://doi.org/10.1371/journal.pone.0056365.s002

(DOC)

Acknowledgments

The authors thank Robert Cormier (University of Minnesota Medical School at Duluth) and Ravi Retnakaran (Mount Sinai Hospital, Toronto) for their help with the manuscript preparation. We are grateful to sample donors for contributing to this research.

Author Contributions

Conceived and designed the experiments: YH FL. Performed the experiments: Y-XZ YL S-JH X-XP Z-XX Q-XL NQ Y-XW FL G-QZ Y-HK Y-XY Q-QG. Analyzed the data: Y-XZ FL Y-XL LZ XH G-XL X-HY WW. Contributed reagents/materials/analysis tools: YL JD JW HW Y-QL K-BJ LW LZ S-FW HZ X-QZ S-HZ JZ H-HW. Wrote the paper: YH HG.

References

  1. 1. Kahn BB (1998) Type 2 diabetes: when insulin secretion fails to compensate for insulin resistance. Cell 92: 593–596.
  2. 2. O’Rahilly S (2009) Human genetics illuminates the paths to metabolic disease. Nature 462: 307–314.
  3. 3. Holmkvist J, Tojjar D, Almgren P, Lyssenko V, Lindgren CM, et al. (2007) Polymorphisms in the gene encoding the voltage-dependent Ca2+ channel CaV2.3 (CACNA1E) are associated with type 2 diabetes and impaired insulin secretion. Diabetologia 50: 2467–2475.
  4. 4. Gloyn AL, Weedon MN, Owen KR, Turner MJ, Knight BA, et al. (2003) Large scale association studies of variants in genes nencoding 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.
  5. 5. Ashcroft FM, Rorsman P (2004) Molecular defects in insulin secretion in type-2 diabetes. Rev Endocr Metab Disord 5: 135–142.
  6. 6. Yasuda K, Miyake K, Horikawa Y, Hara K, Osawa H, et al. (2008) Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus. Nat Genet 40: 1092–1097.
  7. 7. Qi QB, Li HX, Loos RJF, Liu C, Wu Y, et al. (2009) Common variants in KCNQ1 are associated with type 2 diabetes and impaired fasting glucose in a Chinese Han population. Hum Mol Genet 18: 3508–3515.
  8. 8. Hu C, Wang C, Zhang R, Ma X, Wang J, et al. (2009) Variations in KCNQ1 are associated with type 2 diabetes and beta cell function in a Chinese population. Diabetologia 52: 1322–1325.
  9. 9. Unoki H, Takahashi A, Kawaguchi T, Hara K, Horikoshi M, 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.
  10. 10. Rosengren AH, Braun M, Mahdi T, Andersson SA, Travers ME, et al. (2012) Reduced Insulin Exocytosis in Human Pancreatic β-Cells With Gene Variants Linked to Type 2 Diabetes. Diabetes 61: 1726–33.
  11. 11. van Vliet-Ostaptchouk JV, van Haeften TW, Landman GW, Reiling E, Kleefstra N, et al. (2012) Common variants in the type 2 diabetes KCNQ1 Gene are associated with impairments in insulin secretion during hyperglycaemic glucose clamp. PLoS ONE 7: e32148.
  12. 12. Pongs O (2008) Regulation of excitability by potassium channels. Results Probl Cell Differ 44: 145–161.
  13. 13. MacDonald PE, Wheeler MB (2003) Voltage-dependent K(+) channels in pancreatic beta cells: role, regulation and potential as therapeutic targets. Diabetologia 46: 1046–1062.
  14. 14. Jacobson DA, Kuznetsov A, Lopez JP, Kash S, Ammälä CE, et al. (2007) Kv2.1 ablation alters glucose-induced islet electrical activity, enhancing insulin secretion. Cell Metab 6: 229–235.
  15. 15. Wang Q, Curran ME, Splawski I, Burn TC, Millholland JM, et al. (1996) Positional cloning of a novel potassium channel gene: KVLQT1 mutations cause cardiac arrhythmias. Nat Genet 12: 17–23.
  16. 16. Neyroud N, Tesson F, Denjoy I, Leibovici M, Donger C, et al. (1997) A novel mutation in the potassium channel gene KVLQT1 causes the Jervell and Lange-Nielsen cardioauditory syndrome. Nat Genet 15: 186–189.
  17. 17. Gutman GA, Chandy KG, Grissmer S, Lazdunski M, McKinnon D, et al. (2005) International union of pharmacology. LIII. Nomenclature and molecular relationships of voltage-gated potassium channels. Pharmacol Rev 57: 473–508.
  18. 18. Zhuang GQ, Wu W, Liu F, Ma JL, Luo YX, et al. (2009) SNAP-25 (1–180) enhances insulin secretion by blocking Kv2.1 channels in rat pancreatic islet beta-cells. Biochem Biophys Res Commun 379: 812–816.
  19. 19. Leung YM, Kwan EP, Ng B, Kang Y, Gaisano HY (2007) SNAREing voltage-gated K+ and ATP-sensitive K+ channels: tuning beta-cell excitability with syntaxin-1A and other exocytotic proteins. Endocr Rev 28: 653–663.
  20. 20. Alberti KG, Zimmet PZ (1998) Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med 15: 539–553.
  21. 21. Retnakaran R, Qi Y, Goran MI, Hamilton JK (2009) Evaluation of proposed oral disposition index measures in relation to the actual disposition index. Diabet Med 26: 1198–1203.
  22. 22. Herrington J, Sanchez M, Wunderler D, Yan L, Bugianesi RM, et al. (2005) Biophysical and pharmacological properties of the voltage-gated potassium current of human pancreatic beta-cells. J Physiol 567: 159–175.
  23. 23. Schwetz TA, Norring SA, Ednie AR, Bennett ES (2010) Sialic acids attached to O-glycans modulate voltage-gated potassium channel gating. J Biol Chem 286: 4123–4132.
  24. 24. Feinshreiber L, Singer-Lahat D, Ashery U, Lotan I (2009) Voltage-gated potassium channel as a facilitator of exocytosis. Ann N Y Acad Sci 1152: 87–92.
  25. 25. Feinshreiber L, Singer-Lahat D, Friedrich R, Matti U, Sheinin A, et al. (2010) Non-conducting function of the Kv2.1 channel enables it to recruit vesicles for release in neuroendocrine and nerve cells. J Cell Sci 123: 1940–1947.
  26. 26. Ramiez-Ponce MP, Mateos JC, Bellido JA (2003) Human adipose cells have voltage-dependent potassium currents. J Membr Biol 196: 129–134.
  27. 27. Ramírez-Ponce MP, Mateos JC, Bellido JA (2002) Insulin increases the density of potassium channels in white adipocytes: possible role in adipogenesis. J Endocrinol 174: 299–307.
  28. 28. Xiao X, Zhang Y, Wang K (2009) Association of KCNB1 to rheumatoid arthritis via interaction with HLA-DRB1. BMC Proc 3 Suppl 7S134.