25(OH)D Levels in Infancy Is Associated With Celiac Disease Autoimmunity in At-Risk Children: A Case–Control Study

Objectives: An observed variation in the risk of celiac disease, according to the season of birth, suggests that vitamin D may affect the development of the disease. The aim of this study was to investigate if vitamin D concentration is associated with the risk of celiac disease autoimmunity (CDA) in genetically at-risk children.

Study Design: Children prospectively followed in the multinational The Environmental Determinants of Diabetes in the Young study, conducted at six centers in Europe and the US, were selected for a 1-to-3 nested case–control study. In total, 281 case–control sets were identified. CDA was defined as positivity for tissue transglutaminase autoantibodies (tTGA) on two or more consecutive visits. Vitamin D was measured as 25-hydroxyvitamin D [25(OH)D] concentrations in all plasma samples prior to, and including, the first tTGA positive visit. Conditional logistic regression was used to examine the association between 25(OH)D and risk of CDA.

Results: No significant association was seen between 25(OH)D concentrations (per 5 nmol/L increase) and risk for CDA development during early infancy (odds ratio [OR] 0.99, 95% confidence interval [CI] 0.95–1.04) or childhood (OR 1.02, 95% CI 0.97–1.07). When categorizing 25(OH)D concentrations, there was an increased risk of CDA with 25(OH)D concentrations <30 nmol/L (OR 2.23, 95% CI 1.29, 3.84) and >75 nmol/L (OR 2.10, 95% CI 1.28–3.44) in early infancy, as compared with 50–75 nmol/L.

Conclusion: This study indicates that 25(OH)D concentrations <30 nmol/L and >75 nmol/L during early infancy were associated with an increased risk of developing CDA in genetically at-risk children. The non-linear relationship raises the need for more studies on the possible role of 25(OH)D in the relation to celiac disease onset.

Introduction

Vitamin D is essential for bone growth and the effective functioning of innate immunity (1). Most of the circulating vitamin D, 25-hydroxyvitamin D, [25(OH)D] is synthesized by the skin with the help of UVB radiation from sunlight. Populations living in the northern hemisphere often need to compensate for the lower amount of sunlight via fortified foods or dietary supplements (2). Interestingly, vitamin D sufficiency, compared with lower than normal plasma levels, was associated with decreased risk of developing autoimmunity or immune-mediated diseases such as type 1 diabetes (T1D) (3–5), multiple sclerosis, rheumatoid arthritis, and Crohn's disease (6–8). Celiac disease is a chronic enteropathy with autoimmune features caused by an immune-mediated response to dietary gluten, leading to the destruction of intestinal mucosa resulting in malabsorption (9). Celiac disease autoimmunity (CDA) is defined as the presence of tissue transglutaminase autoantibodies (tTGA), indicative of an ongoing gluten-induced inflammatory response, which often precedes small bowel mucosal damage and a celiac disease diagnosis (10, 11). Although celiac disease is strongly associated with human leucocyte antigen (HLA) DQA1^*05:01-DQB1^*02 and DQA1^*03:01-DQB1^*03:02 haplotypes, these risk genes cannot fully explain the risk, suggesting that environmental factors also contribute (9). The most promising promoters so far are gastrointestinal infections (12, 13), high gluten intake (14), or both (15). For reasons yet unknown, children born during the spring and summer months are at the highest disease risk (16–20). This has led to the hypothesis that vitamin D deficiency in early life may predispose to celiac disease due to seasonal differences in UVB exposure and subsequent 25(OH)D concentrations or via dysregulation of the immune response leading to an abnormal intestinal mucosa with increasing permeability (21, 22). Although low 25(OH)D concentrations have been reported at the time of celiac disease diagnosis (23–26), this can be attributed to deranged dietary absorption from a damaged gut epithelium.

The aim of this study was to investigate if the levels of 25(OH)D before the disease onset are associated with increased risk of CDA in a prospective birth cohort of genetically at-risk children.

Materials and Methods

Study Population

The Environmental Determinants of Diabetes in the Young (TEDDY) is a prospective birth cohort consisting of 8,676 genetically at-risk children born between September 2004 and February 2010. Children were enrolled in the study before the age of 4.5 months and followed for 15 years to identify environmental triggers of T1D and celiac disease (27, 28). Children carrying high-risk HLA alleles for T1D and celiac disease were enrolled at six centers, three in the US (Colorado, Washington, and Georgia/Florida) and three in Europe (Finland, Germany, and Sweden). The following HLA-class II genotypes: HLA-DR3/4; HLA-DR4/4, HLA-DR4/8, HLA-DR3/3, and HLA-DR4/4 were the eligibility criteria for enrollment in the study. Children with HLA-DR4/1, HLA-DR4/13, HLA-DR4/9, and HLA-DR3/9 were included if they had a first degree relative (FDR) i.e. having a mother, father, or sibling with T1D.

For all study participants, separate written informed consent was obtained from a parent or primary caretaker for genetic screening and participation in the prospective follow-up beginning at birth. The study was conducted according to the guidelines of the Declaration of Helsinki, and local institutional or regional ethics review boards in all participating countries approved the study.

Study Design

The present study was performed by using children included in two nested case–control cohorts that were constructed with a focus on islet autoimmunity (IA) and T1D, and which aimed to study multiple biomarkers that are expensive to measure in a large birth cohort setting. The design and planning of the TEDDY nested case–control biomarker study have been described in detail elsewhere (29). Cases and controls were identified as of May 31, 2012, and were matched by clinical center, sex, and family history of T1D. All available samples meeting the design criteria by that time were processed in the laboratories chosen for each biomarker analysis. All children in the 1:3 nested case–control studies for vitamin D biomarker that had been screened for CDA as of August 31, 2017 were considered for the present study. Each case–control set included a CDA positive child (“case”) matched up to three controls. All controls were tTGA-negative for at least 6 months from the age of seroconversion of the case. Subjects having at least one 25(OH)D measurement prior or at seroconversion visit of the matched cases were included in the final analysis, resulting in 281 cases and 643 controls (123 cases with three matched controls, 116 cases with two controls, and 42 cases with one control) (Supplementary Figure 1). Maternal and early infant feeding information were collected prospectively at clinic visits every 3 months. Descriptive characteristics of the study population are presented in Table 1.

TABLE 1

Table 1. Descriptive characteristics of subjects in the nested case–control study.

Screening for CDA and Celiac Disease

Annual screening for celiac disease started from the age of 2 years using radiobinding assays (RBA) to measure tTGA in serum. Samples were analyzed at the Bristol Laboratory (University of Bristol, UK) for the European sites and the Barbara Davis Center Laboratory (Aurora, Colorado) for the US sites. In the US, the RBA uses anti-IgA agarose to capture IgA-TGA, whereas, in Bristol, a mixture of both anti-IgA agarose and protein A sepharose is used to assess both IgA-tTG and IgG-tTG. Children from the US sites with a tTGA level > 0.05 were deemed antibody positive. To harmonize the protocol, all sera with tTGA levels > 0.01 at the laboratory in Colorado were reassayed at the laboratory in Bristol for confirmation of being tTGA positive. Children with tTGA levels >1.3 unit in Bristol were deemed antibody positive. All earlier samples of children who were tTGA positive at age 24 months were retrospectively analyzed in the Bristol laboratory to determine the age of seroconversion to tTGA positivity.

Celiac disease autoimmunity was defined as being tTGA positive in two consecutive samples, collected 3–6 months apart. A total of 281 children developed CDA at a median age of the first positive tTGA at 2.8 years of age (IQR; Q1: 2.0, Q3: 3.9 years).

Vitamin D Assessment

Plasma from blood samples was drawn into light-protected tubes (BD Vacutainer®CPT™ cell preparation tubes) and analyzed at the Institute for Health and Welfare, Helsinki, Finland. 25(OH)D concentration was measured using the ARCHITECT 25-OH vitamin D chemiluminescent microparticle immunoassay (3). The laboratory participates in the vitamin D external quality assessment scheme. Plasma samples for 25(OH)D were collected from visits at 3, 6, 9, and 12 months of age and then annually up or until and including the time of seroconversion for CDA-cases. Vitamin D status was categorized according to cut-offs previously used in pediatric populations, which are as follows: <30 nmol/L, 30–50 nmol/L, 50–75 nmol/L, and >75 nmol/L (30, 31).

Genotyping

Genotyping was performed by the University of Virginia using a custom Illumina Infinium ImmunoChip (Illumina, Inc; CA). The ImmunoChip was designed to genotype immune-mediated disease loci identified by genome wide association studies in 12 autoimmune diseases (including celiac disease and T1D). Genes in the vitamin D pathway may modify the efficiency of 25(OH) concentrations; therefore several single nucleotide polymorphisms (SNPs) associated with vitamin D metabolism, such as GC (rs7041), VDR (rs1544410 [Bsml], rs11568820 [Cdx2], rs7975232 [Apal]), CYP27B1 (rs4646536), CYP24A1 (rs4809959, rs2616277), and RXRA (rs3818740, rs10881582) were selected for the analyses as previously described (3).

Statistical Analyses

Conditional logistic regression was used to examine the association between 25(OH)D concentrations and the risk of developing CDA. The magnitudes of the associations were described by odds ratios (ORs) with 95% confidence intervals (CIs). Having an FDR with celiac disease, HLA genotype associated with celiac disease and factors associated with 25(OH)D concentrations including season of drawing blood or season of birth, long-distance protocol (LDP; that is, samples collected locally and shipped to a TEDDY clinic for processing), and IA status (3) were adjusted in the models.

Concentrations of 25(OH)D were examined in early infancy by using the first available sample (i.e., early infancy concentrations) and during childhood, defined as average concentrations of all visits prior and including the seroconversion visit of the case (i.e., childhood concentrations). The early infancy and childhood 25(OH)D concentrations were incorporated into the conditional logistic regression models in two ways, namely: (a) The 25(OH)D concentration as a continuous variable; (b) the 25(OH)D concentration was grouped into four categories as <30, 30–50, 50–75, and >75 nmol/L.

When examining the a priori proposed 25(OH)D vitamin D gene interactions with the 25(OH)D concentration, each SNP was analyzed individually by treating the number of minor alleles as a continuous variable, and an interaction term between the SNP and 25(OH)D concentrations was included in the model. Ancestry (population structure) was adjusted in the interaction analyses by using the two largest principal components from a principal component analysis (PCA) of the ImmunoChip data in the cohort, using KING software (32).

All analyses were performed using SAS (Version 9.4; SAS Institute, Cary, North Carolina) version 9.4. A two-sided value of p < 0.05 was considered statistically significant.

Results

Distribution of 25(OH)D Concentrations

The proportion of children classified as being vitamin D deficient (<30 nmol/L) during early infancy or during childhood were 114 (13.7%) and 49 (5.3%) children, respectively. Early infancy and childhood 25(OH)D concentrations of cases and control subjects according to the matching variables in the study cohort are presented in Supplementary Table 1. Differences in 25(OH)D concentrations were most noticeable between countries and clinical centers, where subjects from Germany and Sweden had the highest concentrations in both early infancy and during childhood. No other major differences in 25(OH)D concentrations for sex and having an FDR with T1D at the two exposure points were observed.

25(OH)D Concentrations and Risk of CDA

There was no association between 25(OH)D concentrations and risk for CDA development during early infancy (OR 0.99, 95% CI 0.95–1.04 5 per nmol/L increase) or childhood (OR 1.02, 95% CI 0.97–1.07 5 per nmol/L increase), after adjusting for HLA, FDR with celiac disease, season of birth, season of drawing blood, being on LDP, and IA status (Table 2). There was no interaction between 25(OH)D concentrations and a clinical center on the risk of CDA (P = 0.30).

TABLE 2

Table 2. Association between 25(OH)D concentrations categories and celiac disease autoimmunity (CDA) in children in the nested case–control study.

When classifying 25(OH)D concentrations into four categories, there was an increased risk of CDA with 25(OH)D concentrations <30 nmol/L (OR 2.23, 95% CI 1.29, 3.84) and >75 nmol/L (OR 2.10, 95% CI 1.28–3.44) in early infancy compared with 25(OH)D concentrations between 50 and 75 nmol/L (Table 2). No significant association was identified for average childhood 25(OH)D concentrations (Table 2). Using conditional logistic regression analysis with smoothing splines (33, 34), a non-linear relationship between vitamin D concentrations in early infancy and the risk of CDA was observed (p = 0.014 for nonlinearity) (Figure 1). The non-linear relationship is consistent with the finding that both <30 nmol/L and >75 nmol/L concentrations in early infancy were associated with a higher risk of CDA compared with 50–75 nmol/L (i.e., a U-shape relationship). Similar non-linear relationships were observed when analyzing the first available sample at 3, 6, and 9 months, respectively (data not shown). In a sensitivity analysis, additional adjustment for gluten consumption prior to seroconversion (14) and being rotavirus vaccinated (12) changed the association for 25(OH)D in early infancy and CDA marginally. Gastrointestinal episodes (12) were not taken into consideration since they could be a mediator rather than a confounder (low vitamin D was associated with higher odds of gastrointestinal infections, p = 0.017).

FIGURE 1

Figure 1. The estimated effect [solid line, with pointwise standard error (dashed lines)] of 25(OH)D concentration in early infancy on the log hazard of CDA from conditional logistic regression analysis with smoothing splines (p = 0.014 for non-linearity).

There was no interaction between any of the selected vitamin D pathway SNPs and 25(OH)D concentration, analyzed as both continuous and categorical variables, on CDA (Table 3).

TABLE 3

Table 3. Test for interaction between vitamin D gene SNPs and 25(OH)D concentrations on the risk of celiac disease autoimmunity in the TEDDY nested case–control study.

To identify potential sociodemographic or other confounding factors explaining the findings, descriptive characteristics of the four 25(OH)D concentration categories are presented in Supplementary Table 3. Among subjects with 25(OH)D concentrations <30 nmol/L, a higher portion of the children were from the sites in Colorado (USA) and Finland, were persistent IA-positive, exclusively breastfed, born during the fall season (September–November), and the blood was drawn during the winter months (December–February). Subjects with 25(OH)D concentrations >75 nmol/L, a higher proportion were from Sweden, reported vitamin D supplementation at the time of drawing the blood, or were on an LDP. Also, only a few subjects were exclusively breastfed and blood drawn was collected during the winter months (December to February).

Discussion

This nested case–control study found an association of low levels of 25(OH)D concentrations in early infancy with an increased risk of CDA in genetically at-risk children. This finding may shed light on the role of vitamin D in the development of celiac disease. The exact mechanism for this association is not clear, but it could be attributed to the active form of vitamin D (1,25-dihydroxyvitamin D) which inhibits the production of proinflammatory cytokines, possibly leading to a shift from an anti-inflammatory state to a more inflammatory state (35). This speculation would be in line with a previous finding of significantly higher levels for IL-6, IL-13, IL-10, IL-1b, and TNF-α in combination with significantly lower 25(OH)D concentrations in screening-detected celiac disease cases compared with healthy controls (24). These differences were no longer present when comparing celiac disease patients on a gluten-free diet compared with controls. A Norwegian study compared 25(OH)D cord plasma levels between children who later developed celiac disease and healthy controls, and found that levels were highly correlated with maternal 25(OH)D at the time of delivery but not related to later development of celiac disease (36). Drawing conclusions on this is difficult as the 25(OH)D concentrations of the infant at birth are highly dependent on maternal vitamin D status. Early life exposures such as vitamin D supplementation may influence the developmental imprinting of the immune system. Vitamin D is able to modulate the innate immune system and enable the system to fight against pathogens, but it may be so that the infants acquired 25(OH)D status during the first year of life may be more important for shaping the adaptive immune system and risk for later autoimmunity (1, 37).

It has been speculated that earlier identified non-linear relationships between 25(OH)D and different disease outcomes may be related to SNPs and their influence on vitamin D status (38, 39). Most of the circulating 25(OH)D is bound to vitamin D binding protein genes. Identified genes directly involved in the vitamin D metabolism and associated with abnormal 25(OH)D concentrations are VDR, GC, CYP2R1, CYP24A1, DHCR7, and RXRA (40). The functions of these genes can be expressed as a lower affinity for 25(OH)D or impaired 25-hydoxylase activity. In the present study, we included the most critical SNPs involved in the vitamin D pathway as previously described (3), but found no interaction between any of the selected vitamin D pathway SNPs and 25(OH)D, analyzed as both continuous and categorical variables. This is in line with findings in a previous study on neonatal vitamin D status in relation to the later development of celiac disease (36).

Many studies have reported an association between the risk of celiac disease and the season of birth (primarily summer births) (17–19, 41). One possible explanation for this association is the difference in received sunlight in pregnant mothers and spring-born infants, leading to lower 25(OH)D concentrations during the second half of infancy (during fall and winter). The lower vitamin D status in these infants may also coincide with a time when there are frequent seasonal infections and the age when gluten is introduced (~6 months after birth). An alternate hypothesis is that the observed seasonal variation is due to the differences in viral infections, but it may also be that low levels of 25(OH)D increase the susceptibility to infections affecting the gut, causing a disruptive barrier to triggering antigens involved in the celiac disease (22, 42).

A more inexplicable finding from the present study was that 25(OH)D concentrations >75 nmol/L during early infancy increased the risk of CDA in childhood. Concentrations exceeding 75 nmol/L are most likely due to frequent vitamin D supplementation. Interestingly, one previous study reported children who received vitamin D supplementation for longer than 3 months to be at an increased risk of developing celiac disease (43). The proposed explanation is that high doses of vitamin D upregulate Th2 cell cytokines associated with immune reaction to external stimuli (44). In the present study, 92% of the children exceeding 25(OH)D concentrations >75 nmol/L reported vitamin D supplementation at the time of drawing the blood compared with 55% among children with low 25(OHD) concentrations. Differences in 25(OH)D concentrations may also be due to strength, frequency, and the duration of vitamin D supplementation. At the time of the study enrollment, there were differences among countries regarding recommendations for vitamin D supplementation to infants (45). In Sweden and Germany, the vitamin D recommendations (10 μg/day, starting within 2 months of birth) did not change during the study period. In Finland and the US, the recommendations changed from 5 to 10 μg during the corresponding period. Results from the full study cohort showed that nearly 80% of the participants (European sites; 97–99%) started vitamin D supplementation within the 1st year of life (45). Of note, the start of vitamin D supplementation may be a confounder in our analyses because information about the length of supplementation before the time of blood draw was not available. Vitamin D supplementation in higher doses may also be a result of parents knowing the genetic risk of their infants for T1D and supplementing their child as an action to prevent disease. However, in TEDDY, less than 3% of the mothers reported giving dietary supplements as an action to prevent T1D before the age of 2 (46).

At present, there is no clear consensus regarding optimal 25(OH)D concentrations in the general pediatric population. Defining appropriate levels of 25(OH)D has been problematic, and recommendations are mostly based on studies on fracture risk and poorly accounted for the risk of immune-mediated diseases in children.

This study has some limitations. First, the study was originally designed to investigate several dietary biomarkers and its association with IA and T1D status. However, the number of IA-positive subjects was well-balanced between cases (29%) and controls (23%) and adjusted for in the analyses. Another limitation is that a major proportion of 25(OH)D samples close to seroconversion, defined as within 1 year before seroconversion to tTGA positivity, were missing. Samples drawn before seroconversion were only available in <50% of the subjects, resulting in lower power in the statistical tests.

It is debated whether circulating total 25(OH)D is the best marker of vitamin D status compared with the biological active form of vitamin D [1.25(OH)₂D] and free (unbound to protein) 25(OH)D (47). However, total 25(OH)D measured in plasma or serum is an accepted indicator of vitamin D status, often used in epidemiological research, and it reflects both endogenous and dietary sources. It has a half-life of 2–3 weeks and is above all mostly used as a measure of vitamin D deficiency. However, the active form of vitamin D has an ability to modulate both innate and adaptive immunity, with pro-and anti-inflammatory actions that may be related to celiac disease (48).

The strengths of this study are the large study cohort, with study participants from six clinical sites in four countries, following the same study protocol to prospectively collect blood samples from birth up until 10 years of age before seroconversion to CDA.

This study adds additional knowledge regarding 25(OH)D status in genetically at-risk infants but raises the need for more studies and preferable randomized controlled trials to further investigate the role of 25(OH)D concentrations in early childhood and its relation to celiac disease onset.

In conclusion, this study shows that 25(OH)D concentrations <30 and >75 nmol/L during early infancy are associated with an increased risk of developing CDA in genetically at-risk children, indicating a possible role of 25(OH)D in the development of celiac disease.

Data Availability Statement

The data presented in the study are deposited in the NIDDK Central Repository at https://repository.niddk.nih.gov/studies/teddy/.

Ethics Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and local institutional or regional ethics review boards in all participating countries approved the study. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

CA: designed the study, drafted the initial manuscript, interpretation of data, and reviewed and revised the manuscript. DA: conceptualized and designed the study and critical review of both intermediary and final versions of the manuscript. WH, MR, J-XS, JT, A-GZ, BA, and JK: coordinated and supervised data collection, interpretation of data, and critical review of both intermediary and final versions of the manuscript. KK, SK, SV, and IE: interpretation of data and critical review of both intermediary and final versions of the manuscript. UU and MB: accusation of data collection, interpretation of data, and critical review of both intermediary and final versions of the manuscript. JN: conceptualized the study, interpretation of data, and critical review of both intermediary and final versions of the manuscript. XL: carried out the statistical analyses, interpretation of data, and reviewed and revised the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work, ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

TEDDY Study Group

Colorado Clinical Center: Marian Rewers, M.D., Ph.D., PI, Aaron Barbour, Kimberly Bautista, Judith Baxter, Daniel Felipe-Morales, Brigitte I. Frohnert, M.D., Marisa Stahl, M.D., Patricia Gesualdo, Rachel Haley, Michelle Hoffman, Rachel Karban, Edwin Liu, M.D., Alondra Munoz, Jill Norris, Ph.D., Stesha Peacock, Hanan Shorrosh, Andrea Steck, M.D., Megan Stern, Kathleen Waugh. University of Colorado, Anschutz Medical Campus, Barbara Davis Center for Childhood Diabetes.

Finland Clinical Center: Jorma Toppari, M.D., Ph.D., PI^¥∧, Olli G. Simell, M.D., Ph.D., Annika Adamsson, Ph.D.^∧, Sanna-Mari Aaltonen^∧, Suvi Ahonen^*±§, Mari Åkerlund^*±§, Leena Hakola^*±, Anne Hekkala, M.D.^μ¤, Henna Holappa^μ¤, Heikki Hyöty, M.D., Ph.D.^*±, Anni Ikonen^μ¤, Jorma Ilonen, M.D., Ph.D.^¥¶, Sanna Jokipuu^∧, Leena Karlsson^∧, Jukka Kero M.D., Ph.D.^¥∧, Miia Kähönen^μ¤, Mikael Knip, M.D., Ph.D.^*±, Minna-Liisa Koivikko^μ¤, Katja Kokkonen^*±, Merja Koskinen^*±, Mirva Koreasalo^*±§, Kalle Kurppa, M.D., Ph.D.^*±, Salla Kuusela, M.D.^*±, Jarita Kytölä^*±, Sinikka Lahtinen^*±, Jutta Laiho, Ph.D.^*, Tiina Latva-aho^μ¤, Laura Leppänen^∧, Katri Lindfors, Ph.D.^*, Maria Lönnrot, M.D., Ph.D.^*±, Elina Mäntymäki^∧, Markus Mattila^*±, Maija Miettinen^§, Katja Multasuo^μ¤, Teija Mykkänen^μ¤, Tiina Niininen^±*, Sari Niinistö^§, Mia Nyblom^*±, Sami Oikarinen, Ph.D.^*±, Paula Ollikainen^μ¤, Zhian Othmani^¥, Sirpa Pohjola ^μ¤, Jenna Rautanen^±§, Anne Riikonen^*±§, Minna Romo^∧, Satu Simell, M.D., Ph.D.^¥, Aino Stenius^μ¤, Päivi Tossavainen, M.D.^μ¤, Mari Vähä-Mäkilä^¥, Eeva Varjonen^∧, Riitta Veijola, M.D., Ph.D.^μ¤, Irene Viinikangas^μ¤, Suvi M. Virtanen, M.D., Ph.D.^*±§. ^¥University of Turku, ^*Tampere University, ^μUniversity of Oulu, ^∧Turku University Hospital, Hospital District of Southwest Finland, ^±Tampere University Hospital,^¤Oulu University Hospital, §Finnish Institute for Health and Welfare, Finland, ^¶University of Kuopio.

Georgia/Florida Clinical Center: Jin-Xiong She, Ph.D., PI, Desmond Schatz, M.D.^*, Diane Hopkins, Leigh Steed, Jennifer Bryant, Katherine Silvis, Michael Haller, M.D.^*, Melissa Gardiner, Richard McIndoe, Ph.D., Ashok Sharma, Stephen W. Anderson, M.D.^∧, Laura Jacobsen, M.D.^*, John Marks, DHSc.^*, P.D. Towe^*. Center for Biotechnology and Genomic Medicine, Augusta University. ^*University of Florida, Pediatric Endocrinology. ^∧Pediatric Endocrine Associates, Atlanta.

Germany Clinical Center: Anette G. Ziegler, M.D., PI, Ezio Bonifacio Ph.D.^*, Cigdem Gezginci, Anja Heublein, Eva Hohoff^¥, Sandra Hummel, Ph.D., Annette Knopff, Charlotte Koch, Sibylle Koletzko, M.D.^¶, Claudia Ramminger, Roswith Roth, Ph.D., Jennifer Schmidt, Marlon Scholz, Joanna Stock, Katharina Warncke, M.D., Lorena Wendel, Christiane Winkler, Ph.D. Forschergruppe Diabetes e.V. and Institute of Diabetes Research, Helmholtz Zentrum München, Forschergruppe Diabetes, and Klinikum rechts der Isar, Technische Universität München. ^*Center for Regenerative Therapies, TU Dresden, ^¶Dr. von Hauner Children's Hospital, Department of Gastroenterology, Ludwig Maximillians University Munich, ^¥University of Bonn, Department of Nutritional Epidemiology.

Sweden Clinical Center: Åke Lernmark, Ph.D., PI, Daniel Agardh, M.D., Ph.D., Carin Andrén Aronsson, Ph.D., Maria Ask, Rasmus Bennet, Corrado Cilio, Ph.D., M.D., Susanne Dahlberg, Malin Goldman Tsubarah, Emelie Ericson-Hallström, Annika Björne Fors, Lina Fransson, Thomas Gard, Monika Hansen, Susanne Hyberg, Berglind Jonsdottir, M.D., Ph.D., Helena Elding Larsson, M.D., Ph.D., Marielle Lindström, Markus Lundgren, M.D., Ph.D., Marlena Maziarz, Ph.D., Maria Månsson Martinez, Jessica Melin, Zeliha Mestan, Caroline Nilsson, Yohanna Nordh, Kobra Rahmati, Anita Ramelius, Falastin Salami, Anette Sjöberg, Carina Törn, Ph.D., Ulrika Ulvenhag, Terese Wiktorsson, Åsa Wimar. Lund University.

Washington Clinical Center: William A. Hagopian, M.D., Ph.D., PI, Michael Killian, Claire Cowen Crouch, Jennifer Skidmore, Christian Chamberlain, Brelon Fairman, Arlene Meyer, Jocelyn Meyer, Denise Mulenga, Nole Powell, Jared Radtke, Shreya Roy, Davey Schmitt, Sarah Zink. Pacific Northwest Research Institute.

Pennsylvania Satellite Center: Dorothy Becker, M.D., Margaret Franciscus, MaryEllen Dalmagro-Elias Smith, Ashi Daftary, M.D., Mary Beth Klein, Chrystal Yates. Children's Hospital of Pittsburgh of UPMC.

Data Coordinating Center: Jeffrey P. Krischer, Ph.D., PI, Rajesh Adusumali, Sarah Austin-Gonzalez, Maryouri Avendano, Sandra Baethke, Brant Burkhardt, Ph.D., Martha Butterworth, Nicholas Cadigan, Joanna Clasen, Kevin Counts, Christopher Eberhard, Steven Fiske, Laura Gandolfo, Jennifer Garmeson, Veena Gowda, Belinda Hsiao, Christina Karges, Qian Li, Ph.D., Shu Liu, Xiang Liu, Ph.D., Kristian Lynch, Ph.D., Jamie Malloy, Cristina McCarthy, Jose Moreno, Hemang M. Parikh, Ph.D., Cassandra Remedios, Chris Shaffer, Susan Smith, Noah Sulman, Ph.D., Roy Tamura, Ph.D., Dena Tewey, Michael Toth, Ulla Uusitalo, Ph.D., Kendra Vehik, Ph.D., Ponni Vijayakandipan, Melissa Wroble, Jimin Yang, Ph.D., R.D., Kenneth Young, Ph.D. Past staff: Michael Abbondondolo, Lori Ballard, Rasheedah Brown, David Cuthbertson, Stephen Dankyi, David Hadley, Ph.D., Kathleen Heyman, Francisco Perez Laras, Hye-Seung Lee, Ph.D., Colleen Maguire, Wendy McLeod, Aubrie Merrell, Steven Meulemans, Ryan Quigley, Laura Smith, Ph.D. University of South Florida.

Project Scientist: Beena Akolkar, Ph.D. National Institutes of Diabetes and Digestive and Kidney Diseases.

Other Contributors: Thomas Briese, Ph.D., Columbia University. Todd Brusko, Ph.D., University of Florida. Suzanne Bennett Johnson, Ph.D., Florida State University. Eoin McKinney, Ph.D., University of Cambridge. Tomi Pastinen, M.D., Ph.D., The Children's Mercy Hospital. Eric Triplett, Ph.D., University of Florida.

Dietary Biomarkers Laboratory: Iris Erlund, Ph.D., Irma Salminen, Jouko Sundvall, Nina Kangas, Petra Arohonka. Finnish Institute for Health and Welfare, Helsinki, Finland.

Funding

The TEDDY study is funded by U01 DK63829, U01 DK63861, U01 DK63821, U01 DK63865, U01 DK63863, U01 DK63836, U01 DK63790, UC4 DK63829, UC4 DK63861, UC4 DK63821, UC4 DK63865, UC4 DK63863, UC4 DK63836, UC4 DK95300, UC4 DK100238, UC4 DK106955, UC4 DK112243, UC4 DK117483, and contract no. HHSN267200700014C from the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), the National Institute of Allergy and Infectious Diseases (NIAID), the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the National Institute of Environmental Health Sciences (NIEHS), the Centers for Disease Control and Prevention (CDC), and Juvenile Diabetes Research Foundation (JDRF). This study was supported in part by the NIH/NCATS Clinical and Translational Science Awards to the University of Florida (UL1 TR000064) and the University of Colorado (UL1 TR002535). The sponsors were represented in The Environmental Determinants of Diabetes in the Young (TEDDY) Steering Committee. The sponsors had a role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; and review and approval to submit for publication. The sponsors did not have the right to veto submission to any journal. The corresponding author had the final say in submitting the manuscript for publication.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2021.720041/full#supplementary-material

Abbreviations

25(OH)D, 25-hydroxyvitamin D; CDA, celiac disease autoimmunity; FDR, first degree relative; HLA, human leukocyte antigen; IA, islet autoimmunity; LDP, long distance protocol; SNP, single nucleotide polymorphisms; TEDDY, the environmental determinants of diabetes in the young; tTGA, tissue transglutaminase autoantibodies; T1D, type 1 diabetes.

References

1. Sassi F, Tamone C, D'Amelio P. Vitamin D: nutrient, hormone, and immunomodulator. Nutrients. (2018) 10:1656. doi: 10.3390/nu10111656

CrossRef Full Text | Google Scholar

2. Misra M, Pacaud D, Petryk A, Collett-Solberg PF, Kappy M. Vitamin D deficiency in children and its management: review of current knowledge and recommendations. Pediatrics. (2008) 122:398–417. doi: 10.1542/peds.2007-1894

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Norris JM, Lee HS, Frederiksen B, Erlund I, Uusitalo U, Yang J, et al. Plasma 25-hydroxyvitamin D concentration and risk of islet autoimmunity. Diabetes. (2018) 67:146–54. doi: 10.2337/db17-0802

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Miettinen ME, Niinistö S, Erlund I, Cuthbertson D, Nucci AM, Honkanen J, et al. Serum 25-hydroxyvitamin D concentration in childhood and risk of islet autoimmunity and type 1 diabetes: the TRIGR nested case-control ancillary study. Diabetologia. (2020) 63:780–7. doi: 10.1007/s00125-019-05077-4

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hypponen E, Laara E, Reunanen A, Jarvelin MR, Virtanen SM. Intake of vitamin D and risk of type 1 diabetes: a birth-cohort study. Lancet. (2001) 358:1500–3. doi: 10.1016/S0140-6736(01)06580-1

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Munger KL, Levin LI, Hollis BW, Howard NS, Ascherio A. Serum 25-hydroxyvitamin D levels and risk of multiple sclerosis. JAMA. (2006) 296:2832–8. doi: 10.1001/jama.296.23.2832

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Lin J, Liu J, Davies ML, Chen W. Serum vitamin D level and rheumatoid arthritis disease activity: review and meta-analysis. PLoS ONE. (2016) 11:e0146351. doi: 10.1371/journal.pone.0146351

PubMed Abstract | CrossRef Full Text | Google Scholar

8. White JH. Vitamin D deficiency and the pathogenesis of Crohn's disease. J Steroid Biochem Mol Biol. (2018) 175:23–8. doi: 10.1016/j.jsbmb.2016.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Lebwohl B, Sanders DS, Green PHR. Coeliac disease. Lancet. (2018) 391:70–81. doi: 10.1016/S0140-6736(17)31796-8

PubMed Abstract | CrossRef Full Text | Google Scholar

10. López Casado MÁ, Lorite P, Ponce de León C, Palomeque T, Torres MI. Celiac Disease Autoimmunity. Arch Immunol Ther Exp. (2018) 66:423–30. doi: 10.1007/s00005-018-0520-z

CrossRef Full Text | Google Scholar

11. Liu E, Lee HS, Aronsson CA, Hagopian WA, Koletzko S, Rewers MJ, et al. Risk of pediatric celiac disease according to HLA haplotype and country. N Engl J Med. (2014) 371:42–9. doi: 10.1056/NEJMoa1313977

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Kemppainen KM, Lynch KF, Liu E, Lonnrot M, Simell V, Briese T, et al. Factors that increase risk of celiac disease autoimmunity after a gastrointestinal infection in early life. Clin Gastroenterol Hepatol. (2017) 15:694–702. e5. doi: 10.1016/j.cgh.2016.10.033

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Stene LC, Honeyman MC, Hoffenberg EJ, Haas JE, Sokol RJ, Emery L, et al. Rotavirus infection frequency and risk of celiac disease autoimmunity in early childhood: a longitudinal study. Am J Gastroenterol. (2006) 101:2333–40. doi: 10.1111/j.1572-0241.2006.00741.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Andren Aronsson C, Lee HS, Hard Af Segerstad EM, Uusitalo U, Yang J, Koletzko S, et al. Association of gluten intake during the first 5 years of life with incidence of celiac disease autoimmunity and celiac disease among children at increased risk. JAMA. (2019) 322:514–23. doi: 10.1001/jama.2019.10329

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lindfors K, Lin J, Lee HS, Hyoty H, Nykter M, Kurppa K, et al. Metagenomics of the faecal virome indicate a cumulative effect of enterovirus and gluten amount on the risk of coeliac disease autoimmunity in genetically at risk children: the TEDDY study. Gut. (2020) 69:1416–22. doi: 10.1136/gutjnl-2019-319809

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Tanpowpong P, Obuch JC, Jiang H, McCarty CE, Katz AJ, Leffler DA, et al. Multicenter study on season of birth and celiac disease: evidence for a new theoretical model of pathogenesis. J Pediatr. (2013) 162:501–4. doi: 10.1016/j.jpeds.2012.08.056

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Ivarsson A, Hernell O, Nystrom L, Persson LA. Children born in the summer have increased risk for coeliac disease. J Epidemiol Commun Health. (2003) 57:36–9. doi: 10.1136/jech.57.1.36

CrossRef Full Text | Google Scholar

18. Namatovu F, Lindkvist M, Olsson C, Ivarsson A, Sandström O. Season and region of birth as risk factors for coeliac disease a key to the aetiology? Arch Dis Child. (2016) 101:1114–8. doi: 10.1136/archdischild-2015-310122

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Lebwohl B, Green PH, Murray JA, Ludvigsson JF. Season of birth in a nationwide cohort of coeliac disease patients. Arch Dis Child. (2013) 98:48–51. doi: 10.1136/archdischild-2012-302360

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Moos C, Duus KS, Frederiksen P, Heitmann BL, Andersen V. Exposure to the danish mandatory vitamin D fortification policy in prenatal life and the risk of developing coeliac disease-the importance of season: a semi ecological study. Nutrients. (2020) 12:1243. doi: 10.3390/nu12051243

CrossRef Full Text | Google Scholar

21. Tanpowpong P, Camargo CA. Early-life vitamin D deficiency and childhood-onset coeliac disease. Public Health Nutr. (2014) 17:823–6. doi: 10.1017/S1368980013003510

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Vici G, Camilletti D, Polzonetti V. Possible role of vitamin D in celiac disease onset. Nutrients. (2020) 12:1051. doi: 10.3390/nu12041051

CrossRef Full Text | Google Scholar

23. Ciacci C, Bilancio G, Russo I, Iovino P, Cavallo P, Santonicola A, et al. 25-hydroxyvitamin D, 1,25-dihydroxyvitamin D, and peripheral bone densitometry in adults with celiac disease. Nutrients. (2020) 12:929. doi: 10.3390/nu12040929

CrossRef Full Text | Google Scholar

24. Björck S, Brundin C, Karlsson M, Agardh D. Reduced bone mineral density in children with screening-detected celiac disease. J Pediatr Gastroenterol Nutr. (2017) 65:526–32. doi: 10.1097/MPG.0000000000001568

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Corazza GR, Di Sario A, Cecchetti L, Tarozzi C, Corrao G, Bernardi M, et al. Bone mass and metabolism in patients with celiac disease. Gastroenterology. (1995) 109:122–8. doi: 10.1016/0016-5085(95)90276-7

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Lionetti E, Galeazzi T, Dominjanni V, Acquaviva I, Catassi GN, Iasevoli M, et al. Lower level of plasma 25-hydroxyvitamin d in children at diagnosis of celiac disease compared with healthy subjects: a case-control study. J Pediatr. (2021) 228:132–137. e1. doi: 10.1016/j.jpeds.2020.08.089

PubMed Abstract | CrossRef Full Text | Google Scholar

27. TEDDY Study Group. The environmental determinants of diabetes in the young (TEDDY) study: study design. Pediatr Diabetes. (2007) 8:286–98. doi: 10.1111/j.1399-5448.2007.00269.x

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Rewers M, Hyöty H, Lernmark Å, Hagopian W, She JX, Schatz D, et al. The environmental determinants of diabetes in the young (TEDDY) study: 2018 update. Curr Diabetes Rep. (2018) 18:136. doi: 10.1007/s11892-018-1113-2

CrossRef Full Text | Google Scholar

29. Lee HS, Burkhardt BR, McLeod W, Smith S, Eberhard C, Lynch K, et al. Biomarker discovery study design for type 1 diabetes in the environmental determinants of diabetes in the young (TEDDY) study. Diabetes Metab Res Rev. (2014) 30:424–34. doi: 10.1002/dmrr.2510

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Hossein-nezhad A, Holick MF. Vitamin D for health: a global perspective. Mayo Clinic Proc. (2013) 88:720–55. doi: 10.1016/j.mayocp.2013.05.011

CrossRef Full Text | Google Scholar

31. Munns CF, Shaw N, Kiely M, Specker BL, Thacher TD, Ozono K, et al. Global consensus recommendations on prevention and management of nutritional rickets. J Clin Endocrinol Metab. (2016) 101:394–415. doi: 10.1210/jc.2015-2175

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen WM. Robust relationship inference in genome-wide association studies. Bioinformatics. (2010) 26:2867–73. doi: 10.1093/bioinformatics/btq559

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Therneau T, Grambsch P. Modeling Survival Data: Extending the Cox Model. New York, NY: Springer (2000).

Google Scholar

34. Therneau T, Lumley T, Atkinson E, Crowson C. Survival: survival analysis. (2020). Available online at: https://cran.r-project.org/package=survival (accessed April 26, 2021).

Google Scholar

35. Yang CY, Leung PS, Adamopoulos IE, Gershwin ME. The implication of vitamin D and autoimmunity: a comprehensive review. Clin Rev Allergy Immunol. (2013) 45:217–26. doi: 10.1007/s12016-013-8361-3

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Mårild K, Tapia G, Haugen M, Dahl SR, Cohen AS, Lundqvist M, et al. Maternal and neonatal vitamin D status, genotype and childhood celiac disease. PLoS ONE. (2017) 12:e0179080. doi: 10.1371/journal.pone.0179080

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Bishop E, Ismailova A, Dimeloe SK, Hewison M, White JH. Vitamin D and immune regulation: antibacterial, antiviral, anti-inflammatory. JBMR Plus. (2020) 5:e10405. doi: 10.1002/jbm4.10405

CrossRef Full Text | Google Scholar

38. McGrath JJ, Saha S, Burne TH, Eyles DW. A systematic review of the association between common single nucleotide polymorphisms and 25-hydroxyvitamin D concentrations. J Steroid Biochem Mol Biol. (2010) 121:471–7. doi: 10.1016/j.jsbmb.2010.03.073

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Grant WB, Karras SN, Bischoff-Ferrari HA, Annweiler C, Boucher BJ, Juzeniene A, et al. Do studies reporting 'U'-shaped serum 25-hydroxyvitamin D-health outcome relationships reflect adverse effects? Dermatoendocrinol. (2016) 8:e1187349. doi: 10.1080/19381980.2016.1187349

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Sepulveda-Villegas M, Elizondo-Montemayor L, Trevino V. Identification and analysis of 35 genes associated with vitamin D deficiency: a systematic review to identify genetic variants. J Steroid Biochem Mol Biol. (2020) 196:105516. doi: 10.1016/j.jsbmb.2019.105516

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Capriati T, Francavilla R, Castellaneta S, Ferretti F, Diamanti A. Impact of the birth's season on the development of celiac disease in Italy. Eur J Pediatr. (2015) 174:1657–63. doi: 10.1007/s00431-015-2589-2

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Aranow C. Vitamin D and the immune system. J Investig Med. (2011) 59:881–6. doi: 10.2310/JIM.0b013e31821b8755

CrossRef Full Text | Google Scholar

43. Bittker SS, Bell KR. Potential risk factors for celiac disease in childhood: a case-control epidemiological survey. Clin Exp Gastroenterol. (2019) 12:303–19. doi: 10.2147/CEG.S210060

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Bittker SS. Elevated levels of 1,25-dihydroxyvitamin D in plasma as a missing risk factor for celiac disease. Clin Exp Gastroenterol. (2020) 13:1–15. doi: 10.2147/CEG.S222353

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Yang J, Tamura RN, Uusitalo UM, Aronsson CA, Silvis K, Riikonen A, et al. Vitamin D and probiotics supplement use in young children with genetic risk for type 1 diabetes. Eur J Clin Nutr. (2017) 71:1449–54. doi: 10.1038/ejcn.2017.140

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Smith LB, Lynch KF, Baxter J, Lernmark B, Roth R, Simell T, et al. Factors associated with maternal-reported actions to prevent type 1 diabetes in the first year of the TEDDY study. Diabetes Care. (2014) 37:325–31. doi: 10.2337/dc13-0449

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Michos ED, Kalyani RR, Segal JB. Why USPSTF Still Finds Insufficient Evidence to Support Screening for Vitamin D Deficiency. JAMA Netw Open. (2021) 4:e213627. doi: 10.1001/jamanetworkopen.2021.3627

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Antonucci R, Locci C, Clemente MG, Chicconi E, Antonucci L. Vitamin D deficiency in childhood: old lessons and current challenges. J Pediatr Endocrinol Metab. (2018) 31:247–60. doi: 10.1515/jpem-2017-0391

PubMed Abstract | CrossRef Full Text | Google Scholar