Gut microbiome markers in subgroups of HLA class II genotyped infants signal future celiac disease in the general population: ABIS study
Patricia L. Milletich1 
Angelica P. Ahrens1 
Jordan T. Russell1 
Joseph R. Petrone1 Meghan A. Berryman1 
Daniel Agardh2 Jonas F. Ludvigsson3,4 
Eric W. Triplett1* 
Johnny Ludvigsson5- 1Department of Microbiology and Cell Science, Institute of Food and Agricultural Sciences, University of Florida, Gainesville, FL, United States
- 2Department of Clinical Sciences, Lund University, Malmö, Sweden
- 3Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, Sweden
- 4Department of Pediatrics, Örebro University Hospital, Örebro, Sweden
- 5Crown Princess Victoria’s Children’s Hospital and Division of Pediatrics, Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
Although gut microbiome dysbiosis has been illustrated in celiac disease (CD), there are disagreements about what constitutes these microbial signatures and the timeline by which they precede diagnosis is largely unknown. The study of high-genetic-risk patients or those already with CD limits our knowledge of dysbiosis that may occur early in life in a generalized population. To explore early gut microbial imbalances correlated with future celiac disease (fCD), we analyzed the stool of 1478 infants aged one year, 26 of whom later acquired CD, with a mean age of diagnosis of 10.96 ± 5.6 years. With a novel iterative control-matching algorithm using the prospective general population cohort, All Babies In Southeast Sweden, we found nine core microbes with prevalence differences and seven differentially abundant bacteria between fCD infants and controls. The differences were validated using 100 separate, iterative permutations of matched controls, which suggests the bacterial signatures are significant in fCD even when accounting for the inherent variability in a general population. This work is the first to our knowledge to demonstrate that gut microbial differences in prevalence and abundance exist in infants aged one year up to 19 years before a diagnosis of CD in a general population.
Introduction
Celiac disease (CD) is an antigen-driven, immune-mediated disorder occurring in 1% of the Western population and increasing in incidence (Lebwohl et al., 2015; Sharma et al., 2019; King et al., 2020). The presence of gluten drives intestinal inflammation and the production of autoantibodies against tissue transglutaminase (tTG) in genetically predisposed individuals (Sollid and Lie, 2005).
Gliadin, a prolamin found in gluten, is rich in glutamine and proline, making it partially resistant to gastrointestinal digestion (Shan et al., 2002). When glutamine-rich gluten peptides reach the lumen, they interact with tTG, deaminating glutamine into glutamate. This deamination causes the peptide to become negatively charged and preferentially binds to the class II human leukocyte antigen (HLA) molecules DQ2 and DQ8 for the presentation to CD4+ T cells (Sollid, 2017). HLA proteins encode the human equivalent of the major histocompatibility complex (MHC) proteins, which display antigens for T-helper cell recognition to elicit tolerance or immune response. Protein structure altering allelic polymorphisms found in DQ2.5, DQ2.2, or DQ8 form regions P4 and P6 for DQ2.5 and DQ2.2, and P1 and P9 for DQ8 that attract negatively charged residues of glutamate making an individual more susceptible to celiac (Moustakas et al., 2000; Tollefsen et al., 2006). HLA is strongly associated with CD, with DQ2/DQ8 found in up to 98% of CD patients (Cecilio and Bonatto, 2015).
It is estimated that approximately 30% of the world population carries either DQ2 or DQ8. Yet only 1% of those with the risk DQs haplotypes are diagnosed with CD, pointing toward other factors that increase the risk for or protect from the disease (Mustalahti et al., 2010; Megiorni and Pizzuti, 2012; Jamnik et al., 2017). Research suggests that intestinal dysbiosis is characteristic of CD etiology, with aberrant microbiome signatures presenting even in individuals who have followed a gluten-free diet (Collado et al., 2009; Wu et al., 2021).
Prospective studies of gut dysbiosis in CD are limited, especially in the years preceding a diagnosis. Environmental and genetic mediators of the infant gut microbiome have been previously described, including mode of delivery, infant feeding (e.g., duration of breastfeeding), gestational age, antibiotic (Penders et al., 2006), and, recently, HLA genetics (Russell et al., 2019). The influence of HLA-DQ and environmental factors on gut colonization continues to be studied. However, these studies focus primarily on infants already at high risk for CD (Palma et al., 2012; Olivares et al., 2015). Restricting the study of gut microbiota to infants with a first-degree relative with CD may lead to conclusions that are not generalizable to the larger population, and cross-sectional study designs limit the identification of early microbial risk marks for future celiac disease (fCD).
The current investigation was derived from the general population cohort, All Babies In Southeast Sweden (ABIS), which sought to identify early risk and protective factors for future disease. Here, ABIS allowed for the study of the contribution of HLA genetics and environmental factors to increased CD risk and for the discovery of microbial differences in infants at one year of age who developed CD up to two decades later.
Materials and methods
Study population
This investigation derives from the ABIS prospective, population-based cohort of 17,055 children born in southeast Sweden from October 1, 1997, to October 1, 1999. The ABIS cohort was designed to identify environmental and genetic factors associated with autoimmune diseases (Ludvigsson et al., 2001). After being provided oral, written, or video information, parents gave their informed consent for their child’s participation. Mothers completed questionnaires and diaries for the first year of life, which included but were not limited to information on gestational age, mode of delivery, antibiotic use, age of introduction or frequency of consumption of certain foods (fruit, vegetables, eggs, meat, gluten, among others), duration of breastfeeding, and environmental factors and toxic exposures. Infants were binned into geographical regions based on their municipality at one year of age to help account for cultural differences possibly contributing to the microbiome structure (Figure S1).
Sample collection
Stool samples were collected from infants’ diapers using provided sterile spatula and tubes. These samples were immediately frozen after collection. Samples collected at home or at the WellBaby Clinic were sent to the Research Laboratory of Division of Pediatrics, Medical Faculty, Linkoping University, where all samples were stored at -80°C. Samples for this study were shipped to the University of Florida in a frozen state and stored at -80°C until extraction (Russell et al., 2019).
Institutional review board approval
The ABIS study was approved by the Research Ethics Committees of the Faculty of Health Science at Linköping University, Sweden, 1997/96287 and 2003/03-092 and the Medical Faculty of Lund University, Sweden (Dnr 99227, Dnr 99321) as described previously (Russell et al., 2019). The microbiome analysis performed at the University of Florida was approved by the University of Florida’s Institutional Review Board as an exempt study IRB201800903.
Diagnoses from the swedish national patient registry
The Swedish Patient Register provided data on autoimmune and neurological disorders in ABIS children, with the most recent update having occurred in 2021 (international classification of disease, ICD-10: K90) (Ludvigsson et al., 2011). For a diagnosis of CD for this investigation, at least two records of CD in the Patient Register were necessary. Up until 2012, Swedish pediatricians required both positive serology and a positive small intestinal biopsy for diagnosis. A validation study found that 67/67 (100%; 95%CI=95–100%) of Swedish pediatric gastroenterologists performed a small intestinal biopsy in >90% of all children with suspected CD prior to diagnosis (Ludvigsson et al., 2011). That same validation study found that in a Swedish setting, 95% of individuals with Marsh III (villous atrophy) had CD (Ludvigsson et al., 2011). In 2012, ESPGHAN offered a non-biopsy pathway in selected children with levels of tTG antibodies ten times the upper limit of normal, and these recommendations were re-iterated by the Swedish working group for pediatric CD (Husby et al., 2012).
HLA genotyping
HLA genotyping was performed with the available blood spots of the ABIS children. HLA-DR-DQ genotypes are defined by the presence of common European HLA-DR-DQ alleles associated with disease risk using sequence-specific hybridization with lanthanide-labeled oligonucleotide probes (Ilonen et al., 2016).
The CD-associated HLA genotypes were categorized as a combination of known risk alleles, HLA-DQA1*05:01-DQB1*02 (DQ2.5), DQA1*03:01-DQB1*03:02 (DQ8), DQA1*02:01-DQB1*02 (DQ2.2), possibly protective HLA-DRB1*01:01/DQA1*01:01-DQB1*05:01 (DR1-DQ5), and CD non-associated HLA-DQs (DQX).
Generation of amplicons for microbiome analysis
Microbial composition analysis was run in stool samples collected at one year of age. DNA extraction, 16S rRNA barcoded PCR, and V3-V4 16S rRNA Illumina sequencing were performed as previously described (Russell et al., 2019). Bacterial quantification through universal 16s rRNA primers was performed as described previously (Russell et al., 2021).
Forward and Reverse fastq reads were merged using default PEAR parameters, and ambiguous primers were removed using Cutadapt using the following parameters: Forward, NNNNCCTACGGGAGGCAGCAG -e 9. Reverse, ATTAGATACCCSBGTAGTCCCC -e 11 (Zhang et al., 2014; Martin, 2011). ASVs were left untruncated and set to a minimum length of 400 nucleotides before filtering using the DADA2 R package (Callahan et al., 2016; R Core Team, 2016). Sequences were filtered to allow no ambiguous nucleotides, a maximum expected error rate of 2, Q-score 10, and PhiX reads removed. Sequences were then “pseudo-pooled” before the chimeras were removed with the default consensus method. Taxonomy was assigned using the SILVA v138 database (Quast et al., 2013).
Population subset and filtering
Infants with fewer than 1000 raw reads were not included in this analysis, nor were those later diagnosed with a neurological disorder or an autoimmune disease other than CD, resulting in 1452 controls and 26 fCD. Other excluded autoimmune diseases include hypothyroidism, type-1 diabetes, juvenile idiopathic arthritis, Crohn’s disease, or ulcerative colitis.
The average read counts per sample were 46627, with a maximum of 632001. In this analysis, 10389 ASVs were filtered to only include those present in at least five infants at five or more reads, resulting in 2376 ASVs that were then binned to 184 unique genera. Relative abundances were calculated using the transform_sample_counts function, and ASVs belonging to the same genus were combined using tax_glom (McMurdie and Holmes, 2013). To calculate total abundance, the relative abundance values were multiplied by sample by the total number of copies of 16s rRNA per gram of stool, as determined through qPCR. Significance p-values in the microbiome analysis were adjusted for false discovery rate (FDR). The statistical program R Foundation for Statistical Computing was used for all computational analysis (R Core Team, 2018).
Early life factors associated with fCD
To determine factors associated with fCD, chi-square tests were conducted to compare the year-one diary responses (n=1478) using Monte Carlo simulation.
Confounders of beta microbial diversity
Confounders of beta diversity were assessed within the total population (n = 1478) at the genus level. Subjects with missing data for each variable were excluded from the calculation. The distance function from the phyloseq package was applied using the binomial method. The distance matrix was tested against the variable of interest using default parameters with the adonis2 R package (Oksanen et al., 2020).
Random permutations for matched case-control iterations
To compensate for the 55-fold increase of controls to fCD infants, 100 random selections of controls were employed to create 100 separate randomized matched controls for our comparisons (Figure 1). The iterative process allowed us to account for the inherent variability within the controls in a general population cohort without the risk of overpowering the results with an imbalance of controls. For each iteration, controls were selected for each future celiac disease infant (fCD) to balance geography (region of Sweden) and the presence of siblings at birth in the comparative groups. These variables were selected because they were identified as confounders of beta diversity in this cohort (ANOVA: padj < 0.05) (Table S1). From the resulting total subset of controls, two infants were randomly selected as matches for each fCD infant. This process was performed 100 times. The set.seed() function was used for each iteration to allow for replication of case-control permutation groups.
Figure 1 Pipeline for the matching iterations algorithm and analysis of future celiac disease (fCD). The differential prevalence (The Prevalence Interval for Microbiome Evaluation, PIME) and abundance (DESeq2) analyses were conducted independently at each for 100 iterations of matched controls with the fCD cases. Matches were randomly selected from subsets based on microbiome confounders of geographic region and siblings at birth. The random selection of controls was applied to investigate bacterial differences in cases and controls while accounting for inherent variability within controls in the general population. Genera that were of most interest to this investigation were defined as those that were shared across at least fifty iterations of randomized case-control groups.
Differentially abundant genera in fCD versus controls
Differential abundances of microbes without a prevalence threshold were assessed with the binomial distribution model R package, DESeq2 (Love et al., 2014). For both the total abundance and relative abundance analyses, the estimateSizeFactor() function was used with “poscount” type, allowing for zeroes. Within DESeq2, the Wald test for calculating p values, with a local fit type, was employed. The additional Cook’s distance filtering was not considered. After each iteration of DESeq2 (within each of the 100 case-control permutations), significance values were adjusted for multiple comparisons across the genera. Genera were deemed significant to this investigation if they appeared in at least 50 of the 100 iterations of case-control matching with a padj < 0.05.
Prevalence filtering of the core microbiome in fCD
Prevalence is a metric that is often overlooked in microbiome studies. Prevalence is defined as a proportion of samples in a specific group that share taxa, irrespective of abundance. The Prevalence Interval for Microbiome Evaluation (PIME) was employed at each of the 100 case-control group comparisons to obtain the core genera representing fCD or controls and to determine whether any bacteria were more likely to be present or absent in fCD (Roesch et al., 2020). The PIME package was designed as a robust filtering methodology for low prevalence, and its reliability has been validated against existing methods with several 16S rRNA independent datasets. PIME determines the set of bacteria that meet the prevalence criteria at increasing increments of 5%, and the optimal prevalence threshold is selected by the user, balancing out of bag error (OOB) rates while maximizing the representation of the bacteria (i.e., number of taxa retained after filtering).
The first ten iterations of case-control comparisons were used as guides for the expected OOB in the remaining ninety iterations. The initial OOB for the first ten runs ranged from 33.3% to 42.3%, with an average of 38.0%. An optimal prevalence threshold of 60% was selected for these analyses, as an average OOB error of 2.82 was observed at this interval, with an average of 35 of the original 184 genera (19%) and 5.02E9 of the original sequences retained. After filtering for prevalence, binomial distances for iterations 25, 50, 75, and 100 were conducted using plot_ordination and stat_ellipse (Wickham, 2016).
To determine which genera consistently differentiated controls and fCD using this approach, the top 30 taxa with the highest mean decrease accuracy (MDA) were considered at each iteration. Taxa were removed with negative MDA values, and the mean abundance for the remaining taxa were determined and binned for all iterations. Taxa were removed if they were significant in fewer than 50 of the 100 iterations. The direction of the effect was assessed for each remaining genus with respect to the average abundance in fCD and controls.
Conditional factors
Kruskal-Wallis or Mann-Whitney U tests were used to determine the effect of environmental and genetic factors against the total abundance and relative abundance of the significant genera identified in the DESeq2 or PIME analyses. These tests were conducted at the level of the complete cohort (n=1478). P-values were adjusted using false-discovery-rate. Factors that were significant in both total abundance and relative abundance are presented in Table S3.
Results
Diagnoses
Twenty-six subjects were diagnosed with CD according to the Swedish National Patient Registry (Walker-Smith and Sanderson, 1999). Female subjects (n=17, 65%) had a mean age of diagnosis of 10.71 ± 5.17 years, while the male age of diagnosis was 11.44 ± 6.65 years (Table 1).
Table 1 Distribution of genetic and environmental factors for infants with future CD.
Core microbiome signatures: Bacteria with differential prevalence in fCD
With a prevalence threshold of 60% (Table S2), distinct clustering of fCD infants (n=26) and controls (n=52) was observed (Figure 2).
Figure 2 Microbiome diversity differences in selected case-control iterations for future celiac disease (fCD). Principal component analysis (PCoA) of the binomial distance of genera between fCD (n=26) and controls (n=52) is presented for Iteration #25, Iteration #50, Iteration #75, and Iteration #100, depicting the unique clustering of the core microbiome after applying a prevalence filter at 60% in PIME. Data are presented for (A) reads/g (total abundance) and (B) relative abundance. The percentages of the variance explained by each of the two components are presented in the axes.
Five genera were higher in fCD infants than controls, including three that were shared in both read/g (total abundance) and relative abundance analyses: Erysipelatoclostridium, Haemophilus, and Lachnospiraceae NK4A136 group. Lachnoclostridium only appeared in the reads/g analysis, while Lachnospira was only present in the relative abundance analysis (Figures 3A, B).
Figure 3 Distribution of the total (reads/g) and relative (%) abundance of genera differentially abundant in future celiac Disease (fCD) or controls. A, C) relative abundance B, D) core bacteria split by trends of abundance. (A, B) Bacteria greater in fCD infants. (C, D) Bacteria greater in controls.
Alternatively, nine genera were found higher in controls than fCD infants, six of which were found in both types of abundance analysis, including Butyriococcus, Collinsella, Monoglobus, Roseburia, Ruminococcus, and Terrisporobacter (Figures 3C, D). Three other genera, Bacteroides, Gemella, and Veillonella, were only present in the reads/g analysis.
Differentially abundant microbes in fCD and controls
Across the iterations of matched case-control, 25 genera were consistently in higher abundance in controls (DESeq: p-adj ≤ 0.05), seven of which were found in both reads/g (total abundance) and relative abundance analyses. These include Anaeroglobus, Barnesiella, Candidatus Soleaferrea, Eubacterium, Monoglobus, Senegalimassilia, and Oscillospiraceae UCG.002 (Figure 4). Candidatus Soleaferrea was found to be completely absent in fCD infants but was present in an average of 21.7% of controls. The remaining 18 genera were found only in the reads/g analysis, five of which were also completely absent in controls (Figure S2).
Figure 4 Differential abundance of genera in future celiac disease (fCD). DESeq Log2FoldChange (LFC) statistics of genera found in at least half of the matching iterations of cases and controls in (A) reads/g (total abundance) and (B) relative abundance. Genera with LFC < 0 are more abundant in controls than in fCD.
Genetic and environmental impact on bacteria
To determine potential causes of the taxonomic differences between fCD and controls, environmental and genetic factors were tested against the abovementioned bacteria that were shared in both relative abundance and read/g (total abundance) of the core microbiome and differential expression analyses (Table 2 and Table S3).
Table 2 Factors associated with prevalence or abundance of gut bacteria enriched or depleted in fCD.
Of the significant core microbes, 66% (n=6) were correlated to the presence or absence of HLA haplotypes; however, only 29% (n=2) of the bacteria identified in the differential expression analysis were associated with HLA. Monoglobus, found in both analyses, was impacted by both DR13-DQ603 and DR3-DQ2.5. Butyricicoccus and Collinsella were also correlated with DR13-DQ603. DR15-DQ602 was associated with Butyricicoccus as well as Lachnospiraceae NK4A136 group. Erysipelatoclostridium was correlated with DR14-DQ5, whereas Terrisporobacter was specifically correlated with DR14-DQ503. Barnesiella, the only non-core microbe impacted by HLA, was associated with DR7-DQ2.2 (Figures S3, S4).
Dietary factors (such as duration of total or exclusive breastfeeding and weekly meals involving eggs or beef) are associated with one (11%) of the core microbes, Collinsella, yet four (57%) of the microbes from differential expression analysis, Anaeroglobus, Eubacterium, Senegalimassilia, and UCG-002. Gastroenteritis in the first year of life was associated with three core microbes, Lachnospiraceae NK4A136 group, Roseburia, and Ruminococcus, and none of the microbes identified in the differential abundance analysis. Biological sex was associated with the greatest number of genera, six total, in both the core microbiome analysis (Lachnospiraceae NK4A136 group and Ruminococcus) and the differential abundance analysis (Candidatus Soleaferrea, Eubacterium, and Senegalimassilia), as well as Monoglobus which was significant in both analyses.
DR3-DQ2.5 was present in a higher proportion of fCD infants compared to all controls, 50% vs. 22% (padj < 0.05), but no other significant factors were found for the selected bacteria (Table S4).
Discussion
This research illustrates through the gut microbiome analysis of 1478 one-year-old infants that microbiome differences exist, both in prevalence and abundance, in the infant’s gut up to 19 years prior to celiac disease diagnosis. The average age of diagnosis of the 26 infants with future CD was 10.96 years, ranging between 2 and 20 years old. Despite the considerable time span between the age of diagnosis and stool collection, we observed microbial differences between healthy controls and fCD infants. Prevalence is a metric that is often overlooked in microbiome studies. In this investigation, we considered both differential abundance of the bacteria, as is customary, and prevalence of the bacteria, which is calculated irrespective of abundance. Whether the bacterium is present or absent may be physiologically meaningful, as may be the architecture of the “core” microbiome.
Core microbiomes were determined using genera prevalent in 60% of fCD cases or controls (n=52). The prevalence threshold was repeated with 100 unique case-control pairings to account for inherent differences in the distribution of controls. We defined genera of interest as those that were significant in 50 or more iterations, which provides further confidence that these core microbiome differences persist, despite the inherent variability of the control constituents.
Three core microbes were higher in fCD infants, Erysipelatoclostridium, Haemophilus, and Lachnospiraceae NK4A136 group. Erysipelatoclostridium and Haemophilus are opportunistic pathogens correlated with disease (Shao et al., 2017; Aranaz et al., 2021; Barberio et al., 2022; Cao et al., 2021; Liu et al., 2021; Milosavljevic et al., 2021). Lachnospiraceae NK4A136 group is a butyrate producer correlated with both health and disease (Ma et al., 2020; Stadlbauer et al., 2020; Park et al., 2022). Since Lachnospiraceae NK4A136 group is registered as an uncultured bacterium in the SILVA database, more research is needed to determine the true characteristics of this genus.
Monoglobus, unique for being the only genus significant in both core and differentially abundant microbes, is understudied but could be a valuable marker for health against CD. M. pectinilyticus, a pectin degrader, is the only current species described (Kim et al., 2019). Further research is needed to determine strain-specific consequences of this genus on host gut health.
Of the differentially abundant microbes, Candidatus Soleaferrea stands out for being absent in fCD infants but present in approximately 20% of the control cohorts. There is evidence that Candidatus Soleaferrea, which was absent in fCD infants, has anti-inflammatory properties, but this genus has also been correlated with several autoimmune diseases including CD (Cai et al., 2020; Cao et al., 2021; Zou et al., 2022). Within each genus, species are diverse and have a range of metabolic features, which can have various consequences on host physiology. Metagenomic sequencing and culture methods are necessary to understand the role of these specific microbes in CD etiology.
Most CD gut microbiome studies focus on high-genetic-risk patients or are conducted with subjects who have already received a CD diagnosis (Verdu and Schuppan, 2021). The focus on high-risk patients may obscure the microbial differences since HLA has been shown to influence the gut microbiome diversity (Russell et al., 2019). Here, we benefit from a prospective general-population cohort that affords a diverse distribution of HLA genetic features.
HLA polymorphisms were associated with the core microbiome composition of ABIS infants. Six of the nine genera (66%) were correlated with the presence or absence of HLA haplotypes: DR13-DQ603, DR14-DQ5[03], DR15-DQ602, DR3-DQ2.5. Interestingly, the correlative HLA are not all associated with risk for or protection from CD. Through differential analysis, Monoglobus and Barnesiella were the only two bacteria with HLA-associated abundance differences. A higher abundance of Barnesiella correlated with the absence of HLA-DR7-DQ2.2, a known risk allele for CD. The molecular mimicry theory suggests that similarities between microorganism antigens resemble self-antigens or cross-react with HLA molecules, which may help explain the mechanism behind HLA’s impact on bacterial abundance (Ogasawara et al., 1986). Core microbes could be impacted by this molecular mimicry mechanism, while differentially abundant bacteria in fCD appeared to be more associated with external environmental factors rather than HLA.
This study was limited by the number of infants who later developed CD due to the nature of the general population study. The average rate of CD in Sweden is approximately 2–3%, with regional differences (Myléus et al., 2009). In the current study, the CD diagnosis incidence is 1.8%. To account for the imbalance of fCD cases and controls, we use a novel algorithm that utilizes iterations of randomized matching to create 100 unique case-control cohorts, then used the genera present in at least half of these iterations.
Altogether, these data illustrate the need to consider a broader range of HLA haplotypes in investigating the role of the gut microbiome in CD, as both genetic and environmental factors bear consequence on prevalence and differential abundance of gut microbes. Without looking into non-risk HLA genotypes in CD studies, significant differences in the gut microbiome could be obscured based on the distribution of HLA genetics in the cohort.
Data availability statement
The forward and reverse 16S raw sequencing data generated in this study is available through the NCBI Sequence Read Archive under BioProject PRJNA854900.HLA genotypes and associated sample metadata used for statistical comparison are available in the source data file. R codes used for statistical comparisons and figure generation are available at: https://github.com/PMilletich/Celiac_Disease_Pipeline.
Author contributions
TM, AA, JT, and JRT performed data analysis. JT performed the laboratory analysis. TM, MRM, JonL, DA, EWT, and JohL did the writing. JohL developed the ABI cohort. ET and JohL conceived of the project. All authors contributed to the article and approved the submitted version.
Funding
ABIS was supported by Barndiabetesfonden (Swedish Child Diabetes Foundation); Swedish Council for Working Life and Social Research, Grant/Award Numbers: FAS2004-1775, FAS2004–1775; Swedish Research Council, Grant/Award Numbers: K2005-72 × -11242-11A and K2008- 69 × -20826-01-4, K2008-69 × -20826-01-4; Östgöta Brandstodsbolag; Medical Research Council of Southeast Sweden (FORSS); JDRF Wallenberg Foundation, Grant/Award Number: K 98-99D-12813-01A; ALF-and LfoU grants from Region Östergötland and Linköping University, Sweden; Joanna Cocozza Foundation.
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/fcimb.2022.920735/full#supplementary-material
References
Aranaz, P., Ramos-Lopez, O., Cuevas-Sierra, A., Martinez, J. A., Milagro, F. I., Riezu-Boj, J. I. (2021). A predictive regression model of the obesity-related inflammatory status based on gut microbiota composition. Int. J. Obes. (Lond). 45 (10), 2261–2268. doi: 10.1038/s41366-021-00904-4
Barberio, B., Facchin, S., Patuzzi, I., Ford, A. C., Massimi, D., Valle, G., et al. (2022). A specific microbiota signature is associated to various degrees of ulcerative colitis as assessed by a machine learning approach. Gut Microbes 14 (1), 2028366. doi: 10.1080/19490976.2022.2028366
Cai, J., Zhou, L., Song, X., Yin, M., Liang, G., Xu, H., et al. (2020). Alteration of intestinal microbiota in 3-deoxyglucosone-induced prediabetic rats. BioMed Res. Int. 2020, 8406846. doi: 10.1155/2020/8406846
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J., Holmes, S. P. (2016). DADA2: High-resolution sample inference from illumina amplicon data. Nat. Methods 13 (7), 581–583. doi: 10.1038/nmeth.3869
Cao, R. R., He, P., Lei, S. F. (2021). Novel microbiota-related gene set enrichment analysis identified osteoporosis associated gut microbiota from autoimmune diseases. J bone miner. Metab. 39 (6), 984–996. doi: 10.1007/s00774-021-01247-w
Cecilio, L. A., Bonatto, M. W. (2015). The prevalence of HLA DQ2 and DQ8 in patients with celiac disease, in family and in general population. Arquivos brasileiros cirurgia digestiva: ABCD = Braz. Arch. Digestive Surg. 28 (3), 183–185. doi: 10.1590/S0102-67202015000300009
Collado, M. C., Donat, E., Ribes-Koninckx, C., Calabuig, M., Sanz, Y. (2009). Specific duodenal and faecal bacterial groups associated with paediatric coeliac disease. J. Clin. Pathol. 62 (3), 264–269. doi: 10.1136/jcp.2008.061366
Husby, S., Koletzko, S., Korponay-Szabó, I. R., Mearin, M. L., Phillips, A., Shamir, R., et al. (2012). European Society for pediatric gastroenterology, hepatology, and nutrition guidelines for the diagnosis of coeliac disease. J. Pediatr. Gastroenterol. Nutr. 54 (1), 136–160. doi: 10.1097/MPG.0b013e31821a23d0
Ilonen, J., Kiviniemi, M., Lempainen, J., Simell, O., Toppari, J., Veijola, R., et al. (2016). Genetic susceptibility to type 1 diabetes in childhood - estimation of HLA class II associated disease risk and class II effect in various phases of islet autoimmunity. Pediatr. Diabetes 17 Suppl 22, 8–16. doi: 10.1111/pedi.12327
Jamnik, J., Villa, C. R., Dhir, S. B., Jenkins, D. J. A., El-Sohemy, A. (2017). Prevalence of positive coeliac disease serology and HLA risk genotypes in a multiethnic population of adults in Canada: a cross-sectional study. BMJ Open 7, e017678. doi: 10.1136/bmjopen-2017-017678
Kassambara, A. (2020) Ggpubr: 'ggplot2' based publication ready plots. r package version 0.4.0. Available at: https://CRAN.R-project.org/package=ggpubr.
Kim, C. C., Healey, G. R., Kelly, W. J., Patchett, M. L., Jordens, Z., Tannock, G. W., et al. (2019). Genomic insights from monoglobus pectinilyticus: a pectin-degrading specialist bacterium in the human colon. ISME J. 13 (6), 1437–1456. doi: 10.1038/s41396-019-0363-6
King, J. A., Jeong, J., Underwood, F. E., Quan, J., Panaccione, N., Windsor, J. W., et al. (2020). Incidence of celiac disease is increasing over time: A systematic review and meta-analysis. Am. J. Gastroenterol. 115 (4), 507–525. doi: 10.14309/ajg.0000000000000523
Lebwohl, B., Ludvigsson, J. F., Green, P. H. (2015). Celiac disease and non-celiac gluten sensitivity. BMJ (Clinical Res. ed.) 351, h4347. doi: 10.1136/bmj.h4347
Liu, Y., Li, W., Yang, H., et al. (2021). Leveraging 16S rRNA microbiome sequencing data to identify bacterial signatures for irritable bowel syndrome. Front. Cell Infect. Microbiol. 11. doi: 10.3389/fcimb.2021.645951
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550. doi: 10.1186/s13059-014-0550-8
Ludvigsson, J. F., Andersson, E., Ekbom, A., Feychting, M., Kim, J. L., Reuterwall, C., et al. (2011). External review and validation of the Swedish national inpatient register. BMC Public Health 11, 450. doi: 10.1186/1471-2458-11-450
Ludvigsson, J., Ludvigsson, M. a .d, Sepa, A. (2001). Screening for prediabetes in the general child population: maternal attitude to participation. Pediatr. Diabetes 2 (4), 170–174. doi: 10.1034/j.1399-5448.2001.20405.x
Ma, L., Ni, Y., Wang, Z., Tu, W., Ni, L., Zhuge, F., et al. (2020). Spermidine improves gut barrier integrity and gut microbiota function in diet-induced obese mice. Gut Microbes 12 (1), 1–19. doi: 10.1080/19490976.2020.1832857
Martin, M.. (2021). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J 17 (1), 10–12. doi: 10.14806/ej.17.1.200
McMurdie, P. J., Holmes, S. (2013). Phyloseq: an r package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8 (4), e61217. doi: 10.1371/journal.pone.0061217
Megiorni, F., Pizzuti, A. (2012). HLA-DQA1 and HLA-DQB1 in celiac disease predisposition: practical implications of the HLA molecular typing. J. Biomed. Sci. 19 (1), 88. doi: 10.1186/1423-0127-19-88
Milosavljevic, M. N., Kostic, M., Milovanovic, J., Zaric, R. Z., Stojadinovic, M., Jankovic, S. M., et al. (2021). Antimicrobial treatment of erysipelatoclostridium ramosum invasive infections: a systematic review. Rev. Inst Med. Trop. Sao Paulo. 63, e30. doi: 10.1590/S1678-9946202163030
Moustakas, A. K., van de Wal, Y., Routsias, J., Kooy, Y. M., van Veelen, P., Drijfhout, J. W., et al. (2000). Structure of celiac disease-associated HLA-DQ8 and non-associated HLA-DQ9 alleles in complex with two disease-specific epitopes. Int. Immunol. 12 (8), 1157–1166. doi: 10.1093/intimm/12.8.1157
Mustalahti, K., Catassi, C., Reunanen, A., Fabiani, E., Heier, M., McMillan, S., et al. (2010). The prevalence of celiac disease in Europe: results of a centralized, international mass screening project. Ann. Med. 42 (8), 587–595. doi: 10.3109/07853890.2010.505931
Myléus, A., Ivarsson, A., Webb, C., Danielsson, L., Hernell, O., Högberg, E., et al. (2009). Celiac disease revealed in 3% of Swedish 12-year-olds born during an epidemic. journal of pediatric. Gastroenterol. Nutr. 49 (2), 170–176. doi: 10.1097/MPG.0b013e31818c52cc
Ogasawara, M., Kono, D. H., Yu, D. T. (1986). Mimicry of human histocompatibility HLA-B27 antigens by klebsiella pneumoniae. Infect Immun. 51 (3), 901–908. doi: 10.1128/iai.51.3.901-908.1986
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2020). Vegan: Community ecology package. R package version 2.5-7. https://CRAN.R-project.org/package=vegan
Olivares, M., Neef, A., Castillejo, G., Palma, G. D., Varea, V., Capilla, A., et al. (2015). The HLA DQ2 genotype selects for early intestinal microbiota composition in infants at high risk of developing coeliac disease. Gut 64 (3), 406–417. doi: 10.1136/gutjnl-2014-306931
Palma, G. D., Capilla, A., Nova, E., Castillejo, G., Varea, V., Pozo, T., et al. (2012). Influence of milk-feeding type and genetic risk of developing coeliac disease on intestinal microbiota of infants: the PROFICEL study. PLoS One 7 (2), e30791. doi: 10.1371/journal.pone.0030791
Park, J. Y., Seo, H., Kang, C. S., Shin, T. S., Kim, J. W., Park, J. M., et al. (2022). Dysbiotic change in gastric microbiome and its functional implication in gastric carcinogenesis. Sci. Rep. 12 (1), 4285. doi: 10.1038/s41598-022-08288-9
Penders, J., Thijs, C., Vink, C., Stelma, F. F., Snijders, B., Kummeling, I., et al. (2006). Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics 118 (2), 511–521. doi: 10.1542/peds.2005-2824
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41 (Database issue), D590–D596. doi: 10.1093/nar/gks1219
R Core Team (2016). : A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing).
R Core Team (2018). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing).
Roesch, L. F. W., Dobbler, P. T., Pylro, V. S., Kolaczkowski, B., Drew, J. C., Triplett, E. W. (2020). Pime: A package for discovery of novel differences among microbial communities. Mol. Ecol. Resour. 20 (2), 415–428. doi: 10.1111/1755-0998.13116
Russell, J. T., Roesch, L., Ördberg, M., Ilonen, J., Atkinson, M. A., Schatz, D. A., et al. (2019). Genetic risk for autoimmunity is associated with distinct changes in the human gut microbiome. Nat. Commun. 10 (1), 3621. doi: 10.1038/s41467-019-11460-x
Russell, J. T., Lauren Ruoss, J., de la Cruz, D., Li, N., Bazacliu, C., Patton, L., et al. (2021). Antibiotics and the developing intestinal microbiome, metabolome and inflammatory environment in a randomized trial of preterm infants. Sci Rep 11 (1), 1943. doi: 10.1038/s41598-021-80982-6
Shan, L., Molberg, Ø., Parrot, I., Hausch, F., Filiz, F., Gray, G. M., et al. (2002). Structural basis for gluten intolerance in celiac sprue. Sci. (New York N.Y.) 297 (5590), 2275–2279. doi: 10.1126/science.1074129
Shao, T., Shao, L., Li, H., Xie, Z., He, Z., Wen, C. (2017). Combined signature of the fecal microbiome and metabolome in patients with gout. Front. Microbiol. 8. doi: 10.3389/fmicb.2017.00268
Sharma, P., Baloda, V., Gahlot, G. P., Singh, A., Mehta, R., Vishnubathla, S., et al. (2019). Clinical, endoscopic, and histological differentiation between celiac disease and tropical sprue: A systematic review. J. Gastroenterol. Hepatol. 34 (1), 74–83. doi: 10.1111/jgh.14403
Sollid, L. M., Lie, B. A. (2005). Celiac disease genetics: current concepts and practical applications. Clin. Gastroenterol. hepatol: Off. Clin. Pract. J. Am. Gastroenterol Assoc. 3 (9), 843–851. doi: 10.1007/s00251-017-0985-7
Stadlbauer, V., Engertsberger, L., Komarova, I., Feldbacher, N., Leber, B., Pichler, , et al. (2020). Dysbiosis, gut barrier dysfunction and inflammation in dementia: a pilot study. BMC Geriatr. 20 (1), 248. doi: 10.1186/s12877-020-01644-2
Tollefsen, S., Arentz-Hansen, H., Fleckenstein, B., Molberg, O., Ráki, M., Kwok, W. W., et al. (2006). HLA-DQ2 and -DQ8 signatures of gluten T cell epitopes in celiac disease. J. Clin. Invest. 116 (8), 2226–2236. doi: 10.1172/JCI27620
Verdu, E. F., Schuppan, D. (2021). Co-Factors, microbes, and immunogenetics in celiac disease to guide novel approaches for diagnosis and treatment. Gastroenterology. 161 (5), 1395–1411.e4. doi: 10.1053/j.gastro.2021.08.016
Walker-Smith, J. A., Sanderson, I. R. (1999). Diagnostic criteria for chronic inflammatory bowel disease in childhood. nestle nutrition workshop series. Clin. Perform. Programme 2, 107–120. doi: 10.1159/000061783
Wickham, H. (2016). ggplot2: Elegant graphics for data analysis (Springer-Verlag New York). Available at: https://CRAN.R-project.org/package=vegan.
Wu, X., Qian, L., Liu, K., Wu, J., Shan, Z. (2021). Gastrointestinal microbiome and gluten in celiac disease. Ann. Med. 53 (1), 1797–1805. doi: 10.1080/07853890.2021.1990392
Zhang, J., Kobert, K., Flouri, T., Stamatakis, A.. (2014). PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 30 (5), 614–620. doi: 10.1093/bioinformatics/btt593
Keywords: celiac disease, autoimmunity, human leucocyte antigen, infant, gut microbiome
Citation: Milletich PL, Ahrens AP, Russell JT, Petrone JR, Berryman MA, Agardh D, Ludvigsson JF, Triplett EW and Ludvigsson J (2022) Gut microbiome markers in subgroups of HLA class II genotyped infants signal future celiac disease in the general population: ABIS study. Front. Cell. Infect. Microbiol. 12:920735. doi: 10.3389/fcimb.2022.920735
Received: 18 April 2022; Accepted: 24 June 2022;
Published: 25 July 2022.
Edited by:
Alberto Caminero, McMaster University, CanadaReviewed by:
Heather Galipeau, McMaster University, CanadaJose Ramon Bilbao, University of the Basque Country, Spain
Copyright © 2022 Milletich, Ahrens, Russell, Petrone, Berryman, Agardh, Ludvigsson, Triplett and Ludvigsson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eric W. Triplett, ewt@ufl.edu