- Open Access
Association of KCNB1 polymorphisms with lipid metabolisms and insulin resistance: a case-control design of population-based cross-sectional study in Chinese Han population
Lipids in Health and Diseasevolume 14, Article number: 112 (2015)
The Erratum to this article has been published in Lipids in Health and Disease 2015 14:144
In our previous study, we had assessed in the Chinese Han population the association of KCNB1 rs1051295 with metabolic traits indicating metabolic syndrome, and showed that KCNB1 rs1051295 genotype TT was associated with increase of waist to hip ratio (WHR), fasting insulin (FINS), triglycerides (TG) and decreased insulin sensitivity at basal condition. Here, we aimed at detecting whether there were associations between other tag SNPs of KCNB1 and favorable or unfavorable metabolic traits.
We conducted a case–control design of population-based cross-sectional study to investigate the association between each of the 22 candidates tag SNPs of KCNB1 and metabolic traits in a population of 733 Chinese Han individuals. The association was assessed by multiple linear regression analysis or unconditional logistic regression analysis.
We found that among the 22 selected tag SNPs, four were associated with an increase (rs3331, rs16994565) or decrease (rs237460, rs802950) in serum cholesterol levels; two of these (rs237460, rs802590) further associated or were associated with reduced serum LDL-cholesterol. Two novel tag SNPs (rs926672, rs1051295) were associated with increased serum TG levels. We also showed that KCNB1 rs926672 associated with insulin resistance by a case–control study, and two tag SNPs (rs2057077and rs4810952) of KCNB1 were associated with the measure of insulin resistance (HOMA-IR) in a cross-section study.
These results indicate that KCNB1 is likely associated with metabolic traits that may either predispose or protect from progression to metabolic syndrome. This study provides initial evidence that the gene variants of KCNB1, encoding Kv2.1 channel, is associated with perturbation of lipid metabolism and insulin resistance in Chinese Han population.
KCNB1, also known as Kv2.1, is located on chromosome 20q13.2, and is mainly expressed in the brain cortex and hippocampus, islet β-cells, cardiac atrium and ventricle, skeletal muscles and some other tissues . Our recently reported case–control study in Chinese Han population showed a strong association of KCNB1 rs1051295 genotype TT with an increased risk of becoming afflicted with metabolic syndrome leading to T2D. Specifically, we showed the genotype TT of this SNP was associated with several metabolic traits of metabolic syndrome, including increased waist to hip ratio (WHR), fasting insulin (FINS) levels, increased serum triglyceride (TG) levels, and decreased insulin sensitivity at basal condition .
Kv2.1 channel’s role in secretion is attributed not only to the regulation of membrane excitability but also to a more direct role in the exoctyotic machinery. Up-regulation of Kv2.1 channel in PC12 , bovine chromaffin cells , rat dorsal root ganglion cells  and rodent and human islet β-cells,  facilitated exocytotic SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) complex assembly initiated by direct interaction of its C-terminus with SNARE protein syntaxin-1A and subsequent further assembly with VAMP2 and SNAP25. In β-cells, these Kv2.1-SNARE complex interactions can enhance insulin secretion . To our knowledge, there have been no in vivo or in vitro studies to demonstrate a possible role of Kv2.1 channel in metabolism to account for the perturbed metabolic traits we had reported . However, in support of this possibility that Kv channels could be involved in body metabolic functions, Chandy and colleagues recently demonstrated that ShK-186, a selective and potent blocker of voltage-gated Kv1.3 channel, reduced weight gain, adiposity, and fatty liver with resultant enhancement of peripheral insulin sensitivity; and also decreased serum levels of cholesterol, glucose, HbA1c, insulin, and leptin .
To delve deeper into the role of KCNB1-encoded Kv2.1 in metabolic syndrome, we examined the associations between the tag SNPs of KCNB1 and different metabolic traits of metabolic syndrome. This study explores the possible role of Kv2.1 and its influence in the body’s metabolism leading to T2D.
General information of the participants
Of the 733 participants, mean age was 39.9 ± 10.4 years old, with males constituting 43 % (315/733) of the total participants. FPG, TC, TG, HDL-C and LDL-C for the participants was 4.9 ± 0.5, 4.6 ± 0.9, 1.1 ± 0.6, 1.7 ± 0.3 and 2.7 ± 0.7 mmol/L, respectively, and the HOMA-IR and HOMA-B% was 1.2 ± 0.4and 108.7 ± 31.9 %, respectively.
SNPs of KCNB1 associated with changes in serum TC or LDL-C
rs237460: Compared with genotype AA, genotype GG (P = 0.02, b = -0.11 (-0.19,-0.02)) and GA (P = 0.02, b = -0.19(-0.34, -0.04)) was associated with decreased serum TC levels. We also found that serum levels of LDL-C of genotype GA (2.60 mmol/L) (P = 0.005, b = -0.18(-0.31, -0.06)) and GG (2.64 mmol/L) (P = 0.07, b = -0.07 (-0.14, 0.004)) was lower or trend to be lower than that in genotype AA (2.78 mmol/L) (Table 1).
rs802590: Genotype AC was associated with decreased serum levels of TC (P = 0.03, b = -0.36(-0.68, -0.05)) and decreased serum levels of LDL-C (P = 0.06, b = -0.12 (-0.24, 0.01)) compared with genotype CC, whereas genotype AA showed no association (for TC: P = 0.17, b = -0.11(-0.26, 0.05); for LDL-C:P = 0.61, b = -0.03(-0.16,0.09)) (Table 1).
rs16994565: Compared with genotype CC, the serum TC levels in genotype TT and TC increased by 22.7 % (P = 0.02, b = 0.36(-0.06, 0.66)) and 18.9 % (P = 0.08, b = 0.58 (-0.06, 1.23)), respectively (Table 1).
rs3331: genotype TT (P = 0.01, b = 0.49 (0.13, 0.84) and TC (P = 0.03, b = 0.83 (0.09, 1.57)) was associated with increased serum levels of TC compared to genotype CC. Although the TG in genotype TT (1.11 mmol/L) increased 15.4 % than genotype CC (0.94 mmol/L), it did not reach a significant level (P = 0.79, b = 0.03(-0.21, 0.27)) (Table 1).
SNPs of KCNB1 associated with changes in serum TG
rs926672: Compared with genotype TT, the serum TG levels in genotype AA (P = 0.004, b = 0.08(0.02, 0.13)) and AT (P = 0.002, b = 0.19(0.07, 0.30)) increased 14.5 and 18.3 %, respectively. In addition, genotype AT was associated with increased TC levels (P = 0.049, b = 0.16(0.00, 0.32)), insulin resistance (P = 0.03, b = 0.74(0.06, 1.42)), and fasting insulin (FINS) (P = 0.02, b = 0.10 (0.01, 0.19)) compared with genotype TT, and showed a trend in increasing islet β-cell function by HOMA-B% (P = 0.09, b = 5.66(-0.86, 12.18)) (Table 1).
rs1051295: Compared with genotype CC, both of genotype TT (P = 0.03, b = 0.06 (0.01, 0.12)) and TC (P = 0.10, b = 0.10(-002, 0.22)) was associated or inclined to be associated with increased serum TG levels, in line with our previous report . We also found that genotype TT tended to be associated with increased BMI (P = 0.09, b = -0.22 (-0.48, 0.04)), HOMA-IR (P = 0.09, b = 0.04(-0.01, 0.09)), FINS (P = 0.11, b = 2.92(-0.06, 0.65)) and HOMA-B (P = 0.09, b = 5.66(-0.86, 12.18)) compared with genotype CC (Table 1).
SNPs of KCNB1 associated with the measure of insulin resistance by HOMA-IR
rs4810952: Compared with genotype TT, genotype TC was associated with increased HOMA-IR (P = 0.01, b = -9.84(-16.65,-3.02)), FINS (P = 0.003, b = 0.86(0.30, 1.43)) and the HOMA-B% (P < 0.001, b = 9.60(4.27, 14.93)) (Table 1).
rs2057077: Compared with genotype TT (1.06 mmol/L and 8.29mUI/L), HOMA-IR and FINS of TC (1.19 mmol/L and 9.26mUI/L) increased 12.3 % (P = 0.008, b = 0.13(0.03, 0.23)) and 11.70 % (P = 0.01, b = 0.92 (0.19, 1.65)), respectively (Table 1).
SNPs of KCNB1 associated with insulin resistance (IR)
According to the definition of insulin resistance, we identified 199 IR subjects and 534 Non-IR subjects from 733 participants, and the prevalence of IR were 27.15 % (199/733) in this study. The clinical anthropometric and biochemical characteristics of the participants are summarized in Table 2. Among 22 available tag SNPs of KCNB1, the rs926672 (P = 0.03) and rs2057077 (P = 0.04) were associated with insulin resistance by univariate analysis. However, only rs926672 was associated with IR (P = 0.03) after being adjusted for age and gender (Table 3). We found that the risk for the carrier of genotype AT suffering from insulin resistance was 1.77 fold to that of genotype TT (P = 0.03, OR = 1.77, (1.07-2.93)). In the dominant or recessive model, neither allele T (genotype TT and AT) (P = 0.50, OR = 1.13(0.79, 1.62)) nor allele A (genotype AA and AT) (P = 0.057, OR = 1.60(0.99-2.60)) was associated with IR. Although we had already found that KCNB1 rs1051295 genotype TT was associated with T2D and decreased insulin sensitivity in our previous study, no association was observed for rs1051295 and insulin resistance in either dominant model (P = 0.72, OR = 1.07 (0.75,1.51)) or recessive model (P = 0.24, OR = 0.75(0.47,1.20)) by univariate or multivariate analysis in this study.
Haplotypes Analysis of KCNB1 with Insulin Resistance
We constructed haplotypes for 22 tagging SNPs to examine whether these SNPs demonstrated any additional evidence of association with IR when tested together by using the Haploview program. We identified 4 haplotypes for the KCNB1 gene and each haplotype represented sites are as follows: Block1 (rs4810952-rs1961192),Block2 (rs6125642-rs16994546-rs2057077), Block3(rs6019852-rs6095516-rs742759), and Block4 (rs802950-rs6019857-rs653070-rs552068) (shown in Fig. 1). We found that the haplotype of Block 3 GCA was significantly different between the IR group and the control group (P = 0.04). Details are shown in Table 4. Rs926672 was not in a haplotype block generated in this study due to weak linkage disequilibrium with other SNPs.
Here, we first report four genetically non-linked KCNB1 tag SNPs (rs237460, rs3331, rs16994565 and rs802950, all r2 < 0.65) to be associated with significant changes in serum TC levels, with two of these (rs237460 and rs802590) further associated or trend to be associated with reduced LDL-C levels. In addition, we also demonstrated that one novel tag SNPs rs926672 was associated with increased serum TG levels, in addition to the previously identified KCNB1 rs1051295 . We also showed KCNB1 rs926672 associated with insulin resistance by a case–control study, and two tag SNPs (rs2057077and rs4810952) of KCNB1 were associated with the measure of insulin resistance as determined by HOMA-IR in a cross-section study. Together, these results indicate that KCNB1 is likely associated with changes in lipid metabolism and insulin sensitivity.
Although the role of Kv2.1 channel in lipid and fat metabolism has not been assessed, previous reports could provide some clues. Besides the simple diffusion manner, fatty acids, previously thought to be transported across the cell membrane by simple diffusion, is now known to also do so by competitive receptor-mediated transporters including fatty acid translocase (FAT), fatty acid binding protein (FABP) and fatty acid transport protein (FATP) [8–10]. Fatty acid transport protein-1 (FATP1) which transports free fatty acids (FFA) is insulin-sensitive and FATP1 expression in adipocytes is influenced by insulin-mediated transcriptional modulation . Recent studies have shown that the interaction between the C-terminus of Kv2.1 channel and SNARE protein syntaxin-1A in islet beta cell is required for efficient insulin exocytosis and glucose-stimulated insulin secretion  . In addition, studies have demonstrated that FAT levels would be improved with adipocyte differentiation [12, 13], while Kv channels Kv2.1 and Kv3.3 were shown to possibly play a role in the differentiation of pre-adipocytes and human mesenchymal stem cells into adipocytes [14, 15]. These postulated mechanisms may underlie the association between SNPs of Kv2.1 channel and lipid metabolism, and a possible reason for these Kv2.1 SNPs being associated with insulin resistance, whereby dysregulated lipid metabolism is one of the risk factors for insulin resistance .
In addition, there was substantial evidence that internal interactions within some loci in gene might have a large impact on the phenotype we observed. Many studies have shown that neither single SNP loci nor multiple SNP loci have brought us greater effectiveness than haplotype-based research in the analysis of complex diseases . Haplotype was likely to have more information about linkage disequilibrium, which made it easier to find and locate causative loci in association analysis [18, 19]. In this study we constructed four haplotypes for the KCNB1 gene and found Block 3 (rs6019852-rs6095516-rs742759) genotype GCA was significantly different distributed between the IR group and the control group (P = 0.04). This result indicated that the haplotype composed of three SNPs above is likely to have association with IR by internal interactions. In addition,we also found that KCNB1 rs926672 was not in any haplotype block generated in this study due to weak linkage disequilibrium with other SNPs, which suggests that rs926672 may be an independent genetic marker of insulin resistance.
KCNB1 SNPs have been showed to be associated with some diseases in human, such as rheumatoid arthritis, left ventricular hypertrophy, long QT syndrome and type 2 diabetes in previous studies [2, 20–22]. In this study, we demonstrated that eight of the 22 selected tag SNPs are associated with changes in serum lipids or insulin resistance, which adds new insight to our previous report showing the association of KCNB1 rs1051295 with type 2 diabetes . Since these SNP are not located in the coding regions of KCNB1, 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 potential mechanisms that these SNPs exert their effects on lipid metabolism and insulin sensitivity.
The limitation of this study is the relative small sample size, thus, the association needs to be verified in a study with larger sample size within the Chinese Han population, and in other non-Chinese populations. Also, because the small sample size, we did not do a multiple testing correction for the significance threshold, although it was conservative for an explore study . In view of lipid metabolism and insulin resistance are central to the etiology of many complex diseases, future research needs to focus on whether these KCNB1 SNPs are associated with these diseases as well.
This study provides the initial evidence for the association between KCNB1 polymorphisms with lipid metabolism and insulin resistance in Chinese Han population. This initial study indicates that Kv2.1 may play a role in the lipid metabolism and thus influence insulin sensitivity. Further studies are needed to examine the precise mechanism underlying the association of KCNB1 polymorphisms with lipid metabolism and insulin resistance.
Materials and methods
This study was approved by the Ethics Committee of Capital Medical University (Beijing, China) and was conducted in accordance with the principles of the Helsinki Declaration II. Written consent was obtained from all participants.
Study participants and data collection
This cross-sectional study involved 733 Chinese Han participants included 315 males and 418 females without overt diabetes. Subjects were consecutively recruited at the Physical Examination Center, Beijing Xuanwu Hospital. Venous blood samples were drawn from all participants after an overnight fast of 8-10 h for DNA preparation and determination of fasting plasma glucose (FPG), fasting insulin (FINS), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C). Height, weight, waist circumference (WC) and hip circumference (HC) of participants were measured according to standard protocols. Body mass index (BMI) was calculated by dividing the weight (in kg) by the height (in m) squared. Waist to hip ratio (WHR) was calculated by dividing the waist circumference by the hip circumference.
Tag SNPs selection and genotyping
Genotyping data was obtained from the CHB (Han Chinese in Beijing) panel of the Build 36 of HapMap Project to choose the 22 tag SNPs within the range of the first and last exon (transcribed region) of KCNB1 on the basis of the following principal criteria: r2 > 0.8(100 % coverage) and minor allele frequency (MAF) > 0.05. Of those tag SNPs, 20 were located in introns, and 2 in the 3’-UTR (rs3331 and rs1051295). Characteristics of the selected tag SNPs, included genomic position, genic position, minor allele frequency, and Hardy-Weinberg equilibrium shown in Table 5.
DNA was isolated from peripheral white blood cells using DNeasy tissue kit (Qiagen) and samples were diluted to a final concentration of 20-50 ng/ml. Genotyping was assessed using Mass Array (Sequenom, San Diego, CA) based on allele-specific MALDI-TOF mass spectrometry. DNA from samples was randomly assigned to 96 well plates, with genotyping performed in blinded manner.
PCR was performed for a 384-well microtiter plate of reactions in which the same assay was applied to different DNAs. Each PCR procedure was carried out in a total volume of 5 μL. The PCR cocktail contained 1 μL of genomic DNA, 1 μL of each primer, 0.5 μL of 10*PCR buffer, 0.1 μL deoxy-ribonucleoside triphosphate, 0.4 μL MgCl2, 0.2 μL Hotstar and 1.8 μL H2O. In the thermal reactor, amplification of the KCNB1 was achieved under the following conditions: denaturation at 94 °C for 4 min followed by 45 amplification cycles of (denaturation at 94 °C for 20 s, annealing at 56 °C for 30 s, extension at 72 °C for 1 min) and a final extension step of 3 min at 72 °C.
The genotyping success rates were > 98.5%and all tag SNPs are consistent with the Hardy-Weinberg equilibrium (P > 0.05) except for rs653070 and rs742759 for which the P value was 0.04 and 0.01, respectively, thus excluded from further analysis. Genotyping of all samples was repeated in 5 % of random samples for verification and quality control; all revealed the genotype data had an error rate of <1 %. The genetic linkage between tag SNPs ranged from ‘none’ (D’ = 0.00, r2 = 0.00) to ‘moderate’ (D’ = 1.00, r2 = 0.79, shown in Fig. 1).
The value of homeostasis model assessment (HOMA) denoting insulin resistance (IR) was termed as HOMA- IR and was calculated as equal to (fasting insulin [in IU/mL]*fasting glucose [in mmol/L]/22.5) . Subjects having HOMA-IR values in the highest quartile (P75) were defined as insulin resistance in our study .
Statistical analysis was performed using SPSS software, version 18.0 (SPSS, Inc, Chicago, IL. purchased by Capital Medical University, China). The quantitative variables of metabolic traits were expressed as mean ± SD. The difference of these variables among the 3 genotypes of tag SNPs or between IR and Non-IR groups was assessed by AVONA or independent t test, respectively. The χ 2 test was used to examine the differences of SNPs or gender between the IR and Non-IR groups. The association between the tag SNPs and the metabolic traits or insulin resistance was finally determined by multiple linear regression analysis or unconditional logistic regression analysis, in which age and gender were adjusted as confounding factors, and b or OR and their 95 % confidence interval (CI) were presented.
The Haploview program (http://www.broad.mit.edu/mpg/haploview/) was used to select tag SNPs, calculate pairwise linkage disequilibrium statistics and examine allelic and haplotype associations with insulin resistance. The Haploview program was also used to draw the figure of pairwise LD pattern. Hardy-Weinberg equilibrium was determined for each SNP distribution. A value of P < 0.05 was considered as significant.
High-density lipoprotein cholesterol
Low-density lipoprotein cholesterol
Fasting plasma glucose
Waist to hip ratio
Body mass index
Minor allele frequency
Gutman GA, Chandy KG, Grissmer S, Lazdunski M, McKinnon D, Pardo LA, et al. International Union of Pharmacology. LIII. Nomenclature and molecular relationships of voltage-gated potassium channels. Pharmacol Rev. 2005;57:473–508.
Zhang Y-XX, Liu Y, Dong J, Wang Y-XX, Wang J, Zhuang G-QQ, et al. An exploratory study of the association between KCNB1 rs1051295 and type 2 diabetes and its related traits in Chinese Han population. PLoS ONE. 2013;8:e56365.
Singer-Lahat D, Sheinin A, Chikvashvili D, Tsuk S, Greitzer D, Friedrich R, et al. K+ channel facilitation of exocytosis by dynamic interaction with syntaxin. J Neurosci. 2007;27:1651–8.
Feinshreiber L, Singer-Lahat D, Ashery U, Lotan I. Voltage-gated potassium channel as a facilitator of exocytosis. Ann N Y Acad Sci. 2009;1152:87–92.
Feinshreiber L, Singer-Lahat D, Friedrich R, Matti U, Sheinin A, Yizhar O, et al. Non-conducting function of the Kv2.1 channel enables it to recruit vesicles for release in neuroendocrine and nerve cells. J Cell Sci. 2010;123(Pt 11):1940–7.
Dai XQ, Manning Fox JE, Chikvashvili D, Casimir M, Plummer G, Hajmrle C, et al. The voltage-dependent potassium channel subunit Kv2.1 regulates insulin secretion from rodent and human islets independently of its electrical function. Diabetologia. 2012;55:1709–20.
Upadhyay SK, Eckel-Mahan KL, Mirbolooki MR, Tjong I, Griffey SM, Schmunk G, et al. Selective Kv1.3 channel blocker as therapeutic for obesity and insulin resistance. Proc Natl Acad Sci U S A. 2013;110:E2239–48.
Lobo S, Wiczer BM, Smith AJ, Hall AM, Bernlohr DA. Fatty acid metabolism in adipocytes: functional analysis of fatty acid transport proteins 1 and 4. J Lipid Res. 2007;48:609–20.
Doege H, Stahl A. Protein-mediated fatty acid uptake: novel insights from in vivo models. Physiology (Bethesda). 2006;21:259–68.
Nickerson JG, Alkhateeb H, Benton CR, Lally J, Nickerson J, Han X-XX, et al. Greater transport efficiencies of the membrane fatty acid transporters FAT/CD36 and FATP4 compared with FABPpm and FATP1 and differential effects on fatty acid esterification and oxidation in rat skeletal muscle. J Biol Chem. 2009;284:16522–30.
Wu Q, Ortegon AM, Tsang B, Doege H, Feingold KR, Stahl A. FATP1 is an insulin-sensitive fatty acid transporter involved in diet-induced obesity. Mol Cell Biol. 2006;26:3455–67.
Hajri T, Abumrad NA. Fatty acid transport across membranes: relevance to nutrition and metabolic pathology. Annu Rev Nutr. 2002;22:383–415.
Laugerette F, Passilly-Degrace P, Patris B, Niot I, Febbraio M, Montmayeur J-PP, et al. CD36 involvement in orosensory detection of dietary lipids, spontaneous fat preference, and digestive secretions. J Clin Invest. 2005;115:3177–84.
You MH, Song MS, Lee SK, Ryu PD, Lee SY, Kim DY. Voltage-gated K+ channels in adipogenic differentiation of bone marrow-derived human mesenchymal stem cells. Acta Pharmacol Sin. 2013;34:129–36.
Ramírez-Ponce MP, Mateos JC, Bellido JA. Insulin increases the density of potassium channels in white adipocytes: possible role in adipogenesis. J Endocrinol. 2002;174:299–307.
Liu J, Jahn LA, Fowler DE, Barrett EJ, Cao W, Liu Z. Free fatty acids induce insulin resistance in both cardiac and skeletal muscle microvasculature in humans. J Clin Endocrinol Metab. 2011;96:438–46.
Chen X, Li Z. Inference of haplotype effects in case-control studies using unphased genotype and environmental data. Biom J. 2008;50:270–82.
Clark AG. The role of haplotypes in candidate gene studies. Genet Epidemiol. 2004;27:321–33.
Morris RW, Kaplan NL. On the advantage of haplotype analysis in the presence of multiple disease susceptibility alleles. Genet Epidemiol. 2002;23:221–33.
Xiao X, Zhang Y, Wang K. Association of KCNB1 to rheumatoid arthritis via interaction with HLA-DRB1. BMC Proc. 2009;7(3 Suppl):S134.
Arnett DK, Li N, Tang W, Rao DC, Devereux RB, Claas SA, et al. Genome-wide association study identifies single-nucleotide polymorphism in KCNB1 associated with left ventricular mass in humans: the HyperGEN Study. BMC Med Genet. 2009;10:43.
Iwasa H, Kurabayashi M, Nagai R, Nakamura Y, Tanaka T. Multiple single-nucleotide polymorphisms (SNPs) in the Japanese population in six candidate genes for long QT syndrome. J Hum Genet. 2001;46:158–62.
William Stafford N. How does multiple testing correction work. Nat Biotechnol. 2009;27((12):1135–7.
Resnick HE, Jones K, Ruotolo G, Jain AK, Henderson J, Lu W, et al. Insulin resistance, the metabolic syndrome, and risk of incident cardiovascular disease in nondiabetic american indians: the Strong Heart Study. Diabetes Care. 2003;26:861–7.
Hanson RL, Imperatore G, Bennett PH, Knowler WC. Components of the “metabolic syndrome” and incidence of type 2 diabetes. Diabetes. 2002;51:3120–7.
This work was supported by grants from the Key Projects in the National Science & Technology Pillar Program (No.2015BAI09B01) and the Natural Science Foundation of Beijing Municipality to Yan He (7132027).
The authors declare that they have no competing interests.
Conception of research idea and Secure funding: YCY, YH, ZWC; Study Design: YCY, YH, ZWC; Implementation (training; Data collection): YCY, JD, RYK, JD, YXZ, FL, YXY, RZ, LLX, XXP, LZ, DH, HYG; Study coordination and management: YCY, YH; Data management: YCY, YH; Data analysis and drafting paper: YCY, YH. All authors read and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/s12944-015-0144-9.