Human Milk Oligosaccharide-Stimulated Bifidobacterium Species Contribute to Prevent Later Respiratory Tract Infections

Editor(s):
Shaillay Kumar Dogra, Francois-Pierre Martin, Dominique Donnicola, Monique Julita †, Bernard Berger and Norbert Sprenger

Abstract

Human milk oligosaccharides (HMOs) may support immune protection, partly via their action on the early-life gut microbiota. Exploratory findings of a randomized placebo- controlled trial associated 2′fucosyllactose (2′FL) and lacto-N-neotetraose (LNnT) formula feeding with reduced risk for reported bronchitis and lower respiratory tract illnesses (LRTI), as well as changes in gut microbiota composition. We sought to identify putative gut microbial mechanisms linked with these clinical observations. (2) Methods: We used stool microbiota composition, metabolites including organic acids and gut health markers in several machine-learning-based classification tools related prospectively to experiencing reported bronchitis or LRTI, as compared to no reported respiratory illness. We performed preclinical epithelial barrier function modelling to add mechanistic insight to these clinical observations. (3) Results: Among the main features discriminant for infants who did not experience any reported bronchitis (n = 80/106) or LRTI (n = 70/103) were the 2-HMO formula containing 2′FL and LNnT, higher acetate, fucosylated glycans and Bifidobacterium, as well as lower succinate, butyrate, propionate and 5-aminovalerate, along with Carnobacteriaceae members and Escherichia. Acetate correlated with several Bifidobacterium species. By univariate analysis, infants experiencing no bronchitis or LRTI, compared with those who did, showed higher acetate (p < 0.007) and B. longum subsp. infantis (p ≤ 0.03). In vitro experiments demonstrate that 2′FL, LNnT and lacto-N-tetraose (LNT) stimulated B. longum subsp. infantis (ATCC15697) metabolic activity. Metabolites in spent culture media, primarily due to acetate, supported epithelial barrier protection. (4) Conclusions: An early-life gut ecology characterized by Bifidobacterium-species- driven metabolic changes partly explains the observed clinical outcomes of reduced risk for bronchitis and LRTI in infants fed a formula with HMOs. (Trial registry number NCT01715246.)

Keywords: Bifidobacterium; acetate; fucosylated glycans; respiratory infections; bronchitis; machine learning based classification models

 

1. Introduction

“Birth seeds, breast milk feeds” essentially captures the very first phase of early-life microbiota development and summarizes differences in microbiota development observed by mode of delivery and diet. The infant gut is progressively and sequentially colonized with distinct microbial populations during infancy and early childhood, from an aerobic to an anaerobic milk-oriented microbiome first, then to a more diverse adult-like microbiome, reaching a stable composition at around 2 to 3 years [1–3]. Typically, during the first few months of life, a high Bifidobacterium dominance is observed before the microbiome diversifies with the introduction of a complementary weaning diet [2,4,5]. Early and frequent antibiotic use, cesarean (C-) section birth and formula feeding alter the earlylife microbiome development pattern during the recommended exclusive breastfeeding period of the first 4–6 months [6,7]. This period of microbiome establishment is increasingly recognized as a critical window related to future health [8–11]. Further insights to support this have been gained from observations related to microbiome disruptors including antibiotic use, C-section birth and no or short breastfeeding that have been related to inappropriate development of immune competence, evidenced by higher risk of infections and allergies [12–15].

Human milk oligosaccharides (HMOs) are one of the major components of breast milk by mass, in order after lactose and fat. Major HMOs such as those composed of Nacetyllactosamine units attached to lactose, those containing fucose and many with sialic acid are basically missing in farmed animal milks [16]. Consequently, infant formula products manufactured from animal milks do not contain those major HMO classes, yet contain small amounts of naturally present galactosyllactoses and sialyllactoses [17,18]. Historically, HMOs were first described following the observation that breastfed and formula- fed infants are colonized with different microbes in early life [19]. In breastfed infants, specific HMOs such as 2′fucosyllactose (2′FL), lacto-N-tetraose (LNT) and lacto-Nneotetraose (LNnT) are supposed to be fermented by the gut microbiota, leading to increased Bifidobacterium abundance and metabolic activity [20–24]. Breastfed infants having Bifidobacterium breve strains able to use FL showed the highest Bifidobacterium abundance and a gut environment with the highest acetate concentration. Equally, LNT and LNnT were shown to specifically boost specific Bifidobacterium species such as B. longum subsp. infantis or B. breve in preclinical models [25–28].

We recently reported exploratory findings describing the effects of 2′FL and LNnT on the early-life gut microbiota of formula-fed infants, and demonstrated that formula supplementation with these two HMOs shifted the microbiota closer to that observed in breastfed infants and was associated with reduced risk of antibiotic use later in life [29,30]. Notably, 2-HMO formula-fed infants had a significantly reduced risk of experiencing reported bronchitis and lower respiratory tract illnesses (LRTI), yet gut microbial mechanisms and markers in relation to these clinical results remain unknown [29]. To this end, we employed machine-learning-based methods [31] to identify microbiota taxa, metabolites and stool gut health indicator proteins measured at 3 months of age that associated with susceptibility to bronchitis or LRTI from 3 to 12 months of infant age. Using this approach, we generated a hypothesis on the possible mechanisms of HMO action, which were further tested.

 

2. Materials and Methods

2.1. Clinical Trial

The trial was registered in www.clinicaltrial.gov with number NCT01715246 and was previously described for the full trial population [29] and for the per protocol formula-fed trial population [30] that we further explored here.

The trial population for the formula-fed groups consisted of healthy, full-term male and female infants from birth to 14 days old at enrollment, who were exclusively formula fed at the time of enrollment. Eligible infants were randomly assigned to one of two study formulas (Control or Test) using mode of delivery (vaginal or cesarean section) and gender as stratification factors (randomization was carried out using a permuted block algorithm with Medidata Balance (New York, NY, USA)). Randomized infants received exclusive feedings with the Test or Control formulas from enrollment through 4 months of age in amounts suitable for their weight, age and appetite. The Test formula was an intact protein cow’s milk protein-based infant formula with 2 HMOs (2′fucosyllactose (2′FL) and lacto-N-neotetraose (LNnT)) at 1.5 g/L and a ratio of 2:1 between the 2 HMOs. The Control formula was the same intact protein cow’s milk protein-based infant formula without the 2 HMOs, but with 1.5 g/L additional lactose.

Randomized infants were allowed the introduction of complementary weaning food and continued feedings with the Control or Test formulas through 6 months of age, and then both groups received the same intact protein cow’s-milk-based follow-up formula for feedings from 6 to 12 months of age.

The primary objective of the randomized, controlled clinical trial was to evaluate the non-inferiority of weight gain from enrollment to 4 months of age in infants fed the Test formula, compared with those fed the Control formula. Secondary objectives included the evaluation of differences between the formula groups in body weight, body length, head circumference, digestive tolerance, formula compliance, stool protein markers for intestinal status, stool microbiota, medication use and morbidity through 12 months of age. Infection- related morbidity was assessed based on an a priori-formulated hypothesis that the 2 HMOs help to prevent infections [30].

Adverse events were reported as described in Puccio et al. [29], then coded and categorized by a single physician, who was not involved in the study conduct using the Medical Dictionary for Regulatory Activities (MedDRA) System, Organ and Class (SOC) categories as well as the Preferred Terms (PTs) within each SOC category. Here, we focused on the per protocol population as described in Berger et al. [30] and on adverse events observed to be reduced with 2-HMO formula feeding [29], namely the reported PT bronchitis (up to 12 months) and the LRTI adverse event cluster (including the PTs: bronchiolitis, bronchitis, pneumonia, LRTI, LRTI viral, RSV (respiratory syncytial virus) bronchiolitis, RSV bronchitis, viral respiratory infection). For PT bronchitis, different values exist based on the incidence of occurrence (0 = 91; 1 = 22; 2 = 4; 3 = 4). Here, we considered any incidence through 12 months of age as a case (Y) and no incidence as a control (N). For PT bronchitis, we had 91 controls and 30 cases. Since we aimed to investigate the 3-month stool markers in a prognostic manner, we removed those cases who only had an incidence of PT bronchitis up to three months, but not from 3 to 12 months. This resulted in a total of 26 cases. Further, three samples from controls were removed in the modeling due to missing values and eight subjects from controls were removed as they did not finish the study, leaving us with 80 controls.

Similarly, for LRTI (up to 12 months) we had some infants with several incidences (0 = 81; 1 = 31; 2 = 5; 3 = 4). Again, we counted any number of incidences as a case (Y) and no incidence as a control (N). Thus, we had for LRTI (up to 12 months) 81 controls and 40 cases. Next, we removed cases with an incidence of LRTI only up to three months, but not from 3 to 12 months, leaving us with 33 cases. Further, three samples from controls were removed in the modeling because of missing values and 4 subjects from controls were removed as they did not finish the study, leaving us with 74 controls.

The adverse event cluster LRTI also comprises PT bronchitis. Hence, we had all 26 bronchitis cases also represented in the LRTI cases, where we had an additional seven cases. These 7 LRTI cases were part of the control PT bronchitis strata. Of the 111 infants in total, we had 78 common controls and 26 common cases for PT bronchitis and LRTI from 3 to 12 months of age.

2.2. Stool Microbiota Composition, Metabolites and Gut Health Markers

Stool samples were processed for microbiota composition analysis by V3 and V4 region of the 16S rRNA gene amplicon sequencing, as previously described [30]. Notably, the employed bioinformatics analysis pipeline allowed for species and subspecies identification of the genus Bifidobacterium [32].

The biochemical composition of the stool samples was analyzed using a well-established metabonomic approach via proton Nuclear Magnetic Resonance (1H NMR) Spectroscopy [33,34]. Quantitative profiles of major metabolites were acquired, including amino acids (phenylalanine, tyrosine, isoleucine), short-chain fatty acids (SCFA) (propionate, butyrate, acetate, valerate, iso-valerate, 5-amino-valerate), organic acids (lactate, succinate) and carbohydrates (fucosyl-glycans).

The gut health markers calprotectin, α1-antitrypsin and elastase were quantified by Biotecon Diagnostics GmbH (Potsdam, Germany). Quantification of α1-antitrypsin was carried out using a stool extraction kit and enzyme-linked immunosorbent assay (ELISA) provided by Immundiagnostic AG (Bensheim, Germany). Elastase was quantified using the Schebo® Quick-PrepTM E1-stool extraction system (Schebo Biotech, Gießen, Germany) and the ScheBo® Pancreatic Elastase 1 assay (ScheBo® Biotech). Calprotectin was quantified upon extraction with the Bühlmann Smart-Prep (Dresden, Germany) stool extraction system with the EK-CAL ELISA kit (Bühlmann, Dresden, Germany).

2.3. In Vitro Fermentation Assay and Analysis

We used as carbohydrate sources the HMOs LNT (Lacto-N-tetraose), LNnT, 2′FL, 6′SL (6′sialyllactose) (all from Glycom A/S, Lyngby, Denmark) together with glucose and lactose (Fluka, Buchs, Switzerland), to ferment using the B. longum subsp. infantis type strain (ATCC 15697). The B. longum subsp. infantis was grown overnight in glucose-free MRS (DeMan, Rogosa and Sharpe) broth supplemented with 1% of each carbohydrate source in anaerobic conditions at 37 °C. Anaerobic conditions were achieved using BD GasPak system (Becton, Dickinson and Company, Franklin Lakes, NJ, USA). The optical density (OD) at 600 nm was measured before pelleting the bacteria by centrifugation and resuspension in PBS. To reach a final OD of 0.1, an aliquot of each preparation was added to 10 mL of DMEM cell culture media containing 0.1% glucose (Gibco™: DMEM, low glucose, GlutaMAX™ Supplement, pyruvate) further supplemented with 1% of each carbohydrate source for one night at 37 °C in anaerobic conditions. The OD and pH were measured and after centrifugation the conditioned media supernatants were sterile filtered using a 0.22 μm filter. These spent culture media were then used for further experiments and analysis.

Organic acids and short-chain fatty acids were quantified in the conditioned DMEM media using an Ultimate 3000 HPLC system (Thermo Fisher Scientific, Dionex, Sunnyvale, CA, USA) equipped with a HiPlex H analytical column (8 μm; 300 × 7.7 mm; Polymer Laboratories (Varian, London, UK)) and an U3000 RS diode array detector (Thermo Fisher Scientific) set to an absorption wavelength of 210 nm. The system was operated isocratically with 5 mM sulfuric acid at a flow rate of 0.6 mL/min and a temperature of 55 °C. Quantification was carried out using external standard curves from 0.5 to 100 mM of the respective authentic organic acids (Fluka (Sigma), Switzerland). Dionex Chromeleon software version 6.8 (Thermo Fisher Scientific) was employed to pilot the HPLC and to quantify the organic acids.

For NFkB stimulation analysis, HT−29 clone 34 cells were seeded in a 96-well plate at 2 × 104 cells/well and used 2 days after 80% confluence was achieved. Culture media were replaced with fresh antibiotic-free media for 1 day and then replaced with the media containing the tested carbohydrates and conditioned media. Cells were then either challenged with 2 ng/well TNFα (Bio-Techne AG, Zug, Switzerland) or Salmonella SL1344 (NCTC 13347) at 106 cfu/well. After 24 h at 37°C at 5% CO2, the NFkB reporter luminescence was measured using the Phospha-light system (Applied Biosystems, Foster City, CA, USA) according to manufacturer’s instructions and expressed as the relative amount compared to the unconditioned glucose control.

Caco-2 (ATCC HTB-37) cells were seeded in 24-well plates at 1 × 105 cells/well and cultured for 14 days at 5% CO2 and 37 °C before performing the adhesion/invasion assay with Salmonella SL 1344 (NCTC 13347). Salmonella SL1344 was grown to logarithmic phase in LB (Luria Bertani) broth and diluted to an OD 0.1 (approximately 1 × 108 cfu/mL) in DMEM at pH 5.8. The day before the experiment, the Caco-2 cell culture media were replaced with DMEM without antibiotics. DMEM media containing the different carbohydrates and the conditioned media were adjusted to pH 5.8 and added to Caco-2 cells together with 1 × 107 cfu/well Salmonella SL1344. Cells were incubated for 60 min. Thereafter, cells were washed 3 times with fresh culture media and incubated for another 60 min with culture media containing 100 μg/mL gentamycin. Cells were again washed 3 times with PBS (phosphate-buffered saline, Sigma, St. Louis, MO, USA) to remove gentamycin, before Caco-2 cells were lysed with sterile MilliQ-water with 1% Triton X-100 (Sigma) and plated on LB agar plates to enumerate the Salmonella that invaded the cells.

2.4. Machine Learning Algorithms and Statistical Analysis

We used sparse partial least square discriminant analysis (sPLS-DA) [35] to investigate the separation of cases and controls based on three-months-of-age stool microbiota composition, metabolites and gut health markers, as well as mode of delivery and formula type (with or without HMOs). The modeling was carried out in R [36] using the package caret [37] and mixOmics [38] using standard cross-validation methods (10-fold with 10 repeats) to avoid over-fitting of the models. The sPLS-DA modeling was run on balanced error rate (BER), which calculates the average fraction of mis-classified samples in each class, weighted by the number of samples in each class. Additionally, we used random down-sampling to balance out the two classes—cases and controls—as we had fewer cases than controls. For statistical analysis, non-parametric tests were performed (Mann–Whitney and Kruskal–Wallis tests) using GraphPad Prism software version 8.3.0.

 

3. Results

3.1. Cohort Description

We sought to further explore microbiota taxa, metabolites and stool gut health indicator proteins from per protocol infants with stool samples at three months of age, and susceptibility to reported bronchitis or LRTI from 3 through 12 months of age [29,30]. Demographics and infant characteristics are shown in Table 1. The flow chart (Figure 1) depicts the number of infants randomized to the trial and the subgroup of infants available for the exploratory study reported here. The grouping of the formula-fed infants by controls and cases is equilibrated for all parameters except for their exposure to the 2-HMOcontaining formula.

 

Fig 1. Flowchart.

Figure 1. Subject flow chart showing the number of per protocol formula-fed infants with stool sample at three months of age, by controls and cases who experienced at least one reported bronchitis or lower respiratory tract infection (LRTI) between three and 12 months of age.

 

Table 1. Subject characteristics

 

# Values are mean ± standard deviation unless otherwise noted; * parameters by Test and Control feeding group were previously reported [29,30]. For categorical data, difference by two-tailed Chi-square test with Fisher exact probability test; for numerical data, by t-test. ns, not significant; Cases, from 3 to 12 months of age, any bronchitis or LRTI; Controls, no bronchitis or LRTI through 12 months of age; LRTI, lower respiratory tract infection adverse events cluster.

3.2. Several Microbiome and Host Features Relate to Lower Respiratory Tract Infections

We used 1H-NMR-identified metabolites, 16S gene-based microbiota composition (at genus and for bifidobacteria at species and subspecies level) and host gut health indicator proteins from infant stool at three months of age along with Test formula and mode of delivery to generate a total of 60 features used to categorize infants who experienced (Yes class) or did not experience (No class) bronchitis or LRTI from 3 to 12 months of age.

The class imbalance of data with controls being about three times greater as compared to cases (bronchitis: 80 controls and 26 cases; LRTI: 70 controls and 33 cases) posed a challenge, particularly for the Random Forest and Support Vector Machine algorithms. These methods required up-sampling to obtain acceptable class separation. Yet, using Random Forest with Boruta feature selection identified relative acetate, butyrate and 5- aminovalerate, together with Lactobacillus and succinate, as the most discriminant features (data not shown). Support vector machine combined with recursive feature elimination identified the same features together with fucosyl-glycans, lactate and Bifidobacterium longum subsp. infantis and other Bifidobacterium species, to name the most prominent (data not shown).

Classification using sPLS-DA with random down-sampling of controls was able to reasonably separate the controls from cases for bronchitis and LRTI, with a median cross validation overall accuracy, from 10-fold with 10 repeats, of about 60 to 70% and 60 to 65% across the 10 repeats for bronchitis and LRTI, respectively, with component 1 alone (Figure S1). We next focused on a more detailed analysis using five independent rounds of random down-sampling of controls for both bronchitis and LRTI. Across the rounds, the cases are always the same and it is controls that are randomly chosen, given that cases are far fewer than controls. The sensitivity–specificity receiver operating characteristic (ROC) curve with component 1 showed an area under the curve (AUC) between 0.74 (p-value 0.004) and 0.82 (p = 6.6 × 10−5) for bronchitis and between 0.71 (p = 0.003) and 0.8 (p = 3.7 × 10−5) for LRTI, across the five runs, respectively (Table S1). Class separation from representative examples of the five independent random down-sampling rounds are shown in Figure 2, together with the features of component 1 used by sPLS-DA to separate the classes for bronchitis and LRTI.

Consistently, separation of controls from bronchitis cases is driven by 2-HMO supplementation followed by fucosylated glycans, relative amount of acetate, lactate and the genus Bifidobacterium, with higher relative amounts in controls. Features picked by the model indicating cases in component 1 are the organic acid succinate, the relative amount of the SCFAs butyrate and 5-amino-valerate and propionate, as well as calprotectin, Csection birth, Trabulsiella, Escherichia and the family of unidentified members of the Carnobacteriaceae.

Fig 2.

Figure 2. Representative sPLS-DA plots and loadings of the five independent rounds of sPLS-DA modeling with downsampling to separate controls and cases for bronchitis (n = 52) (a,b) and LRTI (n = 66) (c,d). Results from all five sPLS-DA modeling rounds for both bronchitis and LRTI are provided in Table S1.

Next, we performed univariate analysis of the sPLS-DA-identified features using the full data set to compare feeding groups. We observed statistically significantly higher relative fucosyl-glycan amounts and Bifidobacterium abundance in feces of 2-HMO (Test) formula- fed infants, and lower abundance of the unidentified Carnobacteriaceae and Escherichia (Figure 3).

Fig 3.

Figure 3. Comparison of the identified main features from the sPLS-DA model by feeding group: (a,b) relative abundance of metabolites; (c) relative abundance of microbial taxa; (d) concentration of calprotectin. (n = 106; * indicate p < 0.05, ** p < 0.01 by non-parametric t-test).

Together, the case–control analysis indicates that an early-life gut ecology characterized by higher relative abundance of Bifidobacterium, together with higher relative acetate and lower relative butyrate, propionate and 5-aminovalerate, associates with protection from reported bronchitis and LRTI from three to 12 months of age. Therefore, we further investigated the relationship between relative acetate and the main Bifidobacterium species identified in our study. We observed a positive correlation between relative acetate amount in stool and the sum of all Bifidobacterium and the sum of the B. longum taxa together (Figure 4a). B. longum subsp. longum and B. longum subsp. infantis. B. breve also showed a positive although less strong correlation (Figure 4a). No statistically significant correlation was seen with B. bifidum, B. dentium and B. adolescentis. Similarly, the linear regression analysis illustrates a statistically significant correlation between the relative Bifidobacterium and acetate abundances (Figure 4b).

Fig 4.

 

Figure 4. Association between relative acetate and Bifidobacterium species in infant stool samples. (a) Spearman correlation; (b) linear regression analysis of the sum of Bifidobacterium species and relative amounts of acetate.

To complete the analysis with the full data set, we investigated in a univariate approach the relative content in SCFAs and Bifidobacterium taxa abundances, in controls and cases for both bronchitis and LRTI from 3 to 12 months of age (Figure 5). As expected from the machine-learning-based classification models, the relative acetate amount was higher in the stool samples of control infants at three months of age as compared with cases who experienced at least one bronchitis or LRTI event thereafter until 1 year of age. For butyrate and 5-aminovalerate, we observed higher relative amounts in the three-month stool samples of infants who experienced at least one bronchitis or LRTI event thereafter. Equally, for propionate, higher relative amounts were observed that reached statistical significance only for LRTI. Interestingly, only B. longum subsp. infantis was consistently higher in the stool of control infants as compared to infants who experienced bronchitis and LRTI. Yet, the B. longum subsp. infantis was only a minor contributor to the total Bifidobacterium population in our study infants. For LRTI, the sum of B. longum showed a slightly higher relative abundance in controls compared to cases.

Fig 5.

Figure 5. Comparison of relative amounts of short-chain fatty acids and Bifidobacterium species in the stool of infants at three months of age between infants who did not (controls) or did (cases) experience any bronchitis and LRTI thereafter from 3 to 12 months of age: (a,b) relative abundance of the short-chain fatty acids for bronchitis or LRTI cases and controls; (c,d) relative abundance of the Bifidobacterium species for bronchitis or LRTI cases and controls (* indicate p < 0.05, ** p < 0.01 by non-parametric t-test).

The identified features from the clinical trial indicate an important role of bifidobacteria metabolism, especially B. longum subsp. infantis and other Bifidobacterium species able to ferment specific HMOs. Therefore, we further explored whether and how HMO-stimulated Bifidobacterium metabolites alter epithelial barrier function upon an inflammatory or pathogen challenge in an experimental model.

3.3. Acetate Produced by HMO Stimulated Bifidobacterium longum subsp. infantis Supports Epithelial Barrier Function

B. longum subsp. infantis (ATCC 15697) were grown in the presence of additional glucose, lactose or the HMOs 2′FL, LNnT, LNT or 6′SL. LNT was assessed due to its structural similarity to LNnT and 6′SL, as it is structurally very different from the other HMOs utilized. B. longum subsp. infantis grew on all substrates, with lactose and LNT generating the highest bacterial growth (Figure 6a). In terms of metabolic activity, we found the highest amounts of acetate in the spent culture media from B. longum subsp. infantis conditioned with LNT and LNnT as substrates (Figure 6b). No changes in pH were observed (Figure 6c). Moreover, spent culture media from B. longum subsp. infantis conditioned with LNT or LNnT strongly reduced TNFα- and Salmonella SL1344 (data not shown)- mediated NFkB activation (Figure 6d). Invasion of Salmonella in Caco-2 cells was significantly reduced by spent culture media from B. longum subsp. infantis conditioned with 2′FL, LNT and LNnT (Figure 6e). Since acetate was particularly increased in spent culture media from B. longum subsp. infantis conditioned with LNT and LNnT, we tested whether acetate alone was sufficient to mediate the observed epithelial barrier protection. Acetate dose dependently reduced Salmonella SL1344 invasion and importantly, the addition of acetate to spent culture media from B. longum subsp. infantis conditioned with glucose led to equivalent protection from Salmonella SL1344 as compared to spent culture media from B. longum subsp. infantis conditioned with LNT (Figure 6f). Hence, acetate is a key Bifidobacterium-produced metabolite stimulated by specific HMOs, such as LNT and LNnT, that promote epithelial barrier protection.

Fig 6.

Figure 6. In vitro experiments with B. longum subsp. infantis conditioned with glucose (Glc), lactose (Lac), 2′fucosyllactose (2′FL), lacto-N-tetraose (LNT), lacto-N-neotetraose (LNnT) and 6′sialyllactose (6′SL): (a) bacterial growth monitoring using optical density at 600 nm (OD600); (b) quantification of acetate in spent culture media from conditioned B. longum subsp. infantis; (c) measurement of conditioned spent culture medium pH; (d) measurement of NFkB activation by TNFα in the presence of control and conditioned (cond) spent culture media; (e) measurement of Salmonella invasion in cultured epithelial cells in the presence of control and conditioned spent culture media; (f) measurement of Salmonella invasion with different amounts of acetate added to fresh or spent culture media in comparison to control and spent culture media. Symbols stand for independent experiments each with 3 replicates (* p < 0.05, ** p < 0.01 by Kruskal–Wallis or non-parametric t-test).

 

4. Discussion

The developing gut microbiome during early life plays a key role in long-term health through its action on mucosal barrier function and immune development [10,22,39]. Yet, the age-specific role of key microbial taxa dominating the early-life gut ecosystem as well as key environmental and nutritional factors influencing their activity remains unclear. Here, we used a clinical intervention trial with 2-HMO and control formula-fed infants to identify early-life infant parameters and microbiome features related to subsequent reported LRTI, including bronchitis. We identified that feeding with a 2-HMO formula, associated with reduced risk for reported bronchitis and LRTI, is characterized by a distinct gut microbial and metabolite phenotype. This phenotype is composed of higher stool acetate, fucosyl-glycans and Bifidobacterium species abundance, along with lower content in butyrate, propionate, 5-aminovalerate and succinate, and lower Escherichia species abundance. Among these features, relative acetate and B. longum subsp. infantis abundance were higher in the stool of infants with no reported bronchitis or LRTI, suggesting an association with protection against respiratory tract infection. Our in vitro experiments further corroborate these findings, through the demonstration of metabolic activity from HMO-conditioned Bifidobacterium longum subsp. infantis, and especially the acetate production that enhanced epithelial barrier function in cell culture models.

The short-chain fatty acid acetate is a key metabolite of bifidobacteria [40], and in our study we observed a significant correlation between these two features. Additionally, by univariate analysis we observed significant differences in relative acetate and B. longum subsp. infantis measured in stool at three months of age when comparing infants who experienced bronchitis or LRTI thereafter or did not. From the correlation matrix between bifidobacteria and acetate, we can assume that other Bifidobacterium species such as B. longum subsp. longum and B. breve also contributed to the higher in vivo production of acetate. All those species are known to have strains that can use specific HMOs as substrate [19,40–42]. Notably, HMO utilization capacity by bifidobacteria was previously shown to lead to higher Bifidobacterium species abundance, with higher acetate concentration and lower stool pH [19].

In our study, B. longum subsp. infantis was associated with protection against respiratory infection in early life. Therefore, we decided to use a B. longum subsp. infantis well adapted to use HMOs [43,44] to further explore causality using a series of in vitro cell culture models. First, we confirmed previous observations, indicating that specific HMOs, mainly LNT, LNnT and 2′FL, but not 6′SL, stimulated B. longum subsp. infantis growth and metabolic activity. Second, we showed that the generated metabolic phenotype of stimulated B. longum subsp. infantis reduced NFkB activation and protected epithelial cells from invasion by Salmonella. Third, we demonstrate that acetate is sufficient to provide the epithelial cell protection. Hence, HMO-stimulated B. longum subsp. infantis type strain exerts epithelial barrier protection. In a series of elegant preclinical animal studies, acetate from bifidobacteria was shown to mediate protection from enterohaemorrhagic E. coli [45]. The authors speculate that acetate improved the epithelial barrier function, thereby providing enhanced protection from the pathogenic E. coli. Recently, acetate was also shown to drive gut secretory IgA production and specificity [46], which provides a mechanistic angle. The epithelial barrier protective effect seen for acetate is not confined to the intestine and provides protection via the gut–lung axis [47]. In adults, colonic-administration 13C-SCFA showed the highest systemic availability for acetate (36%) compared to propionate (9%) and butyrate (2%), indicating that acetate can mediate beneficial effects at other body sites [48]. In a mouse model, acetate was shown to protect from RSV, a major cause of LRTI in infants [49]. Acetate-mediated activation of lung epithelial Gpr43 reduced the cytotoxic effects of RSV and stimulated an interferon-mediated anti-viral response [49]. In another animal model, acetate was further shown to contribute to protection from bacterial superinfection of the lungs following influenza virus infection through the activation of alveolar macrophages [50]. Although bifidobacteria-produced acetate in early life may not be solely responsible, our findings suggest a likely contribution to infant protection from LRTI later in life.

The early-life gut ecology with higher relative acetate and specific bifidobacteria in infants who did not experience any bronchitis or LRTI through 12 months contrasts that of infants who did experience any bronchitis or LRTI, which is characterized by higher relative butyrate, propionate and 5-aminovalerate and succinate. This observation indicates that an age-appropriate gut microbiome development may be important for protection against respiratory infection. In a recent study, Tsukuda et al. reported SCFA and other organic acid profiles in stool samples during the first 2 years of age from a small cohort of mainly breastfed and vaginally born infants [51]. An early succinate peak was paralleled by increasing acetate, while propionate and butyrate increased with cessation of breastfeeding. This is consistent with the general observation of the development of gut microbiota diversity and composition, dominated by bifidobacteria in early life and diversifying rapidly with complementary diet [6,52]. A delay in the bifidobacteria dominance is generally observed with cesarean-section-born infants [5,53], and a recent study linked such a delay together with further alterations in the gut microbiota, leading to higher risk for respiratory infections [12]. Thus, an acetate- and bifidobacteria-dominated gut ecology during the first few months of life seems important for protection against respiratory infection, and specific HMOs support this protection.

While important biological functions leading to enhanced protection can be attributed to acetate, a bifidobacteria-dominated gut ecology associated with other biochemical molecules may also be of importance for host physiology and immune development. An important group of molecules are those derived from the aromatic amino acids tryptophan, phenylalanine and tyrosine by HMO-utilizing Bifidobacterium metabolism; for example, the aromatic lactic acids indolelactate, phenyllactate and 4-hydroxyphenyllactate [22,54]. Indolelactate has been shown to activate the aryl-hydrocarbon receptor. This is of particular interest, as aryl-hydrocarbon receptor signaling is important for intestinal homeostasis, barrier function in gut, skin and lungs and immunity [55,56]. Specifically, B. longum subsp. infantis supplementation to breastfed infants silenced Th2 and Th17 cytokines and induced interferon beta by an indolelactate-mediated upregulation of galectin- 1 [22]. Hence, several metabolites from key HMO-utilizing bifidobacteria act through different molecular routes to strengthen protection against infection.

Using machine learning approaches, we robustly identified other features together with acetate that may contribute to enhanced protection. By our 1H NMR metabolomics approach, the signals integrated as fucosyl-glycans do not allow any exact structural annotation due to the complex biological mixture represented in feces. Yet, through spiking experiments we have confirmed that the 1H NMR signal does not correspond to 2′FL, which was part of the 2-HMO formula ingredients. However, among the identified “fucosyl- glycan” are signals of free fucose and fucose bound to other unidentified glycans. It is possible that part of the fucosyl-glycan signals stem from mucous glycans and may indicate gut mucosal fucosylation changes. Free fucose liberated from 2′FL upon microbial action may have further strengthened mucosal barrier function, possibly through upregulation of the fucosyltransferase 2 activity, as recently elaborated in preclinical models [57].

The exploratory approach we took identified several features indicating different possible early-life-gut-ecology-mediated mechanisms related to later protection from respiratory infection. Yet, study limitations must be emphasized here. Besides the limitations of the trial [1,2], primarily related to the reporting of respiratory illnesses as part of the adverse event reporting, our case–control classification approach suffered from a relatively low number of cases versus controls and from a relatively low number of subjects and features. With such limitations, some machine learning classification algorithms do not perform at their best. Therefore, we initially explored the data using multiple classification methods (Random Forests with Boruta feature selection, Support Vector Machines [58] and sPLS-DA), all based on different algorithms, and observed some features that were consistently selected by the different methods. Among them are the relative acetate amounts, fucosyl-glycans, product and relative butyrate and succinate. For the detailed analysis presented here, we chose to use sPLS-DA. To overcome the class imbalance, we ran the sPLS-DA classification in five independent runs by randomly selecting samples from controls to balance with the cases and found moderate separation, with the primary purpose of hypothesis generation.

 

5. Conclusions

Machine-learning-based classification proved useful to identify possible mechanisms relating early-life gut ecology to subsequent respiratory infection using data from a wellstratified and controlled intervention trial with formula-fed infants. Among the main identified features discriminating between cases and controls for early-life respiratory infections were the 2-HMO formula and the relative stool acetate, butyrate and 5-aminovalerate, as well as succinate and fucosyl-glycans. Together, these indicate the importance of a distinguished early-life gut ecology with higher relative acetate, mediated by HMOutilizing bifidobacteria for protection against respiratory infection. This hypothesis merits further investigation, including additional detailed metabolite analysis and functional assays to assess their contribution to appropriate immune maturation during this critical early-life window.

Supplementary Materials: The following are available online at www.mdpi.com/article/ 10.3390/microorganisms9091939/s1, Figure S1: Representative sPLS-DA model accuracy to separate cases from controls, Table S1: Loadings of 5 independent rounds of sPLS-DA modeling to separate controls and cases for bronchitis (n = 52) and LRTI.

Author Contributions: Conceptualization, N.S., S.K.D., F.-P.M. and B.B.; methodology, D.D., M.J., S.K.D. and F.-P.M.; formal analysis, N.S., S.K.D., F.-P.M., D.D. and B.B.; writing—original draft preparation, N.S.; writing—review and editing, N.S., S.K.D., F.-P.M. and B.B. All authors have read and agreed to the published version of the manuscript.

Funding: This research received no external funding.

Institutional Review Board Statement: The study was conducted according to the guidelines of the Declaration of Helsinki and the International Conference on Harmonization guidelines for Good Clinical Practice between October 2012 and July 2015 in the Dipartimento Materno Infantile AOUP ‘‘Paolo Giaccone,’’ Universita` di Palermo, Palermo, Italy, and in the Department of Paediatrics at Jessa Hospital in Hasselt, Belgium. The study was approved by the Institutional Review Boards of both hospitals (protocol code CT11.26.INF and 22. March 2012 and 16. May 2012). The study was registered at www.clinicaltrials.gov under the number: NCT01715246.

Informed Consent Statement: Informed consent was obtained from the parent or legal guardian of each infant before enrollment.

Data Availability Statement: The microbiota data presented in this study were previously reported [30] and are openly available in SRA under accession number SRP151522.

Acknowledgments: The authors kindly thank all the infants and their families who took part in this study, as well as the clinical investigators P. Alliet, G. Puccio and their teams. For critical reading of the manuscript and support we thank Elizabeth Forbes Blom and Eline van der Beek. This study was funded by Société des Produits Nestlé SA, Switzerland.

Conflicts of Interest: All authors are employees or were previous employees of Société des Produits Nestlé, S.A. Switzerland. The authors declare no other conflict of interest. The employer had no role in the design of the study; the collection, analyses or interpretation of data; the writing of the manuscript, nor in the decision to publish the results.

Abbreviations

HMO, human milk oligosaccharide; 2′FL, 2′fucosyllactose; LNnT, lacto-N-neotetraose; LNT, lacto- N-tetraose; 6′SL, 6′sialyllactose; C-section, cesarean section; LRTI, lower respiratory tract infections; SCFA, short-chain fatty acids; ELISA, enzyme-linked immunosorbent assay; sPLS-DA, sparse partial least square discriminant analysis.

 

References

1. Dogra, S.; Sakwinska, O.; Soh, S.E.; Ngom-Bru, C.; Bruck, W.M.; Berger, B.; Brussow, H.; Lee, Y.S.; Yap, F.; Chong, Y.S.; et al. Dynamics of infant gut microbiota are influenced by delivery mode and gestational duration and are associated with subsequent adiposity. mBio 2015, 6, e02419-14, doi:10.1128/mBio.02419-14.
2. Stewart, C.J.; Ajami, N.J.; O'Brien, J.L.; Hutchinson, D.S.; Smith, D.P.; Wong, M.C.; Ross, M.C.; Lloyd, R.E.; Doddapaneni, H.; Metcalf, G.A.; et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nature 2018, 562, 583–588, doi:10.1038/s41586-018-0617-x.
3. Roswall, J.; Olsson, L.M.; Kovatcheva-Datchary, P.; Nilsson, S.; Tremaroli, V.; Simon, M.C.; Kiilerich, P.; Akrami, R.; Kramer, M.; Uhlen, M.; et al. Developmental trajectory of the healthy human gut microbiota during the first 5 years of life. Cell Host Microbe 2021, 29, 765–776 e763, doi:10.1016/j.chom.2021.02.021.
4. Laursen, M.F.; Bahl, M.I.; Michaelsen, K.F.; Licht, T.R. First Foods and Gut Microbes. Front. Microbiol. 2017, 8, 356, doi:10.3389/fmicb.2017.00356.
5. Korpela, K.; de Vos, W.M. Early life colonization of the human gut: Microbes matter everywhere. Curr. Opin. Microbiol. 2018, 44, 70–78, doi:10.1016/j.mib.2018.06.003.
6. Baumann-Dudenhoeffer, A.M.; D'Souza, A.W.; Tarr, P.I.; Warner, B.B.; Dantas, G. Infant diet and maternal gestational weight gain predict early metabolic maturation of gut microbiomes. Nat. Med. 2018, 24, 1822–1829, doi:10.1038/s41591-018-0216-2.
7. Robertson, R.C.; Manges, A.R.; Finlay, B.B.; Prendergast, A.J. The Human Microbiome and Child Growth—First 1000 Days and Beyond. Trends Microbiol. 2019, 27, 131–147, doi:10.1016/j.tim.2018.09.008.
8. Stinson, L.F. Establishment of the early-life microbiome: A DOHaD perspective. J. Dev. Orig. Heal. Dis. 2020, 11, 201–210, doi:10.1017/S2040174419000588.
9. Ronan, V.; Yeasin, R.; Claud, E.C. Childhood Development and the Microbiome-The Intestinal Microbiota in Maintenance of Health and Development of Disease During Childhood Development. Gastroenterology 2021, 160, 495–506, doi:10.1053/j.gastro.2020.08.065.
10. de Steenhuijsen Piters, W.A.A.; Binkowska, J.; Bogaert, D. Early Life Microbiota and Respiratory Tract Infections. Cell Host Microbe 2020, 28, 223–232, doi:10.1016/j.chom.2020.07.004.
11. Victora, C.G.; Bahl, R.; Barros, A.J.; Franca, G.V.; Horton, S.; Krasevec, J.; Murch, S.; Sankar, M.J.; Walker, N.; Rollins, N.C. Breastfeeding in the 21st century: Epidemiology, mechanisms, and lifelong effect. Lancet 2016, 387, 475–490, doi:10.1016/S0140- 6736(15)01024-7.
12. Reyman, M.; van Houten, M.A.; van Baarle, D.; Bosch, A.; Man, W.H.; Chu, M.; Arp, K.; Watson, R.L.; Sanders, E.A.M.; Fuentes, S.; et al. Impact of delivery mode associated gut microbiota dynamics on health in the first year of life. Nat. Commun. 2019, 10, 4997, doi:10.1038/s41467-019-13014-7.
13. Arrieta, M.C.; Stiemsma, L.T.; Dimitriu, P.A.; Thorson, L.; Russell, S.; Yurist-Doutsch, S.; Kuzeljevic, B.; Gold, M.J.; Britton, H.M.; Lefebvre, D.L.; et al. Early infancy microbial and metabolic alterations affect risk of childhood asthma. Sci. Transl. Med. 2015, 7,
307ra152, doi:10.1126/scitranslmed.aab2271.
14. Feehley, T.; Plunkett, C.H.; Bao, R.; Choi Hong, S.M.; Culleen, E.; Belda-Ferre, P.; Campbell, E.; Aitoro, R.; Nocerino, R.; Paparo, L.; et al. Healthy infants harbor intestinal bacteria that protect against food allergy. Nat. Med. 2019, 25, 448–453, doi:10.1038/s41591-018-0324-z.
15. Samuel, T.M.; Binia, A.; de Castro, C.A.; Thakkar, S.K.; Billeaud, C.; Agosti, M.; Al-Jashi, I.; Costeira, M.J.; Marchini, G.; Martinez-Costa, C.; et al. Impact of maternal characteristics on human milk oligosaccharide composition over the first 4 months of lactation in a cohort of healthy European mothers. Sci. Rep. 2019, 9, 11767, doi:10.1038/s41598-019-48337-4.
16. Barile, D.; Tao, N.; Lebrilla, C.B.; Coisson, J.D.; Arlorio, M.; German, J.B. Permeate from cheese whey ultrafiltration is a source of milk oligosaccharides. Int. Dairy J. 2009, 19, 524–530, doi:10.1016/j.idairyj.2009.03.008.
17. Urashima, T.; Taufik, E.; Fukuda, K.; Asakuma, S. Recent advances in studies on milk oligosaccharides of cows and other domestic farm animals. Biosci. Biotechnol. Biochem. 2013, 77, 455–466, doi:10.1271/bbb.120810.
18. Kunz, C. Historical aspects of human milk oligosaccharides. Adv. Nutr. 2012, 3, 430S–439S, doi:10.3945/an.111.001776.
19. Matsuki, T.; Yahagi, K.; Mori, H.; Matsumoto, H.; Hara, T.; Tajima, S.; Ogawa, E.; Kodama, H.; Yamamoto, K.; Yamada, T.; et al. A key genetic factor for fucosyllactose utilization affects infant gut microbiota development. Nat. Commun. 2016, 7, 11939, doi:10.1038/ncomms11939.
20. Korpela, K.; Salonen, A.; Hickman, B.; Kunz, C.; Sprenger, N.; Kukkonen, K.; Savilahti, E.; Kuitunen, M.; de Vos, W.M. Fucosylated oligosaccharides in mother's milk alleviate the effects of caesarean birth on infant gut microbiota. Sci. Rep. 2018, 8, 13757, doi:10.1038/s41598-018-32037-6.
21. Borewicz, K.; Gu, F.; Saccenti, E.; Hechler, C.; Beijers, R.; de Weerth, C.; van Leeuwen, S.S.; Schols, H.A.; Smidt, H. The association between breastmilk oligosaccharides and faecal microbiota in healthy breastfed infants at two, six, and twelve weeks of age. Sci. Rep. 2020, 10, 4270, doi:10.1038/s41598-020-61024-z.
22. Henrick, B.M.; Rodriguez, L.; Lakshmikanth, T.; Pou, C.; Henckel, E.; Arzoomand, A.; Olin, A.; Wang, J.; Mikes, J.; Tan, Z.; et al. Bifidobacteria-mediated immune system imprinting early in life. Cell 2021, 184, 3884–3898.e11, doi:10.1016/j.cell.2021.05.030.
23. Ioannou, A.; Knol, J.; Belzer, C. Microbial Glycoside Hydrolases in the First Year of Life: An Analysis Review on Their Presence and Importance in Infant Gut. Front. Microbiol. 2021, 12, 631282, doi:10.3389/fmicb.2021.631282.
24. Marcobal, A.; Barboza, M.; Froehlich, J.W.; Block, D.E.; German, J.B.; Lebrilla, C.B.; Mills, D.A. Consumption of human milk oligosaccharides by gut-related microbes. J. Agric. Food Chem. 2010, 58, 5334–5340, doi:10.1021/jf9044205.
25. Yu, Z.T.; Chen, C.; Kling, D.E.; Liu, B.; McCoy, J.M.; Merighi, M.; Heidtman, M.; Newburg, D.S. The principal fucosylated oligosaccharides of human milk exhibit prebiotic properties on cultured infant microbiota. Glycobiology 2013, 23, 169–177, doi:10.1093/glycob/cws138.
26. Yoshida, E.; Sakurama, H.; Kiyohara, M.; Nakajima, M.; Kitaoka, M.; Ashida, H.; Hirose, J.; Katayama, T.; Yamamoto, K.; Kumagai, H. Bifidobacterium longum subsp. infantis uses two different beta-galactosidases for selectively degrading type-1 and type-2 human milk oligosaccharides. Glycobiology 2012, 22, 361–368, doi:10.1093/glycob/cwr116.
27. James, K.; Motherway, M.O.; Bottacini, F.; van Sinderen, D. Bifidobacterium breve UCC2003 metabolises the human milk oligosaccharides lacto-N-tetraose and lacto-N-neo-tetraose through overlapping, yet distinct pathways. Sci. Rep. 2016, 6, 38560, doi:10.1038/srep38560.
28. Ruiz-Moyano, S.; Totten, S.M.; Garrido, D.A.; Smilowitz, J.T.; German, J.B.; Lebrilla, C.B.; Mills, D.A. Variation in consumption of human milk oligosaccharides by infant gut-associated strains of Bifidobacterium breve. Appl. Environ. Microbiol. 2013, 79, 6040–6049, doi:10.1128/AEM.01843-13.
29. Puccio, G.; Alliet, P.; Cajozzo, C.; Janssens, E.; Corsello, G.; Sprenger, N.; Wernimont, S.; Egli, D.; Gosoniu, L.; Steenhout, P. Effects of Infant Formula With Human Milk Oligosaccharides on Growth and Morbidity: A Randomized Multicenter Trial. J. Pediatr. Gastroenterol. Nutr. 2017, 64, 624–631, doi:10.1097/MPG.0000000000001520.
30. Berger, B.; Porta, N.; Foata, F.; Grathwohl, D.; Delley, M.; Moine, D.; Charpagne, A.; Siegwald, L.; Descombes, P.; Alliet, P.; et al. Linking Human Milk Oligosaccharides, Infant Fecal Community Types, and Later Risk To Require Antibiotics. mBio 2020, 11, e03196-19, doi:10.1128/mBio.03196-19.
31. Yu, X.T.; Wang, L.; Zeng, T. Revisit of Machine Learning Supported Biological and Biomedical Studies. Methods Mol. Biol. 2018, 1754, 183–204, doi:10.1007/978-1 4939-7717-8_11.
32. Kieser, S.; Sarker, S.A.; Berger, B.; Sultana, S.; Chisti, M.J.; Islam, S.B.; Foata, F.; Porta, N.; Betrisey, B.; Fournier, C.; et al. Antibiotic Treatment Leads to Fecal Escherichia coli and Coliphage Expansion in Severely Malnourished Diarrhea Patients. Cell Mol. Gastroenterol. Hepatol. 2018, 5, 458–460 e456, doi:10.1016/j.jcmgh.2017.11.014.
33. Moco, S.; Collino, S.; Rezzi, S.; Martin, F.P. Metabolomics perspectives in pediatric research. Pediatr. Res. 2013, 73, 570–576, doi:10.1038/pr.2013.1.
34. Martin, F.P.; Moco, S.; Montoliu, I.; Collino, S.; Da Silva, L.; Rezzi, S.; Prieto, R.; Kussmann, M.; Inostroza, J.; Steenhout, P. Impact of breast-feeding and high- and low-protein formula on the metabolism and growth of infants from overweight and obese mothers. Pediatr. Res. 2014, 75, 535–543, doi:10.1038/pr.2013.250.
35. Le Cao, K.A.; Boitard, S.; Besse, P. Sparse PLS discriminant analysis: Biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinform. 2011, 12, 253, doi:10.1186/1471-2105-12-253.
36. Team, R.C. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; Vienna, Austria, 2018. https://www.R-project.org/.
37. Max Kuhn, J.W.; Weston, S.; Williams, A.; Keefer, C.; Engelhardt, A.; Cooper, T.; Mayer, Z.; Kenkel, B. The R Core Team, Michael Benesty, Reynald Lescarbeau, Andrew Ziem, Luca Scrucca, Yuan Tang, Can Candan, Tyler Hunt. Caret: Classification and Regression Training, 2017. Available online: https://cran.r project.org/web/packages/caret/index.html (accessed on 10. June 2021).
38. Rohart, F.; Gautier, B.; Singh, A.; Le Cao, K.A. mixOmics: An R package for omics feature selection and multiple data integration. PLoS Comput. Biol. 2017, 13, e1005752, doi:10.1371/journal.pcbi.1005752.
39. Dominguez-Bello, M.G.; Godoy-Vitorino, F.; Knight, R.; Blaser, M.J. Role of the microbiome in human development. Gut 2019, 68, 1108–1114, doi:10.1136/gutjnl-2018-317503.
40. Sela, D.A.; Mills, D.A. Nursing our microbiota: Molecular linkages between bifidobacteria and milk oligosaccharides. Trends Microbiol. 2010, 18, 298–307, doi:10.1016/j.tim.2010.03.008.
41. Zuniga, M.; Monedero, V.; Yebra, M.J. Utilization of Host-Derived Glycans by Intestinal Lactobacillus and Bifidobacterium Species. Front. Microbiol. 2018, 9, 1917, doi:10.3389/fmicb.2018.01917.
42. Katoh, T.; Ojima, M.N.; Sakanaka, M.; Ashida, H.; Gotoh, A.; Katayama, T. Enzymatic Adaptation of Bifidobacterium bifidum to Host Glycans, Viewed from Glycoside Hydrolyases and Carbohydrate-Binding Modules. Microorganisms 2020, 8, 481, doi:10.3390/microorganisms8040481.
43. Sela, D.A.; Chapman, J.; Adeuya, A.; Kim, J.H.; Chen, F.; Whitehead, T.R.; Lapidus, A.; Rokhsar, D.S.; Lebrilla, C.B.; German, J.B.; et al. The genome sequence of Bifidobacterium longum subsp. infantis reveals adaptations for milk utilization within the infant microbiome. Proc. Natl. Acad. Sci. USA 2008, 105, 18964–18969, doi:10.1073/pnas.0809584105.
44. Chichlowski, M.; Shah, N.; Wampler, J.L.; Wu, S.S.; Vanderhoof, J.A. Bifidobacterium longum Subspecies infantis (B. infantis) in Pediatric Nutrition: Current State of Knowledge. Nutrients 2020, 12, 1581, doi:10.3390/nu12061581.
45. Fukuda, S.; Toh, H.; Hase, K.; Oshima, K.; Nakanishi, Y.; Yoshimura, K.; Tobe, T.; Clarke, J.M.; Topping, D.L.; Suzuki, T.; et al. Bifidobacteria can protect from enteropathogenic infection through production of acetate. Nature 2011, 469, 543–547, doi:10.1038/nature09646.
46. Takeuchi, T.; Miyauchi, E.; Kanaya, T.; Kato, T.; Nakanishi, Y.; Watanabe, T.; Kitami, T.; Taida, T.; Sasaki, T.; Negishi, H.; et al. Acetate differentially regulates IgA reactivity to commensal bacteria. Nature 2021, 595, 560–564, doi:10.1038/s41586-021-03727-5.
47. Dang, A.T.; Marsland, B.J. Microbes, metabolites, and the gut-lung axis. Mucosal Immunol. 2019, 12, 843–850, doi:10.1038/s41385- 019-0160-6.
48. Boets, E.; Gomand, S.V.; Deroover, L.; Preston, T.; Vermeulen, K.; De Preter, V.; Hamer, H.M.; Van den Mooter, G.; De Vuyst, L.; Courtin, C.M.; et al. Systemic availability and metabolism of colonic-derived short-chain fatty acids in healthy subjects: A stable isotope study. J. Physiol. 2017, 595, 541–555, doi:10.1113/JP272613. Microorganisms 2021, 9, 1939 18 of 18
49. Antunes, K.H.; Fachi, J.L.; de Paula, R.; da Silva, E.F.; Pral, L.P.; Dos Santos, A.A.; Dias, G.B.M.; Vargas, J.E.; Puga, R.; Mayer, F.Q.; et al. Microbiota-derived acetate protects against respiratory syncytial virus infection through a GPR43-type 1 interferon response. Nat. Commun. 2019, 10, 3273, doi:10.1038/s41467-019-11152-6.
50. Sencio, V.; Barthelemy, A.; Tavares, L.P.; Machado, M.G.; Soulard, D.; Cuinat, C.; Queiroz-Junior, C.M.; Noordine, M.L.; Salome- Desnoulez, S.; Deryuter, L.; et al. Gut Dysbiosis during Influenza Contributes to Pulmonary Pneumococcal Superinfection through Altered Short-Chain Fatty Acid Production. Cell Rep. 2020, 30, 2934 2947.e6, doi:10.1016/j.celrep.2020.02.013.
51. Tsukuda, N.; Yahagi, K.; Hara, T.; Watanabe, Y.; Matsumoto, H.; Mori, H.; Higashi, K.; Tsuji, H.; Matsumoto, S.; Kurokawa, K.; et al. Key bacterial taxa and metabolic pathways affecting gut short-chain fatty acid profiles in early life. ISME J. 2021, 15, 2574– 2590, doi:10.1038/s41396-021-00937-7.
52. Laursen, M.F.; Andersen, L.B.; Michaelsen, K.F.; Molgaard, C.; Trolle, E.; Bahl, M.I.; Licht, T.R. Infant Gut Microbiota Development Is Driven by Transition to Family Foods Independent of Maternal Obesity. mSphere 2016, 1, e00069-15, doi:10.1128/mSphere.00069-15.
53. Shao, Y.; Forster, S.C.; Tsaliki, E.; Vervier, K.; Strang, A.; Simpson, N.; Kumar, N.; Stares, M.D.; Rodger, A.; Brocklehurst, P.; et al. Stunted microbiota and opportunistic pathogen colonization in caesarean-section birth. Nature 2019, 574, 117–121, doi:10.1038/s41586-019-1560-1.
54. Martin, F.; Laursen, M.S.; von Burg, N.; Mörbe, U.; Andersen, D.; Moll, J.M.; Pekmez, C.T.; Rivollier, A.; Kim, F.; Mølgaard, M.C.; et al. Breastmilk-Promoted bifidobacteria produce aromatic amino acids in the infant gut. bioRxiv 2020, https://doi.org/10.1101/2020.01.22.914994.
55. Metidji, A.; Omenetti, S.; Crotta, S.; Li, Y.; Nye, E.; Ross, E.; Li, V.; Maradana, M.R.; Schiering, C.; Stockinger, B. The Environmental Sensor AHR Protects from Inflammatory Damage by Maintaining Intestinal Stem Cell Homeostasis and Barrier Integrity. Immunity 2018, 49, 353–362.e5, doi:10.1016/j.immuni.2018.07.010.
56. Stockinger, B.; Di Meglio, P.; Gialitakis, M.; Duarte, J.H. The aryl hydrocarbon receptor: Multitasking in the immune system. Annu. Rev. Immunol. 2014, 32, 403 432, doi:10.1146/annurev-immunol-032713-120245.
57. Li, Y.; Jiang, Y.; Zhang, L.; Qian, W.; Hou, X.; Lin, R. Exogenous l-fucose protects the intestinal mucosal barrier depending on upregulation of FUT2-mediated fucosylation of intestinal epithelial cells. FASEB J. 2021, 35, e21699, doi:10.1096/fj.202002446RRRR.
58. Pasolli, E.; Truong, D.T.; Malik, F.; Waldron, L.; Segata, N. Machine Learning Meta-analysis of Large Metagenomic Datasets: Tools and Biological Insights. PLoS Comput. Biol. 2016, 12, e1004977, doi:10.1371/journal.pcbi.1004977.