- Open Access
Vitamin D receptor polymorphism rs2228570 is significantly associated with risk of dyslipidemia and serum LDL levels in Chinese Han population
Lipids in Health and Disease volume 17, Article number: 193 (2018)
The goal of this study was to determine if vitamin D receptor (VDR) gene polymorphisms underlie susceptibility to dyslipidemia in a Chinese Han population.
Three tag single nucleotide polymorphisms (SNPs) (rs11574129, rs2228570, and rs739837) were genotyped using TaqMan assays to determine VDR SNP associations with dyslipidemia. We genotyped 877 cases of dyslipidemia from a normotensive, non-diabetes mellitus population and 1822 non-dyslipidemia subjects in a stage I study. In a follow-up stage II study, we included a larger sample of 3124 controls and 1679 cases with dyslipidemia. Finally, we explored the potential molecular mechanism for the SNP associations using molecular modeling analysis.
We found a significant association between SNP rs2228570 and dyslipidemia in the additive (adjusted odds ratio (OR) = 1.255, 95% Confidence Interval (CI) = (1.118–1.409), P < 0.001), dominant (OR = 1.384, 95% CI = 1.384 (1.136–1.6), P = 0.001) and recessive models (OR = 1.356, 95%CI = 1.1–1.671, P = 0.004) in stage I. We further established that the rs2228570 variant was significantly associated with dyslipidemia in the additive (adjusted OR = 1.146, 95% CI = 1.053–1247, P = 0.002), dominant (OR = 1.184, 95%CI =1.018–1.376, P = 0.028) and recessive models (OR = 1.209, 95%CI = 1.064–1.374, P = 0.004) in stage II. The TT genotype was significantly higher (4.93 ± 0.75 mmol/L) compared to the TC (4.67 ± 0.47 mmol/L) or CC (4.66 ± 0.44 mmol/L) genotype (P = 0.01) in cases with elevated low-density lipoprotein cholesterol (LDL-C) levels. In contrast, the cases with the TT genotype had significantly lower serum 25(OH)D levels (18.43 ± 5.04 ng/ mL) compared to the TC (26.24 ± 4.16 ng/mL) and CC (36.76 ± 8.10 ng/ mL) genotypes (P < 0.001). Multivariable linear regression analysis indicated that the rs2228750 genotype significantly correlated with serum low-density lipoprotein-C (LDL-C) levels in cases with dyslipidemia. Using molecular modeling analysis, we further found that the rs2228570 variant changed the structure and the stability of VDR and altered the binding energy of its ligand.
The VDR rs2228570 variant may increase susceptibility to dyslipidemia in the Chinese Han population.
The prevalence of lipid abnormalities is increasing in the Chinese Han population, and is considered as a major risk factor for cardiovascular disease. Low vitamin D status has been reported to be correlated with an increased risk of hyperlipidemia [1, 2]. Furthermore, a correlation was found between vitamin D deficiency and increased serum total cholesterol (TC) and low-density lipoprotein-cholesterol (LDL-C) levels, and decreased serum high-density lipoprotein-C (HDL-C) levels [3,4,5]. Nonetheless, the mechanism of the lipid-lowering effect of vitamin D remains unknown. Vitamin D likely increases serum calcium by enhancing intestinal calcium absorption. Subsequently, this elevated calcium level reduces serum triglycerides by suppression of hepatic triglyceride formation and secretion . The second possible mechanism may be through the inhibitory effect of vitamin D on serum PTH concentrations. The lowered level of PTH may reduce serum triglycerides by increasing their peripheral uptake. Another potential mechanism to explain the association between 25(OH) D and triglycerides would be through insulin resistance. In cases of vitamin D deficiency, the risk of insulin resistance increases, which is associated with elevated levels of very-low-density lipoprotein (VLDL) cholesterol and triglycerides .
Dyslipidemia is a multifactorial disease, which can be influenced by both genetic and environmental factors. Vitamin D is involved in lipid metabolism. The actions of vitamin D are mediated by the Vitamin D Receptor (VDR), a nuclear receptor . The VDR gene, which codes for VDR, is located at q11–q13 on chromosome 12. A recent study showed that the CC genotype of VDR, single nucleotide polymorphism (SNP) rs2228570, is a risk factor for diabetic dyslipidemia in elderly males in North China . Another investigation, conducted in Lebanon, suggested that VDR polymorphisms were associated with triglyceride and HDL levels in a healthy young population . VDR FokI polymorphisms appear to be associated with coronary heart disease in Han Chinese. Furthermore, the GG genotype predicts higher HDL-cholesterol levels in adults with coronary heart disease . Considering the lack of studies evaluating the link between VDR and dyslipidemia, we investigated and determined the association between VDR variants and dyslipidemia in a large Chinese Han sample.
The study was approved by the Ethics Committees of the First Affiliated Hospital with Nanjing Medical University, Nanjing, Jiangsu, China (NO2015-SR-032). All participants signed written informed consent.
A two-stage designed study was used to evaluate whether genetic variations in VDR are associated with dyslipidemia in the Chinese Han population. In stage I, we genotyped three SNPs (rs11574129, rs2228570, and rs739837) in 877 cases with dyslipidemia from normotensive and non-diabetes mellitus population and 1822 non-dyslipidemia subjects from a community-based epidemiological survey in Jiangsu Province, China. We followed-up the stage I study to confirm the association between rs2228750 and dyslipidemia using a larger sample, including 1679 cases with dyslipidemia and 3124 age-, sex-, and ethnic origin-matched controls, in a stage II study.
Dyslipidemia was identified according to the Guidelines on Prevention and Treatment of Dyslipidemia in Chinese Adults (2007) : triglyceride (TG) ≥ 2.26 mmol/L as high, total cholesterol (TC) ≥ 6.22 mmol/L as high, low-density lipoprotein cholesterol (LDL-C) ≥ 4.14 mmol/L as high, and high-density lipoprotein cholesterol (HDL-C) < 1.04 mmol/L as low. In our study, high TG, high TC, high LDL-C, or low HDL-C was regarded as dyslipidemia. The exclusion criteria included a history of cancer, stoke, thyroid disorder, kidney, or liver disease. Individuals who had family history of obesity, dyslipidemia, and other metabolic diseases were also excluded. All individuals answered standardized questionnaires inquiring about age, sex, ethnicity, history of smoking, and family medical history. Body-mass index (BMI) was calculated as weight divided by height squared in kg/m2. Blood pressure was measured using a mercury sphygmomanometer. LDL-C, HDL-C, TC, TG, and fasting blood glucose (GLU) serum levels were measured enzymatically using an automatic clinical analyzer (Hitachi Inc., Tokyo, Japan). Serum 25(OH) D concentration was determined using an enzyme-linked immunosorbent assay (ELISA) according to the product specifications (Eagle Biosciences Inc., Nashua, New Hampshire, USA).
Genomic DNA extraction and genotyping
Genomic deoxyribonucleic acid (DNA) was extracted from whole-blood samples using the TIANGEN RelaxGene Blood DNA System (TianGen Biotech Co, Lid, Beijing, China; NO. DP319). Three tagSNPs (rs11574129, rs2228570, and rs739837) were selected using the linkage disequilibrium (LD) method and SNP function prediction software (https://manticore.niehs.nih.gov/snpinfo/snpfunc.html). The rs11574129 and rs739837 variants regulating the microRNA binding were positioned in the 3’untranslated region. The rs2228570 was predicted as a splicing enhancer in the exon-coding region. Genotyping was performed using the TaqMan SNP Genotyping Assay (Applied Biosystems, Foster City, CA, USA). For quality control, 5% of samples were re-genotyped in a blinded fashion.
Molecular Modeling Analysis
The VDR molecular modeling studies using Discovery Studio 3.0 (DS 3.0, Accelrys Inc., Boston, MA, USA) provided insight into two VDR isoforms (wild type and RS2228570 variant). These VDR isoforms were produced by I-TASSER (http://zhanglab.ccmb.med.umich.edu/I-TASSER/). We selected the most appropriate structure from five models to model VDR by FG-MD (http://zhanglab.ccmb.med.umich.edu /FG-MD/). Then, we used DS 3.0 to set minimization (steepest descent, conjugate gradient) and equilibration steps, followed by CHARMM to understand the molecular dynamics.
Data analyses were performed using SPSS 19.0 (SPSS, Inc., Chicago, IL, USA). Quantitative variables were expressed as mean ± standard deviation (SD). Comparisons of the means between the two groups were determined by t-test, and the multiple-group means were compared using one-way ANOVA, with adjustments for age and sex. The Chi-square (χ2) test or Fisher’s exact test was employed for categorical variables. Multiple logistic regression analysis was applied to adjust covariates including age and sex. In addition, multivariable linear regression was utilized to evaluate the association of genotypes with lipid levels. Hardy–Weinberg equilibrium (HWE) was assessed by Fisher’s exact χ2 test in the control. Statistical significance was assumed at P < 0.05.
Study subject characteristics
A total of 436 males and 441 females were included in the case group and 926 males and 956 females were included in the control group in the stage I study. The average levels of BMI, LDL, and TC in the case group were significantly higher than those in the controls. Compared to controls, the case groups had a significantly lower level of HDL. However, there was no significant difference in systolic blood pressure (SBP), diastolic blood pressure (DBP), gender, age, and GLU between the case and control groups (Additional file 1: Table S1).
The clinical characteristics of the cases the in stage II study are presented in Additional file 1: Table S2. The case group consisted of 834 males and 845 females, whereas 1515 males and 1609 females were included in the control group. We observed no significant differences in gender and age between the case and control groups. However, the groups differed significantly in their mean BMI, SBP, DBP, GLU, TC, TG, LDL-C, and HDL-C levels. The hypertension prevalence in the case group was significantly higher than that in the control group.
Association of VDR polymorphisms with dyslipidemia
In the present study, the VDR genotype (rs11574129, rs2228570, and rs739837) frequencies were in the Hardy-Weinberg equilibrium (both P > 0.05). In the stage I study, we found that the SNP rs2228570 statistically correlated with an increased risk of dyslipidemia; after adjusting for confounding factors, the odds ratios (ORs) with 95% confidence interval (CI) of the additive, dominant, and recessive models were 1.255 (1.118–1.409) (P < 0.001),1.384 (1.136–1.6)(P = 0.001) and 1.356 (1.1–1.671)(P = 0.004) (Table 1).
Based on the results of the stage I study, we focused on the relationship between the rs2228570 variant and the risk of dyslipidemia in the stage II study. Subsequently, we conducted a second-stage replication by genotyping a larger case-control population. We found that SNP rs2228570 significantly correlated with an increased risk of dyslipidemia; after adjusting for confounding factors, the ORs with 95% CI of the additive, dominant and recessive models were 1.146 (1.053–1247) (P = 0.002), 1.184(1.018–1.376)(P = 0.028) and 1.209 (1.064–1.374)(P = 0.004) (Table 2). To reduce potential bias in the stage II study, we adjusted for confounding factors, such as blood pressure, smoking, BMI, and GLU. In the GLU-adjusted model, we found a significant association between the rs2228570 variant and increased risk of dyslipidemia. After adjusting for confounding factors, the ORs with 95% CI of the additive, dominant, and recessive models were 1.145 (1.052–1.246) (P = 0.002), 1.177 (1.013–1.369) (P = 0.034), and 1.211(1.065–1.377) (P = 0.003). A similar trend was observed in the blood pressure-, smoking-, and BMI-adjusted models (Table 3).
Serum lipid levels in the cases with different genotypes are shown in Fig. 1, in which it can be seen that the rs2228750 genotype was significantly associated with serum LDL levels in the dyslipidemia cases. The cases with the TT genotype had significantly higher LDL levels (4.93 ± 0.75 mmol/L) than those with the TC (4.67 ± 0.47 mmol/L) or CC (4.66 ± 0.44 mmol/L) genotype (P = 0.01). There were no significant differences among the three genotypes in the levels of serum TC, TG, or HDL (P > 0.05) (Additional file 1: Table S3). The levels of serum 25(OH)D among the different genotype cases with high LDL levels (≥ 4.41 mmol/L) is illustrated in Fig. 2. In the cases with high LDL, the cases with the TT genotype had significantly lower serum 25(OH)D levels (18.43 ± 5.04 ng/mL) than those with the TC (26.24 ± 4.16 ng/mL) or CC (36.76 ± 8.10 ng/mL) genotype (P < 0.001) (Additional file 1: Table S4).
Multiple linear regression analysis revealed that the rs2228570 variant (P = 0.041), age (P = 0.01), and sex (P = 0.021) were independently related to the level of serum LDL (Table 4).
Molecular modeling study
The molecular modeling analysis showed that VDR protein stability was reduced in the rs2228570 variant. The binding energy of VDR to its ligand was also decreased in the rs2228570 variant (Additional file 1: Tables S5, S6). As can be observed in Fig. 3, the conformation and interaction changes are in the polymorphic region (rs2228570), where the replacement of threonine with methionine at amino acid position 1 changed the structure of VDR. This modification resulted in a shift of the mode in which the amino acid residue Tyr236 in the protein VDR was bound to hydrogen in its ligand. The change also provided interactions with neighboring residues.
The active form of vitamin D is 1, 25-dihydroxyvitamin D. VDR tightly regulates vitamin D levels in cells and in the circulation . A previous study  revealed that VDR activation repressed the increased levels of mouse and human CYP7A1 and reduced the cholesterol concentrations exerted by hepatic small heterodimer partner (SHP). Recent research found that the deletion of macrophage VDR promoted insulin resistance and monocyte cholesterol transport, which accelerated atherosclerosis in mice . Based on the findings of the abovementioned studies, we hypothesized that VDR expression is associated with serum lipid levels. Hence, further studies on VDR polymorphisms would improve the understanding of the association of VDR with dyslipidemia.
In our study, we found a significant association of rs2228570 polymorphism with dyslipidemia and serum LDL levels. Moreover, cases with the TT genotype had significantly lower serum 25(OH)D levels compared to cases with other genotypes with high LDL.
As is well known, vitamin D regulates lipid metabolism; it inhibits lipid synthesis and promotes adipopexis in adipocytes . However, the underlying mechanisms remain largely unknown. It is possible that the potential underlying mechanism operates through facilitating calcium homeostasis via vitamin D. It is commonly known that vitamin D and parathyroid hormone (PTH) are the two most important hormones that maintain calcium levels [17, 18]. Therefore, the reduction in lipid concentrations observed in this study may have been caused by increased calcium levels due to a decrease in fatty acid absorption and an increase in fecal fatty acid content, probably arising from the formation of insoluble calcium-fatty soaps in the gut [19, 20]. Such reduced absorption of fat diminishes serum total and LDL cholesterol levels [21, 22]. Calcium can also bind bile acid [23, 24] and increase the conversion of cholesterol to bile acids, thus causing increased cholesterol excretion . Furthermore, microsomal triglyceride transfer protein (MTP) is involved in the secretion of triglyceride-rich VLDL, which leads to hypertriglyceridemia and atherosclerosis . In addition, hepatocellular calcium activates MTP. The increase in hepatocellular calcium could be restrained by increasing serum concentrations . The inhibition of hypertriglyceridemia may be caused by the inhibition of MTP activity via the suppression of hepatocellular calcium. Another potential mechanism is likely to be associated with a decrease in cholesterol absorption . Vitamin D increases intestinal calcium absorption by affecting the rate of fat absorption, which results from the ability of calcium to bind to fatty and bile acids, forming calcium-fatty acid soaps , which are excreted in the feces. Moreover, vitamin D may play a role in insulin secretion and sensitivity, thereby indirectly influencing lipid metabolism .
Using our molecular model, we further established that rs2228570 changed the structure and stability of VDR and altered the binding energy of its ligand. Earlier investigations revealed that an association between low levels of vitamin D and rs2228570 located in VDR’s 5’end. rs2228570 has a T to C variation in the translation initiation codon (AGT) in exon 2 . This change induces the synthesis of a small (49.5 kD) protein with increased biological activity . The reason for the differences between the activities of the two proteins, such as their abilities to bind 1,25-dihydroxy vitamin D3 or induce transcription, still remain to be determined.
In conclusion, our results suggest that the rs2228570 polymorphism correlates with an increased risk of higher LDL and lower serum 25(OH)D levels in Chinese patients. The variation in rs2228750 reduced the stability of VDR and changed the binding energy of its ligand, which may be the underlying molecular mechanism for variant associated dyslipidemia. Our findings provide new insights into the prevention and treatment of dyslipidemia. Nevertheless, further studies are required to verify these results in different ethnic groups. Finally, additional functional investigations of the VDR gene are needed to confirm the present findings.
Low density lipoprotein cholesterol
Single nucleotide polymorphisms
Vitamin D receptor
Botella-Carretero JI, Alvarez-Blasco F, Villafruela JJ, et al. Vitamin D deficiency is associated with the metabolic syndrome in morbid obesity. Clin Nutr. 2007;26:573–80.
Kendrick J, Targher G, Smits G, et al. 25-Hydroxyvitamin D deficiency is independently associated with cardiovascular disease in the Third National Health and Nutrition Examination Survey. Atherosclerosis. 2009;205:255–60.
Forman JP, Curhan GC, Taylor EN. Plasma 25-hydroxyvitamin D levels and risk of incident hypertension among young women. Hypertension. 2008;52:828–32.
Schwartz JB. Effects of vitamin D supplementation in atorvastatin-treated patients: a new drug interaction with an unexpected consequence. Clin Pharmacol Ther. 2009;85:198–203.
Skaaby T, Husemoen LL, Pisinger C, et al. Vitamin D status and changes in cardiovascular risk factors: a prospective study of a general population. Cardiology. 2012;123:62–70.
Martins D, Wolf M, Pan D, et al. Prevalence of cardiovascular risk factors and the serum levels of 25-hydroxyvitamin D in the United States: data from the Third National Health and Nutrition Examination Survey. Arch Intern Med. 2007;167(11):1159–65.
Ginsberg HN, Zhang Y-L, Hernandez-Ono A. Regulation of plasma triglycerides in insulin resistance and diabetes. Arch Med Res. 2005;36(3):232–40.
Jia J, Shen C, Mao L, et al. Vitamin D receptor genetic polymorphism is significantly associated with decreased risk of hypertension in a Chinese Han population. J Clin Hypertens (Greenwich). 2014;16:634–9.
Xia Z, Hu Y, Han Z, et al. Association of vitamin D receptor gene polymorphisms with diabetic dyslipidemia in the elderly male population in North China. Clin Interv Aging. 2017;12:1673–9.
Hajj A, Chedid R, Chouery E, et al. Relationship between vitamin D receptor gene polymorphisms, cardiovascular risk factors and adiponectin in a healthy young population. Pharmacogenomics. 2016;17(15):1675–86.
He L, Wang M. Association of vitamin d receptor-a gene polymorphisms with coronary heart disease in Han Chinese. Int J Clin Exp Med. 2015;8(4):6224.
Joint Committee for Developing Chinese guidelines on Prevention and Treatment of Dyslipidemia in Adults. Chinese guidelines on prevention and treatment of dyslipidemia in adults. Zhonghua Xin Xue Guan Bing Za Zhi. 2007;35(5):390–419.
García-Bailo B, Jamnik J, Da Costa LA, et al. Genetic variation in the vitamin D receptor, plasma 25-hydroxyvitamin D, and biomarkers of cardiometabolic disease in Caucasian young adults. J Nutrigenet Nutrigenomics. 2013;6(4-5):256–67.
Chow EC, Magomedova L, Quach HP, et al. Vitamin D receptor activation down-regulates the small heterodimer partner and increases CYP7A1 to lower cholesterol. Gastroenterology. 2014;146:1048–59.
Oh J, Riek AE, Darwech I, et al. Deletion of macrophage Vitamin D receptor promotes insulin resistance and monocyte cholesterol transport to accelerate atherosclerosis in mice. Cell Rep. 2015;10(11):1872–86.
Jorde R, Grimnes G. Vitamin D and metabolic health with special reference to the effect of vitamin D on serum lipids. Prog Lipid Res. 2011;50:303–12.
Rostand SG, Drueke TB. Parathyroid hormone, vitamin D, and cardiovascular disease in chronic renal failure. Kidney Int. 1999;56:383–92.
Steingrimsdottir L, Gunnarsson O, Indridason OS, et al. Relationship between serum parathyroid hormone levels, vitamin D sufficiency, and calcium intake. JAMA. 2005;294:2336–41.
Denke MA, Fox MM, Schulte MC. Short-term dietary calcium fortification increases fecal saturated fat content and reduces serum lipids in men. J Nutr. 1993;123:1047–53.
Reid IR. Effects of calcium supplementation on circulating lipids: potential pharmacoeconomic implications. Drugs Aging. 2004;21:7–17.
Grundy SM, Denke MA. Dietary influences on serum lipids and lipoproteins. J Lipid Res. 1990;31:1149–72.
Vaskonen T. Dietary minerals and modification of cardiovascular risk factors. J Nutr Biochem. 2003;14:492–506.
Saunders D, Sillery J, Chapman R. Effect of calcium carbonate and aluminum hydroxide on human intestinal function. Dig Dis Sci. 1988;33:409–13.
Van der Meer R, Welberg JW, Kuipers F, et al. Effects of supplemental dietary calcium on the intestinal association of calcium, phosphate, and bile acids. Gastroenterology. 1990;99:1653–9.
Vaskonen T, Mervaala E, Sumuvuori V, et al. Effects of calcium and plant sterols on serum lipids in obese Zucker rats on a low-fat diet. Br J Nutr. 2002;87:239–45.
Cho HJ, Kang HC, Choi SA, et al. The possible role of Ca2+ on the activation of microsomal triglyceride transfer protein in rat hepatocytes. Biol Pharm Bull. 2005;28:1418–23.
Boon N, Hul GB, Stegen JH, et al. An intervention study of the effects of calcium intake on faecal fat excretion, energy metabolism and adipose tissue mRNA expression of lipid-metabolism related proteins. Int J Obes. 2007;31:1704–12.
Kamycheva E, Jorde R, Figenschau Y, et al. Insulin sensitivity in subjects with secondary hyperparathyroidism and the effect of a low serum 25-hydroxyvitamin D level on insulin sensitivity. J Endocrinol Investig. 2007;30:126–32.
Saijo T, Ito M, Takeda E, et al. A unique mutation in the vitamin D receptor gene in three Japanese patients with vitamin D-dependent rickets type II: utility of single-strand conformation polymorphism analysis for heterozygous carrier detection. Am J Hum Genet. 1991;49:668–73.
Arai H, Miyamoto K, Taketani Y, et al. A vitamin D receptor gene polymorphism in the translation initiation codon: effect on protein activity and relation to bone mineral density in Japanese women. J Bone Miner Res. 1997;12:915–21.
We express our sincere gratitude to the Department of Epidemiology and Biostatistics at the School of Public Health of Nanjing Medical University for providing us with the necessary resources and administration support throughout the study.
The work was supported by grants provided by the National Natural Science Foundation of China (No.81570810) and the Technology Project of Jiangsu Province (BL2014085). This work was also supported by a Project funded by Jiangsu Province Health Commission (H201402).
Availability of data and materials
All data generated or analysed during this study are included in this MS and its Additional file 1.
Ethics approval and consent to participate
The experimental protocol was approved by the Institutional Ethics Board of the First Affiliated Hospital with Nanjing Medical University (NO2015-SR-032). All patients signed written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Clinical characteristics of cases and controls (stage I). Table S2. Clinical characteristics of cases and controls (stage II). Table S3. Analysis of rs2228570 association with lipid levels. Table S4. Analysis of the association of rs2228570 and serum 25(OH)D levels in cases with high LDL. Table S5. Changes in VDR stability. Table S6. Changes in binding energy of VDR. (DOC 241 kb)
About this article
Cite this article
Jia, J., Tang, Y., Shen, C. et al. Vitamin D receptor polymorphism rs2228570 is significantly associated with risk of dyslipidemia and serum LDL levels in Chinese Han population. Lipids Health Dis 17, 193 (2018) doi:10.1186/s12944-018-0819-0
- Vitamin D receptor