Skip to main content
  • Research article
  • Open access
  • Published:

Vitamin D receptor gene polymorphisms and the risk of the type 1 diabetes: a meta-regression and updated meta-analysis

Abstract

Background

The association between the polymorphisms in the vitamin D receptor (VDR) gene and the risk of type 1 diabetes mellitus (T1DM) has been evaluated in several studies. However, the findings were inconclusive. Thus, we conducted a meta-analysis to comprehensively evaluate the effect of VDR gene polymorphisms on the risk of T1DM.

Methods

All relevant studies reporting the association between VDR gene polymorphisms and susceptibility to T1DM published up to May 2020 were identified by comprehensive systematic database search in ISI Web of Science, Scopus, and PubMed/MEDLINE. Strength of association were assessed by calculating of pooled odds ratios (ORs) and 95% confidence intervals (CIs). The methodological quality of each study was assessed according to the Newcastle–Ottawa Scale. To find the potential sources of heterogeneity, meta-regression and subgroup analysis were also performed.

Results

A total of 39 case–control studies were included in this meta-analysis. The results of overall population rejected any significant association between VDR gene polymorphisms and T1DM risk. However, the pooled results of subgroup analysis revealed significant negative and positive associations between FokI and BsmI polymorphisms and T1DM in Africans and Americans, respectively.

Conclusions

This meta-analysis suggested a significant association between VDR gene polymorphism and T1DM susceptibility in ethnic-specific analysis.

Peer Review reports

Background

Type 1 diabetes mellitus (T1DM) is a globally-widespread disease that is characterized by a reduction in insulin production or the production of ineffective insulin [1]. It is generally believed that the immune-associated destruction of beta cells of the islets of Langerhans causes the disease, resulting in lower insulin levels (that is called type 1a diabetes mellitus). In a smaller T1DM subset, no evidence of autoimmunity can be found (type 1b) [2]. T1DM constitutes roughly 5 to 10% of all diabetes cases, and its prevalence is still rising [3]. With more than half a million children living with T1DM, and almost 90,000 children diagnosed each year, T1DM inflicts mostly children of under 15 years of age [4]. It is well known that T1DM is a multi-factorial autoimmune disorder caused by interactions between genetic and environmental factors [5].

Vitamin D (VitD) is a steroid molecule that has many roles in the body, such as regulation of the immune cells. In addition to immune responses, VitD is also involved in the etiopathogenesis of several disorders, such as cancer, autoimmune disorders, cardiovascular disorders, asthma, and diabetes [6,7,8,9]. In animal model of T1DM, VitD suppresses the occurrence of diabetes, by regulating the T helper (Th) 1/Th2 cytokine balance in the local pancreatic lesions [10, 11]. Moreover, VitD inhibits T cell activation and secretion of pro-inflammatory cytokines, such as interleukin (IL)-1, IL-6, IL-12, tumor necrosis factor (TNF)-α, and interferon (IFN)-γ, which are involved in the pathogenesis of T1DM [12,13,14]. Mostly, VitD exerts its function through vitamin D receptor (VDR), which is found in the nuclei of target cells, such as lymphocytes, macrophages, and pancreatic cells. VDR is a member of the nuclear hormone receptors superfamily and has been linked to insulin sensitivity and secretion [15].

Four common single nucleotide polymorphisms (SNPs) of VDR gene are FokI (rs2228570), TaqI (rs731236), BsmI (rs1544410), and ApaI (rs7975232). Among them, ApaI, BsmI, and TaqI polymorphisms are located in the 3′-end of VDR gene which lead to silent mutation associated with increased VDR mRNA stability. In contrast, FokI SNP is located in the start codon that produces a protein with shorter size (424 amino acids), which is more active than the long form (427 amino acids) [8, 16, 17]. Over the course of past few decades, the VDR gene polymorphisms have been associated with susceptibility to numerous autoimmune disorders [8, 18, 19].

In recent years, several studies have investigated the association between VDR gene SNPs and T1DM in all over the world, which have yielded conflicting results. The reasons for these discrepancies might be small sample sizes, clinical heterogeneity, and low statistical power. Therefore, a comprehensive meta-analysis might be the best way to solve these problems. Two previous meta-analyses performed by Tizaouia et al. in 2014 [20] and Guo et al. in 2006 [21] reported that VDR gene polymorphisms were not associated with the susceptibility to T1DM. However, Zhang et al. in 2012 [22] demonstrated that BsmI polymorphism was significantly associated with the risk of T1DM. Furthermore, Sahin et al. in 2017 indicated that BsmI and TaqI polymorphisms were associated withT1DM risk in children with less than average 15 years old [23]. Qin et al. in 2014 evaluated the association of only BsmI SNP with T1DM risk and demonstrated its association in the overall analysis, as well as in Asians, Latino, and Africans [24]. In 2014, Wang et al., by including 20 studies, reported that BsmI polymorphism might be a risk factor for susceptibility to T1DM in the East Asian population, and the FokI polymorphism was associated with an increased risk of T1DM in the West Asian population [25].

Since several articles published after the last meta-analysis, here we conducted an updated meta-analysis with the aim of providing a much more reliable conclusion on the significance of the association between VDR gene polymorphisms and T1DM risk.

Methods

This meta-analysis was conducted according to the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, including search strategy, inclusion and exclusion criteria, data extraction and quality assessment, and statistical analysis [26].

Search strategy

Three electronic databases (PubMed/MEDLINE, Scopus, and Web of Science) were systematically searched for studies regarding the association of VDR gene polymorphisms, including FokI (rs2228570) and/or TaqI (rs731236) and/or BsmI (rs1544410) and/or ApaI (rs7975232), and T1DM susceptibility, which were published before May 2020. The following combinations of search terms were used: (“T1D” OR “type 1 diabetes” OR “diabetes”) AND (“VDR” OR “vitamin D receptor”) AND (“polymorphisms” OR “SNP” OR “variation” OR “mutation”). The reference lists of review articles were also manually searched for additional pertinent publications. Original data in English language and human population studies were collected.

Inclusion and exclusion criteria

Eligible studies must meet the following criteria: a) All studies assessing the association of VDR gene polymorphisms and T1DM risk; b) All studies reporting sufficient data to calculate the odds ratio (OR) and its 95% confidence intervals (CIs); c) All studies with distinct case and control groups (case-control and cohort design). The exclusion criteria were: a) studies that their genotype or allele frequency could not be extracted; b) letters, non-English publications, animal studies, case reports, reviews, comments, book chapters, and abstracts; c) duplicate and republished studies. The application of these criteria recognized 39 studies eligible for the quantitative analysis.

Data extraction and quality assessment

According to a standardized extraction form, the following data were independently extracted by two reviewers: the author’s name, journal and year of publication, country of origin, ethnicity, number of case and control for each gender separately, genotype and allele frequencies in cases and healthy groups, mean or range of age, genotyping method, total sample size of cases and controls. The third reviewer finalized the extracted data, and potential discrepancies were resolved by consensus. For quality assessment of the included publications, the Newcastle-Ottawa Scale (NOS) was applied [27]. In this respect, studies with 0–3, 4–6 or 7–9 scores were, respectively, of low, moderate, and high-quality.

Statistical analysis

Deviation from Hardy–Weinberg equilibrium (HWE) for distribution of the allele frequencies was analyzed by χ2-test in control groups. The strength of association between VDR gene polymorphisms and T1DM susceptibility was estimated by calculating pooled OR and its 95% CI. Different comparison model for FokI, TaqI, BsmI, and ApaI were as follows: FokI; dominant model (ff + Ff vs. FF), recessive model (ff vs. Ff + FF), allelic model (f vs. F), homozygote (ff vs. FF), and heterozygote (Ff vs. FF): TaqI; dominant model (tt + Tt vs. TT), recessive model (tt vs. Tt + TT), allelic model (t vs. T), homozygote (tt vs. TT), and heterozygote (Tt vs. TT): BsmI; dominant model (bb + Bb vs. BB), recessive model (bb vs. Bb + BB), allelic model (b vs. B), homozygote (bb vs. BB), and heterozygote (Bb vs. BB): ApaI; dominant model (aa+Aa vs. AA), recessive model (aa vs. Aa+AA), allelic model (a vs. A), homozygote (aa vs. AA), and heterozygote (Aa vs. AA). The heterogeneity among studies was measured by the χ2 test-based Q statistic, and I2 value which quantify the degree of heterogeneity [28]. Accordingly, heterogeneity was considered significant if I2 values exceeded 50% or the Q statistic had a P value of less than 0.1 and random-effects model (DerSimonian–Laird approach) was carried out [29]. Otherwise, the fixed-effects model (Mantel–Haenszel approach) was performed for combination of data [30]. In order to assess the predefined sources of heterogeneity among included studies, subgroup analysis and meta-regression analysis based on year of population, and ethnicity were performed. Stability of our results was assessed by sensitivity analysis. Potential publication bias was estimated by Egger’s linear regression test, and also Begg’s test was employed to estimate the funnel plot asymmetry (P value< 0.05 considered statistically significant) [31, 32]. The data analyses were carried out using STATA (version 14.0; Stata Corporation, College Station, TX) and SPSS (version 23.0; SPSS, Inc. Chicago, IL).

Results

Study characteristics

Regarding to aforementioned keywords, a total of 1116 studies were initially retrieved. Of these studies, 456 publications were duplicate, 559 and 62 publications excluded by title & abstract and full text examination, respectively. Finally, 39 studies qualified for quantitative analysis. It should be noted that while the latest meta-analysis by Tizaouia et al. [20] in 2014 included 23 studies, we performed the updated meta-analysis by adding 16 more articles. Also, no studies were found by hand search (Fig. 1). The eligible studies were published from 1998 to 2019 and had an overall good methodological quality with NOS scores ranging from 6 to 8. Polymerase chain reaction-restriction fragment length polymorphism (PCR- RFLP) and Taq-man were used by majority of included studies as genotyping method. Tables 1 and 2 summarized the characteristics and genotype frequency of the included studies.

Fig. 1
figure 1

Flow diagram of study selection process

Table 1 Characteristics of studies included in meta-analysis of overall T1DM
Table 2 Distribution of genotype and allele among T1DM patients and controls

Quantitative synthesis

Meta-analysis of the association between FokI (rs2228570) polymorphism and T1DM risk

Overall, 29 case-control studies with 3723 cases and 5578 controls were analyzed for assessment of FokI polymorphism and T1DM risk. Of 29 studies, 15 studies were conducted in European countries [15, 34,35,36, 38,39,40,41, 44, 47, 49, 71], 9 studies were in Asian countries [33, 42, 45, 46, 48, 54,55,56,57], 3 studies were in African population [51,52,53] and eventually one study in Australia [50] and one study in American population [43]. Among studies were performed in Europe, Audi et al. [71] conducted an association study in different city of Spain (Barcelona and Navarra) and reported all data separately including genotype and allele frequency; thus we considered each population as a separate study. The pooled results revealed no significant association in overall population across all genotype models, meanwhile subgroup analysis according to ethnicity showed decreased risk of T1DM susceptibility in European population [dominant model (OR = 0.86, 95% CI, 0.74–1.00, P = 0.05) and heterozygote contrast (OR = 0.86, 95% CI, 0.75–0.99, P = 0.04)] and increased risk of T1DM susceptibility in African population under all genotype models; dominant model (OR = 2.06, 95% CI, 1.20–3.53, P = 0.008), recessive model (OR = 2.14, 95% CI, 1.03–4.43, P = 0.04), allelic model (OR = 1.17, 95% CI, 1.06–2.97, P = 0.02), ff vs. FF model (OR = 3.11, 95% CI, 1.44–6.69, P = 0.004), and Ff vs. FF model (OR = 1.81, 95% CI, 1.13–2.91, P = 0.01). Besides, susceptibility to T1DM in Asians compared to Africans and Europeans were not affected by FokI polymorphism (Fig. 2). The results of pooled ORs, heterogeneity tests and publication bias tests in different analysis models are shown in Table 3.

Fig. 2
figure 2

Pooled OR and 95% CI of individual studies and pooled data for the association between ApaI gene polymorphism and T1DM risk in heterozygote contrast (Aa vs. AA)

Table 3 Main results of pooled ORs in meta-analysis of Vitamin D Receptor gene polymorphisms

Meta-analysis of the association between TaqI (rs731236) polymorphism and T1DM risk

There were 20 case-control studies with 1837 cases and 1895 controls concerning TaqI polymorphism and T1DM risk. Studies were performed in different population, 8 studies were in Europeans [15, 34, 35, 38, 41, 44, 59, 60], 8 studies in Asians [42, 46, 48, 57, 58, 62,63,64], 2 studies in Africans [52, 65] and one study each was in Australia [50] and Americans [61]. Meta-analysis rejected any significant association between TaqI SNP and the risk of T1DM susceptibility. Moreover, the results of subgroup analysis by ethnicity were not significant under five genotype models. In subgroup analysis, since there was only one study for the Australians [50], Americans [61], and two studies for Africans [52, 65], these studies were excluded from the analysis. The results of pooled ORs, heterogeneity tests and publication bias tests in different analysis models are shown in Table 3.

Meta-analysis of the association between BsmI (rs1544410) polymorphism and T1DM risk

To examining the association between BsmI polymorphism and T1DM risk, 34 case-control studies with 4826 cases and 7159 controls subjects were included. It was detected that 15 studies with 1938 cases and 4450 controls were performed in European countries [15, 34,35,36, 38, 40, 41, 44, 49, 59, 60, 71] which among these 15 studies, Turpeinen et al. [36] conducted an association study in different city of Finland (Turku, Tampere and Oulu) and reported all data separately, including genotype and allele frequency; thus we considered each population as a separate study. Moreover, 13 studies out of 34 eligible studies were carried out in Asian populations [42, 46, 48, 56,57,58, 62,63,64, 67,68,69,70], 3 studies were in Americans [43, 61, 66] and three studies were in Africans [52, 53, 65]. No significant association between BsmI polymorphism and T1DM risk were found under all genotype models for the overall population. However, pooled results of subgroup analysis indicated markedly significant negative associations between BsmI SNP and the risk of T1DM susceptibility in American populations across all genotype models; dominant model (OR = 0.57, 95% CI, 0.39–0.84, P = 0.004), recessive model (OR = 0.62, 95% CI, 0.41–0.94, P = 0.02), allelic model (OR = 0.66, 95% CI, 0.54–0.81, P < 0.001), bb vs. BB model (OR = 0.52, 95% CI, 0.34–0.80, P = 0.003), except Bb vs. BB model (OR = 0.66, 95% CI, 0.41–1.05, P = 0.08) (Fig. 3). No significant association was detected for European, Asian and African population. The results of pooled ORs, heterogeneity tests and publication bias tests in different analysis models are shown in Table 3.

Fig. 3
figure 3

Pooled odds ratio (OR) and 95% confidence interval of individual studies and pooled data for the association between FokI, BsmI gene polymorphism and T1DM risk in different ethnicity subgroups and overall populations for A; dominant model (FokI), B; Ff vs. FF Model (FokI), and C; Recessive Model (BsmI)

Meta-analysis of the association between ApaI (rs7975232) polymorphism and T1DM risk

Finally, 24 case-control studies with 2436 cases and 4074 controls were identified eligible for quantitative synthesis of the association between ApaI polymorphism and T1DM risk. Overall, 10 studies were conducted in Europe [15, 35, 36, 38, 41, 44, 59, 60], 10 studies were in Asia [42, 46, 48, 54, 55, 57, 58, 62,63,64], 2 studies in Africa [52, 65] and one study each was in Australia [50] and America [61]. Because of limited number of studies performed in Australia, America and Africa these studies were excluded from subgroup analysis. The results demonstrated no significant association between the ApaI polymorphism and risk of T1DM in the overall population and ethnic-specific analysis (Fig. 3). The results of pooled ORs, heterogeneity tests and publication bias tests in different analysis models are shown in Table 3.

Evaluation of heterogeneity and publication bias

During the meta-analysis of VDR gene polymorphism evidence of substantial to moderate heterogeneity was detected. However, partial heterogeneity was resolved while the data were stratified by ethnicity. Publication bias was evaluated by funnel plot, Begg’s test and Egger’s test. There was no obvious evidence of asymmetry from the shapes of the funnel plots (Fig. 4), and all P values of Begg’s test and Egger’s test were > 0.05, which showed no evidences of publication biases.

Fig. 4
figure 4

Begg’s funnel plot for publication bias test. A; dominant model FokI. B; dominant model TaqI. C; dominant model BsmI. D; dominant model ApaI. Each point represents a separate study for the indicated association

Sensitivity analysis

The leave-one-out method was used in the sensitivity analysis to explore the effect of individual data on the pooled ORs. The significance of ORs was not altered through omitting any single study in the dominant model for FokI, TaqI, BsmI and ApaI SNPs, indicating that our results were statistically robust (Fig. 5).

Fig. 5
figure 5

Sensitivity analysis in present meta-analysis investigates the single nucleotide polymorphisms of Vitamin D Receptor contribute to risk for T1DM (A, FokI; B, TaqI; C, BsmI; D, ApaI)

Bayesian meta-regression analysis

Meta-regression and subgroup analyses were performed to explore potential sources of heterogeneity among included studies (Table 4). The findings of meta-regression indicated that ethnicity can be the potential source of heterogeneity, therefore, subgroup analysis was performed to attenuate the effect of these parameters. (Fig. 6).

Table 4 Meta-regression analyses of potential source of heterogeneity
Fig. 6
figure 6

Meta-regression plots of the association between VDR gene polymorphisms and risk of CAD based on; A: Publication year (Dominant model), B: Ethnicity (Recessive model), C: Publication year (Allelic model), C: Ethnicity (aa vs. AA model)

Discussion

In this study, we performed a systematic review and meta-analysis to achieve a vivid and exact approximation of the associations between the VDR gene polymorphisms, including FokI (rs2228570), TaqI (rs731236), BsmI (rs1544410), and ApaI (rs7975232) and susceptibility to T1DM. The findings of meta-analysis on 39 case–control studies, containing 29 studies with 3723 cases and 5578 controls for FokI, 20 studies with 1837 cases and 1895 controls for TaqI, 34 studies with 4826 cases and 7159 controls for BsmI, and 24 studies with 2436 cases and 4074 controls for ApaI, indicated no significant association of VDR gene polymorphisms with T1DM risk in overall population. That notwithstanding, the subgroup analysis resulted in identification of significant associations between FokI and BsmI polymorphism and T1DM in African and American population. Our study provided some beneficial points over previous studies. First, this meta-analysis included further studies with more sample size compared with the previous studies, conferring more conclusive results. Second, we performed subgroup analysis by ethnicity to indicated association of VDR gene polymorphisms with T1DM risk in different ethnical groups.

Over the course of past years, a bulk of studies has addressed the association of VDR gene polymorphisms and risk of T1DM throughout various populations, resulting in conflicting findings [61, 67]. Such discrepancies might stem from diversity in detection methods, differences in diagnostic criterions, clinical heterogeneity, small sample sizes, low statistical power, and interactions between genetic and environmental contributing factors according to variations in the geo-epidemiological factors. As a consequence, three previous meta-analyses by Guo et al. [21] in 2006 [including 11 studies for FokI (1424 cases and 3301 controls), 13 studies for BsmI (1601cases and 4207 controls), 9 studies for ApaI (1101 cases and 2805 controls), and 7 studies for TaqI (681 cases and 781 controls)], Zhang et al. [22] in 2012 [T1DM cases and 4049 controls in 21 studies for BsmI, 2167 T1DM cases and 3402 controls in 17 studies for FokI, 1166 T1DM cases and 2328 controls in 11 studies for ApaI, and 1041 T1DM cases and 1137 controls in 8 studies for TaqI], and Tizaouia et al. [20] in 2014 (13 studies for TaqI, 23 studies for BsmI, 15 studies for ApaI, and 18 studies for FokI) were carried out to resolved the conundrum and attain an exact approximation. They indicated that VDR gene SNPs were not associated with T1DM risk, except than BsmI polymorphism association with T1DM predisposition that was observed in Zhang et al. [22] study. Upon the latest meta-analysis published in 2014, several original association studies evaluated the role of VDR gene polymorphisms with T1DM risk. As a result, the necessity for performing an updated meta-analysis is sensed to come up with resolution of the limitations of individual association studies and to gain a much more valid and comprehensive pooled estimation on the association of VDR gene polymorphisms with T1D risk.

Previous meta-analysis performed by Tizaouia et al. [20] in 2014 reported no significant association of VDR gene FokI polymorphism with risk of T1D. According to our meta-analysis, the pooled results in overall population across all genotype models demonstrated no significant association of VDR gene FokI polymorphism; nonetheless, subgroup analysis according to ethnicity showed a marginally-significant decreased susceptibility to T1DM in European population according to dominant genetic model and heterozygote comparison, while an increased risk of T1DM in African population according to all genotype models. In addition, our meta-analysis did not support any significant association between TaqI SNP and susceptibility to T1DM. Furthermore, the results of subgroup analysis according to ethnicity did not show any significant association in all genetic models. However, in the subgroup analysis, given that there was only one study in the Australian [50] and American [61] populations, and two studies in the African [52, 65] population, the subgroup analysis was not performed in these populations. In line with our findings, previous meta-analysis by Tizaouia et al. [20] also did not show significant association of VDR gene TaqI polymorphism with risk of T1D. According to the previous meta-analysis, BsmI SNP was not the risk factor for T1D susceptibility. However, after excluding one study, a marginal significant (P = 0.051) association was found in the homozygous model. On the other side, our meta-analysis also revealed that BsmI polymorphism was not a risk for T1DM in all genetic models when all of the population were analyzed. Nonetheless, subgroup analysis demonstrated a strong negative significant association between BsmI SNP and the risk of T1DM in American population in all of the genetic model comparisons. Finally neither our meta-analysis nor the previous one by Tizaouia et al. [20] found any significant association of ApaI polymorphism and T1DM risk in overall as well as subgroup analyses. Taken together, although our meta-analysis included further studies compared to the previous study, the overall analysis was almost the same. Nonetheless, our subgroup analysis indicated association of VDR genetic polymorphisms with T1DM risk in different ethnical groups.

In their meta-analysis, Tizaoui et al. [20] indicated in the stratification analysis that publication year, age, gender, estimated VitD levels, and latitude modulated the association between VDR gene polymorphisms and T1D risk. Furthermore, another meta-analysis revealed a relationship between winter ultraviolet radiation (UVR) and VDR gene polymorphisms in T1DM, implying to the influence of the UVR on the association between VDR polymorphisms and T1DM susceptibility [72]. During the four cooler months, it was observed that latitude strongly determines the available levels of VitD producing UV. As latitude increases, the amount of VitD producing UV decreases, which may prevent VitD synthesis in humans [73]. As a result, the latitude of the locations in which the individuals live may impress the susceptibility to develop T1DM.

Despite we tried to conduct best meta-analysis of the VDR gene polymorphisms and susceptibility to RA, there was also a number of limitations that should be taken into account. First, there was significant heterogeneity across studies, which may lessen the certainty of the results. However, we tried to find and attenuate its effect by meta-regression and subgroup analysis. Consequently, heterogeneity was still an unavoidable problem that may influence the accuracy of the overall results. Second, only articles published in the English language were include in this meta-analysis. Third, our meta-analysis was based on crude approximation of the genetic variations regardless of adjusting the analysis by gender, age, VitD intake, and other environmental factors like exposure to sun light, as several studies noted the involvement of these parameters as well as gene-environment and gene-gene interactions in the susceptibility and of RA and we could not analyze it owing to a lack of published well-structured data.

Conclusion

In conclusion, this study was a systematic review and meta-analysis of 40 case–control association studies to come up with the clear estimation of the associations between the VDR gene SNPs [FokI (rs2228570), TaqI (rs731236), BsmI (rs1544410), and ApaI (rs7975232)] and susceptibility to T1DM. The findings of meta-analysis revealed no significant association of VDR gene SNPs with T1DM risk in the overall population. However, the subgroup analysis indicated significant associations between FokI and BsmI polymorphism and T1DM risk in African and American population. As a limitation, we did not evaluate a number of VDR gene SNPs that might act in interaction with environmental factors to determine the fate of T1DM pathogenicity. Further investigations on the VDR, above and beyond the genetic as well as traditional risk factors, may confer a possibility for identification of critical susceptibility factors in the disease development, which might be applicable in the personalized medicine for better and optimized therapy of T1DM patients.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

Abbreviations

T1DM:

Type 1 diabetes mellitus

VDR:

Vitamin D receptor

Vitamin D:

VitD

SNP:

Single nucleotide polymorphisms

IL:

Interleukin

PRISMA:

Preferred Reporting Items for Systematic reviews and Meta-Analyses

NOS:

Newcastle-Ottawa Scale

UVR:

Ultraviolet radiation

Th:

T helper

TNF:

Tumor necrosis factor

IFN:

Interferon

HWE:

Hardy–Weinberg equilibrium

PCR- RFLP:

Polymerase chain reaction-restriction fragment length polymorphism

References

  1. Gupta G, et al. A clinical update on metformin and lung cancer in diabetic patients. Panminerva Med. 2018;60(2):70–5.

    Article  PubMed  Google Scholar 

  2. Atkinson MA, Eisenbarth GS, Michels AW. Type 1 diabetes. Lancet. 2014;383(9911):69–82.

    Article  PubMed  Google Scholar 

  3. Miettinen ME, et al. Genetic determinants of serum 25-hydroxyvitamin D concentration during pregnancy and type 1 diabetes in the child. PLoS One. 2017;12(10):e0184942.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Diaz-Valencia PA, Bougnères P, Valleron A-J. Global epidemiology of type 1 diabetes in young adults and adults: a systematic review. BMC Public Health. 2015;15(1):255.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Todd JA, et al. Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet. 2007;39(7):857.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Vaidya A, Williams JS. The relationship between vitamin D and the renin-angiotensin system in the pathophysiology of hypertension, kidney disease, and diabetes. Metabolism. 2012;61(4):450–8.

    Article  CAS  PubMed  Google Scholar 

  7. Riek AE, et al. Vitamin D suppression of endoplasmic reticulum stress promotes an antiatherogenic monocyte/macrophage phenotype in type 2 diabetic patients. J Biol Chem. 2012;287(46):38482–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Makoui MH, Imani D, Motallebnezhad M, Azimi M, Razi B. Vitamin D receptor gene polymorphism and susceptibility to asthma: Meta-analysis based on 17 case-control studies. Ann Allergy Asthma Immunol. 2020;124(1):57–69.

    Article  CAS  PubMed  Google Scholar 

  9. Bagheri-Hosseinabadi Z, et al. Vitamin D receptor (VDR) gene polymorphism and risk of rheumatoid arthritis (RA): systematic review and meta-analysis. Clin Rheumatol. 2020.

  10. Bhalla AK, et al. 1, 25-Dihydroxyvitamin D3 inhibits antigen-induced T cell activation. J Immunol. 1984;133(4):1748–54.

    CAS  PubMed  Google Scholar 

  11. Lemire JM. Immunomodulatory actions of 1, 25-dihydroxyvitamin D3. J Steroid Biochem Mol Biol. 1995;53(1–6):599–602.

    Article  CAS  PubMed  Google Scholar 

  12. Riachy R, et al. 1, 25-Dihydroxyvitamin D 3 protects human pancreatic islets against cytokine-induced apoptosis via down-regulation of the Fas receptor. Apoptosis. 2006;11(2):151–9.

    Article  CAS  PubMed  Google Scholar 

  13. Trembleau S, et al. The role of IL-12 in the induction of organ-specific autoimmune diseases. Immunol Today. 1995;16(8):383–6.

    Article  CAS  PubMed  Google Scholar 

  14. Uitterlinden AG, et al. Genetics and biology of vitamin D receptor polymorphisms. Gene. 2004;338(2):143–56.

    Article  CAS  PubMed  Google Scholar 

  15. Panierakis C, et al. Vitamin D receptor gene polymorphisms and susceptibility to type 1 diabetes in Crete, Greece. Clin Immunol. 2009;133(2):276–81.

    Article  CAS  PubMed  Google Scholar 

  16. Wang Q, et al. Quantitative assessment of the associations between four polymorphisms (FokI, ApaI, BsmI, TaqI) of vitamin D receptor gene and risk of diabetes mellitus. Mol Biol Rep. 2012;39(10):9405–14.

    Article  CAS  PubMed  Google Scholar 

  17. Ferrari S, et al. Vitamin D receptor gene start codon polymorphisms (FokI) and bone mineral density: interaction with age, dietary calcium, and 3′-end region polymorphisms. J Bone Miner Res. 1998;13(6):925–30.

    Article  CAS  PubMed  Google Scholar 

  18. Imani D, et al. Association between vitamin D receptor (VDR) polymorphisms and the risk of multiple sclerosis (MS): an updated meta-analysis. BMC Neurol. 2019;19(1):339.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Valdivielso JM, Fernandez E. Vitamin D receptor polymorphisms and diseases. Clin Chim Acta. 2006;371(1–2):1–12.

    Article  CAS  PubMed  Google Scholar 

  20. Tizaoui K, et al. Contribution of VDR polymorphisms to type 1 diabetes susceptibility: systematic review of case–control studies and meta-analysis. J Steroid Biochem Mol Biol. 2014;143:240–9.

    Article  CAS  PubMed  Google Scholar 

  21. Guo S-W, et al. Meta-analysis of vitamin D receptor polymorphisms and type 1 diabetes: a HuGE review of genetic association studies. Am J Epidemiol. 2006;164(8):711–24.

    Article  PubMed  Google Scholar 

  22. Zhang J, et al. Polymorphisms in the vitamin D receptor gene and type 1 diabetes mellitus risk: an update by meta-analysis. Mol Cell Endocrinol. 2012;355(1):135–42.

    Article  CAS  PubMed  Google Scholar 

  23. Sahin OA, et al. Association of vitamin D receptor polymorphisms and type 1 diabetes susceptibility in children: a meta-analysis. Endocr Connections. 2017;6(3):159–71.

    Article  CAS  Google Scholar 

  24. Qin W-H, et al. A meta-analysis of association of vitamin D receptor BsmI gene polymorphism with the risk of type 1 diabetes mellitus. J Recept Signal Transduct. 2014;34(5):372–7.

    Article  CAS  Google Scholar 

  25. Wang G, et al. Associations between two polymorphisms (FokI and BsmI) of vitamin D receptor gene and type 1 diabetes mellitus in Asian population: a meta-analysis. PLoS One. 2014;9(3):e89325.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Moher D, et al. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Ann Intern Med. 2009;151(4):264–9.

    Article  PubMed  Google Scholar 

  27. Stang A. Critical evaluation of the Newcastle-Ottawa scale for the assessment of the quality of nonrandomized studies in meta-analyses. Eur J Epidemiol. 2010;25(9):603–5.

    Article  PubMed  Google Scholar 

  28. Huedo-Medina TB, et al. Assessing heterogeneity in meta-analysis: Q statistic or I2 index? Psychol Methods. 2006;11(2):193.

    Article  PubMed  Google Scholar 

  29. DerSimonian R, Laird N. Meta-analysis in clinical trials control. Clin Trials. 1986;7:177–88. Find this article online.

    Article  CAS  Google Scholar 

  30. Mantel N, Haenszel W. Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. 1959;22(4):719–48.

    CAS  PubMed  Google Scholar 

  31. Egger M, et al. Bias in meta-analysis detected by a simple, graphical test. Bmj. 1997;315(7109):629–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Begg CB, Mazumdar M. Operating characteristics of a rank correlation test for publication bias. Biometrics. 1994;50:1088–101.

    Article  CAS  PubMed  Google Scholar 

  33. Ban Y, et al. Vitamin D receptor initiation codon polymorphism influences genetic susceptibility to type 1 diabetes mellitus in the Japanese population. BMC Med Genet. 2001;2(1):7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Fassbender W, et al. VDR gene polymorphisms are overrepresented in German patients with type 1 diabetes compared to healthy controls without effect on biochemical parameters of bone metabolism. Horm Metab Res. 2002;34(06):330–7.

    Article  CAS  PubMed  Google Scholar 

  35. Gyorffy B, et al. Gender-specific association of vitamin D receptor polymorphism combinations with type 1 diabetes mellitus. Eur J Endocrinol. 2002;147(6):803–8.

    Article  CAS  PubMed  Google Scholar 

  36. Turpeinen H, et al. Vitamin D receptor polymorphisms: no association with type 1 diabetes in the Finnish population. Eur J Endocrinol. 2003;149(6):591–6.

    Article  CAS  PubMed  Google Scholar 

  37. Martí G, Audí L, Esteban C, Oyarzábal M, Chueca M, Gussinyé M, Yeste D, Fernández-Cancio M, Andaluz P, Carrascosa A. Asociación de los polimorfismos del gen del receptor de la vitamina D con la diabetes mellitus tipo 1 en dos poblaciones españolas. Med Clín. 2004;123(8):286–90.

    Google Scholar 

  38. Pedro JS, et al. Heterogeneity of vitamin D receptor gene association with celiac disease and type 1 diabetes mellitus. Autoimmunity. 2005;38(6):439–44.

    Article  CAS  Google Scholar 

  39. Zemunik T, et al. FokI polymorphism, vitamin D receptor, and interleukin-1 receptor haplotypes are associated with type 1 diabetes in the Dalmatian population. J Mol Diagn. 2005;7(5):600–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Capoluongo E, et al. Slight association between type 1 diabetes and “ff” VDR FokI genotype in patients from the Italian Lazio region. Lack of association with diabetes complications. Clin Biochem. 2006;39(9):888–92.

    Article  CAS  PubMed  Google Scholar 

  41. Lemos MC, et al. Lack of association of vitamin D receptor gene polymorphisms with susceptibility to type 1 diabetes mellitus in the Portuguese population. Hum Immunol. 2008;69(2):134–8.

    Article  CAS  PubMed  Google Scholar 

  42. Israni N, et al. Interaction of vitamin D receptor with HLA DRB1* 0301 in type 1 diabetes patients from North India. PLoS One. 2009;4(12):e8023.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Mory DB, et al. Prevalence of vitamin D receptor gene polymorphisms FokI and BsmI in Brazilian individuals with type 1 diabetes and their relation to β-cell autoimmunity and to remaining β-cell function. Hum Immunol. 2009;70(6):447–51.

    Article  CAS  PubMed  Google Scholar 

  44. Yavuz DG, et al. Vitamin D receptor gene BsmI, FokI, ApaI, TaqI polymorphisms and bone mineral density in a group of Turkish type 1 diabetic patients. Acta Diabetol. 2011;48(4):329–36.

    Article  CAS  Google Scholar 

  45. Yokota I, et al. Association between vitamin D receptor genotype and age of onset in juvenile Japanese patients with type 1 diabetes. Diabetes Care. 2002;25(7):1244.

    Article  PubMed  Google Scholar 

  46. Bonakdaran, S., et al., Vitamin D receptor gene polymorphisms in type 1 diabetes mellitus: a new pattern from Khorasan province, Islamic Republic of Iran. 2012.

    Google Scholar 

  47. Sahin SB, et al. Fas, Fas ligand, and vitamin D receptor FokI gene polymorphisms in patients with type 1 diabetes mellitus in the Aegean region of Turkey. Genet Test Mol Biomark. 2012;16(10):1179–83.

    Article  CAS  Google Scholar 

  48. Mohammadnejad Z, et al. Association between vitamin D receptor gene polymorphisms and type 1 diabetes mellitus in Iranian population. Mol Biol Rep. 2012;39(2):831–7.

    Article  CAS  PubMed  Google Scholar 

  49. Vedralová M, et al. Polymorphisms in the vitamin D receptor gene and parathyroid hormone gene in the development and progression of diabetes mellitus and its chronic complications, diabetic nephropathy and non-diabetic renal disease. Kidney Blood Press Res. 2012;36(1):1–9.

    Article  PubMed  CAS  Google Scholar 

  50. Greer RM, et al. Serum vitamin D levels are lower in Australian children and adolescents with type 1 diabetes than in children without diabetes. Pediatr Diabetes. 2013;14(1):31–41.

    Article  CAS  PubMed  Google Scholar 

  51. Hamed EO, et al. Vitamin D level and Fok-I vitamin D receptor gene polymorphism in Egyptian patients with type-1 diabetes. Egypt J Immunol. 2013;20(2):1–10.

    PubMed  Google Scholar 

  52. Abd-Allah SH, et al. Vitamin D status and vitamin D receptor gene polymorphisms and susceptibility to type 1 diabetes in Egyptian children. Gene. 2014;536(2):430–4.

    Article  CAS  PubMed  Google Scholar 

  53. El-Kafoury AA, et al. The association of polymorphic sites in some genes with type 1 diabetes mellitus in a sample of Egyptian children. Egypt J Med Hum Genet. 2014;15(3):265–72.

    Article  Google Scholar 

  54. Nasreen M, et al. Serum vitamin D levels and gene polymorphisms (Fok1 and Apa1) in children with type I diabetes and healthy controls. JPMA. 2016;66:1215.

    Google Scholar 

  55. Mukhtar M, et al. Vitamin D receptor gene polymorphisms influence T1D susceptibility among Pakistanis. Int J Genomics. 2017;2017:4171254.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Ali R, et al. Evaluation of vitamin D receptor gene polymorphisms (Fok-I and Bsm-I) in T1DM Saudi children. J Clin Lab Anal. 2018;32(5):e22397.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  57. Rasoul MA, et al. Relationship of four vitamin D receptor gene polymorphisms with type 1 diabetes mellitus susceptibility in Kuwaiti children. BMC Pediatr. 2019;19(1):71.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Chang TJ, et al. Vitamin D receptor gene polymorphisms influence susceptibility to type 1 diabetes mellitus in the Taiwanese population. Clin Endocrinol. 2000;52(5):575–80.

    Article  CAS  Google Scholar 

  59. Škrabić V, et al. Vitamin D receptor polymorphism and susceptibility to type 1 diabetes in the Dalmatian population. Diabetes Res Clin Pract. 2003;59(1):31–5.

    Article  PubMed  Google Scholar 

  60. Bianco M, et al. Vitamin D receptor polymorphisms: are they really associated with type 1 diabetes? Eur J Endocrinol. 2004;151(5):641–2.

    Article  CAS  PubMed  Google Scholar 

  61. García D, et al. VDR polymorphisms influence the immune response in type 1 diabetic children from Santiago, Chile. Diabetes Res Clin Pract. 2007;77(1):134–40.

    Article  PubMed  CAS  Google Scholar 

  62. Cheon CK, et al. Vitamin D receptor gene polymorphisms and type 1 diabetes mellitus in a Korean population. Pediatr Int. 2015;57(5):870–4.

    Article  CAS  PubMed  Google Scholar 

  63. Khalid KE. Vitamin D receptor gene polymorphisms in Sudanese children with type 1 diabetes. AIMS Genet. 2016;3(3):167–76.

    Article  Google Scholar 

  64. Iyer A, et al. Relationship between vitamin D receptor gene polymorphisms and type 1 diabetes mellitus in Saudi patients. Int J Pharmacol. 2017;13(8):1092–7.

    Article  CAS  Google Scholar 

  65. Ahmed AE-A, et al. Vitamin D receptor rs7975232, rs731236 and rs1544410 single nucleotide polymorphisms, and 25-hydroxyvitamin D levels in Egyptian children with type 1 diabetes mellitus: effect of vitamin D co-therapy. Diabetes Metab Syndr Obes. 2019;12:703.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Hauache O, et al. Vitamin D receptor gene polymorphism: correlation with bone mineral density in a Brazilian population with insulin-dependent diabetes mellitus. Osteoporos Int. 1998;8(3):204–10.

    Article  CAS  PubMed  Google Scholar 

  67. Motohashi Y, et al. Vitamin D receptor gene polymorphism affects onset pattern of type 1 diabetes. J Clin Endocrinol Metab. 2003;88(7):3137–40.

    Article  CAS  PubMed  Google Scholar 

  68. Shimada A, et al. Evidence for association between vitamin D receptor BsmI polymorphism and type 1 diabetes in Japanese. J Autoimmun. 2008;30(4):207–11.

    Article  CAS  PubMed  Google Scholar 

  69. Tawfeek M, Habib F, Saultan EM. Vitamin D receptor BsmI gene polymorphisms and gestational diabetes mellitus: a Saudi study. Br J Med Med Res. 2011;1(4):459–68.

    Article  Google Scholar 

  70. Al-Moubarak SH. Gender-specific association of vitamin D receptor polymorphism Bsm-I with type 1 diabetes mellitus. Int J Pharm Sci Rev Res. 2013;21(2):254–7.

    CAS  Google Scholar 

  71. Martí G, et al. Asociación de los polimorfismos del gen del receptor de la vitamina D con la diabetes mellitus tipo 1 en dos poblaciones españolas. Med Clin. 2004;123(8):286–90.

    Article  Google Scholar 

  72. Ponsonby A-L, et al. Variation in associations between allelic variants of the vitamin D receptor gene and onset of type 1 diabetes mellitus by ambient winter ultraviolet radiation levels: a meta-regression analysis. Am J Epidemiol. 2008;168(4):358–65.

    Article  PubMed  Google Scholar 

  73. Kimlin MG, Olds WJ, Moore MR. Location and vitamin D synthesis: is the hypothesis validated by geophysical data? J Photochem Photobiol B Biol. 2007;86(3):234–9.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Mrs. Maryam Izad for all her support.

Disclosure of conflict of interest

Not applicable.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Author information

Authors and Affiliations

Authors

Contributions

NZ participated in study design and manuscript drafting. RB, participated in literature search and contributed to manuscript drafting. MHM analyzed the data and participated in drafting the manuscript. SA analyzed and interpreted the data and participated in manuscript drafting. PM contributed to data analysis and prepared the original draft. BR performed the literature search, analyzed data, and participated in manuscript drafting. DI performed the literature search, developed the main idea, and participated in manuscript drafting. MY performed the literature search and participated in manuscript drafting. HM performed data interpretation and participated in manuscript drafting. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Na Zhai or Haleh Mikaeili.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhai, N., Bidares, R., Makoui, M.H. et al. Vitamin D receptor gene polymorphisms and the risk of the type 1 diabetes: a meta-regression and updated meta-analysis. BMC Endocr Disord 20, 121 (2020). https://doi.org/10.1186/s12902-020-00575-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12902-020-00575-8

Keywords