The role of vitamin D receptor gene polymorphisms in gestational diabetes mellitus susceptibility: a meta-analysis

Abstract

Background

Gestational diabetes mellitus (GDM) is a common disease during pregnancy. The association of vitamin D receptor (VDR) polymorphisms with GDM is still controversial. This study aimed to assess the associations between VDR polymorphisms and GDM risk.

Methods

We searched Cochrane Library, PubMed, and Embase electronic database for all eligible studies published from Jan 1, 1980 to December 31, 2020 to conduct a Meta-analysis. We analyzed four VDR polymorphisms: BsmI (rs1544410), ApaI (rs7975232), TaqI (rs731236), and FokI (rs2228570). Inclusion Criteria: (1) The data can be evaluated; (2) case–control study; and (3) meeting the Hardy–Weinberg’s law. Exclusion criteria: (1) Insufficient or extractable data; (2) Severe publication bias in the data; and (3) duplicate publications. We eventually included 15 studies in seven articles, including 2207 cases and 2706 controls.

Results

We eventually included 15 studies in seven articles, including 2207 cases and 2706 controls. The data showed that ApaI (rs7975232) VDR gene polymorphism was related with the risk of GDM for the comparison of CC vs AA and recessive model in overall population and FokI (rs2228570) VDR gene polymorphism was associated with the risk of GDM for recessive model in overall population. BsmI (rs1544410) polymorphism was not related with the risk of GDM in overall population. However, in the analysis of subgroups grouped by race, BsmI (rs1544410) has certain correlations. And, the data suggested the TaqI (rs731236) polymorphism was not associated with GDM.

Conclusion

Based on the meta-analysis, VDR ApaI (rs7975232) and FokI (rs2228570) polymorphisms increase susceptibility to GDM. In the future, it can be used to diagnose and screen molecular biomarkers for GDM patients.

Background

Gestational diabetes mellitus (GDM) is defined as glucose intolerance diagnosed during pregnancy [1]. GDM is characterized by increased insulin resistance, hyperglycemia, and obesity [2,3,4]. The prevalence of GDM is increasing in decades and floating from 1.7 to 11.6% among populations [5]. Although considerable research effort has been focused on GDM, the pathophysiology of the disease remains incompletely understood. Genetic and environmental factors play an important role in the etiology of GDM [2].

Vitamin D deficiency is associated with diabetes mellitus [6,7,8]. Vitamin D receptor (VDR) gene polymorphisms may contribute to development of diabetes mellitus through calcium metabolism alteration and modulation of insulin secretion [9,10,11]. Three single nucleotide polymorphisms BsmI, ApaI and TaqI of the VDR gene were found in the major untranslated regions that regulate gene expression. FokI is a T > C substitution that results in exon 2 [12, 13]. The above four VDR gene polymorphisms all have a certain effect on insulin production, and secretion plays a role in the pathogenesis of GDM. Therefore, VDR gene polymorphisms may plays a role in the pathogenesis of GDM.

Many studies have researched the role of VDR gene polymorphisms in GDM. It is reported that VDR has four well-characterized di-allelic polymorphisms: BsmI (A > G, rs1544410), ApaI (A > C, rs7975232), TaqI (T > C, rs731236), and FokI (C > T, rs2228570). However, the results of these studies are still uncertain [13,14,15,16,17,18,19]. Different research teams and research designs might lead to differences in results. The objective is to clarify the effect of VDR gene polymorphisms on GDM risk, we conducted a meta-analysis of all eligible case–control studies.

Methods

Search strategy

We identified the keywords “VDR” OR “vitamin D receptor” AND “polymorphism” OR “variant” OR “allele” OR “genotype” OR “gestational diabetes” OR “gestational diabetes mellitus” OR “GDM” to search the articles in Cochrane Library, PubMed, and Embase electronic database. All articles published until December 31, 2020. In addition, manually search the article's reference list for more literature. This article does not collect unpublished data. When multiple articles contain studies of the same population, complete studies were chosen in this study. The language of the publication is limited to English or at least an English abstract.

Inclusion and exclusion criteria

Inclusion Criteria: (1) The data can be evaluated; (2) case–control study; and (3) meeting the Hardy–Weinberg’s law. Exclusion criteria: (1) Insufficient or extractable data; (2) Severe publication bias in the data; and (3) duplicate publications.

Data extraction

The data was independently evaluated by two reviewers according to include and exclude criteria for these documents, discuss whether can be included in the meta-analysis. The difference was not resolved until the consensus of each item was reached. The following information was recorded for each study: author’s name, year of publication, country of origin, racial descent, source of the control population, genotyping methods, matched factors as well as adjusted factors, number of cases and controls.

Statistical analysis

ORs (odds ratios) and 95% CIs were used to estimate the relationships between VDR gene polymorphism and GDM. For heterogeneity detection, we chose the P value to measure. If P < 0.05, we chose the random effect model, otherwise chose the fixed effect model. For publication bias we calculated Egger and Begg’ test, respectively (P < 0.05 was considered representative of statistically significant publication bias). If P < 0.05, it was considered biased. Hardy–Weinberg’s law was detected in all control groups. This meta-analysis was performed using STATA (version 14.0; US).

Results

Study selection

We found 186 records through a full search of the database. After several rounds of screening, 36 articles met our requirements. After two individuals independently evaluated the inclusion and exclusion criteria, 15 case–control studies in a total of seven articles were included in the study [13,14,15,16,17,18,19]. We identified 186 articles from the database, and after excluding irrelevant and duplicate research, 36 articles entered the next step of analysis. According to the inclusion and exclusion criteria, seven articles were included in our study. The specific retrieval process was shown in Fig. 1.

Fig. 1

Flow diagram detailing procedures of selecting eligible studies

Study characteristics

We identified 15 independent studies in seven eligible reports, including 2207 cases and 2706 controls. The main characteristics of all the studies included in our study were shown in Table 1. There were 5 case–control studies on BsmI (rs1544410) [14,15,16,17, 19], 4 case–control studies on TaqI (rs731236) [13,14,15,16], 3 case–control studies on FokI (rs2228570) [14, 15, 17] and 3 case–control studies on ApaI (rs7975232) [13,14,15]. 15 independent studies consisted of 4 Asian [16, 19], 3 African [17] and 8 Caucasian populations [13,14,15, 18].

Table 1 Basic information of the original articles included in this meta-analysis

Publication bias

Funnel plot for comparison of allele models for ApaI (Fig. 2A), FokI (Fig. 2B) and BsmI (Fig. 2C) gene polymorphisms was evaluated to intuitively show the situation of publication bias. We used Begg’s test and Egger’s test to assess publication bias (Table 2). The results of the Egger’s test are P = 0.03 for the contrast of CT vs TT + CC of FokI (rs2228570), while the Begg’s test are P = 0.296. Publication bias was not observed in any other analysis under various other comparative models.

Fig. 2

The funnel plot compared with the allele model for a ApaI (C vs A), b FokI (T vs C) and c BsmI (G vs A) gene polymorphisms to show publication bias

Table 2 Summary ORs (95% CI) of VDR gene polymorphisms and gestational diabetes mellitus risk

ApaI (rs7975232)

The results showed that in the total population of CC vs AA and the recessive model, ApaI (rs7975232) was associated with a higher GDM risk (CC vs AA: OR = 1.974, 95% CI 1.276–3.054, P = 0.002, Fig. 3; CC vs AA + CA: OR = 1.548, 95% CI 1.080–2.217, P = 0.017, Fig. 4). In the subgroup analysis, compared with the CC vs AA and recessive models in the Asian population, it was found to be associated with a higher risk of GDM (CC vs AA: OR = 2.996, 95% CI 1.278–7.022, P = 0.012, Fig. 3; CC vs AA + CA: OR = 3.114, 95% CI 1.403–6.912, P = 0.005, Fig. 4), and CC vs AA comparison among Caucasian populations. (CC vs AA: OR = 1.679, 95% CI 1.006–2.804, P = 0.048, Fig. 3). Table 2 shows other related results of ApaI (rs7975232).

Fig. 3

Fixed-effects meta-analysis on GDM risk and VDR ApaI (rs7975232) polymorphism in overall, Asian and Caucasian population (CC versus AA)

Fig. 4

Fixed-effects meta-analysis on GDM risk and VDR ApaI (rs7975232) polymorphism under recessive model in overall and Asian population (CC vs AA + CA)

FokI (rs2228570)

The results showed that in the recessive model, FokI (rs2228570) was associated with a higher GDM risk in overall population (TT vs CC + CT: OR = 1.454, 95% CI 1.037–2.040, P = 0.030, Fig. 5). In the subgroup, a relationship with a higher GDM risk was found in the Caucasian population under the allele and recessive models (T vs C: OR = 1.631, 95% CI 1.142–2.329, P = 0.007, Fig. 6; TT vs CC + CT: OR = 1.988, 95% CI 1.235–3.200, P = 0.005, Fig. 5). The other related results of FokI (rs2228570) were shown in Table 2.

Fig. 5

Fixed-effects meta-analysis on GDM risk and VDR FokI (rs2228570) polymorphism under recessive model in overall and Caucasian population (TT vs CC + CT)

Fig. 6

Random-effects meta-analysis on GDM risk and VDR FokI (rs2228570) polymorphism under allelic model in Caucasian population (T vs C)

BsmI (rs1544410)

The results showed that BsmI (rs1544410) was not related to GDM risk in the general population. In the subgroup, a relationship with higher GDM risk was found in the Asian population allele model, the comparison of GA vs AA, the dominant model and the over-dominant model. (G vs A: OR = 1.959, 95% CI 1.272–3.017, P = 0.002; GA vs AA: OR = 2.059, 95% CI 1.317–3.217, P = 0.002; GA + GG vs AA: OR = 2.059 95% CI 1.317–3.217, P = 0.002, Fig. 7; GA versus AA + GG: OR = 2.059, 95% CI 1.317–3.217, P = 0.002, Fig. 8). In the subgroup, relationships with lower GDM risk were found in African populations through allele models, GG vs AA, GA vs AA, dominant, recessive and over-dominant models (G vs A: OR = 0.301, 95% CI 0.213–0.427, P < 0.001; GG vs AA: OR = 0.109, 95% CI 0.051–0.232, P < 0.001; GA vs AA: OR = 0.234, 95% CI 0.137–0.401, P < 0.001; GA + GG vs AA: OR = 0.188, 95% CI 0.113–0.312, P < 0.001, Fig. 7; GG vs AA + GA: OR = 0.251, 95% CI 0.126–0.498, P < 0.001; GA vs AA + GG: OR = 0.526, 95% CI 0.329–0.840, P = 0.007, Fig. 8). Other related results of BsmI (rs1544410) were shown in Table 2.

Fig. 7

Random-effects meta-analysis on GDM risk and VDR BsmI (rs1544410) polymorphism under dominant model in Asian and African population (GA + GG vs AA)

Fig. 8

Random-effects meta-analysis on GDM risk and VDR BsmI (rs1544410) polymorphism under over-dominant model in Asian and African population (GA vs AA + GG)

TaqI (rs731236)

The data showed that the TaqI (rs731236) polymorphism of the VDR gene was not related to susceptibility to GDM (Table 2). TaqI (rs731236) was heterogeneous in CT and TT contrast, overt dominant models, and overdominant models in overall population. In the subgroup, CT versus TT showed heterogeneity between the dominant model and the over-dominant model (Table 2).

Sensitivity analyses

One-way sensitivity analysis was performed on the data involved in this meta-analysis. Each study of the meta-analysis was deleted to reflect the overall impact of each data set, and the corresponding combined results did not change substantially.

Discussion

GDM has become major health concern worldwide. Studies suggested that VDR gene polymorphisms might have an impact on GDM risk [14, 16, 18]. However, it is difficult to obtain more accurate results through a single study to determine the relationship between genes and diseases. Meta-analysis can solve the problem of insufficient statistics in a single study, so as to draw more precise conclusions. The association of VDR gene polymorphisms with the incidence of cancer, osteoporosis, and autoimmune thyroid disease has been confirmed in a meta-analysis [20,21,22]. In our study the PICO was shown as follow: P: Gestational diabetes mellitus; I: vitamin D receptor (VDR) polymorphisms; C: control people; O: susceptibility. This study showed ApaI (rs7975232) VDR gene polymorphism was related with GDM for the comparison of CC vs AA and recessive model in overall population and FokI (rs2228570) VDR gene polymorphism was associated with the risk of GDM for recessive model in overall population. The BsmI (rs1544410) and TaqI (rs731236) polymorphisms of the VDR were not related with GDM in overall population.

Due to differences between races, evidence that could cause disease is sometimes not very reliable. This suggests that different races influence genetic background differently [23]. Therefore, based on subgroup analysis of different races, it can be found that the same polymorphisms in disease susceptibility in different populations play different roles. In our study, subgroup analysis suggested that the VDR gene ApaI (rs7975232) polymorphism was significantly associated with GDM for the comparison of CC vs AA and recessive model in Asian population and under the comparison of CC vs AA in Caucasian population. For VDR gene FokI (rs2228570) polymorphism, it was significantly associated with GDM under the comparison of CC vs AA and the recessive model in Asian and under allelic model and the recessive model in Caucasian. However, for VDR gene BsmI (rs1544410) polymorphism, it was significantly associated with GDM under allelic model, the comparison of GA vs AA, dominant model, and over-dominant model in Asian and under allelic model, the comparison of GG vs AA, the comparison of GA vs AA, dominant model, recessive model and over-dominant model in African population. Interestingly, the subgroup analysis in Asia and Africa for BsmI (rs1544410) is the opposite, perhaps because of ethnic differences. Of course, it also may be the difference in results caused by the insufficient number of studies included. We certainly need more and better research to get more reliable results.

However, this meta-analysis has some limitations. Firstly, heterogeneity may influence the results of this meta-analysis. Nonetheless, we use specific research standards to strictly perform data extraction and analysis to minimize this possibility. Secondly, the study only includes published studies, and the existence of results indicating no meaning or negative results may not be published, and this will increase the likelihood of publication bias. Finally, our results have not been adjusted. If you can get more research data, you should be able to analyze it more accurately. We can obtain more accurate results by adjusting other variables, including age and family history, etc. [24,25,26,27]. In addition, an in-depth analysis of these factors provides a more complete understanding of the linkages between these factors and the risks of GDM.

Conclusions

In summary, VDR ApaI (rs7975232) and FokI (rs2228570) polymorphisms increase susceptibility to GDM. In the future, it can be used to diagnose and screen molecular biomarkers for GDM patients. VDR BsmI (rs1544410) polymorphism was associated with GDM in Asian and African population. VDR TaqI (rs731236) polymorphism was not associated with GDM.