Genetic variants at 10q23.33 are associated with plasma lipid levels in a Chinese population

Sijun Liu, Yun Qin, Feng Lu, Meihu Dong, Yudi Lin, Huizhng Li, Chong Shen, Junheng Di, Yue Jing, Gungfu Jin, Zhibin Hu, Hongbing Shen,

aDepartment of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

bDepartment of Public Service Management, Kangda College, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

cDepartment of Chronic Non-communicable Disease Control, Wuxi Center for Disease Control and Prevention, Wuxi , Jiangsu 214023, China.

Genetic variants at 10q23.33 are associated with plasma lipid levels in a Chinese population

Sijun Liua,b,△, Yun Qiana,c,△, Feng Lua, Meihua Dongc, Yudi Linc, Huizhang Lia, Chong Shena, Juncheng Daia, Yue Jianga, Guangfu Jina, Zhibin Hua, Hongbing Shena,

aDepartment of Epidemiology and Biostatistics, School of Public Health, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

bDepartment of Public Service Management, Kangda College, Nanjing Medical University, Nanjing, Jiangsu 210029, China;

cDepartment of Chronic Non-communicable Disease Control, Wuxi Center for Disease Control and Prevention, Wuxi , Jiangsu 214023, China.

Plasma lipid abnormalities are implicated in the pathogenic process of type 2 diabetes. The IDE-KIF11-HHEX gene cluster on chromosome 10q23.33 has been identified as a susceptibility locus for type 2 diabetes. We hypothesized that genetic variants at 10q23.33 may be associated with plasma lipid concentrations. Seven tagging single nucleotide polymorphisms (SNPs: rs7923837, rs2488075, rs947591, rs11187146, rs5015480, rs4646957 and rs1111875) at 10q23.33 were genotyped in 3,281 subjects from a Han Chinese population, using the Taq-Man OpenArray and Sequenom MassARRAY platforms. Multiple linear regression analyses showed that SNP rs7923837 in the 3'-flanking region of HHEX was significantly associated with triglyceride levels (P = 0.019, 0.031 mmol/L average decrease per minor G allele) and that rs2488075 and rs947591 in the downstream region of HHEX were significantly associated with total cholesterol levels (P = 0.041, 0.058 mmol/L average decrease per minor C allele and P = 0.018, 0.063 mmol/L average decrease per minor A allele, respectively). However, the other four SNPs (rs11187146, rs5015480, rs4646957 and rs1111875) were not significantly associated with any plasma lipid concentrations in this Chinese population. Our data suggest that genetic variants in the IDE-KIF11-HHEX gene cluster at 10q23.33 may partially explain the variation of plasma lipid levels in the Han Chinese population. Further studies are required to confirm these findings in other populations.

cholesterol, triglycerides, polymorphism, genetic, iDE-KIF11-HHEX

INTRODUCTION

Plasma lipid abnormalities are associated with risk of type 2 diabetes[1,2]. Some studies have shown that elevated triglyceride (TG) levels and low levels of high-density lipoprotein cholesterol (HDL-C) accelerate the pathogenesis of type 2 diabetes[3-5]. Recently, several studies have investigated the potential effectof genetic factors associated with risk of type 2 diabetes on plasma lipid levels. Onuma et al.[6]analyzed the association of polymorphisms in the glucokinase (hexokinase 4) regulator (GCKR) gene with type 2 diabetes in a case-control study and with fasting blood glucose and TG levels in the general population. The A allele of SNP rs780094 was found to be associated with reduced risk of type 2 diabetes and lower levels of fasting plasma glucose, but higher levels of TG, in a Japanese population. In the study reported by Chen et al.[7], the authors found that subjects with minor alleles of SNPs rs2283228 and rs2237892 in the KQT-like subfamily member 1 (KCNQ1) gene, which were associated with type 2 diabetes, had higher levels of TG. This evidence suggests that genetic variants associated with diabetes risk may also be potential genetic determinants of plasma lipid levels.

A genome-wide association study (GWAS) conducted in a French case-control study identified a novel type 2 diabetes susceptibility locus on chromosome 10q23.33, which is located in a gene cluster including an insulin-degrading enzyme (IDE), a kinesininteracting factor 11 (KIF11), and a hematopoietically expressed homeobox protein (HHEX)[8]. Following this discovery, several studies have confirmed this association in British[9], Finnish[10], Japanese[11,12]and Chinese populations[13-15]. Recently, we have also found that SNPs rs7923837 and rs1111875 in the IDE-KIF11-HHEX locus at 10q23.33 were independently associated with risk of type 2 diabetes in a Chinese population[16]. However, the relationship between IDE-KIF11-HHEX locus and lipid traits in different populations is not clear. Therefore, in an effort to evaluate the influence of the polymorphisms in IDE-KIF11-HHEX locus on plasma lipid concentrations, we performed a fine-mapping study by genotyping seven tagging SNPs at 10q23.33 in 3,281 Han Chinese subjects to examine the associations of these variants with plasma levels of total cholesterol (TC), TG, HDL-C and low density lipoprotein cholesterol (LDL-C).

SUBJECTS AND METHODS

Study subjects

The subjects in the current study were selected from a community-based non-communicable diseases screening program comprised of more than 50,000 participants in Jiangsu Province during 2004 and 2008. All subjects were unrelated, ethnic Han Chinese. Subjects were excluded from the study if they had a history of diabetes, hypertension, coronary heart disease or cancer, or fasting plasma glucose≥5.6 mmol/L. After providing informed consent, all subjects were interviewed face-to-face using a standard questionnaire that included demographic characteristics, risk factors and disease history. Subjects who smoked 1 cigarette per day for over 1 year were defined as smokers, and those who consumed 3 or more alcohol drinks a week for over 6 months were considered as alcohol drinkers. Physical examinations, including measurements of height, weight and blood pressure, as well as laboratory tests to measure TC, TG, HDL-C and fasting plasma glucose concentrations, were performed for each participant. Sitting blood pressure was based on the average of three blood pressure readings measured. Body mass index (BMI) was calculated as weight (in kilograms) divided by the square of height (in meters). Fasting blood samples for routine laboratory examinations were obtained in the early morning after an overnight fast. All biochemical parameters were measured enzymatically on an auto-analyzer (Hitachi 7180 Biochemistry Auto-analyzer, Japan) according to the manufacturer's instructions. For subjects with TG levels ＜ 4.52 mmol/L, LDL-C levels were estimated indirectly using the Friedewald's formula. The Institutional Review Board of Nanjing Medical University approved the study.

SNP selection and genotyping

Based on our previous study[16], we used a blockbased tagging strategy to select tagging SNPs using Haploview 4.2 software according to the HapMap database［http://www.hapmap.org/, phaseII Nov08, on NCBI B36 assembly, dbSNPb126; population: Han Chinese population (CHB) and Japanese population (JPT)］. The criteria included SNPs with minor allele frequency (MAF)≥ 0.10, Hardy-Weinberg equilibrium P≥0.05 and call rate≥95% when using pairwise linkage disequilibrium (r2) of 0.8 as the threshold for each block. Seven tagging SNPs (rs7923837, rs2488075, rs947591, rs11187146, rs5015480, rs4646957 and rs1111875) associated with type 2 diabetes in the IDE-KIF11-HHEX locus at 10q23.33 were included in the current study.

Genomic DNA was isolated from leucocytes of venous blood by proteinase K digestion and phenol/ chloroform extraction. Genotyping was performed using the TaqMan OpenArray Genotyping System (Life Technologies, Carlsbad, CA, USA) and the iPLEX Sequenom MassARRAY platform (Sequenom, Inc.). For quality control, two non-template controls were used in each chip or plate. The overall call rates ranged from 98.8% to 99.8% for all SNPs.

Statistical analysis

Associations between the genotypes and plasma lipidconcentrations were determined by multiple linear regression analysis with adjustment for age, sex, smoking status, drinking status and BMI. The Hardy-Weinberg equilibrium was tested by a goodness-of-fit χ2test to compare the observed genotype frequencies with the expected ones among the 3,281 subjects. All statistical analyses were performed using Statistical Analysis System software version 9.1.3 (SAS Institute, Cary, NC, USA). All tests were two-sided and the significance level was set at P ＜ 0.05.

Table 1 Characteristics of the study population

RESULTS

The demographic and biochemical characteristics of the 3,281 subjects included in this study are shown in Table 1. The mean age was 56.58 (± 9.88) years and the mean BMI value was 22.12 (± 2.63) kg/m2. The mean values of TC, HDL-C, LDL-C and TG were 4.40 (± 0.80) mmol/L, 1.62 (± 0.38) mmol/L, 2.29 (± 0.74) mmol/L and 1.09 (± 0.45) mmol/L, respectively. Among these subjects, 827 subjects (25.38%) were smokers and 635 subjects (19.51%) were drinkers.

The observed genotype frequencies for the seven SNPs were all consistent with the Hardy-Weinberg equilibrium among 3,281 subjects (P ＞ 0.05) (Table 2). We examined the association between each SNP and TC, TG, HDL-C or LDL-C levels, respectively, in an additive model using a linear regression model with adjustment for age, sex, smoking, drinking and BMI (Table 2). We found significant associations between rs7923837 and TG (P = 0.019), between rs2488075 and rs947591 and TC (P = 0.041 and 0.018, respectively). As shown in Table 3, the G allele of rs7923837 was associated with a lower TG levels (0.031 mmol/L average decrease per G allele). Similarly, the C allele of rs2488075 and the A allele of rs947591 were both associated with lower TC levels (0.058 mmol/L average decrease per C allele and 0.063 mmol/L average decrease per A allele, respectively). Conditional analysis indicated that rs2488075 and rs947591 were not significant after adjustment with each other, as the two SNPs were in strong linkage equilibrium (LD) (r2=0.734). However, no significant associations were observed between the other four SNPs (rs11187146, rs5015480, rs4646957 and rs1111875) and blood lipid concentrations.

We then conducted a stratification analysis for rs7923837, rs2488075 and rs947591 by age, sex, BMI, smoking, and drinking status. As shown in Table 4, the associations between rs7923837 and TG levels were more evident among subjects of the low age group (P = 0.009), male subjects (P = 0.003), non-drinkers (P = 0.019), and subjects with low BMI (P = 0.003). The associations between rs2488075, rs947591 and TC levels were more evident among subjects of low age group (P = 0.021 and 0.013, respectively), female subjects (P = 0.048 and 0.032, respectively), non-smokers (P = 0.004 and 0.002, respectively), non-drinkers (P = 0.007 and 0.003, respectively) and subjects with low BMI (P = 0.038 and 0.037, respectively).

DISCUSSION

To the best of our knowledge, this is the first study to investigate the association between IDE-KIF11-HHEX polymorphisms and plasma lipid concentrations in a Chinese population. Of the seven tagging SNPs at the IDE-KIF11-HHEX locus, we found that rs7923837 was associated with plasma concentrations of TG, and rs2488075 and rs947591 were associated with plasma concentrations of TC.

Table 2 Association between lipid concentrations and selected single nucleotide polymorphisms

rs7923837 is located in the 3'-flanking region of the HHEX gene, which encodes a transcription factor that is involved in Wnt signaling and is critical for hepatic and pancreatic development[17,18]. In addition, HHEX may regulate β-cell development and/or function by activating hepatocyte nuclear factor 1α[19]. Several studies reported that the association between rs7923837 and type 2 diabetes is mediated through decreased β-cell secretory capacity or decreased β-cell mass[20-22]. Thus, HHEX is critical for insulin signaling and islet function[23]. HHEX may influence metabolic phenotypes such as TG and TC because insulin is necessary for the regulation of metabolic phenotypes. On chromosome 10q23.33, genetic variants in HHEX have been established as susceptibility loci for type 2 diabetes[8]. In our previous fine-mapping study[16], we reported that rs7923837 and rs1111875 were independently associated with risk of type 2 diabetes in a Chinese population. Several studies have also investigated the relationship between HHEX polymorphisms and other metabolic diseases. Zhao et al.[24]found that the type 2 diabetes risk-associated G allele of rs7923837 was associated with higher pediatric BMI in European American children. Cruz et al.[25]analyzed the association between the HHEX rs5015480 and risk of metabolic syndrome (MS) in a case-control study from Mexico city and found that rs5015480 was significantly associated with MS. Taken together, this evidence suggests that genetic variants in the IDEKIF11-HHEX gene cluster at 10q23.33 may contribute to metabolism-related traits and diseases, including circulating lipid levels and diabetes risk.

Genetic variants associated with the risk of type 2 diabetes have been found to influence plasma TG levels. The variant rs780094 in GCKR was associated with a decreased risk of type 2 diabetes, but with higher TG levels, in a GWAS of European population[26]. Similarly, in the current study, we found that the risk allele of diabetes was associated with lower TG levels. Considering that the subjects included in this study were from a healthy population, the lipid level-related variant might interpret the variation of lipid levels of baseline. The relationships between the same variant and diabetes risk and lipid levels impliesthat low lipid baseline levels in some subjects are, in part, genetically determined, causing them to be more susceptible to type 2 diabetes. This was also supported by the results from the stratification analysis, which showed that the associations between genetic variants at 10q23.33 and lipid levels were more evident among young subjects, non-smokers, non-drinkers and subjects with low BMI. However, the underlying mechanism remains unclear, and further studies are needed to elucidate the roles of genetic variants at 10q23.33 in circulating lipid levels.

Table 3 Effects of rs7923837, rs2488075 and rs947591 on plasma lipid concentrations

Table 4 Stratification analysis for rs7923837, rs2488075, and rs947591 and lipid levels in an additive genetic model

There is a strong relationship between glucose, cholesterol metabolism and type 2 diabetes[27]. The study by Hao et al.[28]suggested that plasma cholesterol plays a direct role in pancreatic islet dysfunction and may be a key factor underlying the progression of type 2 diabetes. Genetic variants that are associated with the risk of type 2 diabetes have also been found to influence plasma TC levels. Sanghera et al.[29]revealed a significant association between rs10885409 in TCF7L2 with type 2 diabetes and TC levels in Asian Indians. Chen et al.[30]found rs2237895 in KCNQ1, which was thought to be a candidate gene of diabetes that influenced plasma TC levels in the Han Chinese population. They argued that this variant might result in an increased expression of KCNQ1 and a subsequent increase in insulin secretion, which could stimulate lipid synthesis. In this study, we also found that rs2488075 and rs947591 were associated with plasma TC levels. rs2488075 and rs947591 are located downstream of HHEX, which may affect metabolic phenotypes. However, the mechanism that allows the CC genotype in rs2488075 and the AA genotype in rs947591 to contribute to lower levels of TC is still unknown.

Our study may be subject to certain limitations. First, the number of subjects in our study was moderate, thus the statistical power was limited. Second, the associations were not strong statistically; none of them passed multiple correction. Third, bias and reverse causation cannot be completely excluded in this observational epidemiological study. A Mendelian randomization approach that uses the random inheritance of genetic variants from parents to offspring may be of benefit in further studies. Thus, larger, well-designed epidemiological studies with ethnically diverse populations are warranted to confirm our findings.

In summary, the results in the current study indicate that genetic variants in the IDE-KIF11-HHEX gene cluster at 10q23.33 are associated with plasma lipid levels in the Chinese population. These findings highlight the important correlation between lipid levels and diabetes development at the genetic level. Further studies are needed to replicate our findings in other populations.

[1] Chaudhary R, Likidlilid A, Peerapatdit T, Tresukosol D, Srisuma S, Ratanamaneechat S, et al. Apolipoprotein E gene polymorphism: effects on plasma lipids and risk of type 2 diabetes and coronary artery disease. Cardiovasc Diabetol 2012; 11: 36.

[2] Cannon CP. Mixed dyslipidemia, metabolic syndrome, diabetes mellitus, and cardiovascular disease: clinical implications. Am J Cardiol 2008; 102: 5L-9L.

[3] Bitzur R, Cohen H, Kamari Y, Shaish A, Harats D.Triglycerides and HDL cholesterol: stars or second leads in diabetes? Diabetes Care 2009; 32(S2): 373-7.

[4] Mooradian AD. Dyslipidemia in type 2 diabetes mellitus. Nat Clin Pract Endocrinol Metab 2009; 5: 150-19.

[5] Wilson PW, Meigs JB, Sullivan L, Fox CS, Nathan DM, D'Agostino RB, Sr. Prediction of incident diabetes mellitus in middle-aged adults: the Framingham Offspring Study. Arch Intern Med 2007; 167: 1068-74.

[6] Onuma H, Tabara Y, Kawamoto R, Shimizu I, Kawamura R, Takata Y, et al. The GCKR rs780094 polymorphism is associated with susceptibility of type 2 diabetes, reduced fasting plasma glucose levels, increased triglycerides levels and lower HOMA-IR in Japanese population. J Hum Genet 2010; 55: 600-64.

[7] Chen Z, Yin Q, Ma G, Qian Q. KCNQ1 gene polymorphisms are associated with lipid parameters in a Chinese Han population. Cardiovasc Diabetol 2010; 9: 35.

[8] Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, et al. A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature 2007; 445: 881-5.

[9] Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, et al. Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science 2007; 316: 1336-41.

[10] Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, et al. A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science 2007; 316: 1341-5.

[11] Tabara Y, Osawa H, Kawamoto R, Onuma H, Shimizu I, Miki T, et al. Replication study of candidate genes associated with type 2 diabetes based on genome-wide screening. Diabetes 2009; 58: 493-8.

[12] Horikawa Y, Miyake K, Yasuda K, Enya M, Hirota Y, Yamagata K, et al. Replication of genome-wide association studies of type 2 diabetes susceptibility in Japan. J Clin Endocrinol Metab 2008; 93: 3136-41.

[13] Lin Y, Li P, Cai L, Zhang B, Tang X, Zhang X et al. Association study of genetic variants in eight genes/loci with type 2 diabetes in a Han Chinese population. BMC Med Genet 2010; 11: 97.

[14] Wu Y, Li H, Loos RJ, Yu Z, Ye X, Chen L, et al. 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 2008; 57: 2834-42.

[15] Zhou DZ, Liu Y, Zhang D, Liu SM, Yu L, Yang YF, et al. Variations in/nearby genes coding for JAZF1, TSPAN8/LGR5 and HHEX-IDE and risk of type 2 diabetes in Han Chinese. J Hum Genet 2010; 55: 810-5.

[16] 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.

[17] Bort R, Martinez-Barbera JP, Beddington RS, Zaret KS. Hex homeobox gene-dependent tissue positioning is required for organogenesis of the ventral pancreas. Development 2004; 131: 797-806.

[18] Hunter MP, Wilson CM, Jiang X, Cong R, Vasavada H, Kaestner KH, et al. The homeobox gene Hhex is essential for proper hepatoblast differentiation and bile duct morphogenesis. Dev Biol 2007; 308: 355-67.

[19] Tanaka H, Yamamoto T, Ban T, Satoh S, Tanaka T, Shimoda M, et al. Hex stimulates the hepatocyte nuclear factor 1alpha-mediated activation of transcription. Arch Biochem Biophys 2005; 442: 117-24.

[20] Pivovarova O, Nikiforova VJ, Pfeiffer AF, Rudovich N. The influence of genetic variations in HHEX gene on insulin metabolism in the German MESYBEPO cohort. Diabetes Metab Res Rev 2009; 25: 156-62.

[21] Pascoe L, Tura A, Patel SK, Ibrahim IM, Ferrannini E, Zeggini E, et al. Common variants of the novel type 2 diabetes genes CDKAL1 and HHEX/IDE are associated with decreased pancreatic beta-cell function. Diabetes 2007; 56: 3101-4.

[22] Cai Y, Yi J, Ma Y, Fu D. Meta-analysis of the effect of HHEX gene polymorphism on the risk of type 2 diabetes. Mutagenesis 2011; 26: 309-14.

[23] Xu P, Che Y, Cao Y, Wu X, Sun H, Liang F, et al. Polymorphisms of TCF7L2 and HHEX genes in Chinese women with polycystic ovary syndrome. J Assist Reprod Genet 2010; 27: 23-8.

[24] Zhao J, Bradfield JP, Zhang H, Annaiah K, Wang K, Kim CE, et al. Examination of all type 2 diabetes GWAS loci reveals HHEX-IDE as a locus influencing pediatric BMI. Diabetes 2010; 59: 751-5.

[25] Cruz M, Valladares-Salgado A, Garcia-Mena J, Ross K, Edwards M, Angeles-Martinez J, et al. Candidate gene association study conditioning on individual ancestry in patients with type 2 diabetes and metabolic syndrome from Mexico City. Diabetes Metab Res Rev 2010; 26: 261-70.

[26] Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, et al. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 2007; 316: 1331-6.

[27] Sanghera DK, Nath SK, Ortega L, Gambarelli M, Kim-Howard X, Singh JR, et al. TCF7L2 polymorphisms are associated with type 2 diabetes in Khatri Sikhs from North India: genetic variation affects lipid levels. Ann Hum Genet 2008; 72: 499-509.

[28] Hao M, Head WS, Gunawardana SC, Hasty AH, Piston DW. Direct effect of cholesterol on insulin secretion: a novel mechanism for pancreatic beta-cell dysfunction. Diabetes 2007; 56: 2328-38.

[29] Sanghera DK, Ortega L, Han S, Singh J, Ralhan SK, Wander GS, et al. Impact of nine common type 2 diabetes risk polymorphisms in Asian Indian Sikhs: PPARG2 (Pro12Ala), IGF2BP2, TCF7L2 and FTO variants confer a significant risk. BMC Med Genet 2008; 9: 59.

[30] Chen XD, Yang YJ, Li SY, Peng QQ, Zheng LJ, Jin L, et al. Several polymorphisms of KCNQ1 gene are associated with plasma lipid levels in general Chinese populations. PloS One 2012; 7: e34229.

Received 12 May 2012, Revised 13 June 2012, Accepted 18 November 2012, Epub 20 December 2012

This work was supported by grants from the Project of National Natural Science Foundation of China (No. 81102180, No. 81072379), Ministry of Health Research Program (No. WKJ2010-2-032), Wuxi Science & Technology Research Program (No. CSE01016) and the Priority Academic Program for the Development of Jiangsu Higher Education Institutions (Public Health and Preventive Medicine).

△These authors contribute equally to this work.

The authors reported no conflict of interests.

10.7555/JBR.27.20120091