BioMed Research International

BioMed Research International / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 8197580 | 9 pages |

Vitamin D Receptor Gene Polymorphisms and Autoimmune Thyroiditis: Are They Associated with Disease Occurrence and Its Features?

Academic Editor: Salvatore Gallone
Received04 Apr 2019
Revised15 Jul 2019
Accepted28 Jul 2019
Published21 Aug 2019


Purpose. Vitamin D, besides its role in calcium-phosphorus metabolism, turned out to play a significant immunomodulating function. Until now four single nucleotide polymorphisms of vitamin D receptor gene (VDR), rs2228570 (FokI), rs1544410 (BsmI), rs7975232 (ApaI), and rs731236 (TaqI), have been studied in autoimmune thyroid disorders, with conflicting results. Another functional polymorphism of the VDR gene, rs11568820 (Cdx2), has been shown to influence the immune system, although it has not been studied for its association with autoimmune thyroiditis to date. Therefore, the study aimed to evaluate the association of these five VDR gene polymorphisms with susceptibility to autoimmune thyroiditis among Caucasian Polish population. A relationship between the studied polymorphisms and selected clinical features of the disease was additionally assessed. Methods. 223 patients with autoimmune thyroiditis and 130 control subjects were enrolled in the study. VDR polymorphisms were studied by PCR-RFLP or TaqMan real-time PCR. Results. Allele and genotype distributions of any of the studied polymorphisms did not differ significantly between patients and controls. Similarly, frequencies of haplotypes derived from rs1544410-rs7975232-rs731236 (BsmI-ApaI-TaqI) polymorphisms were not significantly different in the two studied groups. However, a weak association between rs1544410 (BsmI) or rs7975232 (ApaI) VDR polymorphisms and thyroid volume was found (p = 0.03 and p = 0.04, resp.). Conclusions. Our results suggest that VDR gene is not a major susceptibility factor for autoimmune thyroiditis development, at least in Caucasian Polish population.

1. Introduction

Autoimmune thyroid diseases (AITDs) are common pathologies that affect up to 5% of the population and their incidence is still increasing [1]. Autoimmune thyroiditis (AIT), also known as Hashimoto’s thyroiditis, is the most frequent manifestation of AITDs. The etiology of the disease is complex. It develops as a result of an interaction between predisposing gene variants and environmental triggers. Among the known genetic factors there are HLA class II gene variants and polymorphisms of non-HLA immune-regulating genes (e.g., CTLA-4, FOXP3, and PTPN22) or thyroid-specific genes (TPO) [2]. According to twin studies, genetic factors in general predominate and contribute 70-80 % to the development of AIT. However, these already proven genetic agents are responsible for less than 20% of the genetic susceptibility to the disease [3].

The active form of vitamin D was found to exert an immunomodulating effect. Because of that 1,25(OH)2D may be a meaningful player in the pathogenesis of autoimmune and inflammatory disorders [4, 5]. Vitamin D modifies cytokine secretion pattern, decreases proinflammatory cytokines production, and increases anti-inflammatory signaling [6]. It was also found that vitamin D can alter the expression of MHC class II and costimulatory molecules on the surface of antigen presenting cells (APCs) [6]. T cells can be affected by vitamin D both directly and through its effect on APCs [7]. Vitamin D also inhibits differentiation of B cells, their proliferation, and immunoglobulin production [6]. 1,25(OH)2D downregulates Th17 lymphocytes activity and IL-17 production that are factors potentially affecting the development of autoimmunity [8]. Moreover, an association between vitamin D and regulatory T cells was proven [9]. All these actions of vitamin D seem to induce more tolerogenic phenotype of immune cells.

The biological action of 1,25(OH)2D is mediated by the nuclear receptor (vitamin D receptor (VDR)), which is a ligand-activated transcription factor. It regulates several target genes, including those affecting the immune system [10]. Moreover, VDR is expressed among others in monocytes, dendritic cells, and activated B or T lymphocytes [11]. It is postulated that different polymorphic variants of VDR, through their effect on receptor structure or expression, could influence the risk of autoimmune disorders [4]. To date, four single nucleotide polymorphisms (SNPs) of VDR gene have been most extensively studied in AIT patients: rs2228570 (FokI) in exon 2, rs1544410 (BsmI) and rs7975232 (ApaI) in intron 8, and rs731236 (TaqI) in exon 9. Their role has been confirmed in some autoimmune diseases, for example, type 1 diabetes (T1D), rheumatoid arthritis (RA), or multiple sclerosis (MS) [12]. Rs2228570 is known to be a functional SNP. In the case of rs1544410, rs7975232, and rs731236 functional status is questionable, although these three SNPs reliably represent the largest haplotype block of VDR gene [13]. The results obtained so far in AIT are heterogeneous and there is no clear conclusion about the role of these VDR variants in the development of thyroid autoimmunity [1416]. Another candidate SNP of VDR that has not been yet studied in AIT is rs11568820 (Cdx2). This functional SNP, located in the promoter region of VDR, has been found to affect the immune system [17]. Therefore, the study aimed to assess the relationship between these five VDR SNPs (rs2228570, rs1544410, rs7975232, rs731236, and rs11568820) and AIT susceptibility in the Caucasian Polish population. The association between these SNPs and selected clinical features of the disease was additionally analyzed.

2. Subjects and Methods

2.1. Subjects

223 unrelated adult AIT patients of Caucasian Polish origin were enrolled in the study. They were recruited from endocrinology outpatient clinics over the period of two years, from January 2016 to September 2017. The diagnosis of AIT was made on the basis of standard criteria that include at least three from (1) clinical symptoms of hypothyroidism, (2) hypo- or euthyroidism found by thyroid function tests (TSH), (3) elevated anti-thyroid antibodies level (thyroid peroxidase antibodies (TPOAb) and/or thyroglobulin antibodies (TgAb)), and (4) thyroid ultrasound findings typical for AIT. In a subgroup of patients the diagnosis was confirmed by fine needle aspiration biopsy (FNAB). All patients were diagnosed with hypothyroidism (overt or subclinical dysfunction). However, on enrollment, they were euthyroid (TSH values were within the normal range) on levothyroxine supplementation for at least three months.

Control group consisted of 130 unrelated healthy volunteers of Caucasian Polish origin, predominantly blood donors. They were recruited within the same geographical location as patients. In the control group thyroid gland disorders were excluded by anamnesis, thyroid ultrasonography, and laboratory tests (TSH, TPOAb, and TgAb). None of them had any known family history of AITDs. Basic characteristics of patients and controls are presented in Table 1. History of any other autoimmune, neoplastic, or chronic inflammatory disorder led to exclusion from the study in both patient and control groups. All participants gave written informed consents. The study protocol was approved by the ethics committee of the Poznan University of Medical Sciences (approval number: 443/15).

AITHealthy controls

Age (SD) [years]46.88 (13.50)35.42 (14.96)
Female to male ratio21.3:13.48:1
Thyroid volume (SD) [ml]11.39 (8.80)10.22 (4.09)
TSH (SD) [mIU/l]2.01 (2.08)1.52 (0.49)
Levothyroxine dose (SD) [μg/kg body mass]1.19 (0.50)-
Age at diagnosis (SD) [years]43.37 (14.57)-
Follow-up time (IQR) [years]3 (1-6)-
Elevated TPOAb [%]92.83none
Elevated TgAb [%]64.13none

AIT: autoimmune thyroiditis.
Data are mean (SD, standard deviation), median (IQR, interquartile range), or percentage.

Detailed thyroid ultrasound (US) reports were available in a subgroup of 98 AIT patients. These patients were further subdivided into AIT with nodules (at least one solid lesion ≥ 5 mm in diameter) and AIT without nodules (AIT patients without any solid lesions ≥ 5 mm in diameter). None of the observed thyroid nodules were autonomous. Patients with nodules suspected for malignancy were referred for FNAB, and those with lesions of Bethesda category II or higher were excluded from further analyses. None of the female patients in this group were pregnant or have delivered in the last 12 months.

2.2. Methods

In both patients and controls genomic DNA was isolated from whole blood using spin column-based DNA extraction kits, according to the manufacturer’s instructions (NucleoSpin® Blood L kit, Macherey-Nagel).

Rs2228570, rs1544410, rs7975232, and rs731236 VDR SNPs were genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). The target DNA fragments were amplified by thermal cycling (MJ Mini Thermal Cycler, Bio-Rad Laboratories) with the use of Color OptiTaq PCR Master Mix (EURx, Poland). PCR products were then digested with appropriate restriction endonucleases. Primer sequences, conditions of the PCR reactions, and the size of products yielded are shown in Table 2 [1820]. The products of digestion were separated in a 2% agarose gel stained with SYBR® Safe (Thermo Fisher). SNP variants were indicated by (1) nucleotide present in the polymorphic site or (2) the first letter of the restriction enzyme, with uppercase letters in case of the absence and lowercase letters in case of the presence of the cut site.

SNPForward primerReverse primerPCR conditionsCyclesPCR product size (bp)EnzymeDigested products size (bp)

rs22285705’-AGCTGGCCCTGGCACTGACTCTGGCT-3’5’-ATGGAAACACCTTGCTTCTTCTCCCTC-3’95°C for 5 min; (95°C for 30 s; 63°C for 30 s; 72°C for 30 s); 72°C for 7 min33267FokICC: 267 
CT: 267, 197, and 70 
TT: 197 and 70

rs15444105’-ACTTGCATGAGGAGGAGCATGTC-3’5’-GGAGAGGAGCCTGTGTCCCATTTG-3’95°C for 5 min; (95°C for 45 s; 61°C for 30 s; 72°C for 45 s); 72°C for 7 min36812BsmIGG: 475 and 337 
GA: 812, 475, and 337 
AA: 812

rs7975232 / rs7312365’-GGTGGGATTGAGCAGTG-3’5’-ATGCTGCACTCAGGCTG-3’95°C for 5 min; (95°C for 30 s; 56°C for 30 s; 72°C for 45 s); 72°C for 7 min35280ApaI / TaqIApaI:  
CC: 256 and 24 
CA: 280, 256, and 24 
AA: 280 
TT: 280 
TC: 180 and 100 
CC: 280, 180, and 100

SNP, single nucleotide polymorphism; bp, base pair.

Rs11568820 genotyping was performed by TaqMan real-time PCR with the use of Roche LightCycler 2.0 thermal cycler. The reaction comprised of Probe qPCR Master Mix E0420 (EURx) and a TaqMan® SNP Genotyping Assay ID: C___2880808_10 (Thermo Fisher). The reaction program was as follows: 37°C, 2 min, Uracil-N-glycosylase pretreatment; 95°C, 10 min, predenaturation; 40x: 92°C, 15 s, ramp rate (RR) 5°C/s; 60°C, 1 min, RR 5°C/s with signal acquisition. Figure 1 presents a result sample of rs11568820 genotyping.

To confirm the accuracy of the aforementioned methods randomly selected samples were verified by direct sequencing.

In a subgroup of patients thyroid ultrasonography was performed on enrollment, using a ProSound SSD-3500SX ultrasound system with a 7.5 MHz linear transducer (Hitachi Aloka Medical, Japan). The volume of the thyroid gland [ml] was calculated as the sum of the two lobes (each lobe volume estimated with the formula π/6 x height [cm] x width [cm] x depth [cm]).

2.3. Statistics

The distribution of genotypes in patients and controls was tested for deviations from Hardy-Weinberg equilibrium (HWE) by chi-squared test (χ2) [21]. The comparison of allele frequencies and genotype distributions between groups was performed by χ2, Fisher’s exact, Fisher-Freeman-Halton, or χ2 for trend tests, where appropriate. Odds ratio (OR) with 95% confidence interval (95% CI) was calculated as a measure of the strength of association. To compare the clinical parameters between groups and/or subgroups (e.g., age, duration of the disease, and thyroid volume) we used Student’s t-test or ANOVA for normally distributed variables and Mann–Whitney U test or Kruskal-Wallis test for variables that were not normally distributed. In the case of significant differences between-groups found by multiple comparison tests (ANOVA or Kruskal-Wallis), post-hoc tests were used for pairwise comparisons. Analyses were performed using GraphPad Prism version 7.03 for Windows (GraphPad Software). A p-value of < 0.05 was considered statistically significant.

Linkage disequilibrium (LD) between studied SNPs was calculated to construct haplotype blocks. In the next step, haplotype frequencies were compared between the groups (with the use of Haploview 4.2 version, Broad Institute) [22].

The numerical data were expressed as mean ± standard deviation or median and interquartile range (for data not normally distributed).

3. Results

3.1. Allele and Genotype Frequencies

Genotype distributions were consistent with HWE for the studied VDR polymorphisms in both groups (p > 0.05). When comparing genotype distributions and allele frequencies of five studied VDR SNPs between the AIT group and controls, no statistically significant differences were observed. Allele and genotype frequencies of the studied SNPs and results of the case-control analysis are presented in Table 3.

n=223 (%)
Healthy controls
n=130 (%)
pOR (95% CI)

rs2228570 (FokI)


CC (FF)72 (32.29)38 (29.23)0.71 (0.84)
CT (Ff)98 (43.95)63 (48.46)
TT (ff)53 (23.77)29 (22.31)


C (F)242 (54.26)139 (53.46)0.841.03 (0.76-1.40)
T (f)204 (45.74)121 (46.54)

rs1544410 (BsmI)


GG (bb)86 (39.09)43 (35.83)0.51 (0.31)
GA (bB)99 (45.00)52 (43.33)
AA (BB)35 (15.91)25 (20.83)


G (b)271 (61.59)138 (57.50)0.301.19 (0.86-1.63)
A (B)169 (38.41)102 (42.50)

rs7975232 (ApaI)


CC (aa)53 (24.65)28 (22.40)0.78 (0.94)
CA (aA)112 (52.09)70 (56.00)
AA (AA)50 (23.26)27 (21.60)


C (a)218 (50.70)126 (50.40)0.941.01 (0.74-1.39)
A (A)212 (49.30)124 (49.60)

rs731236 (TaqI)


TT (TT)79 (36.74)42 (33.60)0.63 (0.38)
TC (Tt)106 (49.30)61 (48.80)
CC (tt)30 (13.95)22 (17.60)


T (T)264 (61.40)145 (58.00)0.381.15 (0.84-1.58)
C (t)166 (38.60)105 (42.00)

rs11568820 (Cdx2)


GG118 (71.08)91 (76.47)0.54 (0.39)
GA46 (27.71)26 (21.85)
AA2 (1.20)2 (1.68)


G282 (84.94)208 (87.39)0.411.23 (0.76-2.01)
A50 (15.06)30 (12.61)

AIT, autoimmune thyroiditis; OR, odds ratio; CI, confidence interval; χ2 test for trend.
3.2. Haplotype Analysis

According to D’ values the strongest LD was observed between three SNPs located in the 3’ end of the VDR gene: rs731236, rs7975232, and rs1544410. Instead, the weakest LD was found between rs2228570 SNP and studied 3’ end polymorphisms (detailed results in Table 4). Haplotype block containing three linked SNPs (rs731236, rs7975232, and rs1544410) was constructed, according to Gabriel et al. [20]. Haplotype analysis was possible in 86.10% of patients and 93.85% of controls. In both the patient and the control groups the most frequent haplotype was TCG (Tab) (49.6% and 50.3%, resp.). The second most frequent haplotype was CAA (tAB), again with similar percentages in both groups (37.1% in AIT and 40.1% in controls). Other haplotypes were rare. Haplotype frequencies did not differ significantly between patients and controls (Table 5).

Locus 1Locus 2D’ (CI)LODr2

rs731236 (TaqI)rs7975232 (ApaI)1.0 (0.94-1.0)32.850.73
rs731236 (TaqI)rs1544410 (BsmI)0.95 (0.89-0.98)42.130.90
rs731236 (TaqI)rs2228570 (FokI)0.05 (0.0-0.26)0.050.00
rs731236 (TaqI)rs11568820 (Cdx2)0.50 (0.18-0.73)1.510.05
rs7975232 (ApaI)rs1544410 (BsmI)0.98 (0.9-1.0)30.170.70
rs7975232 (ApaI)rs2228570 (FokI)0.06 (0.0-0.24)0.080.00
rs7975232 (ApaI)rs11568820 (Cdx2)0.59 (0.21-0.8)1.590.05
rs1544410 (BsmI)rs2228570 (FokI)0.06 (0.0-0.26)0.060.00
rs1544410 (BsmI)rs11568820 (Cdx2)0.44 (0.14-0.68)1.290.04
rs2228570 (FokI)rs11568820 (Cdx2)0.12 (0.01-0.52)0.050.00

D’, value of D prime between the loci; LOD, likelihood odds ratio;
r2, correlation coefficient between two loci; CI, confidence interval.

HaplotypeAIT (%)Healthy controls (%)pOR95% CI

TCG (Tab)49.650.30.870.970.71-1.34
CAG (tAb)

AIT, autoimmune thyroiditis.
The table shows haplotypes with frequencies over 1%.
3.3. Association of VDR SNPs with Selected Clinical Features of AIT

The association between VDR SNPs and the ultrasonographic features was analyzed in a subgroup of 98 AIT patients. Patients were divided into two subgroups: AIT with nodules and AIT without nodules (as described in Methods Section). These subgroups did not differ significantly in age or gender distribution and time since diagnosis of the disease (p > 0.05). Thyroid volume (TV) was slightly higher in AIT with nodules, as expected (12.65 ± 9.29 vs. 9.38 ± 7.05, p = 0.06). AIT patients with nodules also had lower daily doses of levothyroxine comparing to AIT without nodules group (75 μg (50-100) vs. 100 μg (75-125), p < 0.01). However, we did not observe any significant differences in genotype distributions or allele frequencies of the studied VDR SNPs between these two AIT subgroups (data shown in Supplementary Table 1S).

The relationship between VDR variants and TV in AIT patients was analyzed. It revealed that the TV differed significantly between different genotypes of rs1544410 and rs7975232 (p = 0.03 and 0.04, resp., Figure 2), but not between rs2228570, rs731236, and rs11568820 genotypes (p = 0.29, 0.12, and 0.54, resp.). In the case of rs1544410 and rs7975232, Dunn’s post-hoc test was performed for pairwise comparison. It showed that differences between GG and AA and between GG and GA genotypes of rs1544410 were of borderline significance (p = 0.07 in both cases). In the case of rs7975232 SNP the difference between CC and AA genotypes showed a trend toward significance (p = 0.06). The subgroups created according to rs1544410 and rs7975232 genotypes did not differ significantly in age, gender, time since diagnosis of the disease, or daily dose of levothyroxine (p > 0.05 in all cases; for details, see Supplementary Tables 2S and 3S).

We also tested the association between studied VDR SNPs and age at AIT diagnosis, although no significant differences were found (data not shown).

4. Discussion

VDR gene polymorphisms have been extensively studied for their potential association with susceptibility to different autoimmune and inflammatory disorders, for example, MS, T1D, RA, systemic lupus erythematosus (SLE), inflammatory bowel diseases, or tuberculosis [12]. Meta-analyses of different studies confirmed some of these associations, such as rs2228570 and rs731236 in RA [23], rs1544410 and rs731236 in T1D [24], rs7975232 in Crohn’s disease, or rs731236 in ulcerative colitis [25]. The interstudy heterogeneity is often reported what may be in part explained by ethnic differences observed in allele and genotype distributions of VDR SNPs [26].

In our case-control study of Polish AIT patients, we failed to find any significant associations between rs2228570, rs1544410, rs7975232, and rs731236 VDR SNPs and the risk of disease development. These results are in accordance with recent European population study from Italy, where no significant difference in VDR SNPs distribution between AIT and controls was also found (rs1544410, rs7975232, and rs731236 were assessed) [27]. On the contrary, there are publications where VDR SNPs were shown to be associated with AIT in Caucasians; rs2228570 was found to influence the disease risk in Serbian patients and rs1544410 in Croatians [28, 29]. Moreover, two studies from Turkey confirmed the association with rs731236 SNP [30, 31]. In Asians rs2228570 SNP was reported as a risk factor of AIT development in several studies [3234], although Meng et al. failed to find such an association among the Chinese [35]. In a recent study of Iranian patients, Zarrin et al. showed that rs2228570 VDR SNP might be weakly associated with the disease risk [36].

In 2013 Feng et al. performed a meta-analysis of available studies on the role of VDR SNPs in AITD pathogenesis. They found a relationship between rs1544410 and rs731236 SNP and AITD risk. In Europeans, similarly to overall results, B allele of rs1544410 (BsmI) and t allele of rs731236 (TaqI) were associated with reduced AITD risk. Unfortunately, the authors did not perform stratification by AITD subtypes [14]. In 2017 Wang et al. assessed studies on VDR SNPs in Hashimoto’s thyroiditis and found rs2228570 SNP as a disease risk marker. However, the analysis of subgroups by ethnicity showed no significant differences between patients and controls among Caucasians [15]. In 2018 Gao et al. analyzed the influence of VDR SNPs on AITD risk again, adding newly published studies. They found an association between rs2228570 or rs731236 SNPs and AITD. This observation was also confirmed for AIT patients in stratification by AITD subtypes. In Caucasian AITD patients the only significant association was found for rs1544410 SNP [16]. It cannot be excluded that the role of VDR gene variants is not the same in two phenotypes of AITD-AIT and Graves’ disease. The impact of VDR SNPs on predisposition to autoimmunity may also vary depending on studied population or ethnicity.

Rs11568820 (Cdx2) VDR SNP was not assessed in AITD patients to date (according to the PubMed database, 31.04.2019). This SNP was previously shown to influence the immune system and modify the risk of some infectious diseases (e.g., tuberculosis and rubella) [17, 37, 38]. However, we did not observe a significant association between rs11568820 SNP and AIT risk. The role of rs11568820 SNP was also studied in different autoimmune diseases, among others in T1D, MS, vitiligo, and psoriasis. Similar to our results, no clear connection with autoimmunity was confirmed [3942].

In our group of unrelated Polish Caucasians, strong LD was present between three VDR SNPs: rs1544410, rs7975232, and rs731236. Rs2228570 and rs11568820 SNPs were not in LD with any of the other studied SNPs. Most frequent haplotypes were TCG (Tab) and CAA (tAB) in both patients and controls, what is in agreement with earlier studies of Caucasian population [27, 43]. We did not confirm the association between VDR haplotypes and AIT. Giovinazzo et al. and Meng et al. also did not observe such associations, contrary to the study by Stefanic et al. where the Tab variant was found to confer susceptibility to AIT [27, 29, 35].

Our study showed that there might be a weak association between rs1544410 and rs7975232 SNPs and volume of thyroid gland in AIT patients. The influence of VDR SNPs on TV has not been already reported. Meng et al. and Inoue et al. studied the relationship between VDR SNPs and goiter degree (evaluated clinically), but they found no associations [32, 35]. Vitamin D level was previously reported to be positively correlated with TV in AIT [44]. One of the possible explanations of the observed associations is that some VDR genotypes contribute to the atrophic variant of AIT, in contrast to goitrous form [45]. It was also proven that vitamin D may influence apoptosis of thyrocytes that could lead to volume changes [46]. As TV may vary according to different genetic and exogenous factors [47], aforementioned results must be interpreted cautiously and need further confirmation.

5. Conclusions

We found no differences in genotype or allele distributions between AIT patients and the control group for any of the studied VDR SNPs. It suggests that VDR gene is probably not a strong susceptibility factor for AIT development, at least in the Polish population. It cannot be excluded that association may be minor in case of Caucasian patients comparing to Asians and larger sample size is needed to confirm VDR role in determining disease risk.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Adam Maciejewski and Michał J. Kowalczyk contributed equally to this work.


This study was supported by Grant from Poznan University of Medical Sciences No. 502-14-02221355-41204.

Supplementary Materials

Table 1S: Allele and genotype frequencies of VDR SNPs in “AIT with nodules” vs. “AIT without nodules.” Table 2S: Characteristic of AIT subgroups according to genotypes of rs1544410 (BsmI) VDR polymorphism. Table 3S: Characteristic of AIT subgroups according to genotypes of rs7975232 (ApaI) VDR polymorphism. (Supplementary Materials)


  1. M. P. J. Vanderpump, “The epidemiology of thyroid disease,” British Medical Bulletin, vol. 99, no. 1, pp. 39–51, 2011. View at: Publisher Site | Google Scholar
  2. Y. Tomer, “Mechanisms of autoimmune thyroid diseases: from genetics to epigenetics,” Annual Review of Pathology: Mechanisms of Disease, vol. 9, no. 1, pp. 147–156, 2014. View at: Publisher Site | Google Scholar
  3. W. M. Wiersinga, “Clinical relevance of environmental factors in the pathogenesis of autoimmune thyroid disease,” Endocrinology and Metabolism Clinics of North America, vol. 31, no. 2, pp. 213–222, 2016. View at: Publisher Site | Google Scholar
  4. B. Altieri, G. Muscogiuri, L. Barrea et al., “Does vitamin D play a role in autoimmune endocrine disorders? A proof of concept,” Reviews in Endocrine and Metabolic Disorders, vol. 18, no. 3, pp. 335–346, 2017. View at: Publisher Site | Google Scholar
  5. R. A. G. Khammissa, J. Fourie, M. H. Motswaledi, R. Ballyram, J. Lemmer, and L. Feller, “The biological activities of vitamin D and Its receptor in relation to calcium and bone homeostasis, cancer, immune and cardiovascular systems, skin biology, and oral health,” BioMed Research International, vol. 2018, Article ID 9276380, 9 pages, 2018. View at: Publisher Site | Google Scholar
  6. B. Prietl, G. Treiber, T. R. Pieber, and K. Amrein, “Vitamin D and immune function,” Nutrients, vol. 5, no. 7, pp. 2502–2521, 2013. View at: Publisher Site | Google Scholar
  7. M. Hewison, “An update on vitamin D and human immunity,” Clinical Endocrinology, vol. 76, no. 3, pp. 315–325, 2012. View at: Publisher Site | Google Scholar
  8. S. Joshi, L.-C. Pantalena, X. K. Liu et al., “1,25-Dihydroxyvitamin D3 ameliorates Th17 autoimmunity via transcriptional modulation of interleukin-17A,” Molecular and Cellular Biology, vol. 31, no. 17, pp. 3653–3669, 2011. View at: Publisher Site | Google Scholar
  9. E. S. Chambers and C. M. Hawrylowicz, “The impact of vitamin D on regulatory T cells,” Current Allergy and Asthma Reports, vol. 11, no. 1, pp. 29–36, 2011. View at: Publisher Site | Google Scholar
  10. G. Heine, U. Niesner, H.-D. Chang et al., “1,25-dihydroxyvitamin D3 promotes IL-10 production in human B cells,” European Journal of Immunology, vol. 38, no. 8, pp. 2210–2218, 2008. View at: Publisher Site | Google Scholar
  11. E. Peelen, S. Knippenberg, A.-H. Muris et al., “Effects of vitamin D on the peripheral adaptive immune system: a review,” Autoimmunity Reviews, vol. 10, no. 12, pp. 733–743, 2011. View at: Publisher Site | Google Scholar
  12. W. Dankers, E. M. Colin, J. P. van Hamburg et al., “Vitamin D in autoimmunity: molecular mechanisms and therapeutic potential,” Frontiers in Immunology, vol. 7, p. 697, 2016. View at: Google Scholar
  13. Y. Fang, J. B. J. van Meurs, A. D'Alesio et al., “Promoter and 3-untranslated-region haplotypes in the vitamin D receptor gene predispose to osteoporotic fracture: the Rotterdam Study,” American Journal of Human Genetics, vol. 77, no. 5, pp. 807–823, 2005. View at: Publisher Site | Google Scholar
  14. M. Feng, H. Li, S. F. Chen, W. F. Li, and F. B. Zhang, “Polymorphisms in the vitamin D receptor gene and risk of autoimmune thyroid diseases: a meta-analysis,” Endocrine Journal, vol. 43, no. 2, pp. 318–326, 2013. View at: Publisher Site | Google Scholar
  15. J. Wang, S. Lv, G. Chen et al., “Meta-analysis of the association between vitamin D and autoimmune thyroid disease,” Nutrients, vol. 7, no. 4, pp. 2485–2498, 2015. View at: Publisher Site | Google Scholar
  16. X. R. Gao and Y. G. Yu, “Meta-analysis of the association between vitamin d receptor polymorphisms and the risk of autoimmune thyroid disease,” International Journal of Endocrinology, Article ID 2846943, 2018. View at: Google Scholar
  17. V. Meyer and L. Bornman, “Cdx-2 polymorphism in the vitamin D receptor gene (VDR) marks VDR expression in monocyte/macrophages through VDR promoter methylation,” Immunogenetics, vol. 70, no. 8, pp. 523–532, 2018. View at: Publisher Site | Google Scholar
  18. S. S. Harris, T. R. Eccleshall, C. Gross, B. Dawson-Hughes, and D. Feldman, “The vitamin D receptor start codon polymorphism (FokI) and bone mineral density in premenopausal American black and white women,” Journal of Bone and Mineral Research, vol. 12, no. 7, pp. 1043–1048, 1997. View at: Publisher Site | Google Scholar
  19. Y. Ozaki, S. Nomura, M. Nagahama, C. Yoshimura, H. Kagawa, and S. Fukuhara, “Vitamin-D receptor genotype and renal disorder in japanese patients with systemic lupus erythematosus,” Nephron, vol. 85, no. 1, pp. 86–91, 2000. View at: Publisher Site | Google Scholar
  20. J.-Y. Oh and E. Barrett-Connor, “Association between vitamin D receptor polymorphism and type 2 diabetes or metabolic syndrome in community-dwelling older adults: the Rancho Bernardo study,” Metabolism, vol. 51, no. 3, pp. 356–359, 2002. View at: Publisher Site | Google Scholar
  21. S. Rodriguez, T. R. Gaunt, and I. N. M. Day, “Hardy-Weinberg equilibrium testing of biological ascertainment for Mendelian randomization studies,” American Journal of Epidemiology, vol. 169, no. 4, pp. 505–514, 2009. View at: Publisher Site | Google Scholar
  22. J. C. Barrett, B. Fry, J. Maller, and M. J. Daly, “Haploview: analysis and visualization of LD and haplotype maps,” Bioinformatics, vol. 21, no. 2, pp. 263–265, 2005. View at: Publisher Site | Google Scholar
  23. K. Tizaoui and K. Hamzaoui, “Association between VDR polymorphisms and rheumatoid arthritis disease: Systematic review and updated meta-analysis of case-control studies,” Immunobiology, vol. 220, no. 6, pp. 807–816, 2015. View at: Publisher Site | Google Scholar
  24. O. A. Sahin, D. Goksen, A. Ozpinar, M. Serdar, and H. Onay, “Association of vitamin D receptor polymorphisms and type 1 diabetes susceptibility in children: a meta-analysis,” Endocrine Connections, vol. 6, no. 3, pp. 159–171, 2017. View at: Publisher Site | Google Scholar
  25. L. Wang, Z. T. Wang, J. J. Hu, R. Fan, J. Zhou, and J. Zhong, “Polymorphisms of the vitamin D receptor gene and the risk of inflammatory bowel disease: a meta-analysis,” Genetics and Molecular Research, vol. 13, no. 2, pp. 2598–2610, 2014. View at: Publisher Site | Google Scholar
  26. K. Köstner, N. Denzer, C. S. L. Müller, R. Klein, W. Tilgen, and J. Reichrath, “The relevance of Vitamin D Receptor (VDR) gene polymorphisms for cancer: a review of the literature,” Anticancer Reseach, vol. 29, no. 9, pp. 3511–3536, 2009. View at: Google Scholar
  27. S. Giovinazzo, T. M. Vicchio, R. Certo et al., “Vitamin D receptor gene polymorphisms/haplotypes and serum 25(OH)D3 levels in Hashimoto’s thyroiditis,” Endocrine Journal, vol. 55, no. 2, pp. 599–606, 2017. View at: Publisher Site | Google Scholar
  28. J. Djurovic, O. Stojkovic, O. Ozdemir et al., “Association between FokI , ApaI and TaqI RFLP polymorphisms in VDR gene and Hashimoto's thyroiditis: preliminary data from female patients in Serbia,” International Journal of Immunogenetics, vol. 42, no. 3, pp. 190–194, 2015. View at: Publisher Site | Google Scholar
  29. M. Stefanic, S. Papic, M. Suver, L. Glavas-Obrovac, and I. Karner, “Association of vitamin D receptor gene 3-variants with Hashimoto's thyroiditis in the Croatian population,” International Journal of Immunogenetics, vol. 35, no. 2, pp. 125–131, 2008. View at: Publisher Site | Google Scholar
  30. D. Yazici, D. Yavuz, O. Tarcin, S. Sancak, O. Deyneli, and S. Akalin, “Vitamin D receptor gene ApaI, TaqI, FokI and BsmI polymorphisms in a group of Turkish patients with Hashimoto's thyroiditis,” Minerva Endocrinologica, vol. 38, no. 2, pp. 195–201, 2013. View at: Google Scholar
  31. B. Guleryuz, F. Akin, M. T. Ata, M. M. Dalyanoğlu, and S. Turgut, “Vitamin-D receptor (VDR) gene polymorphisms (TaqI, FokI) in Turkish patients with Hashimoto’s thyroiditis: Relationship to the levels of Vit-D and Cytokines,” Endocrine, Metabolic & Immune Disorders—Drug Targets, vol. 16, no. 2, pp. 131–139, 2016. View at: Publisher Site | Google Scholar
  32. N. Inoue, M. Watanabe, N. Ishido et al., “The functional polymorphisms of VDR, GC and CYP2R1 are involved in the pathogenesis of autoimmune thyroid diseases,” Clinical & Experimental Immunology, vol. 178, no. 2, pp. 262–269, 2014. View at: Publisher Site | Google Scholar
  33. W. Y. Lin, L. Wan, C. H. Tsai, R. H. Chen, C. C. Lee, and F. J. Tsai, “Vitamin D receptor gene polymorphisms are associated with risk of Hashimoto's thyroiditis in Chinese patients in Taiwan,” Journal of Clinical Laboratory Analysis, vol. 20, no. 3, pp. 109–112, 2006. View at: Publisher Site | Google Scholar
  34. Y. Ban, M. Taniyama, and Y. Ban, “Vitamin D receptor gene polymorphisms in hashimoto's thyroiditis,” Thyroid, vol. 11, no. 6, pp. 607-608, 2001. View at: Publisher Site | Google Scholar
  35. S. Meng, S. T. He, W. J. Jiang et al., “Genetic susceptibility to autoimmune thyroid diseases in a Chinese Han population: Role of vitamin D receptor gene polymorphisms,” Ann Endocrinol (Paris), vol. 76, no. 6, pp. 684–689, 2015. View at: Google Scholar
  36. R. Zarrin, M. Bagheri, A. Mehdizadeh, P. Ayremlou, and A. H. Faghfouri, “The association of FokI and ApaI polymorphisms in vitamin D receptor gene with autoimmune thyroid diseases in the northwest of Iran,” Medical Journal of The Islamic Republic of Iran, vol. 32, no. 1, pp. 18–21, 2018. View at: Publisher Site | Google Scholar
  37. M. Harishankar and P. Selvaraj, “Influence of Cdx2 and TaqI gene variants on vitamin D 3 modulated intracellular chemokine positive T-cell subsets in pulmonary tuberculosis,” Clinical Therapeutics, vol. 39, no. 5, pp. 946–957, 2017. View at: Publisher Site | Google Scholar
  38. I. G. Ovsyannikova, N. Dhiman, I. H. Haralambieva et al., “Rubella vaccine-induced cellular immunity: evidence of associations with polymorphisms in the Toll-like, vitamin A and D receptors, and innate immune response genes,” Human Genetics, vol. 127, no. 2, pp. 207–221, 2010. View at: Publisher Site | Google Scholar
  39. J. Dickinson, D. Perera, A. van der Mei et al., “Past environmental sun exposure and risk of multiple sclerosis: a role for the Cdx-2 Vitamin D receptor variant in this interaction,” Multiple Sclerosis Journal, vol. 15, no. 5, pp. 563–570, 2009. View at: Publisher Site | Google Scholar
  40. B. Frederiksen, E. Liu, J. Romanos et al., “Investigation of the vitamin D receptor gene (VDR) and its interaction with protein tyrosine phosphatase, non-receptor type 2 gene (PTPN2) on risk of islet autoimmunity and type 1 diabetes: the Diabetes autoimmunity study in the young (DAISY),” The Journal of Steroid Biochemistry and Molecular Biology, vol. 133, no. 1, pp. 51–57, 2013. View at: Publisher Site | Google Scholar
  41. I. E. Aydingöz, I. Bingül, S. Dogru-Abbasoglu, P. Vural, and M. Uysal, “Analysis of vitamin D receptor gene polymorphisms in vitiligo,” Dermatology, vol. 224, no. 4, pp. 361–368, 2012. View at: Publisher Site | Google Scholar
  42. X. Zhou, L. Xu, and Y. Li, “The association of polymorphisms of the vitamin D receptor gene with psoriasis in the Han population of northeastern China,” Journal of Dermatological Science, vol. 73, no. 1, pp. 63–66, 2014. View at: Publisher Site | Google Scholar
  43. E. Ramos-Lopez, A. Kurylowicz, T. Bednarczuk, J. Paunkovic, C. Seidl, and K. Badenhoop, “Vitamin D receptor polymorphisms are associated with graves' disease in german and polish but not in serbian patients,” Thyroid, vol. 15, no. 10, pp. 1125–1130, 2005. View at: Publisher Site | Google Scholar
  44. N. C. Bozkurt, B. Karbek, B. Ucan et al., “The association between severity of vitamin D deficiency and Hashimoto's thyroiditis,” Endocrine Practice, vol. 19, no. 3, pp. 479–484, 2013. View at: Publisher Site | Google Scholar
  45. D. Engelbrecht Zantut-Wittmann, L. Persoli, M. A. Tambascia et al., “HLA-DRB1*04 and HLA-DQB1*03 association with the atrophic but not with the goitrous form of chronic autoimmune thyroiditis in a Brazilian Population,” Hormone and Metabolic Research, vol. 36, no. 7, pp. 492–500, 2004. View at: Publisher Site | Google Scholar
  46. E. Kaczmarek, K. Lacka, D. Jarmolowska-Jurczyszyn, A. Sidor, and P. Majewski, “Changes of B and T lymphocytes and selected apopotosis markers in Hashimoto's thyroiditis,” Journal of Clinical Pathology, vol. 64, no. 7, pp. 626–630, 2011. View at: Publisher Site | Google Scholar
  47. P. S. Hansen, T. H. Brix, F. N. Bennedbæk, S. J. Bonnema, K. O. Kyvik, and L. Hegedüs, “Genetic and environmental causes of individual differences in thyroid size: a study of healthy danish twins,” The Journal of Clinical Endocrinology & Metabolism, vol. 89, no. 5, pp. 2071–2077, 2004. View at: Publisher Site | Google Scholar

Copyright © 2019 Adam Maciejewski et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

1171 Views | 364 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.