Association of VDR Polymorphisms and Gene–obesity lnteraction with Type 2 Diabetes: A Case–control Study among Chinese Rural Population＊

XU Ze, ZHANG Dong Dong, LIU Ya Ping, ZHANG Yu Jing, XUE Yuan, GAO Jiao Jiao,ZHAO Tong, HUANG Hao Yue, LI Wen Jie, and LI Xing

In recent years, the prevalence of type 2 diabetes(T2D) has increased at an alarming rate, and the epidemic of T2D mellitus and its complications poses a major global health threat. Globally, the International Diabetes Federation estimated that 10.5% of adults aged 20–79 years (536 million adults) had diabetes mellitus in 2021. If these trends continue, approximately 783 million people will have diabetes by 2045[1]. T2D is a complex metabolic disease caused by a complex interplay of genetic and environmental risk factors. Obesity, a sedentary lifestyle, and energy-dense diets are the best-known environmental factors associated with T2D. The prevalence of overweight and obesity is high in Chinese, and the number of people with obesity accounts for approximately 46% of the adult population.

Accumulating evidence has demonstrated that low circulating vitamin D levels are associated with an increased risk of developing T2D across diverse populations. The vitamin D receptor (VDR) is a member of nuclear receptor superfamily 1, and it is expressed in many cells and tissues, such as pancreatic β-cells, liver, and adipose tissue.VDRmediates the biological effects of vitamin D, and 1,25-dihydroxyvitamin D depends onVDRto exert its biological effects. Most studies on the association betweenVDRpolymorphism and T2D have focused on six single nucleotide polymorphisms (SNPs): Fok1,Bsm1, Taq1, Apa1, Tru91, and 3’UTR[2]. However, the association ofVDR(rs2189480 and rs739837)variants with T2D remains unclear, and the interaction between the loci and obesity on T2D susceptibility has not been explored. Therefore, this case–control study aimed to investigate the relationship betweenVDRSNPs and T2D and further explore the combined effects of obesity andVDRgene polymorphisms on T2D, which may help in predicting the risk of T2D and identifying individuals at high risk. Furthermore, the study may provide novel strategies for T2D prevention and individualized anti-diabetes treatment.

From July to August 2013 and from July to August 2015, we recruited 2,378 Han Chinese from Wuzhi, Xin’an, and Houzhai County in Henan Province. All participants completed a standard questionnaire (general characteristics, lifestyle,medical history, and family history of disease), body measurements, and fasting blood collection. T2D was defined as fasting serum glucose ≥ 7.0 mmol/L and/or taking hypoglycemic medication and/or having a self-reported T2D history[3]. Individuals with a self-reported history of gestational diabetes or other types of diabetes, coronary heart disease,chronic kidney diseases, cancer, or failed DNA extraction were excluded from the study. We performed a 1:2 ratio case–control match by sex and age (within 3 years); finally, a total of 397 T2D and 794 non-T2D cases were included in our study.

All participants provided informed consent before the examination, and the study was approved by the Medical Ethics Committee of Zhengzhou University.

For each participant, height, weight, waist circumference (WC), and hip circumference were measured twice to the nearest 0.5 kg and 0.1 cm,and the body mass index (BMI) was calculated in kg/m2. All participants were asked to assemble on time after an overnight fast (＞ 10 h) for blood sample collection. Oral antidiabetic drugs, insulin, and alcohol should have been avoided for at least 12 h before blood collection. Blood samples were used for the DNA isolation and measurements of glucose,lipid profile, and serum 25-hydroxyvitamin D[25(OH)D] concentrations. 25(OH)D levels were measured by human enzyme-linked immunosorbent assay kits (Sangon Biotech, Shanghai, China). The levels of fasting plasma glucose (FPG), triglyceride(TG), total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and fasting insulin (FINS) were measured by an automatic biochemical analyzer(KHB, Shanghai, China). The Homeostasis Model Assessment Insulin Resistance (HOMA-IR) was applied as a measure of the degree of insulin resistance [HOMA-IR = FPG (mmol/L) × FINS(mIU/L)/22.5]. Genomic DNA was isolated from the peripheral blood following the standard procedures of the DNA extraction kit (Bioteke, Beijing, China).The TaqMan probe method and real-time polymerase chain reaction technology were used to identify the rs2189480 and rs739837 genotypes.

Non-normally distributed variables were shown as median values with interquartile range, and categorical variables were described with frequencies and percentages. The distribution of participants’ characteristics between the case and control groups was compared using the nonparametric Wilcoxon rank-sum test or Chi-square test. Multivariate logistic regression was used to examine the relationship betweenVDRSNPs and T2D, and odds ratio (OR) and corresponding 95%confidence interval (CI) were calculated after adjustment for smoking, alcohol consumption,physical activity, 25(OH)D, and diabetes family history, assuming codominant, recessive, and dominant models of inheritance. Multivariable linear regression was performed to explore the association of BMI and WC with blood glucose-related parameters, including FPG, hemoglobin A1c, HOMAIR, and HOMA-β. Multiple logistic regression analysis was used to examine the interaction betweenVDRSNPs and obesity on T2D risk according to the genotype of each SNP (under the recessive or dominant genotype model) and obesity exposure status. Multiplicative interactions were evaluated using a multiplicative model based on generalized linear regression, with FPG as the dependent variable. Covariance analysis was used to explore the influence ofVDRvariations on cardiovascular metabolic indicators. Unless otherwise stated,statistical analyses were performed using SAS v.9.1(SAS Institute, Cary, NC, USA), and the significance level was set atP＜ 0.05.

Results of the comparisons of the basic characteristics of the two groups are presented in Supplementary Table S1 (available in www.besjournal.com). A total of 1,191 individuals (397 T2D cases/794 controls) were included in this case–control study, and male and female participants accounted for 39.3% and 61.6%,respectively. Compared with the control group,patients with T2D had higher levels of FPG, FINS,HOMA-IR, TC, TG, and LDL-C but lower levels of HOMA-β and HDL-C (P＜ 0.05). As expected, patients with T2D had significantly higher BMI and WC than controls, respectively (P＜ 0.05).

At present, the results of studies on the association betweenVDRpolymorphism and chronic diseases are still contradictory. Our study explored the association between two loci (rs2189480 and rs739837) and T2D risk through multiple logistic regression analysis (Table 1). The genotypic distributions of the two SNPs did not deviate from the Hardy–Weinberg equilibrium (P＞ 0.05). Multiple logistic regression analysis showed that participants with the “AA” genotype of rs2189480 were significantly less likely to develop T2D than those with the “GA+GG” genotype (P= 0.036,OR= 1.30,95%CI1.02–1.66). The rs2189480 locus is one of theVDRSNPs located in the intron and consists of a G-A change. Studies have shown that rs2189480 affects regulatory T cell function and T cells help regulate inflammatory activity[4]. Therefore, the rs2189480 mutation may have a protective effect on T2D by regulating the inflammatory response. Studies have reported the effect of rs739837 variants on diabetes,but no unanimous conclusions have been drawn,and the role of the rs739837 locus is still unclear.Our study results showed that carrying “CC+GC” of rs739837 rather than “GG” was significantly associated with an increased hazard of developing T2D (P= 0.028,OR= 1.32, 95%CI1.03–1.70). The rs739837 is located in the 3’UTR of theVDRgene,which contains many regulatory elements that influence its post-transcriptional regulation. A study showed that variations in rs739837 may influence the stability ofVDRmRNA and its gene transcription and further increase the risk of diabetes[5].

Obesity is believed to be a promoter of T2D.Obesity is strongly associated with the development of insulin resistance, which in turn plays a key role in the pathogenesis of obesity-associated cardiometabolic complications, including metabolic syndrome components, T2D, and cardiovascular diseases. Our study investigated the association between BMI/WC and T2D in the rural Chinese population(Supplementary Table S2, available in www.besjournal.com) and the influence of BMI/WC and blood glucose-related parameters (Supplementary Table S3, available in www.besjournal.com). The results showed that higher BMI or WC was significantly associated with an increased risk of T2D in Chinese rural populations. Moreover, after multivariable adjustment, we found that for each increase in BMI of 1 kg/m2, the level of FINS and HOMA-IR increased by 0.49 mIU/L and 0.19 (P＜ 0.05),respectively. For each increase in WC of 1 cm, the level of FPG, FINS, and HOMA-IR increased by 0.017 mmHg,0.49 mIU/L, and 0.19 (P＜ 0.05), respectively.

To assess the interaction between obesity status and SNPs of theVDRgene on T2D, the participants were divided into four subgroups by genetic models of SNPs and obesity status (Table 2). An interactionbetweenVDRvariants and obesity on T2D risk has been found (Pinteraction= 0.029). Compared with nonobese individuals carrying rs739837 “GG” genotype,those with overweight/obesity carrying the “CC” or“GC” genotype had a higher risk of developing T2D(OR= 3.06, 95%CI2.04–4.61). Further, the results of the generalized linear model showed a significant interaction between rs739837 and BMI (P= 0.041).At higher BMI (＞ 26.96 kg/m2), rs739837 “C” allele carriers had higher FPG levels and showed a stronger effect than individuals with the “GG” genotype(Figure 1).VDRis expressed in pancreatic β cells, and alleles ofVDRregulate β-cell function. Human and animal studies have shown that disruption of theVDRsignaling pathway can affect insulin secretion and decrease insulin mRNA levels[6,7]. The ability of fat cells to store fat is gradually reduced with obesity, and fat is stored in places such as visceral tissue, which adversely affects insulin signaling. TheVDRgene mutation and obesity may have an interactive effect on T2D risk by affecting insulin secretion. Meanwhile, our results found that subjects carrying "GC" or "CC "of rs739837 had higher fasting insulin levels than "GG" carriers,suggesting that mutations in theVDRgene may affect the insulin signaling pathway and lead to insulin resistance in the subjects (Supplementary Table S4, available in www.besjournal.com). On the contrary, studies have reported thatVDRis expressed in adipose tissue and dynamically upregulated during adipocyte differentiation[8]. The T allele of TaqI loci was associated with a higher risk of obesity in Greek individuals (contributing to an elevated BMI of 3 kg/m2per risk allele)[9].VDRplays a key role in fat formation by affecting vitamin D levels.Meanwhile, the expression ofVDRmRNA also changes during adipocyte differentiation[10].Therefore, we speculate that the interaction betweenVDRgene variation and adipocyte differentiation and metabolism may contribute to the risk of T2D. However, the specific mechanism of interaction needs further study.

Table 1. Genotypic distribution of VDR polymorphisms and the association of VDR polymorphisms with T2D risk

Table 2. Interaction analysis of VDR polymorphisms and obesity on T2D

Figure 1. Multiplicative interactions between rs739837 (G ＞ C) and BMI (A) or WC (B) and the associations of these interactions with FPG. The blue trend line represents the correlations between BMI or WC and FPG levels in individuals with risk alleles. The red trend line represents the correlation between BMI or WC and FPG levels in individuals without risk alleles. The point where the two lines cross indicates that for individuals whose BMI or WC is higher than the level value at this point, risk allele carriers have a greater effect on systolic or diastolic blood pressure than non-risk allele carriers. BMI, body mass index;FPG, fasting plasma glucose; WC, waist circumference.

This study explores the interaction betweenVDRgene (rs2189480 and rs739837) variation and obesity on T2D susceptibility in the Chinese rural population. Multiplicative interactions between rs739837 variants and BMI or WC were evaluated using a multiplicative model based on generalized linear regression, with FPG as the dependent variable. This study also employed a 1:2 case–control design, including serum samples from 1191 subjects covering a wide age span, which may be suitable for the exploration of T2D etiology.However, this study has several limitations. First,this was a single study limited to Chinese Han populations, which may limit the generalizability of the findings to other ethnic groups. Second, only one gene was discussed in the study; other regulatory genes may contribute to the pathogenesis of T2D. Furthermore, the study was limited by the lack of data on more informative measures of adiposity in addition to BMI and WC,such as visceral fat index or percentage of body fat.Therefore, further studies are needed to verify the observed interactions in other populations and further investigate the mechanisms.

In conclusion, our study found that the “CC+GC”genotypes of rs739837 in theVDRgene were associated with an increased risk for T2D development. Interactions betweenVDRand obesity were associated with the susceptibility for T2D. The observations from this study constitute evidence that people who carry particular variants of theVDRrs739837, which improved BMI level control, may be beneficial in inhibiting T2D progression.

#Correspondence should be addressed to LI Xing, Email: lixing530@zzu.edu.cn, Tel: 86-371-67781305.

Biographical note of the first author: XU Ze, female,born in 1995, Master, majoring in nutrition and food hygiene.

Received: February 10, 2022;Accepted: September 20, 2022