Next Article in Journal
Study on Water Suitability of Apple Plantations in the Loess Plateau under Climate Change
Next Article in Special Issue
Controversies Surrounding Vitamin D: Focus on Supplementation and Cancer
Previous Article in Journal
The Impact of Psychological Distress on Incident Functional Disability in Elderly Japanese: The Ohsaki Cohort 2006 Study
Previous Article in Special Issue
Knowledge, Attitudes and Perceptions towards Vitamin D in a UK Adult Population: A Cross-Sectional Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vitamin D Status and Analysis of Specific Correlates in Preschool Children: A Cross-Sectional Study in Southern Croatia

1
Teaching Institute of Public Health of Split Dalmatian County, 21000 Split, Croatia
2
Faculty of Kinesiology, University of Split, 21000 Split, Croatia
3
Department of Medical Laboratory Diagnostics, University Hospital of Split, 21000 Split, Croatia
4
School of Medicine, University of Mostar, 88000 Mostar, Bosnia and Herzegovina
5
Department of Pediatrics, University Hospital of Split, 21000 Split, Croatia
6
School of Medicine, University of Split, 21000 Split, Croatia
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(11), 2503; https://doi.org/10.3390/ijerph15112503
Submission received: 30 September 2018 / Revised: 2 November 2018 / Accepted: 6 November 2018 / Published: 8 November 2018
(This article belongs to the Special Issue Vitamin D and Public Health)

Abstract

:
Vitamin D deficiency is a globally important problem, particularly in children, but there is a lack of information regarding this deficiency in preschool children from southeastern Europe. This study aimed to establish the levels of 25-hydroxyvitamin D (25(OH)D) and associations of gender, time spent outdoors, physical activity (PA), and body mass index (as predictors) with the 25(OH)D level (outcome) in healthy preschool children. The participants were preschoolers (all 5–6 years of age) from southern Croatia. All the participants were tested during their mandatory medical examination 6–7 months prior to school enrollment. The PA was obtained using the preschool-age physical activity questionnaire (Pre-PAQ), which categorizes PA into five levels (from sedentary to vigorous PA). The prevalence of 25(OH)D deficiency was high: 58% of the children had 25(OH)D levels of <50 nmol/L (deficiency), and an additional 29% had an insufficient level of 25(OH)D (50–75 nmol/L). Boys had higher levels of 25(OH)D than girls. A multinomial regression using 25(OH)D categories as the outcome and a sufficient level (>75 nmol/L) as the reference value identified gender as the only significant predictor of 25(OH)D status, with boys being at lower risk for 25(OH)D deficiency than girls. These results showed a high prevalence of 25(OH)D deficiency in preschoolers from the southern part of Croatia, which is additionally alarming based on the geographical position of the studied region (42° N) and its high number of sunshine hours (>2600 h per year). Future studies examining other potential correlates of 25(OH)D in the region are warranted.

1. Introduction

Vitamin D is mostly known as a crucial factor in bone health but is also related to hypertension, metabolic syndrome, autoimmune and infectious diseases, and even cancer [1,2,3]. The vitamin D status is based on the serum level of 25-hydroxyvitamin D (25(OH)D), and studies have provided evidence showing that 25(OH)D deficiency differs across regions of the world [4,5,6,7]. Deficiency of 25(OH)D is particularly important in children [8]. More precisely, although it was originally hypothesized that children were not at risk of vitamin D deficiency due to their dietary intake and outdoor activities, more recent studies have identified a high prevalence of 25(OH)D insufficiency in children and adolescents worldwide [6,9,10,11]. Moreover, studies have frequently demonstrated an inverse association between obesity indicators (i.e., body mass index and body fat measures) and 25(OH)D status, with a higher prevalence of 25(OH)D deficiency among obese/overweight children [8,12,13]. In general, it is clear that failure to address a child’s low 25(OH)D status could have serious long-term negative health implications [14]. As a result, an evident increase in the body of knowledge on factors that might be associated with 25(OH)D status in children has been observed in recent years [15,16,17].
Physical activity (PA) is any movement of the body produced by skeletal muscles that results in a greater energy expenditure than that observed at rest levels [18]. The benefits of increased levels of PA are numerous and well documented and include various physiological benefits (i.e., improved cardiovascular health status, lowered blood pressure, and better hormonal function) and psychological improvements (i.e., reduced mental stress and anxiety, improved self-esteem, and socialization) [19,20]. Additionally, the level of PA is reportedly positively correlated with 25(OH)D status. In a study in which the PA level was self-reported, the researchers obtained a positive correlation between PA and 25(OH)D status in U.S. adults (>20 years) [21], and similar conclusions were reached in other studies [22,23,24]. Although outdoor PA is particularly beneficial, even indoor PA has been recognized as beneficial with regard to 25(OH)D status [21,23].
In children, PA is even more important because it positively influences cognitive and emotional functions [25], and appropriate PA in childhood has the potential to reduce the risk of obesity [26]. Recent years have seen an increase in public interest in the relationship between PA and 25(OH)D status in children, mostly because studies completed over the last decade have regularly confirmed a high prevalence of 25(OH)D deficiency in children and adolescents [9,27]. However, the results of studies that specifically investigated a possible correlation between PA and 25(OH)D status in children are not consistent [16,17,28,29].
In brief, in a German study, PA was not correlated with 25(OH)D status in 11- to 13-year-old children [29], which is consistent with the results of a Danish study that observed eight- to 11-year-old children [28]. In contrast, another Danish study reported a positive association between PA and 25(OH)D status in younger children (aged 4–8 years) [16], and similar correlations and potentially positive effects of increased PA on 25(OH)D status were observed in children and adolescents from other world regions [9,17,30,31]. Finally, in a very recent study, Brazilian researchers reported a gender-specific association between PA and 25(OH)D status in 12- to 17-year-old adolescents; specifically, these researchers showed that PA is significantly correlated with 25(OH)D status only in boys [15]. Collectively, it appears that the association between PA and 25(OH)D status in children is not universal and should be identified specifically by (i) gender (i.e., some studies reported evidence of gender-specific associations) [15], (ii) narrow age groups (due to the specific influence of different developmental factors on 25(OH)D status in different age groups) [16,28], and (iii) geographical regions/countries (mainly due to differences in sunshine hours among different geographical locations) [9,18,27].
Croatia is located in the southeastern part of Europe, and during the last decade, the 25(OH)D status has received increasing attention in studies that investigated various populations in Croatia [32,33]. However, these studies mostly examined adolescent or adult populations or very specific samples of children (i.e., children at risk for type 1 diabetes) [34,35], and the 25(OH)D status and its correlates in apparently healthy Croatian children have not been investigated. Indeed, even in the very recent official guidelines, national public health authorities noted a lack of empirical evidence that could provide specific insights into the prevalence and correlates of 25(OH)D status in Croatia [36]. Moreover, few studies have investigated this problem in the region where Croatia is located (i.e., territory of the former Yugoslavia) [37,38,39]. Finally, to the best of our knowledge, no investigation on this territory has examined either the 25(OH)D status specifically in preschool children or the possible associations between PA and 25(OH)D status in this highly vulnerable population.
Thus, this study aimed to investigate the prevalence of 25(OH)D deficiency in healthy preschool children from the southern part of Croatia and the gender-specific associations of time spent outdoors, PA levels, and body mass index (BMI; an indicator of obesity) with 25(OH)D level in this population. Initially, we hypothesized that the PA levels and time spent outdoors would be positively correlated with 25(OH)D status, whereas the BMI would be negatively correlated with 25(OH)D status. We also hypothesized that these associations would be irrespective of gender. Due to the lack of studies on these problems in Croatia and surrounding countries, the results should help raise public awareness of the problem and contribute to a better understanding of the possible influence of the studied correlates on 25(OH)D status in this age group.

2. Materials and Methods

2.1. Participants

The participants in this cross-sectional study were healthy preschool-age children from the southern part of Croatia (age range: 5–6 years; mean age: 6.0 ± 0.4 years; n = 260; 128 females). All the children were tested during their regular and mandatory medical examination prior to school enrollment. Prior to the study, the investigation and protocol were approved by the Ethical Board of University of Split, School of Medicine (EBO: 2181-198-03-04-16-0009), and parent(s) were informed of the study purposes and of the potential risks and benefits of participation. The testing was performed during March and April 2017 after written parental consent was obtained. Although this study included an analysis of PA (see Section 2.2 for details), which is not a regular part of the medical examination, the response rate was 84%. In this study, we included only those children who were identified as healthy with regard to their participation in the regular scholastic system in Croatia (participation in the first year of elementary school starting the following September, approximately six months after testing). Based on (i) an estimated 70% prevalence of 25(OH)D deficiency/insufficiency [37], (ii) a theoretical population sample of approximately 3000 preschool children in the studied region, and (iii) a significance level of p < 0.05, the sample size required for this investigation was calculated to be 250 participants.

2.2. Variables

The variables investigated in this study included the child’s gender (male or female, as indicated by the accompanying parent), body mass index, indicators of PA status, and time spent outdoors, which were treated as independent variables, as well as the 25(OH)D status, which was treated as the dependent variable.
After a 6–8 h fasting period, a 5 mL venous blood sample was collected in a BD Vacutainer® SSTII Advance (BD, Plymouth, UK) container. The sample was centrifuged and stored at −20 °C until its analysis at the accredited laboratory of the University Hospital of Split (Split, Croatia). The 25(OH)D levels were measured using a commercially available Elecsys® Vitamin D total assay with a Cobas e601 analyzer (Roche Diagnostics International Ltd., Rotkreuz, Switzerland), which uses a competitive electrochemiluminescence binding technique. This assay employs a vitamin D-binding protein as the capture protein to bind vitamin D3(25-OH) and vitamin D2(25-OH). The detection range of the test is 7.5–175 nmol/L 25(OH)D, and the sensitivity of the assay is 5 nmol 25(OH)D/L. The intraclass CV was 5.3% at 39 nmol 25(OH)D/L, 5.6% at 67.1 nmol 25(OH)D/L, 6.7% at 165 nmol 25(OH)D/L, 2.2% at 174 nmol 25(OH)D/L, 3.9% at 70.8 nmol 25(OH)D/L, and 5.2% at 39.5 nmol 25(OH)D/L. In addition, the Limit of Blank, Limit of Detection, and Limit of Quantification were 5.0 nmol/L, 7.5 nmol/L and 12.5 nmol/L, respectively. All the samples were analyzed in duplicate. The method was standardized based on international standards [40,41].
The 25(OH)D measurement results are presented as raw values (nmol/L). Moreover, for the interpretation and identification of 25(OH)D levels, the raw data were further categorized into three categories (e.g., vitamin D sufficiency, insufficiency and deficiency). Specifically, although there is no absolute consensus regarding 25(OH)D status, we used the regularly used values of ≥75 nmol/L, 51–75 nmol/L, and <50 nmol/L for the classification of vitamin D sufficiency, insufficiency and deficiency, respectively [42]. Such categorization was also used for some of the statistical analyses (see details in Section 2.3). Additionally, a value of <25 mmol/L was considered to indicate severe vitamin D deficiency [43].
The BMI was calculated based on the measured body height (BH) and body mass (BM) using Equation (1) [44].
BMI = BM   ( kg ) BH   ( m ) 2
The BH and BM were measured using a standardized stadiometer and scale, respectively, and these measurements were performed by experienced technicians. To avoid diurnal variation, all the children were tested at the same time of the day (between 7 and 9 a.m.). The BH was measured to the nearest 1 cm, and BM was reported to the nearest 0.1 kg.
The time spent outdoors was reported by the parents as the number of hours their child spent outdoors during the last week (including the weekend).
The indicators of PA included variables obtained from the preschool-age physical activity questionnaire (Pre-PAQ), which is a reliable and valid measuring tool that aim to evaluate different activity levels of preschool children [45,46,47]. In general, the Pre-PAQ categorizes activity into five progressive levels (L1: stationary with no movement, L2: stationary with limb or trunk movement, L3: slow activity, L4: medium activity, and L5: fast-paced activity). In addition to analyzing all five levels, L1 and L2 (stationary activities) were also combined (L1 + L2) in this study, as previously suggested [45]. The Pre-PAQ (in the Croatian language) was completed by the parent who accompanied the child to the medical examination. The questionnaire was translated to the Croatian language by a native speaker of both English and Croatian (the Pre-PAQ was originally developed in English by Australian authors), and the translation was then checked by three academicians who participated in this study (two medical doctors and one exercise scientist). Although the reliability of the Pre-PAQ was previously proven to be appropriate (as detailed in previous studies [45,48]), for the purpose of this study, the test-retest reliability of the Pre-PAQ was checked by asking 25 parents to complete the questionnaire twice within a time frame of 15 days. The calculated Pearson’s correlation coefficients between the results ranged from 0.73 (for L4) to 0.81 (for L5), indicating that the appropriate to high reliability of the questionnaire.

2.3. Statistics

First, we calculated descriptive statistics for all the variables. Specifically, for parametric variables, we calculated the means and standard deviations, whereas frequencies (counts) and percentages are reported for nonparametric variables. The differences between boys and girls were then analyzed using t-test for independent samples (for parametric variables) and χ2 test (for categorical variables). The associations among the variables were established through Pearson’s product moment correlation coefficients, and the associations between the independent variables and 25(OH)D levels (outcome) were analyzed by multinomial logistic regression. For this statistical analysis, the outcome comprised three results/categories, namely, vitamin D sufficiency, vitamin D insufficiency, and vitamin D deficiency, and vitamin D sufficiency (>75 nmol/L) was used as the reference value. Based on the findings of previous studies that identified gender as a significant covariate of 25(OH)D [15], the analyses were performed for the total sample and then stratified for gender. Statistica for Windows version 13 (Dell Inc., Tulsa, OK, USA) was used for all the calculations, and a p value of <0.05 was considered significant.

3. Results

The descriptive statistics for the measured variables for the total population and the two genders are presented in Table 1. Boys had significantly higher levels of vitamin D (t-value: 2.52, p < 0.01) and higher numbers of PA scores of L2 (t-value: 2.34: p < 0.05) and L5 (t-value: 2.25, p < 0.05) compared with girls.
In the total sample (e.g., not stratified by gender), severe 25(OH)D deficiency was evident in 12% of the participants, and an additional 46% of the participants presented with 25(OH)D levels of 25–50 mmol/L (Figure 1). The prevalence of vitamin D deficiency was higher in girls: 17% of the boys and 8% of the girls had sufficient levels of 25(OH)D (χ2: 9.97, p < 0.02).
The Pearson’s product moment correlation coefficients between the independent variables and 25(OH)D levels (outcome) did not reach statistical significance in either the total sample or the gender-based subsamples (Table 2).
The multinomial regression analysis performed using 25(OH)D status as the outcome (Table 3) confirmed the previously reported χ2 results. In brief, in the total sample of participants, gender was the only significant predictor of the 25(OH)D status, and boys were at lower risk for vitamin D deficiency (odds ratio (OR): 0.33; 95% confidence interval (CI): 0.15–0.74) than girls. Moreover, 25(OH)D level was not influenced by the PA levels, BMI, or time spent outdoors, and this finding was obtained for the total population of participants and for boys (Table 4) and girls (Table 5).

4. Discussion

Our results showed a high prevalence of 25(OH)D deficiency in the studied population of preschool children. Additionally, boys were more physically active and had higher levels of 25(OH)D than girls. Moreover, the PA levels, BMI, and time spent outdoors (predictors) were not significantly correlated with the 25(OH)D level in the studied sample of preschool children. Therefore, our initial study hypotheses on the significant associations between predictors and 25(OH)D were rejected.

4.1. 25(OH)D Level in Croatian Preschool Children

The reported results showed alarming figures for the 25(OH)D status in children from the southern part of Croatia. 25(OH)D deficiency was evident among 58% of the participants, and a higher prevalence of 25(OH)D deficiency was observed in girls (49% and 67% for boys and girls, respectively). Although the high prevalence of 25(OH)D deficiency is disturbing, these results are in accordance with those of previous studies performed in similar samples of participants around the globe. Specifically, a very recent Iranian study reported subnormal 25(OH)D levels in 69% of children aged 30–72 months [11], which is similar to the prevalence reported by U.S. authors for a representative sample of one- to 11-year-old children in the period 2001–2016 [43]. However, a Danish study that investigated four- to eight-year-old children reported a lower prevalence of 25(OH)D deficiency (e.g., approximately 30% children with a 25(OH)D level of <50 nmol/L) [16]. This difference between the results presented here and the results obtained in the Danish study are particularly intriguing due to the latitude difference between Croatia and Denmark (the latitudes of these countries are approximately 42° N and 55° N, respectively). Due to its northern latitude and thus fewer hours of sunshine, it would be expected that Denmark would have a higher prevalence of 25(OH)D deficiency. However, it appears that other factors (i.e., time spent outdoors, sun-protective behavior, and dietary intake) are more significant influencing factors on 25(OH)D level than latitude and the consequent number of sunshine hours [49].
To the best of our knowledge, this study constitutes the first investigation of the 25(OH)D levels specifically in preschool children in Croatia and the surrounding countries (the territory of the former Yugoslavia). Additionally, studies performed on other samples in Croatia (i.e., adults and adolescents) mostly focused on specific populations with severe health conditions, which makes comparison of the results additionally difficult [50,51]. Therefore, we could not specifically compare our results with those previously reported for children, and future studies are needed to confirm our results. However, it appears that 25(OH)D deficiency should be recognized as a problem in this territory. Specifically, although we could not find representative data for Croatia, our colleagues from Bosnia and Herzegovina recently reported alarming figures for 25(OH)D deficiency in a representative sample from their country [37]. In short, although these studies observed a sufficient 25(OH)D level in only 20% of the tested participants, only 34% of the participants younger than 18 years had a sufficient serum level of 25(OH)D [37]. It must be stated that this cited study included participants across all age categories, and therefore, the results are not absolutely comparable to ours because we exclusively investigated preschool children. However, because Croatia, Bosnia and Herzegovina are located within the same territory and share a similar cultural and social background, which results in similar life habits (i.e., time spent outdoors and sun-protection, nutrition), the results are highly indicative and note the necessity of a more profound investigation of the problem of 25(OH)D deficiency in the whole territory.

4.2. Gender Differences in the 25(OH)D and Physical Activity Levels

Our study found that boys had higher levels of 25(OH)D than girls. In general, gender differences in the 25(OH)D status of children are not frequently reported, and most previous studies have shown no significant differences in 25(OH)D levels between boys and girls. For example, investigations of Asian children until early adolescence showed no significant differences in 25(OH)D levels between boys and girls [52,53,54,55], which is consistent with the results of Icelandic and Finnish studies of six- to eight-year-old children [10,56]. Moreover, the previously cited Danish study reported female gender as one of the risk factors for 25(OH)D deficiency in eight- to 11-year-old children [28]. However, it must be stated that in the previously cited studies in which the authors did not confirm gender differences, the male participants exhibited somewhat higher levels of 25(OH)D than the female participants, but these differences did not reach statistical significance [10,52,53,54,55,56].
The differences in PA levels between genders might help explain the relative inconsistency between our results and those of previous studies in which the authors reported no influence of gender on 25(OH)D status. In short, in our study, boys had higher PA levels of L2 and L5 than girls, which could also influence the 25(OH)D levels in the total population of children. Although the present study is the first to use the Pre-PAQ for measuring the PA of children and one of the first to report the PA of preschool children in the region, our results are in accordance with those of previous studies that investigated similarly aged children in other countries, although in some of them, the gender differences appeared to be significant at a somewhat older age (e.g., 10–12 years) [57,58]. Although the interpretation of gender differences in PA is beyond the scope of this paper, it is possible that such differences influence even the differences in 25(OH)D levels between genders discussed previously.

4.3. Lack of Association of Body Mass Index, Physical Activity and Time Spent Outdoors with 25(OH)D Levels

Although sometimes criticized as not being an objective measure of body composition (i.e., both muscle mass and adiposity can increase the BMI numerical value), the BMI is still widely used as an indicator of being overweight and an indirect indicator adiposity, particularly during childhood [59]. The association between adiposity (e.g., high BMI values) and 25(OH)D levels has been frequently confirmed. Specifically, studies have indicated an inverse association between 25(OH)D concentrations and indicators of obesity in children [60,61,62], and this association is generally explained by metabolic pathways (i.e., vitamin D affects lipolysis and adipogenesis in human adipocytes through its role in regulating the intracellular calcium concentrations) [63]. However, in our study, the 25(OH)D level was not correlated with BMI in either preschool boys or preschool girls.
To explain this inconsistency between our findings (e.g., no significant association between BMI and 25(OH)D) and the results of previous studies that found negative correlations between adiposity and 25(OH)D levels, two specific issues should be highlighted. First, our study sample of preschool children had a relatively low proportion of overweight and obese children (i.e., only 1% of the participants had values of >25 kg/m2, which is used as the threshold to indicate being overweight/obese), and this relatively low BMI variability likely decreased the influence of the BMI on 25(OH)D levels. In contrast, previous studies examined samples with higher prevalence of overweight/obesity, and consequently, the association between indicators of adiposity and 25(OH)D levels was more easily observed [60,61,62,64]. Second, prior studies rarely examined preschool children, and in most investigations where preschoolers were included, they formed part of the larger sample (i.e., preschool and school-age children) [64]. In addition, in previous investigations of only preschoolers, indicators of obesity were not related to 25(OH)D levels [65]. Therefore, although the relationship between indicators of adiposity and 25(OH)D status has been frequently confirmed in somewhat older children and adolescents, whether it truly exists in this age group remains questionable [60,61,62,64].
Although some previous studies reported potential beneficial effects of physical activity on 25(OH)D levels in children and adolescents [9,15], our results did not confirm such an association among preschoolers from southern Croatia. To explain our findings, several issues deserve attention. The first explanation for the lack of an association between PA and 25(OH)D levels is specific to the participant sample. Briefly, this study enrolled children in a very narrow age range (5–6 years), and it is possible that some other factors are more important determinants of 25(OH)D levels in such participant samples. Second, the PA level was determined using the Pre-PAQ. As explained in the Methods section, the parents answered the questionnaire for their children, and their participation was voluntary. It is thus possible that the parents did not objectively report their child’s PA either because (i) they were personally unaware of their child’s PA-related behavior or (ii) due to a tendency to provide ocially desirable responses (i.e., because physical activity is “socially desirable behavior”, parents will report a higher level of vigorous physical activity for their child) [66].
The “social desirability bias” is also a probable explanation for the lack of an association between sun-exposure time and 25(OH)D levels. Specifically, until recently, the problem of 25(OH)D deficiency was not particularly discussed in the Croatian scientific and professional literature because of Croatia’s geographical location (42° N) and its high number of sunshine hours (>2600 h per year). However, recent epidemiological data have provided evidence indicating low 25(OH)D levels even in southern countries [67]. Consequently, the importance of sun-exposure time in the prevention of 25(OH)D deficiency has become a widely known issue in our country [36], and an increased level of sun exposure has thus become a certain type of “socially desirable behavior”, similar to physical activity. This trend is particularly apparent in parents with regard to their children’s sun-exposure time. Consequently, we may logically question the objectivity of the parents’ reporting of their children’s sun-exposure time. Additionally, we cannot ignore the fact that time spent outdoors was not reported at the same moment for all participants (i.e., children were tested during March), which might result in differences due to weather conditions between the weeks when children were tested.

4.4. Limitations and Strengths

In addition to the study limitations discussed previously (i.e., social-desirability bias and parental unawareness of the children’s behavior), this study has some additional limitations. First, in this study, we used the children’s BMI as the indicator of anthropometric status, and future studies should observe the association between body composition and 25(OH)D levels using more precise measures of adiposity and body build. Next, knowing the possible influence of dietary intake on 25(OH)D levels, the lack of data on dietary intake is a significant limitation of the study. Additionally, this study observed children from only one region (the Mediterranean part of Croatia). Because of the differences in other factors that might influence 25(OH)D levels in preschool children (i.e., type of diet and 25(OH)D supplementation), the generalization of these results is relatively limited.
This study is one of the first to investigate 25(OH)D levels and some of its correlates specifically in healthy preschool children from the region and is likely the first to address this issue in Croatia. The 25(OH)D levels of all participants were determined within a relatively short time span, and the samples were analyzed in only one (accredited) laboratory. The sample of participants included children in a very narrow age range (5–6 years), and the results may be considered plausible for the studied sample. Finally, the gender-specific approach is one of the important strengths of this investigation. Therefore, despite the previously mentioned limitations, the results and observations presented could be used in future, more profound investigations of this problem.

5. Conclusions

Our results showed that 25(OH)D deficiency is very prevalent in apparently healthy preschool children aged 5–6 years from the southern part of Croatia. Specifically, 25(OH)D deficiency was evident in almost 60% of children (49% of males and 67% of females). Gender was a strong factor influencing 25(OH)D levels, with boys at lower risk for 25(OH)D deficiency than girls. Although the results showed that girls were also less physically active than boys, we could not confirm that PA level is a factor that significantly influences the 25(OH)D levels in the observed children. The latter conclusion could be additionally supported by the fact that the BMI, as an indicator of adiposity, was also not related to the 25(OH)D level. Finally, we did not find an association between time spent outdoors and the 25(OH)D level.
Because this study observed children from one specific geographical region (i.e., the Mediterranean part of Croatia), which is characterized by a mild climate and a high number of sunshine hours, the results should be confirmed in other regions. Additionally, in further studies, special attention should be paid to other potential correlates of 25(OH)D levels, such as sun-protective behavior (i.e., usage of sun-protective creams, lotions, and clothing), type of diet, and 25(OH)D supplementation.

Author Contributions

Data curation, Z.K., D.S.D. and K.E.; Formal analysis, B.G. and D.S.D.; Funding acquisition, O.U., J.M.; Investigation, Z.K., J.M. and K.E.; Methodology, B.G., N.Z. and J.M.; Project administration, Z.K. and J.M.; Resources, Z.S. and O.U.; Writing—original draft, Z.S., N.Z., O.U. and M.P.

Funding

The APC was funded by the University of Split, Split, Croatia.

Acknowledgments

The authors are particularly grateful to all parents and children who voluntarily participated in the study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dobnig, H. A review of the health consequences of the vitamin D deficiency pandemic. J. Neurol. Sci. 2011, 311, 15–18. [Google Scholar] [CrossRef] [PubMed]
  2. Prentice, A.; Schoenmakers, I.; Jones, K.S.; Jarjou, L.M.; Goldberg, G.R. Vitamin D Deficiency and Its Health Consequences in Africa. Clin. Rev. Bone Miner. Metab. 2009, 7, 94–106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Norval, M.; Coussens, A.K.; Wilkinson, R.J.; Bornman, L.; Lucas, R.M.; Wright, C.Y. Vitamin D Status and Its Consequences for Health in South Africa. Int. J. Environ. Res. Public Health 2016, 13, 1019. [Google Scholar] [CrossRef] [PubMed]
  4. Kheiri, B.; Abdalla, A.; Osman, M.; Ahmed, S.; Hassan, M.; Bachuwa, G. Vitamin D deficiency and risk of cardiovascular diseases: A narrative review. Clin. Hypertens. 2018, 24, 9. [Google Scholar] [CrossRef] [PubMed]
  5. Mpandzou, G.; Ait Ben Haddou, E.; Regragui, W.; Benomar, A.; Yahyaoui, M. Vitamin D deficiency and its role in neurological conditions: A review. Rev. Neurol. 2016, 172, 109–122. [Google Scholar] [CrossRef] [PubMed]
  6. Misra, M.; Pacaud, D.; Petryk, A.; Collett-Solberg, P.F.; Kappy, M.; Drug and Therapeutics Committee of the Lawson Wilkins Pediatric Endocrine Society. Vitamin D deficiency in children and its management: Review of current knowledge and recommendations. Pediatrics 2008, 122, 398–417. [Google Scholar] [CrossRef] [PubMed]
  7. Yeum, K.J.; Song, B.C.; Joo, N.S. Impact of Geographic Location on Vitamin D Status and Bone Mineral Density. Int. J. Environ. Res. Public Health 2016, 13, 184. [Google Scholar] [CrossRef] [PubMed]
  8. Alaklabi, A.M.; Alsharairi, N.A. Current Evidence on Vitamin D Deficiency and Metabolic Syndrome in Obese Children: What Does the Evidence from Saudi Arabia Tell Us? Children 2018, 5, 11. [Google Scholar] [CrossRef] [PubMed]
  9. Al-Othman, A.; Al-Musharaf, S.; Al-Daghri, N.M.; Krishnaswamy, S.; Yusuf, D.S.; Alkharfy, K.M.; Al-Saleh, Y.; Al-Attas, O.S.; Alokail, M.S.; Moharram, O.; et al. Effect of physical activity and sun exposure on vitamin D status of Saudi children and adolescents. BMC Pediatr. 2012, 12, 92. [Google Scholar] [CrossRef] [PubMed]
  10. Bjarnadottir, A.; Kristjansdottir, A.G.; Hrafnkelsson, H.; Johannsson, E.; Magnusson, K.T.; Thorsdottir, I. Insufficient autumn vitamin D intake and low vitamin D status in 7-year-old Icelandic children. Public Health Nutr. 2015, 18, 208–217. [Google Scholar] [CrossRef] [PubMed]
  11. Esmaeili Dooki, M.R.; Moslemi, L.; Moghadamnia, A.A.; Alijanpour Aghamaleki, M.; Bijani, A.; Pornasrollah, M.; Ashrafianamiri, H. Vitamin D status in preschool children: Should vitamin D supplementation, preventing vitamin D deficiency be continued in children over 2 years? J. Public Health 2018, fdy147. [Google Scholar] [CrossRef] [PubMed]
  12. Flores, A.; Flores, M.; Macias, N.; Hernández-Barrera, L.; Rivera, M.; Contreras, A.; Villalpando, S. Vitamin D deficiency is common and is associated with overweight in Mexican children aged 1–11 years. Public Health Nutr. 2017, 20, 1807–1815. [Google Scholar] [CrossRef] [PubMed]
  13. Wakayo, T.; Whiting, S.J.; Belachew, T. Vitamin D Deficiency is Associated with Overweight and/or Obesity among Schoolchildren in Central Ethiopia: A Cross-Sectional Study. Nutrients 2016, 8, 190. [Google Scholar] [CrossRef] [PubMed]
  14. Zadka, K.; Palkowska-Gozdzik, E.; Rosolowska-Huszcz, D. The State of Knowledge about Nutrition Sources of Vitamin D, Its Role in the Human Body, and Necessity of Supplementation among Parents in Central Poland. Int. J. Environ. Res. Public Health 2018, 15, 1489. [Google Scholar] [CrossRef] [PubMed]
  15. Da Silva, A.C.M.; Cureau, F.V.; de Oliveira, C.L.; Giannini, D.T.; Bloch, K.V.; Kuschnir, M.C.C.; Dutra, E.S.; Schaan, B.D.; de Carvalho, K.M.B. Physical activity but not sedentary time is associated with vitamin D status in adolescents: Study of cardiovascular risk in adolescents (ERICA). Eur. J. Clin. Nutr. 2018. [Google Scholar] [CrossRef] [PubMed]
  16. Mortensen, C.; Molgaard, C.; Hauger, H.; Kristensen, M.; Damsgaard, C.T. Sun behaviour and physical activity associated with autumn vitamin D status in 4–8-year-old Danish children. Public Health Nutr. 2018, 1–10. [Google Scholar] [CrossRef] [PubMed]
  17. Saki, F.; Dabbaghmanesh, M.H.; Omrani, G.R.; Bakhshayeshkaram, M. Vitamin D deficiency and its associated risk factors in children and adolescents in southern Iran. Public Health Nutr. 2017, 20, 1851–1856. [Google Scholar] [CrossRef] [PubMed]
  18. Fernandes, M.R.; Barreto, W.D.R.J. Association between physical activity and vitamin D: A narrative literature review. Rev. Assoc. Med. Bras. 2017, 63, 550–556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Brown, H.E.; Pearson, N.; Braithwaite, R.E.; Brown, W.J.; Biddle, S.J. Physical activity interventions and depression in children and adolescents: A systematic review and meta-analysis. Sports Med. 2013, 43, 195–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Warburton, D.E.R.; Bredin, S.S.D. Health benefits of physical activity: A systematic review of current systematic reviews. Curr. Opin. Cardiol. 2017, 32, 541–556. [Google Scholar] [CrossRef] [PubMed]
  21. Scragg, R.K.; Camargo, C.A., Jr.; Simpson, R.U. Relation of serum 25-hydroxyvitamin D to heart rate and cardiac work (from the National Health and Nutrition Examination Surveys). Am. J. Cardiol. 2010, 105, 122–128. [Google Scholar] [CrossRef] [PubMed]
  22. Millen, A.E.; Wactawski-Wende, J.; Pettinger, M.; Melamed, M.L.; Tylavsky, F.A.; Liu, S.; Robbins, J.; LaCroix, A.Z.; LeBoff, M.S.; Jackson, R.D. Predictors of serum 25-hydroxyvitamin D concentrations among postmenopausal women: The Women’s Health Initiative Calcium plus Vitamin D clinical trial. Am. J. Clin. Nutr. 2010, 91, 1324–1335. [Google Scholar] [CrossRef] [PubMed]
  23. Kluczynski, M.A.; Lamonte, M.J.; Mares, J.A.; Wactawski-Wende, J.; Smith, A.W.; Engelman, C.D.; Andrews, C.A.; Snetselaar, L.G.; Sarto, G.E.; Millen, A.E. Duration of physical activity and serum 25-hydroxyvitamin D status of postmenopausal women. Ann. Epidemiol. 2011, 21, 440–449. [Google Scholar] [CrossRef] [PubMed]
  24. Tseng, M.; Giri, V.; Bruner, D.W.; Giovannucci, E. Prevalence and correlates of vitamin D status in African American men. BMC Public Health 2009, 9, 191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bidzan-Bluma, I.; Lipowska, M. Physical Activity and Cognitive Functioning of Children: A Systematic Review. Int. J. Environ. Res. Public Health 2018, 15, 800. [Google Scholar] [CrossRef] [PubMed]
  26. Czyz, S.H.; Toriola, A.L.; Starosciak, W.; Lewandowski, M.; Paul, Y.; Oyeyemi, A.L. Physical Fitness, Physical Activity, Sedentary Behavior, or Diet—What Are the Correlates of Obesity in Polish School Children? Int. J. Environ. Res. Public Health 2017, 14, 664. [Google Scholar] [CrossRef] [PubMed]
  27. Maalouf, J.; Nabulsi, M.; Vieth, R.; Kimball, S.; El-Rassi, R.; Mahfoud, Z.; El-Hajj Fuleihan, G. Short- and long-term safety of weekly high-dose vitamin D3 supplementation in school children. J. Clin. Endocrinol. Metab. 2008, 93, 2693–2701. [Google Scholar] [CrossRef] [PubMed]
  28. Petersen, R.A.; Dalskov, S.M.; Sorensen, L.B.; Hjorth, M.F.; Andersen, R.; Tetens, I.; Krarup, H.; Ritz, C.; Astrup, A.; Michaelsen, K.F.; et al. Vitamin D status is associated with cardiometabolic markers in 8–11-year-old children, independently of body fat and physical activity. Br. J. Nutr. 2015, 114, 1647–1655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Lammle, L.; Bergmann, K.; Bos, K.; Koletzko, B. Predictors of differences in vitamin D levels in children and adolescents and their relation to endurance performance. Ann. Nutr. Metab. 2013, 62, 55–62. [Google Scholar] [CrossRef] [PubMed]
  30. Absoud, M.; Cummins, C.; Lim, M.J.; Wassmer, E.; Shaw, N. Prevalence and predictors of vitamin D insufficiency in children: A Great Britain population based study. PLoS ONE 2011, 6, e22179. [Google Scholar] [CrossRef] [PubMed]
  31. Dong, Y.; Pollock, N.; Stallmann-Jorgensen, I.S.; Gutin, B.; Lan, L.; Chen, T.C.; Keeton, D.; Petty, K.; Holick, M.F.; Zhu, H. Low 25-hydroxyvitamin D levels in adolescents: Race, season, adiposity, physical activity, and fitness. Pediatrics 2010, 125, 1104–1111. [Google Scholar] [CrossRef] [PubMed]
  32. Grazio, S.; Naglic, D.B.; Anic, B.; Grubisic, F.; Bobek, D.; Bakula, M.; Kavanagh, H.S.; Kuna, A.T.; Cvijetic, S. Vitamin D serum level, disease activity and functional ability in different rheumatic patients. Am. J. Med. Sci. 2015, 349, 46–49. [Google Scholar] [CrossRef] [PubMed]
  33. Laktasic-Zerjavic, N.; Korsic, M.; Crncevic-Orlic, Z.; Kovac, Z.; Polasek, O.; Soldo-Juresa, D. Vitamin D status, dependence on age, and seasonal variations in the concentration of vitamin D in Croatian postmenopausal women initially screened for osteoporosis. Clin. Rheumatol. 2010, 29, 861–867. [Google Scholar] [CrossRef] [PubMed]
  34. Boraska, V.; Skrabic, V.; Zeggini, E.; Groves, C.J.; Buljubasic, M.; Peruzovic, M.; Zemunik, T. Family-based analysis of vitamin D receptor gene polymorphisms and type 1 diabetes in the population of South Croatia. J. Hum. Genet. 2008, 53, 210–214. [Google Scholar] [CrossRef] [PubMed]
  35. Zemunik, T.; Skrabic, V.; Boraska, V.; Diklic, D.; Terzic, I.M.; Capkun, V.; Peruzovic, M.; Terzic, J. FokI polymorphism, vitamin D receptor, and interleukin-1 receptor haplotypes are associated with type 1 diabetes in the Dalmatian population. J. Mol. Diagn. 2005, 7, 600–604. [Google Scholar] [CrossRef]
  36. Vranesic Bender, D.; Giljevic, Z.; Kusec, V.; Laktasic Zerjavic, N.; Bosnjak Pasic, M.; Vrdoljak, E.; Lkinas Kelecic, D.; Reiner, Z.; Anic, B.; Krznaric, Z. Guidelines for the Prevention, Detection and Therapy of Vitamin D Deficiency in Adults. Lijec. Vjesn. 2016, 138, 121–132. [Google Scholar] [PubMed]
  37. Sokolovic, S.; Alimanovic-Alagic, R.; Dzananovic, L.; Cavaljuga, S.; Beslic, N.; Ferhatbegovic-Opankovic, E. Vitamin D status in Bosnia and Herzegovina: The cross-sectional epidemiological analysis. Osteoporos. Int. 2017, 28, 1021–1025. [Google Scholar] [CrossRef] [PubMed]
  38. Dovnik, A.; Mujezinovic, F.; Treiber, M.; Pecovnik Balon, B.; Gorenjak, M.; Maver, U.; Takac, I. Determinants of maternal vitamin D concentrations in Slovenia: A prospective observational study. Wiener Klinische Wochenschr. 2017, 129, 21–28. [Google Scholar] [CrossRef] [PubMed]
  39. Stokic, E.; Kupusinac, A.; Tomic-Naglic, D.; Zavisic, B.K.; Mitrovic, M.; Smiljenic, D.; Soskic, S.; Isenovic, E. Obesity and vitamin D deficiency: Trends to promote a more proatherogenic cardiometabolic risk profile. Angiology 2015, 66, 237–243. [Google Scholar] [CrossRef] [PubMed]
  40. Sempos, C.T.; Vesper, H.W.; Phinney, K.W.; Thienpont, L.M.; Coates, P.M.; Vitamin, D.S.P. Vitamin D status as an international issue: National surveys and the problem of standardization. Scand. J. Clin. Lab. Investig. 2012, 243, 32–40. [Google Scholar] [CrossRef]
  41. Thienpont, L.M.; Stepman, H.C.; Vesper, H.W. Standardization of measurements of 25-hydroxyvitamin D3 and D2. Scand. J. Clin. Lab. Investig. 2012, 243, 41–49. [Google Scholar] [CrossRef]
  42. Franca Gois, P.H.; Wolley, M.; Ranganathan, D.; Seguro, A.C. Vitamin D Deficiency in Chronic Kidney Disease: Recent Evidence and Controversies. Int. J. Environ. Res. Public Health 2018, 15, 1773. [Google Scholar] [CrossRef] [PubMed]
  43. Mansbach, J.M.; Ginde, A.A.; Camargo, C.A., Jr. Serum 25-hydroxyvitamin D levels among US children aged 1 to 11 years: Do children need more vitamin D? Pediatrics 2009, 124, 1404–1410. [Google Scholar] [CrossRef] [PubMed]
  44. Kondric, M.; Trajkovski, B.; Strbad, M.; Foretic, N.; Zenic, N. Anthropometric influence on physical fitness among preschool children: Gender-specific linear and curvilinear regression models. Coll. Antropol. 2013, 37, 1245–1252. [Google Scholar] [PubMed]
  45. Dwyer, G.M.; Hardy, L.L.; Peat, J.K.; Baur, L.A. The validity and reliability of a home environment preschool-age physical activity questionnaire (Pre-PAQ). Int. J. Behav. Nutr. Phys. Act. 2011, 8, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Vollmer, R.L.; Adamsons, K.; Gorin, A.; Foster, J.S.; Mobley, A.R. Investigating the Relationship of Body Mass Index, Diet Quality, and Physical Activity Level between Fathers and Their Preschool-Aged Children. J. Acad. Nutr. Diet. 2015, 115, 919–926. [Google Scholar] [CrossRef] [PubMed]
  47. Skouteris, H.; Hill, B.; McCabe, M.; Swinburn, B.; Busija, L. A parent-based intervention to promote healthy eating and active behaviours in pre-school children: Evaluation of the MEND 2-4 randomized controlled trial. Pediatr. Obes. 2016, 11, 4–10. [Google Scholar] [CrossRef] [PubMed]
  48. Dwyer, G.M.; Walker, K.; Baur, L.; Badawi, N. Developmental outcomes and physical activity behaviour in children post major surgery: An observational study. BMC Pediatr. 2016, 16, 123. [Google Scholar] [CrossRef] [PubMed]
  49. Akeson, P.K.; Lind, T.; Hernell, O.; Silfverdal, S.A.; Ohlund, I. Serum Vitamin D Depends Less on Latitude Than on Skin Color and Dietary Intake During Early Winter in Northern Europe. J. Pediatr. Gastroenterol. Nutr. 2016, 62, 643–649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Katic, M.; Pirsl, F.; Steinberg, S.M.; Dobbin, M.; Curtis, L.M.; Pulanic, D.; Desnica, L.; Titarenko, I.; Pavletic, S.Z. Vitamin D levels and their associations with survival and major disease outcomes in a large cohort of patients with chronic graft-vs-host disease. Croat. Med. J. 2016, 57, 276–286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Culic, S.; Markic, J.; Petrovic, D.; Konjevoda, P.; Pavelic, J. Serum vitamin D levels in children with newly diagnosed and chronic immune thrombocytopenia. Semin. Hematol. 2016, 53 (Suppl. 1), S67–S69. [Google Scholar] [CrossRef]
  52. Zhu, Z.; Zhan, J.; Shao, J.; Chen, W.; Chen, L.; Li, W.; Ji, C.; Zhao, Z. High prevalence of vitamin D deficiency among children aged 1 month to 16 years in Hangzhou, China. BMC Public Health 2012, 12, 126. [Google Scholar] [CrossRef] [PubMed]
  53. Sanghera, D.K.; Sapkota, B.R.; Aston, C.E.; Blackett, P.R. Vitamin D Status, Gender Differences, and Cardiometabolic Health Disparities. Ann. Nutr. Metab. 2017, 70, 79–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Zhao, X.; Xiao, J.; Liao, X.; Cai, L.; Xu, F.; Chen, D.; Xiang, J.; Fang, R. Vitamin D Status among Young Children Aged 1–3 Years: A Cross-Sectional Study in Wuxi, China. PLoS ONE 2015, 10, e0141595. [Google Scholar] [CrossRef] [PubMed]
  55. Mandlik, R.; Kajale, N.; Ekbote, V.; Patwardhan, V.; Khadilkar, V.; Chiplonkar, S.; Khadilkar, A. Determinants of Vitamin D Status in Indian School-children. Indian J. Endocrinol. Metab. 2018, 22, 244–248. [Google Scholar] [CrossRef] [PubMed]
  56. Soininen, S.; Eloranta, A.M.; Lindi, V.; Venalainen, T.; Zaproudina, N.; Mahonen, A.; Lakka, T.A. Determinants of serum 25-hydroxyvitamin D concentration in Finnish children: The Physical Activity and Nutrition in Children (PANIC) study. Br. J. Nutr. 2016, 115, 1080–1091. [Google Scholar] [CrossRef] [PubMed]
  57. Pate, R.R.; Pfeiffer, K.A.; Trost, S.G.; Ziegler, P.; Dowda, M. Physical activity among children attending preschools. Pediatrics 2004, 114, 1258–1263. [Google Scholar] [CrossRef] [PubMed]
  58. Cleland, V.; Crawford, D.; Baur, L.A.; Hume, C.; Timperio, A.; Salmon, J. A prospective examination of children’s time spent outdoors, objectively measured physical activity and overweight. Int. J. Obes. 2008, 32, 1685–1693. [Google Scholar] [CrossRef] [PubMed]
  59. Kato, N.; Isojima, T.; Yokoya, S.; Tanaka, T.; Ono, A.; Yokomichi, H.; Yamagata, Z.; Tanaka, S.; Matsubara, H.; Ishikuro, M.; et al. Earlier BMI rebound and lower pre-rebound BMI as risk of obesity among Japanese preschool children. Int. J. Obes. 2018, 42, 52–58. [Google Scholar] [CrossRef] [PubMed]
  60. Alemzadeh, R.; Kichler, J.; Babar, G.; Calhoun, M. Hypovitaminosis D in obese children and adolescents: Relationship with adiposity, insulin sensitivity, ethnicity, and season. Metab. Clin. Exp. 2008, 57, 183–191. [Google Scholar] [CrossRef] [PubMed]
  61. Smotkin-Tangorra, M.; Purushothaman, R.; Gupta, A.; Nejati, G.; Anhalt, H.; Ten, S. Prevalence of vitamin D insufficiency in obese children and adolescents. J. Pediatr. Endocrinol. Metab. 2007, 20, 817–823. [Google Scholar] [CrossRef] [PubMed]
  62. Kumar, J.; Muntner, P.; Kaskel, F.J.; Hailpern, S.M.; Melamed, M.L. Prevalence and associations of 25-hydroxyvitamin D deficiency in US children: NHANES 2001–2004. Pediatrics 2009, 124, e362–e370. [Google Scholar] [CrossRef] [PubMed]
  63. Gilbert-Diamond, D.; Baylin, A.; Mora-Plazas, M.; Marin, C.; Arsenault, J.E.; Hughes, M.D.; Willett, W.C.; Villamor, E. Vitamin D deficiency and anthropometric indicators of adiposity in school-age children: A prospective study. Am. J. Clin. Nutr. 2010, 92, 1446–1451. [Google Scholar] [CrossRef] [PubMed]
  64. Motlaghzadeh, Y.; Sayarifard, F.; Allahverdi, B.; Rabbani, A.; Setoodeh, A.; Sayarifard, A.; Abbasi, F.; Haghi-Ashtiani, M.T.; Rahimi-Froushani, A. Assessment of Vitamin D Status and Response to Vitamin D3 in Obese and Non-Obese Iranian Children. J. Trop. Pediatr. 2016, 62, 269–275. [Google Scholar] [CrossRef] [PubMed]
  65. Rodríguez-Dehli, A.C.; Riaño-Galán, I.R.; Fernández-Somoano, A.; Navarrete-Muñoz, E.M.; Espada, M.; Vioque, J.; Tardón, A. Hipovitaminosis D y factores asociados a los 4 años en el norte de España. Anales de Pediatría 2017, 86, 188–196. [Google Scholar] [CrossRef] [PubMed]
  66. Klesges, L.M.; Baranowski, T.; Beech, B.; Cullen, K.; Murray, D.M.; Rochon, J.; Pratt, C. Social desirability bias in self-reported dietary, physical activity and weight concerns measures in 8- to 10-year-old African-American girls: Results from the Girls Health Enrichment Multisite Studies (GEMS). Prev. Med. 2004, 38, S78–S87. [Google Scholar] [CrossRef] [PubMed]
  67. Lips, P.; van Schoor, N. Worldwide vitamin D status. In Vitamin D, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2011; pp. 947–963. [Google Scholar]
Figure 1. 25(OH)D levels in preschool children with differences between genders (Chi square (χ2) and level of significance (p)).
Figure 1. 25(OH)D levels in preschool children with differences between genders (Chi square (χ2) and level of significance (p)).
Ijerph 15 02503 g001
Table 1. Descriptive statistics for the observed variables with differences between genders (t-test for independent samples).
Table 1. Descriptive statistics for the observed variables with differences between genders (t-test for independent samples).
TotalBoysGirlst-Test
(n = 260)(n = 132)(n = 128)
Mean ± SDMean ± SDMean ± SDt-Value (p)
25(OH)D (nmol/L)46.55 ± 20.1649.62 ± 19.6843.38 ± 20.242.52 (0.01)
BMI (kg/m2)15.91 ± 2.1615.91 ± 2.1415.92 ± 2.19−0.04 (0.97)
Pre-PAQ/L1 (score)1.15 ± 0.831.11 ± 0.791.18 ± 0.88−0.66 (0.51)
Pre-PAQ/L2 (score)0.80 ± 0.680.90 ± 0.690.71 ± 0.652.34 (0.02)
Pre-PAQ/L1 + L2 (score)1.95 ± 1.162.02 ± 1.151.89 ± 1.170.87 (0.38)
Pre-PAQ/L3 (score)0.96 ± 0.670.94 ± 0.690.98 ± 0.65−0.52 (0.60)
Pre-PAQ/L4 (score)0.92 ± 0.730.94 ± 0.770.89 ± 0.690.65 (0.52)
Pre-PAQ/L5 (score)0.32 ± 0.450.39 ± 0.490.26 ± 0.42.25 (0.03)
Time spent outdoors (hours/week)17.42 ± 7.1217.94 ± 7.6116.88 ± 6.571.21 (0.23)
Pre-PAQ—preschool-age physical activity questionnaire; Pre-PAQ/L1-L5—levels of physical activity as indicated by the Pre-PAQ; BMI—body mass index; SD—standard deviation.
Table 2. Univariate associations between independent variables and 25(OH)D level.
Table 2. Univariate associations between independent variables and 25(OH)D level.
Independent Variables25(OH)D (nmol/L)
TotalBoysGirls
(n = 260)(n = 132)(n = 128)
rprprp
BMI (kg/m2)−0.040.490.010.85−0.010.25
Pre-PAQ/L1 (score)−0.030.59−0.080.380.020.86
Pre-PAQ/L2 (score)−0.070.25−0.070.42−0.120.17
Pre-PAQ/L1 + L2 (score)−0.060.29−0.090.28−0.060.53
Pre-PAQ/L3 (score)−0.120.06−0.150.08−0.080.35
Pre-PAQ/L4 (score)−0.020.74−0.060.47−0.010.88
Pre-PAQ/L5 (score)−0.040.54−0.070.44−0.050.55
Time spent outdoors (hours/week)−0.050.420.010.90−0.150.09
r—Pearson’s product moment correlation coefficient; p—level of significance.
Table 3. Correlates of the 25(OH)D level in the total sample (n = 260). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
Table 3. Correlates of the 25(OH)D level in the total sample (n = 260). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
25(OH)D Deficiency (<50 nmol/L)25(OH)D Insufficiency (50–75 nmol/L)
Independent VariablesOR95% CIOR95% CI
Gender
Male0.330.15–0.740.600.25–1.43
FemaleREFREF
BMI cont1.080.89–1.311.080.88–1.31
Pre-PAQ/L1 cont1.150.71–1.861.000.59–1.70
Pre-PAQ/L2 cont1.490.81–2.751.140.59–2.20
Pre-PAQ/L1 + L2 cont1.220.86–1.751.050.72–1.54
Pre-PAQ/L3 cont1.400.76–2.580.950.49–1.85
Pre-PAQ/L4 cont1.150.68–1.940.880.49–1.56
Pre-PAQ/L5 cont2.350.84–6.681.920.65–5.69
Time spent outdoors cont1.040.98–1.101.030.97–1.10
OR—odds ratio; CI—confidence interval; cont indicates continuous variables.
Table 4. Correlates of the 25(OH)D level in boys (n = 132). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
Table 4. Correlates of the 25(OH)D level in boys (n = 132). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
25(OH)D Deficiency (<50 nmol/L)25(OH)D Insufficiency (50–75 nmol/L)
Independent VariablesOR95% CIOR95% CI
BMI1.040.82–1.311.130.88–1.44
Pre-PAQ/L11.330.71–2.490.880.44–1.76
Pre-PAQ/L21.210.59–2.471.230.58–2.61
Pre-PAQ/L1 + L21.220.79–1.881.020.64–1.63
Pre-PAQ/L31.450.69–3.040.850.37–1.95
Pre-PAQ/L41.050.57–1.930.830.43–1.62
Pre-PAQ/L51.960.63–6.121.880.57–6.16
Time spent outdoors1.010.94–1.081.020.95–1.09
Table 5. Correlates of the 25(OH)D level in girls (n = 128). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
Table 5. Correlates of the 25(OH)D level in girls (n = 128). 25(OH)D sufficiency (>75 nmol/L) was used as the reference value.
25(OH)D Deficiency (<50 nmol/L)25(OH)D Insufficiency (50–75 nmol/L)
Independent VariablesOR95% CIOR95% CI
BMI1.130.81–1.571.010.71–1.45
Pre-PAQ/L10.990.47–2.131.10.49–2.47
Pre-PAQ/L22.210.62–7.941.070.27–4.35
Pre-PAQ/L1 + L21.250.66–2.361.090.55–2.17
Pre-PAQ/L31.390.46–4.171.110.34–3.68
Pre-PAQ/L41.460.52–4.081.060.35–3.26
Pre-PAQ/L54.290.33–56.452.650.18–39.51
Time spent outdoors1.140.99–1.301.10.95–1.27

Share and Cite

MDPI and ACS Style

Karin, Z.; Gilic, B.; Supe Domic, D.; Sarac, Z.; Ercegovic, K.; Zenic, N.; Uljevic, O.; Peric, M.; Markic, J. Vitamin D Status and Analysis of Specific Correlates in Preschool Children: A Cross-Sectional Study in Southern Croatia. Int. J. Environ. Res. Public Health 2018, 15, 2503. https://doi.org/10.3390/ijerph15112503

AMA Style

Karin Z, Gilic B, Supe Domic D, Sarac Z, Ercegovic K, Zenic N, Uljevic O, Peric M, Markic J. Vitamin D Status and Analysis of Specific Correlates in Preschool Children: A Cross-Sectional Study in Southern Croatia. International Journal of Environmental Research and Public Health. 2018; 15(11):2503. https://doi.org/10.3390/ijerph15112503

Chicago/Turabian Style

Karin, Zeljka, Barbara Gilic, Daniela Supe Domic, Zdenko Sarac, Katarina Ercegovic, Natasa Zenic, Ognjen Uljevic, Mia Peric, and Josko Markic. 2018. "Vitamin D Status and Analysis of Specific Correlates in Preschool Children: A Cross-Sectional Study in Southern Croatia" International Journal of Environmental Research and Public Health 15, no. 11: 2503. https://doi.org/10.3390/ijerph15112503

APA Style

Karin, Z., Gilic, B., Supe Domic, D., Sarac, Z., Ercegovic, K., Zenic, N., Uljevic, O., Peric, M., & Markic, J. (2018). Vitamin D Status and Analysis of Specific Correlates in Preschool Children: A Cross-Sectional Study in Southern Croatia. International Journal of Environmental Research and Public Health, 15(11), 2503. https://doi.org/10.3390/ijerph15112503

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop