Evolution of Vitamin D Status and Vitamin D Receptor Gene Expression Among Professional Handball Athletes During a Competitive Period. Relationship with Body Composition, Calcium, Magnesium and Phosphorous

Introduction: A generalized risk of vitamin D deficiency exists worldwide affecting also professional and elite athletes. This study assesses the evolution of vitamin D status and vitamin D receptor (VDR) gene expression and their relationship with body composition, calcium (Ca), magnesium (Mg) and phosphorous (P) among professional handball athletes during a competitive period. Methods: A total of 26 male subjects were recruited: 13 professional handball athletes and 13 non-athlete controls. An observational follow-up study was conducted in 2 time points over a 16-week period. Nutritional intake, body composition, and routinary biochemical parameters were measured via 24-hours recall, bioimpedance and enzyme immunoassay, respectively. Ca and Mg were measured by flame atomic absorption spectrophotometry and P was determined with the colorimetric method of Fiske-Subbarow. 25-hydroxyvitamin-D (25(OH)D) levels and its forms (i.e., 25(OH)D3 and 25(OH)D2) were measured by liquid chromatography-tandem mass spectrometry (LC-MS/MS), whereas VDR gene expression was measured by quantitative real time-polymerase chain reaction (qRT-PCR). Results: A total of 54% of the athletes showed deficient vitamin D status. Moreover, a prevalence of insufficient vitamin D status in handball players affected 46% at baseline, reaching 61% after 16 weeks. Vitamin D showed no evolution during the competitive period and no differences between groups were observed (all p ≥ 0.05). Handball players increased the VDR expression, enhanced body composition, Ca and Mg levels at 16-weeks follow-up (all p < 0.05). VDR gene expression was positively related with body mass and body mass index at follow-up in athletes (all p ≤ 0.038; r ≥ 0.579) and with Ca at baseline in controls (p = 0.026; r = 0.648). Finally, 25(OH)D2 form was directly associated with P in athletes at 16 weeks of study (p = 0.034; r = 0.588). Conclusion: Players of indoor team sports such as handball would be a population at risk of vitamin D deficiency. The 16-weeks competition improved VDR gene expression, body composition, Ca and Mg levels. The associations observed between VDR gene expression and the variables of the study evidenced the importance of this receptor as a marker involved in health status in handball athletes despite vitamin D − although in a deficient status −, Ca, Mg and P showed no remarkable changes during the competition period.


Introduction
Vitamin D is a pro-hormone which comprises the combination of 25-hydroxyvitamin D 3 (25(OH)D 3 ) and 25-hydroxyvitamin D 2 (25(OH)D 2 ) thus forming 25-hydroxyvitamin D (25(OH)D) [1].Vitamin D status has become a subject of growing interest, existing a generalized deficiency worldwide, especially in at-risk of deficiency groups [2], but also affecting both professional and elite athletes [3,4].Moreover, the novel identified functions of vitamin D during recent years (e.g., optimizing muscle function, maintaining bone health, or minimizing the risk of infection) have underscored its relevance and potential benefits in the world of sports nutrition [4].
The high incidence of poor vitamin D status on athletes became clear in a recent meta-analysis where 56% from a total of 2313 of elite and sub-elite athletes (considering multiple disciplines' across the world) presented vitamin D insufficiency [3].In this respect, previous studies performed in professional sports team's athletes have reported vitamin D levels to be below reference values in approximately 80% of basketball [5] and 70% of football players [6,7].Moreover, the difference is notable between outdoor and indoor athletes, with a prevalence of vitamin D deficiency of 59% in outdoor athletes and 64% in indoor athletes [8].A large heterogeneity exists additionally in the literature regarding vitamin D status according to multitude of variables (i.e., cohorts of athletes, individual characteristics, type and level of sport practice, geographical location and latitude, methodology for vitamin D determination, or clinical criteria to establish vitamin D status' categories) [9].Thus, the evidence about vitamin D status in team sports (such as handball) remains scarce at the present time, highlighting the importance of testing and correcting vitamin D status in athletes, considering the wide variability of sports disciplines and the potential benefits of an adequate vitamin D status [10], especially when a possible protective effect and enhancement of physical performance is considered [11].
Vitamin D exerts its action through binding to the vitamin D receptor (VDR) [12], which is present in almost all body tissues [13].Therefore, not only an adequate level of vitamin D may be of relevance to athletes health, but also VDR could play an important role [14] via modulating immunity [15] or cardiovascular function [16].Previous evidence has suggested increased vitamin D levels and VDR up-regulation in muscle cells to exert a direct influence upon the efficiency of calcium (Ca) binding for muscle fiber contraction [17], being a possible crucial factor for sports athletes and their exercise performance [18].
To the best of our knowledge, only a few studies [19,20] have determined vitamin D status using the "gold standard" based on liquid chromatography-tandem mass spectrometry (LC-MS/MS) [21], and none of them complemented its determination with an evaluation of VDR gene expression in whole blood.Further, methodological and well-designed studies including healthy control groups are lacking, not existing evidence in the case of professional handball players [22,23].Thus, the present study was aimed to: (I) evaluate the status of both circulating vitamin D and its forms through LC-MS/MS in a group of professional handball players, (II) assess VDR gene expression, and (III) determine their relationship with body composition and minerals involved in vitamin D metabolism (i.e., Ca, magnesium (Mg) and phosphorous (P)) during the competitive period.

Study Design and Participants
This observational study was conducted in 26 male subjects divided in 2 groups: handball athletes (n = 13; aged 22.9 ± 2.7 years), and non-athlete healthy controls (n = 13; aged 20.9 ± 2.8 years).All athletes, who participated in the Spanish Professional Handball League, were monitored over a period of 16 weeks, with 2 evaluation timepoints established: baseline (i.e., initial conditions) and followup (i.e., after 16 weeks of training and competition).The study was conducted during the first competitive period of the handball season in which the athletes played an accumulated total of 12 games from September to January.The training duration of the athletes consisted of 10.5 ± 2.1 h/ week of a handball training regimen, in addition to competition matches on weekends.The control group did not meet the current recommendations for physical activity in adults as follows: at least 150-300 min/week of moderate intensity or 75-150 min/week of vigorous intensity aerobic physical activity, or an equivalent combination of moderate-and vigorous-intensity aerobic physical activity [24].
The inclusion criteria for handball athletes were: (I) to be experienced athletes who had been training for 8-12 years, (II) to pass a recruitment medical evaluation consisting of a clinical examination, (III) not to present injuries during and at least 6 months before the study, (IV) non-smoker status, and (V) the absence of medications use.The exclusion criteria were: (I) to consume nutritional supplements of any kind in the 6 weeks before and during the study, and (II) to use sunscreen.The present study was approved by the Ethics Committee of the University of Granada (Granada, Spain), and was conducted in accordance with the last revised guidelines of the Declaration of Helsinki [25].The purpose of the study and its risks and discomforts were explained to the participants before their written consent was obtained.

Anthropometric and Body Composition Analysis
Height (m) was assessed with a stadiometer (Seca, model 213, range 85 to 200 cm; precision: 1 mm; Hamburg, Germany).Body composition was obtained by multi-frequency bioelectrical impedance (Tanita MC-980MA, Barcelona, Spain).The equipment complies with the applicable European standards (93/42 EEC, 90/384 EEC) for use in the medical industry.The participants were informed in advance of the conditions that had to be observed before the measurement: no alcohol consumption (24 h), no vigorous exercise (12 h), no food or drink (3 h), and no urination.All measurements took place in the morning at the same time.The following variables were obtained: body mass (BM) (kg), and fat mass (FM) (%).Body mass index (BMI) was calculated as weight (kg) divided by the square of height (m 2 ).

Dietary Intake
The assessment of energy and micronutrient intake was quantitatively performed by means of a 24-hour dietary recall questionnaire by a qualified dietitian.Data from food intakes were obtained during individual interviews to request information from each participant about the types of foods and serving sizes.Recall accuracy was recorded with a set of photographs of prepared foods and dishes that are commonly consumed in Spain.The validated Nutriber® software package was used to estimate the intake of each nutrient for the individual athletes and non-athlete controls [26].Energy intake was represented both in absolute values and as an energy per weight ratio.Macronutrient intake was expressed as percentage of energy supplied by each nutrient per total energy ingested.Micronutrient intake in turn was stated as nutrient density (i.e., the mass of micronutrient per 1000 kcal).All subjects were required to maintain their dietary habits over the study period.

Biochemical Parameters
All participants' samples were obtained on Monday morning between 8:00 and 10:00 under fasting conditions and after abstaining from physical exercise for at least 12 h.Whole blood was drawn from the antecubital vein and plasma was separated by centrifugation at 4 °C for 15 min at 3000 x g.Samples were frozen until further analysis.

RNA Extraction and VDR Expression
Whole blood samples were collected in PAXgene® Blood RNA Tube for blood collection, stabilization, and transport, and total RNA was extracted using RNA Isolation Kit (PAX-gene®, ref. 762,165, Becton Dickinson, Germany) following the manufacturer's instructions.Quality control of samples was performed by adding 1 µL of RNA in the spectrophotometry (Nanodrop 2000c-Thermo Scientific, Wilmington, Delaware, United States) for nucleic acids.The absorbance was measured at 260 nm (A260) and 280 nm (A280) and the ratio A260/A280 was used to assess the purity of RNA which was accepted when the value was ≥ 1.9.RNA integrity was assessed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, United States of America) and RNA quality was rated according to the RNA integrity number (RIN).The high-quality samples (RIN ≥ 9) and standard quality samples (RIN = 7) were included in the study to secure a more realistic approach.
A quantity of 2 µg of the total RNA with random hexamers was used to cDNA synthesis performed by High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems™ ref. 4,374,966).Thermal cycler conditions were set to run for 2 h at 37 °C and for 10 min at 75 °C, followed by rapid cooling to 4 °C and cDNA samples were stored at − 20 °C until assays.For the determination of gene expression of the VDR, the relative gene expression was determined in triplicate by the Taqman™ qPCR method (Hs00172113_m1, Applied Biosystems; Foster City, California, United States of America), using Quant Studio 12 K Flex Real-Time PCR

Results
The descriptive characteristics of the participants of the study are represented in Table 1.Handball players presented greater height and BM, and lower percentage of FM compared to non-athletes (all p ≤ 0.001; ES ≥ 1.41).Despite energy intake was higher in athletes compared to controls, fat consumption was greater in the control group (all p ≤ 0.038; ES ≥ 0.881).Finally, athletes showed higher plasma creatinine, urea, and triglycerides levels than controls (all p ≤ 0.010; ES ≥ 1.10).
Figure 1 shows the prevalence of vitamin D status (i.e., sufficiency, insufficiency, and deficiency) and its evolution throughout the study period.A total of 54% of the athletes showed deficient vitamin D status, decreasing to 39% at follow-up.The prevalence of insufficient vitamin D status in handball players affected 46% at baseline, reaching 61.0% after 16 weeks.No significant differences were observed between athletes and control group regarding vitamin D status (p ≥ 0.05).
Figure 2 shows the VDR relative mRNA expression in whole blood.The results showed a significant increase in VDR expression in handball players at 16-weeks follow-up (p = 0.019), presenting up-regulation (p < 0.002) in comparison with the controls.No inter-groups differences were observed at baseline (all p ≥ 0.05).
The comparative and evolutionary analysis of body composition, vitamin D and its forms, and related minerals levels in the participants of the study are shown in Table 2. Athletes had a higher BM and lower FM throughout the study (all p ≤ 0.005; ES ≥ 1.26) compared to controls.All the study subjects were within normal ranges for minerals, however, there were between-groups differences, with higher plasma Mg (p = 0.046) at baseline and lower Ca levels (all, p < 0.001) at 16 weeks of study in our athletes relative to non-athletes.Moreover, athletes showed lower levels of 25(OH)D 2 form at baseline compared to non-athletes (p = 0.026).Throughout the study, an increase in both Mg and Ca levels was observed in non-athletes, and a reduction in Mg levels for athletes (all, p ≤ 0.046).Regarding the evolution of vitamin D concentrations in both study groups, no significant differences were found in mean levels of total 25(OH)D and its forms throughout the study period (all p ≥ 0.05).
The relationship among 25(OH)D, 25(OH)D 2 , 25(OH)D 3 levels and VDR gene expression with body composition and the minerals of the study in both non-athletes and handball athletes are represented in Table 3.The relative VDR gene expression was directly associated with BM and BMI (all p ≤ 0.038; r ≥ 0.579) at follow-up in athletes.VDR gene expression was additionally directly related to Ca levels in non-athletes controls at baseline (p = 0.026; r = 0.648).
System (Life Technologies, Carlsbad, California, United States).A 1.2 µL volume of cDNA-or DNase-free water was mixed with 3.8 µL of TaqMan Gene Expression Master Mix (Applied Biosystems; Foster City, California, United States of America).5 µL of PCR Mix was accurately added to each subarray using AccufillTM (Applied Biosystems; Foster City, California, United States of America).The Open Array system conducts each PCR reaction in 33 nL of reagents.Standard cycling conditions were used as recommended by the manufacturer.The selected housekeeping genes were previously validated using qBasePLUS software (Biogazelle, Gent, Belgium) [30].The selected housekeeping genes were: ATP5B (ATP synthase H + transporting, mitochondrial F1 complex, beta polypeptide); LUC7L2 (LUC7 like 2, pre-RNA splicing factor); ARF1 (ADPribosylation factor 1); and PAPOLA (Poly(A) polymerase alpha).The mRNA relative analysis expression was performed through the 2 −∆∆Cq method [31] normalized to the selected housekeeping genes.

Statistical Analysis
Statistical analyses were performed using the SPSS version 25.0 statistical package for MS Windows (SPSS, Inc., Chicago, Illinois, United States of America).GraphPad Prism 9.0 software (GraphPad Software, San Diego, California, United States of America) was used for plotting the graphs.Data were expressed as means (X) and standard deviations (SD).Data normality was checked using the Shapiro-Wilks test.The Student t-test for independent samples was performed to determine group differences between handball athletes and non-athletes.The paired t-test was used to determine the evolution of body composition and biochemical parameters of the study.Moreover, the Mann-Whitney test was utilized to determine the difference in independent group and the Wilcoxon test for dependent samples to verify intragroup differences for VDR gene expression.Group distributions for 25(OH)D were analyzed through chi-square test.Pearson's correlation coefficient was employed: (I) to establish correlations between vitamin D and its forms with the body composition and biochemical variables throughout the study period; (II) to evaluate the relationship between changes in VDR gene expression and changes in the variables of the study over time.Effect sizes (ES) between the non-athletes and athletes were calculated using Cohen's d, and interpreted as small (0.20-0.50), moderate (0.50-0.80) or large (> 0.80) -ES below 0.2 being considered trivial - [32].Differences were considered significant at p-Values < 0.05. of insufficient and deficient vitamin D status, affecting all athletes at baseline and follow-up -25(OH)D levels not improving during the evaluation period -.Moreover, an up-regulation of VDR expression was observed after 16 weeks of study in handball players.Handball players, in comparison to non-athletes, had higher BM, lower FM, decreased Ca and increased Mg levels throughout the study period.Additionally, Ca levels at baseline were associated to VDR gene expression in controls, whereas VDR gene expression in athletes was directly related to BM and BMI after 16-weeks period, also observing a positive correlation between 25(OH)D 2 and P levels at this point.These outcomes suggest a potential role of VDR gene expression as useful biomarker in blood which, despite the non-significant relationships observed for vitamin D levels, could provide additional information of different health variables in indoors handball athletes.Exercise seemed to be an external Finally, 25(OH)D 2 form was positively associated with P levels in athletes at 16 weeks of study (p = 0.034; r = 0.588).The relationship between changes in VDR gene expression and biochemical and body composition variables was further explored throughout the study.However, no significant relationships were observed between them (all p ≥ 0.05) (data not shown).

Discussion
To our knowledge, this is the first study reporting vitamin D status in professional handball players, providing data of VDR gene expression, and determining their relationship with body composition and minerals involved in vitamin D metabolism during the competitive period.The main findings of the present research showed a high prevalence The quantitative variables data were expressed as the mean and standard deviation (SD).The percentage of macronutrients (carbohydrates, proteins, and fats) were represented as percentage of the total energy intake.The nutritional intake of vitamin D, magnesium, calcium and phosphorous were showed as the micronutrient intake per 1000 kcal of energy ingested.The effect sizes (ES) were calculated using the Cohen's comparison of groups of equal size, and the 95% confidence limit (CL) were calculated for all dependent variables.References values: Glycaemia, 70-110 mg/dL; Creatinine, 0.7-1.2mg/dL; Urea, 10-50 mg/dL; Uric acid, 3.4-7 mg/dL; Triglycerides, 50-200 mg/dL; HDL, 40-60 mg/dL; LDL, 70-150 mg/dL; Cholesterol, 110-200 mg/dL; Prealbumin, 20-40 mg/dL; Albumin, 3.5-5.2mg/dL and according to our results, indoor handball athletes present a considerable risk of vitamin D deficiency.Moreover, Bauer et al. [22] observed that vitamin D deficiency is a common finding in professional handball athletes even in summer, suggesting a negative effect on their performance and the need of interventions for reverting that situation.Additionally, the comparative analysis of 25(OH)D forms showed higher 25(OH)D 2 levels in non-athletes at baseline.Although the biological impact of this form has been observed to be less pronounced, assessing their levels could also be of relevance, since previous evidence indicates that the entry of vitamin D 2 into the total body pool of vitamin D dilutes the relative amount of vitamin D 3 [37].Moreover, in this sense, the non-collection of sun exposure precludes us to establish such relationship.
Interestingly, and even though plasma vitamin D levels remained constant throughout the study period, VDR gene expression was upregulated after 16-weeks only in athletes' group -VDR gene expression being directly related to BM and BMI -suggesting a possible response to exercise performance as external stimulus.VDR has been mainly studied in tissues such as skeletal muscle, given its relevance as a target tissue in sports performance [4,38], being described as a robust marker of the hypertrophic response to resistance exercise in humans [39], so these changes could have been a consequence of the adaptations to exercise itself.Nevertheless, it has been postulated that an adequate concentration of vitamin D is necessary to optimize the gene modulation role of VDR in blood [40], affecting a number of physiological stimulus for VDR gene expression, since it was up-regulated in athletes' group at follow-up after 16 weeks.
The prevalence of vitamin D deficiency was around 50% in both groups at baseline -the percentage of insufficiency being around 30-40% at this point -.Interestingly, sufficient vitamin D status was not observed in athletes in any evaluation point, being the majority of control group' participants also below the limits for sufficiency [28,29].Many factors, including insufficient vitamin D intake or sun exposure have been implicated in decreased vitamin D status [9].Although we did not assess qualitatively the vitamin D intake by means of a food frequency questionnaire (FFQ), the lack of reported dietary sources of both vitamin D 3 (e.g., blue fish, meat livers, milk, and eggs) [33] and D 2 (e.g., juices, margarines, and mushrooms) [34] may explain the poor vitamin D status observed throughout the study period.This study was carried out during autumn and early winter, which could partly explain the low vitamin D levels obtained.Sunlight limitations due to geographical location, lifestyle, or the performance of most indoor activities (including sports practice), might have favored such general vitamin D deficient status [35].In European countries, it has been observed higher and lower prevalence rates of vitamin D deficiency in general population during the typical winter (January to March) and summer (July to September) months, respectively [36].Additionally, the prevalence of vitamin D deficiency in athletes from different disciplines and countries is prominent, increasing this risk in periods with fewer hours of sun exposure per day [3].Therefore, Fig. 1 The status of vitamin D in handball athletes (n = 13) compared to non-athletes (n = 13) throughout the study period.Data are represented as frequencies in percentage.The status of vitamin D was classified according to the reference values of vitamin D in plasma: sufficiency > 30.0 ng/mL, insufficiency 20.0-30.0 ng/mL, and deficiency < 20.0 ng/mL for 25(OH)D levels [28,29] In our study, it was observed that both athletes and healthy controls presented biochemical parameters with clinical-nutritional relevance within the reference values for healthy population.In this regard, previous research has suggested that vitamin D deficiency may be associated with lower albumin and/or pre-albumin levels, indicating a possible link between vitamin D status and protein malnutrition [41].Also, circulating Ca, Mg and P were within the reference values for healthy population throughout the study period, indicating that despite the observed high prevalence of insufficiency and deficiency in vitamin D, the analyzed biochemical parameters were not affected and consequently their general health status.However, Mg levels were higher at baseline and Ca levels were lower at 16 weeks of study in athletes than in healthy controls.Despite being in the competitive period, our athletes showed greater stability in Ca, Mg and P concentrations than the non-athlete group − the differences being mainly due to variations in the latter group −.
Vitamin D acts in the body in 2 different ways, through endocrine and autocrine mechanisms, in which the functions important for health, training, and recovery.In fact, it has been suggested that maintaining adequate levels of vitamin D may improve muscle contractility and reduce the risk of muscle injury or may be involved in muscle repair and recovery processes [9].Likewise, it has been also reported that vitamin D deficiency would lead to dysregulation of VDR [4].Therefore, and considering that our population consisted of professional athletes and the analyzed study period was the competitive season, such up-regulation of VDR gene expression could have been related to the vitamin D deficient status or, ultimately, the exercise-promoted adaptations.The absence of a relationship with the percentage of body fat in our athletes and, on the contrary, the association observed with BMI, leads us to consider that these changes could be related to the predominance of muscle mass in our athletes, as BMI depends too on muscle mass.More research is needed to elucidate the role of different variables modulating VDR gene expression and its relationship with other parameters which could affect athletes' health and, therefore, their sport performance.Mann-Whitney U test was used to compare differences between independent groups (Non-Athletes vs. Handball athletes).Wilcoxon test was used to compare intra-group values (Baseline vs. 16 Weeks Follow-up).* Statistical significance was considered for p < 0.05 aspects related to cell cycle [42].Despite the high prevalence of vitamin D insufficiency and deficiency observed in our athletes, no relationships were observed with the analyzed minerals except for P, which could lead us to indicate that exercise may modulate Ca, Mg and P levels, independently of the vitamin D status.The positive association observed between the 25(OH)D 2 form and plasma P concentrations in our handball players may be explained by the effect which vitamin D presents upon (I) increasing P intestinal absorption and (II) facilitating P renal reabsorption, therefore above-mentioned minerals may play a key role.On the one hand, vitamin D is essential for bone growth, density, and remodeling, and low levels of vitamin D may increase the risk of fracture [10].On the other hand, from an endocrine perspective, it is involved in essential bodily processes such as response to gene expression, changes related to protein synthesis, hormone synthesis, immune/inflammatory response, as well as cell turnover and synthesis, in which vitamin D inadequacy availability may limit compromise both the performance of the tissues and their control of  importance of this receptor as a marker involved in health status in handball athletes despite vitamin D − although in a deficient status −, Ca, Mg and P showed no remarkable changes during the competition period.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.
elevating its blood concentrations and assuring an adequate bone mineralization [43].Despite the non-significant correlations observed between total 25(OH)D and plasma Mg or Ca levels, it is important to control these minerals, since an optimal Mg or Ca status may be important for optimizing total 25(OH)D status and its relation to physical sport performance aspects (i.e., ATP regeneration, protein synthesis) [44] and bone health [45].
As a practical application and due to inadequate vitamin D levels can lead to decreased athletic performance, increased risk of injury and poor recovery, it is crucial for coaches and athletes to assess and address vitamin D deficiency.Moreover, athletes should be encouraged to consult with a health professional, such as a physician or sports medicine specialist for a comprehensive assessment, as individual requirements may vary depending on factors such as ethnicity, age and health status.Also, it is advisable to periodically re-evaluate vitamin D levels in athletes, especially in the off-season, during the winter months or in indoor sports when sun exposure may be limited.This helps monitor progress and adjust interventions accordingly, being necessary to consult with a healthcare professional to determine the appropriate dosage based on individual needs and test results.
The present study has a number of limitations including: (I) the hours of sun exposure were not recorded, so we did not control for the confounding effect of this variable upon our results; (II) the small number of participants limited the power of the observed relationships between the variables analyzed in the study; (III) the focus on a single discipline of team sports and in male athletes in turn makes it difficult to extrapolate our results to other sports disciplines or female athletes; (IV) The qualitative dietary intake was not assessed through the utilization of a FFQ, which could have provided elucidation regarding the suboptimal vitamin D status that was observed, potentially resulting from inadequate consumption of vitamin D-rich food sources.Nevertheless, an observational study has been made of professional athletes with a control population; vitamin D was determined with a reference technique such as LC-MS/MS; vitamin D receptor expression was evaluated; and 25(OH)D 2 and 25(OH)D 3 forms were analyzed, allowing us to extract valuable and additional information less explored in athletes [37].

Conclusions
Players of indoor team sports such as handball would be a population at risk of vitamin D deficiency.The 16-weeks competition improved VDR gene expression and body composition.The associations observed between VDR gene expression, and the variables of the study evidenced the

Table 2
Comparative and evolutionary analysis of body composition, vitamin D and its related minerals for non-athletes and handball athletes BMI: Body Mass Index; ES: Cohen's d effect size; 25(OH)D: 25-hydroxyvitamin D; 25(OH)D 2 : 25-hydroxyvitamin D 2 ; 25(OH)D 3 : 25-hydroxyvitamin D 3 .The quantitative variables data were expressed as the mean and standard deviation (SD).Independent t-test analysis was applied for comparing non-athletes vs. handball athletes.Paired t-test analysis was used to compare values for non-athletes and handball athletes at 16-weeks follow-up (p < 0.05 being considered statistically significant).The effect sizes (ES) were calculated using the Cohen's comparison of groups of equal size.References values: Calcium, 8.6-10.2mg/dL; Phosphorous, 2.5-4.5 mg/dL; Magnesium, 1.7-2.60 mg/dL

Table 3
Matrix for correlations coefficients (r) showing the simple linear relationship among 25(OH)D, 25(OH)D 2 , 25(OH)D 3 levels and VDR expression with body composition and vitamin D related minerals in non-athletes and handball athletes