1. Introduction
Vitamin D is a steroid hormone of the fat-soluble vitamin class involved in various biological processes, including cell proliferation, bone metabolism, and cell differentiation [
1]; it is also responsible for a complex multi-step metabolism and acts as a hormone in numerous extra-skeletal targets [
2]. Vitamin D deficiency, which occurs when vitamin D levels in the blood are below 20 ng/mL, is prevalent around the world [
3], mainly affecting countries with little sun exposure due to climatic conditions, high latitudes, and winter regimes, in addition to other conditions (e.g., skin hyperpigmentation and chronic diseases) and affecting eating habits, pregnancy, breastfeeding [
4]. This condition is considered a global health problem, as low vitamin D levels are associated with increased risk for various diseases and metabolic changes [
5]. In Brazil, the Brazilian Society of Endocrinology and Metabology recommends serum levels above 20 ng/mL for the general healthy population and 30–60 ng/mL for risk groups such as the elderly, pregnant women, patients with osteomalacia, rickets, osteoporosis, secondary hyperparathyroidism, pre-bariatric patients, and inflammatory, autoimmune, and chronic kidney diseases [
6].
The circulating effects of 1,25(OH)
2D are mediated by the vitamin D receptor (VDR), which is a member of the superfamily of intracellular nuclear receptors [
7]. Vitamin D also has some genomic effects, including cell apoptosis regulation, differentiation, proliferation, DNA repair, oxidative stress, and cell metabolism, which are driven by transcription factors [
1].
Hence, VDR and vitamin D metabolic enzymes are expressed in all innate and adaptive arms of the immune system, and genomic approaches for gene expression profiling have identified various VDR-regulated genes implicated in regulating innate and adaptive immunity, including CYP and manganese-dependent superoxide dismutase (SOD2) genes [
8]. Genetic variations called single nucleotide polymorphisms (SNP) can reach the human genome and alter different genes’ transcription and translation steps. The VDR gene is located on chromosome 12q13.1 and has several SNPs, more notably BsmI (rs 1544410), which is located in intron 8 and the result of substituting an adenine-guanine (A-G) [
9]. Researchers have associated BsmI with osteoarthritis [
10], breast cancer [
11], melanoma [
12], and system lupus erythematosus [
13], among others, although data on its association with serum vitamin D levels are still conflicting [
13,
14].
Other genes are also known to modulate the expression of the VDR gene. For instance, one study found that the CYP gene, whose expression was constant, is regulated by fasting in the liver, while adipose tissue and the brain are the organs where vitamin D seems to play an important, albeit not fully known, role [
3]. Additionally, the SOD2 gene has an antioxidant function, which may be relevant at low serum vitamin D levels [
15]. Therefore, it is important to shed more light on how the different genotypes of the SNP BsmI of the VDR gene can modulate the expression of the genes of the VDR itself and others that may be in its pathway.
Given the above, this study sought to evaluate the effect of the genotypes of SNP BsmI of the VDR gene on the expression of different genes (VDR, SOD2, and CYP24A1) in individuals with low serum vitamin D levels.
3. Results
Nine-eight individuals with hypovitaminosis D (≤20 ng/dL), with a mean age of 30.5 ± 11 years, 54.1% female, and self-declared white (54.1%), participated in the study.
Table 2 lists the baseline characteristics of the total sample.
Table 3 lists the descriptive analyses of anthropometric and biochemical measurements.
The genotypic and allelic frequencies of the BsmI polymorphism of the VDR gene in the total sample are described in
Table 4. The studied sample presented Hardy-Weinberg equilibrium (χ
2 =2.01).
Comparisons between genotype groups are listed in
Table 5, and no differences were found between the markers analyzed among the different BsmI genotypes.
The VDR, SOD2, and CYP24A1 expressions were performed in each genotype group for the VDR SNP BsmI, and the results are illustrated in
Figure 2; VDR gene expression was significantly lower in the GA and AA genotypes than in the wild-type GG genotype, and AA showed the lowest expression (
Figure 2A). The SOD2 gene was also significantly less expressed in the AA group than GG and GA, which showed no differences (
Figure 2B). In contrast, CYP24A1 overexpression was observed in the GA and AA genotypes (
Figure 2C).
4. Discussion
This study demonstrated, for the first time, the effects of different BsmI genotypes of the VDR gene on the expression of different genes (besides VDR itself) in patients with hypovitaminosis D. BsmI downregulated the VDR gene and SOD2 and overregulated CYP24A1, thus showing the modulating influence of BsmI on other genes in hypovitaminosis D.
The sample represents general characteristics observed in young and active individuals, most of whom were women, as observed in other studies (
Table 2). The higher self-declaration of color (white) is due to southern Brazil (i.e., the study region) presenting a higher prevalence of self-declared white individuals [
20]. As for the level of education, most individuals had incomplete higher education, implying that they could understand the questions of the questionnaire, different from the data found in the rest of the country where the proportion of people aged 25 years or older are those who have completed compulsory basic education (i.e., high school) [
21]. Indeed, another study carried out with individuals with vitamin D insufficiency and assisted by primary health care services reported serum vitamin D levels not being associated with the level of education, age, marital status, and income [
22].
Another important finding is regarding the participants’ being overweight (
Table 3), as the average BMI was within 27.1 kg/m
2. A recent study evaluated serum vitamin D levels in healthy adult women and also showed that the nutritional status of the group evaluated was classified as overweight/obese, although without showing a significant association [
23]. In this study, the distribution of genotypes and alleles of SNP BsmI were in Hardy-Weinberg equilibrium, and there were no significant associations between BsmI genotypes and serum vitamin D levels or the other markers (
Table 5), as previously demonstrated by Retamoso and collaborators [
15]. In other reports, no significant differences between vitamin D levels or genotypic and allelic frequencies of polymorphisms in the VDR gene were found, even though cross-sectional studies have shown that VDR gene polymorphisms can reduce the affinity of the VDR for serum vitamin D levels [
24]. Despite VDR genetic polymorphisms being determinants of vitamin D levels, they have other genetic and environmental factors that are influenced by sun exposure, diet, and even skin pigmentation [
25].
In order to clarify whether the BsmI SNP may influence the expression of other genes related to the physiological role of vitamin D, VDR, CYP24A1, and SOD-2 gene expression was studied in the groups of the three genotypes. Thus, we know that the VDR can modulate the expression of various genes, and its inactivation occurs due to the lack or excess of vitamin D since the almost ubiquitous expression of the VDR gene supports corroborates data from the last 30 years showing that vitamin D not only regulates calcium homeostasis but also promotes immunity, growth, and cell differentiation [
26].
The VDR gene was under-expressed in the GA and AA genotypes compared to GG and AA less than the others (
Figure 2A). In another study that compared the frequency of GG versus AA and AG genotypes, the association with insufficient 25(OH)D concentrations was maintained, suggesting that BsmI, which regulates VDR expression, can modulate vitamin D levels in patients with cognitive disorders [
24]. It is important to emphasize that this study did not evaluate specific pathological conditions but individuals with hypovitaminosis D, which may explain these findings since they were different groups (people with hypovitaminosis D against a population with cognitive impairment) or the degree of genetic mixture of the population studied considering that they are of different ethnic origins and have a high degree of miscegenation among the populations investigated.
Moreover, the SOD-2 gene was significantly less expressed in AA genotype carriers compared to GG and GA (
Figure 2B). This result can be explained by the study of
Dauletbaev and collaborators [28], who investigated the impact of the genome on transcription by 1,25(OH)
2D in carcinogenic cells, in which 1,25(OH)
2D induced genes such as SOD2, IRS2, BIRC3, and DUSP1/5, which are cytoplasmic or mitochondrial signaling molecules that mediate the effects of growth factors and/or cytokine interactions with known anticancer properties. In other words, 1,25(OH)
2D significantly induced mitochondrial expression but not cytosolic SOD2 that converts the free radical O2
•− (superoxide) into H
2 O
2 to defend against free radicals [
8].
Lastly, we evaluated whether the BsmI SNP of the VDR gene would modulate the expression of the CYP24A1 gene, which showed overexpression in carriers of the GA and AA genotypes compared to GG (
Figure 2C), with AA being significantly more expressed than GG and GA, in which the mutated allele “A” possibly increases CYP24A1 expression. This is an unprecedented analysis in the literature, and the expression of most vitamin D target genes is 5x up- or down-regulated, meaning only a few genes respond with significant changes in vitamin D expression [
27], as is the case of CYP24A1, which showed higher expression against BsmI SNP of the VDR gene.
Despite this study not evaluating the individuals by racial group, which was a limiting factor, it was possible to achieve unpublished results in which the SNP BsmI of the VDR gene can modulate CYP24A1 and SOD2 expression, even if it did not influence serum vitamin D levels. Given this context, future studies can seek to shed more light according to skin color.