Vitamin D Metabolism Gene Polymorphisms: Comparison
Please note this is a comparison between Version 1 by Ulises de la Cruz-Mosso and Version 3 by Lily Guo.

Vitamin D is a potent immunonutrient that through its main metabolite calcitriol, regulates the immunomodulation of macrophages, dendritic cells, T and B lymphocytes, which express the vitamin D receptor (VDR), and they produce and respond to calcitriol. Genetic association studies have shown that up to 65% of vitamin D serum variance may be explained due to genetic background. 

  • vitamin D
  • polymorphism
  • calcidiol
  • calcitriol
  • VDBP
  • GC
  • CYP2R1
  • CYP27B1
  • VDR

1. Introduction

The etiology and progression of autoimmune diseases (AIDs) are multifactorial and complex[1]. Genetic and environmental factors such as nutrients have been proposed to partially explain the pathophysiology progression of autoimmunity[2]. Notably, vitamin D regulates the growth and differentiation of various cells of the immune system such as macrophages, dendritic cells, T cells, and B cells, which are able to express the vitamin D receptor (VDR), produce and respond to the active form of vitamin D, calcitriol (1α,25(OH)

The etiology and progression of autoimmune diseases (AIDs) are multifactorial and complex [1]. Genetic and environmental factors such as nutrients have been proposed to partially explain the pathophysiology progression of autoimmunity [2]. Notably, vitamin D regulates the growth and differentiation of various cells of the immune system such as macrophages, dendritic cells, T cells, and B cells, which are able to express the vitamin D receptor (VDR), produce and respond to the active form of vitamin D, calcitriol (1α,25(OH)

2

D

3)[3].

) [3].

In autoimmune diseases, epidemiological studies have reported a high prevalence of vitamin D deficiency by the quantification of calcidiol; this deficiency has been associated with worse disease clinical activity and progression of systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and multiple sclerosis (MS) diseases[4]. Likewise, several clinical trials and murine studies using cholecalciferol supplementation have described vitamin D′s immunomodulatory properties in autoimmune diseases[5]. Lower calcidiol serum levels are a strong predictor for worse MS clinical manifestations [6] as well as RA patients show lower calcidiol levels than control subjects (CS), and they have shown a negative association of calcidiol serum levels with clinical disease activity[7]. Regarding SLE patients, a higher prevalence of calcidiol serum deficiency than the general population has also been reported[8], and the lower serum calcidiol in SLE patients is associated with high clinical disease activity[9].

In autoimmune diseases, epidemiological studies have reported a high prevalence of vitamin D deficiency by the quantification of calcidiol; this deficiency has been associated with worse disease clinical activity and progression of systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and multiple sclerosis (MS) diseases [4]. Likewise, several clinical trials and murine studies using cholecalciferol supplementation have described vitamin D′s immunomodulatory properties in autoimmune diseases [5]. Lower calcidiol serum levels are a strong predictor for worse MS clinical manifestations [6] as well as RA patients show lower calcidiol levels than control subjects (CS), and they have shown a negative association of calcidiol serum levels with clinical disease activity [7]. Regarding SLE patients, a higher prevalence of calcidiol serum deficiency than the general population has also been reported [8], and the lower serum calcidiol in SLE patients is associated with high clinical disease activity [9].

Calcidiol serum deficiency in patients with autoimmune diseases and the general population could be attributed to several factors, including lack of exposure to sunlight, skin pigmentation, sunscreen use, nutrient intake deficiencies, age, use of glucocorticoids, and the genetic background of the populations[10]. The potential roles of 35 genes that could modulate the vitamin D serum levels status have been reported in previous studies, highlighting that multiple single nucleotide polymorphisms (SNPs) in these genes are associated with lower calcidiol serum levels[11], such as the SNPs described in the vitamin D binding protein (VDBP; rs2282679

Calcidiol serum deficiency in patients with autoimmune diseases and the general population could be attributed to several factors, including lack of exposure to sunlight, skin pigmentation, sunscreen use, nutrient intake deficiencies, age, use of glucocorticoids, and the genetic background of the populations [10]. The potential roles of 35 genes that could modulate the vitamin D serum levels status have been reported in previous studies, highlighting that multiple single nucleotide polymorphisms (SNPs) in these genes are associated with lower calcidiol serum levels [11], such as the SNPs described in the vitamin D binding protein (VDBP; rs2282679

GC

), 25-hydroxylase (rs10751657,

CYP2R1

), 1α-hydroxylase (rs10877012,

CYP27B1

) and the vitamin D receptor (

FokI

(rs2228570),

Bsml

(rs1544410),

Apal

(rs7975232), and

Taql

(rs731236)

VDR

). Therefore, the aim of this comprehensive literature review was to discuss the current findings of the functional SNPs in

GC

,

CYP2R1

,

CYP27B1

, and

VDR

related to genetic risk and the most common clinical features of MS, RA, and SLE.

2. Vitamin D Status and Genetic Evidence in the Populations

Regarding vitamin D deficiency, several factors have been described, mainly the lack of exposure to sunlight, latitude, the season of the year, skin pigmentation, and use of sunscreen; other factors involved in vitamin D deficiency are diet, age, pharmacotherapy administered (antiepileptic and glucocorticoids), and particularly, several studies have described that genetic differences between individuals and populations such as genetic polymorphisms could influence the vitamin D deficiencies presented in all populations around the world[10].

Regarding vitamin D deficiency, several factors have been described, mainly the lack of exposure to sunlight, latitude, the season of the year, skin pigmentation, and use of sunscreen; other factors involved in vitamin D deficiency are diet, age, pharmacotherapy administered (antiepileptic and glucocorticoids), and particularly, several studies have described that genetic differences between individuals and populations such as genetic polymorphisms could influence the vitamin D deficiencies presented in all populations around the world [10].
Multiple SNPs in

GC

,

CYP2R1

,

CYP27B1

, and

VDR genes are associated with lower calcidiol serum levels[11]. Moreover, vitamin D serum deficiencies are present in a high frequency in a healthy population, which could be related to the SNPs’ presence in these genes that may modify the response to supplementation of vitamin D in health and disease. A study conducted in healthy Iranian adolescents described the differential effect of the

genes are associated with lower calcidiol serum levels [11]. Moreover, vitamin D serum deficiencies are present in a high frequency in a healthy population, which could be related to the SNPs’ presence in these genes that may modify the response to supplementation of vitamin D in health and disease. A study conducted in healthy Iranian adolescents described the differential effect of the

CYP2R1

(rs10741657) A > G SNP on the supplementation of 50,000 UI of cholecalciferol weekly over 9 weeks, and showed that participants carrying the AA genotype presented 2.5-fold higher calcidiol serum levels in comparison to those that carrying the GG genotype (OR = 2.5 (1.4–4.4);

p = 0.002)[12]. This is evidence of the role of polymorphisms in genes related to vitamin D metabolism in the variation of the response to vitamin D supplementation, even in healthy conditions.

0.002) [54]. This is evidence of the role of polymorphisms in genes related to vitamin D metabolism in the variation of the response to vitamin D supplementation, even in healthy conditions.

3. Polymorphisms in the Main Key Genes Related to Vitamin D Metabolism

The skin produces around 80% of vitamin D in the form of cholecalciferol when exposed to ultraviolet B (UVB) light at a wavelength of 290–320 nm, and the remaining 20% is obtained from the diet as ergocalciferol from mushrooms or as cholecalciferol from fortified dairy products, fish, and eggs, mainly [13].

The skin produces around 80% of vitamin D in the form of cholecalciferol when exposed to ultraviolet B (UVB) light at a wavelength of 290–320 nm, and the remaining 20% is obtained from the diet as ergocalciferol from mushrooms or as cholecalciferol from fortified dairy products, fish, and eggs, mainly [55].
In blood, cholecalciferol and ergocalciferol bind to the vitamin D binding protein (VDBP), a protein encoded by the

GC

gene to be transported to the liver (

Figure 1a). The

2a). The

GC (rs2282679) SNP was been associated with modulation of vitamin D serum levels in several populations[14][15][16]. In these studies, the main vitamin D circulating metabolite in blood, calcidiol, was used to evaluate vitamin D deficiencies [17].

(rs2282679) SNP was been associated with modulation of vitamin D serum levels in several populations [56,57,58]. In these studies, the main vitamin D circulating metabolite in blood, calcidiol, was used to evaluate vitamin D deficiencies [59].

Figure 12.

Polymorphisms in main key enzymes and proteins associated with vitamin D metabolism: localization and functional effects: (

a

) Vitamin D binding protein (VDBP) (encoded by

GC

gene) binds to ergocalciferol/cholecalciferol in order to be transported to the liver;

GC

(rs2282679) single nucleotide polymorphisms (SNP) due to its location in intron may generate a splicing alteration and exon skipping; (

b

) In the liver, 25 hydroxylase (encoded by

CYP2R1

gene) converts ergocalciferol and cholecalciferol to calcidiol and then calcidiol binds to VDBP to be transported to the kidney;

CYP2R1

(rs10741657) SNP located on the 5′ untranslated region (UTR) region may affect the transcript stabilization and the post-transcriptional control; (

c

) In the kidney, calcidiol is converted to calcitriol by the enzyme 1 alpha hydroxylase (encoded by the

CYP27B1

gene);

CYP27B1

(rs10877012) SNP located on 5′ UTR may affect the transcript stabilization and the post-transcriptional control of mRNA; (

d

) After calcitriol enters target cells and binds to vitamin D receptor (VDR) (encoded by

VDR

gene). Then, the VDR-calcitriol complex in the cytosol is translocated to the nucleus, where it binds to retinoid X receptor (RXR) to form a heterodimer, which interacts with vitamin D response element (VDRE) in vitamin D target genes, i.e., in T helper (Th) lymphocytes to suppress

IL-17A

or activate

FOXP3.

Mainly four SNPs have been described in the

VDR

gene: the

FokI

(rs2228570) located on exon 2, which generates a non-synonymous polymorphism with a change of C > T (also called F > f) and this results in a change of threonine to methionine. The presence of the restriction site

FokI

C allele (F allele), generates a new start codon (ATG) 9 bp after of the common starting site, which translate to an shorter truncated VDR protein of 424 amino acids with more transactivation capacity as a transcription factor than the wild type full-length VDR A isoform (VDRA) of 427 amino acids; the

BsmI

(rs1544410) located on intron 8 presents a change of A > G (also called B > b), could affect messenger ribonucleic acid (mRNA) stability and the gene expression of

VDR

, and also it could generate an alteration in the splice sites for mRNA transcription or a change in the intron regulatory elements of

VDR

;

ApaI

(rs7975232) located on intron 8 of

VDR

presents a change of A > C (also called A > a), does not change the amino acid sequence of the VDR protein, therefore could affect mRNA stability and the gene expression of

VDR

;

TaqI

(rs731236) is located on the exon 9 of

VDR, presents a change of C > T (also called T > t) and generates a synonymous change of the isoleucine amino acid in the coding sequence, therefore it does not change the encoded protein, but it could influence the stability of the mRNA. All these SNPs are related to modulating de vitamin D serum status in health and disease. Ile: isoleucine; Thr: threonine; Met: methionine VDBP: vitamin D binding protein; VDR: vitamin D receptor; RXR: retinoid X receptor; VDRE: vitamin D response elements; UTR: untranslated region; THEM4: thioesterase superfamily member 4; Th: T helper lymphocyte; VDRA: wild type full-length VDR A isoform.

, presents a change of C > T (also called T > t) and generates a synonymous change of the isoleucine amino acid in the coding sequence, therefore it does not change the encoded protein, but it could influence the stability of the mRNA. All these SNPs are related to modulating de vitamin D serum status in health and disease. Ile: isoleucine; Thr: threonine; Met: methionine VDBP: vitamin D binding protein; VDR: vitamin D receptor; RXR: retinoid X receptor; VDRE: vitamin D response elements; UTR: untranslated region; THEM4: thioesterase superfamily member 4; Th: T helper lymphocyte; VDRA: wild type full-length VDR A isoform.
In the liver, cholecalciferol and ergocalciferol are converted to calcidiol by the enzyme vitamin D 25-hydroxylase which is encoded by the

CYP2R1

gene (

Figure 1b), the (rs10741657) SNP described in

2b), the (rs10741657) SNP described in

CYP2R1 had also been related to modulate calcidiol serum status in several studies[12][14][15][18]. Then, after the generation of calcidiol in the liver, it binds again to VDBP and subsequently interacts with the enzyme 1-α hydroxylase, encoded by the

had also been related to modulate calcidiol serum status in several studies [54,56,57,60]. Then, after the generation of calcidiol in the liver, it binds again to VDBP and subsequently interacts with the enzyme 1-α hydroxylase, encoded by the

CYP27B1

gene, mainly in the proximal kidney tubule, calcidiol is converted to calcitriol (1α,25 dihydroxyvitamin D), which is the biological and functional active form of vitamin D (

Figure 1c)[13]. Polymorphisms in

2c) [55]. Polymorphisms in

CYP27B1, particularly the rs10877012 SNP, was associated with lower calcidiol serum levels[19].

, particularly the rs10877012 SNP, was associated with lower calcidiol serum levels [61].
Once the generation of calcitriol is performed, it is bound to VDBP and is transported to target cells and tissues where in cytoplasm or membrane, it interacts with the vitamin D receptor (VDR), encoded by the

VDR

gene.

After the calcitriol translocation to the nucleus, this complex forms a heterodimer with the retinoid x receptor (RXR), and thus regulates vitamin D target genes through its binding to specific DNA sequences called vitamin D response elements (VDREs)[20][21]. Since VDR is expressed in different immune cells, such as neutrophils, macrophages, dendritic cells (CDs), T and B lymphocytes, calcitriol may regulate the immune system[8].

After the calcitriol translocation to the nucleus, this complex forms a heterodimer with the retinoid x receptor (RXR), and thus regulates vitamin D target genes through its binding to specific DNA sequences called vitamin D response elements (VDREs) [62,63]. Since VDR is expressed in different immune cells, such as neutrophils, macrophages, dendritic cells (CDs), T and B lymphocytes, calcitriol may regulate the immune system [8].

Calcitriol can be produced by monocytes and macrophages and generate a shift from pro-inflammatory to tolerogenic immune status[22]. Calcitriol promotes M1 phenotype switching to M2 via the nVDR-PPARγ pathway and via the upregulation of the expression of IL-10[23]. Particularly, calcitriol promotes a shift from Th1 and Th17 to Th2 immune profile via suppression of expression of cytokines of the Th1 (IL-2, IFN-γ, TNF-α) and Th17 (IL-17, IL-21) profiles, and induction of the expression of cytokines of Th2 profile (IL-4, IL-5, IL-9, IL-13), in order to limit inflammatory processes and autoimmune reactions [8][22].

Calcitriol can be produced by monocytes and macrophages and generate a shift from pro-inflammatory to tolerogenic immune status [64]. Calcitriol promotes M1 phenotype switching to M2 via the nVDR-PPARγ pathway and via the upregulation of the expression of IL-10 [65]. Particularly, calcitriol promotes a shift from Th1 and Th17 to Th2 immune profile via suppression of expression of cytokines of the Th1 (IL-2, IFN-γ, TNF-α) and Th17 (IL-17, IL-21) profiles, and induction of the expression of cytokines of Th2 profile (IL-4, IL-5, IL-9, IL-13), in order to limit inflammatory processes and autoimmune reactions [8,64].
Several studies have focused on evaluating SNPs described in the

VDR

, such as

FokI

(rs2228570),

BsmI

(rs1544410),

ApaI

(rs7975232), and

TaqI

(rs731236) polymorphisms, which was associated in several populations with lower calcidiol levels (

Figure 1d). Therefore, SNPs′ presence in key enzymes related to vitamin D metabolism may modulate calcidiol levels and calcitriol function; they may also modify the disease activity in MS, RA, and SLE through the vitamin D deficiency, and contribute to genetic susceptibility to autoimmunity. In the following sections, we will describe how polymorphisms in these genes can modulate their expressions.

2d). Therefore, SNPs′ presence in key enzymes related to vitamin D metabolism may modulate calcidiol levels and calcitriol function; they may also modify the disease activity in MS, RA, and SLE through the vitamin D deficiency, and contribute to genetic susceptibility to autoimmunity. In the following sections, we will describe how polymorphisms in these genes can modulate their expressions.

3.1. Vitamin D Binding Protein (VDBP) (SNP rs2282679

3.1. Vitamin D Binding Protein (VDBP) (SNP rs2282679

GC)

Vitamin D binding protein (VDBP) is the main carrier protein for vitamin D, which binds 85 to 90% of the total circulating calcidiol [24]. This protein is encoded by the

)

Vitamin D binding protein (VDBP) is the main carrier protein for vitamin D, which binds 85 to 90% of the total circulating calcidiol [66]. This protein is encoded by the

GC gene located in chromosome 4, position 13.3 (4q13.3), and has 55,136 bp size [25]. The

gene located in chromosome 4, position 13.3 (4q13.3), and has 55,136 bp size [67]. The

GC gene consists of 13 exons and the rs2282679 SNP was described located at the position 71,742, presenting an A ˃ C change (ancestral allele: A) in the intron 12, near to the actin III subdomain[26] and the endonuclease enzyme commonly used to identify this

gene consists of 13 exons and the rs2282679 SNP was described located at the position 71,742, presenting an A ˃ C change (ancestral allele: A) in the intron 12, near to the actin III subdomain [68] and the endonuclease enzyme commonly used to identify this

GC

(rs2282679) SNP is

FokI

(

Flavobacterium okeanokoites) [27]. According to rs2282679 SNP intron location, its hypothetical functional effect may alter genetic mRNA expression due to modification of the splicing process (44). A complete GWAS showed that

) [69]. According to rs2282679 SNP intron location, its hypothetical functional effect may alter genetic mRNA expression due to modification of the splicing process (44). A complete GWAS showed that

GC (rs2282679) C allele was associated with lower serum calcidiol and VDBP levels in a study carried out in approximately 30,000 subjects with European ancestry included in 15 cohorts[28]. Additionally, various studies have shown the association of the

(rs2282679) C allele was associated with lower serum calcidiol and VDBP levels in a study carried out in approximately 30,000 subjects with European ancestry included in 15 cohorts [70]. Additionally, various studies have shown the association of the

GC (rs2282679) SNP with genetic susceptibility and disease modulation in MS, RA, and SLE[29][30][31][32] (

(rs2282679) SNP with genetic susceptibility and disease modulation in MS, RA, and SLE [71,72,73,74] (

Figure 1a).

2a).

3.2. Vitamin D 25-Hydroxylase (SNP rs10741657

3.2. Vitamin D 25-Hydroxylase (SNP rs10741657

CYP2R1)

The enzyme vitamin D 25-hydroxylase is a protein of around 500 amino acids with a molecular weight of 50–55 kilo-daltons (kDa). The liver is the main site where it is synthesized and performs its function of 25-hydroxylation enzyme activity, besides this enzyme activity was also described in the kidney and intestines[33], but cholecalciferol and ergocalciferol are metabolized mainly by the vitamin D 25-hydroxylase to calcidiol in the liver[34].

The enzyme vitamin D 25-hydroxylase is encoded by the

)

The enzyme vitamin D 25-hydroxylase is a protein of around 500 amino acids with a molecular weight of 50–55 kilo-daltons (kDa). The liver is the main site where it is synthesized and performs its function of 25-hydroxylation enzyme activity, besides this enzyme activity was also described in the kidney and intestines [75], but cholecalciferol and ergocalciferol are metabolized mainly by the vitamin D 25-hydroxylase to calcidiol in the liver [76].
The enzyme vitamin D 25-hydroxylase is encoded by the

CYP2R1

gene located in chromosome 11 in the short arm at position 15.2 (11p15.2) and consists of 15,500 bp. The SNP rs10741657

CYP2R1 is positioned at 14,893,332 pb in the gene, it displays an A ˃ G change in the 5′ UTR region (ancestral allele: A)[18].

In order to recognize the presence of this SNP, the endonuclease restriction enzyme

is positioned at 14,893,332 pb in the gene, it displays an A ˃ G change in the 5′ UTR region (ancestral allele: A) [60].
In order to recognize the presence of this SNP, the endonuclease restriction enzyme

MnlI

(

Moraxella nonliquefaciens) is commonly used[35]. According to the SNP location in 5′ UTR, the rs10741657 SNP may regulate gene expression by modifying translation initialization and transcript stabilization of mRNA, and therefore, modulate 25-hydroxylase expression and enzymatic activity rate[36].

A meta-analysis that included 16 articles with 52,417 participants showed that the

) is commonly used [77]. According to the SNP location in 5′ UTR, the rs10741657 SNP may regulate gene expression by modifying translation initialization and transcript stabilization of mRNA, and therefore, modulate 25-hydroxylase expression and enzymatic activity rate [46].
A meta-analysis that included 16 articles with 52,417 participants showed that the

CYP2R1 (rs10741657) GG genotype shows trends of low calcidiol serum levels compared to the AA genotype in the Caucasian and Asian population[18] and when healthy persons are supplemented with cholecalciferol, carriers of the AA genotype showed a 2.5-fold increase in calcidiol levels compared to GG genotype carriers [12], which highlights the role of

(rs10741657) GG genotype shows trends of low calcidiol serum levels compared to the AA genotype in the Caucasian and Asian population [60], and when healthy persons are supplemented with cholecalciferol, carriers of the AA genotype showed a 2.5-fold increase in calcidiol levels compared to GG genotype carriers [54], which highlights the role of

CYP2R1

(rs10741657) SNP in vitamin D metabolism (

Figure 1b).

2b).

3.3. Vitamin D 1-α Hydroxylase (SNP rs10877012

3.3. Vitamin D 1-α Hydroxylase (SNP rs10877012

CYP27B1)

The enzyme vitamin D 1-α hydroxylase corresponds to a P450 protein of 507 amino acids of around 55 kDa[37]. This enzyme metabolizes calcidiol to calcitriol, the active form of vitamin D[34], and it is encoded by the

)

The enzyme vitamin D 1-α hydroxylase corresponds to a P450 protein of 507 amino acids of around 55 kDa [78]. This enzyme metabolizes calcidiol to calcitriol, the active form of vitamin D [76], and it is encoded by the

CYP27B1

gene, which is found in chromosome 12, long arm, position 14.1 (12q14.1), its size is 6653 bp.
The SNP rs10877012

CYP27B1 is located at the non-coding region -1260 and is characterized by a change of G ˃ T (ancestral allele: G)[38][39]. The restriction endonuclease enzyme

is located at the non-coding region -1260 and is characterized by a change of G ˃ T (ancestral allele: G) [79,80]. The restriction endonuclease enzyme

HinfI

(

Haemophilus influenza) is commonly used for its detection [27][35]. Due to its location in the 3′ UTR region, this SNP may alter transcript stabilization regulation and its localization in the cytoplasm [36].

A study carried out in 253 German patients with differentiated thyroid carcinoma showed that patients with the presence of the GG genotype (Referred as CC in this study by its position in the negative DNA strand) was associated with lower calcitriol serum levels than patients carrying the TT genotype (Referred as AA in this study by its position in negative DNA strand) (60 pmol/mL vs. 72 pmol/mL, respectively)[38].

In another study carried out in Caucasian German patients with gestational diabetes, the GG (CC) genotype was also associated with lower calcidiol levels[19]. Likewise, another study carried out in a healthy Caucasian British population demonstrated an association between the

is commonly used for its detection [69,77]. Due to its location in the 3′ UTR region, this SNP may alter transcript stabilization regulation and its localization in the cytoplasm [46].
A study carried out in 253 German patients with differentiated thyroid carcinoma showed that patients with the presence of the GG genotype (Referred as CC in this study by its position in the negative DNA strand) was associated with lower calcitriol serum levels than patients carrying the TT genotype (Referred as AA in this study by its position in negative DNA strand) (60 pmol/mL vs. 72 pmol/mL, respectively) [79].
In another study carried out in Caucasian German patients with gestational diabetes, the GG (CC) genotype was also associated with lower calcidiol levels [61]. Likewise, another study carried out in a healthy Caucasian British population demonstrated an association between the

CYP27B1 (rs10877012) G (C) allele and lower serum calcidiol levels[40].

Besides, in type 1 diabetes in Caucasian patients from Germany, those carrying the GG (CC) genotype had a reduced amount of mRNA from

(rs10877012) G (C) allele and lower serum calcidiol levels [81].
Besides, in type 1 diabetes in Caucasian patients from Germany, those carrying the GG (CC) genotype had a reduced amount of mRNA from

CYP27B1

compared to HS (1.6855 vs. 1.8107, respectively,

p = 0.0220)[19]. Therefore, the SNP rs10877012

= 0.0220) [61]. Therefore, the SNP rs10877012

CYP27B1

may also modulate the vitamin D serum status and the genetic susceptibility or disease modulation in autoimmune diseases such as MS, RA, and SLE (

Figure 1c).

2c).

3.4. Polymorphisms in Vitamin D Receptor (VDR)

Most genetic studies evaluating the potential association of calcidiol serum levels with genetic polymorphisms have focused on evaluating the polymorphisms described in the vitamin D receptor (VDR). The VDR is a member of the steroid/thyroid hormone receptor superfamily; and this receptor is encoded by the gene with the same name,

VDR

, which is located in chromosome 12, position 12q.13.11, comprising a region of approximately 100,000 bp of DNA, and only 4600 bp encode the VDR protein.

Functionally VDR is a transcription factor regulated by ligand binding and possibly by phosphorylation events[41]. It is a soluble 427 amino acid protein located mainly in the nucleus, cell cytoplasm, and cellular membrane, from where it translocates to the nucleus through the microtubule system after interaction with its ligand, calcitriol [42].

Functionally VDR is a transcription factor regulated by ligand binding and possibly by phosphorylation events [82]. It is a soluble 427 amino acid protein located mainly in the nucleus, cell cytoplasm, and cellular membrane, from where it translocates to the nucleus through the microtubule system after interaction with its ligand, calcitriol [83].
VDR is expressed in various organs involved in calcium metabolism, immune cells, and the nervous system (5). Three isoforms of the VDR were described.

The most common is the VDRA isoform of 427 amino acids and 48 kDa, with a start site in exon 2. The second is a long VDRB1 isoform of 477 amino acids and 54 kDa, this isoform presents 50 amino acids more in the N-terminal domain by an ATG start site in the exon 1d, described in the human kidney as well as in intestinal and renal epithelial cell lines[42][43].

The most common is the VDRA isoform of 427 amino acids and 48 kDa, with a start site in exon 2. The second is a long VDRB1 isoform of 477 amino acids and 54 kDa, this isoform presents 50 amino acids more in the N-terminal domain by an ATG start site in the exon 1d, described in the human kidney as well as in intestinal and renal epithelial cell lines [83,84].
The third is a shorter VDRA isoform of 424 amino acids and greater transactivation capacity as a transcription factor, caused by the SNP

FokI in the exon 2[42].

in the exon 2 [83].
More than 14 different polymorphisms were described in human

VDR

, which could influence the modulation of the response to calcitriol by binding to VDR. The four SNPs most frequently studied are:

FokI

(rs2228570),

BsmI

(rs1544410),

ApaI

(rs7975232), and

TaqI (rs731236) [41][44]. These were related to modulating the vitamin D status independently (in alleles and risk genotypes) as well as in haplotypes and haplogenotypes, both in original articles and meta-analyses (

(rs731236) [82,85]. These were related to modulating the vitamin D status independently (in alleles and risk genotypes) as well as in haplotypes and haplogenotypes, both in original articles and meta-analyses (

Figure 1d).

2d).

3.4.1. FokI (rs2228570) VDR SNP

FokI

(rs2228570)

VDR

SNP, also referred as the start codon polymorphism (SCP), was defined using the

FokI

(Flavobacterium okeanokoites) restriction enzyme in a restriction fragment length polymorphism test (RFLP) [45].

) restriction enzyme in a restriction fragment length polymorphism test (RFLP) [86].

FokI (rs2228570) is located in exon 2 and is considered a non-synonymous polymorphism, because the change of C > T (ancestral allele T), also referred as F > f change, which generates a non-synonymous change of threonine to methionine and dictates two potential translation initiation sites[44].

The presence of the restriction site

(rs2228570) is located in exon 2 and is considered a non-synonymous polymorphism, because the change of C > T (ancestral allele T), also referred as F > f change, which generates a non-synonymous change of threonine to methionine and dictates two potential translation initiation sites [85].
The presence of the restriction site

FokI

is when the C allele is presented (also called F allele by the cut of the

FokI restriction enzyme), this C allele generates a new start codon (ATG) 9 bp after the common starting site, which translates to a shorter VDRA protein of 424 amino acids instead of the wild type full-length VDRA isoform of 427 amino acids[45].

restriction enzyme), this C allele generates a new start codon (ATG) 9 bp after the common starting site, which translates to a shorter VDRA protein of 424 amino acids instead of the wild type full-length VDRA isoform of 427 amino acids [86].

FokI VDR SNP (rs2228570) was found to be functional, and the short 424 amino acid VDRA isoform is somewhat more active than the long VDRA isoform of 427 amino acid, in terms of its transactivation capacity as a transcription factor[46][47]. In the absence of the restriction site

SNP (rs2228570) was found to be functional, and the short 424 amino acid VDRA isoform is somewhat more active than the long VDRA isoform of 427 amino acid, in terms of its transactivation capacity as a transcription factor [87,88]. In the absence of the restriction site

FokI, the T allele (also called f allele), translation begins at the first original site at the exon 2, and the VDRA of 427 amino acids is expressed, which is 1.7-fold less active in its transactivation capacity, and presents less stability[48][49][50][51].

, the T allele (also called f allele), translation begins at the first original site at the exon 2, and the VDRA of 427 amino acids is expressed, which is 1.7-fold less active in its transactivation capacity, and presents less stability [89,90,91,92].

3.4.2. BsmI (rs1544410) and ApaI (rs7975232) VDR SNPs

BsmI

(rs1544410) and

ApaI

(rs7975232)

VDR

SNPs were defined using the

BsmI

(

Bacillus stearothermophilus

) and the

ApaI

(

Acetobacter pasteurianus)

restriction enzyme, respectively in a RFLP test.

BsmI

(rs1544410)

VDR SNP, located in the intron 8, presents a change of A > G (also called B > b), and the ancestral allele is the G allele[49]. Regarding its functional effect, it could generate an alteration in the splice sites for mRNA transcription or a change in the intron regulatory elements of

SNP, located in the intron 8, presents a change of A > G (also called B > b), and the ancestral allele is the G allele [90]. Regarding its functional effect, it could generate an alteration in the splice sites for mRNA transcription or a change in the intron regulatory elements of

VDR

.

ApaI (rs7975232) VDR SNP, also is located in the intron 8, presents a change of A > C, (also called A > a), and the ancestral allele is the C allele [21]. Both SNPs are located at the 3′ end of the

SNP, also is located in the intron 8, presents a change of A > C, (also called A > a), and the ancestral allele is the C allele [63]. Both SNPs are located at the 3′ end of the

VDR

and do not change the amino acid sequence of the VDR protein. Therefore, they could affect mRNA stability and the gene expression of

VDR by LD[21][49].

by LD [63,90].

3.4.3.  TaqI (rs731236) VDR SNP

3.4.3. TaqI (rs731236) VDR SNP

TaqI

(rs731236)

VDR

SNP, located in the exon 9, was defined using the

TaqI

(

Thermus aquaticus

) restriction enzyme, and presents a change of C > T, (also called T > t), and the ancestral allele is the C allele.

TaqI generates a synonym change of the coding sequence; therefore, it does not produce an amino acid change of the encoded protein, but it could influence the stability of the mRNA [49].

generates a synonym change of the coding sequence; therefore, it does not produce an amino acid change of the encoded protein, but it could influence the stability of the mRNA [90].
If

TaqI

(rs731236)

VDR

SNP is in high LD with

ApaI

(rs7975232)

VDR SNP, its functional effect is the possible modification of one of the zinc fingers of the nuclear signaling heterodimer that binds to the VDREs located in the target genes[52].

SNP, its functional effect is the possible modification of one of the zinc fingers of the nuclear signaling heterodimer that binds to the VDREs located in the target genes [93].
Due to the closeness of these four

VDR

SNPs, they were studied to determine their LD. In various populations of SLE, RA, and MS patients, the

BsmI

(rs1544410),

ApaI

(rs7975232), and

TaqI

(rs731236)

VDR SNPs were described with a strong LD[53][54], which infers that the alleles of these three polymorphisms could segregate into haplotypes from one generation of persons to another. In the case of

SNPs were described with a strong LD [94,95], which infers that the alleles of these three polymorphisms could segregate into haplotypes from one generation of persons to another. In the case of

FokI

(rs2228570), it was described in a low LD with the other three

VDR

SNPs, which suggests that

FokI

does not segregate in blocks with others downstream

VDR SNPs [53][54]. However, because the genetic recombination points vary between populations,

SNPs [94,95]. However, because the genetic recombination points vary between populations,

FokI

was also studied in haplotypes together with the

BsmI

(rs1544410),

ApaI

(rs7975232), and

TaqI

(rs731236)

VDR SNPs in several populations [55][56][57].

SNPs in several populations [96,97,98].