Research article | Open | Open Peer Review | Published:
Duodenal microbiota composition and mucosal homeostasis in pediatric celiac disease
BMC Gastroenterologyvolume 13, Article number: 113 (2013)
Celiac disease (CD) is an autoimmune disorder of the small intestine which is triggered by dietary gluten in genetically predisposed (HLA-DQ2/DQ8 positive) individuals. Only a fraction of HLA-DQ2/DQ8 positive individuals develop CD indicating that other factors have a role in the disorder. Several studies have addressed intestinal microbiota aberrancies in pediatric CD, but the results are inconsistent. Previously, we demonstrated that pediatric CD patients have lower duodenal expression of TLR2 and higher expression of TLR9 as compared to healthy controls (HC) indicating that microbiota may have a role in CD.
We used bacterial phylogenetic microarray to comprehensively profile the microbiota in duodenal biopsies of CD (n = 10) and HC (n = 9) children. The expression of selected mucosa-associated genes was assessed by qRT-PCR in CD and HC children and in treated CD adults (T-CD, n = 6) on gluten free diet.
The overall composition, diversity and the estimated microbe associated molecular pattern (MAMP) content of microbiota were comparable between CD and HC, but a sub-population profile comprising eight genus-like bacterial groups was found to differ significantly between HC and CD. In HC, increased TLR2 expression was positively correlated with the expression of tight junction protein ZO-1. In CD and T-CD, the expression of IL-10, IFN-g and CXCR6 were higher as co5mpared to HC.
The results suggest that microbiota and altered expression of mucosal receptors have a role in CD. In CD subjects, the increased expression of IL-10 and IFN-g may have partly resulted from the increased TLR9 expression and signaling.
Celiac disease (CD) is a common chronic immune-mediated, inflammatory disorder of the small intestine induced by intolerance to gluten-containing dietary products [1, 2]. When a CD patient consumes gluten, an inflammatory cascade occurs in the small intestinal mucosa, eventually resulting in an active disease that is characterized by villous atrophy, crypt hyperplasia and increased numbers of lymphocytes in the lamina propria [1, 3]. Untreated CD is manifested by gastrointestinal symptoms, malabsorption and even malnutrition [1, 3] and also by extra-intestinal symptoms such as dermititis herpetiformis (skin rash) in some individuals . Both genetic predisposition and environmental factors are considered to be involved in the development of CD . Individuals who carry the alleles human leukocyte antigen (HLA)-DQ2 or HLA-DQ8 have an increased risk of developing the disease, but only less than one tenth of them eventually get CD, indicating that other genetic factors and/or environmental factors are also important in the pathogenesis . Recently, 13 new CD risk loci were identified, bringing the number of known CD loci to 40 and giving a refined picture of the genetic risk of CD . Most of these loci contain candidate genes of immunological function, but the pathways leading from genetic predisposition to an affected person are poorly understood [6, 7]. The intestinal microbiota is a major dictator of the antigen milieu of enterocytes, and it may have a role in the CD pathogenesis.
Gut colonization starts immediately after birth, and considerable microbiota maturation takes place during the first years of life followed by a gradual microbiota succession until young adulthood [8–10]. A balanced commensal microbiota contributes to the physiological development of the gut and the maturation of the immune system; thereby, alterations in the intestinal microbiota could play a role in the onset of different diseases, including CD [3, 11]. Increased bacterial diversity and changes in several bacterial groups in the microbiota of pediatric CD patients have been reported in several studies [12–16]. However, other recent studies have failed to show major microbiota differences between children with and without CD [17–19]. Two recent studies have addressed the microbiota in infants with a genetic predisposition to CD [4, 11]. Both of these reported the microbiota of predisposed infants to be different from that of non-predisposed, but while Sellitto et al. reported a reduction or lack of Bacteroides in predisposed infants , De Palma et al. found that Bacteroides fragilis and staphylococci were increased and bifidobacteria were reduced in genetically susceptible infants . Although the idea that the microbiota is involved in the etiology of CD has been addressed in numerous studies, the results on specific CD-associated microbiota changes remain inconclusive.
Recently, we demonstrated that pediatric CD patients have decreased duodenal expression of Toll-like receptor 2 (TLR2) and the negative regulator of Toll-receptor signaling (Tollip), and increased expression of TLR9 and interleukin 8 (IL-8), which is a marker of intestinal inflammation . TLRs are a family of pattern recognition receptors, which recognize conserved microbe-associated motifs such as lipopolysaccharides (LPS, ligand for TLR4), lipoproteins and lipoteichoeic acids (LTA, ligand for TLR2), flagellin (TLR5) and nucleic acid motifs (TLR3, 7 and 9). Intestinal epithelial homeostasis is dependent on the activation of TLRs at adequate levels in order to keep defense against microbes balanced and to avoid an excessive inflammatory response to gut commensals [20, 21]. Further, TLR signaling is known to affect epithelial barrier function by having a bearing on the expression of tight junction proteins, mucus protein mucin 2 and antimicrobial peptides such as RegIII-γ [20, 21]. Thus, our previous results suggested that microbiota may play a role in the etiology of CD, but we were unable to reveal aberrancies in the duodenal mucosal microbiota of pediatric CD by targeting selected bacterial groups by quantitative PCR (qPCR) .
The objectives of the present study were to comprehensively characterize the total duodenal mucosal microbiota and to re-evaluate the possible microbiota differences in pediatric CD patients and healthy controls by using a high-throughput bacterial phylogenetic microarray (HITChip). Further, the expression of nine mucosa-associated genes, IL-10, interferon-gamma (IFN-γ), tumor necrosis factor alpha (TNF-α), a tight junction protein zonula occludens-1 (ZO-1), a gap junction protein connexin-43 (Cx43), a mucus protein mucin 2 (MUC2), an antimicrobial peptide RegIIIγ, a chemokine CXCL16 and its receptor CXCR6 were measured by using a quantitative reverse transcription-PCR (qRT-PCR).
Study subjects and DNA and RNA extraction
Duodenal biopsy samples were collected from 20 Finnish children and 6 Finnish adults: 10 children with newly diagnosed CD (median age 9.5 ± 4.1 y; 3-14 y, 4 males and 6 females) before the implementation of gluten-free diet (GFD), 10 healthy control (HC) children (median age 8.5 ± 3.8 y; 4-16y, 4 males and 6 females) with gastrointestinal complaints or other reasons for esophagogastroduodenoscopy (abdominal pain - 4 patients, gastroesophageal reflux disease – 2 patients, growth retardation – 2 patients esophagitis -1 patient, achalasia – 1 patient) and 6 adults (median age 46 ± 11.4 y; 30-60 y) with CD who had been on a GFD at least for a year (treated CD, T-CD). All HC and T-CD subjects had both negative celiac serology and normal small intestinal mucosa (Marsh 0 lesions), while CD patients had both positive celiac serology markers (anti-tissue transglutaminase antibodies and/or anti-endomysium antibodies) and villous atrophy and crypt hyperplasia (Marsh III lesions) in duodenal biopsy. The study subjects were the same as in our previous study  and the same DNA and RNA preparations were used in this study. The DNA and RNA extraction and purification protocols were described in detail earlier .
Quantitative reverse-transcriptase-PCR (qRT-PCR)
Gene expression assays were performed using comparative Ct (threshold cycle)-method with ABI 7300 Real Time PCR System. (Applied Biosystems/Life Technologies Corporation, Carlsbad, CA). Taqman Gene Expression Assays (Applied Biosystems) used in analyses were: IL-10, assay ID: Hs00174086_m1; IFN-g, assay ID: Hs00174143_m1; TNF-α assay ID: Hs00174128; ZO-1, assay ID: Hs01551861_m1; Cx43, assay ID: Hs00748445_s1; MUC2, assay ID: Hs03005103_g1; RegIIIγ, assay ID: Hs01595405_m1; CXCL16, assay ID: Hs00222859_m1; CXCR6, assay ID: Hs01890898_s1. Gene expression assays were performed according to kit’s protocol. Reactions were run in three replicates in a total volume of 20 μl or 50 μl with 25 ng of cDNA in each. Thermal cycler conditions used were 1) 50°C for 2 min, 2) 95°C for 10 min, 3) 95°C for 15 sec, 4) 60°C for 1 min. Steps 3 and 4 were repeated 40 times. Gene expression of 18S rRNA was used as an endogenous control (a house-keeping gene) to normalize the gene expression. Negative control and Universe Human Reference RNA (Agilent Technologies, Santa Clara, CA) as a control RNA were included in every PCR run. Results were analyzed with RQ-Study program (Applied Biosystems) to receive Ct values for all samples. The relative expression of genes was then calculated as described in detail previously .
Microbiota analysis by phylogenetic microarray
The microbiota was analyzed with the bacterial phylogenetic microarray –the HITChip (Human Intestinal Tract Chip) [22–24]. This microarray consists of over 4 800 oligonucleotide probes targeting the V1 and V6 hypervariable regions of the 16S rRNA gene of 1038 intestinal phylotypes [22–24]. It allows a comprehensive and high-resolution analysis of the microbiota composition at different taxonomic levels. Genus-like level (L2) taxa correspond to bacteria having 90% or higher similarity in their 16S rRNA gene, whereas L1 taxa correspond to a phylum-like level .
The HITChip measurements were performed as previously described [22–24]. In brief, amplification of 16S rRNA gene was carried out from 100-200 ng of DNA with primer T7 prom-Bact-27 F and Bact-1369R , followed by in vitro transcription, dye labeling, fragmentation, and hybridization, as described earlier . The arrays were scanned with Agilent DNA Microarray Scanner G2505C (Agilent, USA) and the intensity values for each image were extracted from the generated images using Agilent Extraction Software, version 10.7.3.1. Normalization and quality control of HITChip array data were performed with scripts in R statistical software, as described earlier [22, 24, 26]. The technical replicates having a Pearson correlation over 0.94 were selected for further analysis and the replicates were averaged. The between sample normalization was carried out with min-max algorithm [27, 28]. Signal intensity threshold was applied to reduce experimental and possible cross-hybridization noise, as described previously .
Estimation of the microbe associated molecular pattern (MAMP) content
The relative content of selected MAMPs was estimated from the HITChip profiles summarizing the abundance of Gram-positive, Gram–negative or flagellated genus-like bacterial groups. The abundance of Gram-positive bacteria is assumed to reflect the load of LTA, i.e. TLR2 ligands, the abundance of Gram-negative bacteria the load of LPS i.e. TLR4 ligands, and the abundance of flagellated bacteria the load of TLR5 ligands. The general GC% of the microbiome was estimated based on the genomic GC-content of representative species from each of the genus-like group (Additional file 1: Table S1), weighted by the relative abundance of each of these groups in the total microbiota. The GC% is taken to reflect the load of unmethylated CpG motifs, i.e. TLR9 ligands, as a high GC content of a bacterial genome correlates with a higher number of CpG motifs in the genome (Kant R, de Vos WM, Palva A, Satokari R, unpublished results) [29, 30].
The data analysis was performed in R version 2.15.1 (R Development CT 2012). The sum of signal intensities for probes targeting a genus-like group was used as a quantitative measure of the abundance of the group in a sample . When computing the signal at higher level taxa, the probe intensities were divided by the number of known target phylotypes per probe. The signal intensity above threshold was log10-transformed.
The diversity of the microbial profiles was computed by Simpson’s reciprocal index of diversity (1/D) and Shannon indices on probe-level data [31, 32]. Principal component analysis (PCA) and Redundancy analysis (RDA) [33, 34] were computed using R packages ‘stats’ and ‘vegan’. The significance of separation in RDA was assessed with a permutation test  using 50000 permutations.
Bacterial groups that are associated with the health status were selected with 9-fold cross-validation (CV). The data was split into 9 subsets (folds) of equal size with stratification such that both CD and HC samples were present in all folds. In the CV procedure, eight subsets were used for training, and one subset in turn was left out for testing. Within each CV fold, six random forests were learned to predict the study group using a set of bacteria, selected by: i) choosing the 2, 4, 6, 8, or 10 bacteria having smallest p-values from two-sample t-test carried out within the training set, or ii) selecting the bacterial groups with high mean decrease in gini score from random forest using all bacteria . The prediction error of the random forests was then estimated with the test data. The feature selection method with lowest prediction error rate from cross-validation (that is, 8 bacteria with smallest p-values) was then applied for the full data, and its prediction error was estimated with 9-fold cross validation. Significance of the prediction was tested with permutation test using 10,000 permutations.
Differentially present bacterial groups were analyzed with Student’s t-test assuming two tailed distribution . False discovery rate (FDR) correction of p-values was carried out using Benjamini-Hochberg (BH) correction . The average relative abundance of each phylum- or genus-like group was estimated by first calculating the percentage of signal intensity within each sample, and then computing the average of percentage of signal intensities within CD and HC groups.
For relative gene expression, the normality of the data was tested by Shapiro-Wilk normality test function. When normality assumption applied, ANOVA and Tukey Honest Significant Differences post-hoc analysis was carried out. If the data residuals differed significantly from normality, Kruskal-Wallis and associated posthoc analysis were carried out . For all analyses, FDR corrected p-values below 0.05 were considered significant. The correlations between bacterial groups, gene expressions and MAMP content were estimated by Spearman coefficient, followed by FDR correction (BH) of p-value.
The study was accepted by the ethical committee of the Hospital District of Southwest Finland. Written informed consent was obtained from all of the study patients or their parents. Pediatric T-CD subjects could not be included as a control group in the study, because in Finland follow-up biopsies are not taken from children if celiac serology turns negative within two years after the implementation of GFD.
Composition of the duodenal mucosal microbiota
HITChip microbiota profiles were obtained from 10 CD subjects (4 males and 6 females) and 9 healthy controls (HC, median age 9 ± 4.1 years, 4 males and 5 females) (no significant difference in median age, one HC was not included due to a low reproducibility of HITChip profile, Additional file 2: Table S2).
At the bacterial phylum-like level, representatives from 13 groups were detected in the duodenal mucosa of the pediatric subjects (Table 1). The phylum-like level microbiota profiles were found to be individual-specific with large inter-individual variation (Figure 1). Proteobacteria, Bacilli and Bacteroidetes were found to be the major bacterial groups in the duodenal mucosa of both CD and HC (Figure 2). Proteobacteria and Bacilli constituted each approximately one third of the community while Bacteroidetes amounted to around 14% (Figure 2). There were no significant differences in the abundance of bacterial phylum-like groups between CD and HC. Similarly, the bacterial diversity was comparable between the HC and CD groups (Figure 3).
Altogether 65 genus-like groups were detected in the duodenal biopsies (Table 1). HC and CD shared the same predominant genus-like groups, whose average proportion was above 5% of the total signal. The predominant groups were Sutterella wadsworthensis et rel., Streptococcus mitis et rel., Aquabacterium, Streptococcus bovis et rel., Streptococcus intermedius et rel., and Prevotella melaninogenica et rel. (Table 1). Among these, the most abundant bacterial groups in both HC and CD were Sutterella wadsworthensis et rel. and Streptococcus mitis et rel., both with an average abundance of 14 to 19% (Table 1). No single genus-like bacterial group abundance differed significantly between HC and CD (Table 1).
The general MAMP content of the duodenal microbiota
The general MAMP content was estimated based on the microbiota profiles (Table 1). No significant difference between the abundance of Gram-positive or Gram-negative bacteria carrying LTA or LPS respectively was found, nor in the abundance of potentially flagellated bacteria between HC and CD. Thus, the ligand load for TLR2, 4 and 5 seems comparable between HC and CD. In addition, the average GC% content of the total microbiota and the relative abundance of high-GC% bacteria (genomic GC% > 58%) were comparable between the groups suggesting a similar load of unmethylated CpG motifs, which are ligands for the TLR9.
Celiac disease associated microbiota profile
In PCA and RDA plots, HC vs. CD subjects did not cluster separately, showing that the groups do not differ in terms of the total microbiota profile (Additional file 3: Figure S1, Additional file 4: Figure S2). Moreover, as already mentioned above, no single genus-like bacterial group’s abundance differed significantly between HC and CD (Table 1).
Further analysis was performed with random forests using feature selection as preprocessing step, where the best method was chosen from six candidates using cross-validation. The aim was to identify a possible health status-related bacterial sub-population within the total microbiota. The feature selection method resulting in lowest cross-validated prediction error was to choose a set of eight genus-like bacteria having lowest p-values from t-test (Figure 4). Thus, although no single genus-like group individually differed significantly between HC and CD, a random forest learned with a profile of selected eight genus-like groups predicted the health status with error rate of 31.6%. This is significantly better than random guess, having a median error rate of 53%, and 95% confidence intervals of [32%, 74%].
As can be seen from Table 1, P. melaninogenica et rel. has a relative abundance of 5.3% and 6.5% in HC and CD, respectively and it is among the predominant genus-like groups contributing to the separation between HC and CD. Also Haemophilus ssp. and Serratia ssp. had relatively higher abundance in CD, whereas the other five bacterial groups were higher in HC, see Table 1.
Human host gene expression
The gene expression levels of ZO-1, CXCL16, CXCR6, IL-10, IFN-γ, TNF-α, Cx43, MUC2, and RegIIIγ in duodenal biopsies were assessed by relative quantitative reverse transcription-PCR in the HC and CD children. Biopsies from six T-CD adults were included for comparison, because biopsies from T-CD children were not available (see ethical considerations above). The expression of ZO-1, CXCL16, CXCR6 could be assessed only from 8 HC and 9 CD subjects, because RNA of some of the samples did not suffice after the other qRT-PCR analysis (Figures 5 and 6).
The gene expression of ZO-1 was decreased in CD as compared to T-CD (Figure 5). The expression of CXCR6 was higher in CD and T-CD as compared to HC (Figure 6A). The expressions of IL-10 and IFN-γ were higher, whereas the ratio of IL-10 to IFN-γ was significantly reduced in both CD and T-CD as compared to HC (Figure 7). The gene expressions of TNF-α, Cx43, MUC2, RegIIIγ (Table A2), and CXCL16 (Figure 6B) were found to be comparable between the groups.
When correlating of TLR2 expression  with the expression of proteins related to the physical barrier function or the production of antimicrobial peptides (Cx43, ZO-1, MUC2 and RegIIIγ), the ZO-1 expression correlated positively with the TLR2 expression in the HC group (p = 0.03, rho = 0.88). The expression of Cx43 negatively correlated with the relative abundance of bacteria related to Solobacterium moorei (previously Bulleidia moorei) et rel. (p = 0.03, rho = -0.82). No other correlations were found between the relative abundance of genus-like groups or the estimated MAMP content of the microbiota and the gene expressions (data not shown).
Human small intestine, especially duodenum, is a distinctive environment for microbial life because of the excretion of digestive enzymes and bile, and it harbors a microbiota that is different from the one in the large intestine . The HITChip analysis showed that representatives of the genus Streptococcus constitute 29% and 32% of the total duodenal microbiota in HC and CD subjects, respectively, and that the signal obtained for all Bacilli (phylum-like group) predominantly resulted from streptococci. This is in line with the previous findings [39, 40]. In addition to Bacilli, Proteobacteria, Bacteroidetes and Clostridium cluster XIVa were found to be the most abundant phylum-like groups in the duodenum. These bacteria have previously been found to be dominant groups in the ileum, distal duodenum and proximal jejunum [17, 39, 41] and in the duodenum . Interestingly, within Proteobacteria a single most abundant genus-like group was found to be Sutterella wadsworthensis et rel. both in CD and HC subjects, which suggests that it belongs to the normal duodenal microbiota. Similarly, Mukhopadhya et al. detected S. wadsworthensis with an equally high frequency of 84 to 86% in the colonic biopsies of healthy adults (n = 64) or ulcerative colitis patients (n = 69), showing that the species belongs to the normal microbiota of intestinal mucosa and is unlikely to have a role in IBD . Genus Aquabacterium belonging to the Proteobacteria also had a relatively high abundance (~9.5%, prevalence 100% in both CD and HC) in the duodenum of both groups of children. Previously, Aquabacterium has been detected in human colonic mucosal biopsy , but to our knowledge, this is the first time that Aquabacterium has been described as an abundant inhabitant of the human small intestine.
In the PCA and RDA analysis, HC and CD subjects did not cluster separately. Furthermore, none of the 65 genus-like bacterial groups was found to be significantly different in abundance between HC and CD. Thus, the overall duodenal microbiota composition seems comparable between HC and CD, which is in line with the results obtained by Ou et al  and Nistal et al  using small intestinal biopsies. Moreover, the bacterial diversity was also found to be comparable between the study groups. Previously, bacterial diversity assessed by PCR-T/DGGE has been observed to be increased in children with CD [13, 16]. It is noteworthy that PCR-D/TGGE analysis detects only the most abundant bacteria and therefore may strongly underestimate microbiota diversity in complex communities. Our results from high-throughput microbiota profiling, like the results by Nistal et al , give a more in depth view of the duodenal microbiota regarding the bacterial groups inhabiting duodenal mucosa and the overall diversity. In addition, several studies have found differences in specific bacterial groups between HC and CD [12, 14, 15, 44]. The most consistent findings from these studies were that children with CD have increased counts of Bacteroides and reduced counts of bifidobacteria either in feces or duodenal epithelium [12, 14, 15] which were not found to differ between CD and HC in this study.
As the overall microbiota profile, diversity or individual genus-like groups did not show a significant difference between HC and CD, random forest was used to explore whether a sub-population bacterial profile could be associated with the health status. A profile of eight bacterial groups was found to distinguish HC from CD. The abundance of Prevotella melaninogenica, and the total abundance of Prevotella spp (Table 1) were found to be higher in CD, which is in line with the results of Ou et al. Haemophilus et rel. were also found to be enriched, although insignificantly as an individual bacterial group, in CD children by Nistal et al. . Serratia spp. was also found to be present in higher abundance in CD. The phylogenetic microarray targets mainly S. marcescens, which is considered as an opportunistic pathogen able to cause invasive infections (sepsis, meningitis, pneumonitis) in neonates [45, 46]. These bacteria may impair the intestinal integrity, but their possible role in CD remains to be elucidated.
In the distinctive profile detected by random forest, P. oralis, R. bromii, P. cinnamivorans, Proteus and C. stercorarium groups are increased in the HC group. R. bromii is of particular interest, because it acts as key species for fiber/resistant starch degradation in the intestine  feeding butyrate producing bacteria including P. cinnamivorans. Butyrate is a major source of energy to the enterocytes and acts as a regulator of gene expression, inflammation and differentiation in host cells . Previously, R. bromii has been detected at increased level in healthy subjects as compared to Crohn’s disease patients , indicating its potential role for benefiting the gut mucosal homeostasis.
Finally, it should be taken into account that the highly individual-specific microbiota-profiles may have a strong impact on the results within small study groups. Moreover, children in the HC group had healthy duodenal mucosa, but had gastrointestinal complaints or other reasons for gastroscopy and the possibility of microbiota alterations in these children as compared to children devoid of any symptoms can´t be excluded. Therefore, the profile comprising eight genus-like bacterial groups, which showed significant difference between CD and HC should be verified in future studies.
Similar average proportion of Gram-positive and Gram-negative bacteria, which are the carriers of LTA and LPS, respectively, was found in CD and HC. Our previous findings from the same cohort showed that the expression of TLR2 is higher in HC . As the LTA content is comparable between HC and CD, the net signaling through TLR2 is presumably higher in HC. In vitro, TLR2 stimulation of intestinal epithelial cells (IECs) has been shown to increase Cx43 synthesis, the apical reorganization of ZO-1 and trans-epithelial resistance, which reflects the strength of tight junctions between IECs and barrier function [51–53]. Although no significant difference was found in the expression ZO-1 and Cx43 between HC and CD, the expression of ZO-1 was significantly lower in CD as compared to T-CD. Moreover, there was a significant positive correlation between the expressions of ZO-1 and TLR2 in HC further supporting earlier findings, which demonstrated more permeable epithelial barrier in CD as compared to HC due to the decreased expression of tight junction proteins . Negative correlation was found between the abundance of S. moorei et rel. and the expression of Cx43. Previously this bacterium has been associated with oral cavity diseases  and its possible down-regulatory effect of Cx43 may facilitate invasion and colonization. In CD, however, S. moorei et rel. seems to be irrelevant as it was found to be as abundantly present in HC and CD.
CXCL16 has been shown to function as a scavenger receptor in antigen-presenting cells where it mediates adhesion and phagocytosis of both Gram-positive and Gram-negative bacteria . In addition, it works as a chemokine for CXCR6-expressing cells such as natural killer T cells and T helper 1 (Th1) - polarized CD4 T cells . Since both microbiota alterations and Th1-polarized inflammation have been linked to the pathogenesis of celiac disease, we decided to evaluate the expression of both CXCL16 and its receptor CXCR6 in this study. We found that the expression of cytokine CXCL16 was found to be comparable in CD and HC and T-CD, whereas the expression of its receptor CXCR6 was higher in CD and T-CD as compared to HC. CXCL16 has dual functions as a transmembrane adhesion molecule and a soluble chemokine . Both the membrane-bound form of CXCL16 and its receptor CXCR6 have been found to be expressed not only by dendritic cells/macrophages and T cells respectively, but also by IECs [58–60]. Previously, increased intestinal CXCL16 expression has been observed in the colonic biopsies of Crohn´s disease patients due to immune cell infiltration . In this study, an increased expression of CXCR6 was not only observed in the inflamed mucosa of CD but also in T-CD without immune cell infiltration suggesting an altered expression in epithelial cells. However, increased expression of CXCR6 in mucosa associated immune cells cannot be excluded when studying whole biopsies with different cell types. Diegelmann et al observed that in vitro stimulation of CXCR6 activates several distinct signaling pathways in IECs and they suggested that the CXCL16-CXCR6 chemokine-receptor system contributes to the integrity of epithelium and the regulation of mucosal innate and adaptive immune systems . To our knowledge our data report for the first time an increased duodenal expression of CXCR6 in CD subjects in whom it may have an important role in the mucosal immunity.
The average genomic GC content (GC%) of the total microbiota and the proportion of high-GC% bacteria was found to be comparable between the HC and CD groups. High GC content of a bacterial genome correlates with a higher number of potentially immunostimulatory CpG motifs in the genome (Kant R, de Vos WM, Palva A, Satokari R, unpublished results) [29, 30] and thus, the GC% was taken to reflect the load of TLR9 ligands in the microbiota. Consequently, HC and CD are likely to harbor similar load of TLR9 ligands. Our previous results from the same cohort showed that the expression of TLR9 in duodenum is higher in CD . Thus, the signaling through TLR9 is presumably higher in CD as compared to HC due to the equally high load of TLR9 ligands and higher expression of TLR9. In immune cells, TLR9 stimulation is known to trigger Th1 type immune responses [61–63]. Further, a recent in vitro study demonstrated that apical TLR9 stimulation of IECs leads to an increased expression of IFN-γ and IL-10, but not TNF-α from peripheral blood mononuclear cells (PBMCs) on the basolateral side . In CD patients, increased mucosal expression of IL-10 and IFN-γ is well documented and has been associated with the activity of gliadin-reactive T-cells upon the gluten/gliadin stimulation either in vivo or in vitro[65–69]. We found that the expression of IL-10 and IFN-γ in the duodenal mucosa was increased not only in CD but also in T-CD with gluten-free diet, suggesting for an additional gluten/gliadin-independent route of stimulation. Similarly to previous studies [65, 68, 69], TNF-α expression was found to be unaffected in CD and T-CD as compared to HC. The ratio of IL-10 to IFN-γ transcripts was strongly reduced in both CD and T-CD as compared to HC suggesting that although IL-10 expression was also increased there is a significant inclination towards a Th1 response both in untreated and treated CD. Since both CD and T-CD had increased expression of TLR9 , we hypothesize that the increased TLR9 signaling in the small intestine may contribute to the persistent activation of Th1 (IFN-γ) signaling pathway markers in the small intestine found in CD children despite gluten-free diet treatment . Due to the limitations in studying whole biopsies with different cell types we cannot conclude which of the proposed signaling routes i.e. the direct TLR9 stimulation of immune cells or the stimulation of IECs with subsequent triggering of immune cells would possibly encompass more into the observed cytokine profile.
Our results suggest that intestinal microbiota and host-microbe cross-talk play a role in the CD. While the overall microbiota composition in the duodenal mucosa was comparable between the CD and healthy children, a sub-population profile comprising eight genus-like bacterial groups was found to differ significantly between the study groups. The sub-population of bacteria has potentially a specific role in e.g. maintaining gut homeostasis in the healthy individuals or in compromising the epithelial function in CD. Secondly, these and our previous results from the same cohort  suggest that altered expression of mucosal receptors and epithelial host-microbe cross-talk have a role in CD. We hypothesize that increased TLR9 signaling in the duodenum may contribute to the Th1 response (increased IFN-γ) found in the small intestinal mucosa of CD subjects even after the implementation of GFD.
Chemokine (C-X-C motif) receptor 6
Chemokine (C-X-C motif) ligand 16
- et rel:
False discovery rate
Human intestinal tract chip
Human leukocyte antigen
Intestinal epithelial cells
Microbe associated molecular pattern
Peripheral blood mononuclear cells
Principal component analysis
Treated celiac disease
Tumor necrosis factor alpha
Green P, Jabri B: Celiac disease. Annu Rev Med. 2006, 57: 207-221. 10.1146/annurev.med.57.051804.122404.
Papista C, Gerakopoulos V, Kourelis A, Sounidaki M, Kontana A, Berthelot L, Moura IC, Monteiro RC, Yiangou M: Gluten induces coeliac-like disease in sensitised mice involving IgA, CD71 and transglutaminase 2 interactions that are prevented by probiotics. Lab Invest. 2012, 92: 625-635. 10.1038/labinvest.2012.13.
Soni S, Badawy SZA: Celiac disease and its effect on human reproduction a review. J Reprod Med. 2010, 55: 3-8.
Sellitto M, Bai G, Serena G, Fricke WF, Sturgeon C, Gajer P, White JR, Koenig SSK, Sakamoto J, Boothe D, Gicquelais R, Kryszak D, Puppa E, Catassi C, Ravel J, Fasano A: Proof of concept of microbiome-metabolome analysis and delayed gluten exposure on celiac disease autoimmunity in genetically at-risk infants. PLoS One. 2012, 7: e33387-10.1371/journal.pone.0033387.
Trynka G, Wijmenga C, van Heel DA: A genetic perspective on coeliac disease. Trends Mol Med. 2010, 16: 537-550. 10.1016/j.molmed.2010.09.003.
Trynka G, Hunt KA, Bockett NA, Romanos J, Mistry V, Szperl A, Bakker SF, Bardella MT, Bhaw-Rosun L, Castillejo G, de la Concha EG, de Almeida RC, Dias KRM, van Diemen CC, Dubois PCA, Duerr RH, Edkins S, Franke L, Fransen K, Gutierrez J, Heap GAR, Hrdlickova B, Hunt S, PlazaIzurieta L, Izzo V, Joosten LAB, Langford C, Mazzilli MC, Mein CA, Midah V, et al: Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nat Genet. 2011, 43: 1193-1201. 10.1038/ng.998.
Sperandeo MP, Tosco A, Izzo V, Tucci F, Troncone R, Auricchio R, Romanos J, Trynka F, Auricchio S, Jabri B, Greco L: Potential celiac patients: a model of celiac disease pathogenesis. PLoS One. 2011, 6: e21281-10.1371/journal.pone.0021281.
Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5: e177-10.1371/journal.pbio.0050177.
Nylund L, Satokari R, Nikkilä J, Rajilic-Stojanovic M, Kalliomäki M, Isolauri E, Salminen S, de Vos WM: Microarray analysis reveals marked intestinal microbiota aberrancy in infants having eczema compared to healthy children in at-risk for atopic disease. BMC Microbiol. 2013, 13: 12-10.1186/1471-2180-13-12.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Gordon JI: Human gut microbiome viewed across age and geography. Nature. 2012, 486: 222-227.
De Palma G, Capilla A, Nova E, Castillejo G, Varea V, Pozo T, Antonio GJ, Polanco I, Lopez A, Ribes-Koninckx C, Marcos A, Dolores GNM, Calvo C, Ortigosa L, Pena-Quintana L, Palau F, Sanz Y: Influence of milk-feeding type and genetic risk of developing coeliac disease on intestinal microbiota of infants: the PROFICEL study. PLoS One. 2012, 7: e30791-10.1371/journal.pone.0030791.
Nadal I, Donant E, Ribes-Koninckx C, Calabuig M, Sanz Y: Imbalance in the composition of the duodenal microbiota of children with coeliac disease. J Med Microbiol. 2007, 56: 1669-1674. 10.1099/jmm.0.47410-0.
Sanz Y, Sanchez E, Marzotto M, Calabuig M, Torriani S, Dellaglio F: Differences in faecal bacterial communities in coeliac and healthy children as detected by PCR and denaturing gradient gel electrophoresis. FEMS Immunol Med Microbiol. 2007, 51: 562-568. 10.1111/j.1574-695X.2007.00337.x.
Collado MC, Donat E, Ribes-Koninckx C, Calabuig M, Sanz Y: Specific duodenal and faecal bacterial groups associated with paediatric coeliac disease. J Clin Pathol. 2009, 62: 264-269. 10.1136/jcp.2008.061366.
De Palma G, Nadal I, Medina M, Donat E, Ribes-Koninckx C, Calabuig M, Sanz Y: Intestinal dysbiosis and reduced immunoglobulin-coated bacteria associated with coeliac disease in children. BMC Microbiol. 2010, 10: 63-10.1186/1471-2180-10-63.
Schippa S, Iebba V, Barbato M, Di Nardo G, Totino V, Checchi MP, Longhi C, Maiella G, Cucchiara S, Conte MP: A distinctive 'microbial signature' in celiac pediatric patients. BMC Microbiol. 2010, 10: 175-10.1186/1471-2180-10-175.
Ou G, Hedberg M, Hörstedt P, Baranov V, Forsberg G, Drobni M, Sandström O, Wai SN, Johansson I, Hammarström ML, Hernell O, Hammarström S: Proximal small intestinal microbiota and identification of rod-shaped bacteria associated with childhood celiac disease. Am J Gastroenterol. 2009, 104: 3058-3067. 10.1038/ajg.2009.524.
Kalliomäki M, Satokari R, Lähteenoja H, Vähämiko S, Grönlund J, Routi T, Salminen S: Expression of microbiota, Toll-like receptors, and their regulators in the small intestinal mucosa in celiac disease. J Pediatr Gastroenterol Nutr. 2012, 54: 727-732. 10.1097/MPG.0b013e318241cfa8.
Nistal E, Caminero A, Herran AR, Arias L, Vivas S, Ruiz de Morales JM, Jose M, Calleja S, Saenz De M, Luis E, Arroyo P, Casqueiro J: Differences of small intestinal bacteria populations in adults and children with/without celiac disease: effect of age, gluten diet, and disease. Inflamm Bowel Dis. 2012, 18: 649-656. 10.1002/ibd.21830.
Maynard CL, Elson CO, Hatton RD, Weaver CT: Reciprocal interactions of the intestinal microbiota and immune system. Nature. 2012, 489: 231-241. 10.1038/nature11551.
Abreu MT: Toll-like receptor signalling in the intestinal epithelium: how bacterial recognition shapes intestinal function. Nat Rev Immunol. 2010, 10: 131-144. 10.1038/nri2707.
Rajilic-Stojanovic M, Heilig HGHJ, Molenaar D, Kajander K, Surakka A, Smidt H, de Vos WM: Development and application of the human intestinal tract chip, a phylogenetic microarray: analysis of universally conserved phylotypes in the abundant microbiota of young and elderly adults. Environ Microbiol. 2009, 11: 1736-1751. 10.1111/j.1462-2920.2009.01900.x.
Biagi E, Nylund L, Candela M, Ostan R, Bucci L, Pini E, Nikk J, Monti D, Satokari R, Franceschi C, Brigidi P, De Vos W: Through ageing, and beyond: gut microbiota and inflammatory status in seniors and centenarians. PLoS One. 2010, 5: e10667-10.1371/journal.pone.0010667.
Jalanka-Tuovinen J, Salonen A, Nikkilä J, Immonen O, Kekkonen R, Lahti L, Palva A, de Vos WM: Intestinal Microbiota in Healthy Adults: Temporal Analysis Reveals Individual and Common Core and Relation to Intestinal Symptoms. PLoS One. 2011, 6: e23035-10.1371/journal.pone.0023035.
Iwamoto T, Tani K, Nakamura K, Suzuki Y, Kitaqawa M, Equchi M, Nasu M: Monitoring impact of in situ biostimulation treatment on groundwater bacterial community by DGGE. FEMS Microbiol Ecol. 2000, 32: 129-141. 10.1111/j.1574-6941.2000.tb00707.x.
Salonen A, Nikk J, Jalanka-Tuovinen J, Immonen O, Rajilic-Stojanovic M, Kekkonen RA, Palva A, de Vos WM: Comparative analysis of fecal DNA extraction methods with phylogenetic microarray: effective recovery of bacterial and archaeal DNA using mechanical cell lysis. J Microbiol Methods. 2010, 81: 127-134. 10.1016/j.mimet.2010.02.007.
Bolstad B, Irizarry R, Astrand M, Speed T: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.
Han JKM: Data mining: concepts and techniques. 2006, San Francisco, CA: Morgan Kaufmann Publishers
Menard O, Gafa V, Kapel N, Rodriguez B, Butel M, Waligora-Dupriet A: Characterization of immunostimulatory CpG-rich sequences from different Bifidobacterium species. Appl Environ Microbiol. 2010, 76: 2846-2855. 10.1128/AEM.01714-09.
Lee KW, Jung JW, Lee Y, Kim TY, Choi SY, Park J, Kim DS, Kwon HJ: Immunostimulatory oligodeoxynucleotide isolated from genome wide screening of Mycobacterium bovis chromosomal DNA. Mol Immunol. 2006, 43: 2107-2118. 10.1016/j.molimm.2005.12.004.
Simpson E: Measurement of diversity. Nature. 1949, 163: 688-688-
Hill MO: Diversity and evenness - a unifying notation and its consequences. Ecology. 1973, 54: 427-432. 10.2307/1934352.
Legendre PLL: Numerical Ecology. 1998, Amsterdam, NL: Elservier Science BV
Venables WNRB: Modern applied statistics with S. 2002, New York, NY: Springer
Breiman L: Random forests. Mach Learn. 2001, 45 (1): 5-32. 10.1023/A:1010933404324.
Chambers JM, Hastie TJ: Statistical models in S. 1992, Pacific Grove, CA: Wasdsworth & Brooks/Cole
Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125: 279-284. 10.1016/S0166-4328(01)00297-2.
Siegel S, Castellan NJ: Non parametric statistics for the behavioural sciences. 1988, New York, NY: MacGraw-Hill
Zoetendal EG, Raes J, van de Bogert B, Arumugam M, Booijink CC, Troost FJ, Bork P, Wels M, de Vos WM, Kleerebezem M: The human small intestinal microbiota is driven by rapid uptake and conversion of simple carbohydrates. ISME J. 2012, 6: 1415-1426. 10.1038/ismej.2011.212.
Sullivan A, Tornblom H, Lindberg G, Hammarlund B, Palmgren AC, Einarsson C, Nord CE: The micro-flora of the small bowel in health and disease. Anaerobe. 2003, 9: 11-14. 10.1016/S1075-9964(03)00060-X.
Ricanek P, Lothe SM, Frye SA, Rydning A, Catn MH, Tønjum T: Gut bacterial profile in patients newly diagnosed with treatment-naive Crohn's disease. Clin Exp Gastroenterol. 2012, 5: 173-186.
Mukhopadhya I, Hansen R, Nicholl CE, Alhaidan YA, Thomson JM, Berry SH, Pattinson C, Stead DA, Russell RK, EI-Omar EM, Hold GL: A comprehensive evaluation of colonic mucosal isolates of Sutterella wadsworthensis from inflammatory bowel disease. PLoS One. 2011, 6: e27076-10.1371/journal.pone.0027076.
Wang X, Heazlewood SP, Krause DO, Florin TH: Molecular characterization of the microbial species that colonize human ileal and colonic mucosa by using 16S rDNA sequence analysis. J Appl Microbiol. 2003, 95: 508-520. 10.1046/j.1365-2672.2003.02005.x.
Sanchez E, Donat E, Ribes-Koninckx C, Calabuig M, Sanz Y: Intestinal Bacteroides species associated with coeliac disease. J Clin Pathol. 2010, 63: 1105-1111. 10.1136/jcp.2010.076950.
Polilli E, Parruti G, Fazii P, D'Antonio D, Palmieri D, D'Incecco C, Mangifesta A, Garofalo G, Del Duca L, D'Amario C, Scimia M: Rapidly controlled outbreak of Serratia marcescens infection/colonisations in a neonatal intensive care unit, Pescara General Hospital, Pescara, Italy. Eurosurveillance. 2011, 16: 16-18.
Adamson V, Mitt P, Pisarev H, Metsvaht T, Telling K, Naaber P, Maimets M: Prolonged outbreak of Serratia marcescens in Tartu University Hospital: a case–control study. BMC Infect Dis. 2012, 12: 281-10.1186/1471-2334-12-281.
Ze X, Duncan SH, Louis P, Flint HJ: Ruminococcus bromii is a keystone species for the degradation of resistant starch in the human colon. ISME J. 2012, 6: 1535-1543. 10.1038/ismej.2012.4.
Defnoun S, Labat M, Ambrosio M, Garcia JL, Patel BKC: Papillibacter cinnamivorans gen. nov., sp nov., a cinnamate-transforming bacterium from a shea cake digester. Int J Syst Evol Microbiol. 2000, 50: 1221-1228. 10.1099/00207713-50-3-1221.
Flint HJ, Scott KP, Louis P, Duncan SH: The role of the gut microbiota in nutrition and health. Nat Rev Gastroenterol Hepatol. 2012, 9: 577-589. 10.1038/nrgastro.2012.156.
Mondot S, Kang S, Furet JP, de Carcer DA, McSweeney C, Morrison M, Marteau P, Dore J, Leclerc M: Highlighting new phylogenetic specificities of Crohn's disease microbiota. Inflamm Bowel Dis. 2011, 17: 185-192. 10.1002/ibd.21436.
Cario E, Gerken G, Podolsky DK: Toll-like receptor 2 enhances ZO-1-associated intestinal epithelial barrier integrity via protein kinase C. Gastroenterology. 2004, 127: 224-238. 10.1053/j.gastro.2004.04.015.
Cario E, Gerken G, Podolsky DK: Toll-like receptor 2 controls mucosal inflammation by regulating epithelial barrier function. Gastroenterology. 2007, 132: 1359-1374. 10.1053/j.gastro.2007.02.056.
Ey B, Eyking A, Gerken G, Podolsky DK, Cario E: TLR2 mediates gap junctional intercellular communication through connexin-43 in intestinal epithelial barrier injury. J Biol Chem. 2009, 284: 22332-22343. 10.1074/jbc.M901619200.
Pizzuti D, Bortolami M, Mazzon E, Buda A, Guariso G, D'Odorico A, Chiarelli S, D'Inca R, De Lazzari F, Martines D: Transcriptional downregulation of tight junction protein ZO-1 in active coeliac disease is reversed after a gluten-free diet. Dig Liver Dis. 2004, 36: 337-341. 10.1016/j.dld.2004.01.013.
Pedersen RM, Holt HM, Justesen US: Solobacterium moorei bacteremia: identification, antimicrobial susceptibility, and clinical characteristics. J Clin Microbiol. 2011, 49: 2766-2768. 10.1128/JCM.02525-10.
Shimaoka T, Nakayama T, Kume N, Takahashi S, Yamaguchi J, Minami M, Hayashida K, Kita T, Ohsumi J, Yoshie O, Yonehara S: Cutting edge: SR-PSOX/CXC chemokine ligand 16 mediates bacterial phagocytosis by APCs through its chemokine domain. J Immunol. 2003, 171: 1647-1651.
Uza N, Nakase H, Yamamoto S, Yoshino T, Takeda Y, Ueno S, Inoue S, Mikami S, Matsuura M, Shimaoka T, Kume N, Minami M, Yonehara S, Ikeuchi H, Chiba T: SR-PSOX/CXCL16 plays a critical role in the progression of colonic inflammation. Gut. 2011, 60: 1494-1505. 10.1136/gut.2010.221879.
Abel S, Hundhausen C, Mentlein R, Schulte A, Berkhout TA, Broadway N, Hartmann D, Sedlacek R, Dietrich S, Muetze B, Schuster B, Kallen KJ, Saftig P, Rose-John S, Ludwig A: The transmembrane CXC-chemokine ligand 16 is induced by IFN-gamma and TNF-alpha and shed by the activity of the disintegrin-like metalloproteinase ADAM10. J Immunol. 2004, 172: 6362-6372.
Hase K, Murakami T, Takatsu H, Shimaoka T, Iimura M, Hamura K, Kawano K, Ohshima S, Chihara R, Itoh K, Yonehara S, Ohno H: The membrane-bound chemokine CXCL16 expressed on follicle-associated epithelium and m cells mediates lympho-epithelial interaction in GALT. J Immunol. 2006, 176: 43-51.
Diegelmann J, Seiderer J, Niess J, Haller D, Goeke B, Reinecker H, Brand S: Expression and regulation of the chemokine CXCL16 in Crohn's disease and models of intestinal inflammation. Inflamm Bowel Dis. 2010, 16: 1871-1881. 10.1002/ibd.21306.
Hemmi H, Takeuchi O, Kawai T, Kaisho T, Sato S, Sanjo H, Matsumoto M, Hoshino K, Wagner H, Takeda K, Akira S: A Toll-like receptor recognizes bacterial DNA. Nature. 2000, 408: 740-745. 10.1038/35047123.
Krieg A: CpG motifs in bacterial DNA and their immune effects. Annu Rev Immunol. 2002, 20: 709-760. 10.1146/annurev.immunol.20.100301.064842.
Yu D, Putta MR, Bhagat L, Li Y, Zhu F, Wang D, Tang JX, Kandimalla ER, Agrawal S: Agonists of Toll-like receptor 9 containing synthetic dinucleotide motifs. J Med Chem. 2007, 50: 6411-6418. 10.1021/jm070881l.
de Kivit S, van Hoffen E, Korthagen N, Garssen J, Willemsen LEM: Apical TLR ligation of intestinal epithelial cells drives a T(h)1-polarized regulatory or inflammatory type effector response in vitro. Immunobiology. 2011, 216: 518-527. 10.1016/j.imbio.2010.08.005.
Forsberg G, Hernell O, Melgar S, Israelsson A, Hammarstrom S, Hammarstrom ML: Paradoxical coexpression of proinflammatory and down-regulatory cytokines in intestinal T cells in childhood celiac disease. Gastroenterology. 2002, 123: 667-678. 10.1053/gast.2002.35355.
Salvati VM, Mazzarella G, Gianfrani C, Levings MK, Stefanile R, De Giulio B, Iaquinto G, Giardullo N, Auricchio S, Roncarolo MG, Troncone R: Recombinant human interleukin 10 suppresses gliadin dependent T cell activation in ex vivo cultured coeliac intestinal mucosa. Gut. 2005, 54: 46-53. 10.1136/gut.2003.023150.
Gianfrani C, Levings MK, Sartirana C, Mazzarella G, Barba G, Zanzi D, Camarca A, Iaquinto G, Giardullo N, Auricchio S, Troncone R, Roncarolo MG: Gliadin-specific type 1 regulatory T cells from the intestinal mucosa of treated celiac patients inhibit pathogenic T cells. J Immunol. 2006, 177: 4178-4186.
Nilsen EM, Lundin KEA, Krajci P, Scott H, Sollid LM, Brandtzaeg P: Gluten specific, Hla-Dq restricted T-cells from celiac mucosa produce cytokines with Thl or Th0 profile dominated by interferon-gamma. Gut. 1995, 37: 766-776. 10.1136/gut.37.6.766.
Nilsen EM, Jahnsen FL, Lundin KEA, Johansen FE, Fausa O, Sollid LM, Jahnsen J, Scott H, Brandtzaeg P: Gluten induces an intestinal cytokine response strongly dominated by interferon gamma in patients with celiac disease. Gastroenterology. 1998, 115: 551-563. 10.1016/S0016-5085(98)70134-9.
Lahdenperä A, Ludvigsson J, Falth-Magnusson K, Högberg L, Vaarala O: The effect of gluten-free diet on Th1-Th2-Th3-associated intestinal immune responses in celiac disease. Scand J Gastroenterol. 2011, 46: 538-549. 10.3109/00365521.2011.551888.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-230X/13/113/prepub
Riikka Lankinen is acknowledged for technical assistance. The financial support from the Academy of Finland (grant nr. 138902), the Finnish Medical Foundation, the Paediatric Research Foundation, the Hospital District of Southwest Finland, the Finnish Coeliac Society, and the Center of Excellence in Microbial Food Safety (CoE-MiFoSa) research program of the Academy of Finland are gratefully acknowledged.
The authors declare that they have no competing interest.
Study concept and design: MK, WdV, RS; collection of clinical samples: MK, HL; acquisition of data: JC, MK, HH, HL, RS; analysis and interpretation of data: JC, MK, HH, JS, RS; statistical analysis: JC, JS; drafting of the manuscript: JC, MK, RS; critical reading and revision of the manuscript: JC, MK, AP, WdV, JS, RS; obtained funding: MK, AP, RS; administrative: MK, AP, RS; study supervision: MK, JS, RS. All authors read and approved the final version for publication.