Polymorphisms of the Vitamin D Receptor Gene and Sex-Differential Associations with Lipid Profiles in Chinese Han Adults＊

CHEN Yan Mei , XU Ping , WANG Zhou Tian , ZHU Yu Mei , GONG Chun Mei ,HUANG Chang Hua , LIU Xiao Li, and ZHOU Ji Chang,,4,#

1.School of Public Health (Shenzhen),Sun Yat-sen University,Shenzhen 518100,Guangdong,China;2.Shenzhen Center for Chronic Disease Control,Shenzhen 518020,Guangdong,China;3.Shenzhen Qilinshan Sanatorium,Shenzhen 518055,Guangdong,China;4.Guangdong Province Engineering Laboratory for Nutrition Translation,Guangzhou 510080,Guangdong,China

Abstract Objective To explore the association of single nucleotide polymorphisms (SNPs) of the vitamin D receptor gene (VDR) with circulating lipids considering gender differences.Methods Of the Han Chinese adults recruited from a health examination center for inclusion in the study,the circulating lipids,25-hydroxyvitamin D (25OHD),and other parameters were measured.The VDR SNPs of Cdx2 (rs11568820),Fok1 (rs2228570),Apa1 (rs7975232),and Taq1 (rs731236) were genotyped with a qPCR test using blood DNA samples,and their associations with lipids were analyzed using logistic regression.Results In the female participants (n=236 with dyslipidemia and 888 without dyslipidemia),multiple genotype models of Fok1 indicated a positive correlation of B (not A) alleles with LDLC level (P ＜ 0.05).In the male participants (n=299 with dyslipidemia and 564 without dyslipidemia),the recessive model of Cdx2 and the additive and recessive models of Fok1 differed (P ＜ 0.05) between the HDLC-classified subgroups,respectively,and Fok1 BB and Cdx2 TT presented interactions with 25OHD in the negative associations with HDLC (P ＜ 0.05).Conclusion In the Chinese Han adults included in the study,the Fok1 B-allele of VDR was associated with higher LDLC in females,and the Fok1 B-allele and the Cdx2 T-allele of VDR were associated with lower HDLC in males.The interaction of VD and Fok1 BB or Cdx2 TT in males synergistically decreased HDLC levels.

Key words:Vitamin D;Vitamin D receptor;Gene polymorphism;Lipid;Dyslipidemia

INTRODUCTION

Vitamin D (VD) is a micronutrient for human health,which is predominantly derived from 7-dehydrocholesterol in the skin under sunlight exposure and is also absorbed from food as a secondary source.Through two steps of hydroxylation,VD is converted to 25-hydroxyvitamin D (25OHD) and then 1,25-dihydroxyvitamin D [1,25(OH)2D].The circulating concentration of 25OHD is generally used as a biomarker of VD status,and 1,25(OH)2D is the functional metabolite working in the manner of a hormone in the osseous and extra-osseous metabolisms[1].Previous studies suggested that VD deficiency was associated with abdominal obesity and/or dyslipidemia[2],and many studies supported the opinion that VD supplementation was beneficial for ameliorating lipid profiles[3].However,confounders in various populations led to inconsistency for each of the lipid parameters,namely triglyceride (TG),total cholesterol (TC),lowdensity lipoprotein cholesterol (LDLC),and highdensity lipoprotein cholesterol (HDLC),as suggested by human studies and meta-analyses[3,4].Researchers have to conduct a close examination of those potential confounders and pay special attention to the functioning of the VD receptor (VDR),which is the most important mediator of the VD function[5].

VDRgene variants have the potential to influence biological outcomes[5,6].It has been reported that VDR has an effect on skeletal muscles independent of VD,and the decreased expression of VDR protein is related to various disease states and aging[7].In the VDR knockout animal model,a decreased size and muscle fiber strength were observed,and VDR expression is also required in myoblasts[8,9].Aside from skeletal health status,it was also reported that the single nucleotide polymorphisms (SNPs) ofVDRwere associated with risks to non-skeletal health[6],such as myocardial infarction,cancer,and death[10].Moreover,these kinds of associations aroused concerns about dyslipidemia since expressions of some lipid metabolism-related genes are regulated by calcitriolstimulation[11],and some may even have VDresponsive elements (VDREs) to bind the 1,25(OH)2Dactivated VDR[12].Thus,both the VD status and the VDR function have the potential to affect the expressions of lipid metabolic genes,ultimately affecting the physiological activities of cells[1].Moreover,by involving biological (in terms of reproductive function,sexual hormones,expressions of heterosome-located genes,etc.) and lifestyle (e.g.,smoking,drinking,physical activity,etc.) factors[13],sexual dimorphism should be taken into careful consideration when conducting medical research[14].

Located on chromosome 12,theVDRgene spans at least 75 kb in length.From the 5’-3’ direction of its sense (forward) strand,Cdx2(rs11568820,C ＞ T),Fok1[rs2228570,A ＞ B (degenerate base standing for C,G,and T,i.e.,not A)],Apa1(rs7975232,C ＞ A),andTaq1(rs731236,A ＞ G) are the generally the relevant SNPs in various ethnicity groups[15-17].These SNPs may affect mRNA expression,amino acid sequence,and protein activity.Particularly,VDRSNPs have been reported to have effects on 25OHD levels in circulation[18,19].For example,Santos et al.[20]studied a group of Brazilian girls aged 7-18 years and found that theVDRwild-types ofBsm1,Apa1,andTaq1were associated with lower 25OHD levels.Therefore,it is necessary to further explore the relationship betweenVDRSNPs and specific lipid parameters,as well as VD status in both females and males.In the present study,we focused on theVDRSNPs ofApa1,Cdx2,Fok1,andTaq1in this regard,which have been studied in genotype-phenotype associations,in order to avoid weak evidence or unknown biological significance for less studied or new SNPs regarding their associations with lipids.

METHODS

Subjects and Sample Collection

The sample size (N) was estimated with the formulaN≥ (Z1-α/2/δ)2×p× (1 -p),whereinZwas 1.96 for the two-sided 95% confidential intervals(CIs),Pwas the dyslipidemia prevalence of 34.64% in Shenzhen adults as known in 2012 but unpublished at the design stage of our study[21],δwas 2.2%,representing tolerable error,and the calculatedNwas 1,798.

After signing written consent forms to be subjects in the study,adult volunteers visiting a health examination center in Shenzhen City in Guangdong Province,China,were recruited between July 2013 and January 2014.A questionnaire survey administeredviaface-to-face interview was conducted to collect the basic health information of each participant.The inclusion criteria were:1) Han Chinese aged ≥ 20 years old;2) living in Shenzhen for＞ 2 years;3) free of liver diseases,renal diseases,and any cancers in the past six months;and 4) not pregnant (for women).The excluded subjects were those who 1) had severe organic diseases,2) had acute infection symptoms,allergic diseases,and malignant tumors,3) had a family history of genetic diseases including familial dyslipidemia,4) had taken VD supplements in the past six months,5) had taken medicines to control lipid levels within the past 12 hours,or 6) had taken diuretics,engaged in strenuous exercise,or had overeaten within the 24 hours before the test.

From ulnar veins of the subjects,who had fasted overnight,blood samples were collected with vacuum tubes,and within two hours post-collection,supernatant and blood cells were separated by centrifugation of 3,000 ×gat 4 °C for 10 min.The body mass index (BMI) was calculated as body weight (kg)/ height (m)2.The above protocol for the cross-sectional study was approved by the Ethics Committee of the Shenzhen Center for Chronic Disease Control and was in accordance with the Declaration of Helsinki.The study was registered at ClinicalTrials.gov(registration ID number NCT04707612).

Index Measurement

Parameters of lipids,glucose,hemogram,etc.,were assayed right on the day,and the rest of the samples were aliquoted with EP tubes and stored at-80 °C for later use.The systolic blood pressure (SBP)and diastolic blood pressure (DBP) of the subjects were measured by a standard mercury sphygmomanometer.Glycated hemoglobin (HbA1c)was detected with a Japan Arkray Instruments analyzer (ion chromatography method).The serum levels of TG,TC,LDLC,HDLC,and fasting plasma glucose (FPG) were tested by the Beckman-LX20 automatic biochemical analyzer.The plasma 25OHD concentration was detected with a commercial ELISA kit (Cat.#:AC-57F1,IDS Ltd.,UK).According to the Chinese guidelines for dyslipidemia management[22],TG ≥ 2.3 mmol/L,TC ≥ 6.2 mmol/L,LDLC ≥4.1 mmol/L,HDLC ＜ 1.0 mmol/L,or previously diagnosed dyslipidemia,were the criteria used to define the dyslipidemia group (DL),and those subjects who did not match the DL criteria were assigned to the non-dyslipidemia group (ND).

In the logistic analyses,from the viewpoint of preventive medicine,as well as to avoid the small subject numbers for the genotypes of low minor allele frequency (MAF)[23],the available marginallyelevated cut-off values for TG (≥ 1.7 mmol/L),TC (≥5.2 mmol/L),and LDLC (≥ 3.4 mmol/L),as well as the only available cut-off value for HDLC (＜ 1.0 mmol/L),were used to define the abnormal subgroups for the corresponding lipid parameters.

DNA Preparation and SNP Genotyping

The protocols for DNA preparation and SNP genotyping described in our previous studies[23,24]were also adopted in the present study.Briefly,an asymmetric amplification with molecular beaconbased real-time quantitative PCR (MB-qPCR)followed by a melting step was performed to determine the variation of the SNP loci.Before its application to the large sample size analysis,the method was verified with the gold standard of Sanger sequencing (ThermoFisher,Shanghai,China).The information of primers and molecular beacons for the genotyping experiment with qPCR (Roche 480II,Singapore) is presented in Supplementary Table S1,available in www.besjournal.com.

The nucleotide bases used in our study were consistent with those in dbSNP of NCBI,i.e.,C ＞ T for rs11568820 (Cdx2),A ＞ B (not A) for rs2228570(Fok1),C ＞ A for rs7975232 (Apa1),and A ＞ G for rs731236 (Taq1).If the complementary bases [G ＞ A forCdx2,T ＞ V (standing for A,C,and G) forFok1,G ＞ T forApa1,and T ＞ C forTaq1],or the letter standing for the DNA-digestibility with the corresponding restriction enzyme (f＞FforFok1,a＞AforApa1,andT＞tforTaq1),were used in the references,the original usages were provided as bracketed annotations in our article.

Data Analysis

The statistical analyses were performed using SPSS for Windows version 25 (IBM Corp.,Armonk,NY,USA).Clinical data were presented as medians and interquartile ranges (25% to 75%) and compared with rank-sum tests (Wilcoxon rank test or Kruskal-Wallis rank test).The genotypes of the four SNPs were tested with Hardy-Weinberg equilibrium (HWE)analyses for sampling representation.For genotypic comparisons,differences in allele and genotype frequencies were evaluated using the chi-square (χ2)test.The homozygous genotypes of CC forCdx2,AA forFok1,CC forApa1,and AA forTaq1were used as reference genotypes,respectively.For both genders,the additive,dominant,recessive,homozygous,and allelic models for each of the SNPs were entered into the logistic regression analyses for odds ratios (ORs)and 95%CIs with adjustment for age,BMI,FPG,HbA1c,25OHD,SBP,and DBP.In order to analyze the interaction between VD nutritional levels andVDRSNPs,the interaction factor calculated by 25OHD concentrations andVDRSNPs were entered into the logistic regression with adjustment for those confounders mentioned above.A simple linear regression was used to assess the association between circulating 25OHD and lipid profiles.APvalue of less than 0.05 was considered to be statistically significant.

RESULTS

Clinical Profiles of the Participants

A total of 1,987 adults were included in the analysis,aged from 20 to 81 years old.Among them,there were 1,124 females aged 39 (31-49) (median and interquartile range) years old and 863 males aged 36 (30-44) years old.The clinical profiles of the participants were divided into DL and ND for each gender,as summarized in Table 1.It was indicated that several metabolic or metabolism-related parameters were statistically different between DL and ND adults (P＜ 0.05),such as age,BMI,FPG,HbA1c,SBP,DBP,TG,TC,LDLC,HDLC,etc.In particular,25OHD concentrations were statistically different between DL and ND groups in men but not in women.

Genotypic Frequencies of VDR

The results of the genotyping experiment on the SNPs with MB-qPCR are shown in Supplementary Figure S1,available in www.besjournal.com.The HWE test showed that the population sampled from the Health Examination Center was representative of the population at large (Supplementary Table S2,available in www.besjournal.com).The genotypic and allelic frequencies ofVDRSNPs between DL and ND groups are summarized in Supplementary Table S3,available in www.besjournal.com.In the female participants,both genotypic and allelic frequencies ofFok1showed differences between DL and ND groups (P＜ 0.05);while in the male participants,these differences were not observed.The genotypic and allelic frequencies of the other three SNPs,namelyCdx2,Apa1,andTaq1,did not display significant differences between the DL and ND groups for either gender (P≥ 0.05),and neither did their different genetic models (Supplementary Tables S4-S6,available in www.besjournal.com).

Logistic Regression Analyses of the Relationship between VDR SNPs and Dyslipidemia

Logistic regression analyses were performed with adjustment for age and BMI since the adjustment for 25OHD did not change the overall correlations (see Supplementary Tables S7-S8). For female participants (see Table 2),between the subgroups of LDLC (≥ 3.4vs.＜ 3.4 mmol/L),VDR Fok1presented significance (P＜ 0.05) in its additive (BBvs.ABvs.AA,P=0.03,OR=1.28,95%CI:1.03-1.59),recessive(BBvs.AB+AA,P=0.03,OR=1.44,95%CI:1.04-2.00),homozygous (BBvs.AA,P=0.02,OR=2.89,95%CI:1.18-7.05),and allelic (Bvs.A,P=0.04,OR=1.26,95%CI:1.01-1.59) models,but between the subgroups of TC (≥ 5.2vs. ＜5.2 mmol/L) and HDLC (＜ 1.0vs.≥ 1.0 mmol/L),no significance (P≥ 0.05) was found in all models of the four SNPs.

Table 3 displays the logistic regression analyses of the relationships between the lipids and the SNPs ofVDRin the male participants.Between the subgroups of HDLC (＜ 1.0vs.≥ 1.0 mmol/L),the data presented significance (P＜ 0.05) in the recessive model ofCdx2(TTvs.CT+CC,P=0.04,OR=2.70,95%CI:1.08-6.80)and the additive (BBvs.ABvs.AA,P=0.02,OR=2.25,95%CI:1.12-4.52) and recessive (BBvs.AB+AA,P=0.02,OR=2.73,95%CI:1.11-6.70) models ofFok1,while other models of the genotype showed no significance (P≥ 0.05).Differences were not found forApa1andTaq1between the subgroups defined by any of the lipid profiles (P≥ 0.05).

Table 1. Clinical profiles of study subjects (medians and interquartile ranges)

Comparison of 25OHD Concentration

No differences were found in the comparison of25OHD concentrations across genotype models ofVDRgene polymorphisms in both genders(see Supplementary Table S9,available in www.besjournal.com). When the plasma 25OHD concentrations in subjects with abnormal and normal lipid profiles were compared,it was shown that women with low HDLC had higher 25OHD while men with high TG had lower 25OHD than the subjects with normal lipid parameters (P＜ 0.05,Supplementary Table S10,available in www.besjournal.com).

Table 2. Logistic regression analyses of lipids with SNPs of VDR in adult females with adjustment for age and body mass index

Table 3. Logistic regression analyses of lipids with SNPs of VDR in adult males with adjustment for age and body mass index

Relationship between Circulating 25OHD and Lipid Profiles

In men,lg TG (β=-0.003,P＜ 0.001),LDLC (β=-0.003,P=0.006),and HDLC (β=-0.001,P=0.025)were inversely associated with the 25OHD concentration while in women,none of the lipid profiles showed a significant correlation with the 25OHD concentration (P≥ 0.05,Supplementary Table S11,available in www.besjournal.com).

Interaction between VD Nutritional Status and VDR SNPs on the Occurrence of Dyslipidemia

According to the logistic regression results (see Tables 2 and 3) and the simple linear regression results (see Supplementary Table S11),VD andVDRSNPs had some parameters that were related to the results of dyslipidemia. Therefore,we further explored whether VD andVDRSNPs were synergistic or antagonistic to the results of dyslipidemia.The results (see Table 4) further showed that in men,the combined effect of VD,Fok1BB (P=0.02,OR=1.03,95%CI:1.01-1.06),andCdx2TT (P=0.02,OR=1.03,95%CI:1.01-1.06) synergistically led to a decrease in HDLC levels while in women,there was no significant difference in the interaction between VD andFok1orCdx2(OR＜ 1 andP≥ 0.05) on the lipid profiles.

Supplementary Table S1. The information of molecular beacon probes and primers for genotyping the four single nucleotide polymorphisms (SNPs) of vitamin D receptor gene (VDR)

Supplementary Table S2. Hardy-Weinberg equilibrium tests on the vitamin D receptor gene polymorphisms,Cdx2 (rs11568820),Fok1 (rs2228570),Apa1 (rs7975232),and Taq1 (rs731236)

Supplementary Table S3.Comparisons of the genotypic and allelic frequencies of vitamin D receptor gene polymorphisms,Cdx2 (rs11568820),Fok1 (rs2228570),Apa1 (rs7975232),and Taq1 (rs731236) between dyslipidemia (DL) and non-dyslipidemia (ND) adults

Supplementary Table S4.Genetic model comparison of Cdx2 (rs11568820) in vitamin D receptor gene between dyslipidemia (DL) and non-dyslipidemia (ND) adults

Supplementary Table S5.Genetic model comparison of Apa1 (rs7975232) in vitamin D receptor gene between dyslipidemia (DL) and non-dyslipidemia (ND) adults

Supplementary Table S6.Genetic model comparison of Taq1 (rs731236) in vitamin D receptor gene between dyslipidemia (DL) and non-dyslipidemia (ND) adults

Supplementary Table S7. Logistic regression analyses of lipids with SNPs of VDR in adult females with adjustment for age,body mass index,dietary 25-hydroxyvitamin D,and the sunshine time per week

Supplementary Table S8. Logistic regression analyses of lipids with SNPs of VDR in adult males with adjustment for age,body mass index,dietary 25-hydroxyvitamin D,and the sunshine time per week

Supplementary Table S9. Comparison of 25-hydroxyvitamin D concentration (nmol/L) across genotype models of VDR gene polymorphisms in adult males and females

Continued

Supplementary Table S10.Comparison of 25-hydroxyvitamin D concentration (nmol/L) in abnormal and normal lipid groups of adult males and females

Supplementary Table S11.Simple linear regressions between serum 25-hydroxyvitamin D and lipid profiles

Supplementary Figure S1. Melting curve analyses on the polymorphisms of vitamin D receptor gene(VDR).Panel A:Melting curve of VDR Cdx2 (rs11568820,A).Panel B:Melting curve of Fok1 (rs2228570,B).Panel C:Melting curve of VDR Apa1 (rs7975232,C).Panel D:Melting curve of VDR Taq1 (rs731236,D).

Table 4.Interactions between plasma 25OHD concentrations and SNPs of VDR,Cdx2 (rs11568820) and Fok1(rs2228570) on the occurrence of dyslipidemia

DISCUSSION

Cellular and animal studies have demonstrated that both VD and VDR play important roles in adipocyte differentiation[25,26].VD activates VDR to inhibit the differentiation of pre-adipocytes into mature adipocytesin vitroat an early stage[27,28],and gene expression data showed that the abundance ofVDRmRNA was changed in the process of adipocyte differentiation[29-31].InVdrknockout mice,VD could not block the expression of peroxisome proliferatoractivated receptor γ (PPARγ) and the corresponding adipocyte differentiation,which indicated that VDR was a key mediator of VD in adipocyte differentiation[27,28].The TC level ofVdrknockout mice was more than 20% higher than that of wildtype mice[32].Therefore,the existence of VDR is very important for the healthy regulation of serum lipid profiles,and it is speculated that VDR has the function of lowering TC levels.Knockout ofVdrin macrophages will cause insulin resistance and promote the movement of cholesterol,accelerating atherosclerosis in mice[33].The above all suggest that the perfect function of VDR has a positive effect on the normal lipid metabolism process.

As SNPs of a gene may impact gene expression,as well as the structure and function of the encoded protein and the consequent phenotypes,VDRpolymorphisms have been found to modify the risk of metabolic diseases[34].Among the SNPs ofVDR,Fok1is the most widely studied.Some studies have found that the AA (complementary to TT,i.e.,ff)carriers ofFok1had higher TC in Arab adults[35],lower HDLC and higher TG in Moroccans[36],andhigher LDLC in Chinese individuals[37].Moreover,Fok1A allele carriers had a higher risk of having coronary heart disease,and patients with a BB (orFF,e.g.,GG) genotype had a higher level of HDLC[38].The genetic benefits of AA (ff) or risk of BB (FF) or their alleles were also reported for other lipid parameters,or even the same lipid parameters,in different studies.In adult male Poles,the AA (ff)genotype carriers had lower fasting insulin and higher HDLC[39].Similarly,in the men in our study,the AA (ff) genotype carriers had higher HDLC.In the Han Chinese population,the risk of dyslipidemia was associated with the BB (e.g.,CC) genotype in the elderly male T2D patients of Beijing[40]and with the B(e.g.,C) allele in a community-based population in Nanjing[37].Largely consistent with these findings,we found that the B allele was positively associated with LDLC in women and that the BB genotype was negatively associated with HDLC in men.

The other threeVDRSNPs,Cdx2,Apa1,andTaq1,were also reported to be associated with lipid profiles. For example,regarding theCdx2polymorphism ofVDRin a Lebanese student cohort,CC (GG) or CT (AG) carriers of both genders had higher LDLC than the TT (AA) carriers[41].Higher TC and LDLC levels were also observed in theApa1A allele (T orA) (Russian)[42]andApa1CA (GT)(Chinese)[43]carriers than in the CC (GG) carriers.TheVDRTaq1AA (TT) genotype was associated with increased TG and HDLC levels compared to the genotypes of GG (CC) and GA (CT) in obese Greek subjects[44].However,Taq1GA and GG genotypes were associated with higher TC and LDLC levels in T2D patients (Arab)[45],and the GA genotype exhibited a higher TC concentration compared with the AA genotypes (Arab)[46].However,Karonova and colleagues[42]noted that no difference was found regarding the association betweenTaq1and any lipid parameter.These findings suggested that the associations betweenVDRSNPs and the individual lipid parameters differed among various races or health conditions.

In general,gender differences and sexual dimorphism also affect the genotype and phenotype association.Studies from Lebanon,Poland,and Brazil concluded that men withFok1B (F) had lower HDLC and higher TG levels[39,41,47]while female students with the A allele (f) had higher TC and LDLC[41].By analyzing genotypic frequencies of the fourVDRSNPs between the DL and ND groups in our study,we found thatFok1showed significant differences between the two groups.Adjusting for age and BMI,we found thatVDR Fok1B (C,G,and T)alleles were related to the occurrence of high LDLC in the female group and low HDLC in the male group.Moreover,VDR Cdx2T was related to low HDLC in the male group. Gender-specific differences betweenVDRSNPs and various cardiovascular risk factors and adiponectin[41]also suggest thatVDRSNPs may be a predictor of cardiovascular diseases.Thus,based on our findings,Fok1B in both sexes andCdx2T in men may serve as predictors for dyslipidemia and related diseases,including cardiovascular diseases.

Though it has been widely acknowledged that plasma or serum 25OHD levels are inversely correlated with TG,TC,and/or LDLC[14,48-53],the relationship between VD and HDLC is complicated[3].We previously found that oral supplementation of VD (50,000 IU/week for eight weeks) induced a decrease in HDLC in male adults[54].In interpreting the discrepancy in gender,age,and physiological status that might be involved,we found that sunlight exposure also seemed to affect the association between 25OHD and HDLC.Specifically,at lower sunlight exposure (＜ 1 h/d),a positive relationship was observed,while at higher sunlight exposure (＞ 2 h/d),25OHD was negatively associated with HDLC[55].In our present study,the men in the DL group had a lower level of VD than the men in the ND group,and plasma 25OHD was inversely associated with lg TG,LDLC,and HDLC in men.There was no significant difference in the distribution of VD concentration among the variousVDRSNPs genotypes,but we noticed that the women in the DL group were much older than the women in the other groups,which may affect the VD levels among the groups.Other possible reasons,such as the population selection,dietary VD,time spent outdoors for sun exposure,menopause,etc.,may also have complicated the results.

In addition,VDRpolymorphisms may affect the responses of metabolic parameters to VD supplementation.After VD supplementation,Apa1CC (aa) as compared toApa1AA (AA) orCdx2TT(AA) as compared toCdx2CC (GG) was associated with a greater decrease in plasma LDLC[56].A previous study conducted in New Zealand found that after VD supplementation in women with VD deficiency,theVDRSNP genotypes ofApa1AA (AA)andTaq1GG (tt) were predicted to show greater improvement in insulin resistance[57].Thus,it was easy to understand that different genotypes ofVDRmay have an impact on VD levels and biochemical metabolic indicators[58-60].For example,theFok1AA (TT,i.e.,ff) genotype was associated with lower 25OHD levels and higher LDLC levels in the aforementioned Chinese population[37].Nevertheless,we failed to find such differences in the plasma 25OHD concentration between the genotypes for each of the four SNPs in our examination center-based population.According to the logistic regression and simple linear regression results,some parameters of VD andVDRSNPs were related to results of dyslipidemia.Therefore,we further explored whether VD andVDRSNPs were synergistic or antagonistic to the results of dyslipidemia.The results showed that in men,the combined effect of VD andFok1BB orCdx2TT did not optimize HDLC parameters but synergistically decreased HDLC levels,which tend to be unhealthy.This suggests that,in men,in particular those who carry theFok1BB orCdx2TT genotypes,supplementation of VD may not help improve their low-HDLC dyslipidemia;however,further studies are needed to confirm this.

Aside from the inconsistent results of the epidemiological evidence due to complicated confounders existing in real life,the fundamental biological aspects of some SNPs ofVDRhave been interpreted.TheFok1restriction enzyme site as GGATGN(9)^ at the 5'-end ofVDRin exon 2 has two ATG triplets that can serve as the initial codon for the long or short isoforms of the VDR protein,respectively. The impaired nuclear uptake of 1,25(OH)2D in fibroblasts was related to the conversion of the A (complementary to T in the first ATG triplet ofVDR) allele to B (e.g.,C) for theFok1polymorphism[61],which initiated the translation from the downstream ATG for the smaller VDR protein (49.5 kD),three amino acid residues shorter than the longer form.Further analysis of the VDR molecular model found that theFok1B-allele changed the structure and stability of VDR and its binding energy with its ligand[37]but is more effective in activating the VD target genes[62].However,the expression and migration efficiency of the two VDR isoforms have not yet been determined.TheCdx2polymorphism is located in the promoter region of theVDRgene,and its minor allele T (or A on its complementary strand) resulted in a higher transcriptional activity ofVDR[63]. But these mechanisms remain to be studied for the polymorphisms ofApa1located in the 8thintron andTaq1in the 9thexon,causing a nonsense variation.

LIMITATIONS AND STRENGTHS

This study is the first to propose that different genotypes ofVDRSNPs affect serum lipid profiles with gender differences.This can be useful for improving preventive strategies for dyslipidemia and related diseases in populations with a specific gender and genetic makeup.However,attention should be paid to our findings considering at least the following aspects:1) the participants were recruited from a health examination center where the dyslipidemia prevalence was lower than that indicated in an epidemiological study of the entire Shenzhen population (27%vs.34.6%),and larger population and multi-center studies are needed to further confirm the relationship betweenVDRSNPs and different lipid parameters;2) identification and quantification of confounding factors,especially sex,age,and lifestyle factors,need more comprehensive consideration in interpreting the relationship between dyslipidemia andVDRpolymorphisms,as well as serum VD levels.

CONCLUSIONS

Cdx2,Fok1polymorphisms ofVDRwere sexdifferentially associated with lipid profiles in Han Chinese adults.In women,Fok1B was positively correlated with LDLC,indicating thatFok1B is a risk factor of high-LDLC-dyslipidemia,while in men,Cdx2T andFok1B were negatively correlated with HDLC,suggesting thatCdx2T andFok1B are risk factors of low-HDLC-dyslipidemia.In men,the VD status interacted withCdx2TT andFok1BB in association with low-HDLC levels.

CONFLICTS OF INTEREST

None of the authors have any conflicts of interest to declare.