Association of vitamin D receptor polymorphisms and type 1 diabetes susceptibility in children: a meta-analysis

in Endocrine Connections
Authors:
Ozlem Atan Sahin Department of Pediatrics, Acıbadem University School of Medicine, Atasehir, Istanbul, Turkey

Search for other papers by Ozlem Atan Sahin in
Current site
Google Scholar
PubMed
Close
,
Damla Goksen Department of Pediatric Endocrinology, Faculty of Medicine, Ege University, Bornova, Izmir, Turkey

Search for other papers by Damla Goksen in
Current site
Google Scholar
PubMed
Close
,
Aysel Ozpinar Department of Biochemistry, Acıbadem University, School of Medicine, Atasehir, Istanbul, Turkey

Search for other papers by Aysel Ozpinar in
Current site
Google Scholar
PubMed
Close
,
Muhittin Serdar Department of Biochemistry, Acıbadem University, School of Medicine, Atasehir, Istanbul, Turkey

Search for other papers by Muhittin Serdar in
Current site
Google Scholar
PubMed
Close
, and
Huseyin Onay Department of Medical Genetics, Faculty of Medicine, Ege University, Bornova, Izmir, Turkey

Search for other papers by Huseyin Onay in
Current site
Google Scholar
PubMed
Close

Open access

Sign up for journal news

Background

There have been studies focused on FokI, BsmI, ApaI and TaqI polymorphisms of the vitamin D receptor (VDR) gene and susceptibility to type 1 diabetes mellitus with controversial results.

Methods

This present study is a meta-analysis investigating the association between FokI, ApaI, TaqI and BsmI polymorphisms of VDR gene and type 1 DM in children. A literature search was performed using Medline, EMBASE, Cochrane and PubMed. Any study was considered eligible for inclusion if at least one of FokI, ApaI, TaqI and BsmI polymorphisms was determined, and outcome was type 1 DM at pediatric age.

Results

A total of 9 studies comprising 1053 patients and 1017 controls met the study inclusion criteria. The pooled odds ratios (ORs) of the FokI, ApaI, TaqI and BsmI polymorphisms were combined and calculated. Forest plots and funnel plots of the OR value distributions were drawn. Our meta-analysis has demonstrated statistically significant associations between DM1 and VDR genotypes, BsmIBB (P < 0.05), BsmIBb, (P < 0.05), BsmIbb (P < 0.05), TaqITT (P < 0.05) and TaqItt (P < 0.05) in children.

Conclusion

The results indicated that BsmIBB, BsmIBb and TaqItt polymorphisms were associated with an increased risk of type 1 DM, whereas BsmIbb and TaqITT had protective effect for type 1 DM in children.

Abstract

Background

There have been studies focused on FokI, BsmI, ApaI and TaqI polymorphisms of the vitamin D receptor (VDR) gene and susceptibility to type 1 diabetes mellitus with controversial results.

Methods

This present study is a meta-analysis investigating the association between FokI, ApaI, TaqI and BsmI polymorphisms of VDR gene and type 1 DM in children. A literature search was performed using Medline, EMBASE, Cochrane and PubMed. Any study was considered eligible for inclusion if at least one of FokI, ApaI, TaqI and BsmI polymorphisms was determined, and outcome was type 1 DM at pediatric age.

Results

A total of 9 studies comprising 1053 patients and 1017 controls met the study inclusion criteria. The pooled odds ratios (ORs) of the FokI, ApaI, TaqI and BsmI polymorphisms were combined and calculated. Forest plots and funnel plots of the OR value distributions were drawn. Our meta-analysis has demonstrated statistically significant associations between DM1 and VDR genotypes, BsmIBB (P < 0.05), BsmIBb, (P < 0.05), BsmIbb (P < 0.05), TaqITT (P < 0.05) and TaqItt (P < 0.05) in children.

Conclusion

The results indicated that BsmIBB, BsmIBb and TaqItt polymorphisms were associated with an increased risk of type 1 DM, whereas BsmIbb and TaqITT had protective effect for type 1 DM in children.

Introduction

Type 1 diabetes (DM1) is a complex disease characterized by the autoimmune destruction of pancreatic β cells. Vitamin D is an immune regulatory hormone that exerts its effects through highly polymorphic VDR that belongs to steroid-receptor superfamily, and it is expressed in many cell types such as lymphocytes and antigen-presenting cells (APCs) (1). During the last decade, VDR gene polymorphisms have been shown to be associated with autoimmune pathologies (2). Vitamin D seems to downregulate type 1 helper (Th-1) cells, by decreasing their proliferation and inhibiting the production of cytokines such as IL-2, TNF-α and interferon-γ (3, 4). For many years, the strongest genetic contribution to DM1 susceptibility had been attributed to the presence of human leukocyte antigen region (HLA) on chromosome 6 (5, 6). Recently, single nucleotide polymorphisms (SNPs) in the VDR gene have been investigated namely FokI F > f (rs10735810), BsmI B > b (rs1544410), ApaI A > a (rs7975232) and TaqI t T > t (rs731236). Allele F of the FokI SNP creates an alternative ATG initiation codon in exon 2 leading to a three amino acid longer VDR protein. The ApaI, BsmI and TaqI polymorphisms take place near the 3′ end of the VDR gene; BsmI and ApaI SNPs are located in intron 8, and the TaqI is a silent SNP in exon 9. Several studies with small data sets that suggested an association between these SNPs and type 1 DM had inconsistent results (7, 8, 9, 10). These inconclusive results had been attributed to ethnic diversities among populations or due to the environmental factors involved in type 1 DM pathogenesis. BsmI is strongly linked with 3 poly(A) microsatellite repeat in the 3′ untranslated region, which may influence VDR mRNA stability (11). BsmI, ApaI and TaqI polymorphisms were also designated as polymorphisms without functional effects in several studies (12, 13). There are some discrepancies in the literature with small data sets, regarding the relative importance of polymorphisms of VDR genes ApaI, BsmI, FokI and TaqI. In our meta-analysis, we included a total of 9 studies comprising 1053 patients and 1017 controls that indicated statistically significant associations between DM1 and VDR genotypes, BsmIBB (P < 0.05), BsmIBb (P < 0.05), BsmIbb (P < 0.05), TaqITT (P < 0.05) and TaqItt (P < 0.05) in children. This present meta-analysis involves a large data set to investigate the associations between type 1 DM and VDR gene polymorphisms ApaI, BsmI, FokI and TaqI, accordingly demonstrate more reliable statistical results to rule out genotype–phenotype correlation of type 1 DM in children.

Methods

Search strategy criteria

For meta-analysis, all published studies evaluating the associations between type 1 DM and FokI, ApaI, TaqI and BsmI polymorphisms that are investigated in patients diagnosed as DM1 at pediatric age are included. A literature search for the MeSH terms ‘type 1 Diabetes mellitus’ or ‘DM1’ was performed by O A, D G and M S Medline, Cochrane and PubMed abstracts were reviewed for relevance. No language and date of study restriction were applied to search strategy. Search to include the eligible studies ended on 05/14/2016. Any study was considered to be eligible for inclusion if it met the following criteria: (1) the publication was an association study of the case control type, (2) at least one of the FokI, ApaI, TaqI and BsmI polymorphism was determined, (3) the outcome was DM in children and (4) there was at least one unrelated control group.

Data extraction

Study selection and data extraction were performed independently by three authors (O A, D G and M S) based on a customized database for extraction. For each study, the following information was collected: first author, year and location of the study, average age at the time of diagnosis, ethnicity, number of participants, number of cases and controls and number of the genotypes in cases and controls. The disagreements were resolved between the reviewers by consensus. For quality assessment, six domains were assessed. Those were representativeness of classes, representativeness of the controls, ascertainment of DM1, genotypic examination and association of assessment. The primary outcome considered in the meta-analysis was the association between DM1 and the presence of FokI, ApaI, TaqI or BsmI polymorphism at pediatric age. For the primary analysis and to allow appropriate comparison of all studies, cases and controls were classified based on FokI, ApaI, TaqI and BsmI genotypes.

Statistical analysis

The odds ratios (OR) with 95% confidence intervals, representativeness of controls, ascertainment of DM1, ascertainment of controls, genotypic examination and association assessments were done. The primary outcome considered in the meta-analysis was the association between DM1 and the presence of FokI, ApaI, TaqI or BsmI polymorphisms. MedCalc Software Acacialaan 22, 8400 (Ostend, Belgium) was used to perform meta-analysis. The odds ratios (OR) of the genetic polymorphisms were combined and calculated, and the funnel plots were drawn. All of the four studied SNPs (FokI, ApaI, TaqI and BsmI) were diallelic, and we calculated summary odds ratios incorporating both within- and between-study variation using a random effects model proposed by DerSimonian and Laird (14).

Results

Our search yielded a total of 50 references. After screening the titles and abstracts, 41 studies were excluded because they were not considered relevant to the study topic, leaving 9 potentially eligible studies (Fig. 1) (15, 16, 17, 18, 19, 20, 21, 22, 23).

Figure 1
Figure 1

Flow chart for identification and selection of studies.

Citation: Endocrine Connections 6, 3; 10.1530/EC-16-0110

In the 9 published papers included in the meta-analysis, ApaI, BsmI, FokI and TaqI polymorphisms were investigated in pediatric population as case–control studies (Table 1).

Table 1

Characteristics features of studies included in the meta-analysis of ApaI, BsmI, FokI and TaqI polymorphisms in the vitamin D receptor gene.

First author Year Region Mean age of cases/diagnosis (years) Cases Source of controls Mean age of controls (years) Controls
Diego Garcia 2007 Santiago-Chile 9.3 ± 4.2 216 Unrelated children 10.3 ± 2.5 203
J I San Pedro 2005 Bilbao-Spain 14.5 ± 9.9 71 Healthy blood donors 88
Tatijana Semunik 2005 Split-Croatia 8.6 ± 4.3 132 Unrelated children 8.2 ± 4.9 232
Vaselin Scrabic 2003 Split-Croatia 8.6 ± 4.3 134 Unrelated children 8.24 ± 4.9 132
Balazs Gyorffy 2002 Budapest-Hungary 5.8 ± 3.2 107 Healthy blood donors 103
Tien-Jyung Chang 2000 Han Chinese-Taiwan 8.8 ± 5.6 157 Healthy subjects 248
Charalambos Panierakis 2009 Crete-Greece Children 100 Unrelated children 96
Greear R M 2013 Brisbane-Australia <15 55 Healthy subjects <15 50
Chong-Kun Cheon 2015 Pusan-South Korea 10.28 ± 3.23 81 Healthy children 9.98 ± 3.56 113

Eight studies on the ApaI-type 1 diabetes association recruited 921 cases/patients and 1033 controls, whereas seven studies on the BsmI polymorphism recruited 866 cases and 983 controls. For the FokI polymorphism, five studies included 465 cases and 569 controls, whereas eight studies on the TagI polymorphism included 921 cases and 1033 controls. Individual and pooled odds ratio estimates of four single-nucleotide polymorphisms in the vitamin receptor gene, P values testing Hardy–Weinberg proportion, test for heterogeneity (Tables 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 and 13) and funnels plots (Figs 2 and 3) are documented for BsmI and TaqI, respectively.

Figure 2
Figure 2

(A, B and C) Forest plots showing individual and pooled odds ratio estimates of BsmI BB, BsmI Bb, BsmIbb polymorphisms, respectively.

Citation: Endocrine Connections 6, 3; 10.1530/EC-16-0110

Figure 3
Figure 3

(A and B) Forest plots showing individual and pooled odds ratios of TaqITT and TaqItt polymorphisms, respectively.

Citation: Endocrine Connections 6, 3; 10.1530/EC-16-0110

Table 2

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for ApaI AA polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 16/157 13/248 2.051 0.958–4.391
Balazs Gyorffy 2002 33/107 23/103 1.551 0.835–2.881
Vaselin Scrabic 2003 66/134 51/132 1.542 0.947–2.509
J J San Pedro 2005 15/71 28/88 0.574 0.278–1.185
Diego Garcia 2007 54/216 43/203 1.240 0.786–1.957
C Panierakis 2009 23/100 37/96 0.476 0.256–0.886
Greear R M 2013 15/55 12/50 1.187 0.493–2.861
Chung Cheon 2015 5/81 9/113 0.760 0.245–2.359
Total (fixed effects) 227/921 216/1033 1.113 0.893–1.388 0.954 0.340
Total (random effects) 227/921 216/1033 1.081 0.755–1.547 0.425 0.671
Q 16.3471
DF 7
Significance level P = 0.0221
I2 (inconsistency) 57.18%
95% CI for I2 5.93–80.51

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.142857 P = 0.5484 (low power); Egger: bias = −1.369742 (95% CI = −6.646958 to 3.907474) P = 0.5488; Harbord–Egger: bias = −1.179831 (92.5% CI = −5.949539 to 3.589877) P = 0.6138.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 3

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for ApaI Aa polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 76/157 105/248 1.278 0.855–1.910
Balazs Gyorffy 2002 45/107 54/103 0.659 0.382–1.136
Vaselin Scrabic 2003 52/134 66/132 0.634 0.390–1.032
J J San Pedro 2005 37/71 43/88 1.139 0.609–2.129
Diego Garcia 2007 115/216 125/203 0.710 0.481–1.048
C Panierakis 2009 58/100 57/96 0.945 0.535–1.669
Greear R M 2013 24/55 32/50 0.435 0.198–0.956
Chung Cheon 2015 32/81 34/113 1.517 0.833–2.765
Total (fixed effects) 439/921 516/1033 0.873 0.729–1.046 −1.473 0.141
Total (random effects) 439/921 516/1033 0.866 0.664–1.131 −1.054 0.292
Q 14.2512
DF 7
Significance level P = 0.0469
I2 (inconsistency) 50.88%
95% CI for I2 0.00–78.02

Bias indicators: Begg–Mazumdar: Kendall’s Tau = 0 P = 0.9049 (low power); Egger: bias = −0.928241 (95% CI = −7.113964 to 5.257482) P = 0.7261; Harbord–Egger: bias = −0.806491 (92.5% CI = −6.321591 to 4.708608) P = 0.7638.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 4

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for ApaI aa polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 65/157 130/248 0.641 0.428–0.960
Balazs Gyorffy 2002 27/107 26/103 1.000 0.536–1.863
Vaselin Scrabic 2003 16/134 15/132 1.058 0.500–2.238
J J San Pedro 2005 19/71 17/88 1.526 0.724–3.217
Diego Garcia 2007 44/216 35/203 1.228 0.751–2.009
C Panierakis 2009 15/100 6/96 2.647 0.982–7.139
Greear R M 2013 11/55 11/50 0.886 0.346–2.270
Chung Cheon 2015 44/81 70/113 0.731 0.409–1.304
Total (fixed effects) 241/921 310/1033 0.956 0.772–1.182 −0.419 0.675
Total (random effects) 241/921 310/1033 1.005 0.754–1.339 0.0349 0.972
Q 11.2531
DF 7
Significance level P = 0.1280
I2 (inconsistency) 37.79%
95% CI for I2 0.00–72.53

Bias indicators: Begg–Mazumdar: Kendall’s Tau = 0.428571 P = 0.1789 (low power); Egger: bias = 2.766246 (95% CI = −0.351565 to 5.884058) P = 0.073; Harbord–Egger: bias = 2.78392 (92.5% CI = −0.118976 to 5.686817) P = 0.0847.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 5

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for BsmI BB polymorphism, respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 4/157 1/248 6.458 0.715–58.314
Balazs Gyorffy 2002 17/107 19/103 0.835 0.407–1.713
Vaselin Scrabic 2003 24/134 17/132 1.476 0.752–2.896
J J San Pedro 2005 15/71 17/88 1.119 0.514–2.435
Diego Garcia 2007 21/216 14/203 1.454 0.718–2.943
C Panierakis 2009 38/100 23/96 1.945 1.048–3.611
Chung Cheon 2015 0/81 1/113 0.460 0.0185–11.440
Total (fixed effects) 119/866 92/983 1.397 1.034–1.888 2.174 0.030
Total (random effects) 119/866 92/983 1.386 1.021–1.880 2.092 0.036
Q 5.7387
DF 6
Significance level P = 0.4531
I2 (inconsistency) 0.00%
95% CI for I2 0.00–69.98

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.333333 P = 0.2389 (low power); Egger: bias = 0.081309 (95% CI = −2.562044 to 2.724662) P = 0.94; Harbord–Egger: bias = 0.057432 (92.5% CI = −2.538737 to 2.653601) P = 0.9624.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 6

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for BsmI Bb polymorphism, respectively.

Study Intervention Controls Odds ratio 95% CI z P
Chang 2000 16/157 16/248 1.645 0.798–3.393
Balazs Gyorffy 2002 53/107 44/103 1.316 0.764–2.268
Vaselin Scrabic 2003 58/134 74/132 0.598 0.368–0.971
J J San Pedro 2005 40/71 44/88 1.290 0.688–2.418
Diego Garcia 2007 110/216 74/203 1.809 1.224–2.674
C Panierakis 2009 43/100 23/96 2.394 1.296–4.422
Chung Cheon 2015 13/81 4/113 5.210 1.632–16.633
Total (fixed effects) 333/866 279/983 1.430 1.160–1.762 3.347 0.001
Total (random effects) 333/866 279/983 1.534 1.001–2.350 1.966 0.049
Q 21.6238
DF 6
Significance level P = 0.0014
I2 (inconsistency) 72.25%
95% CI for I2 40.02–87.16

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.333333 P = 0.3813 (low power); Egger: bias = 2.518064 (95% CI = −4.133965 to 9.170093) P = 0.3752; Harbord–Egger: bias = 2.595692 (92.5% CI = −3.668787 to 8.860172) P = 0.3955.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 7

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for BsmI bb polymorphism, respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 137/157 231/248 0.504 0.255–0.995
Balazs Gyorffy 2002 35/107 40/103 0.766 0.435–1.348
Vaselin Scrabic 2003 52/134 41/132 1.407 0.848–2.336
J J San Pedro 2005 16/71 27/88 0.657 0.321–1.347
Diego Garcia 2007 77/216 115/203 0.424 0.286–0.628
C Panierakis 2009 15/100 20/96 0.671 0.321–1.402
Chung Cheon 2015 68/81 108/113 0.242 0.0826–0.710
Total (fixed effects) 400/866 582/983 0.632 0.508–0.786 −4.114 <0.001
Total (random effects) 400/866 582/983 0.624 0.418–0.933 −2.298 0.022
Q 17.5208
DF 6
Significance level P = 0.0075
I2 (inconsistency) 65.76%
95% CI for I2 23.26–84.72

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.142857 P = 0.5619 (low power); Egger: bias = −0.656941 (95% CI = −7.053883 to 5.74) P = 0.8023; Harbord–Egger: bias = 0.561504 (92.5% CI = −6.312289 to 5.189281) P = 0.8354.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 8

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for FokI FF polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
Balazs Gyorffy 2002 36/107 29/103 1.294 0.719–2.328
Tatijana Semunik 2005 42/132 73/232 1.016 0.642–1.609
J J San Pedro 2005 31/71 41/88 0.888 0.474–1.666
C Panierakis 2009 64/100 50/96 1.636 0.923–2.898
Greear R M 2013 21/55 28/50 0.485 0.223–1.058
Total (fixed effects) 194/465 221/569 1.057 0.817–1.367 0.420 0.675
Total (random effects) 194/465 221/569 1.036 0.733–1.465 0.202 0.840
Q 6.8448
DF 4
Significance level P = 0.1443
I2 (inconsistency) 41.56%
95% CI for I2 0.00–78.48

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.6 P = 0.00833 (low power); Egger: bias = −3.653382 (95% CI = −14.60785 to 7.301086) P = 0.3664; Harbord–Egger: bias = −3.796013 (92.5% CI = −13.262673 to 5.670648) P = 0.3612.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 9

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for FokI Ff polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
Balazs Gyorffy 2002 49/107 56/103 0.709 0.412–1.221
Tatijana Semunik 2005 63/132 136/232 0.645 0.419–0.991
J J San Pedro 2005 35/71 39/88 1.222 0.652–2.287
C Panierakis 2009 31/100 43/96 0.554 0.309–0.993
Greear R M 2013 21/55 22/50 0.786 0.361–1.714
Total (fixed effects) 199/465 296/569 0.724 0.564–0.929 −2.538 0.011
Total (random effects) 199/465 296/569 0.723 0.563–0.929 −2.535 0.011
Q 3.8098
DF 4
Significance level P = 0.4324
I2 (inconsistency) 0.00%
95% CI for I2 0.00–79.45

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.4 P = 0.4833 (low power); Egger: bias = 1.95378 (95% CI = 5.613844–9.521404) P = 0.4715; Harbord–Egger: bias = 2.01803 (92.5% CI = −4.313428 to 8.349488) P = 0.4557.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 10

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for FokI ff polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
Balazs Gyorffy 2002 22/107 18/103 1.222 0.612–2.441
Tatijana Semunik 2005 29/132 23/232 2.558 1.410–4.643
J J San Pedro 2005 5/71 8/88 0.758 0.237–2.426
C Panierakis 2009 1/100 7/96 0.128 0.0155–1.064
Greear R M 2013 7/55 5/50 1.312 0.388–4.435
Total (fixed effects) 64/465 61/569 1.374 0.943–2.003 1.656 0.098
Total (random effects) 64/465 61/569 1.159 0.573––2.344 0.411 0.681
Q 10.1246
DF 4
Significance level P = 0.0384
I2 (inconsistency) 60.49%
95% CI for I2 0.00–85.20

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.6 P = 0.0833 (low power); Egger: bias = −3.173487 (95% CI = −6.536026 to 0.189052) P = 0.575; Harbord–Egger: bias = −4.109035 (92.5% CI = −8.5417147 to 0.323644) P = 0.0889.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 11

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for TaqITT polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
Chang 2000 142/157 233/248 0.609 0.289–1.284
Balazs Gyorffy 2002 44/107 42/103 1.014 0.585–1.759
Vaselin Scrabic 2003 54/134 48/132 1.181 0.720–1.938
J J San Pedro 2005 24/71 31/88 0.939 0.486–1.813
Diego Garcia 2007 115/216 121/203 0.772 0.524–1.137
C Panierakis 2009 10/100 34/96 0.203 0.0933–0.440
Greear R M 2013 18/55 26/50 0.449 0.204–0.990
Chung Cheon 2015 66/81 105/113 0.335 0.135–0.834
Total (fixed effects) 473/921 640/1033 0.713 0.580–0.876 −3.213 0.001
Total (random effects) 473/921 640/1033 0.644 0.440–0.942 −2.270 0.023
Q 20.6282
DF 7
Significance level P = 0.0044
I2 (inconsistency) 66.07%
95% CI for I2 28.03–84.00

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.642857 P = 0.0141 (low power); Egger: bias = −3.773452 (95% CI = −8.197852 to 0.650947) P = 0.0819; Harbord–Egger: bias = −3.522136 (92.5% CI = −8.048273 to 1.004) P = 0.1452.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 12

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for TaqITt polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
Chang 2000 15/157 14/248 1.766 0.828–3.766
Balazs Gyorffy 2002 28/107 33/103 0.752 0.414–1.367
Vaselin Scrabic 2003 55/134 72/132 0.580 0.357–0.943
J J San Pedro 2005 36/71 43/88 1.076 0.576–2.012
Diego Garcia 2007 79/216 69/203 1.120 0.750–1.673
C Panierakis 2009 64/100 59/96 1.115 0.625–1.990
Greear R M 2013 26/55 24/50 0.971 0.451–2.091
Chung Cheon 2015 15/81 8/113 2.983 1.199–7.423
Total (fixed effects) 318/921 322/1033 1.017 0.829–1.248 0.165 0.869
Total (random effects) 318/921 322/1033 1.062 0.785–1.438 0.390 0.697
Q 13.8658
DF 7
Significance level P = 0.0536
I2 (inconsistency) 49.52%
95% CI for I2 0.00–77.47

Bias indicators: Begg–Mazumdar: Kendall’s Tau = 0.047619 P  > 0.9999 (low power); Egger: bias = 1.307065 (95% CI = −4.196619 to 6.810748) P = 0.5682; Harbord–Egger: bias = 1.410342 (92.5% CI = −3.284733 to 6.105417) P = 0.5305.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Table 13

P values testing Hardy–Weinberg proportion and test for heterogeneity of studies included in the meta-analysis for TaqItt polymorphism respectively.

Study Intervention Controls Odds ratio 95% CI z P
T J Chang 2000 0/157 1/248 0.524 0.0212–12.939
Balazs Gyorffy 2002 33/107 26/103 1.321 0.721–2.418
Vaselin Scrabic 2003 25/134 11/132 2.523 1.186–5.367
J J San Pedro 2005 11/71 14/88 0.969 0.410–2.290
Diego Garcia 2007 22/216 13/203 1.657 0.811–3.386
C Panierakis 2009 22/100 7/96 3.586 1.453–8.849
Greear R M 2013 6/55 5/50 1.102 0.314–3.862
Chung Cheon 2015 0/81 0/113
Total (fixed effects) 119/921 77/1033 1.677 1.225–2.295 3.230 0.001
Total (random effects) 119/921 77/1033 1.655 1.163–2.356 2.798 0.005
Q 6.9433
DF 6
Significance level P = 0.3261
I2 (inconsistency) 13.59%
95% CI for I2 0.00–75.19

Bias indicators: Begg–Mazumdar: Kendall’s Tau = −0.047619 P  > 0.9999 (low power); Egger: bias = −6.680722 (95% CI = −3.88097 to 2.519525) P = 0.608; Harbord–Egger: bias = −1.122847 (92.5% CI = −3.907446 to 1.661601) P = 0.4074.

DF, degree of freedom; Q, heterogeneity in meta analysis.

Of the articles included in the study, investigators of all studies included in the meta-analysis specifically looked for the presence of autoantibodies to diagnose type 1 diabetes and all fulfilled World Health Organization and the American Diabetes Association criteria (24).

Selection of controls varied across studies. Groups of controls included healthy blood donors and unrelated children.

ApaI, BsmI, FokI and TaqI polymorphisms and the risk for type 1 diabetes

For ApaI-AA, the odds ratio ranged from 0.476 to 2.051 (Table 2). The random-effects model yielded a pooled odds ratio of 1.081 (95 percent confidence interval (CI): 0.755–1.547). There was indication of heterogeneitiy (P = 0.0221).

For ApaI-Aa, the odds ratio ranged from 0.435 to 1.517 (Table 3). The random-effects model yielded a pooled odds ratio of 0.866 (95 percent confidence interval (CI): 0.664–1.131). There was indication of heterogeneity (P = 0.0469).

For ApaI-aa, the odds ratio ranged from 0.641 to 2.647 (Table 4). The fixed-effects model yielded a pooled odds ratio of 0.956 (95 percent confidence interval (CI): 0.772–1.182). There was indication of homogeneity (P = 0.1280). In view of these estimates, there is no evidence that any of the three alleles alone is associated with type 1 diabetes.

For BsmI-BB, the odds ratio ranged from 0.460 to 6.458 (Table 5). The fixed-effects model yielded a pooled odds ratio of 1.397 (95 percent confidence interval (CI): 1.034–1.888, P = 0.030). There was indication of homogeneity (P = 0.4531).

For BsmI-Bb, the odds ratio ranged from 0.598 to 5.210 (Table 6). The random-effects model yielded a pooled odds ratio of 1.534 (95 percent confidence interval (CI): 1.001–2.350, P = 0.049). There was indication of heterogeneity (P = 0.0014).

For BsmI-bb, the odds ratio ranged from 0.242 to 1.407 (Table 7). The random-effects model yielded a pooled odds ratio of 0.624 (95 percent confidence interval (CI): 0.418–0.933, P = 0.022). There was indication of heterogeneity (P = 0.0075).

Forest plots are shown in Fig. 2A, B and C for BsmI BB, Bb and bb alleles, respectively. Individual and pooled odds ratio estimates for the BsmI alleles are represented as squares and diamonds. In view of these estimates, there is evidence that BsmI-Bb or BsmI-bb alone is associated with type 1 diabetes.

For FokI-FF, the odds ratio ranged from 0.485 to 1.636 (Table 8). The fixed-effects model yielded a pooled odds ratio of 1.057 (95 percent confidence interval (CI): 0.817–1.367, P = 0.675). There was indication of homogeneity (P = 0.1443).

For FokI-Ff, the odds ratio ranged from 0.645 to 1.222 (Table 9). The fixed-effects model yielded a pooled odds ratio of 0.724 (95 percent confidence interval (CI): 0.564–0.929, P = 0.011). There was indication of homogeneity (P = 0.4324).

For FokI-ff, the odds ratio ranged from 0.128 to 2.558 (Table 10). The random-effects model yielded a pooled odds ratio of 1.159 (95 percent confidence interval (CI): 0.573–2.344, P = 0.681). There was indication of heterogeneity (P = 0.0384).

For TagI-TT, the odds ratio ranged from 0.203 to 1.181 (Table 11). The random-effects model yielded a pooled odds ratio of 0.644 (95 percent confidence interval (CI): 0.440–0.942, P = 0.023). There was indication of heterogeneity (P = 0.0044).

For TaqI-Tt, the odds ratio ranged from 0.580 to 2.983 (Table 12). The random-effects model yielded a pooled odds ratio of 1.062 (95 percent confidence interval (CI): 0.785–1.438, P = 0.697). There was some indication of heterogeneity (P = 0.0536).

For TaqI-tt, the odds ratio ranged from 0.524 to 3.586 (Table 13). The fixed-effects model yielded a pooled odds ratio of 1.655 (95 percent confidence interval (CI): 1.677–2.295, P = 0.001). There was indication of heterogeneity (P = 0.3261).

Forest plots are shown in Fig. 3 and B for TaqI-TT and tt alleles, respectively. Individual and pooled odds ratio estimates for the TaqI alleles are represented as squares and diamonds. In view of these estimates, there is evidence that TaqI-TT and TaqI-tt alone is associated with type 1 diabetes.

Discussion

There are a number of reports on FokI, ApaI, TaqI and BsmI polymorphisms of the VDR gene in diabetic patients, there have not been conclusive evidence that any of these polymorphisms has causative association with type 1 DM in children. In a 2006 meta-analysis that focused on VDR polymorphisms, FokI, ApaI, TaqI, BsmI and type 1 DM association included mainly adult samples and did not reveal any specific association (25). However, out of 19 published papers included in this meta-analysis, authors of only five papers specifically looked for the presence of autoantibodies to distinguish type 1 diabetes from type 2. In other 14 studies included in this meta-analysis, investigators used only one criteria (e.g., ketosis, early requirement of insulin). This may be one of the main reasons for different statistical results other than ethnic diversities when compared to our results. DM1 is mainly a disease of pediatric age: considering the qualitative assessment of study inclusion criteria, choosing studies where diagnosis is at pediatric age with age matching control samples, and/or taking American Diabetes Association criteria would end up with more reliable meta-analysis results. In our study, mean age of control samples are in pediatric range. In another meta-analysis involving Chinese adult samples, authors concluded that BsmI polymorphisms in the VDR region would increase the risk of type 1 DM in East Asians (26). In the study of Zhang J, Asian samples with BsmI polymorphism was found to have a significant association with increased risk of type 1 DM (27). The study of Qin WH demonstrated a significant relationship among BsmI B allele and BB genotype and increased risk for type 1 DM in Asians, whereas this study included Latino and African adult samples and authors also found another specific association with BsmIbb genotype and type 1 DM in overall populations (28). The novel finding in our study was the presence of an increased risk of type 1 DM in carriers of BsmIBB, BsmIBb and TaqItt polymorphisms and decreased risk of type 1 DM in children with BsmIbb and TaqITT polymorphisms. There are GWAS studies that widen our approach to vitamin D receptor polymorphisms and DM1.

Meta-analysis may be more reliable when evaluating genotype frequencies in certain diseases because in a way it may reduce the effect of biased sampling or nonrandom mating in individual study population. Results of studies so far, regarding VDR polymorphisms and DM1 susceptibility, are conflicting. In the study of Garcia and coworkers, an association was found between BsmI polymorphism and DM1 (15). The frequency of genotype bb was found to be significantly lower in the cases than that in controls. Among five prevalent haplotypes, BAT has been found to be statistically more frequent in study group in the same study. Among genotype combinations, AabbTT was found to be higher in controls. In our study population, genotype combination frequencies were not assessed because of unavailable data, and this may be one of the limitations of our study.

In the study of San-Pedro JI, an association of an haplotype ‘fBAt’ and risk of type 1 DM in Basque population has been identified (16). In the study of Skrabic V, BBAAtt genotype combination was found to be associated with type 1 DM in Dalmatian population of southern Croatia, with the ‘tt’ genotype preferentially presented in the affected individuals (17). This was also noticed in previous studies that focused on association of VDR gene polymorphisms with increased susceptibility to T1DM in Taiwanese, Japanese, South Asian (Indian) and German populations (29, 30, 31, 32). TaqI polymorphism among type 1 DM patients and control subjects differed significantly, with the VDR tt genotype occurring more frequently in T1DM patients. No difference was noticed in the genotype frequencies of BsmI and ApaI polymorphisms in cases and controls. We evaluated TaqItt polymorphism frequency and demonstrated a significant increase in diabetic children in our study. In the study of Zemunik and coworkers, some evidence of association of Tru91–BsmI haplotype and type 1 DM in population of South Croatia was found (17). In our meta-analysis, we have included two studies from Croatia. One of the limitations of this study was that its sample size was small, it only included nine studies and the power of this study is not high. In the study of Panierakis and coworkers, homogeneous southern European population with low incidence of type 1 DM was included in the study group, and they found an association of T1DM and FokI, BsmI, ApaI and TaqI polymorphisms. In this study, FokIFF genotype and F allele and BsmIBB genotype and B allele were less frequent in individuals with T1DM (21). In the same study, ApaIAA genotype and A allele, as well as TaqITT genotype and T allele were more frequent in individuals. Greear and coworkers also studied the association of Taq1, Fok1, Apa1 and type 1 DM and found no significant difference in distribution of VDR polymorphisms in diabetic patients, whereas diabetic patients had significantly decreased levels of vitamin D levels than healthy controls (22). In the study of Cheon CK and coworkers, the frequency of bb and TT genotype has been found to be significantly increased among carriers demonstrating a protective effect in Korean subjects (23). Gyorffy and coworkers suggested a strong linkage disequilibrium between the ‘b’ and ‘a’ alleles in his study. The close loci of these polymorphisms might be an explanation for the stability of linkage, but the background of these combinations need further investigation (19). In the same study, a strong association has been found between carrier state of the ‘b’ + ‘a’ + ‘u’ alleles and the presence of type 1 DM in females. There are other reports as well that point out to a gender-specific association and consequence of gene polymorphism (33). A number of studies have shown that patients with 1 DM have low levels of vitamin D although other studies have conflicting results (34, 35, 36). In 2010, GWAS study in approximately 30,000 individuals from European descent identifies variants at four loci that were associated with 25(OH)D levels: GC rs2282679, dhcr7 RS 12785878 and CYP2R1 (37). A second GWAS of 25(OH)D levels confirmed the findings with GC, DHCR7 and CYP2R1 (38). These variants are located within or near genes involved in vitamin D transport (GC), cholesterol synthesis (DHCR7) and hydroxylation (CYP2R1 and CTP24A1) (37). Cooper and coworkers recently tested genetic variants influencing 25(OH)D metabolism for an association with both circulating 25(OH)D metabolism for an association with both circulating 25(OH)D concentrations and T1D. They replicated the associations found in the GWAS of the four vitamin D metabolism genes (GC, DHCR7, CYP2R1 and CYP24A1) with 25(OH)D in control subjects and found that CYP27B1, DHCR7 and CYP2R1 were associated with type 1 diabetes (39). The Fok1 polymorphism of the VDR, which increases the transcriptional activity of VDR, has been suggested as an influencing factor for susceptibility to T1DM. It affects insulin secretion and sensitivity and has been found to be a susceptibility factor for the development of diabetic retinopathy (40, 41). In addition, vitamin D-binding protein gene polymorphisms were found to be associated with diabetes-associated antibody insulinoma antigen 2 and with T1DM (42). In the study of Greear and coworkers, low vitamin D levels in the diabetic children have been attributed to inflammatory or other pathologic processes, mainly as a consequence of the disease rather than being a risk factor, as previously stated in DAISY study (43). In the study of Chang and coworkers, the allele frequency of the BsmI differed between Taiwanese patients and controls significantly (20). There are some limitations in this meta-analysis. The power of this study should further be increased by additional studies, and this meta-analysis involves only nine studies. Some of the studies contained small number of cases, and background of the patients varied across included studies. However our meta-analysis employed a random-effects model designed to encounter these variations and found significant effect of polymorphisms on type 1 DM susceptibility.

As a conclusion, our meta-analysis of accessible published data has demonstrated statistically significant association between BsmIBB, BsmIBb, BsmIbb, TaqITT and Taqtt polymorphisms and susceptibility to type 1 DM in children; however, influence of vitamin D receptor gene polymorphisms on susceptibility to type 1 diabetes deserves further investigations. Meta-analysis includes larger data sets and accordingly may demonstrate more reliable statistical results to rule out genotype–phenotype correlations of diseases.

Declaration of interest

The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.

Funding

This work did not receive any specific grant from any funding agency in the public, commercial, or not-for-profit sector.

References

  • 1

    Todd JA, Bell JI & McDevitt HO. HLA DQ beta gene contributes to susceptibility and resistence to insulin dependent diabetes mellitus. Nature 1987 329 599604. (doi:10.1038/329599a0)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 2

    Ban Y, Taniyama M, Yanagawa T, Yamada S, Maruyama T, Kasuga A & Ban Y. Vitamin D receptor initiation codon polymorphism influences susceptibility to type 1 dibetes mellitus in Japanese population. BMC Medical Genetics 2001 2 713. (doi:10.1186/1471-2350-2-7)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 3

    Lemire JM, Adams JS, Sakai R & Jordan SC. 1alpha,25-Dihydroxyvitamin D3 suppresses proliferation and immunoglobulin production by normal human peripheral blood mononuclear cells. Journal of Clinical Investigation 1984 74 657661. (doi:10.1172/JCI111465)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 4

    Yang S, Smith C, Prahl JM, Luo X & DeLuca HF. Vitamin D deficiency suppresses cell-mediated immunity in vivo. Archives of Biochemistry and Biophysics 1993 1 98106. (doi:10.1006/abbi.1993.1260)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 5

    Pugliese A. Genetics of type 1 diabetes. Endocrinology and Metabolism Clinics of North America 2004 33 116. (doi:10.1016/S0889-8529(03)00082-3)

  • 6

    Vyse TJ & Todd JA. Genetic analysis of autoimmune disease. Cell 1996 85 311318. (doi:10.1016/S0092-8674(00)81110-1)

  • 7

    McDermott MF, Ramachandran A, Ogunkolade BW, Aganna E, Curtis D, Boucher BJ, Snehalatha C & Hitman GA. Allelic variation in the vitamin D receptor influences susceptibility to IDDM in Indian Asians. Diabetologia 1997 8 971975. (doi:10.1007/s001250050776)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 8

    Motohashi Y, Yamada S & Yanagawa T. Vitamin D receptor gene polymorphism affects onset pattern of type 1 diabetes. Journal of Clinical Endocrinology and Metabolism 2003 88 31373140. (doi:10.1210/jc.2002-021881)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 9

    Turpeinen H, Hermann R, Vaara S, Laine AP, Simell O, Knip M, Veijola R & Ilonen J. Vitamin D receptor polymorphisms: no association with type 1 diabetes in the Finnish population. European Journal of Endocrinology 2003 6 591596. (doi:10.1530/eje.0.1490591)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 10

    Lemos MC, Fagulha A, Coutinho E, Gomes L, Bastos M, Barros L, Carrilho F, Geraldes E, Regateiro FJ & Carvalheiro M. Lack of association of vitamin D receptor polymorphisms with susceptibility to type 1 diabetes mellitus in the Portugeese polpulation. Human Immunology 2008 2 134138. (doi:10.1016/j.humimm.2008.01.008)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 11

    Obi-Tabot ET, Tian XQ, Chen TC & Holick MF. A human skin equivalent model that mimics the photoproduction of vitamin D3 in human skin. In Vitro Cellular and Developmental Biology 2000 3 201204. (doi:10.1290/1071-2690(2000)036<0201:AHSEMT>2.0.CO;2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 12

    San-Pedro JI, Bilbao JR, Perez de Nanclares G. Heterogeneity of vitamin D receptor gene association with celiac disease and type1 diabetes mellitus, Autoimmunity 2005 38 439444. (doi:10.1080/08916930500288455)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 13

    Nejentsev S, Cooper JD, Godfrey L, Howson JM, Rance H, Nutland S, Walker NM, Guja C, Ionescu-Tirgovişte C & Savage DA et al. Analysis of the vitamin D receptor gene sequence variants in type 1 diabetes. Diabetes 2004 10 641642.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 14

    DerSimonian R & Laird N. Meta-analysis in clinical trials. Controlled Clinical Trials 1986 3 177188. (doi:10.1016/0197-2456(86)90046-2)

  • 15

    Garcia D, Angel B, Carrasco E, Albala C, Santos JL & Pérez-Bravo F. VDR polymorphisms influence the immune response in type 1 diabetic children from Santiago, Chile. Diabetes Research and Clinical Practice 2007 1 134140. (doi:10.1016/j.diabres.2006.10.018)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 16

    San Pedro JI, Bilbao JR, Perez de Nanclares G, Vitoria JC, Martul P, Castaño L.Heterogeneity of vitamin D receptor gene association with celiac disease and type 1 diabetes mellitus. Autoimmunity 2005 6 439444. (doi:10.1080/08916930500288455)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 17

    Zemunik T, Scrabic V, Boraska V, Diklic D, Terzic IM, Capkun V, Peruzovic M & Terzic J. FokI polymorphism, vitamin D receptor, and interleukin-1 receptor haplotypes are associated with type 1 diabetes in the dalmatian population. Journal of Molecular Diagnostics 2005 5 600604. (doi:10.1016/S1525-1578(10)60593-4)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 18

    Skrabic V, Zemunik T, Situm M & Terzić J. Vitamin D receptor polymorphism and susceptibility to type 1 diabetes in the Dalmatian population. Diabetes Research and Clinical Practice 2003 1 3135. (doi:10.1016/S0168-8227(02)00195-X)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 19

    Gyorffy B, Vasarhelyi B, Krikovszky D, Madácsy L, Tordai A, Tulassay T & Szabó A. Gender-spesific association of vitamin D receptor polymorphism combinations with type 1 diabetes mellitus. European Journal of Endocrinology 2002 6 803808. (doi:10.1530/eje.0.1470803)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 20

    Chang TJ, Lei HH, Yeh JI, Chiu KC, Lee KC, Chen MC, Tai TY & Chuang LM. Vitamin D receptor gene polymorphisms influence susceptibility to type 1 diabetes mellitus in the Taiwanese population. Clinical Endocrinology 2000 5 575580. (doi:10.1046/j.1365-2265.2000.00985.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 21

    Panierakis C, Goulielmos G, Mamoulakis D, Petraki E, Papavasiliou E & Galanakis E. Vitamin D receptor gene polymorphisms and susceptibility to type 1 diabetes in Crete, Greece. Clinical Immunology 2009 2 276281. (doi:10.1016/j.clim.2009.08.004)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 22

    Greear RM, Portelli SL, Hung BS, Cleghorn GJ, McMahon SK, Batch JA & Conwell LS. Serum vitamin D levels are lower in Australian children and adolescents with type 1 diabetes than in children without diabetes. Pediatric Diabetes 2013 1 3141. (doi:10.1111/j.1399-5448.2012.00890.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 23

    Cheon CK, Nam HK, Lee KH, Kim SY, Song JS & Kim C. Vitamin D receptor gene polymorphisms and type 1 diabetes mellitus in a Korean population. Pediatrics International 2015 5 870874. (doi:10.1111/ped.12634)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 24

    American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care 2010 33 (Supplement) 62S69S. (doi:10.2337/dc11-s062)

  • 25

    Guo SW, Magnuson VL, Schiller JJ, Wang X, Wu Y & Ghosh S. Meta-analysis of vitamin D receptor polymorphisms and type 1 diabetes: a HuGE Review of Genetic Association Studies. American Journal of Epidemiology 2006 8 711724. (doi:10.1093/aje/kwj278)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 26

    Wang Q, Xi B, Reilly KH, Liu M, Fu M. Quantitative assessment of the associations between four polymorphisms (FokI, ApaI, BsmaI, TaqI) of vitamin D receptor gene and risk of diabetes mellitus. Molecular Biology Reports 2012 10 94059414. (doi:10.1007/s11033-012-1805-7)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 27

    Zhang J, Li W, Liu J, Wu W, Ouyang H, Zhang Q, Wang Y, Liu L, Yang R & Liu X et al. Polymorphisms in the vitamin D receptor gene and type 1 diabetes mellitus risk: an update by meta-analysis. Molecular and Cellular Endocrinology 2012 1 135142. (doi:10.1016/j.mce.2012.02.003)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 28

    Qin WH, Wang XH, Qiu LJ, Huang XB, Huang Y, Wu NR & Liang HS. A meta-analysis of association of vitamin D receptor BsmI gene polymorphisms with the risk of type 1 diabetes mellitus. Journal of Receptor and Signal Transduction Research 2014 5 372377. (doi:10.3109/10799893.2014.903420)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 29

    Yokota I, Satomura S, Kitamura S, Taki Y, Naito E, Ito M, Nisisho K & Kuroda Y. Association between vitamin D receptor genotype and age of onset in juvenile Japanese patients with type 1 diabetes. Diabetes Care 2002 7 1244. (doi:10.2337/diacare.25.7.1244)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 30

    McDermott MF, Ramachandran A, Ogunkolde BW, Aganna E, Curtis D, Boucher BJ, Snehalatha C & Hitman GA. Allelic variations in the vitamin D receptor influences susceptibility to IDDM in Indian Asians. Diabetologia 1997 8 971975. (doi:10.1007/s001250050776)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 31

    Pani MA, Knapp M, Donner H, Braun J, Baur MP, Usadel KH & Badenhoop K. Vitamin D receptor allele combinations influence genetic susceptibility to type I Diabetes in Germans. Diabetes 2000 3 504507. (doi:10.2337/diabetes.49.3.504)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 32

    Ban Y, Taniyama M, Yanagawa T, Yamada S, Maruyama T, Kasuga A & Ban Y. Vitamin D receptor initiation codon polymorphism influences genetic susceptibility to type 1 diabetes mellitus in the Japanese population. BMC Medical Genetics 2001 7 27.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 33

    Kitagawa I, Kitagawa Y, Kawase Y, Nagaya T & Tokudome S. Advanced onset of menarche and higher bone mineral density depending on itamin D receptor gene polymorphism. European Journal of Endocrinology 1998 139 522527. (doi:10.1530/eje.0.1390522)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 34

    Baumgartl HJ, Standl E, Schmidt-Gayk H, Kolb HJ, Janka HU & Ziegler AG. Changes of vitamin D3 serum concentrations at the onset of immune-mediated type 1 (insulin dependent) diabetes mellitus. Diabetes Research 1991 16 145148.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 35

    Holick MF. Vitamin D: importance in the prevention of cancers, type 1 diabetes, heart disease and osteoporosis. American Journal of Clinical Nutrition 2004 37 680683.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 36

    Pozzili P, Manfrini S, Crino A, Picardi A, Leomanni C, Cherubini V, Valente L, Khazrai M, Visalli NIMDIAB Group. Low levels of 25-hydroxyvitamin D3 and 1,25-dihydroxyvitamin D3 in patients with nely diagnosed type 1 diabetes. Hormone and Metabolic Research 2005 3 680683. (doi:10.1055/s-2005-870578)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 37

    Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, Kiel DP, Streeten EA, Ohlsson C & Koller DL et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet 2010 376 180188. (doi:10.1016/S0140-6736(10)60588-0)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 38

    Ahn J, Yu K, Stolzenberg-Solomon R, Simon KC, McCullough ML, Gallicchio L, Jacobs EJ, Ascherio A, Helzlsouer K & Jacobs KB et al. Genome-wide association study of circulating vitamin D levels. Human Molecular Genetics 2010 19 27392745. (doi:10.1093/hmg/ddq155)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 39

    Cooper JD, Smyth DJ, Walker NM, Stevens H, Burren OS, Wallace C, Greissl C, Ramos-Lopez E, Hyppönen E & Dunger DB et al. Inherited variation in vitamin D genes is associated with predisposition to autoimmune disease type 1 diabetes. Diabetes 2011 60 16241631. (doi:10.2337/db10-1656)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 40

    Taverna MJ, Selam JL & Slama G. Association between a protein polymorphism in the start codon of the vitamin D receptor gene and severe diabetic retinopathy in C-peptide-negative type 1 diabetes. Journal of Clinical Endocrinology and Metabolism 2005 90 48034808. (doi:10.1210/jc.2004-2407)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 41

    Taverna MJ, Sola A, Guyot-Argenton C, Pacher N, Bruzzo F, Slama G, Reach G & Selam JL. TaqI polymorphism of the vitamin D receptor and risk of severe diabetic retinopathy. Diabetologia 2002 45 436442. (doi:10.1007/s00125-001-0769-2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 42

    Ongagna JC, Pinget M, Belcourt A. Vitamin D-binding protein gene polymorphism association with IA-2 utoantibodies in type 1 diabetes. Clinical Biochemistry 2005 38 415419. (doi:10.1016/j.clinbiochem.2004.12.013)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 43

    Simpson M, Brady H, Yin X, Seifert J, Barriga K, Hoffman M, Bugawan T, Barón AE, Sokol RJ & Eisenbarth G et al. No association of vitamin D intake or 25-hydroxyvitamin D levels in childhood with risk of islet autoimmunity and type 1 diabetes: the Daibetes Autoimmunity Study in the Young (DAISY). Diabetologia 2011 54 27792788. (doi:10.1007/s00125-011-2278-2)

    • PubMed
    • Search Google Scholar
    • Export Citation

 

  • Collapse
  • Expand
  • Flow chart for identification and selection of studies.

  • (A, B and C) Forest plots showing individual and pooled odds ratio estimates of BsmI BB, BsmI Bb, BsmIbb polymorphisms, respectively.

  • (A and B) Forest plots showing individual and pooled odds ratios of TaqITT and TaqItt polymorphisms, respectively.

  • 1

    Todd JA, Bell JI & McDevitt HO. HLA DQ beta gene contributes to susceptibility and resistence to insulin dependent diabetes mellitus. Nature 1987 329 599604. (doi:10.1038/329599a0)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 2

    Ban Y, Taniyama M, Yanagawa T, Yamada S, Maruyama T, Kasuga A & Ban Y. Vitamin D receptor initiation codon polymorphism influences susceptibility to type 1 dibetes mellitus in Japanese population. BMC Medical Genetics 2001 2 713. (doi:10.1186/1471-2350-2-7)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 3

    Lemire JM, Adams JS, Sakai R & Jordan SC. 1alpha,25-Dihydroxyvitamin D3 suppresses proliferation and immunoglobulin production by normal human peripheral blood mononuclear cells. Journal of Clinical Investigation 1984 74 657661. (doi:10.1172/JCI111465)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 4

    Yang S, Smith C, Prahl JM, Luo X & DeLuca HF. Vitamin D deficiency suppresses cell-mediated immunity in vivo. Archives of Biochemistry and Biophysics 1993 1 98106. (doi:10.1006/abbi.1993.1260)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 5

    Pugliese A. Genetics of type 1 diabetes. Endocrinology and Metabolism Clinics of North America 2004 33 116. (doi:10.1016/S0889-8529(03)00082-3)

  • 6

    Vyse TJ & Todd JA. Genetic analysis of autoimmune disease. Cell 1996 85 311318. (doi:10.1016/S0092-8674(00)81110-1)

  • 7

    McDermott MF, Ramachandran A, Ogunkolade BW, Aganna E, Curtis D, Boucher BJ, Snehalatha C & Hitman GA. Allelic variation in the vitamin D receptor influences susceptibility to IDDM in Indian Asians. Diabetologia 1997 8 971975. (doi:10.1007/s001250050776)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 8

    Motohashi Y, Yamada S & Yanagawa T. Vitamin D receptor gene polymorphism affects onset pattern of type 1 diabetes. Journal of Clinical Endocrinology and Metabolism 2003 88 31373140. (doi:10.1210/jc.2002-021881)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 9

    Turpeinen H, Hermann R, Vaara S, Laine AP, Simell O, Knip M, Veijola R & Ilonen J. Vitamin D receptor polymorphisms: no association with type 1 diabetes in the Finnish population. European Journal of Endocrinology 2003 6 591596. (doi:10.1530/eje.0.1490591)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 10

    Lemos MC, Fagulha A, Coutinho E, Gomes L, Bastos M, Barros L, Carrilho F, Geraldes E, Regateiro FJ & Carvalheiro M. Lack of association of vitamin D receptor polymorphisms with susceptibility to type 1 diabetes mellitus in the Portugeese polpulation. Human Immunology 2008 2 134138. (doi:10.1016/j.humimm.2008.01.008)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 11

    Obi-Tabot ET, Tian XQ, Chen TC & Holick MF. A human skin equivalent model that mimics the photoproduction of vitamin D3 in human skin. In Vitro Cellular and Developmental Biology 2000 3 201204. (doi:10.1290/1071-2690(2000)036<0201:AHSEMT>2.0.CO;2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 12

    San-Pedro JI, Bilbao JR, Perez de Nanclares G. Heterogeneity of vitamin D receptor gene association with celiac disease and type1 diabetes mellitus, Autoimmunity 2005 38 439444. (doi:10.1080/08916930500288455)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 13

    Nejentsev S, Cooper JD, Godfrey L, Howson JM, Rance H, Nutland S, Walker NM, Guja C, Ionescu-Tirgovişte C & Savage DA et al. Analysis of the vitamin D receptor gene sequence variants in type 1 diabetes. Diabetes 2004 10 641642.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 14

    DerSimonian R & Laird N. Meta-analysis in clinical trials. Controlled Clinical Trials 1986 3 177188. (doi:10.1016/0197-2456(86)90046-2)

  • 15

    Garcia D, Angel B, Carrasco E, Albala C, Santos JL & Pérez-Bravo F. VDR polymorphisms influence the immune response in type 1 diabetic children from Santiago, Chile. Diabetes Research and Clinical Practice 2007 1 134140. (doi:10.1016/j.diabres.2006.10.018)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 16

    San Pedro JI, Bilbao JR, Perez de Nanclares G, Vitoria JC, Martul P, Castaño L.Heterogeneity of vitamin D receptor gene association with celiac disease and type 1 diabetes mellitus. Autoimmunity 2005 6 439444. (doi:10.1080/08916930500288455)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 17

    Zemunik T, Scrabic V, Boraska V, Diklic D, Terzic IM, Capkun V, Peruzovic M & Terzic J. FokI polymorphism, vitamin D receptor, and interleukin-1 receptor haplotypes are associated with type 1 diabetes in the dalmatian population. Journal of Molecular Diagnostics 2005 5 600604. (doi:10.1016/S1525-1578(10)60593-4)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 18

    Skrabic V, Zemunik T, Situm M & Terzić J. Vitamin D receptor polymorphism and susceptibility to type 1 diabetes in the Dalmatian population. Diabetes Research and Clinical Practice 2003 1 3135. (doi:10.1016/S0168-8227(02)00195-X)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 19

    Gyorffy B, Vasarhelyi B, Krikovszky D, Madácsy L, Tordai A, Tulassay T & Szabó A. Gender-spesific association of vitamin D receptor polymorphism combinations with type 1 diabetes mellitus. European Journal of Endocrinology 2002 6 803808. (doi:10.1530/eje.0.1470803)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 20

    Chang TJ, Lei HH, Yeh JI, Chiu KC, Lee KC, Chen MC, Tai TY & Chuang LM. Vitamin D receptor gene polymorphisms influence susceptibility to type 1 diabetes mellitus in the Taiwanese population. Clinical Endocrinology 2000 5 575580. (doi:10.1046/j.1365-2265.2000.00985.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 21

    Panierakis C, Goulielmos G, Mamoulakis D, Petraki E, Papavasiliou E & Galanakis E. Vitamin D receptor gene polymorphisms and susceptibility to type 1 diabetes in Crete, Greece. Clinical Immunology 2009 2 276281. (doi:10.1016/j.clim.2009.08.004)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 22

    Greear RM, Portelli SL, Hung BS, Cleghorn GJ, McMahon SK, Batch JA & Conwell LS. Serum vitamin D levels are lower in Australian children and adolescents with type 1 diabetes than in children without diabetes. Pediatric Diabetes 2013 1 3141. (doi:10.1111/j.1399-5448.2012.00890.x)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 23

    Cheon CK, Nam HK, Lee KH, Kim SY, Song JS & Kim C. Vitamin D receptor gene polymorphisms and type 1 diabetes mellitus in a Korean population. Pediatrics International 2015 5 870874. (doi:10.1111/ped.12634)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 24

    American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care 2010 33 (Supplement) 62S69S. (doi:10.2337/dc11-s062)

  • 25

    Guo SW, Magnuson VL, Schiller JJ, Wang X, Wu Y & Ghosh S. Meta-analysis of vitamin D receptor polymorphisms and type 1 diabetes: a HuGE Review of Genetic Association Studies. American Journal of Epidemiology 2006 8 711724. (doi:10.1093/aje/kwj278)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 26

    Wang Q, Xi B, Reilly KH, Liu M, Fu M. Quantitative assessment of the associations between four polymorphisms (FokI, ApaI, BsmaI, TaqI) of vitamin D receptor gene and risk of diabetes mellitus. Molecular Biology Reports 2012 10 94059414. (doi:10.1007/s11033-012-1805-7)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 27

    Zhang J, Li W, Liu J, Wu W, Ouyang H, Zhang Q, Wang Y, Liu L, Yang R & Liu X et al. Polymorphisms in the vitamin D receptor gene and type 1 diabetes mellitus risk: an update by meta-analysis. Molecular and Cellular Endocrinology 2012 1 135142. (doi:10.1016/j.mce.2012.02.003)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 28

    Qin WH, Wang XH, Qiu LJ, Huang XB, Huang Y, Wu NR & Liang HS. A meta-analysis of association of vitamin D receptor BsmI gene polymorphisms with the risk of type 1 diabetes mellitus. Journal of Receptor and Signal Transduction Research 2014 5 372377. (doi:10.3109/10799893.2014.903420)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 29

    Yokota I, Satomura S, Kitamura S, Taki Y, Naito E, Ito M, Nisisho K & Kuroda Y. Association between vitamin D receptor genotype and age of onset in juvenile Japanese patients with type 1 diabetes. Diabetes Care 2002 7 1244. (doi:10.2337/diacare.25.7.1244)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 30

    McDermott MF, Ramachandran A, Ogunkolde BW, Aganna E, Curtis D, Boucher BJ, Snehalatha C & Hitman GA. Allelic variations in the vitamin D receptor influences susceptibility to IDDM in Indian Asians. Diabetologia 1997 8 971975. (doi:10.1007/s001250050776)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 31

    Pani MA, Knapp M, Donner H, Braun J, Baur MP, Usadel KH & Badenhoop K. Vitamin D receptor allele combinations influence genetic susceptibility to type I Diabetes in Germans. Diabetes 2000 3 504507. (doi:10.2337/diabetes.49.3.504)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 32

    Ban Y, Taniyama M, Yanagawa T, Yamada S, Maruyama T, Kasuga A & Ban Y. Vitamin D receptor initiation codon polymorphism influences genetic susceptibility to type 1 diabetes mellitus in the Japanese population. BMC Medical Genetics 2001 7 27.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 33

    Kitagawa I, Kitagawa Y, Kawase Y, Nagaya T & Tokudome S. Advanced onset of menarche and higher bone mineral density depending on itamin D receptor gene polymorphism. European Journal of Endocrinology 1998 139 522527. (doi:10.1530/eje.0.1390522)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 34

    Baumgartl HJ, Standl E, Schmidt-Gayk H, Kolb HJ, Janka HU & Ziegler AG. Changes of vitamin D3 serum concentrations at the onset of immune-mediated type 1 (insulin dependent) diabetes mellitus. Diabetes Research 1991 16 145148.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 35

    Holick MF. Vitamin D: importance in the prevention of cancers, type 1 diabetes, heart disease and osteoporosis. American Journal of Clinical Nutrition 2004 37 680683.

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 36

    Pozzili P, Manfrini S, Crino A, Picardi A, Leomanni C, Cherubini V, Valente L, Khazrai M, Visalli NIMDIAB Group. Low levels of 25-hydroxyvitamin D3 and 1,25-dihydroxyvitamin D3 in patients with nely diagnosed type 1 diabetes. Hormone and Metabolic Research 2005 3 680683. (doi:10.1055/s-2005-870578)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 37

    Wang TJ, Zhang F, Richards JB, Kestenbaum B, van Meurs JB, Berry D, Kiel DP, Streeten EA, Ohlsson C & Koller DL et al. Common genetic determinants of vitamin D insufficiency: a genome-wide association study. Lancet 2010 376 180188. (doi:10.1016/S0140-6736(10)60588-0)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 38

    Ahn J, Yu K, Stolzenberg-Solomon R, Simon KC, McCullough ML, Gallicchio L, Jacobs EJ, Ascherio A, Helzlsouer K & Jacobs KB et al. Genome-wide association study of circulating vitamin D levels. Human Molecular Genetics 2010 19 27392745. (doi:10.1093/hmg/ddq155)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 39

    Cooper JD, Smyth DJ, Walker NM, Stevens H, Burren OS, Wallace C, Greissl C, Ramos-Lopez E, Hyppönen E & Dunger DB et al. Inherited variation in vitamin D genes is associated with predisposition to autoimmune disease type 1 diabetes. Diabetes 2011 60 16241631. (doi:10.2337/db10-1656)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 40

    Taverna MJ, Selam JL & Slama G. Association between a protein polymorphism in the start codon of the vitamin D receptor gene and severe diabetic retinopathy in C-peptide-negative type 1 diabetes. Journal of Clinical Endocrinology and Metabolism 2005 90 48034808. (doi:10.1210/jc.2004-2407)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 41

    Taverna MJ, Sola A, Guyot-Argenton C, Pacher N, Bruzzo F, Slama G, Reach G & Selam JL. TaqI polymorphism of the vitamin D receptor and risk of severe diabetic retinopathy. Diabetologia 2002 45 436442. (doi:10.1007/s00125-001-0769-2)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 42

    Ongagna JC, Pinget M, Belcourt A. Vitamin D-binding protein gene polymorphism association with IA-2 utoantibodies in type 1 diabetes. Clinical Biochemistry 2005 38 415419. (doi:10.1016/j.clinbiochem.2004.12.013)

    • PubMed
    • Search Google Scholar
    • Export Citation
  • 43

    Simpson M, Brady H, Yin X, Seifert J, Barriga K, Hoffman M, Bugawan T, Barón AE, Sokol RJ & Eisenbarth G et al. No association of vitamin D intake or 25-hydroxyvitamin D levels in childhood with risk of islet autoimmunity and type 1 diabetes: the Daibetes Autoimmunity Study in the Young (DAISY). Diabetologia 2011 54 27792788. (doi:10.1007/s00125-011-2278-2)

    • PubMed
    • Search Google Scholar
    • Export Citation