Vitamin D Abnormalities: More Common in US Children with T1D than in Healthy Children ()
1. Introduction
Adults living with type 1 diabetes mellitus (T1D) have decreased bone mineral density and are more likely to have osteoporosis compared with healthy controls [1] . This is significant for overall health as low bone mineral density is associated with an increased risk of fragility fractures, non-cause specific mortality, and cardiovascu- lar mortality [2] . The decreased bone mineral density associated with T1D may be due to a variety of factors, in- cluding earlier age at onset of diabetes [3] and coexisting conditions such as celiac disease [4] or thyroid disease [5] . Similarly, low serum 25-hydroxyvitamin D levels are associated with poor bone mineral density, and pa- tients with T1D may have low 25-hydroxyvitamin D levels, which may be most significant in patients with poor diabetes control [6] .
The majority of bone mineral density is accrued during childhood and adolescence, and failure to achieve adequate peak bone mass during childhood may lead to bone abnormalities and increased fracture risk during adulthood [7] . As many as 70% of healthy children in the United States have vitamin D abnormalities, and many may therefore not achieve adequate peak bone mass [8] . This risk may be even more pronounced in children with T1D, as there have been several recent studies reporting the high prevalence of vitamin D abnormalities in such children [9] - [11] . Additionally, low levels of serum 25-hydroxyvitamin D has recently been correlated with non-skeletal effects in children and adolescents with T1D. It has been reported that 25-hydroxyvitamin D defi- ciency is associated with a significantly higher insulin requirement [12] and an increased prevalence of retino- pathy [13] .
To better understand possible disparities, several recent studies have offered a comparison between subjects with T1D and control subjects without T1D, with mixed results. A study in Qatari youth with age-matched con- trols reported increased 25-hydroxyvitamin D deficiency/insufficiency with T1D [14] ; however this population may not be generalizable to US children due to a much higher overall prevalence of vitamin D deficiency/insuf- ficiency reported amongst the Qatari (90.6% in subjects with T1D compared to 85.3% in control subjects). In Australia, Greer et al. also reported a significant difference between groups [15] . In contrast, it has been reported in Denmark that no difference was detected between children with newly diagnosed T1D and their healthy sibl- ings [16] .
As serum 25-hydroxyvitamin D levels are reflective of dietary intake as well as sun exposure, geographical distribution of 25-hydryoxyvitamin D abnormalities varies significantly; even within North America, results have varied. A study performed in Florida reported no significant difference between T1D and control groups [17] ; however, this study included children, adolescents, and adults, thus preventing accurate comparison. In contrast, a more recent study evaluating Canadian adolescents and teenagers reported 32.8% of individuals with established T1D have 25-hydroxyvitamin D deficiency, compared to only 11.8% of the similar-aged healthy Canadian population [18] . Given these conflicting results within North American populations, additional clari- fication is needed to better understand the relationship between T1D and 25-hydroxyvitamin D in this region.
Therefore, we explored whether pediatric subjects with T1D differ from a control population of pediatric sub- jects living in middle Tennessee, within the United States. We hypothesize that the prevalence of 25-hydroxy- vitamin D deficiency/insufficiency in children with T1D is higher than in healthy children and that diabetes du- ration, HbA1c, and insulin dose/kg is inversely associated with 25-hydroxyvitamin D levels.
2. Methods
2.1. Study Design
The inclusion criteria for this cross-sectional study were patients who were less than 21 years of age diagnosed with T1D who receive care at Monroe Carell Jr. Children’s Hospital at Vanderbilt University Medical Center. Only patients of diabetes providers who screen all patients’ serum 25-hydroxyvitamin D levels were included, to alleviate bias. No patients were excluded from this study. Data was collected from the medical record, and the following characteristics were obtained for each patient within the past 2 years: height, weight, age, gender, race/ethnicity, date of diabetes diagnosis, and the presence of additional diagnoses including: thyroid disease, celiac disease, Turner Syndrome, Trisomy 21, and Addison’s disease. All medications were recorded including type of insulin utilized, total insulin dose/kg, method of insulin delivery (injection versus pump) and use of sup- plementation with vitamin D or multivitamin. Laboratory values included a serum 25-hydroxyvitamin D level and all HbA1c values from the two years prior to 25-hydroxyvitamin D measurement.
The control population data was collected from a de-identified database of 100 healthy pediatric patients who live in similar middle-Tennessee geography, as previously reported [19] . This population included individuals between the ages of 1 - 21 years who live in similar middle-Tennessee geography and who were having blood- work obtained at an outpatient clinic at Vanderbilt Children’s Hospital. Control subjects were excluded if they had the following diagnoses: osteopenia, osteoporosis, osteogenesis imperfecta, previously diagnosed vitamin D insufficiency/deficiency, more than two bone fractures during the prior year, chronic or current glucocorticoid use, thyrotoxicosis, gastrointestinal disease such as celiac disease causing possible malabsorption, diabetes mel- litus, history of or current treatment for malignancy, or were of non-weight bearing status. No patient approach- ed about the study declined to participate. This dataset contains age, gender, race/ethnicity, height, weight, use of vitamin supplementation, and a single serum 25-hydroxyvitamin D level for each patient.
All data were stored in a secured, de-identified password-protected REDcap database [20] . The study was ap- proved by the Institutional Review Board.
2.2. Anthropometric Evaluation
Body mass index was calculated based on the formula: weight (kg)/height (m2). BMI Z-scores were calculated using age and gender-standardized growth population norms based on the Centers for Disease Control and Pre- vention’s Year 2000 Growth Charts.
2.3. Laboratory Evaluation
The vitamin D measurement utilized for both groups was 25-hydroxyvitamin D, which is the standard indicator of vitamin D status, and levels were measured via liquid chromatography/ tandem mass spectrometry (Mayo Medical Laboratories, Rochester, MN). The total 25-hydroxyvitamin D concentrations were calculated by sum- ming the measured values of 25-hydroxyvitamin D2 and 25-hydroxyvitamin D3. Intra- and inter-assay coeffi- cients of variation have previously been reported to be <7% [21] . The criteria used to determine 25-hydroxyvi- tamin D deficiency, insufficiency, and sufficiency were 25-hydroxyvitamin D levels of <20 ng/mL, 20 - 29 ng/mL, and ≥30 ng/mL respectively [22] . HbA1c was measured by point-of-care immunoassay using the DCA Vantage Analyzer (Siemens Healthcare Diagnostics, Deerfield, IL, USA).
2.4. UV Index
Ambient ultraviolet (UV) light was examined by using the reported UV index for Nashville, TN. Records of the UV index for Nashville were obtained from the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center (CPC) website [23] . Reported UV indices issued for the 30 days preceding serum 25- hydroxyvitamin D measurement were averaged, and this value is referred to as the UV30.
2.5. Statistical Methods
Descriptive statistics were provided for both subject groups. A contingency table was created with a corres- ponding Chi-square test of association to compare the prevalence of 25-hydroxyvitamin D insufficiency and de- ficiency among T1D and normal controls. The alpha-value employed was 0.05. Multiple linear regression was performed to regress 25-hydroxyvitamin D levels on groups, adjusting for race, age, gender, UV index, BMI Z-score, and multivitamin supplementation. A similar multiple linear regression model was used to assess the association of the covariates with the 25-hydroxyvitamin D levels within the T1D patient group. Restricted cu- bic splines were initially considered for all continuous variables (age, UV index, BMI Z-score, HbA1c values, duration of diabetes and insulin dose/kg) to allow for non-linear associations with 25-hydroxyvitamin D levels, but the splines were reduced to linear coefficients if there was not strong evidence of non-linearity.
3. Results
A total of 276 subjects with T1D and 100 healthy pediatric subjects without T1D were evaluated in this study. Of the patients with T1D, there were 268 with established T1D, while 8 were considered new-onset (less than 90 days from diagnosis). T1D patients were analyzed collectively. Participant characteristics are presented in Table 1. Subjects with T1D were more likely than subjects without T1D to be older (median of 14 vs. 11.2 years, p < 0.001), Caucasian (85% vs. 67%, p < 0.001), and to have a lower UV30 (5.7 vs. 6.8, p = 0.021). Subjects with T1D were less likely to report regular supplementation with a multivitamin compared to those without T1D (11% vs. 28%, p < 0.001).
Within the T1D group, the median length of time since diagnosis was 4.9 years (IQR: 2.8, 9.2 years). The median HbA1c during the 2 years prior to the time of 25-hydroxyvitamin D measurement and the median HbA1c value at the time of 25-hydroxyvitamin D measurement were the same [8.6% (IQR 7.9%, 9.4%) vs. 8.6% (IQR: 7.9%, 9.8%), respectively]; 48% of subjects used injections while the remaining 52% used an insulin pump for delivery. Additional diagnosis included: thyroid disease (n = 39), celiac disease (n = 17), Addison’s disease (n = 3), Turner syndrome (n = 2), and Trisomy 21 (n = 1).
The T1D population had a significantly lower 25-hydroxyvitamin D level than the control population, 24 ng/mL versus 31 ng/mL (p < 0.001, Table 2) univariately. The T1D subjects also had significantly higher per- centages of 25-hydroxyvitamin D deficiency and insufficiency compared to healthy controls. In the T1D group, 32% were 25-hydroxyvitamin D deficient, 36% insufficient and 33% sufficient. In the control group, 14% were deficient, 30% insufficient, and 56% sufficient (p < 0.001, Table 2, Figure 1).
Results shown as median and [IQR] interquartile range. *UV30: A measurement of the average UV index (amount of UV light) in the patient’s geo- graphic location during the 30 days prior to the serum 25-hydroxyvitamin D evaluation.
![]()
Table 2. Serum 25-hydroxyvitamin D status: Pediatric TID patients compared with healthy pediatric control patients.
Results shown as median and [IQR] interquartile range.
![]()
Figure 1. Descriptive chart showing serum 25-hydroxyvitamin D status between children and adolescents with T1D (n = 276) versus pediatric control po- pulation (n = 100). Categorization based on Endocrine Society Clinical Practice Guidelines: Deficient is defined as 25-hydroxyvitamin D level < 20 ng/mL; insufficient as 20 - 29 ng/ mL; sufficient as ≥30 ng/ mL.
In the multivariable linear model adjusting for covariates, the 25-hydroxyvitamin D levels for T1D patients were on average 6.67 ng/mL lower than those of normal controls, with a 95% confidence interval of (4.13, 9.22) and p value of <0.001 (Table 3). Caucasian race and UV30 were positively associated with 25-hydroxyvitamin D levels (p < 0.001 and p = 0.0006, respectively). There were significant linear and non-linear associations of 25-hydroxyvitamin D with both age [p = 0.002 and p = 0.002 (CI −0.14, 0.00), Figure 2] and BMI Z-score [p = 0.005 and p = 0.002 (CI −5.99, −1.62), Figure 3]. No significant association between 25-hydroxyvitamin D le- vels and gender or reported multivitamin/vitamin D supplementation was detected.
For the subset analysis of T1D subjects, none of the continuous explanatory variables showed evidence of non-linear associations with 25-hydroxyvitamin D level and hence were modeled with linear coefficients. After adjusting for other variables in the model, Caucasians had elevated 25-hydroxyvitamin D levels compared to African Americans and other races (p = 0.001) (Table 4). No statistically significant association was found be- tween serum 25-hydroxyvitamin D levels and gender, body mass index, UV30, duration of diabetes, HbA1c value, use of injection versus pump, total insulin dose/kg, or supplementation with a multivitamin.
An additional subset analysis was performed comparing patients with T1D and additional comorbid autoim- mune diagnoses and those with only T1D. The median 25-hydroxyvitamin D level in patients with T1D and other autoimmune diagnoses (hypothyroidism, n = 39; celiac disease, n = 17; and Addison’s disease, n = 3) was 24.0 ng/mL (20.0, 32.0) and the median 25-hydroxyvitamin D level in patients with T1D but without other au- toimmune diagnoses (n = 224) was 23.0 ng/mL (18.8, 33.0), p = 0.73. Of the 52 patients with additional au- toimmune diagnoses, four had both hypothyroidism and celiac disease and all three patients who had Addison’s disease also had hypothyroidism. No patients had both Addison’s disease and celiac disease.
4. Discussion
As 25-hydroxyvitamin D levels are generally dependent upon both dietary intake and ultraviolet exposure, geo- graphic location has a significant effect on 25-hydroxyvitamin D level. It is therefore imperative to evaluate the prevalence of 25-hydroxyvitamin D deficiency/insufficiency based on regional or local data when available. When sub-populations are considered, it is preferably with comparison to a control population from the same area. We report that children and adolescents with T1D have a significantly higher prevalence of 25-hydrox- yvitamin D deficiency and insufficiency compared with a similar-aged population without diabetes living in the United States, in middle Tennessee (68% versus 44%).
This finding refutes the previous US based study performed in Florida, where Bierschenk et al. reported no particular association between reduced serum 25-hydroxyvitamin D levels and T1D [17] . Bierschenk reported a median 25-hydroxyvitamin D level of 23.2 ng/mL, consistent with our reported median of 24 ng/mL. The dif-
![]()
Figure 2. Predicted serum 25-hydroxyvitamin D level by age. Predicted means plot for age: Adjusted to: case = T1D, BMI z-score 0.6533, Caucasian race, male sex, UV30 = 5.648, no reported vitamin D supplement use.
![]()
Table 3. All patients: Multivariable analysis of serum 25-hydroxyvitamin D.
*UV30: A measurement of the average UV index (amount of UV light) in the patient’s geographic location during the 30 days prior to the serum 25-hydroxyvitamin D evaluation; **These are significant non-linear associations; see Figure 2 for predicted means plot for age and Figure 3 for pre- dicted means plot for BMI z-score.
![]()
Table 4. T1D patients: Multivariable analysis of serum 25-hydroxyvitamin D.
*UV30: A measurement of the average UV index (amount of UV light) in the patient’s geographic location during the 30 days prior to the serum 25-hydroxyvitamin D evaluation.
![]()
Figure 3. Predicted serum 25-hydroxyvitamin D level for BMI z-score. Predicted means plot for BMI z-score: Adjusted to: case = T1D, 156 months, Caucasian race, male sex, UV30 = 5.648, no reported vitamin D supplement use.
ference was found in healthy controls, with a median of 20.1 ng/mL versus 31 ng/mL, respectively. This discre- pancy is likely explained by the age difference in the study reported by Bierschenk; the control population held a median age of 22 years (range 5 - 65.1 years) versus 16 years (range 5.1 - 62.6 years) for patients with estab- lished T1D. Our study includes only pediatric patients, allowing a more accurate comparison. Thus, social and recreational habits influenced by age likely explain the different outcomes.
Yeshayahu et al. report 25-hydroxyvitamin D deficiency in 32.8% of pediatric patients with T1D living in Canada [18] , which is similar to the prevalence in the current study in middle Tennessee (32%). However, the middle Tennessee population with T1D who are sufficient in 25-hydroxyvitamin D is different than the Cana- dian population: 33% in Tennessee versus 14% in Canada. This disparity was equally true for healthy popula- tions also: 56% versus 35.2%. This is likely due to the significant difference in availability of vitamin D from sunlight between Canada and Tennessee.
In our total group analysis, UV30 was significantly associated with serum 25-hydroxyvitamin D. However, in the multivariable analysis of only patients with T1D, there was no association between UV30 and serum 25-hy- droxyvitamin D after adjusting for other covariates. A lack of association between UV30 and 25-hydroxyvita- min D was similarly reported by Greer et al. [15] . In this study, it could indicate that children and adolescents with T1D may have decreased time outdoors, as the UV30 is only a measure of ambient UV light and does not necessarily represent any individual’s time spent outside or amount of skin exposed to sunlight. Alternatively, the smaller sample size of the subset analysis of T1D subjects resulted in wider 95% confidence intervals, the- reby not excluding the possibility of a clinically meaningful association in the population. In the majority of stu- dies of 25-hydroxyvitamin D abnormalities, investigators have typically adjusted for seasonality effects using three-month intervals (e.g. summer, fall, winter, spring). In this study, we have chosen to use UV30 instead of seasons, with the justification that seasons can be viewed as a crude categorization of the ultraviolet index, and hence lose information relative to the UV30.
Consistent with other studies, we report that serum 25-hydroxyvitamin D levels are higher amongst Caucasian children as compared to African American and non-white children [8] [9] [24] . We found a similar association between BMI z-score and vitamin D status as reported in previous studies; both underweight [22] and over- weight [8] [25] [26] status are associated with lower serum 25-hydroxyvitamin D levels. Additionally, our data suggest that vitamin D levels may have a bimodal distribution and are higher in the youngest and oldest children, which has been reported previously [27] . However, this is in contrast to other studies which have reported that vitamin D levels are lowest amongst youngest children and teenagers [8] [9] [24] [28] .
We did not find a significant association between supplementation with vitamin D or multivitamin and 25- hydroxyvitamin D levels, in contrast to some studies [8] [29] [30] . This is likely due, in part, to the low number of subjects reporting vitamin D supplementation in both groups as well as the variability of vitamin D dosing between patients. In addition, we were unable to ascertain with certainty the frequency of patient supplementa- tion. A treatment study with varying dosing regimens and careful evaluation of compliance should be performed to better evaluate this possibility. Although the frequency of multivitamin or other vitamin D supplementation was low in both groups, supplementation was significantly lower amongst those with T1D. As patients with T1D have to adhere to a complex insulin treatment regimen, less concern may be given to other aspects of health perceived by parents as non-essential.
The T1D population had numerous co-morbid conditions, as there were 39 patients with a diagnosis of thyroid disease (14.1%), 17 with celiac disease (6.2%), 3 with Addison’s disease (1.1%), 2 with Turner syndrome (0.7%), and 1 with Trisomy 21 (0.4%). It has been well-described that patients with T1D are at increased risk for other such genetic and autoimmune disorders [31] -[33] . Our findings are similar to the recently reported coexistence of each of the aforementioned diagnoses within a population of young adults with T1D as follows: thyroid dis- ease (9.5%), celiac disease (7%), Addison’s disease (0.7%), Turner syndrome (3.5%), and Trisomy 21 (0.7%) [34] . Several of these disorders are also associated with decreased bone mineral density due to either disease course or therapy [4] [5] [35] [36] , warranting special concern in such patients. An additional multivariable analysis of the subjects with T1D and either thyroid disease, celiac disease, Addison’s disease, Turner Syndrome, or Trisomy 21 revealed no significant change in results. Therefore, these additional diagnoses are not associated with lower levels of 25-hydroxyvitamin D in T1D subjects.
The relationship between 25-hydroxyvitamin D and T1D is complex and not fully understood. Vitamin D may be protective against development of T1D, as a recent meta-analysis demonstrated decreased risk for T1D in subjects who received vitamin D supplementation during early life compared with those who did not receive vitamin D supplementation (pooled OR 0.71, p < 0.001) [37] . Moreover, it has been suggested that there is a dose-dependent positive correlation between vitamin D supplementation at a young age and decreased risk for T1D [38] . Evaluation of this relationship in US military personnel revealed that non-Hispanic white individuals with 25-hydroxyvitamin D levels of ≥100 nmol/L had a 44% lower risk of developing T1D than those with le- vels < 75 nmol/L (p = 0.03) [39] . Similarly, a recent study reported that 25-hydroxyvitamin D levels are lower in children with pre-T1D as compared to autoantibody-negative children [40] . Thus, there is compelling evidence that low 25-hydroxyvitamin D may reflect the primary immunological disorder. 25-hydroxyvitamin D may have a role in regulation of the immune system including T and B lymphocytes, promoting a more tolerant immuno- logical status [37] [41] . Additionally, 25-hydroxyvitamin D may have direct effects on B-cells, such as improv- ing insulin secretion and improving islet morphology [41] [42] .
In conclusion, children and adolescents with T1D have significantly increased prevalence of 25-hydroxyvi- tamin D deficiency/insufficiency compared to a similarly aged control population in middle-Tennessee, and this is not moderated by control, duration, or insulin management of diabetes. This may increase risk of complica- tions such as diabetic retinopathy, and have a significant impact upon bone health in patients with T1D. These youth may not attain adequate peak bone mass, the effects of which may contribute to increased risk for osteo- porosis, fragility fractures, non-cause specific mortality, and cardiovascular mortality as adults [2] . Further re- search should be pursued to determine whether decreased serum 25-hydroxyvitamin D can be alleviated via die- tary or behavioral modifications in this population.