- Open Access
An analysis of the association between a polymorphism of KCNJ11 and diabetic retinopathy in a Chinese Han population
European Journal of Medical Research volume 20, Article number: 3 (2015)
Genome-wide association studies (GWAS) have reported that the polymorphism rs5219 of the potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11) is associated with type 2 diabetes mellitus (T2DM). Given that diabetic retinopathy (DR) is one of the most common microvascular complications of T2DM, GWAS have identified a number of potential susceptibility genes for DR. However, only a fraction of them have been replicated in different studies and show consistent genetic associations with the occurrence of DR. The aim of the present study was to investigate whether common variants of KCNJ11 confer DR in a cohort of the Chinese Han population.
A case–control study of 580 T2DM patients, including 105 T2DM with DR and 475 T2DM without DR was performed. A single nucleotide polymorphism (SNP) of KCNJ11 (rs5219) was genotyped, and its association with DR was explored using a dominant genetic model. Genotyping was performed by iPLEX technology. Univariate and multivariate logistic regression (MLR) analysis controlling for confounders was conducted to evaluate the association between rs5219 and DR.
The A allele frequency of rs5219 was significantly higher in DR patients than that in the patients without DR (49.01% versus 38.68%, P <0.05). We found the minor A allele could increase the risk to develop DR (ORint = 1.58, 95% CI: 1.139 to 2.192 for allele and P = 0.006, ORint = 1.607, 95% CI: 1.267 to 2.038 for genotype and P <0.001) in the Chinese Han population.
Our findings provided evidence that KCNJ11 was associated with DR in Chinese Han patients with T2DM.
Type 2 diabetes mellitus (T2DM) is a multifactorial metabolic syndrome, characterized by insulin resistance and/or pancreatic β-cell dysfunction resulting from both genetic and environmental factors [1,2]. Its prevalence has increased rapidly over the last decade owing to human longevity and a surge of obesity in a lot of countries , which causes great economic impact on individuals, families, and health systems worldwide.
Diabetic retinopathy (DR), one of the most common microvascular complications of T2DM, is a leading cause of acquired visual loss in working-age individuals [4,5]. Blindness develops primarily from either increased DME (diabetic macular edema) or proliferation of new retinal vessels . Major risk factors, such as prolonged duration diabetes, poor control of blood glucose and high blood pressure have been considered to be responsible for the onset and progression of DR [7,8]. Nevertheless, these factors do not account completely for the risk of DR. Familial aggregation and the various incidence attributed to the diverse ethnic backgrounds suggest the existence of an heritable component in DR, independent of glycemic control and duration of T2DM .
Recently, four small-scale GWAS with modest sample sizes conducted in Mexican-American, Chinese, and Caucasian populations indicated that several SNPs or loci had borderline or weak associations with DR in either type 1 or type 2 diabetes mellitus [10-13]. Although several candidate gene studies had reported promising genes of DR, only few of them had been replicated and showed evident genetic associations.
The potassium inwardly rectifying channel, subfamily J, member 11 (KCNJ11) is a pore-forming subunit of ATP-sensitive potassium (KATP) channels in pancreatic β-cells  and plays a critical role in the regulation of insulin secretion . The KATP channels consist of four subunits of the inwardly rectifying potassium channel Kir6.2 and four subunits of the sulfonylurea receptor 1 (SUR1), which are targets of sulfonylurea drugs. Variations of KCNJ11 could contribute to the decreased sensitivity of the ion channel to ATP, leading to more ATP consumption, which further contributes to insulin-release impairment. E23K (rs5219) in the KCNJ11 gene, substituting glutamate for lysine at position 23, is identified as a SNP associated with T2DM susceptibility . In addition, it has been reported that KCNJ11 is associated with the therapeutic response to sulfonylureas because of its regulation function of insulin secretion . Based on the reported association between KCNJ11 polymorphisms and T2DM, we hypothesized that this gene could also be related to the risk of DR.
Up to now, none of the risk genes of T2DM contributing to the susceptibility of DR had been reported. The aim of the present study was to investigate whether the genetic variant (rs5219) of KCNJ11 confers DR in a Chinese Han population with T2DM.
Our study involved 580 Chinese patients who resided in the metropolitan area of Shanghai and had been diagnosed with T2DM. T2DM was diagnosed on the basis of the WHO criteria (1999) . Known subtypes of diabetes were excluded based on antibody measurements and inheritance. Patients with diabetic ketoacidosis or ketonuria also were excluded. All the patients underwent digital non-mydriatic fundus photography, and two certified ophthalmologists diagnosed DR. Patients without DR were selected as controls. All study participants registered in the analysis were recruited from the Endocrinology and Metabolism outpatient clinics at Fudan University Huashan Hospital in Shanghai. All subjects provided written informed consent for participation in the study and donation of samples. The Ethics Committee of Huashan Hospital affiliated with Fudan University approved this protocol.
All participants were interviewed for the documentation of medical histories, medications, regular physical examinations, and laboratory assessment of T2DM risk factors. Physician-obtained systolic and diastolic blood pressure (BP) values were taken on the left arm of the seated participants. All participants underwent a complete hematological examination while fasting, including serum total cholesterol (TC), triglyceride (TG), blood urea nitrogen (BUN), uric acid (UA), and C-peptide (CP) levels that were measured by an enzymatic method with a chemical analyzer (Hitachi 7600–020, Tokyo, Japan). Postprandial plasma glucose (PPG) was measured 2 h after eating. Fasting plasma glucose(FPG) was measured in fasting state. The blood was centrifuged at 3,000 rpm for 10 min for plasma separation and immediately used to measure biomarkers. Serum creatinine (Cr) was measured by radioimmunoassay (Beijing Atom High-Tech Co. Ltd.). FPG and PPG were quantified by the glucose oxidase-peroxidase procedure. Glycated hemoglobin (HbA1c) was estimated by high-pressure liquid chromatography using an analyzer (HLC-723G7, Tosoh Corporation, Japan). The day-to-day and interassay coefficients of variation at the central laboratory in our hospital for all analyses were between 1% and 3%.
Diabetes was defined as a self-reported history of physician-diagnosed T2DM or according to 1999 WHO criteria , which consisted of one of the following: fasting plasma glucose FPG ≥7.0 mmol/L, plasma glucose ≥11.1 mmol/L 2 hours after an oral glucose tolerance test (OGTT), or random plasma glucose ≥11.1 mmol/L. All the T2DM patients were tested for DR using digital non-mydriatic fundus photography and image analysis. Fundus photography was performed at each site following a standardized protocol. Both eyes of each participant were photographed with a 45-degree 6.3-megapixel digital non-mydriatic camera (Canon CR6-45NM, Lake Success, NY), repeated only once if necessary. DR was determined by two independent retinal specialists without knowledge of patient clinical details. The patients were classified according to the presence or absence of DR, regardless of its degree of severity. The duration was defined as the interval between the first diagnosis of diabetes and the time of enrollment in the present study. Age of onset year was the age at which an individual was diagnosed with T2DM for the first time. The clinical characteristics of participants are summarized in Table 1.
Peripheral venous blood samples were collected from all participants and the genomic DNA was isolated from venous blood leukocytes by the conventional phenol/chloroform method. The genetic variant (rs5219) of KCNJ11 was genotyped using iPLEX (Sequenom, San Diego, CA, USA) with detection by the matrix-assisted laser desorption/ionization time-of-flight mass spectrometry platform. The with-DR and without-DR groups were mixed for genotyping. There was a 99.9% genotype concordance rate when duplicated samples were compared across plates.
Continuous variables were detected by the Kolmogorov-Smirnov Test according to whether they followed normal distribution. Variables that were not normally distributed were log-transformed to approximate normal distribution for analysis. Results are described as mean ± SD or median unless stated otherwise. Differences in variables between the with-DR group and the with-DR group were determined by unpaired t-test. Between groups, differences in properties were accessed by χ 2 analysis. Univariate logistic regression was performed to determine variables associated with DR and to estimate confounding factors possibly disturbing the relation of genetic variants to DR. Multivariable logistic regression (MLR) was carried out to control potential confounders for determining independent contribution of variables to DR. In order to better investigate interaction between DR and genetic variants of KCNJ11, we performed two analyses according to the variables of allele and genotype of KCNJ11, respectively. Odds ratios (OR) with 95% confidence intervals (CI) were calculated for the relative risk of genetic variants of KCNJ11 with DR. Results were analyzed using the Statistical Package for Social Sciences for Windows version 16.0 (SPSS, Chicago, IL, USA). Tests were two-sided and a P value of <0.05 was considered significant.
Clinical characteristics of subjects
The baseline clinical characteristics of the 580 subjects were listed in Table 1. There are 176 males and 299 females (mean age, 64.77 ± 10.9 years) in control and 42 males and 63 females (mean age, 64.54 ± 10.33 years) in case. DR patients had more weight than non-DR patients. The variables of FPG, PPG, HbAlc, BUN and Cr levels were significantly higher in cases than in controls (P <0.05 for all). In addition, there was a significantly longer duration of T2DM and earlier onset of T2DM in with-DR group as compared with the without-DR group (P <0.05 for all). Other variables of age, sex, height, weight, SBP, DBP, CP, TC, TG, UA was similar between the two groups (P >0.05 for all). The minor allele (A) frequency of rs5219 was 38.68% and 49.01% in controls and cases, respectively. The prevalence of DR was 14.67% and 23.21% in diabetic patients with G and A allele, respectively (Figure 1). The prevalence of DR was 14.47%, 16.85% and 28.57% in diabetic patients with GG, GA and GG genotype, respectively (Figure 2).
Univariate logistic regression analysis for diabetes
To estimate the association of various clinical factors and DR, univariate logistic regression models were developed to include age, sex, BMI, hypertension, lipid profiles, renal function parameters and SNP (rs5219) (Table 2). The univariate logistic analyses indicated that FPG, PPG, HbAlc, BUN, Cr, duration of T2DM, onset age of T2DM were significantly associated with DR (P <0.05 for all). In subjects with MAF of rs5219 in KCNJ11, the OR for DR was 1.524 (95% CI: 1.123 to 2.069, P = 0.007).
Multiple logistic regression analysis for diabetes
MLR demonstrated that the genetic variant (rs5219) of KCNJ11 remained significantly different between case and control after adjustment for potential confounders (P = 0.006 for allele analysis and P <0.001 for genotype analysis, Table 3). After adjusting for confounding factors, in subjects with minor allele frequency (MAF) of rs5219 in KCNJ11, the OR for DR was 1.580 (95% CI: 1.139 to 2.192).
We conducted a case–control study to evaluate the association between KCNJ11 (rs5219) and T2DM in a cohort of a Chinese Han population with T2DM. To our knowledge, this study is the first to find that the T2DM-susceptibility gene polymorphism (KCNJ11 rs5219) is associated with an increased risk of DR (P = 0.006 for allele analysis and P <0.001 for genotype analysis).
T2DM is a complex and multifactorial disease resulting from a complex interaction between genes and environmental factors, whose microvascular and macrovascular complications are highly prevalent among T2DM patients. DR is one of the most common and specific microvascular complications of T2DM, which is considered to be influenced by hereditary [19,20] and environmental factors [21,22]. Given that DR is one of the leading risk factors and causes of blindness worldwide , the identification of risk factors for DR could lead to a decrease in vision loss associated with T2DM.
Attempting to identify risk genes of DR, numerous studies have been conducted over the past few decades . Using candidate gene analysis and GWAS, investigators have identified a number of potential susceptibility genes of DR. Several biochemical pathways, including formation of advanced glycation end products (AGEs), polyol accumulation, activation of protein kinase C, oxidative stress and so on, have been proposed to mediate the pathogenesis of DR . At present, three candidate genes of DR that have been extensively studied are vascular endothelial growth factor (VEGF), aldose reductase (ALR), and the receptor for advanced glycation end products (RAGE) [24,25].
Protein encoded by KCNJ11 has a crucial function in insulin secretion, thus making it a potential susceptibility gene of T2DM. The rs5219 variant of KCNJ11 has been widely reported to be associated with T2DM in various ethnic populations [26-28]. This gene encodes for the Kir6.2 subunit, which forms the KATP channels with the regulatory subunit sulfonylurea receptor 1 (SUR1). These channels are present in many tissues, including heart, vascular smooth muscles, and vascular endothelial cells  and play important roles in responses to stress , blood pressure regulation  and the physiology and pathophysiology of various tissues by coupling the metabolic state of the cells with cellular electrical activity .
On this basis, interestingly, a strong association is found between the common variant of KCNJ11 (rs5219) with DR in our study (P = 0.006 for allele analysis and P <0.001 for genotype analysis), which has not been reported in previous studies. The associations between variants of KCNJ11 and hypertension were demonstrated in Chinese  and Korean populations , suggesting a possible role of KCNJ11 in macrovascular diseases. As DR is one of microvascular complications of T2DM, a hypothesis that KCNJ11 may have a relationship with DR comes up. Our finding implies that KCNJ11 (rs5219) may disturb DR susceptibility. However, the molecular mechanisms are yet to be expounded.
To date, there were no consistent and definitive genetic associations with DR in the literature, likely due to the complexity and multifactorial etiology of DR: both genetic and environmental factors play critical roles in its development and progression. Most published studies were based on small patient samples, and their case definitions were not based on the widely accepted and standardized criteria. In our study, we used the method of non-mydriatic digital stereoscopic retinal imaging (NMDSRI) for screening and diagnosing DR. NMDSRI is a sensitive and specific method for the screening and diagnosis of diabetic retinopathy, which may help improve compliance with the standards of eye care for patients with T2DM .
Several limitations of this study deserve comment. First, the subjects who participated in the study were recruited from Shanghai, so they may not have been representative of China as a whole. Second, we need a larger sample to further validate the effect of the polymorphism on DR risk among different ethnic populations. In addition, our study focused primary analyses on a comparison of patients with DR to patients without DR rather than comparing those with advanced or proliferative DR to those without advanced DR. The possibility of an overlap with other causes of retinopathy may have limited the proper identification of DR genes in mild stages of the disease.
In summary, we identified the AA genotype or A allele of KCNJ11 (rs5219) as a genetic risk factor for DR in T2DM patients. Further confirmatory studies on the functional effect of KCNJ11 variants should be conducted to reveal its real contribution to T2DM.
advanced glycation end products
blood urea nitrogen
0 h C-peptide
diastolic blood pressure
diabetic macular edema
fasting plasma glucose
Genome- wide association studies
minor allele frequency
multivariate logistic regression
non-mydriatic digital stereoscopic retinal imaging
oral glucose tolerance test
postprandial plasma glucose
receptor for advanced glycation end products
systolic blood pressure
single nucleotide polymorphism
Type 2 diabetes mellitus
ubunit sulfonylurea receptor 1
vascular endothelial growth factor
Qian Y, Lu F, Dong M, Lin Y, Li H, Chen J, et al. Genetic variants of IDE-KIF11-HHEX at 10q23.33 associated with type 2 diabetes risk: a fine-mapping study in Chinese population. PLoS One. 2012;7:e35060. doi:10.1371/journal.pone.0035060.
Uno S, Zembutsu H, Hirasawa A, Takahashi A, Kubo M, Akahane T, et al. A genome-wide association study identifies genetic variants in the CDKN2BAS locus associated with endometriosis in Japanese. Nat Genet. 2010;42:707–10.
Shaw JE, Sicree RA, Zimmet PZ. Global estimates of the prevalence of diabetes for 2010 and 2030. Diabetes Res Clin Pract. 2010;87:4–14.
Cheung N, Mitchell P, Wong TY. Diabetic retinopathy. Lancet. 2010;376:124–36.
Resnikoff S, Pascolini D, Etya’ale D, Kocur I, Pararajasegaram R, Pokharel GP, et al. Global data on visual impairment in the year 2002. Bull World Health Organ. 2004;82:844–51.
Paine SK, Basu A, Mondal LK, Sen A, Choudhuri S, Chowdhury IH, et al. Association of vascular endothelial growth factor, transforming growth factor beta, and interferon gamma gene polymorphisms with proliferative diabetic retinopathy in patients with type 2 diabetes. Mol Vis. 2012;18:2749–57.
Hoerger TJ, Segel JE, Gregg EW, Saaddine JB. Is glycemic control improving in U.S. adults? Diabetes Care. 2008;31:81–6.
Brownlee M. Biochemistry and molecular cell biology of diabetic complications. Nature. 2001;414:813–20.
Hallman DM, Huber Jr JC, Gonzalez VH, Klein BE, Klein R, Hanis CL. Familial aggregation of severity of diabetic retinopathy in Mexican Americans from Starr County, Texas. Diabetes Care. 2005;28:1163–8.
Fu YP, Hallman DM, Gonzalez VH, Klein BE, Klein R, Hayes MG, Cox NJ, Bell GI, Hanis CL: Identification of diabetic retinopathy genes through a genome-wide association study among Mexican-Americans from Starr County, Texas. J Ophthalmol, 2010. doi:10.1155/2010/861291
Grassi MA, Tikhomirov A, Ramalingam S, Below JE, Cox NJ, Nicolae DL. Genome-wide meta-analysis for severe diabetic retinopathy. Hum Mol Genet. 2011;20:2472–81.
Huang YC, Lin JM, Lin HJ, Chen CC, Chen SY, Tsai CH, et al. Genome-wide association study of diabetic retinopathy in a Taiwanese population. Ophthalmology. 2011;118:642–8.
Sheu WH, Kuo JZ, Lee IT, Hung YJ, Lee WJ, Tsai HY, et al. Genome-wide association study in a Chinese population with diabetic retinopathy. Hum Mol Genet. 2013;22:3165–73.
Li Q, Chen M, Zhang R, Wang J, Zhou J, Bao Y, et al. KCNJ11 E23K variant is associated with therapeutic effect of sulfonylureas in Chinese type 2 diabetic patients. Clin Exp Pharmacol Physiol. 2014;41:748–54.
Bell GI, Polonsky KS. Diabetes mellitus and genetically programmed defects in beta-cell function. Nature. 2001;414:788–91.
Gloyn AL, Weedon MN, Owen KR, Turner MJ, Knight BA, Hitman G, et al. 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. 2003;52:568–72.
Javorsky M, Klimcakova L, Schroner Z, Zidzik J, Babjakova E, Fabianova M, et al. KCNJ11 gene E23K variant and therapeutic response to sulfonylureas. Eur J Intern Med. 2012;23:245–9.
Alberti KG, Zimmet PZ. 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. 1998;15:539–53.
Chun MY, Hwang HS, Cho HY, Chun HJ, Woo JT, Lee KW, et al. Association of vascular endothelial growth factor polymorphisms with nonproliferative and proliferative diabetic retinopathy. J Clin Endocrinol Metab. 2010;95:3547–51.
Galetovic D, Karlica D, Bucan K, Znaor L. The role of clinical and metabolic factors in the genesis and development of diabetic retinopathy. Acta Med Croatica. 2008;62:23–7.
Zhong ZL, Han M, Chen S. Risk factors associated with retinal neovascularization of diabetic retinopathy in type 2 diabetes mellitus. Int J Ophthalmol. 2011;4:182–5.
Esteves J, Laranjeira AF, Roggia MF, Dalpizol M, Scocco C, Kramer CK, et al. Diabetic retinopathy risk factors. Arq Bras Endocrinol Metabol. 2008;52:431–41.
Shah CA. Diabetic retinopathy: A comprehensive review. Indian J Med Sci. 2008;62:500–19.
Balasubbu S, Sundaresan P, Rajendran A, Ramasamy K, Govindarajan G, Perumalsamy N, et al. Association analysis of nine candidate gene polymorphisms in Indian patients with type 2 diabetic retinopathy. BMC Med Genet. 2010;11:158.
Fong DS, Aiello LP, Ferris 3rd FL, Klein R. Diabetic retinopathy. Diabetes Care. 2004;27:2540–53.
Zhou D, Zhang D, Liu Y, Zhao T, Chen Z, Liu Z, et al. The E23K variation in the KCNJ11 gene is associated with type 2 diabetes in Chinese and East Asian population. J Hum Genet. 2009;54:433–5.
Koo BK, Cho YM, Park BL, Cheong HS, Shin HD, Jang HC, et al. Polymorphisms of KCNJ11 (Kir6.2 gene) are associated with Type 2 diabetes and hypertension in the Korean population. Diabet Med. 2007;24:178–86.
Nielsen EM, Hansen L, Carstensen B, Echwald SM, Drivsholm T, Glumer C, et al. The E23K variant of Kir6.2 associates with impaired post-OGTT serum insulin response and increased risk of type 2 diabetes. Diabetes. 2003;52:573–7.
Haider S, Antcliff JF, Proks P, Sansom MS, Ashcroft FM. Focus on Kir6.2: a key component of the ATP-sensitive potassium channel. J Mol Cell Cardiol. 2005;38:927–36.
Kane GC, Behfar A, Dyer RB, O'Cochlain DF, Liu XK, Hodgson DM, et al. KCNJ11 gene knockout of the Kir6.2 KATP channel causes maladaptive remodeling and heart failure in hypertension. Hum Mol Genet. 2006;15:2285–97.
Seino S, Miki T. Physiological and pathophysiological roles of ATP-sensitive K+ channels. Prog Biophys Mol Biol. 2003;81:133–76.
Li JY, Li ZB, Zhu M, Liu YQ, Li Y, Wang SW, et al. Mutational analysis of KCNJ11 in Chinese elderly essential hypertensive patients. J Geriatr Cardiol. 2012;9:153–7.
Tong Z, Yang Z, Patel S, Chen H, Gibbs D, Yang X, et al. Promoter polymorphism of the erythropoietin gene in severe diabetic eye and kidney complications. Proc Natl Acad Sci U S A. 2008;105:6998–7003.
We thank the grants from the National Natural Science Foundation of China (No. 81270903, 30900501, 81370884) and the Natural Science Foundation of Shanghai City (No.10ZR1405500) for support of the study.
The authors declare that they have no competing interests.
N-JL analyzed data and wrote the manuscript. H-HW, ZY, X-MT, Y-PD, X-CW, BL and Z-YZ contributed samples, reagents and analysis tools. JW and R-MH conceived and designed the study. All authors read and approved the final manuscript.
Nai-Jia Liu and Hui-Hui Wu equal contributors.
About this article
Cite this article
Liu, N., Wu, H., Li, Y. et al. An analysis of the association between a polymorphism of KCNJ11 and diabetic retinopathy in a Chinese Han population. Eur J Med Res 20, 3 (2015). https://doi.org/10.1186/s40001-014-0075-3
- Diabetic retinopathy
- Chinese Han population