In the present study, a high prevalence of vitamin D deficiency or insufficiency was found in 7 to 18-year old girls from South Brazil. In addition, VDR wild-type genotypes BsmI, ApaI and TaqI were associated with lower 25(OH)D levels both individually and as the GGT haplotype, even after adjustment for age and season at the time of blood collection, and were independent positive predictors of 25(OH)D levels. To the best of our knowledge, this is the first report relating VDR gene polymorphisms and haplotypes to 25(OH)D levels in healthy children and adolescent girls from the general population.
Previous studies have reported low levels of 25(OH)D in different pediatric and adolescent populations [
2,
3]. In our sample, only 9.4% of the participants had adequate circulating concentrations of 25(OH)D, which is far below the 39.0% reported for a 16–20

year old population from a rural town in the state of São Paulo, Brazil (latitude −23°) [
17]. Since the latitude and solar radiation are similar in the cities included in those previous studies and in the cities studied by us, diverse dietary habits and the age of participants may explain the discrepancies in the frequency of sufficient 25(OH)D status. In this sense, one limitation of the present study is the lack of control for some confounders, such as dietary vitamin D intake and duration of sun exposure. However, no associations were found between 25(OH)D levels and season at the time of blood collection. One potential explanation for the lack of seasonal difference is that below a latitude of approximately 35°, UVB radiation is sufficient for vitamin D synthesis all year round [
18].
Concerning young populations, several studies have shown that vitamin D levels are lower in the presence of overweight/obesity, and that low levels of 25(OH)D could influence the risk of developing metabolic disorders and cardiovascular disease in pediatric and adolescent populations [
4-
6]. Obesity associated-vitamin D insufficiency is likely caused by deposition of skin and dietary vitamin D3 in body fat compartments, resulting in decreased bioavailability [
19].
In addition, studies have shown that obese girls with vitamin D deficiency presented lower insulin sensitivity [
4], and adolescents with reduced levels of vitamin D had increased risk of metabolic syndrome [
5]. Moreover, a recent study with children aged 8 to 18

years, in which 47% were obese, reported a negative association between BMI, percentage of body fat, subcutaneous and visceral adipose tissue and 25(OH)D levels in white and black subjects [
3]. Despite the age similarity, the frequency of obesity was very low in our sample, and therefore we were unable to assess these associations.
Only a few studies are available in the literature relating VDR polymorphisms and 25(OH)D levels, especially in young and healthy populations. In a study with adolescent girls, the BsmI polymorphism was marginally, but not significantly, associated with 25(OH)D levels [
20]. Regarding adult populations, the results are conflicting. In India, one study reported that TaqI polymorphism was associated with higher 25(OH)D levels in healthy adults aged 25 to 60

years [
13], while another study analyzing BsmI and TaqI SNPs did not confirm this association [
21]. In addition, other SNPs in different regions of the VDR gene have been associated with levels of 25(OH)D in adolescent girls [
20] and other populations [
14].
In contrast, studies focusing on specific conditions, such as the overweight Insulin Resistance Atherosclerosis family study, in which only the BsmI SNP [
22] was analyzed, or the Metabolites in Multiple Sclerosis study, covering ApaI and TaqI SNPs [
23], did not find an association between VDR polymorphism and 25(OH)D levels. However, these studies were carried out with disease-affected adult populations, whereas our results refer to younger, healthy subjects. In studies with older overweight subjects, BsmI, ApaI, and TaqI SNPs of VDR gene were not associated with 25(OH)D levels. Nevertheless, the lack of association may be explained, at least in part, by the fact that most of the participants in that study were using vitamin D supplements [
24,
25].
VDR gene polymorphisms have also been studied in some other conditions. BsmI, ApaI, and TaqI polymorphic genotypes were associated with suppressed cytokine response in pulmonary tuberculosis [
26] and appear to be related with higher glucose levels [
12] and diabetes [
27,
28]. ApaI gene variant has been linked to higher bone metabolism and bone mineral density [
10]. We also have recently described that the wild genotype of the ApaI polymorphism is more frequent in girls with precocious pubarche than in controls with normal pubertal development [
8].
The BsmI, ApaI and TaqI polymorphism are located at the 3’ untranslated region (3′ UTR) of the VDR gene. The BsmI and ApaI SNPs are located in intron 8, and TaqI is a silent SNP in exon 9. Although 3’ UTR has been recognized as a region involved in the modulation of gene expression, especially through the regulation of mRNA stability and efficiency of protein translation [
29], the functional role of the BsmI, ApaI, and TaqI VDR gene polymorphisms has not been established. Moreover, the most common haplotype for the VDR gene in the present study (Ht1) is the same found in Caucasians and Asians, followed by Ht7 and Ht3, as found for Caucasians [
30]. A strong linkage disequilibrium between these three SNPs has been shown in the present and previous studies in different populations [
31].
In the present study, girls with the wild genotype for BsmI and TaqI SNPs or presenting the Ht1 haplotype were found to be at higher risk for 25(OH)D deficiency. In turn, Ht7 was regarded as the protective haplotype for circulating 25(OH)D levels. Therefore, it is possible to hypothesize that the VDR genotype influences the susceptibility to vitamin D deficiency in children and adolescents. Thus, the present results open new possibilities for further studies in this field, especially on the genotype-related mechanisms of calcium and vitamin D metabolism and prevention of 25(OH)D deficiency in young subjects.