Multi-omics analysis reveals the influence of genetic and environmental risk factors on developing gut microbiota in infants at risk of celiac disease

doi:10.21203/rs.2.24237/v1

Download PDF

Research

Multi-omics analysis reveals the influence of genetic and environmental risk factors on developing gut microbiota in infants at risk of celiac disease

https://doi.org/10.21203/rs.2.24237/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 11 Sep, 2020

Read the published version in Microbiome →

You are reading this older preprint version

Read the latest preprint version →

Background : Celiac disease (CD) is an autoimmune digestive disorder that occurs in genetically susceptible individuals in response to ingesting gluten, a protein found in wheat, rye, and barley. Research shows that genetic predisposition and exposure to gluten are necessary but not sufficient to trigger the development of CD. This suggests that exposure to other environmental stimuli early in life, e.g., cesarean section delivery, exposure to antibiotics or formula feeding, may also play a key role in CD pathogenesis through yet unknown mechanisms. Here, we use multi-omics analysis to investigate how genetic and early environmental risk factors alter the development of the gut microbiota in infants at risk of CD.

Results : Toward this end, we took advantage of a large-scale prospective birth cohort study of infants with a first-degree relative with CD. We then performed rigorous multivariate association, cross-sectional and longitudinal analyses using metagenomic and metabolomic data collected at birth, three months and six months of age to explore the impact of genetic predisposition and environmental risk factors on the gut microbiota composition, function and metabolome prior to the introduction of trigger (gluten). These analyses revealed several microbial species, functional pathways and metabolites that are associated with each genetic and environmental risk factor or that are differentially abundant between environmentally exposed and non-exposed infants or between time points. Among our significant findings, we found that cesarean section delivery is associated with a decreased abundance of Bacteroides vulgatus and Bacteroides dorei and of folate biosynthesis pathway, and with an increased abundance of hydroxyphenylacetic acid, alterations that are implicated in immune system dysfunction and inflammatory conditions. Additionally, longitudinal analysis revealed that, in infants not exposed to any environmental risk factor, the abundances of Bacteroides uniforms and of metabolite 3-3-hydroxyphenylproprionic acid increase over time, while those for lipoic acid and methane metabolism pathways decrease, patterns that are linked to beneficial immunomodulatory and anti-inflammatory effects.

Conclusions : Overall, our study provides unprecedented insights into major taxonomic and functional shifts in the developing gut microbiota of infants at risk of CD ultimately linking genetic and environmental risk factors to detrimental immunomodulatory and inflammatory effects.

Keywords: Microbiota; Celiac disease; Multi-omics analysis, gut microbiome

General Microbiology

Microbiota

Celiac disease

Multi-omics analysis

gut microbiome

Celiac disease (CD) is an autoimmune enteropathy, which affects three million Americans and 1% of the population worldwide [1]. CD occurs in genetically predisposed individuals that have specific variants of the human leukocyte antigen (HLA) DQ2 and DQ8 genes in response to ingesting gluten, a protein found in wheat, rye, and barley [2]. Notably, CD is the only autoimmune disorder for which the environmental trigger is known [3] (exposure to gluten). Given that the timing of exposure to gluten and the dose of gluten ingested can be carefully monitored, and since gluten removal results in the resolution of symptoms and enteropathy for most patients [4–8], CD can serve as a tunable model of chronic immune-based disorders [9]. This allows for insights into its pathogenesis to be applied not only to individuals with CD but those with other autoimmune diseases as well.

Globally, the incidence of autoimmune diseases including CD is expected to triple by 2050 [10, 11], yet the genes associated with CD (HLA DQ2) and DQ8, and the trigger (gluten) have not changed. Research shows that more than 30% of the population carry the predisposing gene(s) and are exposed to the trigger yet only 2–3% of these individuals develop CD in their lifetime thus suggesting a critical role for environmental factors [12]. Mode of delivery, infant feeding type, timing of gluten introduction into the diet, occurrence of viral infections and early exposure to antibiotics are just a few of the many environmental factors suggested to influence the development of chronic inflammatory diseases such as CD [13]. When evaluating these factors independently, case-control studies and meta-analyses have found that cesarean section delivery [14, 15], lack of breast-feeding [16, 17], timing of gluten introduction [17, 18] and exposure to antibiotics [19] increase the risk of developing CD. However, two independent double blind placebo controlled prospective studies in Europe involving infants with compatible HLA genetics and a first-degree relative with CD, and therefore, at high risk of developing CD, found that vaginal delivery, breast-feeding, and timing of gluten introduction were not protective against developing CD [20, 21].

Accumulating evidence suggests that the gut microbiota may be involved in several immune based disorders [13] such as inflammatory bowel disease [22], type 1 diabetes(T1D) [23] and multiple sclerosis [24]. A limited number of studies have also started to explore the link between the gut microbiota and CD development [25–30]. Initial studies focused on the contribution of HLA genetics to the developing microbiota. In particular, two studies analyzed exclusively breastmilk-fed infants up to 4 months of age with a first-degree relative with CD and found that Bacteroides-Prevotella group [25], Firmicutes, Proteobacteria, and Bifidobacterium [26] were more abundant in infants at high genetic risk for CD (those with two copies of HLA DQ2). Additionally, in a preliminary prospective study, we used 16S rRNA amplicon sequencing to examine the microbiota from 16 infants with a first-degree relative with CD and with a compatible HLA type and found a lower abundance of Bacteroides and a higher abundance of Firmicutes in these subjects compared to controls [27]. Other studies of the gut microbiota and CD have assessed changes, within one year of age, in the microbiota composition of individuals who later developed CD compared to controls [29, 30]. For example, Olivares et al [29] identified increases in the abundances of Firmicutes, Enterococcaceae, and Peptostreptococcaceae in controls from 4 to 6 months but no differences over time were observed in cases [29]. While the link between environmental factors and alterations in the gut microbiota of at-risk subjects has been recently explored for a number of chronic immune based disorders [31, 32], studies addressing this question for CD are scarce [28]. The only study in this direction is the work of Pozo-Rubio et al [28], where they found associations between a limited number of pre-selected specific fecal microbial taxa in subjects at risk of CD and delivery mode, infant feeding type, antibiotic exposure and rotavirus vaccine administration [28].

While these studies have provided valuable insights into the development of the gut microbiota early in life in subjects at risk of CD, solid food has already been introduced into the infants’ diet in many of these studies without accounting for its impact on the microbiota. In addition, to the best of our knowledge, no microbiome-wide study of the effect of environmental risk factors for CD currently exists. More importantly, existing studies are primarily based on 16S rRNA amplicon sequencing, which are not capable of addressing how the functional characterization of the microbiota will affect CD onset. To mitigate these limitations, here, we utilize a large-scale prospective cohort study called the Celiac Disease Genomic, Environmental, Microbiome and Metabolome study (CDGEMM) [33], where we have been following over 400 infants with a first-degree relative with CD who are thus at a high risk of developing CD. In this study, we present multivariate association as well as inter-subject and intra-subject analyses using metagenomic and metabolomic data collected over the first six months after birth to investigate the impact of both genetic and environmental risk factors on the developing gut microbiota of infants at risk of CD prior to the introduction of solid foods.

We selected 31 children recruited into the CDGEMM cohort for whom stool samples were available at birth, 3 months, and 4–6 months for this study (see Table 1 and see Additional file 1: Table S1 for more detailed metadata). None of these infants consumed solid foods before 6 months, which makes them ideal for studying the effect of genetic and environmental risk factors on the gut microbiota in the absence of gluten as a confounder. Twenty-six of these infants were genetically susceptible to developing CD out of which 19 were either heterozygous for DQ2 or DQ8 or carried both DQ2 and DQ8 (referred to as “standard genetic risk” hereafter) and seven were homozygous for DQ2 (referred to as “high genetic risk” hereafter). Additionally, 19 infants who were genetically predisposed to CD and that have been exposed to at least one environmental risk factor are referred to as “environmentally exposed” infants throughout the rest of manuscript. This means that these infants were born via cesarean section or exposed to antibiotics at or during birth (i.e., antibiotics administered to the mother during delivery) or were not exclusively breastmilk-fed (i.e., formula-fed or both formula- and breastmilk-fed). Seven infants who were genetically susceptible and that were not exposed to any of these environmental risk factors, i.e., those born vaginally and not exposed to antibiotics at or during delivery and exclusively breastmilk-fed, are referred to as “environmentally non-exposed” hereafter.

Table 1

Study cohort metadata and genotype. This study cohort was extracted from the larger CDGEMM prospective longitudinal birth cohort study [33].
	USA (n = 18)	Italy (n = 13)	Total (n = 31)
Gender (%)
Male	11 (61.1)	7 (53.8)	18 (58.0)
Female	7 (38.9)	6 (46.2)	13 (42.0)
Mode of Delivery (%)
Vaginal	11 (61.1)	7 (53.8)	18 (58.0)
C-section	7 (38.9)	6 (46.2)	13 (42.0)
Feeding Type (4–6 months of age) (%)
Breastmilk only	12 (66.7)	4 (30.7)	16 (51.6)
Formula only	5 (27.8)	6 (46.2)	11 (35.5)
Both	1 (5.5)	3 (23.1)	4 (12.9)
Antibiotic Exposure (%)
At delivery (mother)	7 (38.9)	2 (15.4)	9 (29.0)
At birth (infant)	2 (11.1)	2 (15.4)	4 (12.9)
Before 6 months of age (infant)	0 (0.0)	4 (30.8)	4 (12.9)
Genotype (%)
DQ2 Homozygous	6 (33.3)	1 (7.7)	7 (22.6)
DQ2 Heterozygous	6 (33.3)	6 (46.2)	12 (38.7)
DQ2/DQ8	3 (16.7)	2 (15.4)	5 (16.1)
DQ8	1 (5.5)	1 (7.7)	2 (6.5)
Negative	2 (11.1)	3 (23.1)	5 (16.1)
Relative with CD
Mother	15 (83.3)	7 (53.8)	22 (70.9)
Father	1 (5.5)	1 (7.7)	2 (6.5)
Sibling	2 (11.1)	5 (38.4)	7 (22.6)

Collected stool samples underwent shotgun metagenomic sequencing and metabolomic profiling. We analyzed metagenomic sequencing reads (see Methods) to profile microbial taxa at species-level resolution (see Additional file 2: S2) and functional pathways encoded by metagenomes (see Additional file3: S3). Additionally, stool samples underwent metabolomic profiling and were analyzed to identify metabolites present in each stool sample (see Additional file 4: S4). The identified microbial taxa, functional pathways and metabolites were then analyzed rigorously to explore how genetic and environmental risk factors influence the development of the gut microbiota as outlined below.

Associations between genetic and environmental risk factors and various microbiota features

We used the MaAslin procedure [34] to investigate how various microbiome features including microbial species, functional pathways and metabolites at each time point are associated with genetic risk for developing CD and three key environmental risk factors including mode of delivery, exposure to antibiotics and infant feeding type (see Fig. 1).

Genetic risk

We found that both high and standard genetic risk to develop CD are associated with a decreased abundance of several species of Streptococcus and Coprococcus at 4–6 months of age compared to those lacking genetic compatibility (Fig. 1A; FDR-adjusted p-value < 0.05). Notably, a decreased abundance of Coprococcus has been previously reported in the gut of individuals who carry a genetic risk to develop autoimmune conditions including CD [35]. Standard and high genetic risk for developing CD are also associated with an increased abundance of Bacteroides and Enterococcus species, respectively, at enrollment compared to no genetic risk. These observations are in agreement with previous studies [25, 26], however, an association between genetic risk and increased abundance of Bifidobacterium or Proteobacteria, which were reported before [25, 26] were not observed here. Among other significant associations, we found a decreased abundance of Veillonella, Parabacteroides and Clostridium perfringens at 4–6 months after birth in infants with high and standard genetic compatibility. This observation is contrary to case-control studies that report an increased abundance of these microbes in autoimmune conditions such as autoimmune liver disease [36], Bechet’s disease [37] and neuromyelitis optica [38].

In addition to association with microbial species, we found that a high genetic risk of developing CD is associated with a decreased abundance of a number of functional pathways at 4–6 months of age (Fig. 1B; p-value < 0.05). These pathways include amino acid metabolism, biosynthesis of secondary metabolites and metabolism of cofactors including ubiquinone and other terpenoid-quinone biosynthesis. Furthermore, we identified an association between high genetic risk and a number of metabolites, e.g., an increased abundance of butanoic acid and a decreased abundance of dihydroxyacteone at 3 and 4–6 months of age (Fig. 1C, p-value < 0.05).

Mode of delivery

We found that cesarean section delivery is associated with a decreased abundance of several species of Bacteroides and Parabacteroides at all time points and with an increased abundance of Enterococcus faecalis (at 3 months after birth) compared to vaginal delivery (Fig. 1A; FDR-adjusted p-value < 0.05) in agreement with previous work [23, 39–41]. For example, we found associations between cesarean section delivery and a decreased abundance of beneficial species Bacteroides vulgatus and Bacteroides dorei. An increased abundance of these species has been reported to lead to a decreased gut microbial production of lipopolysaccharide, which will improve immune function through mechanisms such as major histocompatibility production and T cell activation, among others [42]. Analysis of pathways shows also an association between cesarean section delivery and decreased riboflavin metabolism and folate biosynthesis at 4–6 months after birth and an increase in the abundance of glycerolipid metabolism at 3 and 4–6 months (Fig. 1B; p-value < 0.05). Of note, defects in folate biosynthesis have been linked to an impaired immune response to viral infections and reduced natural killer cell response possibly contributing to T1D onset [43]. Finally, metabolites analysis unveiled an association between cesarean section delivery and an increase in the abundance of a number of metabolites such as butanoic acid (at 3 and 4–6 months), glycolic acid, oxalic acid, and hydroxyphenlacetic acid (at 4–6 months) and a decrease in that of valine, serine, and arabinoic acid among others (at 4–6 months) (Fig. 1C, p-value < 0.05). An increased abundance of hydroxyphenlacetic acid in the serum has been associated with ulcerative colitis in a previous study [44], however, no clear links between the level of metabolites in the gut and those in the serum have been established yet. Additionally, serine, which is decreased in cesarean section delivery, has been reported to be required for effector T cell expansion and thus for modulating the adaptive immune response [45].

Infant feeding type

We examined three infant feeding types in this study including exclusive breastmilk feeding, exclusive formula feeding and both breastmilk and formula feeding, the last two of which were considered as environmental risk factors. Previous work shows an association between infant feeding type and distinct species of Bifidobacterium [23, 46]. Consistent with these reports, we observed that exposure to both breastmilk and formula is associated with a decreased abundance of Bifidobacterium breve (at 4–6 months) while exclusive formula feeding is associated with an increased abundance of Bifidobacterium adolescentis compared to exclusive breastmilk feeding (Fig. 1A; FDR-adjusted p-value < 0.05). We also found that exclusive formula feeding is associated with a decreased abundance of Staphylococcus epidermis (at enrollment) consistent with previous work [47], and with an increased abundance of Ruminococcus gnavus and Lachnospiraceae bacterium (at 3 and 4–6 months), which have been linked to allergic disease [48], diabetes [49] and colonic inflammation [50]. Pathway analysis shows that exposure to formula only or both breastmilk and formula is associated with an increased abundance of pathways for lipids, amino acids and terpendoids metabolism and xenobiotic degradation, and with a decreased abundance of pathways for carbohydrate and energy metabolism (Fig. 1B; p-value < 0.05). Additionally, metabolomic analysis uncovered an association between both breastmilk and formula feeding with a decreased abundance of homoserine, alpha-D-glucopyranoside, and hydrocinnamic acid (at 4–6 months) (Fig. 1C; p-value < 0.05). Exclusive formula feeding is also associated with an increase in sucrose and threonine and a decrease in oxalic acid and dihydroxyacetone abundances, among others (at 4–6 months).

Antibiotic use

We found an association between antibiotic exposure (as an environmental risk factor) and an increased abundance of Bacteroides thetaiotaomicron (at 4–6 months of age) (Fig. 1A; FDR-adjusted p-value < 0.05). This is corroborated with previous work suggesting that this species, which is an important metabolizer of polysaccharides, increases in abundance in response to amoxicillin exposure [51]. Other identified associations for antibiotic exposure not previously reported include an increased Propionibacterium, Subdoligranulum species and a decreased abundance of Bifidobacterium merycicum and Strepotococcus lutetiensis (at 4–6 months). Pathway analysis also revealed an association between antibiotic exposure and a decreased abundance of phenylalanine metabolism and an increased abundance of cyanoamino acid (3 and 4–6 months) and galactose metabolism (4–6 months) (Fig. 1B; p-value < 0.05). Analysis of metabolites showed associations between antibiotic exposure and a number of metabolites including decreased sucrose abundance (at 4–6 months) (Fig. 1C; p-value < 0.05).

Changes in the microbiota of environmentally exposed vs. non-exposed infants

Here, we performed a cross-sectional (inter-subject) analysis to explore how various features of the gut microbiota (microbes, pathways and metabolites) change between genetically predisposed infants who were exposed to at least one environmental risk factor noted before (environmentally exposed infants) vs. those who were not (environmentally non-exposed infants). This analysis did not identify any microbial species whose abundance is significantly different between the environmentally exposed and non-exposed infants (FDR-adjusted p-value < 0.05). Pathways analysis, however, revealed that environmentally exposed infants have a higher abundance of pathways for xenobiotic degradation, fatty acid metabolism, and lipid metabolism among others (at enrollment) and of pathways such as toluene and xylene and biphenyl degradation (at 4–6 months) (Fig. 2A; p-value < 0.05). Furthermore, metabolomic analysis identified a decrease in homoserine and 2-ketobutryic acid among others and an increase in ribose and tyramine at all study time points in environmentally exposed infants (Fig. 2B; p-value < 0.05).

Longitudinal changes in the microbiota of environmentally exposed and non-exposed infants

Given the unique prospective study design of our cohort, we were able to perform a longitudinal (intra-subject) analysis to gain additional insights beyond a cross-sectional analysis by identifying dynamic alterations in the gut microbiota composition, function and metabolome in the first six months after birth. To this end, we explored changes in the microbiota features noted above between all pairs of time points that are observed exclusively in environmentally exposed or exclusively in environmentally non-exposed infants (Fig. 3).

By longitudinal analysis of microbial species, we found that the abundance of a number of species increases over time in the environmentally exposed infants (Fig. 3A; FDR-adjusted p-value < 0.05). For example, the abundance of Anaerostipes caccae monotonically increases during the study period and that of Klebsiella species and Erysipelotrichaceae bacterium increases from enrollment to 4–6 months. Among these, Klebsiella, has been associated with the autoimmune condition ankylosing spondylitis [52]. When examining environmentally non-exposed infants, we observe that the abundance of Bacteroides uniformis monotonically increases during the first 6 months after birth, a pattern which has previously been reported in breastmilk-fed infants [53]. In addition, work in mice found that Bacteroides uniformis improves immune defense mechanisms, which are impaired in obesity, by decreasing TNF-α production and increasing IL-10 production [54]. In our study, we also observed a decrease in the abundance of Veillonella species from enrollment to 4–6 months in non-exposed infants. An increased abundance of Veillonella species has been associated with autoimmune hepatitis [36].

Longitudinal pathway analysis revealed that the abundance of ether lipid metabolism increases from 3 to 4–6 months of age in environmentally exposed infants (Fig. 3B; p-value < 0.05). Notably, a decreased abundance of ether lipids in the serum of children with T1D compared to healthy controls has been observed, [55] although the relationship between the abundance of microbial pathways for ether lipid metabolism in the gut and the level of ether lipids in the serum are yet to be explored. For the non-exposed infants, we observe a decrease in the abundance of sulfur metabolism and lipoic acid metabolism at 3 and 4–6 months, and of methane metabolism and biotin metabolism at 4–6 months compared to enrollment. These patterns are consistent with previous reports [34, 56–62]. For example, increased sulfur metabolism is associated with the development of type 1 diabetes [56] and linked to IBD [34]. Additionally, lipoic acid is an antioxidant that has been suggested to have beneficial immunomodulatory effects on the innate and adaptive immune systems in autoimmune diseases [57]. Methane has also been shown to have an anti-inflammatory effect, promoting immune tolerance in the intestine when tested in animal models [58, 59]. Furthermore, biotin is known to enhance innate [60] and adaptive immune responses [61] and biotin deficiency has been associated with immune disorders and inflammation[63, 64]. A previous study also found that high dose of biotin may be useful in treating multiple sclerosis [62].

Metabolomic analysis revealed a monotonic increase in erythritol abundances during the study period and a decrease in propionic acid abundance from enrollment to 4–6 months in environmentally exposed infants (Fig. 3C; p-value < 0.05). Propionic acid produced in the colon via bacterial fermentation of fiber promotes regulatory T cell generation [65]. Additionally, increased serum levels of erythritol have been associated with central obesity and weight gain [66], though the link between metabolite levels in the gut and those in the serum is not clear. In environmentally non-exposed infants, we observed an increased abundance of uracil, 3-3-hydroxyphenylpropionic acid and dihydroxyacetone from enrollment to 4–6 months. Previous work suggests that 3-hydroxyphenylproprionic acid acts as an anti-inflammatory and antioxidant agent [67].

Several studies have linked exposure to a variety of genetic and environmental risk factors with the onset of non-infective chronic inflammatory diseases. This link has been typically based on the results obtained from either clinical case-control studies [14, 15, 19, 68] or metanalyses [69–72] in which cause-effect relationship cannot always conclusively be determined. Since host genetics and environmental factors are known to influence the gut microbiota composition and function, researchers have started to explore alterations in the gut microbiota of infants at risk of autoimmune conditions such as inflammatory bowel disease [31] or T1D [32]. However, to date, there is no systematic study of how protective or detrimental genetic and environmental factors may change the gut microbiota engraftment and its maturation during the first months of life in CD. In an effort to fill this gap, in this study we used metagenomic and metabolomic data collected in the first 6 months after birth to associate individual risk factors (HLA DQ2/DQ8 genetics, cesarean section delivery, antibiotic use and partial or exclusive formula feeding) with microbial species, pathways, or metabolites in the gut. Additionally, we performed cross-sectional analysis to identify microbes, pathways, or metabolites that are differentially abundant between infants exposed to at least one environmental risk factor and infants who were not, as well as longitudinal analysis to identify dynamic changes in the gut microbiota in the first six months of life.

Many microbes, pathways, or metabolites that we identified in these analyses are well supported in the literature to be associated with inflammation, autoimmune disease, or immune system dysfunction, thereby suggesting that they may have similar effects in CD. For example, we found that high risk HLA genetics and formula feeding are both associated with an increased abundance of Ruminococcus gnavus and Lachnospiraceae bacterium, which are linked to allergic diseases [48] and diabetes [49], respectively. Among other significant findings are associations between cesarean section delivery and a decreased abundance of Bacteroides vulgatus and Bacteroides dorei and folate biosynthesis pathway along with an increased abundance of hydroxyphenlacetic acid. All of these patterns have been reported to be associated with impaired immune function [42] and inflammatory conditions such as T1D and ulcerative colitis [43, 44] suggesting that could also predispose infants to develop CD.

While the cross-sectional analysis did not identify any microbial species whose abundance significantly changes between the environmentally exposed and non-exposed infants at any given time point, our longitudinal analysis yielded significant results further stressing the power of intra-subject analysis. This allows us to prospectively evaluate the impact of risk factors on the dynamics of the gut microbiota development and to link dynamics to increased susceptibility to inflammation. For example, environmentally exposed infants show an increasing abundance over time of Klebsiella, a microbe linked to autoimmune disease [52] and a decreasing abundance over time of propionic acid, a metabolite that promotes innate and adaptive immunity[65]. In contrast, in infants not exposed to environmental risk factors, we observe patterns associated with beneficial immunomodulatory effects and protection against immune system activation and inflammation such as increasing abundance of Bacteroides uniformis over time and decreasing abundance of lipoic acid and methane metabolism [54, 57–59, 67]. Notably, during our analyses, we identified a number of metabolites and pathways with altered abundances in the gut (including hydroxyphenlacetic acid, erythritol and ether lipid metabolism) for which similar variations in the serum are reported to be associated with autoimmune conditions. While the importance of gut-blood axis has been realized fairly recently [73, 74], further investigations are needed to better understand the relationship between different features of the gut microbiota and host- or microbially-derived metabolites in the blood.

Unlike previous microbiome studies for CD that are often based on 16S rRNA amplicon sequencing, here, we use shotgun metagenomic sequencing, which is amenable to functional characterization of the microbiota. This is particularly important as previous studies have shown that functional characterization is a more robust descriptor of the status of the microbiota compared to taxonomic composition alone [75, 76]. Furthermore, unlike typical case-control studies, where disease symptoms have already emerged in cases, our prospective birth cohort provides the opportunity to mechanistically link major shifts in the gut microbiota early in life, due to genetic risk factors and environmental exposures, to increased susceptibility to CD. Nevertheless, our data should be considered exploratory given the relatively small sample size. This limitation can be mitigated through ongoing recruitment into our CDGEMM cohort, which will allow us to validate our findings using a much larger number of subjects in the future.

In this paper, we utilized an ongoing prospective study and multi-omics analysis to perform an in-depth analysis of the impact of genetic and environmental risk factors on the longitudinal development of the gut microbiota in infants at risk for CD, before solid foods (including the trigger of CD, gluten) is introduced. These analyses revealed several microbial species, functional pathways and metabolites that have been previously linked to inflammation or immune system dysfunction as well as several new ones that have not been reported before and could be specific to CD. In this study, we restricted our analysis to the first 6 months of life and particularly prior to the introduction of solid foods in order to proactively “regress out” the effect of gluten on the gut microbiota as a major confounder when analyzing the effect of genetic and environmental risk factors. Further work is required not only to investigate alterations in the gut microbiota over a longer period of time, but also to consider other environmental factors such as viral infections, timing of solid food (gluten) introduction, amount of gluten ingested and household exposures, e.g. family size and contact with pets, which have been reported to be associated with altered microbiomes [77], or with protection against autoimmune conditions such as asthma [78] and T1D [79]. These investigations warrant future studies, which can utilize this longitudinal study design and multi-omics analysis to connect alterations in the gut microbiota early in life to the loss of tolerance to gluten and the development of CD.

Subjects, sampling, and factors of interest

Thirty one healthy infants from the United States (n = 18) and Italy (n = 13) with a first-degree relative with CD participating in the CDGEMM prospective birth cohort study [33] were included in our analysis. These subjects consist of all infants from CDGEMM with available stool samples collected before the introduction of solid foods at 7–15 days (enrollment), 3 months and 4–6 months after birth. Parents answered a detailed questionnaire at enrollment that addressed pregnancy, delivery, family history, household factors, and many other factors related to the infants’ environment before birth and at delivery. Parents also filled out monthly diaries, which addressed infant food intake and any exposure to antibiotics. Infant feeding type was determined according to the reported exclusive feeding type for at least two of the three sample time point collections. Infants who received both breastmilk and formula for at least two of three sample collection points were classified as “both breastmilk and formula fed.” HLA genetic type was determined from whole blood at time of birth (cord blood) or 12 months of age using the DQ-CD Typing Plus (BioDiagne, Palermo, Italy) per the manufacturer’s instructions. Written informed consent was obtained from the parents of infants included in the study according to the standards outlined and approved by the Partners Human Research Committee Institutional Review Board.

DNA extraction

All fecal samples included in the metagenomic analysis were stored and processed centrally in the United States. Total DNA from each sample was extracted using the Qiagen Power soil DNA extraction kit (Qiagen, Hilden, Germany).

Taxonomic profiling

Metagenomic sequencing reads were analyzed by using the CosmosID’s (CosmosID Inc., Rockville, MD) commercial metagenomic analysis platform (formerly knowns as GENIUS; https://app.cosmosid.com/login) [80, 81] to reveal the underlying microbial community composition up to the species-level resolution (see Additional file 5 for a description of this platform).

Functional profiling

After trimming the raw sequencing reads using BBDuk (https://jgi.doe.gov/data-and-tools/bbtools/), we used the SPAdes tool [82] for the assembly of metagenomes and subsequently Prodigal [83] to identify protein coding sequences in the assembled metagenomes. We then utilized InterProScan [84] to annotate the identified genes with biochemical functions based on the KEGG pathways [85]. The relative abundance of each gene was computed as

where, G is fragments per kilobase per million (FPKM) for each gene, L is the length of the gene, C is the coverage of contig in which the gene is identified, R is the read length and K is the k-mer size [86]. The relative abundance of each KEGG pathway was then quantified by summing the relative abundances of all the genes associated to that pathway.

Metabolomic profiling

All stool samples for metabolomics were stored and processed in Italy. The metabolome extraction, purification and derivatization were carried by the MetaboPrep GC kit (Theoreo, Montecorvino Pugliano, Italy) according to manufacturer instructions. Instrumental analyses were performed with a GC-MS system (GC-2010 Plus gas chromatograph and QP2010 Plus mass spectrometer; Shimadzu Corp., Kyoto, Japan). Sample analysis was performed in triplicate. Additional information related to the extraction, purification, derivatization, GC-MS analysis, and data preprocessing can be found in the Supplementary Text (see Additional file 5). The molecular identity of metabolites was determined by analysis of the corresponding mass spectrum in the chromatogram, setting the linear index difference max tolerance to 10. These identified metabolites were further confirmed using external standards according to level 1 Metabolomics Standards Initiative (MSI) [87].

Identifying associations between genetic and environmental risk factors and microbiome features

We used the widely-used multivariate statistical framework, MaAsLin, [22] to determine associations between microbial species, functional pathways or metabolites and genetic and environmental risk factors including HLA genetics, delivery mode, infant feeding type and antibiotic exposure at each time point. No genetic risk, vaginal delivery, exclusive breastmilk-feeding and no antibiotic exposure were considered as the reference levels for HLA genetics, delivery mode, infant feeding type and antibiotic exposure, respectively. All metadata variables were forced simultaneously to control for confounders. We chose a false discovery rate-adjusted p-value of 0.01 for reporting significant microbial species and a p-value of 0.05 for functional pathways and metabolites.

Cross-sectional and longitudinal analysis

For the cross-sectional analysis, we performed the Mann-Whitney U (Wilcoxon Rank Sum) test to compare the abundance of microbial species, pathways and metabolites at each time point between the environmentally exposed and non-exposed groups (using a FDR-adjusted p-value threshold of 0.01 for microbial species and a p-value threshold of 0.05 for pathways and metabolite to report significant results). For the longitudinal analysis, we performed the paired Wilcoxon (Wilcoxon Signed Rank) test to compare the abundances of microbial species, pathways and metabolites between each pair of time points using the same (adjusted) p-value thresholds noted above to report the significant results. All analyses were performed in Python.

Celiac disease (CD); Human leukocyte antigen (HLA); Type 1 diabetes (T1D); False Discovery Rate (FDR).

Ethics approval and consent to participate:

Written informed consent was obtained from the parents of infants included in the study according the standards outlined and approved by the Partners Human Research Committee Institutional Review Board.

Clinical Trial Registration

This study is registered at clinicaltrials.gov with the identifier NCT02061306.

Availability of data and material

The datasets supporting the conclusions of this article are submitted to the NCBI Short Read Archive (SRA) repository, under BioProjectID: PRJNA486782 and SRA accession number SRP158417. Additional data from the analyses presented in this paper are available in the Supplementary Material.

Competing interests

AF is a stockholder at Alba Therapeutics, serves as a consultant for Inova Diagnostics and Innovate Biopharmaceuticals, is an advisory board member for Axial Biotherapeutics and Ubiome, and has a speaker agreement with Mead Johnson Nutrition. MML serves as a consultant to HealthMode and Anokion, has a speaker agreement with Takeda Pharmaceuticals, and performs sponsored research with Glutenostics LLC. HK is a former employee, BF is a current employee, PS is a consultant and RRC and NAH are stockholders at CosmosID Inc. Other authors have declared no competing interests exist.

Authors' contributions

MML designed the study, analyzed the data analysis results and drafted the manuscript. HK performed taxonomic and functional analysis. PS contributed to taxonomic profiling and analysis. MP, JT, and AC designed metabolomic studies and performed metabolomic data analysis. GS isolated the DNA. SC and VK recruited participants and coordinated the US study. FV recruited and supervised the Italian sites and secured sample collection. SC, MM, BM, RF, and LE recruited and coordinated the study, secured sample collection, and subject participation. NAH and RRC oversee CosmosID metagenomic sequencing and data analysis platforms. BF and NAH performed taxonomic profiling of metagenomic data analysis. ARZ designed the computational studies, analyzed the data analysis results and drafted the manuscript. AF conceived the study, analyzed the data analysis results and aided in manuscript preparation. All authors read and approved the manuscript.

Acknowledgements:

The authors would like to thank the families that participate in this study and whose contribution was instrumental to the findings described in this manuscript and the CDGEMM team including Pasqua Piemontese, Angela Calvi, Mariella Baldassarre, Lorenzo Norsa, Chiara Maria Trovato, Celeste Lidia Raguseo, Tiziana Passaro, Paola Roggero, Marco Crocco, Annalisa Morelli, Michela Perrone, Marcello Chieppa, Giovanni Scala,Maria Elena Lionetti, Carlo Catassi, Adelaide Serretiello, Corrado Vecchi, and Gemma Castillejo de Villsante.

Funding

This work was partially supported by funding from the NIH NIDDK; DK104344 to AF, DK109620 and K23DK122127 to MML, funding from Nutrition Obesity Research Center at Harvard (P30-DK040561) and the Thrasher Research Fund to MML, and the faculty start-up funding by Mucosal Immunology and Biology Research Center at Massachusetts General Hospital to ARZ, and through the generous support of Joyce and Hugh McCormick.

Additional file 1. Table S1. Clinical metadata for the subjects in this study.

Additional file 2. Table S2. Results of the taxonomic profiling of metagenomic samples at each time point.

Additional file 3. Table S3. Results of the functional profiling of metagenomic samples at each time point.

Additional file 4. Table S4. Results of the metabolomic profiling of stool samples at each time point.

Additional file 5. Text S5. Supplementary text describing details of data analysis methods.

Lionetti E, Gatti S, Pulvirenti A, Catassi C: Celiac disease from a global perspective. Best Pract Res Clin Gastroenterol 2015, 29(3):365-379.
Schuppan D: Current concepts of celiac disease pathogenesis. Gastroenterology 2000, 119(1):234-242.
Green PH, Cellier C: Celiac disease. N Engl J Med 2007, 357(17):1731-1743.
Vecsei E, Steinwendner S, Kogler H, Innerhofer A, Hammer K, Haas OA, Amann G, Chott A, Vogelsang H, Schoenlechner R et al: Follow-up of pediatric celiac disease: value of antibodies in predicting mucosal healing, a prospective cohort study. BMC Gastroenterol 2014, 14:28.
Leonard MM, Weir DC, DeGroote M, Mitchell PD, Singh P, Silvester JA, Leichtner AM, Fasano A: Value of IgA tTG in Predicting Mucosal Recovery in Children With Celiac Disease on a Gluten-Free Diet. J Pediatr Gastroenterol Nutr 2017, 64(2):286-291.
Ciacci C, Cirillo M, Cavallaro R, Mazzacca G: Long-term follow-up of celiac adults on gluten-free diet: prevalence and correlates of intestinal damage. Digestion 2002, 66(3):178-185.
Rubio-Tapia A, Rahim MW, See JA, Lahr BD, Wu TT, Murray JA: Mucosal recovery and mortality in adults with celiac disease after treatment with a gluten-free diet. Am J Gastroenterol 2010, 105(6):1412-1420.
Valitutti F, Trovato CM, Montuori M, Cucchiara S: Pediatric Celiac Disease: Follow-Up in the Spotlight. Adv Nutr 2017, 8(2):356-361.
Valitutti F, Fasano A: Breaking Down Barriers: How Understanding Celiac Disease Pathogenesis Informed the Development of Novel Treatments. Dig Dis Sci 2019, 64(7):1748-1758.
West J, Fleming KM, Tata LJ, Card TR, Crooks CJ: Incidence and prevalence of celiac disease and dermatitis herpetiformis in the UK over two decades: population-based study. The American journal of gastroenterology 2014, 109(5):757.
Catassi C, Kryszak D, Bhatti B, Sturgeon C, Helzlsouer K, Clipp SL, Gelfond D, Puppa E, Sferruzza A, Fasano A: Natural history of celiac disease autoimmunity in a USA cohort followed since 1974. Ann Med 2010, 42(7):530-538.
Ricano-Ponce I, Wijmenga C, Gutierrez-Achury J: Genetics of celiac disease. Best Pract Res Clin Gastroenterol 2015, 29(3):399-412.
Tamburini S, Shen N, Wu HC, Clemente JC: The microbiome in early life: implications for health outcomes. Nat Med 2016, 22(7):713-722.
Decker E, Engelmann G, Findeisen A, Gerner P, Laass M, Ney D, Posovszky C, Hoy L, Hornef MW: Cesarean delivery is associated with celiac disease but not inflammatory bowel disease in children. Pediatrics 2010, 125(6):e1433-1440.
Marild K, Stephansson O, Montgomery S, Murray JA, Ludvigsson JF: Pregnancy outcome and risk of celiac disease in offspring: a nationwide case-control study. Gastroenterology 2012, 142(1):39-45 e33.
Akobeng AK, Ramanan AV, Buchan I, Heller RF: Effect of breast feeding on risk of coeliac disease: a systematic review and meta-analysis of observational studies. Arch Dis Child 2006, 91(1):39-43.
Szajewska H, Chmielewska A, Piescik-Lech M, Ivarsson A, Kolacek S, Koletzko S, Mearin ML, Shamir R, Auricchio R, Troncone R et al: Systematic review: early infant feeding and the prevention of coeliac disease. Aliment Pharmacol Ther 2012, 36(7):607-618.
Norris JM, Barriga K, Hoffenberg EJ, Taki I, Miao D, Haas JE, Emery LM, Sokol RJ, Erlich HA, Eisenbarth GS et al: Risk of celiac disease autoimmunity and timing of gluten introduction in the diet of infants at increased risk of disease. JAMA 2005, 293(19):2343-2351.
Marild K, Ye W, Lebwohl B, Green PH, Blaser MJ, Card T, Ludvigsson JF: Antibiotic exposure and the development of coeliac disease: a nationwide case-control study. BMC Gastroenterol 2013, 13:109.
Lionetti E, Castellaneta S, Francavilla R, Pulvirenti A, Tonutti E, Amarri S, Barbato M, Barbera C, Barera G, Bellantoni A et al: Introduction of gluten, HLA status, and the risk of celiac disease in children. N Engl J Med 2014, 371(14):1295-1303.
Vriezinga SL, Auricchio R, Bravi E, Castillejo G, Chmielewska A, Crespo Escobar P, Kolacek S, Koletzko S, Korponay-Szabo IR, Mummert E et al: Randomized feeding intervention in infants at high risk for celiac disease. N Engl J Med 2014, 371(14):1304-1315.
Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, Reyes JA, Shah SA, LeLeiko N, Snapper SB et al: Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol 2012, 13(9):R79.
Stewart CJ, Ajami NJ, O’Brien JL, Hutchinson DS, Smith DP, Wong MC, Ross MC, Lloyd RE, Doddapaneni H, Metcalf GA: Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature 2018, 562(7728):583.
Jangi S, Gandhi R, Cox LM, Li N, Von Glehn F, Yan R, Patel B, Mazzola MA, Liu S, Glanz BL: Alterations of the human gut microbiome in multiple sclerosis. Nature communications 2016, 7:12015.
De Palma G, Capilla A, Nadal I, Nova E, Pozo T, Varea V, Polanco I, Castillejo G, Lopez A, Garrote JA et al: Interplay between human leukocyte antigen genes and the microbial colonization process of the newborn intestine. Curr Issues Mol Biol 2010, 12(1):1-10.
Olivares M, Neef A, Castillejo G, Palma GD, Varea V, Capilla A, Palau F, Nova E, Marcos A, Polanco I et al: The HLA-DQ2 genotype selects for early intestinal microbiota composition in infants at high risk of developing coeliac disease. Gut 2015, 64(3):406-417.
Sellitto M, Bai G, Serena G, Fricke WF, Sturgeon C, Gajer P, White JR, Koenig SS, Sakamoto J, Boothe D et al: Proof of concept of microbiome-metabolome analysis and delayed gluten exposure on celiac disease autoimmunity in genetically at-risk infants. PLoS One 2012, 7(3):e33387.
Pozo-Rubio T, de Palma G, Mujico JR, Olivares M, Marcos A, Acuna MD, Polanco I, Sanz Y, Nova E: Influence of early environmental factors on lymphocyte subsets and gut microbiota in infants at risk of celiac disease; the PROFICEL study. Nutr Hosp 2013, 28(2):464-473.
Olivares M, Walker AW, Capilla A, Benitez-Paez A, Palau F, Parkhill J, Castillejo G, Sanz Y: Gut microbiota trajectory in early life may predict development of celiac disease. Microbiome 2018, 6(1):36.
Rintala A, Riikonen I, Toivonen A, Pietila S, Munukka E, Pursiheimo JP, Elo LL, Arikoski P, Luopajarvi K, Schwab U et al: Early fecal microbiota composition in children who later develop celiac disease and associated autoimmunity. Scand J Gastroenterol 2018:1-7.
Torres J, Hu J, Seki A, Eisele C, Nair N, Huang R, Tarassishin L, Jharap B, Cote-Daigneault J, Mao Q et al: Infants born to mothers with IBD present with altered gut microbiome that transfers abnormalities of the adaptive immune system to germ-free mice. Gut 2020, 69(1):42-51.
Vatanen T, Franzosa EA, Schwager R, Tripathi S, Arthur TD, Vehik K, Lernmark Å, Hagopian WA, Rewers MJ, She JX et al: The human gut microbiome in early-onset type 1 diabetes from the TEDDY study. Nature 2018, 562(7728):589-594.
Leonard MM, Camhi S, Huedo-Medina TB, Fasano A: Celiac Disease Genomic, Environmental, Microbiome, and Metabolomic (CDGEMM) Study Design: Approach to the Future of Personalized Prevention of Celiac Disease. Nutrients 2015, 7(11):9325-9336.
Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, Reyes JA, Shah SA, LeLeiko N, Snapper SB: Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome biology 2012, 13(9):R79.
Hov JR, Zhong H, Qin B, Anmarkrud JA, Holm K, Franke A, Lie BA, Karlsen TH: The influence of the autoimmunity-associated ancestral HLA haplotype AH8. 1 on the human gut microbiota: a cross-sectional study. PLoS One 2015, 10(7):e0133804.
Wei Y, Li Y, Yan L, Sun C, Miao Q, Wang Q, Xiao X, Lian M, Li B, Chen Y et al: Alterations of gut microbiome in autoimmune hepatitis. Gut 2019.
Ye Z, Zhang N, Wu C, Zhang X, Wang Q, Huang X, Du L, Cao Q, Tang J, Zhou C et al: A metagenomic study of the gut microbiome in Behcet's disease. Microbiome 2018, 6(1):135.
Cree BA, Spencer CM, Varrin-Doyer M, Baranzini SE, Zamvil SS: Gut microbiome analysis in neuromyelitis optica reveals overabundance of Clostridium perfringens. Ann Neurol 2016, 80(3):443-447.
Shao Y, Forster SC, Tsaliki E, Vervier K, Strang A, Simpson N, Kumar N, Stares MD, Rodger A, Brocklehurst P: Stunted microbiota and opportunistic pathogen colonization in caesarean-section birth. Nature 2019:1-5.
Bokulich NA, Chung J, Battaglia T, Henderson N, Jay M, Li H, Lieber AD, Wu F, Perez-Perez GI, Chen Y: Antibiotics, birth mode, and diet shape microbiome maturation during early life. Science translational medicine 2016, 8(343):343ra382-343ra382.
Wampach L, Heintz-Buschart A, Fritz JV, Ramiro-Garcia J, Habier J, Herold M, Narayanasamy S, Kaysen A, Hogan AH, Bindl L: Birth mode is associated with earliest strain-conferred gut microbiome functions and immunostimulatory potential. Nature communications 2018, 9(1):5091.
Yoshida N, Emoto T, Yamashita T, Watanabe H, Hayashi T, Tabata T, Hoshi N, Hatano N, Ozawa G, Sasaki N: Bacteroides vulgatus and Bacteroides dorei reduce gut microbial lipopolysaccharide production and inhibit atherosclerosis. Circulation 2018, 138(22):2486-2498.
Bayer AL, Fraker CA: The Folate Cycle As a Cause of Natural Killer Cell Dysfunction and Viral Etiology in Type 1 Diabetes. Front Endocrinol (Lausanne) 2017, 8:315.
Sitkin SI, Tkachenko EI, Vakhitov T, Oreshko LS, Zhigalova TN: [Serum metabolome by gas chromatography-mass spectrometry (GC-MS) in patients with ulcerative colitis and celiac disease]. Eksp Klin Gastroenterol 2013(12):44-57.
Ma EH, Bantug G, Griss T, Condotta S, Johnson RM, Samborska B, Mainolfi N, Suri V, Guak H, Balmer ML et al: Serine Is an Essential Metabolite for Effector T Cell Expansion. Cell Metab 2017, 25(2):482.
Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, Li Y, Xia Y, Xie H, Zhong H: Dynamics and stabilization of the human gut microbiome during the first year of life. Cell host & microbe 2015, 17(5):690-703.
Lundequist B, Nord CE, Winberg J: The composition of the faecal microflora in breastfed and bottle fed infants from birth to eight weeks. Acta Paediatr Scand 1985, 74(1):45-51.
Chua HH, Chou HC, Tung YL, Chiang BL, Liao CC, Liu HH, Ni YH: Intestinal Dysbiosis Featuring Abundance of Ruminococcus gnavus Associates With Allergic Diseases in Infants. Gastroenterology 2018, 154(1):154-167.
Kameyama K, Itoh K: Intestinal colonization by a Lachnospiraceae bacterium contributes to the development of diabetes in obese mice. Microbes and environments 2014:ME14054.
Zeng H, Ishaq SL, Zhao F-Q, Wright A-DG: Colonic inflammation accompanies an increase of β-catenin signaling and Lachnospiraceae/Streptococcaceae bacteria in the hind gut of high-fat diet-fed mice. The Journal of nutritional biochemistry 2016, 35:30-36.
Cabral DJ, Penumutchu S, Reinhart EM, Zhang C, Korry BJ, Wurster JI, Nilson R, Guang A, Sano WH, Rowan-Nash AD: Microbial metabolism modulates antibiotic susceptibility within the murine gut microbiome. Cell metabolism 2019, 30(4):800-823. e807.
Wilson C, Tiwana H, Ebringer A: Molecular mimicry between HLA-DR alleles associated with rheumatoid arthritis and Proteus mirabilis as the Aetiological basis for autoimmunity. Microbes Infect 2000, 2(12):1489-1496.
Sanchez E, De Palma G, Capilla A, Nova E, Pozo T, Castillejo G, Varea V, Marcos A, Garrote JA, Polanco I et al: Influence of environmental and genetic factors linked to celiac disease risk on infant gut colonization by Bacteroides species. Appl Environ Microbiol 2011, 77(15):5316-5323.
Cano PG, Santacruz A, Moya Á, Sanz Y: Bacteroides uniformis CECT 7771 ameliorates metabolic and immunological dysfunction in mice with high-fat-diet induced obesity. PloS one 2012, 7(7):e41079.
Orešič M, Simell S, Sysi-Aho M, Näntö-Salonen K, Seppänen-Laakso T, Parikka V, Katajamaa M, Hekkala A, Mattila I, Keskinen P: Dysregulation of lipid and amino acid metabolism precedes islet autoimmunity in children who later progress to type 1 diabetes. Journal of Experimental Medicine 2008, 205(13):2975-2984.
Brown CT, Davis-Richardson AG, Giongo A, Gano KA, Crabb DB, Mukherjee N, Casella G, Drew JC, Ilonen J, Knip M: Gut microbiome metagenomics analysis suggests a functional model for the development of autoimmunity for type 1 diabetes. PloS one 2011, 6(10):e25792.
Liu W, Shi L-j, Li S-g: The Immunomodulatory Effect of Alpha-Lipoic Acid in Autoimmune Diseases. BioMed research international 2019, 2019.
Boros M, Ghyczy M, Érces D, Varga G, Tokés T, Kupai K, Torday C, Kaszaki J: The anti-inflammatory effects of methane. Critical care medicine 2012, 40(4):1269-1278.
Zhang X, Li N, Shao H, Meng Y, Wang L, Wu Q, Yao Y, Li J, Bian J, Zhang Y: Methane limit LPS-induced NF-κB/MAPKs signal in macrophages and suppress immune response in mice by enhancing PI3K/AKT/GSK-3β-mediated IL-10 expression. Scientific reports 2016, 6:29359.
Agrawal S, Agrawal A, Said HM: Biotin deficiency enhances the inflammatory response of human dendritic cells. American Journal of Physiology-Cell Physiology 2016, 311(3):C386-C391.
Kung JT, Mackenzie CG, Talmage DW: The requirement for biotin and fatty acids in the cytotoxic T-cell response. Cellular immunology 1979, 48(1):100-110.
Sedel F, Bernard D, Mock DM, Tourbah A: Targeting demyelination and virtual hypoxia with high-dose biotin as a treatment for progressive multiple sclerosis. Neuropharmacology 2016, 110:644-653.
Abad-Lacruz A, Fernandez-Banares F, Cabre E, Gil A, Esteve M, Gonzalez-Huix F, Xiol X, Gassull M: The effect of total enteral tube feeding on the vitamin status of malnourished patients with inflammatory bowel disease. International journal for vitamin and nutrition research Internationale Zeitschrift fur Vitamin-und Ernahrungsforschung Journal international de vitaminologie et de nutrition 1988, 58(4):428-435.
Fernandez-Banares F, Abad-Lacruz A, Xiol X, Gine J, Dolz C, Cabre E, Esteve M, Gonzalez-Huix F, Gassull M: Vitamin status in patients with inflammatory bowel disease. American Journal of Gastroenterology 1989, 84(7).
Arpaia N, Campbell C, Fan X, Dikiy S, van der Veeken J, deRoos P, Liu H, Cross JR, Pfeffer K, Coffer PJ et al: Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature 2013, 504(7480):451-455.
Hootman KC, Trezzi J-P, Kraemer L, Burwell LS, Dong X, Guertin KA, Jaeger C, Stover PJ, Hiller K, Cassano PA: Erythritol is a pentose-phosphate pathway metabolite and associated with adiposity gain in young adults. Proceedings of the National Academy of Sciences 2017, 114(21):E4233-E4240.
Fan FY, Sang LX, Jiang M: Catechins and Their Therapeutic Benefits to Inflammatory Bowel Disease. Molecules 2017, 22(3).
Baron S, Turck D, Leplat C, Merle V, Gower-Rousseau C, Marti R, Yzet T, Lerebours E, Dupas JL, Debeugny S et al: Environmental risk factors in paediatric inflammatory bowel diseases: a population based case control study. Gut 2005, 54(3):357-363.
Xu L, Lochhead P, Ko Y, Claggett B, Leong RW, Ananthakrishnan AN: Systematic review with meta-analysis: breastfeeding and the risk of Crohn's disease and ulcerative colitis. Aliment Pharmacol Ther 2017, 46(9):780-789.
Ungaro R, Bernstein CN, Gearry R, Hviid A, Kolho KL, Kronman MP, Shaw S, Van Kruiningen H, Colombel JF, Atreja A: Antibiotics associated with increased risk of new-onset Crohn's disease but not ulcerative colitis: a meta-analysis. Am J Gastroenterol 2014, 109(11):1728-1738.
Costenbader KH, Kim DJ, Peerzada J, Lockman S, Nobles-Knight D, Petri M, Karlson EW: Cigarette smoking and the risk of systemic lupus erythematosus: a meta-analysis. Arthritis Rheum 2004, 50(3):849-857.
McCormic ZD, Khuder SS, Aryal BK, Ames AL, Khuder SA: Occupational silica exposure as a risk factor for scleroderma: a meta-analysis. Int Arch Occup Environ Health 2010, 83(7):763-769.
Wilmanski T, Rappaport N, Earls JC, Magis AT, Manor O, Lovejoy J, Omenn GS, Hood L, Gibbons SM, Price ND: Blood metabolome predicts gut microbiome α-diversity in humans. Nat Biotechnol 2019, 37(10):1217-1228.
Analysis of blood and fecal microbiome profile in patients with celiac disease. Human <strong data-auto="strong_text" style="margin: 0px; padding: 0px; border: 0px; font-size: 16200000762939453px; vertical-align: baseline; caret-color: rgb(51, 51, 51); color: rgb(51, 51, 51); font-family: Helvetica, Arial, sans-serif;" xmlns:translation="urn:EBSCO-Translation">Microbiome Journal 2019, 11.
Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, Rodriguez-Mueller B, Zucker J, Thiagarajan M, Henrissat B et al: Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol 2012, 8(6):e1002358.
Human Microbiome Project C: Structure, function and diversity of the healthy human microbiome. Nature 2012, 486(7402):207-214.
Sjogren YM, Jenmalm MC, Bottcher MF, Bjorksten B, Sverremark-Ekstrom E: Altered early infant gut microbiota in children developing allergy up to 5 years of age. Clin Exp Allergy 2009, 39(4):518-526.
Ownby DR, Johnson CC, Peterson EL: Exposure to dogs and cats in the first year of life and risk of allergic sensitization at 6 to 7 years of age. JAMA 2002, 288(8):963-972.
Virtanen SM, Takkinen HM, Nwaru BI, Kaila M, Ahonen S, Nevalainen J, Niinisto S, Siljander H, Simell O, Ilonen J et al: Microbial exposure in infancy and subsequent appearance of type 1 diabetes mellitus-associated autoantibodies: a cohort study. JAMA Pediatr 2014, 168(8):755-763.
Hasan NA, Young BA, Minard-Smith AT, Saeed K, Li H, Heizer EM, McMillan NJ, Isom R, Abdullah AS, Bornman DM et al: Microbial community profiling of human saliva using shotgun metagenomic sequencing. PLoS One 2014, 9(5):e97699.
Ponnusamy D, Kozlova EV, Sha J, Erova TE, Azar SR, Fitts EC, Kirtley ML, Tiner BL, Andersson JA, Grim CJ et al: Cross-talk among flesh-eating Aeromonas hydrophila strains in mixed infection leading to necrotizing fasciitis. Proc Natl Acad Sci U S A 2016, 113(3):722-727.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD: SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. Journal of computational biology 2012, 19(5):455-477.
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ: Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC bioinformatics 2010, 11(1):119.
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, McWilliam H, Maslen J, Mitchell A, Nuka G: InterProScan 5: genome-scale protein function classification. Bioinformatics 2014, 30(9):1236-1240.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research 2000, 28(1):27-30.
Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn. Genome Research 2004.
Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, Fan TW-M, Fiehn O, Goodacre R, Griffin JL: Proposed minimum reporting standards for chemical analysis. Metabolomics 2007, 3(3):211-221.

Download PDF

Journal Publication

published 11 Sep, 2020

Read the published version in Microbiome →

Review #1 received at journal
16 May, 2020
Editorial decision: Minor revision
16 May, 2020
Review #2 received at journal
23 Apr, 2020
Reviewer #2 agreed at journal
20 Apr, 2020
Reviewer #1 agreed at journal
15 Apr, 2020
Reviewers invited by journal
14 Apr, 2020
Editor assigned by journal
19 Feb, 2020
Submission checks completed at journal
18 Feb, 2020
Editor invited by journal
18 Feb, 2020
First submitted to journal
14 Feb, 2020

You are reading this older preprint version

Read the latest preprint version →

Multi-omics analysis reveals the influence of genetic and environmental risk factors on developing gut microbiota in infants at risk of celiac disease

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Results

Discussion

Conclusions

Methods

List of abbreviations

Declarations

Ethics approval and consent to participate:

Clinical Trial Registration

Availability of data and material

Competing interests

Authors' contributions

Acknowledgements:

Funding

Supplementary Information

References

Supplementary Files

Status:

Journal Publication

Version 1