Next Article in Journal
Correction: Wist et al. Phenotypic and Genotypic Traits of Vancomycin-Resistant Enterococci from Healthy Food-Producing Animals. Microorganisms 2020, 8, 261
Next Article in Special Issue
Bovine Milk Oligosaccharides and Human Milk Oligosaccharides Modulate the Gut Microbiota Composition and Volatile Fatty Acid Concentrations in a Preclinical Neonatal Model
Previous Article in Journal
Impact of Chronic Exposure to Sublethal Doses of Glyphosate on Honey Bee Immunity, Gut Microbiota and Infection by Pathogens
Previous Article in Special Issue
The Effects of Genetic Relatedness on the Preterm Infant Gut Microbiota
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 4
10.3390/microorganisms9040846
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mediation Analysis to Identify Links between Gut Bacteria and Memory in Context of Human Milk Oligosaccharides

by
Stephen A. Fleming
1,*,
Jonas Hauser
2,
Jian Yan
3,
Sharon M. Donovan
4,5,
Mei Wang
4 and
Ryan N. Dilger
1,5,6,7
1
Traverse Science, Inc., Champaign, IL 61820, USA
2
Société des Produits Nestlé SA, 1000 Lausanne, Switzerland
3
Nestlé Product Technology Center Nutrition, CH-1800 Vevey, Switzerland
4
Department of Food Science and Human Nutrition, University of Illinois, Urbana, IL 61801, USA
5
Division of Nutritional Sciences, University of Illinois, Urbana, IL 61801, USA
6
Piglet Nutrition and Cognition Laboratory, Department of Animal Sciences, University of Illinois, Urbana, IL 61801, USA
7
Neuroscience Program, University of Illinois, Urbana, IL 61801, USA
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(4), 846; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040846
Submission received: 18 March 2021 / Revised: 12 April 2021 / Accepted: 14 April 2021 / Published: 15 April 2021
(This article belongs to the Special Issue Novel Insights into the Infant Gut Microbiome and Host–Microbe Interactions)
Browse Figures
Review Reports Versions Notes

Abstract

:
Elucidating relationships between the gut and brain is of intense research focus. Multiple studies have demonstrated that modulation of the intestinal environment via prebiotics or probiotics can induce cognitively beneficial effects, such as improved memory or reduced anxiety. However, the mechanisms by which either act remain largely unknown. We previously demonstrated that different types of oligosaccharides affected short- and long-term memory in distinct ways. Given that the oligosaccharide content of human milk is highly variable, and that formula-fed infants typically do not consume similar amounts or types of oligosaccharides, their potential effects on brain development warrant investigation. Herein, a mediation analysis was performed on existing datasets, including relative abundance of bacterial genera, gene expression, brain volume, and cognition in young pigs. Analyses revealed that numerous bacterial genera in both the colon and feces were related to short- and/or long-term memory. Relationships between genera and memory appeared to differ between diets. Mediating variables frequently included GABAergic and glutamatergic hippocampal gene expression. Other mediating variables included genes related to myelination, transcription factors, brain volume, and exploratory behavior. Overall, this analysis identified multiple pathways between the gut and brain, with a focus on genes related to excitatory/inhibitory neurotransmission.

1. Introduction

During early life, neurodevelopment is strongly impacted by nutrition [1,2,3,4,5], among other factors. Of the components in human milk, oligosaccharides provide little to no direct nutritional value to the infant yet are the third most concentrated solid in human milk [6]. Human milk contains a greater variety and higher concentration of oligosaccharides than other mammals [7] and their concentrations are highly variable between mothers [8]. Their functions are multifaceted, impacting immunity [9], mucosal physiology [10], gut bacteria [11], and potentially brain development [12,13,14,15].
Regarding the impact on neurodevelopment, oligosaccharides of multiple origins have shown the capacity to alter human and animal behavior. Rats provided 2′-fucosyllactose (2′-FL) demonstrated vagal-dependent improved learning and memory [16,17], mice consuming chitosan-derived oligosaccharides showed attenuations in neurodegenerative-associated phenotypes [18], and a recent paper associated human milk levels of 2′-FL at 1 month with cognitive function at 24 months in human infants [19]. The potential mechanisms linking the gut and brain are myriad, and include alteration of the central neurotransmitter systems, alteration of the microbiota or microbial-derived metabolites, and direct vagal mediation. For example, the hippocampal serotonergic system is altered in germ-free animals and partially restored by colonization [20], and fecal microbial transplants (FMT) from schizophrenic patients to germ-free mice altered gamma-aminobutyric acid (GABA)-related metabolism in the brain and behavior [21]. Mice that demonstrated depressive-like behavior after a FMT from a non-obese diabetic donor mouse had increased fecal cresol, which has been shown to impair myelination in culture [22]. Short-chain fatty acids from microbial fermentation have been linked to anxiety and aggression [23,24], or conversely proposed as potential therapeutics [25]. Lastly, the vagus nerve mediates improvements to behavior after ingestion of oligosaccharides or probiotics [17,26], as vagotomy prevents pre-/probiotic-related improvements to behavior.
Desire to modulate the microbiota-gut-brain axis has pervaded numerous psychological disciplines: cognitive development [27], depression [28], autism [29,30], schizophrenia [21], and others. The development of the microbiome and its relationship to neurodevelopment is of special interest as the neonatal period represents a critical window when the both the microbiome and brain are sensitive to both distress and eustress, potentially shaping long-term development [31].
We previously demonstrated that consumption of human and non-human oligosaccharides significantly improved performance in a test of recognition memory in pigs [32,33], but such a benefit was specific to the type of oligosaccharide ingested and the type of memory tested (i.e., short or long term). Furthermore, oligosaccharide intake altered average brain volume and hippocampal gene expression [32,33]. Surprisingly, minimal differences between diets in relative abundance of bacterial genera in the colon were found [34]. We identified potential mechanistic relationships between neurotransmitter-related genes and recognition memory; however, these relationships were not consistently found in pigs fed different oligosaccharides. The objective of this study was to perform a secondary analysis on the same datasets to determine whether a relationship exists between the bacteria present in the ascending colon and feces and recognition memory, and, if so, whether that relationship was mediated by gene expression or magnetic resonance imaging outcomes in the context of oligosaccharide intake.

2. Materials and Methods

Procedures have been reviewed in detail in [32,33,34] and are summarized below.

2.1. Animals and Housing

All animal care and experimental procedures were in accordance with the National Research Council Guide for Care and Use of Laboratory Animals and approved by the University of Illinois at Urbana-Champaign Institutional Animal Care and Use Committee. Seventy-two intact male pigs (Camborough® breed, Line 2 and Line 3) were artificially-reared from postnatal day (PND) 2 until PND 33 across six independent cohorts (n = 12 per cohort, each cohort separated in time). Pigs were randomized to each group such that litter representation and initial bodyweight (BW) were counterbalanced within each cohort and between each dietary group, respectively. All pigs were housed in master caging units that contained six individual stainless-steel cages (L × W × H of 87.6 × 88.9 × 50.8 cm) with clear, polycarbonate facades on three sides of the cage and vinyl-coated, expanded-metal flooring (Tenderfoot®, Minneapolis, MN, USA). The flooring allowed fecal and urinary output to drop to a collection system below the pan, thereby reducing the opportunity for fecal contamination of the diets or coprophagy. Pigs were able to see, smell, hear, and minimally touch one another. A towel and toy were included in each cage to provide enrichment, and all pigs were removed from cages and allowed to socialize with each other for approximately 30 min each day. From subjective observation, minimal transfer of fecal material between groups occurred during socialization.
All pigs were reared in the same room with ambient temperature maintained between 27 and 29 °C and a 12 h light/dark cycle maintained from 600 to 1800 h. Prior to placement in the artificial rearing system, pigs were administered 5.0 mL of Clostridium perfingens antitoxin C + D per the manufacturer’s recommendations (Colorado Serum Company, Denver, CO, USA), a standard veterinary procedure to prevent enterotoxemia [35]. At study conclusion (PND 33), pigs were anesthetized using a telazol:ketamine:xylazine solution (50.0 mg tiletamine plus 50.0 mg of zolazepam reconstituted with 2.50 mL ketamine [100 g/L] and 2.50 mL xylazine [100 g/L]; Fort Dodge Animal Health) by intramuscular injection at 0.03 mL/kg BW. After anesthetic induction, pigs were euthanized via intracardiac administration of sodium pentobarbital (86.0 mg/kg of body weight; Euthasol, Virbac Animal Health, Fort Worth, TX, USA). Two pigs (from the BMOS + HMO group) were removed from study due to failure-to-thrive (i.e., exhibited very low growth).

2.2. Dietary Treatments

Pigs (n = 12 per diet) were provided milk replacers reconstituted at 200 g of dry powder per 800 g of water. Reconstituted diets were analyzed to contain approximately 0 g/L oligosaccharide (OS) (control [CON], ProNurse® Specialty Milk Replacer, Purina Animal Nutrition, Gray Summit, MO, USA), 5.79 g/L bovine milk oligosaccharides (BMOS, Nestlé Product & Technology center, Konolfingen, Switzerland), 1.23 g/L HMO (HMO, 0.81 g/L of 2′-fucosyllactose [2′-FL] + 0.42 g/L of Lacto-N-neotetraose [LNnT], Glycom, Hørsholm, Denmark), both bovine and human milk oligosaccharides (BMOS + HMO; 5.75 g/L of BMOS + 1.00 g/L of 2′-FL + 0.53 g/L of LNnT), 3.62 g/L oligofructose (OF, Orafti® P95, Beneo-Orafti, Tiene, Belgium), or 3.41 g/L OF + 1.12 g/L 2′-FL (OF + 2′-FL). Briefly, the BMOS product was derived from bovine whey and composed primarily of galactooligosaccharide and trace amounts of 3′- and 6′-sialyllactose. All diets were formulated with the addition of lactose, such that each diet contained the same amount of total carbohydrate (See Table 1). The nutritional composition of the base formula and oligosaccharide content have been previously reported [32,33]. Though concentrations of each oligosaccharide differed, confounding the ability to differentiate between dose and type of oligosaccharide, the concentrations were chosen to remain consistent with previously conducted clinical trials on BMOS [36,37,38], 2′-FL [39,40], and OF [41,42]. Ultimately, the difference in lactose between diets minimally contributed to the metabolizable energy and growth was equivalent between groups [32,33].
Pigs received approximately 500 mL of experimental diets on the day of arrival to the rearing facility and were fed at a rate of 285 mL and 325 mL of reconstituted milk replacer per kg BW from PND 3-6 and PND 7-33, respectively. Bodyweight was recorded daily to accurately dose meals, which were administered 10 times per day, approximately every 100 min, between 1000 h and 0400 h using an automated feeding system. Feed refusals were not quantified. All pigs were allowed ad libitum access to water at all times.

2.3. Behavior

Pigs were tested on the novel object recognition (NOR) task using two different delays to assess short- and long-term recognition memory. Methods used were adapted from previous studies using this task in pigs from other labs [43,44,45,46], our own lab [47,48,49], and are described from this study [32,33]. Testing consisted of a habituation phase, a sample phase, and a test phase. During the habituation phase, each pig was placed in an empty testing arena for 10 min each day for two days leading up to the sample phase. In the sample phase, the pig was placed in the arena containing two identical objects and given 5 min for exploration. After a delay of 1 or 48 h (representing a short or long delay, respectively), the pig was returned to the arena for the test phase of the NOR task. During the test phase, the pig was placed in the arena containing one object from the sample phase and a novel object and allowed to explore for 5 min. Habituation trials began at PND 22 and testing on the sample phase began on PND 24. Recognition index, or the proportion of time spent with the novel object compared to total exploration of both objects, was used to measure recognition memory.

2.4. Magnetic Resonance Imaging (MRI)

All pigs underwent MRI procedures at PND 32 at the Beckman Institute for Advanced Science and Technology Biomedical Imaging Center using a Siemens MAGNETOM Trio 3T equipment with a Siemens 32-channel head coil. Methods described here are adapted from previous studies using MRI in pigs [50,51,52] and are published [32,33]. Each pig underwent imaging protocols only once. The pig neuroimaging protocol included three magnetization prepared rapid gradient-echo (MPRAGE) sequences and diffusion tensor imaging (DTI) to assess brain macrostructure and microstructure, respectively, as well as magnetic resonance spectroscopy (MRS) to obtain brain metabolite concentrations. In preparation for MRI procedures, anesthesia was induced using an intramuscular injection of telazol (50.0 mg of tiletamine plus 50.0 mg of zolazepam reconstituted with 5.0 DI water; Zoetis, Florham Park, NJ, USA) administered at 0.07 mL/kg BW, and maintained with inhalation of isoflurane (98% O2, 2% isoflurane). Pigs were immobilized during all MRI procedures. Visual observation of each pig’s well-being, as well as observations of heart rate, PO2 and percent of isoflurane were recorded every 5 min during the procedure and every 10 min post-procedure until animals recovered. Total scan time for each pig was approximately 60 min.

2.5. Hippocampal Gene Expression

As previously described [32,33], relative mRNA copy numbers on 93 genes in the hippocampus were quantified using the NanoString nCounter™ system (NanoString Technologies Inc., Seattle, WA, USA) according to the manufacturer’s instructions using 100 ng of RNA as the starting material. Using nSolver software (Version 4.0, NanoString Technologies Inc., Seattle, WA, USA), background subtraction using the median of all eight negative controls was followed by positive control normalization using the geometric mean of six positive controls and housekeeping normalization using the geometric mean of six housekeeping genes (Ribosomal Protein L19,RPL-19; Ribosomal Protein S18, RPS18; Actin β, ACTB; Tubulin β Class 1, TUBB; Glyceraldehyde-3-Phosphate Dehydrogenase, GAPDH; β-glucoronidase, GUS).

2.6. DNA Extraction

DNA was extracted from ascending colon (AC, henceforth interchangeably referred to as colonic) contents and feces by the QIAamp Fast DNA Stool Mini Kit (Qiagen, Valencia, CA, USA) in combination with bead beating on the FastPrep-24 System (MP Biomedicals, Carlsbad, CA, USA) as previously described [53]. DNA concentration was measured with a NanoDrop 1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA).

2.7. PCR Amplification and Sequencing of 16S rRNA Genes

PCR amplification and sequencing of 16S rRNA genes were performed at the DNA Services Lab at the University of Illinois at Urbana-Champaign according to methods previously described [54]. Briefly, the V3-V4 region of bacterial 16S rRNA genes were amplified using primers V3f (5′-CCTACGGGAGGCAGCAG-3′) and V4r (5′-GGACTACHVGGGTWTCTAAT-3′) as previously described. The amplicons were mixed in equimolar concentrations and sequenced on an Illumina MiSeq flow cell for 251 cycles from each end of the fragments using a MiSeq 500-cycle sequencing kit version 2 (2 × 250 nt paired end reads).

2.8. Sequence Processing

Sequences were demultiplexed at the sequencing facility with the bcl2fastq v2.17.1.14 Conversion Software (Illumina, San Diego, CA, USA), and allowed 0 mismatches in the barcode sequences. De-multiplexed forward (read 1) and reverse reads (read 2) were processed using the QIIME software package as described previously [54,55]. Briefly, the paired-end reads were merged, quality filtered, and split into libraries at Phred quality score ≥ 25. Operational taxonomic unit (OTU) assignment, representative sequence picking, and chimera removal were performed as described by Monaco et al., [54]. An OTU table was created and further filtered to remove non-aligned and chimeric OTUs and singletons [54]. The representative sequence of each OTU was assigned to different taxonomic levels based on Greengenes taxonomy (gg_13_8) with the use of the Ribosomal Database Project naïve Bayesian rRNA Classifier at 80% confidence level [56].

2.9. Statistical Analysis

2.9.1. Data Cleaning and Variable Selection

Genes expressed below threshold (defined as the mean 2 standard deviations below the negative controls) were omitted from analysis. Microbial genera with a median relative abundance of less than 0.05% were omitted from analysis on a group-wise basis (i.e., if genera median abundance was above 0.05% in a dietary group, it was retained for that group only). Given extremely high multicollinearity, DTI measures of axial, medial, and radial diffusivity were omitted from the analysis, with only fractional anisotropy retained as the primary DTI outcome. All variables analyzed are shown in Table S1. Missing values were imputed with the median of that group. Values of zero or observations that were not quantified (e.g., non-compliance during behavioral testing) were not considered as “missing” and were not imputed. Observations that were not quantified were not included in the analysis. After cleaning, a total of 70 pigs and 289 variables remained. The dataset contained 19,186 observations, of which 46 were not quantified (0.24%), and 998 were missing (5.2%). For a breakdown of sample size and missingness per group and variable, see Table S2.

2.9.2. Mediation Analysis

All data cleaning and statistics were conducting using R software [57]. Ordinary Least Squares regression coefficients were constructed using the stats package. Mediation analyses used the R package “Mediation” [58], with a bootstrap sample size of 2000 and 95% confidence interval estimates constructed using the percentile method. Predictor variables included microbial genera from both the ascending colon (62 variables) and fecal samples (63 variables). Here, the same genera from a different sampling location were included as separate variables. Mediating variables included gene expression (73 variables), MRI (60 variables), and behavioral variables (28 variables) of exploration during the novel object recognition task. The predicted variable was the recognition index from either a short- or a long-delay. All mediations analyses were performed on each diet group by using a sub-group analysis. All variables were rank transformed prior to mediation. A depiction of the model and equations are shown in Figure 1. All heatmap visualizations were created using the package “ComplexHeatmap” [59].

3. Results

A descriptive analysis of bacterial composition at the phylum and genus level is provided in Figure 2. Given this is a secondary analysis, reports of the impact of the diet on 16S sequencing outcomes, gene expression, MRI, and behavior are reported elsewhere [32,33,34]. Due to the volume of data generated, only mediations where the confidence interval of the indirect effect did not include zero are discussed.

3.1. Short-Term Memory

No variables were identified as mediators within the CON or BMOS groups. Mediating variables in the HMO group were dominated by two types: GABA-related (GAD, glutamate decarboxylase; SLC32A1, vesicular GABA transporter [VGAT]) and cortical volume (Figure 3). With the exception of ascending colon Lactobacillus relative abundance, all genera were from fecal samples. Mediating variables in the BMOS + HMO group were mostly volume of the pons or novel object visit time, which were related to the relative abundance of Prevotella, Lactobacillus, and Desulfovibrio in the ascending colon. In the OF group, whole brain volume mediated the relationship between fecal Prevotella and short-term memory.
In pigs fed OF + 2′-FL, mediators included several different neurotransmitter-related genes (5HTR2, Serotonin Receptor 2A; CHRNA7, Cholinergic Receptor Nicotinic Alpha 7 Subunit; GABRA5, GABA receptor Subunit Alpha-5; GRIA2, Glutamate Ionotropic Receptor AMPA Type Subunit 2; GRIN2D, NMDAR2D), and mediated relationships between colonic Prevotella, Anaerotruncus, Megasphaera, and short-term recognition memory.

3.2. Long-Term Memory

No mediating variables between bacterial genera and long-term memory were identified in the CON and OF + 2′-FL groups (Figure 4). In the other groups, a few variables tended to mediate multiple pathways. In the BMOS group, myelin-related genes (MAG, myelin-associated glycoprotein; PLP, proteolipid protein; MBP, myelin basic protein) mediated the relationship between fecal Ruminococcus, Escherichia, colonic Butyricicoccus, and long-term memory. Transcription-related genes (CREBBP, cyclic AMP response element binding protein [CREB] binding protein; NR4A2, nuclear receptor-related 1 protein [NURR1]) acted as mediators between long-term memory and colonic Alistipes, and fecal Ruminococcus, Synergistes, Blautia, unclassified Veillonellaceae, and Megasphaera.
In the HMO group, total distance moved in the novel object recognition task acted as a mediator for more predictors than any other variable. Predictor variables tended to be colonic genera (eight variables), rather than fecal genera (three variables). Additionally, GABRD (GABA type A receptor delta subunit) and relative volume of the corpus callosum mediated the relationship between colonic Eubacterium and fecal unclassified Paraprevotellaceae.
In pigs fed BMOS + HMO, mediating variables were dominated by GABA/Glu-related genes (GABRA2, GABA type A receptor alpha 2 subunit; GRIN2B, NMDARD2B/GluN2B; SLC17A6, vesicular glutamate transporter 2 [VGLUT2]; SLC17A8, vesicular glutamate transporter 3 [VGLUT3]), or mean novel object visit time. In the OF group, absolute volume of the caudate mediated the relationship between colonic Bilophila and long-term memory.

3.3. Suppressor Variables

Numerous mediating variables acted as suppressor variables, as their inclusion in the model increased the strength of the relationship between host bacteria and memory (Figure S1). These variables themselves were weakly associated with memory. They included measures of exploratory behavior (mean visit time, total distance moved), gene expression (mostly related to GABA/Glu or myelination), fractional anisotropy, magnetic resonance spectroscopy, or brain volume. These variables were present across a range of dietary groups for both short- and long-term memory, and nearly half of them were related to MRI outcomes.

4. Discussion

The goal of the present analysis was to uncover the complexity of the relationships between the previously observed findings linking dietary oligosaccharides and neurocognitive development [32,33], with a focus on gut bacteria. The relationships between each of these variables were not clear, however. For example, the human milk oligosaccharides 2′-FL and LNnT improved short-term, but not long-term memory, and likewise a combination of human and bovine milk oligosaccharides improved long-term, but not short-term memory [33]. Furthermore, the effects of these oligosaccharides on brain volume did not follow an easily discernible pattern. However, we did find that the relationship between gene expression and memory was dependent on the type of oligosaccharide consumed. Given that oligosaccharides have been shown to impact the composition and/or metabolic activity of the gut microbiota, often simultaneously with an improvement in behavioral outcomes, we hypothesized that gut bacteria may impact cognition by acting through a mediating variable.
Here, we combined datasets containing relative abundance of colonic/fecal bacteria, brain volume and chemistry, gene expression, and exploratory behavior using a mediation analysis to identify potential paths between the host microbiota and cognition. At a high level, we found that mediating variables varied with oligosaccharide intake, but tended to appear in similar patterns. Importantly, even though the median relative abundance of genera were similar between groups, the relationships these genera exhibited with mediating and predictor variables appeared to differ between groups.
Prior to discussing the biological importance of (or lack thereof) these mediations, several limitations to the present analysis should be addressed. Given the presence of high skew and outliers in several of the datasets, especially 16S sequencing data, we chose to rank transform all variables prior to regression. Hauke and Kossowski provide a helpful explanation on the consequences of rank transformation in context of Pearson and Spearman correlation [60]. While the use of Spearman correlation (and by extension linear regression on ranked variables) can lessen the impact of outliers and skewness while simultaneously allowing detection of monotonic relationships between 2 variables, it does have several limitations. Primarily, loss of information regarding variability and the spacing between variables: the distance between each observation is equal in a ranked variable, whereas in reality the distance between observations may be extremely variable. Furthermore, the data become less reproducible. In other words, future research will be unable to replicate our exact findings as the rank of a specific observation is only relevant within context of the raw data collected in this study. However, we felt the benefit of being able to identify both linear and non-linear relationships, reduced weighting of outliers, and reduced impact on skewness outweighed these drawbacks.
Herein, we used a sub-group analysis as an alternative to moderated mediation to identify the impact of a multicategorical variable (oligosaccharide intake across 6 different diets) on a mediation. Although a multicategorical, moderated mediation would be preferable to model, not all variables were equally present between groups (e.g., some microbial genera were completely absent in one group and abundant in another). Furthermore, attempts to create a conditional process model [61] were precluded by the inadvertent creation of rank-deficient matrices or singular inverse matrices. Though not clearly identified, these issues may have arisen due to high multicollinearity and interactions between predictor variables in the model that centering was unable to overcome. To avoid these issues a sub-group analysis was used, which presents its own limitations. As described by Hayes [62], such a method should not be used when: (1) the hypothesis is specific to moderation of a specified path; (2) a causal steps approach is used (i.e., mediation is only estimated after a relationship between X and Y, X and M, or M and Y are established); (3) sample size between groups differs; and (4) the moderating variable is not categorical. Here, we did not state an a priori hypothesis regarding which path would be moderated, did not perform a causal steps analysis, and the moderator chosen (diet) was categorical. Thus, we violated the third criterion, as not all groups contained equal sample size, creating an inequality in power to detect a mediation between groups. Therefore, the presence or absence of a mediation in one group compared to another should not be interpreted definitively that the mediation was moderated at the group level.
We chose not to interpret mediations as “partial” or “full”, due in part to the inequality in power between groups, as smaller sample sizes are more likely to result in the detection of full rather than partial mediations, artificially inflating the probability that a full mediation would be observed. Furthermore, despite conducting numerous mediations, we chose not to use a false discovery rate correction. Our reasoning was as follows: (1) small sample size resulted in lower power, increasing the probability of detecting false negatives (i.e., high likelihood that true mediating variables were undetected); and (2) adjusting for multiple comparisons would further increase the probability of detecting false negatives. Rather, our opinion is that in a high-level exploratory analysis, significant results should be interpreted in context of the literature and followed by replication and greater statistical power to confirm whether the observed effects are likely to have been false. With these limitations in mind the results should be interpreted in context of what can reasonably (i.e., theoretically and practically) be inferred from the present analysis.

4.1. Evidence of Common Paths to Cognition

4.1.1. GABA and Glutamatergic Mediators

Perhaps one of the more surprising findings was the convergence of multiple bacterial genera on a few classes of mediators. For example, in the HMO group, the relationship between bacterial genera and short-term memory was mediated by two GABAergic genes (GAD and SLC32A1) and both the absolute and relative volumes of the left and right cortices (Figure 3). Genera (mostly fecal, some colonic) such as Lactobacillus, Oscillospira, and unclassified Ruminococcaceae, Christensenellaceae, and Elusimicrobiaceae converged onto a combination of these two types of mediators. Furthermore, the direction of the relationships was consistent between models. All relationships with cortical volumes (paths a, b, and c) were positive. Genera mediated by GABA-related genes were all positively related to short-term memory and inversely related to GAD or SLC32A1, which themselves were inversely related to memory. In short, downregulation of GABA neurotransmission-related genes and increased cortical volume mediated the positive relationship between short-term memory and relative abundance of Lactobacillus, Oscillospira, and unclassified Ruminococcaceae, Christensenellaceae, and Elusimicrobiaceae. A similar phenomenon was found for the BMOS + HMO group regarding long-term memory.
In the BMOS + HMO group, mediating variables included GABRA2, GRIN2B, SLC17A6, SLC17A8, and average visit time to the novel object. Here, all mediators except for SLC17A6 positively predicted long-term recognition memory. Each of these mediated more than one bacterial genus. Fecal Blautia was mediated by GABRA2, GRIN2B, and average visit time to the novel object. Fecal Dorea was mediated by GABRA2 and SLC17A8. Fecal Acidaminococcus was mediated by GABRA2 and SLC17A6. Here, lower relative amounts of Blautia and higher relative amounts of Dorea and Acidaminococcus were related to improved long-term memory.
It was unsurprising that GABRA2, GRIN2B, SLC17A6, and SLC17A8 would be related to recognition memory. Each of these has been shown in some fashion to play a role in learning and memory. The expression of GABRA2 in the frontal cortex is associated with cognitive decline [63]. Inhibition of GRIN2B impairs memory [64], reduces spine density [65], and it’s expression decreases with age [66]. The role of SLC17A6 is less clear, however, SLC17A6 is highly expressed during early development in subcortical structures [67,68], and age-impaired rats demonstrate increased levels of SLC17A6 as compared to younger animals, with increased levels of SLC17A6 related to cognitive deficit and lower expression of glutamatergic receptors [69]. Lastly, SLC17A8 deletion has been shown to induce anxiety related behavior [70] and mild impairments in learning and memory [71]. The relationships of genera to these mediating variables suggests that the excitatory/inhibitory balance is an important nexus through which host bacteria can influence behavior.
These results echo those shown by Bravo et al., [26] who demonstrated that supplementation with L. rhamnosus reduced expression of GABA receptors in subcortical regions such as the hippocampus and amygdala simultaneous with decreases in anxiety- and depression-like behaviors. Yet, overall, there are little data linking microbial genera to glutamatergic neurotransmission in a neurotypical context, highlighting the novelty of these findings. Blautia has been shown to be inversely associated with sociability and positively correlated with repetitive and anxiety-like behaviors in a Spearman correlation that pooled control mice in a model of autism spectrum disorder (BTBR T+ Itpr3tf/J mice) [29]. In that model, a reduction in abundance of Blautia was found, concurrent with impairments in bile acid and tryptophan metabolism (despite correlation analysis revealing that reduced Blautia was related with improved sociability and anxiety outcomes). Another study investigating autism spectrum disorder (ASD) in children found reductions in Dorea formicigenerans and Blautia luti [30], but these genera were not related to serotonin, tryptophan, or cytokine-related outcomes. In a Shank3 knock-out (KO) mouse model of ASD, Shank3 KO mice demonstrated reduced hippocampal gene expression of GABRA1, GABRA2, GABRB1 and abundance of L. reuteri. A positive relationship between these GABAergic genes and L. reuteri was found. In a follow-up, Shank3 KO mice supplemented with L. reuteri showed increased hippocampal and prefrontal cortex gene expression of GABAergic receptors in males and females [72]. Clinical data from a large cohort study (Flemish Gut Flora Project, n = 1054) demonstrated that microbial metabolic pathways related to GABA synthesis and glutamate degradation were increased and decreased in subjects with depression, respectively. However, pathways most related to quality of life included DOPAC synthesis, isovaleric acid synthesis, and histamine synthesis [28].
In addition to glutamatergic-related genes, we also observed multiple receptors (GABRA5, 5HTR2, CHRNA7, GRIA2, GRIN2D) for different neurotransmitters (GABA, serotonin, acetylcholine, glutamate) that were all inversely related to short-term memory, and mediated relationships with genera in the colon such as Prevotella, unclassified Prevotellaceae, Anaerotruncus, and Megasphaera in the OF + 2′-FL group. Here, the positive relationship between colonic Prevotella and short-term memory was mediated by its inverse relationship to all of the aforementioned receptors. In contrast to a “convergence” upon similar mediators, here Prevotella diverged upon multiple different genes encoding receptors with a variety of functions, yet all bore the same inverse relationship to short-term memory.

4.1.2. Myelination and Transcription Factors

Up to this point, the discussion has centered around GABA/Glu-related gene expression. However, a combination of myelin- and transcription-related genes were repeatedly found as mediators of multiple genera and long-term memory in pigs fed BMOS. Of the mediators identified (MAG, MBP, PLP, NR4A2, and CREBBP, all of which were inversely related to memory), Ruminococcus and long-term memory were mediated by all of them. There appear to be relatively little data linking myelination and gut bacteria, as compared to data linking gut bacteria to various neurotransmitter systems. Germ-free mice exhibit hyper-myelinated axons in the prefrontal cortex [73], whereas FMT from nonobese diabetic mice to C57BL/6 mice reduced prefrontal cortex myelination [22]. A human/animal trial has provided some mechanistic evidence linking the gut microbiome with myelination. Pre-term infants were categorized as “low-growth” or “high-growth” according to body weight and FMT from these infants were conducted on eight- to nine-week-old germ-free pregnant mice. Pups from germ-free mice colonized with feces from “low-growth” preterm infants demonstrated reduced NeuN, neurofilament-L, and MBP in cortical homogenates as compared to those born from dams colonized with feces from “high-growth” infants [74]. Reduction of these proteins was age dependent, as MBP was impaired at 2 weeks, but NeuN and NFL were only impaired at 4 weeks. At 4 weeks, germ-free and “low-growth” mice tended to show increased expression of glutamatergic-, GABAergic-, and ion channel-related genes, with downregulation of serotonergic and dopaminergic genes in total brain homogenate. In part, these data corroborate our own findings. We found that reduced GABAergic- and myelin-related genes, but increased expression of glutamatergic genes, was related to improved recognition memory in a diet-dependent manner. While the direction of change (whether increasing GABAergic/myelin gene expression is “good”) differs, both studies have centered around these candidates.
Even less is known regarding transcriptional regulation in the brain and host microbiome. We found CREBBP and NR4A2 to mediate relationships between several genera and long-term memory but are unaware of existing data supporting this finding.

4.1.3. Exploratory Behavior

With respect to long-term memory in the HMO group, total distance moved during the task was a mediator for 12 colonic and fecal genera, four of which were negatively related to long-term memory, eight of which were positively related to recognition memory (Figure 4). Otherwise, the only other mediators included GABRD and relative volume of the corpus callosum, which mediated the relationships between colonic Eubacterium and fecal unclassified Paraprevotellaceae, respectively. Thus, these genera likely did not directly affect recognition memory, but affected general exploratory behavior, which itself is closely related to memory. Conversely, it may be likely that the reverse is true, where highly exploratory animals were more likely to demonstrate long-term memory and have higher/lower abundance of these genera.
Several studies have shown relationships between the gut microbiome and locomotor activity (typically measured as distance moved during a behavioral task). Mice provided an oral antibiotic show slightly less locomotor activity compared to baseline [22], germ-free mice or mice humanized via FMT from schizophrenic patients show increased locomotor activity and reduced anxiety-like behavior [21,75], zebrafish given a probiotic demonstrate increased distance moved during shoaling-related behaviors [76], and pigs fed HMO demonstrated reduced distance moved during the habitation trial of the NOR task [33]. Clearly, both insults to and developmental support of the microbiome result in various effects on locomotor activity across species, which must be interpreted in context of the present study. Here, although pigs fed HMO demonstrated the least movement [33], within that group pigs with the highest locomotor activity demonstrated improved long-term recognition memory.
Average time per visit to the novel object was also a mediator between genera in the BMOS + HMO group for both short- and long-term memory. Here, visits that lasted longer tended to correlate with improved performance. The genera this variable mediated differed between short- and long-term memory and were both negatively and positively related to average visit time. Previously, we found that BMOS + HMO improved long- but not short-term memory [33].

4.1.4. Brain Volume

Volume of the brain (and its subregions) was rarely a mediating a variable between genera and long-term memory (with the exception of caudate volume in pigs fed OF and relative volume of the corpus callosum in pigs fed HMO). For short-term memory, volume of the cortices, pons, and whole brain were mediators in the HMO, BMOS + HMO, and OF groups, respectively. Genera within these models included Lactobacillus, Prevotella, and Desulfovibrio. It was clear, however, that cortical volume (both absolute and relative) was a frequent mediator for multiple genera including Oscillospira, Lactobacillus, and unclassified genera of Ruminococcaceae, Christensenellaceae, and Elusimicrobiaceae in the HMO group. For each mediation, the strength of the mediating variable on memory was moderate (for every rank increase in brain volume, short-term memory increased by approximately 0.4–1 ranks). This aligns with the observation that pigs fed HMO (in the HMO or BMOS + HMO group) demonstrated increased relative volumes of the left and right cortex. Despite this finding, it is unclear whether structural volume as a mediator would persist into later life.
Carlson et al., [27] investigated the relationship between infant microbiome, cognition, and neuroimaging outcomes. Infants were grouped into three clusters based on cluster analysis of the microbiome. These groups tended to cluster around relative abundance of Faecalibacterium, Bacteroides, and an unclassified genus of Ruminoccocaceae. At 2 years of age, scores on the Mullen Scales of Early Learning were different between clusters and differed by receptive language and expressive language, with the group clusters around Bacteroides performing the best. Some associations were found between alpha diversity and regional gray matter volume, but most regions measured were not different between clustered groups. The authors suggest that the gut microbiome may have minimal effects on regional brain volume by 1–2 years of age. When comparing germ-free to specific pathogen free (SPF) mice, SPF mice tend to show higher regional volume, fractional anisotropy, macromolecular proton fraction (marker of myelination), and myelin staining in numerous regions [77]. Yet, these differences largely disappeared between 4 and 12 weeks, while behavioral differences persisted, further confirming the hypothesis by Carlson et al., [27] that structure of the brain is minimally impacted by the gut microbiome. Furthermore, we observed no evidence of spectroscopic or fractional anisotropy measures as mediating variables. Based on the number of mediating variables that included genes, it would appear that molecular function better describes the relationship between cognition and gut bacteria than does structure or metabolic composition of the brain.

4.2. Potential Mechanisms of Action of Prebiotics

The objective of this analysis was to identify potential pathways between gut bacteria and cognition and likely mechanistic paths between the two. Ultimately, the goal was to identify a mechanistic path that might explain the differential effects that various types of oligosaccharides have on their host. However, while new data were generated, as usual in most scientific research, more questions were raised than answered. The clearest finding however was that the action of microbial genera on memory was frequently associated with GABAergic and glutamatergic genes.
Several mechanisms have been put forth linking the gut and brain. These include vagal activity, the neuroimmune axis, and circulation of microbial metabolites [78]. The present data do not suggest any of these routes are inaccurate. However, it may be that glutamatergic activity is modulated by a combination of these factors. It is likely that numerous mediators may be identified between the gut and brain, however the question remains as to which are the most likely candidates. Our findings support the hypothesis that the human milk oligosaccharides used in this study act via alteration of hippocampal expression of GABAergic and glutamatergic genes, most likely in turn impacting long-term potentiation to improve memory, and long-term potentiation itself is what may be mediated by another pathway [17,26].
Beyond shedding light on potential biomarkers, the present analysis highlights the importance of including groups in statistical modelling (whether as a sub-group analysis or as categorical predictors in the model). Given the presence of mediators in some groups and not others, it is likely that they would not have been identified if correlation or regression analyses were conducted irrespective of treatment level. Furthermore, the identification of multiple suppressor variables suggests that simple correlation matrices might result in a significant number of false negatives, or at least underestimate the true relationship between predictor and dependent variables.

4.3. Limitations

Important limitations of the current analysis should be taken into consideration. First, all data were ranked prior to analysis. As described previously, this sacrifices reproducibility across studies but provides the benefit of attenuating non-normality, skew, and the weight of outliers. Second, a sub-group analysis was performed, rather than performing a moderated mediation analysis. Combined with the lack of equal sample size across models a sub-group analysis may lead to the inaccurate assessment that a lack of mediators in one group versus another represents a true biological difference. The data should be viewed in context of these limitations, with the understanding that other methods of modelling the data may reveal different insights. Lastly, 16S sequencing, mRNA copy number, and structural MRI outcomes do not necessarily equate to relevant biological activity. Similarly, mRNA copy number and MRI outcomes do not directly equate to biological activity However, these data present compelling evidence that may be used to justify further, more rigorous, and well-controlled examinations into the gut–brain axis.

5. Conclusions

We conducted a sub-group mediation analysis between gut bacteria, brain outcomes, and cognition. Gene expression of GABA and glutamate-related neurotransmission frequently mediated the relationship between gut bacteria and both short- and long-term memory. Additional mediators included myelination-related genes, transcription factors, brain volume, and exploratory behavior. Importantly, mediating variables were not equally present in all dietary groups. These data provide a step in the understanding of pathways mediating the association between dietary intake of prebiotics in early life and memory function. In particular, these data provide further support for a key role of hippocampal GABAergic and glutamatergic genes in the link between gut and brain development. These data should help identify potential biomarkers and mechanistic pathways for future research.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms9040846/s1, Figure S1: Suppressor variables between gut bacteria and memory, Table S1: Variables included in the mediation analysis, Table S2: Sample size and missingness.

Author Contributions

conceptualization, S.A.F., J.H., J.Y. and R.N.D.; methodology, S.A.F.; software, S.A.F.; validation, S.A.F.; formal analysis, S.A.F.; investigation, S.A.F., S.M.D., M.W. and R.N.D.; resources, S.M.D. and R.N.D.; data curation, S.A.F. and M.W.; writing—original draft preparation, S.A.F.; writing—review and editing, all authors; visualization, S.A.F.; supervision, J.H., S.M.D. and R.N.D.; project administration, S.A.F., J.H. and R.N.D.; funding acquisition, J.H. and J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Société des Produits Nestlé.

Data Availability Statement

Data may be provided upon request.

Conflicts of Interest

This study was funded by the Société des Produits Nestlé, of which the authors Jonas Hauser and Jian Yan were employees of at the time the work was conducted. Stephen Fleming, Ryan Dilger, and Sharon Donovan have consulted for and received grant funding from the Société des Produits Nestlé. In addition, both Stephen Fleming and Ryan Dilger share ownership in Traverse Science, a consulting firm. Mei Wang received consulting funding from Traverse Science to assist with the interpretation of the manuscript.

References

  1. Georgieff, M.K. Long-term Brain and Behavioral Consequences of Early Iron Deficiency. Nutr. Rev. 2011, 69, S43–S48. [Google Scholar] [CrossRef] [Green Version]
  2. Zeisel, S.H. Choline: Needed for Normal Development of Memory. J. Am. Coll. Nutr. 2000, 19, 528S–531S. [Google Scholar] [CrossRef]
  3. Mudd, A.T.; Dilger, R.N. Early-life nutrition and neurodevelopment: Use of the piglet as a translational model. Adv. Nutr. 2017, 8, 92–104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Prado, E.L.; Dewey, K.G. Nutrition and brain development in early life. Nutr. Rev. 2014, 72, 267–284. [Google Scholar] [CrossRef] [Green Version]
  5. Schneider, N.; Garcia-Rodenas, C.L. Early nutritional interventions for brain and cognitive development in preterm infants: A review of the literature. Nutrients 2017, 9, 187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bode, L. Human milk oligosaccharides: Every baby needs a sugar mama. Glycobiology 2012, 22, 1147–1162. [Google Scholar] [CrossRef] [Green Version]
  7. Urashima, T.; Saito, T.; Nakamura, T.; Messer, M. Oligosaccharides of milk and colostrum in non-human mammals. Glycoconj. J. 2001, 18, 357–371. [Google Scholar] [CrossRef] [PubMed]
  8. McGuire, M.K.; Meehan, C.L.; McGuire, M.A.; Williams, J.E.; Foster, J.; Sellen, D.W.; Kamau-Mbuthia, E.W.; Kamundia, E.W.; Mbugua, S.; Moore, S.E.; et al. What’s normal? Oligosaccharide concentrations and profiles in milk produced by healthy women vary geographically. Am. J. Clin. Nutr. 2017, 105, 1086–1100. [Google Scholar] [CrossRef]
  9. Jeurink, P.V.; van Esch, B.C.; Rijnierse, A.; Garssen, J.; Knippels, L.M. Mechanisms underlying immune effects of dietary oligosaccharides. Am. J. Clin. Nutr. 2013, 98, 572S–577S. [Google Scholar] [CrossRef] [Green Version]
  10. Donovan, S.M.; Comstock, S.S. Human Milk Oligosaccharides Influence Neonatal Mucosal and Systemic Immunity. Ann. Nutr. Metab. 2016, 69 (Suppl. 2), 42–51. [Google Scholar] [CrossRef] [PubMed]
  11. Asakuma, S.; Hatakeyama, E.; Urashima, T.; Yoshida, E.; Katayama, T.; Yamamoto, K.; Kumagai, H.; Ashida, H.; Hirose, J.; Kitaoka, M. Physiology of consumption of human milk oligosaccharides by infant gut-associated bifidobacteria. J. Biol. Chem. 2011, 286, 34583–34592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Mudd, A.T.; Fleming, S.A.; Labhart, B.; Chichlowski, M.; Berg, B.M.; Donovan, S.M.; Dilger, R.N. Dietary Sialyllactose Influences Sialic Acid Concentrations in the Prefrontal Cortex and Magnetic Resonance Imaging Measures in Corpus Callosum of Young Pigs. Nutrients 2017, 9, 1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Jacobi, S.K.; Yatsunenko, T.; Li, D.; Dasgupta, S.; Yu, R.K.; Berg, B.M.; Chichlowski, M.; Odle, J. Dietary Isomers of Sialyllactose Increase Ganglioside Sialic Acid Concentrations in the Corpus Callosum and Cerebellum and Modulate the Colonic Microbiota of Formula-Fed Piglets. J. Nutr. 2016, 146, 200–208. [Google Scholar] [CrossRef] [Green Version]
  14. Tarr, A.J.; Galley, J.D.; Fisher, S.E.; Chichlowski, M.; Berg, B.M.; Bailey, M.T. The prebiotics 3′Sialyllactose and 6′Sialyllactose diminish stressor-induced anxiety-like behavior and colonic microbiota alterations: Evidence for effects on the gut–brain axis. Brain. Behav. Immun. 2015, 50, 166–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Oliveros, E.; Ramirez, M.; Vazquez, E.; Barranco, A.; Gruart, A.; Delgado-Garcia, J.M.; Buck, R.; Rueda, R.; Martin, M.J. Oral supplementation of 2′-fucosyllactose during lactation improves memory and learning in rats. J. Nutr. Biochem. 2016, 31, 20–27. [Google Scholar] [CrossRef]
  16. Vázquez, E.; Barranco, A.; Ramírez, M.; Gruart, A.; Delgado-García, J.M.; Martínez-Lara, E.; Blanco, S.; Martín, M.J.; Castanys, E.; Buck, R.; et al. Effects of a human milk oligosaccharide, 2′-fucosyllactose, on hippocampal long-term potentiation and learning capabilities in rodents. J. Nutr. Biochem. 2015, 26, 455–465. [Google Scholar] [CrossRef] [PubMed]
  17. Vazquez, E.; Barranco, A.; Ramirez, M.; Gruart, A.; Delgado-Garcia, J.M.; Jimenez, M.L.; Buck, R.; Rueda, R. Dietary 2′-Fucosyllactose Enhances Operant Conditioning and Long-Term Potentiation via Gut-Brain Communication through the Vagus Nerve in Rodents. PLoS ONE 2016, 11, e0166070. [Google Scholar] [CrossRef]
  18. Jia, S.; Lu, Z.; Gao, Z.; An, J.; Wu, X.; Li, X.; Dai, X.; Zheng, Q.; Sun, Y. Chitosan oligosaccharides alleviate cognitive deficits in an amyloid-β1–42-induced rat model of Alzheimer’s disease. Int. J. Biol. Macromol. 2016, 83, 416–425. [Google Scholar] [CrossRef]
  19. Berger, P.K.; Plows, J.F.; Jones, R.B.; Alderete, T.L.; Yonemitsu, C.; Poulsen, M.; Ryoo, J.H.; Peterson, B.S.; Bode, L.; Goran, M.I. Human milk oligosaccharide 2′-fucosyllactose links feedings at 1 month to cognitive development at 24 months in infants of normal and overweight mothers. PLoS ONE 2020, 15, e0228323. [Google Scholar] [CrossRef]
  20. Clarke, G.; Grenham, S.; Scully, P.; Fitzgerald, P.; Moloney, R.; Shanahan, F.; Dinan, T.; Cryan, J. The microbiome-gut-brain axis during early life regulates the hippocampal serotonergic system in a sex-dependent manner. Mol. Psychiatry 2013, 18, 666–673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Zheng, P.; Zeng, B.; Liu, M.; Chen, J.; Pan, J.; Han, Y.; Liu, Y.; Cheng, K.; Zhou, C.; Wang, H.; et al. The gut microbiome from patients with schizophrenia modulates the glutamate-glutamine-GABA cycle and schizophrenia-relevant behaviors in mice. Sci. Adv. 2019, 5, 1–11. [Google Scholar] [CrossRef] [Green Version]
  22. Gacias, M.; Gaspari, S.; Santos, P.-M.G.; Tamburini, S.; Andrade, M.; Zhang, F.; Shen, N.; Tolstikov, V.; Kiebish, M.A.; Dupree, J.L.; et al. Microbiota-driven transcriptional changes in prefrontal cortex override genetic differences in social behavior. Elife 2016, 5. [Google Scholar] [CrossRef] [PubMed]
  23. Hanstock, T.L.; Clayton, E.H.; Li, K.M.; Mallet, P.E. Anxiety and aggression associated with the fermentation of carbohydrates in the hindgut of rats. Physiol. Behav. 2004, 82, 357–368. [Google Scholar] [CrossRef] [PubMed]
  24. Hanstock, T.L.; Mallet, P.E.; Clayton, E.H. Increased plasma d-lactic acid associated with impaired memory in rats. Physiol. Behav. 2010, 101, 653–659. [Google Scholar] [CrossRef] [PubMed]
  25. Stilling, R.M.; van de Wouw, M.; Clarke, G.; Stanton, C.; Dinan, T.G.; Cryan, J.F. The neuropharmacology of butyrate: The bread and butter of the microbiota-gut-brain axis? Neurochem. Int. 2016, 99, 110–132. [Google Scholar] [CrossRef] [PubMed]
  26. Bravo, J.A.; Forsythe, P.; Chew, M.V.; Escaravage, E.; Savignac, H.M.; Dinan, T.G.; Bienenstock, J.; Cryan, J.F. Ingestion of Lactobacillus strain regulates emotional behavior and central GABA receptor expression in a mouse via the vagus nerve. Proc. Natl. Acad. Sci. USA 2011, 108, 16050–16055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Carlson, A.L.; Xia, K.; Azcarate-Peril, M.A.; Goldman, B.D.; Ahn, M.; Styner, M.A.; Thompson, A.L.; Geng, X.; Gilmore, J.H.; Knickmeyer, R.C. Infant Gut Microbiome Associated With Cognitive Development. Biol. Psychiatry 2018, 83, 148–159. [Google Scholar] [CrossRef] [PubMed]
  28. Valles-Colomer, M.; Falony, G.; Darzi, Y.; Tigchelaar, E.F.; Wang, J.; Tito, R.Y.; Schiweck, C.; Kurilshikov, A.; Joossens, M.; Wijmenga, C.; et al. The neuroactive potential of the human gut microbiota in quality of life and depression. Nat. Microbiol. 2019, 4, 623–632. [Google Scholar] [CrossRef]
  29. Golubeva, A.V.; Joyce, S.A.; Moloney, G.; Burokas, A.; Sherwin, E.; Arboleya, S.; Flynn, I.; Khochanskiy, D.; Moya-Pérez, A.; Peterson, V.; et al. Microbiota-related Changes in Bile Acid & Tryptophan Metabolism are Associated with Gastrointestinal Dysfunction in a Mouse Model of Autism. EBioMedicine 2017, 24, 166–178. [Google Scholar] [CrossRef] [Green Version]
  30. Luna, R.A.; Oezguen, N.; Balderas, M.; Venkatachalam, A.; Runge, J.K.; Versalovic, J.; Veenstra-VanderWeele, J.; Anderson, G.M.; Savidge, T.; Williams, K.C. Distinct Microbiome-Neuroimmune Signatures Correlate With Functional Abdominal Pain in Children With Autism Spectrum Disorder. CMGH 2017, 3, 218–230. [Google Scholar] [CrossRef] [Green Version]
  31. Borre, Y.E.; O’Keeffe, G.W.; Clarke, G.; Stanton, C.; Dinan, T.G.; Cryan, J.F. Microbiota and neurodevelopmental windows: Implications for brain disorders. Trends Mol. Med. 2014, 20, 509–518. [Google Scholar] [CrossRef]
  32. Fleming, S.A.; Mudd, A.T.; Hauser, J.; Yan, J.; Metairon, S.; Steiner, P.; Donovan, S.M.; Dilger, R.N. Dietary Oligofructose Alone or in Combination with 2′-Fucosyllactose Differentially Improves Recognition Memory and Hippocampal mRNA Expression. Nutrients 2020, 12, 2131. [Google Scholar] [CrossRef] [PubMed]
  33. Fleming, S.A.; Mudd, A.T.; Hauser, J.; Yan, J.; Metairon, S.; Steiner, P.; Donovan, S.M.; Dilger, R.N. Human and Bovine Milk Oligosaccharides Elicit Improved Recognition Memory Concurrent with Alterations in Regional Brain Volumes and Hippocampal mRNA Expression. Front. Neurosci. 2020, 14, 1–14. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, M.; Monaco, M.H.; Hauser, J.; Yan, J.; Dilger, R.N.; Donovan, S.M. Bovine Milk Oligosaccharides and Human Milk Oligosaccharides Modulate the Gut Microbiota Composition and Volatile Fatty Acid Concentrations in a Neonatal Model. Microorganisms 2021, 2, 331–351. [Google Scholar]
  35. Niilo, L. Clostridium perfringens Type C Enterotoxemia. Can. Vet. J. 1988, 29, 658–664. [Google Scholar] [PubMed]
  36. Meli, F.; Puccio, G.; Cajozzo, C.; Ricottone, G.L.; Pecquet, S.; Sprenger, N.; Steenhout, P. Growth and safety evaluation of infant formulae containing oligosaccharides derived from bovine milk: A randomized, double-blind, noninferiority trial. BMC Pediatr. 2014, 14, 306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cooper, P.; Bolton, K.D.; Velaphi, S.; De Groot, N.; Emady-Azar, S.; Pecquet, S.; Steenhout, P. Early Benefits of a Starter Formula Enriched in Prebiotics and Probiotics on the Gut Microbiota of Healthy Infants Born to HIV+ Mothers: A Randomized Double-Blind Controlled Trial. Clin. Med. Insights Pediatr. 2016, 10, 119–130. [Google Scholar] [CrossRef] [Green Version]
  38. Radke, M.; Picaud, J.C.; Loui, A.; Cambonie, G.; Faas, D.; Lafeber, H.N.; De Groot, N.; Pecquet, S.S.; Steenhout, P.G.; Hascoet, J.M. Starter formula enriched in prebiotics and probiotics ensures normal growth of infants and promotes gut health: A randomized clinical trial. Pediatr. Res. 2017, 81, 622–631. [Google Scholar] [CrossRef] [Green Version]
  39. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Marriage, B.J.; Buck, R.H.; Goehring, K.C.; Oliver, J.S.; Williams, J.A. Infants Fed a Lower Calorie Formula with 2′FL Show Growth and 2′FL Uptake Like Breast-Fed Infants. J. Pediatr. Gastroenterol. Nutr. 2015, 61, 649–658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Yao, M.; Lien, E.L.; Capeding, M.R.Z.; Fitzgerald, M.; Ramanujam, K.; Yuhas, R.; Northington, R.; Lebumfacil, J.; Wang, L.; DeRusso, P.A. Effects of Term Infant Formulas Containing High sn-2 Palmitate With and Without Oligofructose on Stool Composition, Stool Characteristics, and Bifidogenicity. J. Pediatr. Gastroenterol. Nutr. 2014, 59, 440–448. [Google Scholar] [CrossRef] [Green Version]
  42. Wernimont, S.; Northington, R.; Kullen, M.J.; Yao, M.; Bettler, J. Effect of an α-lactalbumin-enriched infant formula supplemented with oligofructose on fecal microbiota, stool characteristics, and hydration status: A randomized, double-blind, controlled trial. Clin. Pediatr. 2015, 54, 359–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Moustgaard, A.; Lind, N.M.; Hemmingsen, R.; Hansen, A.K. Spontaneous Object Recognition in the Göttingen Minipig. Neural Plast. 2002, 9, 255–259. [Google Scholar] [CrossRef]
  44. Kornum, B.R.; Thygesen, K.S.; Nielsen, T.R.; Knudsen, G.M.; Lind, N.M. The effect of the inter-phase delay interval in the spontaneous object recognition test for pigs. Behav. Brain Res. 2007, 181, 210–217. [Google Scholar] [CrossRef] [PubMed]
  45. Gifford, A.K.; Cloutier, S.; Newberry, R.C. Objects as enrichment: Effects of object exposure time and delay interval on object recognition memory of the domestic pig. Appl. Anim. Behav. Sci. 2007, 107, 206–217. [Google Scholar] [CrossRef]
  46. Kouwenberg, A.-L.; Walsh, C.J.; Morgan, B.E.; Martin, G.M. Episodic-like memory in crossbred Yucatan minipigs (Sus scrofa). Appl. Anim. Behav. Sci. 2009, 117, 165–172. [Google Scholar] [CrossRef]
  47. Fleming, S.A.; Dilger, R.N. Young pigs exhibit differential exploratory behavior during novelty preference tasks in response to age, sex, and delay. Behav. Brain Res. 2017, 321, 50–60. [Google Scholar] [CrossRef]
  48. Fleming, S.A.; Monaikul, S.; Patsavas, A.J.; Waworuntu, R.V.; Berg, B.M.; Dilger, R.N. Dietary polydextrose and galactooligosaccharide increase exploratory behavior, improve recognition memory, and alter neurochemistry in the young pig. Nutr. Neurosci. 2017, 22, 499–512. [Google Scholar] [CrossRef]
  49. Fleming, S.; Chichlowski, M.; Berg, B.; Donovan, S.; Dilger, R. Dietary Sialyllactose Does Not Influence Measures of Recognition Memory or Diurnal Activity in the Young Pig. Nutrients 2018, 10, 395. [Google Scholar] [CrossRef] [Green Version]
  50. Jacob, R.M.; Mudd, A.T.; Alexander, L.S.; Lai, C.-S.; Dilger, R.N. Comparison of Brain Development in Sow-Reared and Artificially Reared Piglets. Front. Pediatr. 2016, 4, 95. [Google Scholar] [CrossRef] [Green Version]
  51. Radlowski, E.C.; Conrad, M.S.; Lezmi, S.; Dilger, R.N.; Sutton, B.; Larsen, R.; Johnson, R.W. A neonatal piglet model for investigating brain and cognitive development in small for gestational age human infants. PLoS ONE 2014, 9, e91951. [Google Scholar] [CrossRef] [Green Version]
  52. Mudd, A.T.; Getty, C.M.; Sutton, B.P.; Dilger, R.N. Perinatal choline deficiency delays brain development and alters metabolite concentrations in the young pig. Nutr. Neurosci. 2016, 19, 425–433. [Google Scholar] [CrossRef]
  53. Reznikov, E.A.; Comstock, S.S.; Hoeflinger, J.L.; Wang, M.; Miller, M.J.; Donovan, S.M. Dietary Bovine Lactoferrin Reduces Staphylococcus aureus in the Tissues and Modulates the Immune Response in Piglets Systemically Infected with S. aureus. Curr. Dev. Nutr. 2018, 2, nzy001. [Google Scholar] [CrossRef]
  54. Monaco, M.H.; Wang, M.; Pan, X.; Li, Q.; Richards, J.D.; Chichlowski, M.; Berg, B.M.; Dilger, R.N.; Donovan, S.M. Evaluation of Sialyllactose Supplementation of a Prebiotic-Containing Formula on Growth, Intestinal Development, and Bacterial Colonization in the Neonatal Piglet. Curr. Dev. Nutr. 2018, 2, nzy067. [Google Scholar] [CrossRef]
  55. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  56. Cole, J.R.; Wang, Q.; Cardenas, E.; Fish, J.; Chai, B.; Farris, R.J.; Kulam-Syed-Mohideen, A.S.; McGarrell, D.M.; Marsh, T.; Garrity, G.M.; et al. The Ribosomal Database Project: Improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37, D141–D145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. R Core Team R. A Language and Environment for Statistical Computing; Vienna, Austria. 2019. Available online: https://www.R-project.org/ (accessed on 1 January 2021).
  58. Tingley, D.; Yamamoto, T.; Hirose, K.; Keele, L.; Imai, K. mediation: R package for Causal Mediation Analysis. J. Stat. Softw. 2014, 59, 1–38. [Google Scholar] [CrossRef] [Green Version]
  59. Gu, Z.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [Green Version]
  60. Hauke, J.; Kossowski, T.M. Comparison of values of pearson’s and spearman’s correlation coefficients on the same sets of data. Quaest. Geogr. 2011, 30, 87–93. [Google Scholar] [CrossRef] [Green Version]
  61. Hayes, A.F.; Rockwood, N.J. Regression-based statistical mediation and moderation analysis in clinical research: Observations, recommendations, and implementation. Behav. Res. Ther. 2017, 98, 39–57. [Google Scholar] [CrossRef]
  62. Hayes, A.F. Introduction to Mediation, Moderation, and Conditional Process Analysis, 2nd ed.; Little, T., Ed.; The Guilford Press: New York, NY, USA, 2018. [Google Scholar]
  63. Hines, R.M.; Hines, D.J.; Houston, C.M.; Mukherjee, J.; Haydon, P.G.; Tretter, V.; Smart, T.G.; Moss, S.J. Disrupting the clustering of GABAA receptor α2 subunits in the frontal cortex leads to reduced γ-power and cognitive deficits. Proc. Natl. Acad. Sci. USA 2013, 110, 16628–16633. [Google Scholar] [CrossRef] [Green Version]
  64. Ge, Y.; Dong, Z.; Bagot, R.C.; Howland, J.G.; Phillips, A.G.; Wong, T.P.; Wang, Y.T. Hippocampal long-term depression is required for the consolidation of spatial memory. Proc. Natl. Acad. Sci. USA 2010, 107, 16697–16702. [Google Scholar] [CrossRef] [Green Version]
  65. Brigman, J.L.; Wright, T.; Talani, G.; Prasad-Mulcare, S.; Jinde, S.; Seabold, G.K.; Mathur, P.; Davis, M.I.; Bock, R.; Gustin, R.M.; et al. Loss of GluN2B-Containing NMDA Receptors in CA1 Hippocampus and Cortex Impairs Long-Term Depression, Reduces Dendritic Spine Density, and Disrupts Learning. J. Neurosci. 2010, 30, 4590–4600. [Google Scholar] [CrossRef]
  66. Zhao, X.; Rosenke, R.; Kronemann, D.; Brim, B.; Das, S.R.; Dunah, A.W.; Magnusson, K.R. The effects of aging on N-methyl-D-aspartate receptor subunits in the synaptic membrane and relationships to long-term spatial memory. Neuroscience 2009, 162, 933–945. [Google Scholar] [CrossRef] [Green Version]
  67. Li, H.; Santos, M.S.; Park, C.K.; Dobry, Y.; Voglmaier, S.M. VGLUT2 Trafficking Is Differentially Regulated by Adaptor Proteins AP-1 and AP-3. Front. Cell. Neurosci. 2017, 11, 324. [Google Scholar] [CrossRef]
  68. Birgner, C.; Nordenankar, K.; Lundblad, M.; Mendez, J.A.; Smith, C.; le Grevès, M.; Galter, D.; Olson, L.; Fredriksson, A.; Trudeau, L.-E.; et al. VGLUT2 in dopamine neurons is required for psychostimulant-induced behavioral activation. Proc. Natl. Acad. Sci. USA 2010, 107, 389–394. [Google Scholar] [CrossRef] [Green Version]
  69. Ménard, C.; Quirion, R.; Vigneault, E.; Bouchard, S.; Ferland, G.; El Mestikawy, S.; Gaudreau, P. Glutamate presynaptic vesicular transporter and postsynaptic receptor levels correlate with spatial memory status in aging rat models. Neurobiol. Aging 2015, 36, 1471–1482. [Google Scholar] [CrossRef] [Green Version]
  70. Amilhon, B.; Lepicard, E.; Renoir, T.; Mongeau, R.; Popa, D.; Poirel, O.; Miot, S.; Gras, C.; Gardier, A.M.; Gallego, J.; et al. VGLUT3 (Vesicular Glutamate Transporter Type 3) Contribution to the Regulation of Serotonergic Transmission and Anxiety. J. Neurosci. 2010, 30, 2198–2210. [Google Scholar] [CrossRef]
  71. Fazekas, C.L.; Balázsfi, D.; Horváth, H.R.; Balogh, Z.; Aliczki, M.; Puhova, A.; Balagova, L.; Chmelova, M.; Jezova, D.; Haller, J.; et al. Consequences of VGluT3 deficiency on learning and memory in mice. Physiol. Behav. 2019, 112688. [Google Scholar] [CrossRef]
  72. Tabouy, L.; Getselter, D.; Ziv, O.; Karpuj, M.; Tabouy, T.; Lukic, I.; Maayouf, R.; Werbner, N.; Ben-Amram, H.; Nuriel-Ohayon, M.; et al. Dysbiosis of microbiome and probiotic treatment in a genetic model of autism spectrum disorders. Brain. Behav. Immun. 2018, 73, 310–319. [Google Scholar] [CrossRef]
  73. Hoban, A.E.; Stilling, R.M.; Ryan, F.J.; Shanahan, F.; Dinan, T.G.; Claesson, M.J.; Clarke, G.; Cryan, J.F. Regulation of prefrontal cortex myelination by the microbiota. Transl. Psychiatry 2016, 6, e774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Lu, J.; Lu, L.; Yu, Y.; Cluette-Brown, J.; Martin, C.R.; Claud, E.C. Effects of Intestinal Microbiota on Brain Development in Humanized Gnotobiotic Mice. Sci. Rep. 2018, 8, 5443. [Google Scholar] [CrossRef] [Green Version]
  75. Heijtz, R.D.; Wang, S.; Anuar, F.; Qian, Y.; Björkholm, B.; Samuelsson, A.; Hibberd, M.L.; Forssberg, H.; Pettersson, S. Normal gut microbiota modulates brain development and behavior. Proc. Natl. Acad. Sci. USA 2011, 108, 3047–3052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Borrelli, L.; Aceto, S.; Agnisola, C.; De Paolo, S.; Dipineto, L.; Stilling, R.M.; Dinan, T.G.; Cryan, J.F.; Menna, L.F.; Fioretti, A. Probiotic modulation of the microbiota-gut-brain axis and behaviour in zebrafish. Sci. Rep. 2016, 6, 1–9. [Google Scholar] [CrossRef] [Green Version]
  77. Lu, J.; Synowiec, S.; Lu, L.; Yu, Y.; Bretherick, T.; Takada, S.; Yarnykh, V.; Caplan, J.; Caplan, M.; Claud, E.C.; et al. Microbiota influence the development of the brain and behaviors in C57BL/6J mice. PLoS ONE 2018, 13, e0201829. [Google Scholar] [CrossRef]
  78. Sampson, T.R.; Mazmanian, S.K. Control of brain development, function, and behavior by the microbiome. Cell Host Microbe 2015, 17, 565–576. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 09 00846 g001 550
Figure 1. Schematic depicting the mediation models and regression paths. The objective of the mediation analysis was to determine if the indirect effect (paths a*b, which is equal to paths c-c’) was different from zero, suggesting the mediating variable altered the strength of the relationship between X and Y. Three models were used to assess the mediation: X regressed on Y, X regressed on M, and both X and M regressed on Y. The letters a, b, c, and c’ refer to estimates the beta coefficients for each respective model, with e representing the error term.
Figure 1. Schematic depicting the mediation models and regression paths. The objective of the mediation analysis was to determine if the indirect effect (paths a*b, which is equal to paths c-c’) was different from zero, suggesting the mediating variable altered the strength of the relationship between X and Y. Three models were used to assess the mediation: X regressed on Y, X regressed on M, and both X and M regressed on Y. The letters a, b, c, and c’ refer to estimates the beta coefficients for each respective model, with e representing the error term.
Microorganisms 09 00846 g001
Microorganisms 09 00846 g002 550
Figure 2. (A) Relative abundance of taxa at the phylum level in the ascending colon and feces. Each bar represents an individual pig. (B) Relative abundance of common and uncommon taxa at the genera level. Groups with a median abundance less than 0.05% for any given taxa are not shown. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Figure 2. (A) Relative abundance of taxa at the phylum level in the ascending colon and feces. Each bar represents an individual pig. (B) Relative abundance of common and uncommon taxa at the genera level. Groups with a median abundance less than 0.05% for any given taxa are not shown. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Microorganisms 09 00846 g002
Microorganisms 09 00846 g003 550
Figure 3. Heatmaps demonstrating the relationships between predictor (X), mediating (M), and predicted variables (Y, short-term recognition memory). Only mediations whose confidence intervals of the indirect effect do not include zero are shown. The right-most box represents the estimates of the indirect effect (c-c’) bounded by its lower and upper 95% confidence intervals. Paths a, bym, b, c, and c’ are the regression coefficients for X → M, M → Y, X → Y, and X + M → Y, respectively. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Figure 3. Heatmaps demonstrating the relationships between predictor (X), mediating (M), and predicted variables (Y, short-term recognition memory). Only mediations whose confidence intervals of the indirect effect do not include zero are shown. The right-most box represents the estimates of the indirect effect (c-c’) bounded by its lower and upper 95% confidence intervals. Paths a, bym, b, c, and c’ are the regression coefficients for X → M, M → Y, X → Y, and X + M → Y, respectively. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Microorganisms 09 00846 g003
Microorganisms 09 00846 g004 550
Figure 4. Heatmaps demonstrating the relationships between predictor (X), mediating (M), and predicted variables (Y, long-term recognition memory). Only mediations whose confidence intervals of the indirect effect do not include zero are shown. The right-most box represents the estimates of the indirect effect (c-c’) bounded by its lower and upper 95% confidence intervals. Paths a, bym, b, c, and c’ are the regression coefficients for X → M, M → Y, X → Y, and X + M → Y, respectively. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Figure 4. Heatmaps demonstrating the relationships between predictor (X), mediating (M), and predicted variables (Y, long-term recognition memory). Only mediations whose confidence intervals of the indirect effect do not include zero are shown. The right-most box represents the estimates of the indirect effect (c-c’) bounded by its lower and upper 95% confidence intervals. Paths a, bym, b, c, and c’ are the regression coefficients for X → M, M → Y, X → Y, and X + M → Y, respectively. Abbreviations: Uncl., unclassified; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS, pigs fed bovine milk oligosaccharides; BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; OF, pigs feds oligofructose; OF + 2′-FL, pigs fed oligofructose and 2′-fucosyllactose.
Microorganisms 09 00846 g004
Table
Table 1. Carbohydrate content of the diets 1.
Table 1. Carbohydrate content of the diets 1.
Bovine Milk
Oligosaccharides 2		2′Fucosyllactose 3Lacto-N-Neotetraose 3Oligofructose 4Additional LactoseTotal OSTotal LactoseTotal
Carbohydrate
DietPer kgPer LPer kgPer LPer kgPer LPer kgPer LPer kgPer LPer kgPer LPer kgPer LPer kgPer L
Formulated
CON0000000069.3013.8600412.7982.56412.7982.56
BMOS61.9812.4000000061.981.4761.9812.40350.8270.16412.7982.56
HMO004.920.982.410.48007.3312.407.331.47405.4781.09412.7982.56
BMOS + HMO61.9812.404.920.982.410.48000069.3013.86343.4968.70412.7982.56
OF00000026.085.2261.988.6426.085.22386.7177.34412.7982.56
OF + 2′-FL004.920.980026.085.2238.307.6631.006.20381.7976.36412.7982.56
Analyzed
CON00000000NQNQ00NQNQNQNQ
BMOS28.955.79000000NQNQ28.955.79NQNQNQNQ
HMO004.050.812.100.4200NQNQ6.151.23NQNQNQNQ
BMOS + HMO28.755.755.001.002.650.5300NQNQ36.407.28NQNQNQNQ
OF00000018.103.62NQNQ18.103.62NQNQNQNQ
OF + 2′-FL005.601.120017.053.41NQNQ22.654.53NQNQNQNQ
1 Abbreviations: OS, oligosaccharide; CON, control group; HMO, pigs fed human milk oligosaccharides; BMOS; pigs fed bovine milk oligosaccharides, BMOS + HMO, pigs fed both human and bovine milk oligosaccharides; NQ, not quantified. 2 Nestlé Product & Technology center, Konolfingen, Switzerland. 3 Glycom, Hørsholm, Denmark. 4 Orafti® P95; Beneo-Orafti, Tienen, Belgium.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fleming, S.A.; Hauser, J.; Yan, J.; Donovan, S.M.; Wang, M.; Dilger, R.N. A Mediation Analysis to Identify Links between Gut Bacteria and Memory in Context of Human Milk Oligosaccharides. Microorganisms 2021, 9, 846. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040846

AMA Style

Fleming SA, Hauser J, Yan J, Donovan SM, Wang M, Dilger RN. A Mediation Analysis to Identify Links between Gut Bacteria and Memory in Context of Human Milk Oligosaccharides. Microorganisms. 2021; 9(4):846. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040846

Chicago/Turabian Style

Fleming, Stephen A., Jonas Hauser, Jian Yan, Sharon M. Donovan, Mei Wang, and Ryan N. Dilger. 2021. "A Mediation Analysis to Identify Links between Gut Bacteria and Memory in Context of Human Milk Oligosaccharides" Microorganisms 9, no. 4: 846. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms9040846

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop