Next Article in Journal
Mediterranean Diet Adherence, Physical Activity, and Advanced Glycation End Products in Complex PTSD: A Comprehensive Examination of Lifestyle and Cardiovascular Risk in War Veterans
Previous Article in Journal
What Can the Gut Microbiota of Animals Teach Us about the Relationship between Nutrition and Burden of Lifestyle Diseases?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Prebiotic Diet Containing Galactooligosaccharides and Polydextrose Produces Dynamic and Reproducible Changes in the Gut Microbial Ecosystem in Male Rats

by
Robert S. Thompson
1,2,
Samuel J. Bowers
3,
Fernando Vargas
4,
Shelby Hopkins
1,2,
Tel Kelley
1,
Antonio Gonzalez
5,
Christopher A. Lowry
1,2,
Pieter C. Dorrestein
4,
Martha Hotz Vitaterna
3,
Fred W. Turek
3,
Rob Knight
5,6,7,
Kenneth P. Wright, Jr.
1,2 and
Monika Fleshner
1,2,*
1
Department of Integrative Physiology, University of Colorado Boulder, Boulder, CO 80309, USA
2
Center for Neuroscience, University of Colorado Boulder, Boulder, CO 80309, USA
3
Department of Neurobiology, Northwestern University, Center for Sleep and Circadian Biology, Evanston, IL 60208, USA
4
Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California San Diego, San Diego, CA 92093, USA
5
Department of Pediatrics, University of California San Diego, San Diego, CA 92093, USA
6
Department of Computer Science and Engineering, University of California San Diego, San Diego, CA 92093, USA
7
Center for Microbiome Innovation, University of California San Diego, San Diego, CA 92093, USA
*
Author to whom correspondence should be addressed.
Nutrients 2024, 16(11), 1790; https://doi.org/10.3390/nu16111790
Submission received: 7 May 2024 / Revised: 30 May 2024 / Accepted: 4 June 2024 / Published: 6 June 2024

Abstract

:
Despite substantial evidence supporting the efficacy of prebiotics for promoting host health and stress resilience, few experiments present evidence documenting the dynamic changes in microbial ecology and fecal microbially modified metabolites over time. Furthermore, the literature reports a lack of reproducible effects of prebiotics on specific bacteria and bacterial-modified metabolites. The current experiments examined whether consumption of diets enriched in prebiotics (galactooligosaccharides (GOS) and polydextrose (PDX)), compared to a control diet, would consistently impact the gut microbiome and microbially modified bile acids over time and between two research sites. Male Sprague Dawley rats were fed control or prebiotic diets for several weeks, and their gut microbiomes and metabolomes were examined using 16S rRNA gene sequencing and untargeted LC–MS/MS analysis. Dietary prebiotics altered the beta diversity, relative abundance of bacterial genera, and microbially modified bile acids over time. PICRUSt2 analyses identified four inferred functional metabolic pathways modified by the prebiotic diet. Correlational network analyses between inferred metabolic pathways and microbially modified bile acids revealed deoxycholic acid as a potential network hub. All these reported effects were consistent between the two research sites, supporting the conclusion that dietary prebiotics robustly changed the gut microbial ecosystem. Consistent with our previous work demonstrating that GOS/PDX reduces the negative impacts of stressor exposure, we propose that ingesting a diet enriched in prebiotics facilitates the development of a health-promoting gut microbial ecosystem.

Graphical Abstract

1. Introduction

The gut microbiome is a diverse ecosystem that consists of bacteria, archaea, eukaryotes, fungi, and viruses that live in the host’s digestive tract [1,2,3]. Microorganisms residing in the digestive tract comprise a micro-ecosystem displaying established principles of ecosystem dynamics [4]. Like any ecosystem, gut microbes both compete and cooperate for limited resources [5]. Dietary macronutrients [6,7] and micronutrients [8] can rapidly change the gut’s microbial composition. Non-digestible complex carbohydrates and types of fermentable fiber, for example, are dietary substrates selectively utilized by host microorganisms that can rapidly alter the gut microbiome and the fecal metabolome and positively impact host health. In 2017, the International Scientific Association for Dietary Probiotics and Prebiotics released a consensus report defining prebiotics as substrates selectively utilized by host microorganisms conferring health benefits [9]. The complex changing dynamics in the gut’s microbial composition after introducing dietary prebiotics have seldom been examined, because doing so requires repeated sampling over time and costly sequencing of large numbers of samples.
Despite substantial evidence supporting the efficacy of prebiotics for promoting host health, there is a paucity of literature replicating prebiotic impacts on bacteria and bacterially modified metabolites over time [10,11,12,13]. The failure to reproduce findings could be due, in part, to multiple bacterial taxonomy databases, ongoing taxonomic revisions, as well as differences in sample storage, DNA extraction and sequencing, and analytic pipelines [14]. In addition, commonly overlooked are the influences of environmental factors on the gut microbiome, including geographic location and elevation [15,16,17,18], and animal source [19,20].
GOS and PDX increase the relative abundance of the bacterial species Parabacteroides distasonis and Clostridium leptum [21,22], decrease microbially modified secondary bile acids like deoxycholic and lithocholic acid [23,24], and reduce the adverse effects of stress exposure on host sleep physiology [25,26,27]. To make progress towards elucidating the mechanisms for the stress-protective impact on host sleep physiology [27,28], the effects of prebiotics on the gut microbiota and metabolome must be sufficiently robust to resist any potential environmental and methodological influences.
Here, we present the results from two dietary prebiotic animal studies conducted at Northwestern University (NW) in Evanston, Illinois, and the University of Colorado Boulder (CU) in Boulder, Colorado. The two sites have several environmental differences, including different research personnel, vivarium facilities, elevations (182 m vs. 1624 m), and animal sources (Envigo vs. Harlan). To reduce the impact of other factors, NW and CU adhered to standardized fecal sample collection and storage protocols, DNA extraction and sequencing, as well as untargeted LC–MS/MS metabolomics protocols.
The first goal of this project was to determine whether consumption of the same dietary prebiotic formulation tested at different universities, in different locations across the country, and at different times of the year, would produce similar dynamic changes in the gut microbial composition and microbially modified bile acids. The second goal of the study was to explore the potential functional metabolic pathways and networks impacted by the prebiotic diet. We hypothesize that the consumption of GOS/PDX by rats at NW and CU produces robust changes over time in the gut microbiome, fecal metabolome, functional metabolic pathways, and networks.

2. Materials and Methods

2.1. Animals

Male Sprague Dawley rats were tested. Female rats were not tested in these experiments because this study was supported by funding from the Office of Naval Research (ONR MURI N00014-15-1-2809), and ~80–90% of submariners are male, making males a priority for the ONR’s limited funding. Some data presented here were included in previously published work from the more extensive ONR study, which demonstrated that diets enriched in prebiotics (GOS and PDX) facilitate host sleep/circadian recovery both during and after stressor exposure [21,22].

2.1.1. Northwestern (NW) Study

The animals (N = 64, Envigo Laboratories, Madison, WI, USA) were singly housed in a controlled temperature (23 ± 2 °C) and humidity. All the protocols were approved by the Northwestern Institutional Animal Care and Use Committee, as previously described [21]. The animals weighed 40–50 g upon arrival at postnatal day (PND) 23 and were maintained in a 12:12 h light/dark cycle. On arrival, all the rats were housed in Nalgene Plexiglas cages (Thermo Fisher Scientific, Waltham, MA, USA) and were placed on a control or prebiotic diet (ad libitum).

2.1.2. University of Colorado Boulder (CU) Study

The animals (N = 82, Harlan Laboratories, Indianapolis, IN, USA) were singly housed with a controlled temperature (23 ± 2 °C) and humidity. All procedures were approved by the University of Colorado Boulder Institutional Animal Care and Use Committee, as previously described [22]. Briefly, the animals weighed 40–50 g upon arrival at PND 23 and were maintained in a 12:12 h light/dark cycle. On arrival, all the rats were housed in Nalgene Plexiglas cages (Thermo Fisher Scientific, Waltham, MA, USA) and were placed on a control or prebiotic diet (ad libitum).

2.2. Experimental Design

The rats arrived at NW on PND 23 and were randomly placed on either the control or prebiotic diet for the duration of the study (Figure 1). Animal numbers for the NW microbiome data were control (n = 30) and prebiotic diet (n = 32), while the animal numbers for the NW metabolome data were control (n = 31) and prebiotic diet (n = 32). The rats arrived at CU on PND 23 and were immediately placed on either the control or prebiotic diet for the duration of the study (Figure 1). Animal numbers for the CU microbiome data were control (n = 37) and prebiotic diet (n = 37), while the animal numbers for the CU metabolome data were control (n = 40) and prebiotic diet (n = 42). Only samples present for all time points with viable data (i.e., useable fecal samples, high-quality sequencing, quality feature detection, etc.) were included in the final analysis.
At NW, fecal samples were collected on experimental days 0, 28, 42, and 51. At CU, fecal samples were collected on experimental days 2, 33, 75, and 94 (Figure 1). The days chosen for fecal collection differed between the sites due to the other goals of the larger ONR project. The repeated sample collection over time allows one to capture any changes in the microbiome and metabolome due to aging (i.e., adolescence to young adulthood). The two experiments were conducted on Sprague Dawley rats eating identical diets and, thus, give us unique insights into how the gut microbiome and gut metabolome change from adolescence to young adulthood between study sites in response to a prebiotic diet (Figure 1).

2.3. Diets

Rats at both facilities had ad libitum access to control or prebiotic diets immediately upon arrival on PND 23. The control and prebiotic diets fed to rats at NW and CU were the same formulation. The diets were initially formulated by Mead Johnson Nutrition (MJN, Evansville, IN, USA) based on AIN-93G specifications, were custom made by Envigo Teklad (TD.110889; now Inotiv, Lafayette, IN, USA), and were isocaloric, with similar carbohydrate, protein, fat, vitamin, and mineral levels, the details of which have been previously published [27,29]. The prebiotic diet contained the following prebiotic substrates, which were absent from the control diet: galactooligosaccharides (GOS, 24.14 g/kg (7.0 g active); FrieslandCampina, Zwolle, The Netherlands), and polydextrose (PDX, 7.69 g/kg (7.0 g active); Danisco, Terre Haute, IN, USA).

2.4. Fecal Sample Collection Procedures

Fecal samples were collected and prepared as previously described [30] and were collected after cage change. Sterile forceps (100% ethanol) were used to obtain each sample, which were then placed in 1.5 mL sterile screw cap tubes (USA Scientific, Ocala, FL, USA) and put in liquid nitrogen. The samples were then transferred and stored at −80 °C for analyses later. Weekly fecal samples were collected during the light cycle (~900–1100 h) shortly after cage changes. Investigators collected the rat fecal samples immediately after the rats defecated in the new bedding, i.e., within ~10–30 min. At each collection time point, duplicate samples of bedding, water, food, and blank tubes were also collected to control for potential environmental influences on the microbiome and metabolome data. For both study sites, the fecal samples were cut in half lengthwise to ensure each animal’s microbiome and metabolomics data were generated from the same fecal pellet [31].

2.5. The 16S rRNA Gene Sequencing

For both study sites, DNA was extracted from fecal samples and the V4 region of the 16S rRNA gene was amplified using the 515f/806r primer pair with the barcode on the forward read [32], and sequenced as previously described [33]. The samples were purified and precipitated to remove polymerase chain reaction (PCR) artifacts; the samples were sequenced in multiplex using an Illumina HiSeq 2000 (San Diego, CA, USA). All the target gene sequence processing was conducted with Quantitative Insights Into Microbial Ecology (QIIME2) [34] via Qiita. The raw sequencing data were trimmed and demultiplexed at 150 bases. Amplicon sequence variants (ASVs) were generated using the deblur algorithm. Phylogeny was created via SEPP within the QIIME2 fragment insertion plugin, using default parameters. Taxonomy classification was conducted via the QIIME2 feature classifier plugin and based on SILVA [35]. The resulting ASV table was filtered to remove mislabeled samples with a probability above 0.20 using the sample type field, as described in the Human Microbiome Project [36]. The resulting table was then rarefied at 10,000 sequences/sample to correct for an uneven sequencing depth due to amplification differences between the samples.
Beta diversity was examined with principal coordinate analysis (PCoA) using unweighted UniFrac distances (sensitive to rarer taxa) and weighted UniFrac distances (sensitive to abundances of taxa), which are the best ways to visualize the microbiome between treatments as a whole [37]. For analysis, PERMANOVA was used on each time point in QIIME2. Alpha diversity is a within-samples measure and was examined using evenness, observed OTUs, and Faith’s phylogenetic diversity [38]. Differential abundance was assessed in regard to the ASVs, using analysis of the composition of microbiomes (ANCOM) [39], as implemented in QIIME2 and matched with the SILVA database. Consistent with current recommended best practices [40], we refer to the taxonomy assignments as they are designated in the SILVA database since it is updated annually [35,41], and is based on ASVs, not the construction of molecular operational taxonomic units (OTUs) [40].
PICRUSt2 (https://github.com/picrust/picrust2, accessed on 30 June 2022) was performed in the conda environment for both studies, to identify functionally enriched signaling pathways due to prebiotic diet consumption [42].
The 16S rRNA gene sequencing data were uploaded to Qiita, are publicly available, and can be found at https://qiita.ucsd.edu/study/description/11697 (accessed on 18 November 2021) for the NW study and at https://qiita.ucsd.edu/study/description/11525 (accessed on 18 November 2021) for the CU study.

2.6. LC–MS/MS Metabolomics

Fecal and environmental samples were transferred overnight via dry ice to the University of California San Diego and processed for metabolomic analysis. The fecal samples were stored in 1.5 mL centrifuge tubes at −80 °C prior to extraction. Sample IDs were uploaded into an electronic spreadsheet and subsequently used to assign filenames during LC–MS/MS data acquisition. All solvents used for the metabolomic analysis were of LC–MS grade.
This method was adapted from a previously published protocol [43]. Fecal pellets were weighed at 50.0 ± 2 mg wet weight and transferred to 2.0 mL round bottom microcentrifuge tubes (Qiagen Catalog# 990381, Hilden, Germany) for metabolite extraction. A clean stainless-steel bead (Qiagen Catalog# 69989) and 1.5 mL of chilled extraction solvent (50% MeOH) were added to each sample. The samples were then homogenized for 5 min at 25 Hz using a TissueLyser II system (Qiagen Catalog# 85300) and incubated for 20 min at −20 °C. The fecal homogenates were centrifuged at 14,000 rpm for 15 min at 4 °C. Then, 1.2 mL aliquots were transferred into a Nunc 2.0 mL DeepWell plate (Thermo Catalog# 278743) and frozen at −80 °C, before lyophilization using a FreeZone 4.5 L Benchtop Freeze Dryer with Centrivap Concentrator (Labconco, Kansas City, MO, USA). The wells were resuspended with 200 µL of resuspension solvent (50% MeOH spiked with 2.0 µM sulfadimethoxine), vortexed for 30 s, and centrifuged at 2000 rpm for 15 min at 4 °C. Then, 150 µL of the supernatant was transferred into a 96-well plate and maintained at 4 °C, before LC–MS analysis. A resuspension solvent QC and a six standard mix QC (50% MeOH spiked with 1.0 µM sulfamethazine, 1.0 µM sulfamethizole, 1.0 µM sulfachloropyridazine, 1.0 µM amitriptyline, and 1.0 µM coumarin 314) was run every 12th sample to assess the sample background, carry over, chromatography behavior, peak picking, and plate effects.
The fecal extracts were analyzed using an ultra-high performance liquid chromatography system (Vanquish, Thermo Fisher Scientific, Waltham, MA, USA), coupled to a hybrid quadrupole-Orbitrap mass spectrometer (Q-Exactive, Thermo), fitted with a HESI probe. Reverse phase chromatographic separation was achieved using a Kinetex C18 1.7 µm, 100 Å, 50 × 2.1 mm column (Phenomenex, Torrance, CA, USA) held at 40 °C, with a 0.5 mL/min flow rate. Moreover, 5.0 µL aliquots were injected per sample/QC. The mobile phase used was: (A) 0.1% formic acid in water and (B) 0.1% formic acid in acetonitrile. The elution gradient was: 5% B for 1 min, increased to 100% B in the next 8 min, held at 100% B for 2 min, returned to 5.0% B in 0.5 min, and equilibrated at 5.0% B for 2 min. The positive electrospray ionization parameters were: a sheath gas flow rate of 52 (arb. units), an aux gas flow rate of 14 (arb. units), a sweep gas flow rate of 3 (arb. units), a spray voltage of 3.5 kV, a capillary temperature of 270 °C, an S-Lens RF level of 50 (arb. units), and an aux gas heater temperature of 435 °C. The negative electrospray ionization parameters were: a sheath gas flow rate of 52 (arb. units), an aux gas flow rate of 14 (arb. units), a sweep gas flow rate of 3 (arb. units), a spray voltage of 2.5 kV, a capillary temperature of 270 °C, an S-Lens RF level of 50 (arb. units), and an aux gas heater temperature of 435 °C. MS data were acquired using a data dependent acquisition method, with a resolution of 35,000 in MS1 and 17,000 in MS2. An MS1 scan from 100–1500 m/z was followed by an MS2 scan, produced by collision-induced disassociation, of the five most abundant ions from the prior MS1 scan.
Feature tables were generated for the control and prebiotic diet samples. To annotate features with a level 1 metabolome standard initiative (MSI) level of confidence, the mass and retention time were aligned and the MS/MS fragmentation pattern was compared between the features and 20 purified bile acid reference standards, as previously described in detail [30,44]. Primary, secondary, conjugated, and unconjugated bile acids were purchased (Cayman Chemical, Ann Arbor, MI, USA) and used to identify level 1 bile acid identification in fecal metabolomics samples. The samples were solubilized to a final concentration of 10 μM in 50% MeOH, before LC–MS/MS injection.
All untargeted mass spectrometry data can be found in the online mass spectrometry repository, Massive (http://massive.ucsd.edu, accessed on 11 May 2022), using the following accession numbers for NW, MSV000083073, and for CU, MSV000080628.

2.7. Statistical Analysis

The data were analyzed using R statistics version 4.2.2 GUI 1.79 Big Sur ARM build (8160). The data depicted in the figures were made in Prism (version 9.3.1). For the gut microbiome analysis of the UniFrac distance matrices, permutation multivariate analysis of variance (PERMANOVA) was used at each time point [45,46]. Measures of alpha diversity were analyzed separately using repeated measures ANOVA. To investigate differential abundance of genera level taxa between the control and prebiotic diets, a first-level analysis of the composition of the microbiome (ANCOM) was performed on the ASVs [39] to reveal reliable changes. ANCOM analysis will correct for multiple comparisons of ASVs identified in the sequencing data. The ASVs that were undefined/unclassified at the genera level were excluded from the final analysis. Once the taxonomy was assigned, we performed a second level of analysis on genus-level taxonomy assignments using the Nonparametric Tests for Repeated Measures Data in Factorial Designs (nparLD) package version 2.2. Importantly, only genera that were significantly changed by the prebiotic diet based on the ANCOM analysis are presented in this manuscript. Lower relative abundance genera were non-normally distributed; therefore, these data were analyzed using the nparLD package. The bile acid data were log transformed, as previously described [22,44], and analyzed using the nparLD package. Multiple significant p-values in the bile acid data were adjusted using the Holm method. The pathways output from PICRUSt2 was analyzed via DESeq2 version 1.14.1, using the Bioconductor R package, as previously described [21], and volcano plot analysis by the time point. The pathways affected by the prebiotic diet between the study sites and over time were analyzed using nparLD. Tukey’s post hoc analysis was used when appropriate using the nparcomp, the nonparametric relative contrast effects (nparcomp) package version 3.0, for relative abundance of the genera, bile acids, and pathway data. Network analyses examining the relationships between functionally significant pathways and bile acids were performed using the corrr package version 0.4.4. The two-tailed alpha level was set at p < 0.05.

3. Results

3.1. Microbiome

A prebiotic diet significantly changed the beta diversity of the gut microbiome at both study sites (Figure 2). Table 1 denotes the significant effects of a prebiotic diet on weighted and unweighted UniFrac distances. The prebiotic diet had no effect at 0 days on the diet on the weighted or unweighted UniFrac distance (NW). There was a significant effect of the prebiotic diet at 2 days on the diet on the weighted UniFrac distance (CU). The prebiotic diet significantly impacted both the weighted and unweighted UniFrac distance on all the remaining days on the diet at both study sites (Figure 2, Table 1).
The prebiotic diet altered the two main phyla (Firmicutes and Bacteroidetes) in the rat gut microbiome (Supplemental Figure S1). The main significant main effects of the prebiotic diet were on the Firmicutes at NW (F(1,2.77) = 9.72; p = 0.002; Supplemental Figure S1A) and CU (F(1,2.85) = 20.94; p = 0.0000078; Supplemental Figure S1A). The effect of the prebiotic diet changed over time at NW (time-by-diet interaction, F(1,59.99) = 4.61; p = 0.004), but not CU. The prebiotic diet also impacted the Bacteroidetes at both NW (F(1,2.73) = 6.01; p = 0.014; Supplemental Figure S1B) and CU (F(1,2.83) = 20.94; p = 0.0000047; Supplemental Figure S1B). Finally, there were time-by-diet interactions in regard to the Bacteroidetes that changed over time at both NW (F(1,59.76) = 3.33; p = 0.022; Supplemental Figure S1B) and CU (F(1,71.01) = 3.39; p = 0.019; Supplemental Figure S1B). While these phyla changes are important, it is of more interest to examine the taxonomic changes in greater detail at the genera level.
The top nine most abundant genera increased by the prebiotic diet, when compared to the control diet, are shown in Figure 3. The prebiotic diet increased the relative abundance of the genera Bacteroides (Figure 3A), Parabacteroides (Figure 3B), Clostridia_UCG_014 (Figure 3C), Incertae_Sedis (formerly known as Ruminiclostridium V) (Figure 3E), Parasutterella (Figure 3F), Ruminococcus_gauvreauii_group (Figure 3G), and UCG-007 (Figure 3H), at both study sites. See Table 2 and Figure 3 for complete statistical and post hoc analysis results. There were differing effects between study sites in regard to Christensenellaceae_R-7_group (Figure 3D) and Lachnospiraceae_UCG-006 (Figure 3I), such that the prebiotic diet increased these genera in the CU study, but not in the NW study. However, there was a diet-by-time interaction in regard to Lachnospiraceae_UCG-006 at NW (see Figure 3I for the results of the post hoc analysis).
The top six (out of nine) most abundant genera that were consistently lower in the prebiotic diet group, when compared to the control diet group, are shown in Figure 4. Lachnospiraceae_NK4A136_group (Figure 4A), Eubacterium_fissicatena_group (Figure 4E), Eubacterium_ruminantium_group (Figure 4F), GCA-900066575 (Figure 4G), Rosburia (Figure 4H), and Rikenellaceae_RC9_gut_group (Figure 4I) were consistently lower in the prebiotic diet group, when compared to the control diet groups, at each study site (see Table 2 for the statistical analysis; see Figure 4 for the results of the post hoc tests). Colidextribacter (Figure 4D) and UCG-005 (Figure 4C) were lower in the prebiotic diet group in the NW study, but not in the CU study (Table 2). The prebiotic diet had no significant effects on Eubacterium_coprostanoligenes_group at either study site (Figure 4B).
Supplemental Figure S2 depicts the additional ultra-low abundance genera (less than 1% relative abundance). Supplemental Table S1 contains the corresponding statistics. Overall, the impact of a prebiotic diet on ultra-low relative abundance (<1%) genera was less consistent between the study sites, except for the genus Tuzzerella. Tuzzerella was lower in the prebiotic diet groups at both study sites (Supplemental Figure S2, Supplemental Table S1). Finally, the prebiotic diet impacted six genera with opposite effects between study sites, although these were not statistically significant once corrected via Tukey’s post hoc analysis (Supplemental Figure S3, Supplemental Table S2). One interesting finding was the difference in the relative abundance in the genus Muribaculaceae between the study sites. The relative abundance of this genus was higher in the NW versus the CU microbiome data. This large inherent environmental difference present in the genus Muribaculaceae may have played a role in the contrasting effects of the prebiotic diet in the alpha diversity results between the study sites.
Dietary prebiotics inconsistently altered the alpha diversity at both study sites. In the NW study, there was no effect of the prebiotic diet on evenness (Figure 5A); however, in the CU study, the prebiotic diet increased the overall evenness (F(1,3) = 12.27; p = 0.00084; Figure 5A). In the NW study, there was a significant decrease in Faith’s phylogenetic diversity due to the prebiotic diet (F(1,3) = 5.82; p = 0.021), while no effects were found due to the prebiotic diet in Faith’s phylogenetic diversity in the CU study (Figure 5B). The observed features were lower due to the prebiotic diet at NW (F(1,3) = 6.25; p = 0.017), and there was a significant diet-by-time interaction (F(1,3) = 2.95; p = 0.035), but the prebiotic diet did not affect the observed features in the CU study (Figure 5C).

3.2. Metabolome—Bile Acids

Overall, the relative abundance of several bile acids was lower in the prebiotic diet groups compared to the control diet groups and these results were directionally consistent across the study sites (Figure 6; see Table 3 with statistics for all bile acids identified). Specifically, the primary bile acid muricholic beta was lower in the prebiotic diet groups in both studies (Figure 6A; significant main effect at CU, significant interaction at NW). The secondary bile acids, deoxycholic acid (Figure 6B; significant main effects) and lithocholic acid (Figure 6C; significant main effect at CU, significant interaction at NW), were lower in the prebiotic diet groups between the study sites. In the CU study, ursodeoxycholic acid was also lower in the prebiotic diet group (significant main effect) but was unaffected in the NW study (Figure 6D; see Table 3). Finally, the secondary conjugated bile acid, glycodeoxycholic acid, was impacted by the prebiotic diet in the NW study and unaffected in the CU study (see Figure 6E for the results of the post hoc analyses). Table 3 lists the bile acids not affected by the prebiotic diet when corrected for multiple comparisons.

3.3. PICRUSt2—Pathways

In both studies, the prebiotic diet consistently affected the inferred functional metabolic pathways PWY-7332, PWY-7090, PWY-6572, and PWY-6545 over time (Figure 7). The superpathway UDP-N-acetylglucosamine-derived O-antigen building blocks biosynthesis or PWY-7332 was significantly higher in the prebiotic diet in the NW study (F(2.71,59.49) = 60.04; p < 0.0001; Figure 7A) and in the CU study (F(1,2.62) = 182.60; p < 0.0001; Figure 7A). The UDP-2,3-diaetamido-2,3-dideoxy-α-D-mannuronate biosynthesis or PWY-7090 was significantly higher in the prebiotic diet in the NW study (F(1,2.81) = 71.76; p < 0.0001; Figure 7B) and in the CU study (F(1,2.47) = 132.21; p < 0.0001; Figure 7B). The chondroitin sulfate degradation I (bacterial) or PWY-6572 was also significantly higher in the prebiotic diet in the NW study (F(1,2.79) = 56.76; p < 0.0001; Figure 7C) and in the CU study (F(1,2.80) = 43.33; p < 0.0001; Figure 7C). Finally, the prebiotic diet increased the pyrimidine deoxyribonucleotides de novo biosynthesis III or PWY-6545 in the NW study (F(1,2.87) = 17.29; p < 0.0001; Figure 7D) and in the CU study (F(1,2.74) = 30.10; p < 0.0001; Figure 7D). There was also significant diet-by-time interactions for PWY-7332 in the NW study (F(2.71,59.49) = 11.90; p < 0.0001; see Figure 7A for the results of the post hoc comparisons) and the CU study (F(2.62,69.98) = 16.90; p < 0.0001; see Figure 7A for the results of the post hoc comparisons); for PWY-7090 in the NW study (F(2.81,57.58) = 9.13; p < 0.0001; see Figure 7B for the results of the post hoc comparisons) and the CU study (F(2.47,62.56) = 15.61; p < 0.0001; see Figure 7B for the results of the post hoc comparisons); for PWY-6572 in the NW study (F(2.79,58.32) = 4.05; p = 0.008; see Figure 7C for the results of the post hoc comparisons) and the CU study (F(2.80,67.23) = 12.36; p < 0.0001; see Figure 7C for the results of the post hoc comparisons); and for PWY-6545 in the NW study (F(2.87,57.55) = 3.47; p = 0.017; see Figure 7D for the results of the post hoc comparisons) and the CU study (F(2.74,69.77) = 16.90; p < 0.0001; see Figure 7D for the results of the post hoc comparisons).

3.4. Correlation Network Analysis

Correlation network analysis was performed to examine the similarities in the prebiotic diet effects between the study sites with output from the network analyses, as shown in Figure 8. The input into the networks were bile acids (Figure 6) and inferred pathways (Figure 7), which were significantly affected by the prebiotic diets between the study sites. There were no consistent correlations between the pathways and bile acids in the control diets across the studies (Figure 8A,B). In contrast, there were consistent correlation networks between the inferred pathways and the bile acid data in the prebiotic diet groups (Figure 8C,D). The prebiotic diet groups had consistent negative correlations between deoxycholic acid and the four inferred pathways (Figure 8). There was also a consistent positive correlation between lithocholic acid and beta muricholic acid beta between the study sites in the prebiotic diet groups. One difference, however, was a negative correlation between deoxycholic acid and lithocholic acid at NW (Figure 8C), but a positive correlation between these two bile acids at the CU study site (Figure 8D).

4. Discussion

The ingestion of a diet enriched in GOS/PDX produces dynamic and robust changes in the gut microbial composition and microbially dependent bile acids. Despite differences in research personnel, animal facilities, geographic locations, elevations, and animal sources, the temporal pattern of changes in the microbial community structure, microbially dependent metabolites, and functional metabolic pathways, was replicated between the study sites. The prebiotic diet also modulated the relative abundance of several genera, reduced microbially modified bile acids, and altered the networks between inferred functional microbial pathways and microbially modified gut bile acids. Importantly, these changes were sufficiently robust to overcome potential environmental differences between the studies.
Based on measures of β-diversity (UniFrac distance), which take into account phylogenetic relationships [47], dietary prebiotics changed both the weighted and unweighted UniFrac distance at both study sites. In the CU study, the weighted UniFrac distance was altered after 2 days on a prebiotic diet, suggesting the rapid growth of higher abundance genera. Dietary prebiotics produced significant compositional changes in the α-diversity metrics (evenness, Faith’s phylogenetic diversity, observed features) at both study sites; however, the metrics of the induced changes were different. In the NW study, prebiotics reduced Faith’s phylogenetic diversity and observed species, whereas in the CU study, prebiotics increased the evenness. These variable impacts of a prebiotic diet on α-diversity between the study sites could reflect inherent differences in the starting microbiomes between the study sites.
The consumption of a diet enriched in GOS/PDX at NW and CU increased the relative abundance of the Bacteroides genus. Based on the ASV and prior shotgun sequencing data from a subset of these samples, Bacteroides uniformis, a member of the Bacteroides genus, was also significantly increased (p = 0.0003) by GOS/PDX [22]. The ingestion of Bacteroides uniformis produces metabolic, immune, and exercise endurance benefits [48,49]. These studies support the idea that an increased relative abundance of specific taxa within the Bacteroides genus may be health promoting.
The consumption of a diet enriched in GOS/PDX also increased the relative abundance of the Parabacteroides genus. The Parabacteroides genus has been shown to be decreased with a high-fat diet and increased with exercise [50]. Parabacteroides distasonis is a species within the Parabacteroides genus. GOS/PDX supplementation increases Parabacteroides distasonis and restores disturbed sleep and circadian rhythm [21,22]. Based on these studies, increases in the relative abundance of specific taxa within the Parabacteroides genus may be health promoting. Importantly, however, Parabacteroides growth left unchecked or not kept in balance could be detrimental to the gut microbial ecosystem [51].
Additional changes to the gut microbial ecosystem include prebiotic-induced increases in the specific taxa within the genera Incertae_Sedis (formerly known as Ruminiclostridium V based on ASV) and the Ruminococcus gauvreauii group. Increases in the levels of Ruminiclostridium V subsequent to the administration of PDX are associated with improved cognitive performance [52]. And, in contrast, low levels of Ruminiclostridium V have been reported for people with kidney stones [53] and rats with acute necrotizing pancreatitis [54]. Consistent with our data, the genera Ruminococcus gauvreauii group is increased by fructooligosaccharides [24], and this genus is lower in individuals with obesity [55], coronary artery disease [56] and Parkinson’s disease [57]. These findings taken together, therefore, suggest that the genera Incertae_Sedis (formerly known as Ruminiclostridium V based on ASV) and Ruminococcus gauvreauii group may be health promoting.
The genus UCG-007 was also increased over time similarly between studies, but little is known about it other than that it varies seasonally [58]. The genera Clostridia_UCG-014, Christensenellaceae_R-7_group, Parasutterella, and Lachnospiraceae_UCG-006 were also all elevated due to the prebiotic diet, but the temporal effects on these genera were less consistent between the study sites.
In addition to increases in the relative abundance of health-promoting genera, several genera were reduced by the prebiotic diet. Most notably, the genus Lachnospiraceae_NK4A136_group was consistently lower in the prebiotic diet groups at both study sites and has recently been implicated in gut mucous membrane function [59]. The genus UCG-005, within the Oscillospiraceae family, was lower in the prebiotic diet groups. This lower relative abundance of UCG-005 may be health promoting given that UCG-005 is elevated in diabetes patients and is associated with elevated uric acid [60]. The genus Eubacterium_fissicatena_group was lower in the prebiotic diet groups and is potentially harmful to bone mineral density [61] and correlates with obesity in a high-fat diet model [62]. The prebiotic diet also lowered Eubacterium_ruminantium_group, GCA-900066575, and Rikenellaceae_RC9_gut_group. Less is known about how and if these genera are related to host health.
Not only did prebiotics change the microbial composition of the gut microbiome, but they also impacted specific features of the gut metabolome. The sequencing data were analyzed using PICRUSt2 and annotated with the MetaCyc metabolic pathway database. These analyses identified four inferred functional metabolic pathways that were changed by the prebiotic diet. Importantly, the prebiotic diet impacted the same pathways between the study sites, with remarkably similar time courses. The first pathway, the UDP-sugar superpathway (PWY-7332), is involved in building the O-antigen polysaccharide for gram-negative bacteria, including Parabacteroides distasonis, which is a component of lipopolysaccharide. The second pathway, the UDP mannuronate biosynthesis pathway (PWY-7090), was identified for both study sites and is involved in UDP-sugar metabolism. Clearly, the consumption of dietary prebiotics affected the UDP-sugar pathway. The third pathway affected by prebiotics, chondroitin sulfate degradation I (PWY-6572), is involved in the degradation of chondroitin sulfate, which is a sulfated glycosaminoglycan that can affect the gut microbiome composition [63] and increase fecal butyrate levels in stressed mice [64]. The fourth pathway was pyrimidine DNA biosynthesis III (PWY-6545), which is involved in the biosynthesis of the activated precursors of DNA/RNA.
While the significance of how the gut metabolome and host physiology are affected by changes in these inferred pathways cannot be deduced from the PICRUSt2 analysis, there is evidence that consumption of GOS/PDX facilitated host sleep/circadian recovery after stressor exposure [21,22]. Clearly, the consumption of GOS/PDX consistently affected these four functional metabolic pathways between the study sites similarly over time. These findings support the idea that dietary prebiotics consistently and similarly altered the micro-ecosystem of the gut microbiome.
A prebiotic diet changes specific gut metabolites with bioactive potential, including microbially modified secondary bile acids [21,22]. Prebiotic diet consumption produced similar decreases in fecal deoxycholic acid and lithocholic acid between the study sites. It has been demonstrated that the consumption of a diet enriched in isomaltulose [23] and fructooligosaccharide [24] prebiotics also reduces fecal lithocholic and deoxycholic acid. In contrast, a high-fat diet increases both fecal deoxycholic acid and intestinal inflammation [65]. Here, we report that GOS/PDX reduces fecal deoxycholic acid and lithocholic acid, and this finding was consistent at both study sites. The current data and prior studies support the conclusion that the consumption of a prebiotic diet reduces fecal bile acids and changes the micro-ecosystem of the gut, similarly.
Our findings indicate that the consumption of a prebiotic diet consistently affects functional metabolic pathways and fecal bile acid profiles. We conducted network correlation analyses between functional metabolic pathways and fecal bile acids to determine whether these changes are related. Correlational networks between pathways and bile acids were not observed in the control diet groups. However, network correlations were found in the prebiotic groups. Specifically, bile acids were significantly correlated with the functional metabolic pathways. The network correlations in both the prebiotic diet groups were remarkably similar between the study sites, with what appears to be a network hub related to deoxycholic acid. Based on these findings and previous work, we hypothesize that decreases in deoxycholic acid may be a key metabolic feature underlying the potential health-promoting effects of GOS/PDX. Deoxycholic acid can bind to the Takeda G protein-coupled receptor 5 (TGR5), which is specific to bile acids and is known to activate several intracellular signaling pathways [66,67].

5. Conclusions

We demonstrate that dietary GOS/PDX produces robust and reproducible changes in the microbial composition of the gut micro-ecosystem, sufficient to overcome unforeseen environmental impacts, addressing a gap in the literature [10,11,12,13]. Although some variations between the NW study and the CU study exist, the consistent pattern of taxonomic changes over time and impacts on functional metabolic pathways are similar. We identified consistent correlational networks associating the changes in bile acids and functional pathways, which supports the robust nature of the effects. Notably, the networks were found in the prebiotic groups and not the control diet groups, supporting the conclusion that the changes are driven by prebiotics. Finally, these key findings were reproduced at both study sites. Overall, a prebiotic diet increases and decreases the relative abundance of several genera, which may support a health-promoting gut micro-ecosystem.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu16111790/s1, Figure S1: Phylum-level relative abundance data demonstrating consistent effects of prebiotic diet on the two main phyla: Firmicutes and Bacteroidetes. (A) There was a significant main effect of prebiotic diet on Firmicutes at NW (F(1,2.77) = 9.72; p = 0.002) and a significant time-by-diet interaction (F(1,59.99) = 4.61; p = 0.004). There was also a significant main effect of prebiotic diet on Firmicutes at CU (F(1,2.85) = 20.94; p = 0.0000078) but there was not a significant time-by-diet interaction. (B) There was also a significant main effect of prebiotic diet at on Bacteroidetes at NW (F(1,2.73) = 6.01; p = 0.014) and a significant time-by-diet interaction (F(1,59.76) = 3.33; p = 0.022). At CU, there was also a significant main effect of prebiotic diet on Bacteroidetes (F(1,2.83) = 20.94; p = 0.0000047) and there was a significant time-by-diet interaction (F(1,71.01) = 3.39; p = 0.019). Notably, in both studies, the relative abundance of Firmicutes increased and Bacteroidetes decreased across time regardless of diet. * p < 0.05 when compared to control diet; Figure S2: Data demonstrating the effects of dietary prebiotics on lower abundance genera between study sites across time. * p < 0.05 when compared to control diet; Figure S3: Data demonstrating inconsistent effects on six different genera between study sites across time. * p < 0.05 when compared to control diet; Table S1: PERMANOVA table demonstrating significant effects of prebiotic diet by time point at both study sites. Numbers represent days on diet; Table S2: Nonparametric longitudinal data (nparLD) Table: ANOVA-type statistic (ATS) showing the similar significant effects of a prebiotic diet across time on genera identified through ANCOM between study sites.

Author Contributions

Conceptualization, R.S.T., C.A.L., P.C.D., M.H.V., F.W.T., R.K., K.P.W.J. and M.F.; Methodology, S.J.B., S.H., T.K., A.G., P.C.D., F.W.T., R.K. and M.F.; Formal analysis, R.S.T. and S.J.B.; Investigation, S.J.B., F.V., S.H., T.K. and A.G.; Resources, K.P.W.J.; Data curation, R.S.T., F.V. and A.G.; Writing—original draft, R.S.T.; Writing—review & editing, S.J.B., F.V., S.H., T.K., A.G., C.A.L., P.C.D., M.H.V., F.W.T., R.K., K.P.W.J. and M.F.; Supervision, R.S.T. and M.F.; Project administration, R.S.T., S.J.B., F.V., M.H.V. and M.F.; Funding acquisition, C.A.L., P.C.D., M.H.V., F.W.T., R.K., K.P.W.J. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

Department of the Navy, Office of Naval Research Multidisciplinary University Research Initiative (MURI) Award, award number N00014-15-1-2809.

Institutional Review Board Statement

CU ONR protocol—Protocol Number: 2361-14MAY2018-DT. Approval Date: 20 November 2017. CU Animal Welfare Assurance Number: D16-00388. NW ONR protocol—Protocol Number: IS00001718. Approval Date: 15 June 2015. NW Animal Welfare Assurance Number: A3283-01.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

Pieter C. Dorrestein is an advisor and holds equity in Sirenas and Cybele, consulted for MSD animal health in 2023. He is a co-founder, scientific advisor, and holds equity in Ometa Labs, Arome, and Enveda with prior approval by UC San Diego. Rob Knight is a scientific advisory board member, and consultant for BiomeSense, Inc., has equity and receives income. He is a scientific advisory board member and has equity in GenCirq. He is a consultant and scientific advisory board member for DayTwo, and receives income. He has equity in and acts as a consultant for Cybele. He is a co-founder of Biota, Inc., and has equity. He is a cofounder of Micronoma, and has equity, and is a scientific advisory board member. Christopher A. Lowry is a co-founder, board member, and Chief Scientific Officer of Mycobacteria Therapeutics Corporation. The remaining authors have no known competing financial interests.

References

  1. Zhang, Y.; Wang, R. The human gut phageome: Composition, development, and alterations in disease. Front. Microbiol. 2023, 14, 1213625. [Google Scholar] [CrossRef]
  2. Zhang, X.; Meng, H.; Hu, X.; Yuan, Z. Diversity and functional profile of gut symbiotic bacteria between Lysinibacillus sphaericus C(3)-41 susceptible and resistant Culex quinquefasciatus Say as revealed by 16S rRNA gene high-throughput sequencing. Front. Microbiol. 2022, 13, 991105. [Google Scholar] [CrossRef]
  3. Naya-Catala, F.; Piazzon, M.C.; Calduch-Giner, J.A.; Sitja-Bobadilla, A.; Perez-Sanchez, J. Diet and Host Genetics Drive the Bacterial and Fungal Intestinal Metatranscriptome of Gilthead Sea Bream. Front. Microbiol. 2022, 13, 883738. [Google Scholar] [CrossRef]
  4. Vincenot, C.E.; Giannino, F.; Rietkerk, M.; Moriya, K.; Mazzoleni, S. Theoretical considerations on the combined use of System Dynamics and individual-based modeling in ecology. Ecol. Model. 2011, 222, 210–218. [Google Scholar] [CrossRef]
  5. Molly, K.; Vandewoestyne, M.; Desmet, I.; Verstraete, W. Validation of the Simulator of the Human Intestinal Microbial Ecosystem (Shime) Reactor Using Microorganism-Associated Activities. Microb. Ecol. Health D 1994, 7, 191–200. [Google Scholar] [CrossRef]
  6. Liu, S.; Qin, P.; Wang, J. High-Fat Diet Alters the Intestinal Microbiota in Streptozotocin-Induced Type 2 Diabetic Mice. Microorganisms 2019, 7, 176. [Google Scholar] [CrossRef]
  7. Liu, X.; Blouin, J.M.; Santacruz, A.; Lan, A.; Andriamihaja, M.; Wilkanowicz, S.; Benetti, P.H.; Tome, D.; Sanz, Y.; Blachier, F.; et al. High-protein diet modifies colonic microbiota and luminal environment but not colonocyte metabolism in the rat model: The increased luminal bulk connection. Am. J. Physiol. Gastrointest. Liver Physiol. 2014, 307, G459–G470. [Google Scholar] [CrossRef]
  8. Barra, N.G.; Anhe, F.F.; Cavallari, J.F.; Singh, A.M.; Chan, D.Y.; Schertzer, J.D. Micronutrients impact the gut microbiota and blood glucose. J. Endocrinol. 2021, 250, R1–R21. [Google Scholar] [CrossRef]
  9. Gibson, G.R.; Hutkins, R.; Sanders, M.E.; Prescott, S.L.; Reimer, R.A.; Salminen, S.J.; Scott, K.; Stanton, C.; Swanson, K.S.; Cani, P.D.; et al. Expert consensus document: The International Scientific Association for Probiotics and Prebiotics (ISAPP) consensus statement on the definition and scope of prebiotics. Nat. Rev. Gastroenterol. Hepatol. 2017, 14, 491–502. [Google Scholar] [CrossRef]
  10. Schloss, P.D. Identifying and Overcoming Threats to Reproducibility, Replicability, Robustness, and Generalizability in Microbiome Research. mBio 2018, 9, e00525-18. [Google Scholar] [CrossRef]
  11. Vitek, J.; Kalibera, T. R-3—Repeatability, Reproducibility and Rigor. Acm Sigplan Not. 2012, 47, 30–36. [Google Scholar] [CrossRef]
  12. Eaton, K.; Pirani, A.; Snitkin, E.S.; Reproducibility Project: Cancer Biology. Replication Study: Intestinal inflammation targets cancer-inducing activity of the microbiota. eLife 2018, 7, e34364. [Google Scholar] [CrossRef]
  13. Repass, J.; Reproducibility Project: Cancer Biology. Replication Study: Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. eLife 2018, 7, e25801. [Google Scholar] [CrossRef]
  14. Panek, M.; Cipcic Paljetak, H.; Baresic, A.; Peric, M.; Matijasic, M.; Lojkic, I.; Vranesic Bender, D.; Krznaric, Z.; Verbanac, D. Methodology challenges in studying human gut microbiota—Effects of collection, storage, DNA extraction and next generation sequencing technologies. Sci. Rep. 2018, 8, 5143. [Google Scholar] [CrossRef]
  15. De Filippo, C.; Cavalieri, D.; Di Paola, M.; Ramazzotti, M.; Poullet, J.B.; Massart, S.; Collini, S.; Pieraccini, G.; Lionetti, P. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. USA 2010, 107, 14691–14696. [Google Scholar] [CrossRef]
  16. Gupta, V.K.; Paul, S.; Dutta, C. Geography, Ethnicity or Subsistence-Specific Variations in Human Microbiome Composition and Diversity. Front. Microbiol. 2017, 8, 1162. [Google Scholar] [CrossRef]
  17. Kumari, M.; Bhushan, B.; Eslavath, M.R.; Srivastava, A.K.; Meena, R.C.; Varshney, R.; Ganju, L. Impact of high altitude on composition and functional profiling of oral microbiome in Indian male population. Sci. Rep. 2023, 13, 4038. [Google Scholar] [CrossRef]
  18. Dong, W.; Ma, L.; Huang, Q.; Yang, X.; Mei, Z.; Kong, M.; Sun, Z.; Zhang, Z.; Li, J.; Zou, J.; et al. Gut microbiome alterations in pulmonary hypertension in highlanders and lowlanders. ERJ Open Res. 2023, 9, 00617-2022. [Google Scholar] [CrossRef]
  19. Rothschild, D.; Weissbrod, O.; Barkan, E.; Kurilshikov, A.; Korem, T.; Zeevi, D.; Costea, P.I.; Godneva, A.; Kalka, I.N.; Bar, N.; et al. Environment dominates over host genetics in shaping human gut microbiota. Nature 2018, 555, 210–215. [Google Scholar] [CrossRef]
  20. Ericsson, A.C.; Franklin, C.L. The gut microbiome of laboratory mice: Considerations and best practices for translational research. Mamm. Genome 2021, 32, 239–250. [Google Scholar] [CrossRef]
  21. Bowers, S.J.; Summa, K.C.; Thompson, R.S.; Gonzalez, A.; Vargas, F.; Olker, C.; Jiang, P.; Lowry, C.A.; Dorrestein, P.C.; Knight, R.; et al. A Prebiotic Diet Alters the Fecal Microbiome and Improves Sleep in Response to Sleep Disruption in Rats. Front. Neurosci. 2022, 16, 889211. [Google Scholar] [CrossRef] [PubMed]
  22. Thompson, R.S.; Gaffney, M.; Hopkins, S.; Kelley, T.; Gonzalez, A.; Bowers, S.J.; Vitaterna, M.H.; Turek, F.W.; Foxx, C.L.; Lowry, C.A.; et al. Ruminiclostridium 5, Parabacteroides distasonis, and bile acid profile are modulated by prebiotic diet and associate with facilitated sleep/clock realignment after chronic disruption of rhythms. Brain Behav. Immun. 2021, 97, 150–166. [Google Scholar] [CrossRef] [PubMed]
  23. Yang, Z.D.; Guo, Y.S.; Huang, J.S.; Gao, Y.F.; Peng, F.; Xu, R.Y.; Su, H.H.; Zhang, P.J. Isomaltulose Exhibits Prebiotic Activity, and Modulates Gut Microbiota, the Production of Short Chain Fatty Acids, and Secondary Bile Acids in Rats. Molecules 2021, 26, 2464. [Google Scholar] [CrossRef] [PubMed]
  24. Phungviwatnikul, T.; Lee, A.H.; Belchik, S.E.; Suchodolski, J.S.; Swanson, K.S. Weight loss and high-protein, high-fiber diet consumption impact blood metabolite profiles, body composition, voluntary physical activity, fecal microbiota, and fecal metabolites of adult dogs. J. Anim. Sci. 2022, 100, skab379. [Google Scholar] [CrossRef]
  25. McMillin, M.; DeMorrow, S. Effects of bile acids on neurological function and disease. FASEB J. 2016, 30, 3658–3668. [Google Scholar] [CrossRef] [PubMed]
  26. Perino, A.; Demagny, H.; Velazquez-Villegas, L.A.; Schoonjans, K. Molecular Physiology of Bile Acid Signaling in Health, Disease and Aging. Physiol. Rev. 2020, 101, 683–731. [Google Scholar] [CrossRef] [PubMed]
  27. Thompson, R.S.; Roller, R.; Mika, A.; Greenwood, B.N.; Knight, R.; Chichlowski, M.; Berg, B.M.; Fleshner, M. Dietary Prebiotics and Bioactive Milk Fractions Improve NREM Sleep, Enhance REM Sleep Rebound and Attenuate the Stress-Induced Decrease in Diurnal Temperature and Gut Microbial Alpha Diversity. Front. Behav. Neurosci. 2017, 10, 240. [Google Scholar] [CrossRef] [PubMed]
  28. Mika, A.; Gaffney, M.; Roller, R.; Hills, A.; Bouchet, C.A.; Hulen, K.A.; Thompson, R.S.; Chichlowski, M.; Berg, B.M.; Fleshner, M. Feeding the developing brain: Juvenile rats fed diet rich in prebiotics and bioactive milk fractions exhibit reduced anxiety-related behavior and modified gene expression in emotion circuits. Neurosci. Lett. 2018, 677, 103–109. [Google Scholar] [CrossRef] [PubMed]
  29. Mika, A.; Day, H.E.; Martinez, A.; Rumian, N.L.; Greenwood, B.N.; Chichlowski, M.; Berg, B.M.; Fleshner, M. Early life diets with prebiotics and bioactive milk fractions attenuate the impact of stress on learned helplessness behaviours and alter gene expression within neural circuits important for stress resistance. Eur. J. Neurosci. 2016, 45, 342–357. [Google Scholar] [CrossRef]
  30. Bowers, S.J.; Vargas, F.; Gonzalez, A.; He, S.; Jiang, P.; Dorrestein, P.C.; Knight, R.; Wright, K.P., Jr.; Lowry, C.A.; Fleshner, M.; et al. Repeated sleep disruption in mice leads to persistent shifts in the fecal microbiome and metabolome. PLoS ONE 2020, 15, e0229001. [Google Scholar] [CrossRef]
  31. Krauth, S.J.; Coulibaly, J.T.; Knopp, S.; Traore, M.; N’Goran, E.K.; Utzinger, J. An in-depth analysis of a piece of shit: Distribution of Schistosoma mansoni and hookworm eggs in human stool. PLoS Negl. Trop. Dis. 2012, 6, e1969. [Google Scholar] [CrossRef]
  32. Apprill, A.; McNally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef]
  33. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. [Google Scholar] [CrossRef]
  34. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857, Corrected in Nat. Biotechnol. 2019, 37, 1091. [Google Scholar] [CrossRef] [PubMed]
  35. Yilmaz, P.; Parfrey, L.W.; Yarza, P.; Gerken, J.; Pruesse, E.; Quast, C.; Schweer, T.; Peplies, J.; Ludwig, W.; Glockner, F.O. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014, 42, D643–D648. [Google Scholar] [CrossRef] [PubMed]
  36. The Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature 2012, 486, 207–214. [Google Scholar] [CrossRef]
  37. Lozupone, C.; Lladser, M.E.; Knights, D.; Stombaugh, J.; Knight, R. UniFrac: An effective distance metric for microbial community comparison. ISME J. 2011, 5, 169–172. [Google Scholar] [CrossRef]
  38. Faith, D.P. Phylogenetic pattern and the quantification of organismal biodiversity. Philos. Trans. R. Soc. Lond. B Biol. Sci. 1994, 345, 45–58. [Google Scholar] [CrossRef] [PubMed]
  39. Mandal, S.; Van Treuren, W.; White, R.A.; Eggesbo, M.; Knight, R.; Peddada, S.D. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health Dis. 2015, 26, 27663. [Google Scholar] [CrossRef]
  40. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017, 11, 2639–2643. [Google Scholar] [CrossRef]
  41. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glockner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef] [PubMed]
  42. Douglas, G.M.; Maffei, V.J.; Zaneveld, J.R.; Yurgel, S.N.; Brown, J.R.; Taylor, C.M.; Huttenhower, C.; Langille, M.G.I. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 2020, 38, 685–688. [Google Scholar] [CrossRef] [PubMed]
  43. Melnik, A.V.; da Silva, R.R.; Hyde, E.R.; Aksenov, A.A.; Vargas, F.; Bouslimani, A.; Protsyuk, I.; Jarmusch, A.K.; Tripathi, A.; Alexandrov, T.; et al. Coupling Targeted and Untargeted Mass Spectrometry for Metabolome-Microbiome-Wide Association Studies of Human Fecal Samples. Anal. Chem. 2017, 89, 7549–7559. [Google Scholar] [CrossRef] [PubMed]
  44. Thompson, R.S.; Vargas, F.; Dorrestein, P.C.; Chichlowski, M.; Berg, B.M.; Fleshner, M. Dietary prebiotics alter novel microbial dependent fecal metabolites that improve sleep. Sci. Rep. 2020, 10, 3848. [Google Scholar] [CrossRef] [PubMed]
  45. Kelly, B.J.; Gross, R.; Bittinger, K.; Sherrill-Mix, S.; Lewis, J.D.; Collman, R.G.; Bushman, F.D.; Li, H. Power and sample-size estimation for microbiome studies using pairwise distances and PERMANOVA. Bioinformatics 2015, 31, 2461–2468. [Google Scholar] [CrossRef] [PubMed]
  46. Tang, Z.Z.; Chen, G.; Alekseyenko, A.V. PERMANOVA-S: Association test for microbial community composition that accommodates confounders and multiple distances. Bioinformatics 2016, 32, 2618–2625. [Google Scholar] [CrossRef] [PubMed]
  47. Lozupone, C.; Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [PubMed]
  48. Lopez-Almela, I.; Romani-Perez, M.; Bullich-Vilarrubias, C.; Benitez-Paez, A.; Gomez Del Pulgar, E.M.; Frances, R.; Liebisch, G.; Sanz, Y. Bacteroides uniformis combined with fiber amplifies metabolic and immune benefits in obese mice. Gut Microbes 2021, 13, 1–20. [Google Scholar] [CrossRef] [PubMed]
  49. Morita, H.; Kano, C.; Ishii, C.; Kagata, N.; Ishikawa, T.; Hirayama, A.; Uchiyama, Y.; Hara, S.; Nakamura, T.; Fukuda, S. Bacteroides uniformis and its preferred substrate, alpha-cyclodextrin, enhance endurance exercise performance in mice and human males. Sci. Adv. 2023, 9, eadd2120. [Google Scholar] [CrossRef]
  50. Carbajo-Pescador, S.; Porras, D.; Garcia-Mediavilla, M.V.; Martinez-Florez, S.; Juarez-Fernandez, M.; Cuevas, M.J.; Mauriz, J.L.; Gonzalez-Gallego, J.; Nistal, E.; Sanchez-Campos, S. Beneficial effects of exercise on gut microbiota functionality and barrier integrity, and gut-liver crosstalk in an in vivo model of early obesity and non-alcoholic fatty liver disease. Dis. Model. Mech. 2019, 12, dmm039206. [Google Scholar] [CrossRef]
  51. Ezeji, J.C.; Sarikonda, D.K.; Hopperton, A.; Erkkila, H.L.; Cohen, D.E.; Martinez, S.P.; Cominelli, F.; Kuwahara, T.; Dichosa, A.E.K.; Good, C.E.; et al. Parabacteroides distasonis: Intriguing aerotolerant gut anaerobe with emerging antimicrobial resistance and pathogenic and probiotic roles in human health. Gut Microbes 2021, 13, 1922241. [Google Scholar] [CrossRef] [PubMed]
  52. Berding, K.; Long-Smith, C.M.; Carbia, C.; Bastiaanssen, T.F.S.; van de Wouw, M.; Wiley, N.; Strain, C.R.; Fouhy, F.; Stanton, C.; Cryan, J.F.; et al. A specific dietary fibre supplementation improves cognitive performance-an exploratory randomised, placebo-controlled, crossover study. Psychopharmacology 2021, 238, 149–163. [Google Scholar] [CrossRef] [PubMed]
  53. Tang, R.; Jiang, Y.; Tan, A.; Ye, J.; Xian, X.; Xie, Y.; Wang, Q.; Yao, Z.; Mo, Z. 16S rRNA gene sequencing reveals altered composition of gut microbiota in individuals with kidney stones. Urolithiasis 2018, 46, 503–514. [Google Scholar] [CrossRef] [PubMed]
  54. Chen, J.; Huang, C.; Wang, J.; Zhou, H.; Lu, Y.; Lou, L.; Zheng, J.; Tian, L.; Wang, X.; Cao, Z.; et al. Dysbiosis of intestinal microbiota and decrease in paneth cell antimicrobial peptide level during acute necrotizing pancreatitis in rats. PLoS ONE 2017, 12, e0176583. [Google Scholar] [CrossRef]
  55. Verheggen, R.; Konstanti, P.; Smidt, H.; Hermus, A.; Thijssen, D.H.J.; Hopman, M.T.E. Eight-week exercise training in humans with obesity: Marked improvements in insulin sensitivity and modest changes in gut microbiome. Obesity 2021, 29, 1615–1624. [Google Scholar] [CrossRef] [PubMed]
  56. Toya, T.; Corban, M.T.; Marrietta, E.; Horwath, I.E.; Lerman, L.O.; Murray, J.A.; Lerman, A. Coronary artery disease is associated with an altered gut microbiome composition. PLoS ONE 2020, 15, e0227147. [Google Scholar] [CrossRef] [PubMed]
  57. Babacan Yildiz, G.; Kayacan, Z.C.; Karacan, I.; Sumbul, B.; Elibol, B.; Gelisin, O.; Akgul, O. Altered gut microbiota in patients with idiopathic Parkinson’s disease: An age-sex matched case-control study. Acta Neurol. Belg. 2023, 123, 999–1009. [Google Scholar] [CrossRef] [PubMed]
  58. Yildirim, E.; Ilina, L.; Laptev, G.; Filippova, V.; Brazhnik, E.; Dunyashev, T.; Dubrovin, A.; Novikova, N.; Tiurina, D.; Tarlavin, N.; et al. The structure and functional profile of ruminal microbiota in young and adult reindeers (Rangifer tarandus) consuming natural winter-spring and summer-autumn seasonal diets. PeerJ 2021, 9, e12389. [Google Scholar] [CrossRef] [PubMed]
  59. Ma, L.; Ni, Y.; Wang, Z.; Tu, W.; Ni, L.; Zhuge, F.; Zheng, A.; Hu, L.; Zhao, Y.; Zheng, L.; et al. Spermidine improves gut barrier integrity and gut microbiota function in diet-induced obese mice. Gut Microbes 2020, 12, 1–19. [Google Scholar] [CrossRef]
  60. Zhang, W.; Wang, T.; Guo, R.; Cui, W.; Yu, W.; Wang, Z.; Jiang, Y.; Jiang, M.; Wang, X.; Liu, C.; et al. Variation of Serum Uric Acid Is Associated With Gut Microbiota in Patients With Diabetes Mellitus. Front. Cell Infect. Microbiol. 2021, 11, 761757. [Google Scholar] [CrossRef]
  61. Wang, Y.; Zhang, X.; Tang, G.; Deng, P.; Qin, Y.; Han, J.; Wang, S.; Sun, X.; Li, D.; Chen, Z. The causal relationship between gut microbiota and bone mineral density: A Mendelian randomization study. Front. Microbiol. 2023, 14, 1268935. [Google Scholar] [CrossRef] [PubMed]
  62. Song, Y.; Shen, H.; Liu, T.; Pan, B.; De Alwis, S.; Zhang, W.; Luo, X.; Li, Z.; Wang, N.; Ma, W.; et al. Effects of three different mannans on obesity and gut microbiota in high-fat diet-fed C57BL/6J mice. Food Funct. 2021, 12, 4606–4620. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, K.; Qi, L.; Zhao, L.; Liu, J.; Guo, Y.; Zhang, C. Degradation of chondroitin sulfate: Mechanism of degradation, influence factors, structure-bioactivity relationship and application. Carbohydr. Polym. 2023, 301, 120361. [Google Scholar] [CrossRef] [PubMed]
  64. Liu, F.; Zhang, N.; Li, Z.; Wang, X.; Shi, H.; Xue, C.; Li, R.W.; Tang, Q. Chondroitin sulfate disaccharides modified the structure and function of the murine gut microbiome under healthy and stressed conditions. Sci. Rep. 2017, 7, 6783. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, L.; Gong, Z.; Zhang, X.; Zhu, F.; Liu, Y.; Jin, C.; Du, X.; Xu, C.; Chen, Y.; Cai, W.; et al. Gut microbial bile acid metabolite skews macrophage polarization and contributes to high-fat diet-induced colonic inflammation. Gut Microbes 2020, 12, 1–20. [Google Scholar] [CrossRef] [PubMed]
  66. Duboc, H.; Tache, Y.; Hofmann, A.F. The bile acid TGR5 membrane receptor: From basic research to clinical application. Dig. Liver Dis. 2014, 46, 302–312. [Google Scholar] [CrossRef]
  67. Pols, T.W.; Noriega, L.G.; Nomura, M.; Auwerx, J.; Schoonjans, K. The bile acid membrane receptor TGR5: A valuable metabolic target. Dig. Dis. 2011, 29, 37–44. [Google Scholar] [CrossRef]
Figure 1. Experimental timeline detailing methods and fecal sampling events. In both studies, animals arrived on postnatal day 23 and were immediately placed on either the control diet or prebiotic diet. In the Northwestern study, fecal samples were taken on experimental (postnatal) days 0 (23), 28 (51), 42 (65), and 51 (74), while in the CU study, fecal samples were taken on experimental days 2 (25), 33 (58), 75 (100), and 94 (119).
Figure 1. Experimental timeline detailing methods and fecal sampling events. In both studies, animals arrived on postnatal day 23 and were immediately placed on either the control diet or prebiotic diet. In the Northwestern study, fecal samples were taken on experimental (postnatal) days 0 (23), 28 (51), 42 (65), and 51 (74), while in the CU study, fecal samples were taken on experimental days 2 (25), 33 (58), 75 (100), and 94 (119).
Nutrients 16 01790 g001
Figure 2. Unweighted and weighted UniFrac distance examining β-diversity of the fecal microbiome between studies. (A) In the NW study, unweighted UniFrac distance at experimental day 0 was not different between the control and prebiotic diets, but was different on subsequent days 28, 42, and 51. (B) In the CU study, unweighted UniFrac distance at experimental day 2 was not different between the control and prebiotic diets, but was different on subsequent days 33, 75, and 94. (C) In the NW study, weighted UniFrac distance was not different on day 0 between the control and prebiotic diets, but was different on the remaining days examined. (D) In the CU study, weighted UniFrac distance was significantly different on day 2 between the control and prebiotic diets, an effect that persisted for days 33, 75, and 94.
Figure 2. Unweighted and weighted UniFrac distance examining β-diversity of the fecal microbiome between studies. (A) In the NW study, unweighted UniFrac distance at experimental day 0 was not different between the control and prebiotic diets, but was different on subsequent days 28, 42, and 51. (B) In the CU study, unweighted UniFrac distance at experimental day 2 was not different between the control and prebiotic diets, but was different on subsequent days 33, 75, and 94. (C) In the NW study, weighted UniFrac distance was not different on day 0 between the control and prebiotic diets, but was different on the remaining days examined. (D) In the CU study, weighted UniFrac distance was significantly different on day 2 between the control and prebiotic diets, an effect that persisted for days 33, 75, and 94.
Nutrients 16 01790 g002
Figure 3. Consumption of the prebiotic diet produced increases in 9 higher abundance genera between studies. There were consistent increases over time due to the prebiotic diet in: (A) Bacteroides, (B) Parabacteroides, (E) Incertae_Sedis (Ruminiclostridium V), (G) Ruminococcus_gauvreauii_group, and (H) UCG-007. While there were prebiotic diet-induced increases in (C) Clostridia_UCG-014, (D) Christensenellaceae_R-7_group, (F) Parasutterella, and (I) Lachnospiraceae_UCG-006, these genera had less consistent increases over time between studies. * p < 0.05 when compared to control diet.
Figure 3. Consumption of the prebiotic diet produced increases in 9 higher abundance genera between studies. There were consistent increases over time due to the prebiotic diet in: (A) Bacteroides, (B) Parabacteroides, (E) Incertae_Sedis (Ruminiclostridium V), (G) Ruminococcus_gauvreauii_group, and (H) UCG-007. While there were prebiotic diet-induced increases in (C) Clostridia_UCG-014, (D) Christensenellaceae_R-7_group, (F) Parasutterella, and (I) Lachnospiraceae_UCG-006, these genera had less consistent increases over time between studies. * p < 0.05 when compared to control diet.
Nutrients 16 01790 g003
Figure 4. Consumption of prebiotic diet led to decreases in 6 higher abundance genera between studies. There were consistent decreases over time due to prebiotic diet consumption in: (A) Lachnospiraceae_NK4A136_group, (C) UCG-005, (E) Eubacterium_fissicatena_group, (F) Eubacterium_ruminantium_group, (G) GCA-900066575, and (I) Rikenellaceae_R9-gut_group. There were less consistent effects due to diet between studies in: (B) Eubacterium_ coprostanoligenes_group, (D) Colidextribacter, and (H) Roseburia. * p < 0.05 when compared to prebiotic diet.
Figure 4. Consumption of prebiotic diet led to decreases in 6 higher abundance genera between studies. There were consistent decreases over time due to prebiotic diet consumption in: (A) Lachnospiraceae_NK4A136_group, (C) UCG-005, (E) Eubacterium_fissicatena_group, (F) Eubacterium_ruminantium_group, (G) GCA-900066575, and (I) Rikenellaceae_R9-gut_group. There were less consistent effects due to diet between studies in: (B) Eubacterium_ coprostanoligenes_group, (D) Colidextribacter, and (H) Roseburia. * p < 0.05 when compared to prebiotic diet.
Nutrients 16 01790 g004
Figure 5. There was a significant main effect of the prebiotic diet, increasing (A) the evenness of the alpha diversity in the CU study. In contrast, the main significant effects of the prebiotic diet involved decreases in both (B) Faith’s phylogenetic diversity and (C) the observed features of the alpha diversity in the NW study. There were no significant time-by-diet interactions in regard to the measures of alpha diversity, except at NW in observed features. * p < 0.05 effect of diet.
Figure 5. There was a significant main effect of the prebiotic diet, increasing (A) the evenness of the alpha diversity in the CU study. In contrast, the main significant effects of the prebiotic diet involved decreases in both (B) Faith’s phylogenetic diversity and (C) the observed features of the alpha diversity in the NW study. There were no significant time-by-diet interactions in regard to the measures of alpha diversity, except at NW in observed features. * p < 0.05 effect of diet.
Nutrients 16 01790 g005
Figure 6. Consumption of dietary prebiotics affected fecal bile acids between studies, including: (A) muricholic acid beta, (B) deoxycholic acid, and (C) lithocholic acid. Moreover, (D) ursodeoxycholic acid was decreased in the CU study, and (E) glycodeoxycholic acid was decreased in the NW study.
Figure 6. Consumption of dietary prebiotics affected fecal bile acids between studies, including: (A) muricholic acid beta, (B) deoxycholic acid, and (C) lithocholic acid. Moreover, (D) ursodeoxycholic acid was decreased in the CU study, and (E) glycodeoxycholic acid was decreased in the NW study.
Nutrients 16 01790 g006
Figure 7. Functional metabolic pathways affected by prebiotic diet, annotated with the MetaCyc metabolic pathway database. Consumption of dietary prebiotics altered the: (A) superpathway UDP-N-acetylglucosamin-derived O-antigen building blocks biosynthesis (PWY-7332) or the UDP-sugar superpathway, the (B) UDP-2,3-diaetamido-2,3-dideoxy-α-D-mannuronate biosynthesis (PWY-7090) or UDP mannuronate pathway, the (C) chondroitin sulfate degradation I (bacterial) pathway (PWY-6572), and the (D) pyrimidine deoxyribonucleotides de novo biosynthesis III pathway (PWY-6545), when compared to the control diet. These effects were consistent between the study sites and over time. * p < 0.05 when compared to control diet.
Figure 7. Functional metabolic pathways affected by prebiotic diet, annotated with the MetaCyc metabolic pathway database. Consumption of dietary prebiotics altered the: (A) superpathway UDP-N-acetylglucosamin-derived O-antigen building blocks biosynthesis (PWY-7332) or the UDP-sugar superpathway, the (B) UDP-2,3-diaetamido-2,3-dideoxy-α-D-mannuronate biosynthesis (PWY-7090) or UDP mannuronate pathway, the (C) chondroitin sulfate degradation I (bacterial) pathway (PWY-6572), and the (D) pyrimidine deoxyribonucleotides de novo biosynthesis III pathway (PWY-6545), when compared to the control diet. These effects were consistent between the study sites and over time. * p < 0.05 when compared to control diet.
Nutrients 16 01790 g007
Figure 8. Network correlations from both study sites, demonstrating consistent networks between inferred functional metabolic pathways and bile acids in prebiotic diet groups. There were no consistent correlation networks present in the control diet groups between the study sites (A,B). The consistent correlation networks in the prebiotic diet groups (C) at NW and (D) at CU imply that the microbially modified secondary bile acid, deoxycholic acid, could be an important component underlying the beneficial effects of dietary prebiotics.
Figure 8. Network correlations from both study sites, demonstrating consistent networks between inferred functional metabolic pathways and bile acids in prebiotic diet groups. There were no consistent correlation networks present in the control diet groups between the study sites (A,B). The consistent correlation networks in the prebiotic diet groups (C) at NW and (D) at CU imply that the microbially modified secondary bile acid, deoxycholic acid, could be an important component underlying the beneficial effects of dietary prebiotics.
Nutrients 16 01790 g008
Table 1. PERMANOVA table demonstrating significant effects of prebiotic diet by time point at both study sites. Numbers represent days on the diet.
Table 1. PERMANOVA table demonstrating significant effects of prebiotic diet by time point at both study sites. Numbers represent days on the diet.
PERMANOVAs (Pseudo-F)
Northwestern
0284251
UnweightedF(2,68) = 1.24; p = 0.154F(2,80) = 7.68; p = 0.001F(2,66) = 5.60; p = 0.001F(2,69) = 5.87; p = 0.001
WeightedF(2,68) = 2.19; p = 0.053F(2,80) = 9.31; p = 0.001F(2,66) = 4.26; p = 0.001F(2,66) = 4.34; p = 0.001
University of Colorado Boulder
2337594
UnweightedF(2,48) = 1.31; p = 0.053F(2,78) = 4.89; p = 0.001F(2,83) = 3.84; p = 0.001F(2,84) = 4.16; p = 0.001
WeightedF(2,48) = 3.97; p = 0.006F(2,78) = 10.99; p = 0.001F(2,83) = 7.39; p = 0.001F(2,84) = 3.93; p = 0.001
Table 2. Nonparametric longitudinal data (nparLD) table: ANOVA-type statistics (ATSs) showing the similar significant effects of a prebiotic diet over time on genera identified through ANCOM, between study sites.
Table 2. Nonparametric longitudinal data (nparLD) table: ANOVA-type statistics (ATSs) showing the similar significant effects of a prebiotic diet over time on genera identified through ANCOM, between study sites.
Nonparametric Longitudinal Data (nparLD) Table: ANOVA-Type Statistics (ATSs)
Genera (Relative Abundance)
Diet—F-Value; p-ValueTime—F-Value; p-ValueDiet × Time
Higher in Prebiotic Diet, color indicates consistent effect across study site
High Relative Abundance (2–20%)
Bacteroides (Figure 3A)NWF(1,2.74) = 12.62p = 0.00038F(2.74,59.91) = 41.26p = 1.47 × 1024F(2.74,59.91) = 2.78p = 0.044
CUF(1,2.76) = 7.42p = 0.0064F(2.76,70.42) = 17.71p = 9.39 × 1011F(2.76,70.42) = 9.90p = 3.63 × 106
Clostridia_UCG-014 (Figure 3C)NWF(1,2.65) = 4.72p = 0.029F(2.65,55.97) = 38.09p = 4.85 × 1022F(2.65,55.97) = 1.01p = 0.379
CUF(1,2.71) = 35.93p = 2.05 × 109F(2.71,67.05) = 12.88p = 8.44 × 108F(2.71,67.05) = 2.18p = 0.095
Christensenellaceae_R7_group (Figure 3D)NWF(1,2.82) = 3.70p = 0.054F(2.82,53.48) = 79.63p = 1.30 × 1048F(2.82,53.48) = 2.59p = 0.055
CUF(1,2.81) = 48.65p = 3.06 × 1012F(2.81,69.23) = 197.98p = 2.63 × 10120F(2.81,69.23) = 9.16p = 8.53 × 106
Incertae_Sedis (Ruminiclostridium V) (Figure 3E)NWF(1,2.81) = 76.70p = 1.99 × 1018F(2.81,56.25) = 46.03p = 4.62 × 1028F(2.81,56.25) = 16.85p = 2.09 × 1010
CUF(1,2.51) = 210.79p = 9.26 × 1048F(2.51,71.88) = 85.26p = 1.03 × 1046F(2.51,71.88) = 19.26p = 7.86 × 1011
Parabacteroides (Figure 3B)NWF(1,2.74) = 158.1p = 2.96 × 1036F(2.74,59.18) = 44.69p = 1.32 × 1026F(2.74,59.18) = 21.04p = 1.19 × 1012
CUF(1,2.74) = 467.75p = 9.91 × 10104F(2.74,71.88) = 71.88p = 6.98 × 1032F(2.74,71.88) = 4.99p = 0.00258
Low Relative Abundance (1–2%)
Parasutterella (Figure 3F)NWF(1,78) = 29.19p =6.57 × 108F(2.78,59.45) = 40.76p = 1.46 × 1024F(2.78,59.45) = 1.09p = 0.127
CUF(1,2.78) = 9.15p = 0.0025F(2.78,71.92) = 63.79p = 1.78 × 1038F(2.78,71.92) = 2.25p = 0.052
Ruminococcus_gauvreauii_group (Figure 3G)NWF(1,2.71) = 104.03p = 1.99 × 1024F(2.71,59.31) = 27.24p = 3.73 × 1016F(2.71,59.31) = 17.95p = 9.61 × 1011
CUF(1,2.33) = 16.93p = 0.000039F(2.33,69.45) = 19.51p = 2.48 × 1010F(2.33,69.45) = 4.79p = 0.00018
UCG-007 (Figure 3H)NWF(1,2.84) = 289.83p = 5.42 × 1065F(2.84,55.73) = 40.50p = 7.18 × 1025F(2.84,55.73) = 31.13p = 3.74 × 1019
CUF(1,2.66) = 140.28p = 2.31 × 1032F(2.66,57.78) = 32.24p = 9.19 × 1019F(2.66,57.78) = 10.11p = 3.89 × 106
Lachnospiraceae_UCG-006 (Figure 3I)NWF(1,2.77) = 1.76p = 0.184F(2.77,59.93) = 42.81p = 9.83 × 1026F(2.77,59.93) = 4.89p = 0.0028
CUF(1,2.77) = 6.33p = 0.0118F(2.77,65.94) = 51.97p = 3.48 × 1031F(2.77,65.94) = 4.61p = 0.00410
Higher in Control Diet, color indicates consistent effect across study site
High Relative Abundance (2–20%)
Lachnospiraceae_NK4A136_group (Figure 4A)NWF(1,2.81) = 36.70p = 1.38 × 109F(2.81,59.99) = 20.53p = 1.34 × 1012F(2.81,59.99) = 2.81p = 0.020
CUF(1,2.83) = 13.13p = 0.0003F(2.83,71.45) = 2.30p = 0.079F(2.83,71.45) = 3.99p = 0.0087
Eubacterium_coprostanoligenes_group (Figure 4B)NWF(1,2.48) = 1.34p = 0.247F(2.48,57.50) = 30.14p = 1.56 × 1016F(2.48,57.50) = 1.66p = 0.183
CUF(1,2.87) = 3.64p = 0.056F(2.87,71.52) = 56.35p = 5.72 × 1035F(2.87,71.52) = 0.55p = 0.638
UCG-005 (Figure 4C)NWF(1,2.82) = 8.07p = 0.0451F(2.82,59.10) = 21.62p = 2.59 × 1013F(2.82,59.10) = 1.62p = 0.184
CUF(1,2.37) = 0.841p = 0.359F(2.37,71.81) = 18.66p = 4.57 × 1010F(2.37,71.81) = 1.73p = 0.171
Low Relative Abundance (1–2%)
Colidextribacter (Figure 4D)NWF(1,2.73) = 13.95p = 0.00019F(2.73,59.62) = 18.55p = 3.78 × 1011F(2.73,59.62) = 0.816p = 0.475
CUF(1,2.78) = 0.013p = 0.911F(2.78,71.96) = 32.26p = 1.64 × 108F(2.78,71.96) = 0.328p = 0.790
Eubacterium_fissicatena_group (Figure 4E)NWF(1,2.41) = 9.51p = 0.002F(2.41,54.79) = 25.73p = 7.31 × 1014F(2.41,54.79) = 7.08p = 0.00034
CUF(1,2.26) = 4.64p = 0.031F(2.26,68.41) = 18.53p = 1.32 × 109F(2.26,68.41) = 3.09p = 0.039
Eubacterium_ruminantium_group (Figure 4F)NWF(1,2.62) = 17.80p = 0.00002F(2.62,38.97) = 7.97p = 0.00006F(2.62,38.97) = 11.57p = 6.82 × 107
CUF(1,2.44) = 6.22p = 0.013F(2.44,63.83) = 8.63p = 0.00005F(2.44,63.83) = 6.31p = 0.0008
GCA-900066575 (Figure 4G)NWF(1,2.92) = 20.93p = 0.000005F(2.92,58.16) = 24.09p = 3.24 × 1015F(2.92,58.16) = 5.18p = 0.0016
CUF(1,2.91) = 9.67p = 0.0019F(2.91,71.91) = 29.78p = 49.98 × 1019F(2.91,71.91) = 0.937p = 0.420
Roseburia (Figure 4H)NWF(1,2.80) = 6.48p = 0.0109F(2.80,59.66) = 2.90p = 0.037F(2.80,59.66) = 1.73p = 0.161
CUF(1,2.74) = 4.79p = 0.029F(2.74,72.00) = 8.71p = 0.000019F(2.74,72.00) = 0.776p = 0.50
Rikenellaceae_RC9_gut_group (Figure 4I)NWF(1,2.74) = 25.70p = 3.99 × 107F(2.74,59.55) = 9.20p = 0.00006F(2.74,59.55) = 5.55p = 0.0012
CUF(1,2.72) = 10.90p = 0.00096F(2.72,70.925) = 2.88p = 0.040F(2.72,70.925) = 1.42p = 0.236
Table 3. Nonparametric longitudinal data (nparLD) table: ANOVA-type statistics (ATSs) showing significant effects of a prebiotic diet over time on all identified bile acids, between study sites.
Table 3. Nonparametric longitudinal data (nparLD) table: ANOVA-type statistics (ATSs) showing significant effects of a prebiotic diet over time on all identified bile acids, between study sites.
Nonparametric Longitudinal Data (naprLD) Table: ANOVA-Type Statistics (ATSs)
Bile Acids
Diet—F-Value; p-Valuep-adj. (Holm)Time—F-Value; p-Valuep-adj. (Holm)Diet × Timep-adj. (Holm)
Color indicates consistent effect across study site
Primary Bile Acids
Cholic AcidNWF(1,2.846) = 0.190p = 0.663n/aF(2.846,58.12) = 9.534p = 4.493 × 106p = 8.98 × 106F(2.846,58.12) = 0.921p = 0.426ns
CUF(1,2.605) = 4.0759p = 0.0435p = 0.136F(2.605,36.68) = 4.934p = 0.0033p = 0.0066 F(2.605,36.68) = 0.2756p = 0.815n/a
Muricholic_alphaNWF(1,2.72) = 0.188p = 0.665n/aF(2.72,58.508) = 35.817p = 3.20 × 1021p = 1.92 × 1020F(2.72,58.508) = 0.755p = 0.507ns
CUF(1,2.81) = 2.24p = 0.135nsF(2.81,77.408) = 57.30p = 7.157 × 1035p = 1.0024 × 1033F(2.81,77.408) = 0.397p = 0.742n/a
Muricholic_beta (Figure 6A)NWF(1,2.911) = 2.623p = 0.105n/aF(2.911,54.878,) = 36.129p = 1.011 × 1022p = 8.08 × 1022F(2.911,54.878,) = 2.706p = 0.0453p = 0.0453
CUF(1,2.68) = 9.452p = 0.0021p = 0.019F(2.68,78.32) = 81.99p = 1.00 × 1047p = 1.9 × 1046F(2.68,78.32) = 0.272p = 0.823n/a
Conjugated Bile Acids
Glycochenodeoxycholic AcidNWF(1,2.819) = 0.578p = 0.447n/aF(2.819,60.355) = 59.90p = 1.503 × 1036p = 1.95 × 1035F(2.819,60.355) = 1.784p = 0.151ns
CUF(1,2.917) = 2.459p = 0.116nsF(2.917,78.25) = 17.47p = 4.35 × 1011p = 3.48 × 1010F(2.917,78.25) = 0.508p = 0.671n/a
Glycocholic AcidNWF(1,2.627) = 0.146p = 0.701n/aF(2.627,60.03) = 108.142p = 1.084 × 1061p = 1.728 × 1060F(2.627,60.03) = 0.257p = 0.831ns
CUF(1,2.744) = 4.479p = 0.0343p = 0.136F(2.744,75.63) = 41.22p = 1.38 × 1024p = 1.794 × 1023F(2.744,75.63) = 0.109p = 0.274n/a
Glycohyocholic AcidNWF(1,2.913) = 0.092p = 0.762n/aF(2.913,60.943) = 28.238p = 8.523 × 1018p = 4.26 × 1017F(2.913,60.943) = 0.514p = 0.667ns
CUF(1,2.706) = 0.543p = 0.4611nsF(2.706,79.81) = 29.44p = 2.083 × 1017p = 1.872 × 1016F(2.706,79.81) = 0.146p = 0.917n/a
Taurochenodeoxycholic AcidNWF(1,2.57) = 0.453p = 0.501n/aF(2.57,58.52) = 43.784p = 1.290 × 1024p = 1.161 × 1023F(2.57,58.52) = 1.378p = 0.250ns
CUF(1,2.82) = 3.133p = 0.0688nsF(2.82,76.881) = 34.930p = 2.46 × 1021p = 2.706 × 1020F(2.82,76.881) = 0.567p = 0.625n/a
Taurocholic AcidNWF(1,2.83) = 0.417p = 0.518n/aF(2.83,60.603) = 15.834p = 7.724 × 1010p = 2.316 × 109F(2.83,60.603) = 0.743p = 0.519ns
CUF(1,2.773) = 6.388p = 0.0115p = 0.069F(2.773,36.7228) = 1.232p = 0.296nsF(2.773,36.7228) = 0.644p = 0.575n/a
Taurohyocholic AcidNWF(1,2.581) = 2.896p = 0.0878n/aF(2.581,58.664) = 57.948p = 1.301 × 1032p = 1.56 × 1031F(2.581,58.664) = 1.256p = 0.288ns
CUF(1,2.381) = 4.492p = 0.0341nsF(2.381,79.613) = 42.977p = 1.39 × 1022p = 1.668 × 1021F(2.381,79.613) = 1.348p = 0.259n/a
Secondary Bile Acids
Deoxycholic Acid (Figure 6B)NWF(1,2.2670) = 5.557 p = 0.0184n/aF(2.2.267,28.994) = 84.80p = 3.18 × 1042p = 4.77 × 1041F(2.267,28.994) = 2.19p = 0.104ns
CUF(1,2.79) = 12.219p = 0.00047p = 0.005F(2.79,79.83) = 62.44p = 8.56 × 1038p = 1.3696 × 1036 F(2.79,79.83) = 2.188p = 0.0918n/a
Lithocholic Acid (Figure 6C)NWF(1,2.832) = 0.240p = 0.624n/aF(2.832,60.296) = 123.84p = 6.77 × 1076p = 1.2186 × 1074F(2.832,60.296) = 3.374p = 0.0196p = 0.0392
CUF(1,2.90) = 10.84p = 0.0009p = 0.010F(2.90,79.89) = 12.19p = 8.72 × 108p = 3.49 × 107F(2.90,79.89) = 1.476p = 0.220n/a
Ursodeoxycholic Acid (Figure 6D)NWF(1,2.539) = 2.465p = 0.164n/aF(2.539,60.446) = 115.64p = 7.468 × 1064p = 1.2699 × 1062F(2.539,60.446) = 2.228p = 0.0935ns
CUF(1,2.532) = 9.188p = 0.00243p = 0.019F(2.532,78.672) = 4.966p = 0.00349p = 0.0066 F(2.532,78.672) = 1.098p = 0.343n/a
Secondary Conjugated Bile Acids
Glycodeoxycholic Acid (Figure 6E)NWF(1,2.818) = 0.485p = 0.486n/aF(2.818,60.903) = 48.193p = 2.045 × 1029p = 2.255 × 1028F(2.818,60.903) = 5.310p = 0.0015p = 0.0045
CUF(1,2.916) = 5.013p = 0.0252p = 0.126F(2.916,79.05) = 31.25p = 1.064 × 1019p = 1.06 × 1018 F(2.916,79.05) = 0.972p = 0.403n/a
Glycolithocholic AcidNWF(1,2.513) = 1.268p = 0.260n/aF(2.513,59.330) = 72.0p = 1.784 × 1039p = 2.492 × 1038F(2.513,59.330) = 0.534p = 0.627ns
CUF(1,2.513) = 0.009p = 0.923nsF(2.513,78.64) = 76.44p = 6.81 × 1042p = 1.1577 × 1040 F(2.513,78.64) = 0.150p = 0.903n/a
Glycoursodeoxycholic AcidNWF(1,2.488) = 0.898p = 0.343n/aF(2.488,56.613) = 48.827p = 1.26 × 1026p = 1.26 × 1025F(2.488,56.613) = 0.194p = 0.920ns
CUF(1,2.911) = 1.24p = 0.265nsF(2.911,79.71) = 74.60p = 6.64 × 1047p = 1.1952 × 1045 F(2.911,79.71) = 0.214p = 0.881n/a
Taurodeoxycholic AcidNWF(1,2.501) = 0.262p = 0.609n/aF(2.501,59.216) = 40.489p = 3.021 × 1022p = 2.114 × 1021F(2.501,59.216) = 0.835p = 0.456ns
CUF(1,2.768) = 4.484p = 0.0342nsF(2.768,78.280) = 16.797p = 3.07 × 1010p = 1.84 × 109F(2.768,78.280) = 0.511p = 0.659n/a
Taurohyodeoxycholic AcidNWF(1,2.817) = 1.459p = 0.227n/aF(2.817,60.123) = 150.64p = 7.138 × 1092p = 1.4994 × 1090F(2.817,60.123) = 1.111p = 0.341ns
CUF(1,2.785) = 7.212p = 0.00724p = 0.050F(2.785,79.746) = 13.787p = 1.68 × 108p = 8.40 × 108F(2.785,79.746) = 0.473p = 0.687n/a
Taurolithocholic AcidNWF(1,2.895) = 0.001p = 0.974n/aF(2.895,60.534) = 18.164p = 1.887 × 1011p = 7.56 × 1011F(2.895,60.534) = 2.139p = 0.095ns
CUF(1,2.85) = 0.006p = 0.937nsF(2.85,78.972) = 17.603p = 5.79 × 1011p = 4.05 × 1010F(2.85,78.972) = 0.784p = 0.497n/a
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Thompson, R.S.; Bowers, S.J.; Vargas, F.; Hopkins, S.; Kelley, T.; Gonzalez, A.; Lowry, C.A.; Dorrestein, P.C.; Vitaterna, M.H.; Turek, F.W.; et al. A Prebiotic Diet Containing Galactooligosaccharides and Polydextrose Produces Dynamic and Reproducible Changes in the Gut Microbial Ecosystem in Male Rats. Nutrients 2024, 16, 1790. https://doi.org/10.3390/nu16111790

AMA Style

Thompson RS, Bowers SJ, Vargas F, Hopkins S, Kelley T, Gonzalez A, Lowry CA, Dorrestein PC, Vitaterna MH, Turek FW, et al. A Prebiotic Diet Containing Galactooligosaccharides and Polydextrose Produces Dynamic and Reproducible Changes in the Gut Microbial Ecosystem in Male Rats. Nutrients. 2024; 16(11):1790. https://doi.org/10.3390/nu16111790

Chicago/Turabian Style

Thompson, Robert S., Samuel J. Bowers, Fernando Vargas, Shelby Hopkins, Tel Kelley, Antonio Gonzalez, Christopher A. Lowry, Pieter C. Dorrestein, Martha Hotz Vitaterna, Fred W. Turek, and et al. 2024. "A Prebiotic Diet Containing Galactooligosaccharides and Polydextrose Produces Dynamic and Reproducible Changes in the Gut Microbial Ecosystem in Male Rats" Nutrients 16, no. 11: 1790. https://doi.org/10.3390/nu16111790

APA Style

Thompson, R. S., Bowers, S. J., Vargas, F., Hopkins, S., Kelley, T., Gonzalez, A., Lowry, C. A., Dorrestein, P. C., Vitaterna, M. H., Turek, F. W., Knight, R., Wright, K. P., Jr., & Fleshner, M. (2024). A Prebiotic Diet Containing Galactooligosaccharides and Polydextrose Produces Dynamic and Reproducible Changes in the Gut Microbial Ecosystem in Male Rats. Nutrients, 16(11), 1790. https://doi.org/10.3390/nu16111790

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