Background: Vitamin D deficiency is common in children and may adversely affect linear growth, yet regional data from coastal western India are limited.Objectives: To estimate the prevalence of vitamin D deficiency and examine its association with growth parameters among children aged 1–10 years in Sindhudurg, Maharashtra.Methods: A hospital-based cross-sectional study was conducted over 6 months at SSPM Medical College and Lifetime Hospital, Padve. In 400 children (1–10 years), serum 25-hydroxyvitamin D [25(OH)D] was measured and categorized as deficient (<20 ng/mL), insufficient (20–29 ng/mL), or sufficient (≥30 ng/mL). Anthropometry was used to derive HAZ, WAZ, and BAZ; determinants and growth associations were assessed using multivariable regression.Results: Vitamin D deficiency was present in 40.3% (161/400), insufficiency in 25.3%, and sufficiency in 34.5%. Deficient children had lower mean HAZ than sufficient children (−1.06 ± 1.07 vs −0.73 ± 1.17; p=0.028). In adjusted linear models, deficiency remained associated with lower HAZ (β −0.37, 95% CI −0.64 to −0.10; p=0.008) and WAZ (β −0.29, 95% CI −0.55 to −0.04; p=0.024), with no association with BAZ. Higher sun exposure was independently protective against deficiency (aOR 0.77 per hour/day, 95% CI 0.60–0.99; p=0.042). Conclusion: Nearly two-thirds of children had suboptimal vitamin D status, and deficiency was independently associated with lower linear growth indices. Improving effective sun exposure and appropriate supplementation may reduce deficiency burden and support child growth.
Vitamin D plays a central role in calcium–phosphate homeostasis and skeletal mineralization during childhood, and deficiency has been increasingly recognized as a common nutritional problem even among apparently healthy pediatric populations. Indian studies have consistently reported a substantial burden of hypovitaminosis D across diverse pediatric age groups, underscoring that deficiency is not limited to children with overt illness or severe malnutrition. In a study from North India, Angurana et al. reported a high prevalence of vitamin D deficiency in apparently healthy children, highlighting the hidden nature of this micronutrient deficit and the likelihood of under-detection in routine care settings [1].
Despite wide use of serum 25-hydroxyvitamin D (25[OH]D) as the biomarker of vitamin D status, clinically meaningful thresholds for deficiency remain debated across populations and outcomes. Some studies have attempted to derive cutoffs using metabolic risk markers rather than bone endpoints alone. For example, Sharifi et al. proposed a deficiency cutoff anchored to insulin resistance metrics in children, illustrating that vitamin D status may have implications beyond skeletal growth and that threshold selection can influence estimated prevalence [2]. Such variability in definitions, alongside differences in geography, sunlight exposure, diet, and supplementation practices, contributes to heterogeneity in reported prevalence and complicates comparisons between regions.
Beyond prevalence estimation, an important public health question is whether vitamin D deficiency is associated with impaired growth in children. Evidence from outside India suggests that low vitamin D status may be linked to adverse anthropometric outcomes. Mokhtar et al. demonstrated an association between vitamin D status and underweight and stunting among young children in the Ecuadorian Andes, providing biologic and epidemiologic support for evaluating growth correlates of vitamin D in pediatric populations [3]. However, the strength and pattern of association may vary by setting due to differences in baseline nutrition, infection burden, socioeconomic gradients, and sunlight availability.
Within India, multicentric data have reinforced that vitamin D deficiency is widespread and shaped by multiple determinants. Khadilkar et al., in a large multicentre study of Indian children and adolescents, described both the distribution of vitamin D status and key determinants such as demographic and environmental factors, emphasizing that risk is multifactorial and not explained by a single exposure [4]. These findings further support the need for region-specific estimates, particularly in underrepresented settings such as coastal and semi-rural districts, where sunlight patterns, clothing practices, dietary habits, and socioeconomic profiles may differ from large metropolitan cohorts. The same multicentre evidence base also underscores the importance of simultaneously assessing determinants (e.g., sun exposure, diet, supplementation, socioeconomic position) when examining growth correlates, to reduce confounding and improve interpretability of observed associations [5].
Sindhudurg district in coastal Maharashtra has distinctive environmental and lifestyle features that may influence pediatric vitamin D status and growth. Yet, there remains limited local evidence quantifying deficiency and evaluating its relationship with growth parameters among children in the 1–10 year age range. Therefore, this cross-sectional study was designed to estimate the prevalence of vitamin D deficiency in children aged 1–10 years attending SSPM Medical College and Lifetime Hospital, Padve, and to examine the association of vitamin D status with anthropometric growth indicators, while accounting for key demographic and exposure-related determinants.
Objectives
Aim
To determine the prevalence of vitamin D deficiency and evaluate its association with growth parameters among children aged 1–10 years attending SSPM Medical College and Lifetime Hospital, Padve, Sindhudurg, Maharashtra.
Specific objectives
Study design and setting A hospital-based cross-sectional study was conducted over 6 months at SSPM Medical College and Lifetime Hospital, Padve, Sindhudurg, Maharashtra. Participants and eligibility Children aged 1–10 years attending the pediatric outpatient department and/or admitted for non-critical indications during the study period were eligible. Written informed consent was obtained from parents/guardians prior to enrollment. Children were excluded if they had conditions known to affect vitamin D metabolism or growth, including chronic renal or liver disease, malabsorption syndromes, endocrine disorders, genetic syndromes influencing growth, current or recent therapeutic (high-dose) vitamin D treatment, or critical illness requiring emergency resuscitation. Sample size and sampling A total sample of 400 children was enrolled using consecutive sampling over the study period. Data collection A structured case record form was used to collect demographic and exposure information, including age, sex, residence (rural/urban), socioeconomic status, season of sampling, average daily sun exposure (hours/day), outdoor activity level, dietary frequency (milk/egg/fish intake expressed as days per week), and history of vitamin D supplementation in the preceding three months. Anthropometry and growth assessment Weight was measured using a calibrated digital weighing scale with children in light clothing and without footwear. Standing height was measured using a stadiometer; where age-appropriate, recumbent length was measured using a length board. Body mass index (BMI) was calculated as weight (kg) divided by height squared (m²). Age- and sex-standardized indices—height-for-age Z-score (HAZ), weight-for-age Z-score (WAZ), and BMI-for-age Z-score (BAZ)—were derived using standard reference algorithms. Nutritional status was categorized as stunting (HAZ < −2), underweight (WAZ < −2), thinness (BAZ < −2), overweight (BAZ > +1), and obesity (BAZ > +2). Laboratory assessment and definitions Venous blood samples were collected under aseptic precautions. Serum 25-hydroxyvitamin D [25(OH)D] was measured in the institutional laboratory using a standardized assay with routine internal quality checks. Vitamin D status was defined as deficiency (<20 ng/mL), insufficiency (20–29 ng/mL), and sufficiency (≥30 ng/mL). Where available from routine evaluation, serum calcium and alkaline phosphatase were also recorded. Outcomes The primary outcome was the prevalence of vitamin D deficiency (<20 ng/mL). Secondary outcomes included associations between vitamin D status and growth parameters (height, weight, BMI and derived HAZ/WAZ/BAZ), associations with nutritional status categories (stunting, underweight, thinness, overweight/obesity), and determinants of vitamin D deficiency. Statistical analysis Continuous variables were summarized as mean ± standard deviation and categorical variables as number (percentage). Prevalence estimates were reported with 95% confidence intervals. Comparisons across vitamin D categories were performed using one-way ANOVA for continuous variables and chi-square tests for categorical variables. Multivariable linear regression was used to evaluate associations between vitamin D status and Z-scores (HAZ, WAZ, BAZ), while multivariable logistic regression was used to assess determinants of vitamin D deficiency and associations with binary nutritional outcomes. Models were adjusted for prespecified covariates including age, sex, residence, socioeconomic status, season, sun exposure, outdoor activity, dietary frequency variables, and recent vitamin D supplementation. A two-sided p value <0.05 was considered statistically significant. Ethics The study protocol was approved by the Institutional Ethics Committee. Confidentiality was maintained using de-identified study codes.
1.Participant characteristics (N=400)
A total of 400 children aged 1–10 years were included in the analysis. Participant demographics, sampling season, sun exposure profile, dietary frequency indicators, and supplementation history are summarized in Table 1.
Table1. Participant characteristics (N=400)
|
Characteristic |
Overall (N=400) |
|
Demographics |
|
|
Age (years), mean ± SD |
4.88 ± 2.69 |
|
Age group 1–5 years, n (%) |
237 (59.2) |
|
Age group 6–10 years, n (%) |
163 (40.8) |
|
Male sex, n (%) |
199 (49.8) |
|
Socio-demographic profile |
|
|
Rural residence, n (%) |
288 (72.0) |
|
Socioeconomic status, n (%) |
|
|
Lower |
126 (31.5) |
|
Lower-middle |
132 (33.0) |
|
Middle |
94 (23.5) |
|
Upper-middle/Upper |
48 (12.0) |
|
Season and sunlight exposure |
|
|
Season of sampling, n (%) |
|
|
Monsoon |
196 (49.0) |
|
Post-monsoon |
99 (24.8) |
|
Winter |
105 (26.2) |
|
Sun exposure (hours/day), mean ± SD |
1.66 ± 0.87 |
|
Sun exposure category (hours/day), n (%) |
|
|
<0.5 |
37 (9.2) |
|
0.5–1.5 |
142 (35.5) |
|
1.6–3.0 |
203 (50.7) |
|
>3.0 |
18 (4.5) |
|
Lifestyle and diet |
|
|
Outdoor activity, n (%) |
|
|
Low |
117 (29.2) |
|
Moderate |
205 (51.2) |
|
High |
78 (19.5) |
|
Vitamin D supplementation (last 3 months), n (%) |
82 (20.5) |
|
Milk intake (days/week), mean ± SD |
4.76 ± 1.86 |
|
Egg intake (days/week), mean ± SD |
1.84 ± 1.38 |
|
Fish intake (days/week), mean ± SD |
1.25 ± 1.03 |
|
Note: |
|
|
Values are n (%) unless stated otherwise. Sun exposure categories are in hours/day; diet variables are days/week. |
|
Overall, 161/400 children (40.3%) were vitamin D deficient (25[OH]D <20 ng/mL; 95% CI 35.6–45.1). An additional 101/400 (25.3%) had insufficiency (20–29 ng/mL; 95% CI 21.2–29.7), while 138/400 (34.5%) were sufficient (≥30 ng/mL; 95% CI 30.0–39.3). The distribution of vitamin D status is presented in Table 2and shown graphically in Figure 1
Table 2. Prevalence of vitamin D status categories among children aged 1–10 years (N=400).
(Serum 25[OH]D classified as deficient <20 ng/mL, insufficient 20–29 ng/mL, and sufficient ≥30 ng/mL; values shown as n, % and 95% confidence intervals.)
|
Vitamin D status |
n |
Percent |
95% CI (low) |
95% CI (high) |
|
Deficient (<20) |
161 |
40.2 |
35.6 |
45.1 |
|
Insufficient (20–29) |
101 |
25.2 |
21.2 |
29.7 |
|
Sufficient (≥30) |
138 |
34.5 |
30.0 |
39.3 |
Figure 1. Prevalence of vitamin D deficiency, insufficiency, and sufficiency among children aged 1–10 years (N=400).
3.Baseline characteristics by vitamin D status
Baseline characteristics stratified by vitamin D category are summarized in Table 3. Mean daily sun exposure differed across vitamin D groups, with children who were vitamin D deficient reporting lower sun exposure (1.56 ± 0.85 h/day) compared with those who were vitamin D sufficient (1.82 ± 0.82 h/day; p = 0.027). Age distribution, sex, residence, socioeconomic status, season of sampling, and reported dietary frequency (milk/egg/fish) were broadly comparable across vitamin D categories (all p > 0.05). Milk intake frequency showed a borderline trend toward lower intake in deficient children (4.56 ± 1.94 days/week) versus sufficient children (5.05 ± 1.72 days/week; p = 0.065), and outdoor activity pattern also showed a borderline difference (p = 0.064). Reported vitamin D supplementation in the preceding 3 months did not differ significantly by vitamin D status (p = 0.144).
Table 3. Baseline characteristics of children aged 1–10 years stratified by vitamin D status (N=400).
|
Characteristic |
Deficient (<20) |
Insufficient (20–29) |
Sufficient (≥30) |
p-value |
|
Age (years), mean ± SD |
5.18 ± 2.63 |
4.54 ± 2.72 |
4.78 ± 2.73 |
0.154 |
|
Sun exposure (hours/day), mean ± SD |
1.56 ± 0.85 |
1.62 ± 0.92 |
1.82 ± 0.82 |
0.027 |
|
Milk intake (days/week), mean ± SD |
4.56 ± 1.94 |
4.67 ± 1.89 |
5.05 ± 1.72 |
0.065 |
|
Egg intake (days/week), mean ± SD |
1.88 ± 1.43 |
1.78 ± 1.29 |
1.85 ± 1.39 |
0.850 |
|
Fish intake (days/week), mean ± SD |
1.22 ± 1.02 |
1.26 ± 1.07 |
1.31 ± 1.01 |
0.733 |
|
Sex, n (%) |
|
|
|
0.752 |
|
— Male |
83 (51.6) |
48 (47.5) |
68 (49.3) |
|
|
— Female |
78 (48.4) |
53 (52.5) |
70 (50.7) |
|
|
Residence, n (%) |
|
|
|
0.572 |
|
— Rural |
119 (73.9) |
74 (73.3) |
95 (68.8) |
|
|
— Urban |
42 (26.1) |
27 (26.7) |
43 (31.2) |
|
|
Socioeconomic status, n (%) |
|
|
|
0.315 |
|
— Lower |
53 (32.9) |
37 (36.6) |
36 (26.1) |
|
|
— Lower-middle |
59 (36.6) |
31 (30.7) |
42 (30.4) |
|
|
— Middle |
34 (21.1) |
22 (21.8) |
38 (27.5) |
|
|
— Upper-middle/Upper |
15 (9.3) |
11 (10.9) |
22 (15.9) |
|
|
Season of sampling, n (%) |
|
|
|
0.072 |
|
— Monsoon |
83 (51.6) |
56 (55.4) |
57 (41.3) |
|
|
— Post-monsoon |
36 (22.4) |
24 (23.8) |
39 (28.3) |
|
|
— Winter |
42 (26.1) |
21 (20.8) |
42 (30.4) |
|
|
Outdoor activity, n (%) |
|
|
|
0.064 |
|
— Low |
55 (34.2) |
33 (32.7) |
29 (21.0) |
|
|
— Moderate |
73 (45.3) |
45 (44.6) |
87 (63.0) |
|
|
— High |
33 (20.5) |
23 (22.8) |
22 (15.9) |
|
|
Vitamin D supplementation in last 3 months, n (%) |
25 (15.5) |
25 (24.8) |
32 (23.2) |
0.144 |
|
Overweight (BAZ > +1), n (%) |
9 (5.6) |
10 (9.9) |
7 (5.1) |
0.272 |
Note: Values are mean ± SD or n (%). p-values from one-way ANOVA (continuous variables) or χ² test (categorical variables).
4.Growth parameters and nutritional status by vitamin D status
Growth parameters stratified by vitamin D category are shown in Table 4. Mean HAZ differed significantly across categories, with children who were vitamin D deficient showing lower HAZ than those who were sufficient (−1.06 ± 1.07 vs −0.73 ± 1.17; p = 0.028). Mean WAZ showed a similar direction but did not reach conventional statistical significance (p = 0.058), while BAZ did not differ (p = 0.747). The proportions of stunting, underweight, thinness, overweight, and obesity were not significantly different across vitamin D categories (all p > 0.05). Biochemical correlates showed expected variation: mean serum calcium differed across groups (p < 0.001) and alkaline phosphatase was higher in the deficient group (p < 0.001) (Table 4).
Table 4. Growth parameters and nutritional status stratified by vitamin D status (N=400)
|
Parameter / outcome |
Deficient (<20) (n=161) |
Insufficient (20–29) (n=101) |
Sufficient (≥30) (n=138) |
p-value |
|
25(OH)D (ng/mL), mean ± SD |
14.28 ± 3.49 |
24.98 ± 3.25 |
36.37 ± 6.29 |
<0.001 |
|
Height (cm), mean ± SD |
105.87 ± 17.39 |
103.00 ± 18.13 |
104.97 ± 18.30 |
0.449 |
|
Weight (kg), mean ± SD |
17.32 ± 6.28 |
16.68 ± 6.69 |
17.54 ± 7.07 |
0.598 |
|
BMI (kg/m²), mean ± SD |
15.05 ± 2.62 |
15.27 ± 2.81 |
15.35 ± 2.70 |
0.597 |
|
HAZ, mean ± SD |
−1.06 ± 1.07 |
−0.78 ± 1.10 |
−0.73 ± 1.17 |
0.028 |
|
WAZ, mean ± SD |
−1.11 ± 0.97 |
−0.86 ± 1.16 |
−0.83 ± 1.12 |
0.058 |
|
BAZ, mean ± SD |
−0.51 ± 1.02 |
−0.43 ± 1.04 |
−0.43 ± 0.99 |
0.747 |
|
Stunting (HAZ < −2), n (%) |
32 (19.9) |
13 (12.9) |
20 (14.5) |
0.257 |
|
Underweight (WAZ < −2), n (%) |
35 (21.7) |
21 (20.8) |
20 (14.5) |
0.244 |
|
Thinness (BAZ < −2), n (%) |
8 (5.0) |
6 (5.9) |
4 (2.9) |
0.498 |
|
Overweight (BAZ > +1), n (%) |
9 (5.6) |
10 (9.9) |
7 (5.1) |
0.272 |
|
Obesity (BAZ > +2), n (%) |
2 (1.2) |
0 (0.0) |
1 (0.7) |
0.525 |
|
Calcium (mg/dL), mean ± SD |
9.24 ± 0.37 |
9.45 ± 0.41 |
9.35 ± 0.40 |
<0.001 |
|
ALP (IU/L), mean ± SD |
271.92 ± 63.94 |
240.50 ± 73.75 |
246.60 ± 73.70 |
<0.001 |
Note: Continuous variables compared by one-way ANOVA; categorical outcomes compared by χ² test. Values are mean ± SD or n (%).
Abbreviations: HAZ (Height-for-age Z-score), WAZ (Weight-for-age Z-score), BAZ (BMI-for-age Z-score)
5.Adjusted association between vitamin D status and growth outcomes
In multivariable linear regression (adjusted for age, sex, residence, socioeconomic status, season, sun exposure, outdoor activity, diet frequency, and recent supplementation), vitamin D deficiency was independently associated with lower HAZ and lower WAZ compared with vitamin D sufficiency. Specifically, deficient children had a 0.37 SD lower HAZ (adjusted β −0.37, 95% CI −0.64 to −0.10; p=0.008) and a 0.29 SD lower WAZ (adjusted β −0.29, 95% CI −0.55 to −0.04; p=0.024). There was no adjusted association with BAZ (p>0.05). Vitamin D insufficiency was not significantly associated with HAZ, WAZ, or BAZ after adjustment (all p>0.05).
In adjusted logistic regression models, vitamin D deficiency showed higher odds of stunting, underweight, and thinness compared with sufficiency; however, these associations did not reach statistical significance and confidence intervals were wide (Table 5).
Table 5. Adjusted association of vitamin D status with growth indicators and undernutrition outcomes (N=400)
Table 5A. Adjusted linear regression (reference: Sufficient ≥30 ng/mL)
Outcome = Z-score (HAZ/WAZ/BAZ); values are adjusted β (95% CI)
|
Outcome |
Deficient (<20) vs Sufficient |
p-value |
Insufficient (20–29) vs Sufficient |
p-value |
|
HAZ |
−0.37 (−0.64 to −0.10) |
0.008 |
−0.08 (−0.37 to 0.22) |
0.617 |
|
WAZ |
−0.29 (−0.55 to −0.04) |
0.024 |
−0.04 (−0.34 to 0.26) |
0.795 |
|
BAZ |
−0.05 (−0.29 to 0.19) |
0.694 |
−0.01 (−0.27 to 0.26) |
0.963 |
Table 5B. Adjusted logistic regression (reference: Sufficient ≥30 ng/mL)
Outcome = binary nutritional status; values are adjusted OR (95% CI)
|
Outcome |
Deficient (<20) vs Sufficient |
p-value |
Insufficient (20–29) vs Sufficient |
p-value |
|
Stunting (HAZ < −2) |
1.55 (0.81 to 2.97) |
0.187 |
0.92 (0.42 to 1.99) |
0.829 |
|
Underweight (WAZ < −2) |
1.69 (0.90 to 3.17) |
0.104 |
1.56 (0.78 to 3.12) |
0.208 |
|
Thinness (BAZ < −2) |
1.80 (0.44 to 7.32) |
0.409 |
3.26 (0.72 to 14.67) |
0.124 |
Note: Models adjusted for age, sex, residence, SES, season, sun exposure, outdoor activity, milk/egg/fish intake frequency, and recent vitamin D supplementation.
6.Determinants of vitamin D deficiency
In a multivariable logistic regression model with vitamin D deficiency (<20 ng/mL) as the outcome (adjusted for age, sex, residence, socioeconomic status, season, sun exposure, outdoor activity, diet frequency, and recent supplementation), greater daily sun exposure was independently associated with lower odds of deficiency (aOR 0.77 per 1 hour/day increase; 95% CI 0.60–0.99; p=0.042). Compared with children reporting high outdoor activity, those with moderate outdoor activity had lower odds of deficiency (aOR 0.56; 95% CI 0.32–0.96; p=0.036). Recent vitamin D supplementation showed a protective direction but did not meet conventional significance (aOR 0.59; 95% CI 0.35–1.02; p=0.057). Other demographic variables (age, sex, residence, SES, season) and dietary frequency measures did not show statistically significant independent associations (Table 6).
Table 6. Multivariable determinants of vitamin D deficiency (<20 ng/mL) among children aged 1–10 years (N=400)
|
Predictor |
Adjusted OR (aOR) |
95% CI |
p-value |
|
Age (per 1-year increase) |
1.06 |
0.98–1.14 |
0.146 |
|
Male (vs Female) |
1.00 |
0.66–1.53 |
0.982 |
|
Urban residence (vs Rural) |
0.96 |
0.60–1.53 |
0.863 |
|
SES: Lower-middle (vs Lower) |
1.24 |
0.74–2.09 |
0.421 |
|
SES: Middle (vs Lower) |
0.85 |
0.48–1.52 |
0.588 |
|
SES: Upper-middle/Upper (vs Lower) |
0.77 |
0.38–1.58 |
0.479 |
|
Season: Post-monsoon (vs Monsoon) |
0.82 |
0.50–1.36 |
0.455 |
|
Season: Winter (vs Monsoon) |
0.70 |
0.42–1.18 |
0.184 |
|
Outdoor activity: Low (vs High) |
0.65 |
0.35–1.20 |
0.170 |
|
Outdoor activity: Moderate (vs High) |
0.56 |
0.32–0.96 |
0.036 |
|
Sun exposure (per 1 hour/day increase) |
0.77 |
0.60–0.99 |
0.042 |
|
Vitamin D supplementation in last 3 months (Yes vs No) |
0.59 |
0.35–1.02 |
0.057 |
|
Milk intake (per 1 day/week increase) |
0.91 |
0.81–1.02 |
0.090 |
|
Egg intake (per 1 day/week increase) |
0.96 |
0.83–1.11 |
0.606 |
|
Fish intake (per 1 day/week increase) |
0.94 |
0.77–1.16 |
0.573 |
Note: Outcome = vitamin D deficiency (<20 ng/mL). Reference categories: Female, Rural, Lower SES, Monsoon, High outdoor activity.
Vitamin D deficiency was common in children aged 1–10 years in coastal Sindhudurg, with 40.3% deficient and 65.6% below sufficiency. Deficiency was independently associated with lower linear growth indices (HAZ and WAZ), and greater sun exposure was protective against deficiency. These findings support strengthening routine screening/targeted testing in high-risk children and emphasizing pragmatic exposure and supplementation strategies alongside broader nutrition programs, while larger community-based longitudinal studies are needed to confirm growth impacts and guide region-specific policy.