==== Front Neuronal Signal Neuronal Signal ns Neuronal Signaling 2059-6553 Portland Press Ltd. 10.1042/NS20200007 NS20200007 Gastrointestinal, Renal & Hepatic Systems Host-Microbe Interactions Neuroscience Research Articles Lactobacillus rhamnosus GG soluble mediators ameliorate early life stress-induced visceral hypersensitivity and changes in spinal cord gene expression McVey Neufeld Karen-Anne InvestigationMethodologyWriting - review & editingWriting - original draft Strain Conall R. MethodologyWriting - review & editing Pusceddu Matteo M. Methodology Waworuntu Rosaline V. ConceptualizationWriting - review & editing Manurung Sarmauli ConceptualizationWriting - review & editing Gross Gabriele ConceptualizationWriting - review & editing M. Moloney Gerry Methodology Hoban Alan E. Methodology Murphy Kiera MethodologyWriting - review & editing Stanton Catherine Methodology Dinan Timothy G. ConceptualizationFunding acquisitionWriting - review & editing Cryan John F. ConceptualizationFunding acquisitionWriting - review & editing O'Mahony Siobhain M. somahony@ucc.ieConceptualizationResourcesSupervisionFunding acquisitionWriting - original draftWriting - review & editing 1 APC Microbiome Ireland, University College Cork, Cork, Ireland 2 Teagasc Food Research Centre, Moorepark, Co. Cork 3 Mead Johnson Pediatric Nutrition Institute, Evansville, IN, USA 4 Mead Johnson Pediatric Nutrition Institute, Nijmegen, The Netherlands 5 Department of Psychiatry and Neurobehavioural Science, University College Cork, Cork, Ireland 6 Department of Anatomy and Neuroscience, University College Cork, Cork, Ireland Correspondence: Siobhain M. O’Mahony (somahony@ucc.ie) 12 2020 23 11 2020 4 4 NS2020000716 4 2020 26 10 2020 10 11 2020 11 11 2020 © 2020 The Author(s).2020This is an open access article published by Portland Press Limited on behalf of the Biochemical Society and distributed under the .Abstract Visceral hypersensitivity is a hallmark of many functional and stress-related gastrointestinal disorders, and there is growing evidence that the gut microbiota may play a role in its pathophysiology. It has previously been shown that early life stress-induced visceral sensitivity is reduced by various probiotic strains of bacteria (including Lactobacillus rhamnosus GG (LGG)) alone or in combination with prebiotic fibres in rat models. However, the exact mechanisms underpinning such effects remain unresolved. Here, we investigated if soluble mediators derived from LGG can mimic the bacteria’s effects on visceral hypersensitivity and the microbiota–gut–brain axis. Rats were exposed to maternal separation (MS) from postnatal days 2–12. From weaning onwards both non-separated (NS) and MS offspring were provided drinking water with or without supplementation of standardized preparations of the LGG soluble mediators (LSM). Our results show that MS led to increased visceral sensitivity and exaggerated corticosterone plasma levels following restraint stress in adulthood, and both of these effects were ameliorated through LSM supplementation. Differential regulation of various genes in the spinal cord of MS versus NS rats was observed, 41 of which were reversed by LSM supplementation. At the microbiota composition level MS led to changes in beta diversity and abundance of specific bacteria including parabacteroides, which were ameliorated by LSM. These findings support probiotic soluble mediators as potential interventions in the reduction of symptoms of visceral hypersensitivity. Early life stressLactobacillusMicrobiome-Gut-Brain axisProbioticsVisceral pain ==== Body Introduction The occurrence of stress during critical windows of development has extensive and long-term effects on behavior and physiology, which have been observed both preclinically and clinically [1–6]. Observed alterations include, but are not limited to, changes in behavior, brain and gut function and the intestinal microbial profile [7–10]. Early life stressful events are considered a risk factor for a number of psychiatric and gastrointestinal disorders, some of which are highly comorbid [11], indicating a role for the gut–brain axis in the pathogenesis of these disorders. Maternal separation (MS) of rat pups is a robust and reliable model of early life stress that induces long-term alterations to behavior, brain neurochemistry and gut physiology [12]. Affected adult offspring demonstrate stress hyperresponsivity, depressive-like behavior, increased visceral sensitivity, intestinal epithelial barrier dysfunction and reduced gut microbial diversity [12,13], and probiotic feedings have been shown to ameliorate some of these stress-induced alterations [14–17]. We have previously shown that the probiotic Lactobacillus rhamnosus GG (LGG) has potential beneficial effects in this model of early life stress when administered in the drinking water [18]. MS-induced visceral hypersensitivity to colorectal distension (CRD) is a particularly robust outcome that is sensitive to reversal by a variety of probiotic strains (including LGG) alone or in combination with prebiotic fibres in rat models [10,19–23]. In a different animal model, nociceptive changes induced by chronic stress were partially explained by neuroinflammation in the spinal cord [24]. Specifically, hypersensitivity related to the colorectal region of the gastrointestinal tract is associated with changes in gene expression in the lumbosacral segment of the spinal cord [25,26]. This is due to the fact that the colorectum receives predominant innervation from the lumbosacral spinal cord, via the sacral pelvic nerves [27]. There is a growing emphasis on potential novel probiotic-related therapeutic approaches across different fields of biomedicine [28]. Is it possible that products of probiotic metabolism are involved in biological effects, and therefore could deliver benefits without requiring the presence of the live organism per se. Interestingly, supernatants derived from probiotic LGG have been shown to exert a number of effects including immunomodulatory activity, and feeding neonatal mice these soluble mediators ameliorates inflammatory responses in a model of allergic airway inflammation [29]. Additionally, in a model of short-bowel syndrome, supplementing rats with LGG soluble mediators improved intestinal barrier disruption and inflammation, and in fact was more effective than viable LGG supplementation [30]. While a number of studies have shown effects of soluble mediators derived from LGG and other bacterial strains [31–36] in other cases only live bacterial cells exerted bioactivity on host physiology compared with supernatant [37]. Thus, there is a current interest in better understanding mechanisms of action underlying beneficial effects of soluble mediators of probiotics on gut brain axis function. To this end, we assessed the potential long-term benefits of feeding a specific preparation of soluble mediators derived from the probiotic LGG to rats that had been previously subjected to early life MS in terms of visceral pain sensitivity, and whether this observation could partially be explained by spinal gene expression and modulation of the gut microbiota profile. Experimental Animals Fifteen pairs of male and female Sprague-Dawley rats were purchased from Harlan, U.K. and habituated to the Preclinical Research Facility, University College Cork for 1 week. Rats were mated and all 15 females became pregnant. Animals were kept under controlled conditions (21 ± 1°C) on a 12-h light/dark cycle (lights on 7:00 am) and fed ad libitum. Dams with litters and later weaned offspring were maintained in animal cages, dimensions (42 × 25 × 13 cm). All experiments were conducted in accordance with the European Directive 2010/63/EEC, the Recommendation of the S.I No 543 of 2012, and approved by the Animal Experimentation Ethics Committee of the University College Cork. Please see timing of each experimental procedure on Supplementary Figure S1. Maternal separation Day of birth was designated as postnatal day (PND) 0. MS protocol was adapted from that described in [38]. Briefly, at PND0 litters were randomly assigned to undergo MS or to remain as non-separated (NS) controls. At PND2 the litters assigned to MS were removed from the main colony room to an adjacent room where the entire litter was placed into a separate, smaller cage maintained at 30–33°C. Mothers were returned to the home cage and the main colony room so as to prevent communication via ultrasonic vocalization. After a 3 h separation the mothers were again brought into the adjacent room and reunited with their respective litters. Control litters were left undisturbed in the mothers’ home cage except for routine cage cleaning performed once a week. This procedure was repeated from PND2 to PND12 inclusive. The period of separation was carried out between 9:00 am and 12:00 pm daily. Dietary intervention and preparations At PND21 pups were weaned, females were culled and only male rats were used for the remainder of the experiment. Male offspring were randomly group-housed by intervention (three per cage) with no littermates housed together. NS and MS male offspring were randomized into different groups and provided drinking water with or without supplementation of LGG soluble mediators (LSM) from 3 weeks of age until approximately 13 weeks old (end of the experiment). This yielded four separate groups. LSM was produced as previously described (Wu et al., 2018): LGG was cultured under standardized conditions, and spent medium was collected at a defined time point during the late exponential growth phase. This supernatant was desalted, sterile filtered and lyophilized. This specific preparation protocol intended to remove lactic acid that is a nonspecific product of lactic acid bacteria fermentation, as well as any bacterial cells. Reconstituted material was daily freshly prepared in de-ionised drinking water with mixing of the drinking bottle 8 to 10 h after first placing it in the cage. The preparation was at a dose equivalent to 1.3–4.3 × 108 CFU viable LGG/animal/day) of the original LGG culture. All animals received a modified American Institute of Nutrition 93-G control diet (Envigo, WI U.S.A.). Interventions continued throughout behavioural testing. Body weight was measured weekly, water intake daily, and food intake 3 times/week. Twelve animals/group were assessed throughout. Colorectal distension (CRD) CRD occurred when the animals were 11 weeks of age and was carried out during the light cycle of the day. This test was conducted as previously described [39]. Animals were fasted for 24 h and then lightly anaesthetised using isofluorane. A 6 cm polyethylene balloon with a connecting catheter was inserted in the distal colon 1 cm proximal to the anus and fixed to the tail with tape to avoid displacement. Animals were allowed to recover from anaesthesia for at least 15 min before beginning colorectal distension. A customized barostat was used to control balloon inflation and pressure during colorectal distension. The distension paradigm used was an ascending phasic distension from 0 to 80 mmHg. A trained and blinded observer scored each rat for the threshold, or pressure when the first pain behaviour was noted as well as the number of pain behaviours exhibited by each animal. Determination of plasma corticosterone At 12 weeks of age and one week following CRD, animals were removed from home cages at 0830 h and placed in clear Perspex rat restrainers with tails exposed. The tip of each rat’s tail was then cut at an angle with a scalpel blade and blood collected in a 1.5 ml collection tube (approximately 200 μls/time point). The initial blood collection was marked T0 and corresponded to a baseline collection. Animals were then singly housed in restrainers for 30 min. Following the stressor, the animal was removed from the restrainer and a T30 blood collection occurred. The animal was placed back into the home cage after a brief recovery period. Blood was stored on wet ice until spun in a centrifuge and the plasma collected. Plasma was stored in −80°C freezer until an ELISA assay for corticosterone was performed. Tissue extraction and RNA sequencing At 13 weeks of age all animals were culled and samples collected. RNA from rat lumbosacral spinal cord was extracted using commercially available mirVana™ Total Rna Extraction Kit (Ambion/Life Technologies) and DNase treated (Turbo DNA-free, Ambion/Life Technologies) according to manufacturer’s protocol. RNA was high quality with an RNA integrity number value for all samples above 8 (Bioanalyzer, Agilent). Spinal cord samples were then randomly pooled within each group by combining equal amounts of RNA from two to three animals resulting in a final four samples per group sent for RNA sequencing at GATC (Konstanz, Germany), where library preparation, paired-end sequencing (2 × 125 bp) and Fastq-file generation was conducted on an Illumina HiSeq sequencer. Differential gene expression and GO-term enrichment analyses Ensembl release 78 was used for read counting for each gene using HTSeq-Count (v0.6.0) with the following non-default parameters: -s: no; -r: pos; -q –f bam –m intersection-nonempty. Differential gene expression was determined using the DESeq2 R-package (v1.6.2) with default parameters on pairwise comparisons of all possible group combinations. An adjusted P-value ≤ 0.1 (Benjamini–Hochberg method) was considered significantly differentially regulated [40–42]. Differentially expressed genes were then analysed for enrichment of GO-Terms using the DAVID Bioinformatic Resources (v6.8). Microbiota analysis Fresh caecal samples were collected from individual rats at the end of the study. At least 20 mg of fresh caecal material were placed in a microcentrifuge tube and kept on ice until storage at −80°C. DNA was extracted using the DNA Fast Stool DNA Extraction Kit (Qiagen) using the protocol for Gram positive bacteria and including an additional bead beating step at the beginning of the procedure. The microbiota composition of the samples was established by amplicon sequencing of a ∼380 bp fragment of the V3-V4 hypervariable region of the bacterial 16S rRNA gene as outlined in the Illumina 16S Metagenomic Sequencing Library preparation guide (Illumina). Amplicons were quantified, normalised and pooled using the Qubit® dsDNA HS Assay Kit (Life Technologies). Samples were sequenced on the MiSeq sequencing platform at Clinical Microbiomics, using a 2 × 300 cycle kit, following standard Illumina guidelines. Microbiota bioinformatics Forward and reverse reads were merged using USEARCH requiring a merged minimum length of 300 bp, maximum 600 bp and minimum 100 bp overlap. A maximum of one expected error was allowed per read. Primers were trimmed, discarding any sequences that do not have a perfect match to both forward and reverse primers. Sequences were dereplicated, discarding clusters smaller than 5 and clustered at 97% sequence similarity using USEARCH [43]. Suspected chimeras were removed using UCHIME [44]. Taxonomy was assigned to each operational taxonomic unit (OTU) using the RDP Classifier [45]. To estimate alpha diversity, the data set was rarefied to the same number of sequences and the Shannon diversity index was calculated using the diversity function from the vegan package in R [46]. Beta diversity was also calculated in R using Bray–Curtis distances and Principal Coordinate Analysis (PCoA) was performed using the vegan package. Statistical analysis Data was analysed using a repeated measures two-way ANOVA, a two-way ANOVA and LSD’s post-hoc test where appropriate using the statistical software package SPSS 22.0 (IBM). A P-value of 0.05 was selected as the threshold of statistical significance. The rank Kruskal–Wallis test was used for comparison of taxa followed by Dunn’s test with false discovery rate P-value adjustment for multiple comparisons. The Mann U Whitney test was used to assess for differences between the various groups. The Benjamini–Hochberg for false discovery rate was used to calculate q (Q was set at 10%). A two-way ANOVA was carried out on alpha diversity indices for intervention and separation and the interaction between intervention and separation and the independent t test was employed for pairwise comparisons. Beta diversity was compared between groups using the adonis function in the R vegan package at the genus level, using Bray–Curtis as distance measurement. Results Effect of LSM supplementation on animal weights and food and water intake A repeated measures two-way ANOVA showed that all rats gained weight over time (F(8,352) = 4960, P<0.0001) and post hoc analyses showed no specific differences between groups at any timepoint (Supplementary Figure S2A). A repeated measures two-way ANOVA showed that all animals increased their food intake over time (F(6,72) = 471.8, P<0.0001) and a main effect for intervention on food intake (F(3,12) = 5.38, P=0.01), but again post hoc analyses showed no specific differences between groups in terms of amount of food ingested at a particular timepoint (Supplementary Figure S2B). There was a main effect for fluid intake over time (F(7,84) = 131, P<0.0001), but no main effect for LSM on fluid intake (F(3,12) = 2.73, P=0.09) and no specific differences between groups in terms of fluid intake at a particular timepoint (Supplementary Figure S2C). Effect of feeding LSM on visceral sensitivity A significant main effect for the intervention was observed in the CRD test for visceral sensitivity (F(1,44) = 7.1, P=0.01), but no main effect for early life stress (F(1,44) = 2.66, P=0.11). There was no significant interaction (F(1,44) = 0.8267, P=0.3683). Fisher’s LSD post hoc showed a specific difference between MS control rats versus NS control (P=0.04) with MS rats showing lower pain thresholds. There was also a significant difference between MS LSM rats versus MS controls (P=0.01) in that LSM feeding resulted in higher pain thresholds in the early life stressed groups (see Figure 1). There was no difference between MS control and MS rats fed bacterial free unconditioned medium (data not shown) nor was there any differences noted for the number of pain behaviors. Figure 1 The threshold of visceral sensitivity in colorectal distension Early life stress caused lower pain thresholds as measured by CRD in MS versus NS control rats that was ameliorated by supplementation with LSM. *P=0.04 versus NS control; #P=0.001 versus MS control. Data are represented as mean ± SEM, n=12/grp. Effects of feeding LSM on stress responsivity Measures of plasma corticosterone were made for both the T0 (baseline) and T30 (stressed) timepoints and a Δcorticosterone value was calculated from these for each animal. A significant early life and intervention interaction was observed (F(1,41) = 8.49, P=0.006), but no main effect for either early life stress (F(1,41) = 0.006, P=0.93) or supplementation (F(1,41) = 1.32, P=0.26). Fisher’s LSD post hoc showed a significantly higher Δcorticosterone for MS versus NS controls (P=0.04) that was attenuated after LSM intervention (P=0.007) (see Figure 2). Figure 2 Stress reactivity as measured by plasma corticosterone Early life stress caused a significant increase in plasma corticosterone values following a 30 min restraint stress, which was attenuated by LSM supplementation. *P=0.04 versus NS control; #P=0.007 versus MS control. Data represented as mean ± SEM, n=12/grp. Effects of feeding LSM on mRNA expression in the spinal cord Unbiased deep sequencing of mRNA was carried out on spinal cord tissue. Analysis of gene expression by pairwise comparisons revealed a number of differentially regulated genes across experimental groups (P<0.1). First, when comparing mRNA expression profile of NS and MS control animals a total of 86 differentially expressed genes were identified (23 down-regulated and 63 up-regulated) (Figure 3A). Next, we investigated the impact of LSM supplementation on spinal cord gene expression in MS rats compared with MS control animals. A total of 171 genes were identified (119 down-regulated and 52 up-regulated). Interestingly, we found 41 genes that were up-regulated in MS control animals compared with NS, and were down-regulated after administration of LSM (Figure 3A and Table 1). We also determined that amongst the total down-regulated genes (78+41 = 119) in MS animals after LSM supplementation, there was a significant enrichment in Wnt family of signaling peptides which are key players in pain (Dkk2, Tcf7, Nxn, Lef1, Nfatc4, Cdh1, Col1a1, Apcdd1 and Cpz) (Figure 3A). Within this subset of overlapping genes, we investigated biological function linked to these genes. GO analysis showed significant enrichment for biological processes specifically associated with stress (Figure 3B). All gene changes are now included on volcano plots with the horizontal dashed line being fold change and vertical dashed line is 0.1 which is the cut off for a P value when it is logged (Supplementary Figure S3A,B). Figure 3 mRNA expression in spinal cord and functional enrichment analysis of the overlapping genes showed that 12 of these were associated with pain, stress or both (A) Numbers of gene differentially expressed between groups including GO enrichment of overlapping genes particularly associated with the Wnt signaling. (B) Biological processes linked to differentially expressed genes. Of the 63 spinal cord genes that were up-regulated by MS, 41 of these were down-regulated with LSM administration. Spinal cord samples from three animals within each group were randomly pooled resulting in a final four samples per group sent for RNA sequencing. Table 1 Spinal cord genes noted to be up-regulated due to maternal separation and down-regulated by the administration of LSM Description Associated gene name NS versus MS Log2 (fold change) MS versus MS + LSM Log2 (fold change) Annexin A1 Anxa1 −0.89775 −0.822068 Annexin A8 Anxa8 −1.23281 −2.17727 Asialoglycoprotein receptor 2 Asgr2 −0.51856 −0.481867 Biglycan Bgn −0.26825 −0.419479 Complement component 7 C7 −0.46183 −0.467474 Cyclin F Ccnf −0.69414 −0.638514 Cholinergic receptor nicotinic beta 3 subunit Chrnb3 −0.81698 −0.819548 Collagen, type XVIII, alpha 1 Col18a1 −0.38027 −0.549118 Collagen, type I, alpha 1 Col1a1 −0.70161 −0.904032 Collagen, type VI, alpha 1 Col6a1 −0.40196 −0.590088 Collagen, type VI, alpha 2 Col6a2 −0.43194 −0.516948 Carboxypeptidase X (M14 family), member 2 Cpxm2 −0.38027 −0.907097 Carboxypeptidase Z Cpz −0.9267 −0.823511 Cytokine like 1 Cytl1 −0.73819 −1.21628 Elastin Eln −0.78315 −0.946259 Endothelial cell-specific molecule 1 Esm1 −1.14174 −1.53702 Filamin A Flna −0.29314 −0.716614 Glutathione S-Transferase M2 Gstm2 −0.6405 −0.43816 Hydroxycarboxylic acid receptor 1 Hcar1 −0.8035 −0.66941 Interferon induced transmembrane protein 2 Ifitm2 −0.56305 −0.39145 Insulin-like growth factor binding protein 3 Igfbp3 −0.80728 −0.78094 Insulin-like growth factor binding protein 4 Igfbp4 −0.47974 −0.50599 Keratin 10 Krt10 −0.94947 −0.43816 Keratin 13 Krt13 −3.88679 −2.58248 Lectin, galactoside-binding, soluble, 7 Lgals7 −2.47941 −2.52878 Leukocyte immunoglobulin-like receptor, subfamily B (with TM and ITIM domains), member 3-like Lilrb3l −1.26735 −0.87364 Lysyl oxidase Lox −0.47596 −1.04757 Microtubule-associated protein, RP/EB family, member 2 Mapre2 −0.60424 −0.80124 Mannose receptor, C type 1 Mrc1 −0.46494 −0.56762 Myosin, heavy chain 11, smooth muscle Myh11 −0.46892 −0.33294 Pappalysin 2 Pappa2 −0.76301 −1.38181 Procollagen C-endopeptidase enhancer Pcolce −0.63158 −1.34621 Periostin Postn −0.87747 −0.41055 Polymerase I and transcript release factor Ptrf −0.31889 −0.36173 RT1 class II, locus Bb RT1-Bb −0.43733 Serpin family F member 1 Serpinf1 −0.71682 −0.52703 Thrombospondin 2 Thbs2 −0.29513 −0.29096 Effects of feeding LSM on caecal microbiota Alpha diversity Alpha diversity of the samples was measured by observed species and Shannon diversity. The observed species index measures the number of different bacterial species per sample that is defined as ‘richness’. However, regarding diversity, not only the qualitative amount of species, but also the abundance of the species must be taken into account. The relative abundances of the different species making up the samples’ richness are defined as ‘evenness’. The Shannon-diversity index relates both, OTU richness and evenness. A two-way ANOVA for LSM and stress and the interaction between these on alpha diversity indices was carried out. There was a significant effect of LSM on number of observed species F(1,43) = 5.57, P=0.02 (Figure 4A), with no significant interactions between stress and LSM. There was a significant interaction between stress and LSM for Shannon F(1,43) = 6.97, P=0.02. There was no significant effect of stress on the alpha diversity indexes (Figure 4B). Figure 4 Alpha diversity of the caecal microbiome (A) Observed species/ richness was increased on the MS LSM group. (B) Shannon index/species diversity was increased in the MS LSM group also. *P<0.05 versus MS water group; **P<0.01 versus MS water group. Taxonomic analysis NS control versus MS control The effect of MS was assessed at the three taxonomic levels by comparing NS unsupplemented control rats with MS control rats. At the phylum level, there was a decrease in the relative abundance of Proteobacteria in MS rats (P=0.04; q=0.01), while at the family level, Cryomorphaceae (P=0.02; q=0.003) abundances were significantly lower and Prevotellaceae (P=0.02; q=0.005) significantly higher in MS animals compared with NS. MS significantly lowered the relative abundance of 6 genera (Clostridium cluster XIVa (P=0.002; q=0.001); Bilophia (P=0.006; q=0.002); Wandonia (P=0.02; q=0.004); Acetivibrio (P=0.03; q=0.005); Parabacteroides (P=0.04; q=0.006) and Anaerosinus (P=0.04; q=0.007) and increased Alloprevotella (P=0.02; q=0.003) populations (see Table 1 and Supplementary Figure S4). NS control versus NS LSM LSM supplemented NS animals had significantly lower Proteobacteria (P=0.01; q=0.01) and Desulfovibrionaceae (P=0.03; q=0.005) relative abundances compared with NS control rats. While relative abundance of Ruminococcaceae was significantly increased in LSM NS animals (P=0.002; q=0.003) as was Faecalibacterium (P=0.03; q=0.002). Acetitomaculum (P=0.03; q=0.001) relative abundances was noted to be lower than in NS control rats (see Table 1 and Supplementary Figure S4). MS control versus MS LSM In MS rats, LSM treated animals had significantly greater Peptostreptococcaceae (P=0.006; q=0.006), Barnesiella (P=0.03; q=0.006), Parabacteroides (P=0.019; q=0.005), Oribacterium (P=0.02; q=0.004) and Clostridium cluster XI (P=0.006; q=0.003) relative abundances compared with MS control rats. In contrast, relative abundance of Sutterellaceae (P=0.004; q=0.003), Tannerella (P<0.001; q=0.001) and Parasutterella (P=0.004; q=0.002) in LSM MS rats was lower compared with control MS animals (see Table 2 and Supplementary Figure S4). Table 2 The significant differences seen in the caecal microbiota of all groups at phyla, family and genus level NS Control versus MS Control Taxa NS Control MS Control P value q value Phylum Proteobacteria 2.49% 1.83% 0.039 0.0125 Family Cryomorphaceae 0.04% 0.01% 0.02 0.0026 Prevotellaceae 0.97% 2.22% 0.024 0.0053 Genus Clostridium cluster XIVa 7.66% 3.60% 0.002 0.0010 Bilophia 0.22% 0.07% 0.006 0.0019 Alloprevotella 0.93% 2.21% 0.02 0.0029 Wandonia 0.04% 0.01% 0.02 0.0039 Acetivibrio 0.03% 0.11% 0.028 0.0049 Parabacteroides 0.58% 0.44% 0.039 0.0058 Anaerosinus 0.04% 0.02% 0.039 0.0068 NS control versus NS LSM Taxa NS Control NS LSM P value q value Phylum Proteobacteria 2.49% 1.78% 0.017 0.0125 Family Ruminococcaceae 22.80% 31.29% 0.002 0.0026 Desulfovibrionaceae 1.95% 1.34% 0.033 0.0053 Genus Faecalibacterium 3.95% 9.13% 0.028 0.0019 Acetitomaculum 0.03% 0.21% 0.033 0.0010 MS Control versus MS LSM Taxa MS Water MS LSM P value q value Family Sutterellaceae 0.09% 0.01% 0.004 0.0026 Peptostreptococcaceae 1.84% 3.37% 0.006 0.0053 Genus Tannerella 18.50% 9.32% 0.001 0.0010 Parasutterella 0.09% 0.01% 0.004 0.0019 Clostridium cluster XI 1.84% 3.37% 0.006 0.0029 Oribacterium 0.18% 0.67% 0.016 0.0039 Parabacteroides 0.44% 0.62% 0.019 0.0049 Barnesiella 19.53% 24.02% 0.032 0.0058 % refers to percentage of relative abundance of bacteria. Q value is the adjusted P value for false discovery rates Beta diversity of caecal microbiome Bray–Curtis plots were generated from genus files. The Adonis Permanova test was carried out to assess whether caecal microbiome had significantly dissimilar beta diversity due to early life stress or administration of LSM. Early life stress significantly impacted on the beta diversity, seen as differences between the MS control group receiving water only and the NS control group (P<0.05). There was no significant difference between the NS groups receiving LSM or water (P=0.18). In contrast, MS LSM rats had a significantly dissimilar caecal microbiome to control counterparts (P=0.001) (see Figure 5). Figure 5 Beta diversity of the caecal microbiome (A) Bray–Curtis plot of NS Water and MS Water rats caecal microbiomes shows significant impact of early life stress on beta diversity. (B) NS LSM and NS Water rats caecal microbiomes with no differences in beta diversity observed (overlapping area). (C) Bray–Curtis plot of MS LSM and MS Water rats caecal microbiomes shows significant impact of LSM on beta diversity in the stressed group. Discussion There is a growing interest on the potential of probiotic strains to ameliorate the effects of stress. Within the mechanisms underpinning such effects, the production of neuroactive metabolites is emerging as a key player [47]. The present study implicates soluble mediators derived from the probiotic LGG to exert beneficial effects on stress induced visceral hypersensitivity. In the present study, as previously shown early life stress led to an increased sensitivity to colorectal distension [26,48–51]. Interestingly, we noted that feeding with LSM from weaning onwards resulted in an attenuation of this visceral hypersensitivity, indicating that this probiotic-derived intervention ameliorates early-life stress-induced dysfunction. Given the co-morbidity between gastrointestinal and stress-related disorders we also investigated the impact of LSM on hypothalamic–pituitary–adrenal (HPA) axis responsivity by measuring the response to restraint stress. We noted that an increased HPA axis activity induced by early life stress was also amenable to reversal by LSM. We have previously reported that early life MS in rats results in specific changes in gene expression in the lumbosacral spinal cord that were associated with increased visceral pain sensitivity [25,26]. Here using RNA sequencing we expand on our previous data to demonstrate that the expression levels of a total of 86 genes differed between the NS and MS groups. Of interest, 63 of these genes were up-regulated as a result of early life MS and from this subset, expression levels were reversed in 41 of these after LSM administration in MS animals. These 41 overlapping genes were further examined by a functional enrichment analysis and 12 of these were found to be associated with pain circuitry, stress, or both. These findings of changes in spinal cord gene expression are thus consistent with both the CRD and HPA results we observed. One of the genes found to be up-regulated by MS and then down-regulated following LSM intervention coded for the protein annexin A1. Annexin A1 is found in abundance in cells involved in neuroendocrine communication, and its synthesis and release is dependent on glucocorticoids. It is also thought to play a significant role in the perinatal programming of the HPA axis [52]. Of the total down-regulated genes seen in MS animals after LSM supplementation, there was a significant involvement in Wnt signaling. Previous studies has indicated a role for WNT signaling in a rat model of neuropathic pain where a rapid-onset and long-lasting expression of WNTs in the primary sensory neurons, the spinal dorsal horn neurons, and astrocytes were noted [53]. Moreover, spinal blockade of WNT signaling pathways prevented the induction of neuropathic pain, the accompanying neurochemical changes and inflammation without affecting normal pain sensitivity and locomotor activity [53]. While further research is needed to better understand the mechanisms regulating the molecular interactions among stress, visceral sensitivity and Wnt signaling, these pathways with regard to the brain–gut axis communication may provide a novel mechanistic insight and effective approach into the possibility to ameliorate the impact of early life stress. At the microbiota level, alpha diversity of the caecal samples was measured by observed species and Shannon diversity. While early life stress did not induce changes in alpha diversity per se, LSM increased both richness and evenness of the microbiota in the rats exposed to the early life stress perhaps indicating that the microbiota population in the MS rats was more amenable to change. Beta diversity was assessed using Bray–Curtis dissimilarity based on abundance level and we noted that the NS and MS control groups had different species abundances. Furthermore, we also noted that the MS group receiving LSM group when compared with its control had a different beta diversity. This highlights the fact that both early life stress and LSM induce significant changes to the microbiota at a species abundance level. Furthermore, the impact of LSM on specific caecal microbial populations was not consistent between MS and NS animals suggesting that the effect of LSM on the microbiome could be through altering signaling rather than a direct metabolic/antimicrobial effect. MS altered the microbiome at various taxonomic levels, notably significantly reduced clostridium cluster XIVa relative abundances at the genus level. The clostridium cluster XIVa contains many butyrate producing species and are considered important microbial mediators in gut health and homeostasis [54]. Many species belonging to this cluster have putative anti-inflammatory activities including induction of colonic regulatory T cells [55]. Furthermore, clostridium cluster XIVa have been shown to be reduced in several disease states such as Crohns disease [56] and cystic fibrosis [57]. Drilling down into the microbiota data confirms that LSM administration led to a significantly different microbiota profile yet was unable to reverse the effects of early life stress per se. LSM had a differential impact on the gut microbiota of the NS and MS rats. More specifically in the MS rats, LSM increased the abundance of Barnesiella species. Barnesiella is a genus that has been shown to protect the gut against infection in rats and may be a key regulator of gut homeostasis [58]. LSM also significantly reduced the abundance of Tannerella species. Tannerella species are considered a major oral pathogen [59] and its presence has been linked to increased risk of esophageal cancer [60]. Of particular note is the decrease of the genus Parabacteroides in the MS compared with the NS control group that was reversed by LSM supplementation. Some species within this genus, such as Parabacteroides distasonis has been shown to significantly decrease the severity of gastrointestinal inflammation in murine models of acute and chronic colitis induced by dextran sulphate sodium [61] and also to be negatively correlated with multiple sclerosis, non-alcoholic fatty liver disease and obesity [62,63]. Furthermore, in vitro cultivation of Parabacteroides distasonis demonstrated its capacity to transform bile acids and production of succinate [64]. While more mechanistic investigations are necessary, it is possible that changes in the abundance of Parabacteroides may play a role in the beneficial impact of dietary interventions on behavior. This is in line with our previous findings related to the effects of a blend of prebiotics (polydextrose/galacto-oligosaccharide prebiotic blend) positively influencing behaviour and brain gene expression while also increasing the abundance of parabacteroides [65]. The significant and differential increase in Faecalibacterium genera observed NS animals supplemented with LSM compared with NS control animals are also notable as this genera is attributed to have anti-inflammatory properties [66] including antinociceptive effects in anti-inflammatory disease models [37] and has been shown to be protective against pathological gut inflammation in a murine model [67]. LSM was also found to increase Ruminococcaceae populations, another family noted for butyrate production and therefore may play a role in promoting gut health [68]. Taken together, these results suggest that LSM ameliorates stress induced visceral hypersensitivity, can promote gut microbiome diversity in rats that have undergone MS stress and can impact gut microbial composition as well as spinal cord gene expression. In this model, the effects of LSM have not been dependent on administration of viable probiotic bacteria but might have been mediated by LGG metabolites and fermentation products present in the preparation such as short-chain fatty acids, bacteriocins, or immunomodulatory components of LGG such as toll-like receptor ligands. Besides potential direct effects on the host, these might have also exerted indirect effects by influencing the abundance of bacterial species in the intestine. Further studies are needed to address the exact mechanisms. Overall, these findings intriguingly support probiotic soluble mediators as potential interventions in the management of stress induced gastrointestinal symptoms or disorders. Supplementary Material Supplementary Figures S1-S4 Click here for additional data file. Data Availability All supporting are included in the main article and its supplementary files. Competing Interests The authors declare that there are no competing interests associated with the manuscript. Funding This work was funded by Mead Johnson Nutrition. CRediT Author Contribution Karen-Anne McVey Neufeld: Investigation, Methodology, Writing - review & editing, Writing - original draft. Conall R. Strain: Methodology, Writing - review & editing. Matteo M. Pusceddu: Methodology. Rosaline V. Waworuntu: Conceptualization, Writing - review & editing. Sarmauli Manurung: Conceptualization, Writing - review & editing. Gabriele Gross: Conceptualization, Writing - review & editing. Gerry Moloney: Methodology. Alan E. Hoban: Methodology. Kiera Murphy: Methodology, Writing - review & editing. Catherine Stanton: Methodology. Timothy G. Dinan: Conceptualization, Funding acquisition, Writing - review & editing. John F. Cryan: Conceptualization, Funding acquisition, Writing - review & editing. Siobhain M. O’Mahony: Conceptualization, Resources, Funding acquisition, Writing - original draft, Writing - review & editing. Abbreviations CRDcolorectal distension HPAhypothalamic–pituitary–adrenal LGGLactobacillus rhamnosus GG LSMLGG soluble mediators MSmaternal separation NSnon-separated PNDpostnatal day ==== Refs References 1. O'Mahony S.M. , Clarke G. , Dinan T.G. and Cryan J.F. (2017 ) Early-life adversity and brain development: Is the microbiome a missing piece of the puzzle? Neuroscience 342 , 37 –54 10.1016/j.neuroscience.2015.09.068 26432952 2. Cryan J.F. and Dinan T.G. (2013 ) Unraveling the longstanding scars of early neurodevelopmental stress . Biol. Psychiatry 74 , 788 –789 10.1016/j.biopsych.2013.10.004 24188696 3. McVey Neufeld K.A. , Luczynski P. , Dinan T.G. and Cryan J.F. (2016 ) Reframing the Teenage Wasteland: Adolescent Microbiota-Gut-Brain Axis . Can. J. Psychiatry 61 , 214 –221 10.1177/0706743716635536 27254413 4. Maccari S. , Krugers H.J. , Morley-Fletcher S. , Szyf M. and Brunton P.J. (2014 ) The consequences of early-life adversity: neurobiological, behavioural and epigenetic adaptations . J. Neuroendocrinol. 26 , 707 –723 10.1111/jne.12175 25039443 5. Kaplow J.B. and Widom C.S. (2007 ) Age of onset of child maltreatment predicts long-term mental health outcomes . J. Abnorm. Psychol. 116 , 176 –187 10.1037/0021-843X.116.1.176 17324028 6. Coplan J.D. , Andrews M.W. , Rosenblum L.A. , Owens M.J. , Friedman S. , Gorman J.M. et al. (1996 ) Persistent elevations of cerebrospinal fluid concentrations of corticotropin-releasing factor in adult nonhuman primates exposed to early-life stressors: implications for the pathophysiology of mood and anxiety disorders . Proc. Natl. Acad. Sci. U.S.A. 93 , 1619 –1623 10.1073/pnas.93.4.1619 8643680 7. Syed S.A. and Nemeroff C.B. (2017 ) Early Life Stress, Mood, and Anxiety Disorders . Chronic Stress (Thousand Oaks) 1 , 10.1177/2470547017694461 28649671 8. Park S.H. , Videlock E.J. , Shih W. , Presson A.P. , Mayer E.A. and Chang L. (2016 ) Adverse childhood experiences are associated with irritable bowel syndrome and gastrointestinal symptom severity . Neurogastroenterol. Motil. 28 , 1252 –1260 10.1111/nmo.12826 27061107 9. Moloney R.D. , O'Leary O.F. , Felice D. , Bettler B. , Dinan T.G. and Cryan J.F. (2012 ) Early-life stress induces visceral hypersensitivity in mice . Neurosci. Lett. 512 , 99 –102 10.1016/j.neulet.2012.01.066 22326388 10. O'Mahony S.M. , Marchesi J.R. , Scully P. , Codling C. , Ceolho A.M. , Quigley E.M. et al. (2009 ) Early life stress alters behavior, immunity, and microbiota in rats: implications for irritable bowel syndrome and psychiatric illnesses . Biol. Psychiatry 65 , 263 –267 10.1016/j.biopsych.2008.06.026 18723164 11. Whitehead W.E. , Palsson O. and Jones K.R. (2002 ) Systematic review of the comorbidity of irritable bowel syndrome with other disorders: what are the causes and implications? Gastroenterology 122 , 1140 –1156 10.1053/gast.2002.32392 11910364 12. O'Mahony S.M. , Hyland N.P. , Dinan T.G. and Cryan J.F. (2011 ) Maternal separation as a model of brain-gut axis dysfunction . Psychopharmacology (Berl.) 214 , 71 –88 10.1007/s00213-010-2010-9 20886335 13. Soderholm J.D. , Yates D.A. , Gareau M.G. , Yang P.C. , MacQueen G. and Perdue M.H. (2002 ) Neonatal maternal separation predisposes adult rats to colonic barrier dysfunction in response to mild stress . Am. J. Physiol. Gastrointest. Liver Physiol. 283 , G1257 –G1263 10.1152/ajpgi.00314.2002 12388189 14. Vanhaecke T. , Aubert P. , Grohard P.A. , Durand T. , Hulin P. , Paul-Gilloteaux P. et al. (2017 ) L. fermentum CECT 5716 prevents stress-induced intestinal barrier dysfunction in newborn rats . Neurogastroenterol. Motil. 29 , 10.1111/nmo.13069 28370715 15. Cowan C.S. , Callaghan B.L. and Richardson R. (2016 ) The effects of a probiotic formulation (Lactobacillus rhamnosus and L. helveticus) on developmental trajectories of emotional learning in stressed infant rats . Transl. Psychiatry 6 , e823 10.1038/tp.2016.94 27244232 16. Desbonnet L. , Garrett L. , Clarke G. , Kiely B. , Cryan J.F. and Dinan T.G. (2010 ) Effects of the probiotic Bifidobacterium infantis in the maternal separation model of depression . Neuroscience 170 , 1179 –1188 10.1016/j.neuroscience.2010.08.005 20696216 17. Gareau M.G. , Jury J. , MacQueen G. , Sherman P.M. and Perdue M.H. (2007 ) Probiotic treatment of rat pups normalises corticosterone release and ameliorates colonic dysfunction induced by maternal separation . Gut 56 , 1522 –1528 10.1136/gut.2006.117176 17339238 18. McVey Neufeld K.A. , O'Mahony S.M. , Hoban A.E. , Waworuntu R.V. , Berg B.M. , Dinan T.G. et al. (2019 ) Neurobehavioural effects of Lactobacillus rhamnosus GG alone and in combination with prebiotics polydextrose and galactooligosaccharide in male rats exposed to early-life stress . Nutr. Neurosci. 22 , 425 –434 10.1080/1028415X.2017.1397875 29173065 19. Ait-Belgnaoui A. , Payard I. , Rolland C. , Harkat C. , Braniste V. , Theodorou V. et al. (2018 ) Bifidobacterium longum and Lactobacillus helveticus Synergistically Suppress Stress-related Visceral Hypersensitivity Through Hypothalamic-Pituitary-Adrenal Axis Modulation . J. Neurogastroenterol. Motil. 24 , 138 –146 10.5056/jnm16167 29291614 20. Kannampalli P. , Pochiraju S. , Chichlowski M. , Berg B.M. , Rudolph C. , Bruckert M. et al. (2014 ) Probiotic Lactobacillus rhamnosus GG (LGG) and prebiotic prevent neonatal inflammation-induced visceral hypersensitivity in adult rats . Neurogastroenterol. Motil. 26 , 1694 –1704 10.1111/nmo.12450 25298006 21. Distrutti E. , Cipriani S. , Mencarelli A. , Renga B. and Fiorucci S. (2013 ) Probiotics VSL#3 protect against development of visceral pain in murine model of irritable bowel syndrome . PLoS ONE 8 , e63893 23691109 22. McKernan D.P. , Fitzgerald P. , Dinan T.G. and Cryan J.F. (2010 ) The probiotic Bifidobacterium infantis 35624 displays visceral antinociceptive effects in the rat . Neurogastroenterol. Motil. 22 , 1029e268 –1035e268 10.1111/j.1365-2982.2010.01520.x 20518856 23. Verdu E.F. , Bercik P. , Verma-Gandhu M. , Huang X.X. , Blennerhassett P. , Jackson W. et al. (2006 ) Specific probiotic therapy attenuates antibiotic induced visceral hypersensitivity in mice . Gut 55 , 182 –190 10.1136/gut.2005.066100 16105890 24. Bradesi S. , Karagiannides I. , Bakirtzi K. , Joshi S.M. , Koukos G. , Iliopoulos D. et al. (2015 ) Identification of Spinal Cord MicroRNA and Gene Signatures in a Model of Chronic Stress-Induced Visceral Hyperalgesia in Rat . PLoS ONE 10 , e0130938 10.1371/journal.pone.0130938 26222740 25. Moloney R.D. , Sajjad J. , Foley T. , Felice V.D. , Dinan T.G. , Cryan J.F. et al. (2016 ) Estrous cycle influences excitatory amino acid transport and visceral pain sensitivity in the rat: effects of early-life stress . Biol. Sex Differ 7 , 33 10.1186/s13293-016-0086-6 27429736 26. Moloney R.D. , Stilling R.M. , Dinan T.G. and Cryan J.F. (2015 ) Early-life stress-induced visceral hypersensitivity and anxiety behavior is reversed by histone deacetylase inhibition . Neurogastroenterol. Motil. 27 , 1831 –1836 10.1111/nmo.12675 26403543 27. Brierley S.M. , Hibberd T.J. and Spencer N.J. (2018 ) Spinal Afferent Innervation of the Colon and Rectum . Front. Cell Neurosci. 12 , 467 10.3389/fncel.2018.00467 30564102 28. Aguilar-Toalá J.E. , Garcia-Varela R. , Garcia H.S. , Mata-Haro V. , González-Córdova A.F. , Vallejo-Cordoba B. et al. (2019 ) Postbiotics: An evolving term within the functional foods field . Trends Food Sci. Technol. 75 , 105 –114 10.1016/j.tifs.2018.03.009 29. Harb H. , van Tol E.A. , Heine H. , Braaksma M. , Gross G. , Overkamp K. et al. (2013 ) Neonatal supplementation of processed supernatant from Lactobacillus rhamnosus GG improves allergic airway inflammation in mice later in life . Clin. Exp. Allergy 43 , 353 –364 10.1111/cea.12047 23414544 30. Wu J. , Yang K. , Wu W. , Tang Q. , Zhong Y. , Gross G. et al. (2018 ) Soluble Mediators From Lactobacillus rhamnosus Gorbach-Goldin Support Intestinal Barrier Function in Rats After Massive Small-Bowel Resection . JPEN J. Parenter. Enteral Nutr. 42 , 1026 –1034 10.1002/jpen.1044 30133842 31. Bermudez-Brito M. , Munoz-Quezada S. , Gomez-Llorente C. , Matencio E. , Bernal M.J. , Romero F. et al. (2013 ) Cell-free culture supernatant of Bifidobacterium breve CNCM I-4035 decreases pro-inflammatory cytokines in human dendritic cells challenged with Salmonella typhi through TLR activation . PLoS ONE 8 , e59370 10.1371/journal.pone.0059370 23555025 32. Ditu L.M. , Chifiriuc M.C. , Bezirtzoglou E. , Marutescu L. , Bleotu C. , Pelinescu D. et al. (2014 ) Immunomodulatory effect of non-viable components of probiotic culture stimulated with heat-inactivated Escherichia coli and Bacillus cereus on holoxenic mice . Microb. Ecol. Health Dis. 25 , 10.3402/mehd.v25.23239 25317114 33. Fong F.L.Y. , Kirjavainen P.V. and El-Nezami H. (2016 ) Immunomodulation of Lactobacillus rhamnosus GG (LGG)-derived soluble factors on antigen-presenting cells of healthy blood donors . Sci. Rep. 6 , 22845 10.1038/srep22845 26961406 34. Gao J. , Li Y. , Wan Y. , Hu T. , Liu L. , Yang S. et al. (2019 ) A Novel Postbiotic From Lactobacillus rhamnosus GG With a Beneficial Effect on Intestinal Barrier Function . Front. Microbiol. 10 , 477 10.3389/fmicb.2019.00477 30923519 35. Ludwig I.S. , Broere F. , Manurung S. , Lambers T.T. , van der Zee R. and van Eden W. (2018 ) Lactobacillus rhamnosus GG-Derived Soluble Mediators Modulate Adaptive Immune Cells . Front. Immunol. 9 , 1546 10.3389/fimmu.2018.01546 30042761 36. Yan F. , Cao H. , Cover T.L. , Whitehead R. , Washington M.K. and Polk D.B. (2007 ) Soluble proteins produced by probiotic bacteria regulate intestinal epithelial cell survival and growth . Gastroenterology 132 , 562 –575 10.1053/j.gastro.2006.11.022 17258729 37. Miquel S. , Martin R. , Lashermes A. , Gillet M. , Meleine M. , Gelot A. et al. (2016 ) Anti-nociceptive effect of Faecalibacterium prausnitzii in non-inflammatory IBS-like models . Sci. Rep. 6 , 19399 10.1038/srep19399 26775847 38. Neumann I.D. , Wigger A. , Kromer S. , Frank E. , Landgraf R. and Bosch O.J. (2005 ) Differential effects of periodic maternal separation on adult stress coping in a rat model of extremes in trait anxiety . Neuroscience 132 , 867 –877 10.1016/j.neuroscience.2005.01.034 15837146 39. O'Mahony S.M. , Bulmer D.C. , Coelho A.M. , Fitzgerald P. , Bongiovanni C. , Lee K. et al. (2010 ) 5-HT(2B) receptors modulate visceral hypersensitivity in a stress-sensitive animal model of brain-gut axis dysfunction . Neurogastroenterol. Motil. 22 , 573e124 –578e124 20003079 40. Hoban A.E. , Stilling R.M. , Ryan F.J. , Shanahan F. , Dinan T.G. , Claesson M.J. et al. (2016 ) Regulation of prefrontal cortex myelination by the microbiota . Transl. Psychiatry 6 , e774 10.1038/tp.2016.42 27045844 41. Hoban A.E. , Stilling R.M. , Moloney G. , Shanahan F. , Dinan T.G. , Clarke G. et al. (2018 ) The microbiome regulates amygdala-dependent fear recall . Mol. Psychiatry 23 , 1134 –1144 10.1038/mp.2017.100 28507320 42. Love M.I. , Huber W. and Anders S. (2014 ) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 . Genome Biol. 15 , 550 10.1186/s13059-014-0550-8 25516281 43. Edgar R.C. (2010 ) Search and clustering orders of magnitude faster than BLAST . Bioinformatics 26 , 2460 –2461 10.1093/bioinformatics/btq461 20709691 44. Edgar R.C. , Haas B.J. , Clemente J.C. , Quince C. and Knight R. (2011 ) UCHIME improves sensitivity and speed of chimera detection . Bioinformatics 27 , 2194 –2200 10.1093/bioinformatics/btr381 21700674 45. Wang Q. , Garrity G.M. , Tiedje J.M. and Cole J.R. (2007 ) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy . Appl. Environ. Microbiol. 73 , 5261 –5267 10.1128/AEM.00062-07 17586664 46. Oksanen J. , Blanchet G. , Kindt R. , Legendre P. , O'Hara R. , Simpson G. et al. (2011 ) vegan: Community Ecology Package 47. Cryan J.F. , O'Riordan K.J. , Cowan C.S.M. , Sandhu K.V. , Bastiaanssen T.F.S. , Boehme M. et al. (2019 ) The Microbiota-Gut-Brain Axis . Physiol. Rev. 99 , 1877 –2013 10.1152/physrev.00018.2018 31460832 48. Hyland N.P. , O'Mahony S.M. , O'Malley D. , O'Mahony C.M. , Dinan T.G. and Cryan J.F. (2015 ) Early-life stress selectively affects gastrointestinal but not behavioral responses in a genetic model of brain-gut axis dysfunction . Neurogastroenterol. Motil. 27 , 105 –113 10.1111/nmo.12486 25443141 49. Felice V.D. , Gibney S.M. , Gosselin R.D. , Dinan T.G. , O'Mahony S.M. and Cryan J.F. (2014 ) Differential activation of the prefrontal cortex and amygdala following psychological stress and colorectal distension in the maternally separated rat . Neuroscience 267 , 252 –262 10.1016/j.neuroscience.2014.01.064 24513388 50. Coutinho S.V. , Plotsky P.M. , Sablad M. , Miller J.C. , Zhou H. , Bayati A.I. et al. (2002 ) Neonatal maternal separation alters stress-induced responses to viscerosomatic nociceptive stimuli in rat . Am. J. Physiol. Gastrointest. Liver Physiol. 282 , G307 –G316 10.1152/ajpgi.00240.2001 11804852 51. O' Mahony S.M. , Coelho A.M. , Fitzgerald P. , Lee K. , Winchester W. , Dinan T.G. et al. (2011 ) The effects of gabapentin in two animal models of co-morbid anxiety and visceral hypersensitivity . Eur. J. Pharmacol. 667 , 169 –174 10.1016/j.ejphar.2011.05.055 21645509 52. Buckingham J.C. , John C.D. , Solito E. , Tierney T. , Flower R.J. , Christian H. et al. (2006 ) Annexin 1, glucocorticoids, and the neuroendocrine-immune interface . Ann. N.Y. Acad. Sci. 1088 , 396 –409 10.1196/annals.1366.002 17192583 53. Zhang Y.K. , Huang Z.J. , Liu S. , Liu Y.P. , Song A.A. and Song X.J. (2013 ) WNT signaling underlies the pathogenesis of neuropathic pain in rodents . J. Clin. Invest. 123 , 2268 –2286 10.1172/JCI65364 23585476 54. Louis P. and Flint H.J. (2009 ) Diversity, metabolism and microbial ecology of butyrate-producing bacteria from the human large intestine . FEMS Microbiol. Lett. 294 , 1 –8 10.1111/j.1574-6968.2009.01514.x 19222573 55. Atarashi K. , Tanoue T. , Shima T. , Imaoka A. , Kuwahara T. , Momose Y. et al. (2011 ) Induction of colonic regulatory T cells by indigenous Clostridium species . Science 331 , 337 –341 10.1126/science.1198469 21205640 56. Andoh A. , Tsujikawa T. , Sasaki M. , Mitsuyama K. , Suzuki Y. , Matsui T. et al. (2009 ) Faecal microbiota profile of Crohn's disease determined by terminal restriction fragment length polymorphism analysis . Aliment. Pharmacol. Ther. 29 , 75 –82 10.1111/j.1365-2036.2008.03860.x 18945264 57. Duytschaever G. , Huys G. , Bekaert M. , Boulanger L. , De Boeck K. and Vandamme P. (2013 ) Dysbiosis of bifidobacteria and Clostridium cluster XIVa in the cystic fibrosis fecal microbiota . J. Cyst. Fibros. 12 , 206 –215 10.1016/j.jcf.2012.10.003 23151540 58. Ubeda C. , Bucci V. , Caballero S. , Djukovic A. , Toussaint N.C. , Equinda M. et al. (2013 ) Intestinal microbiota containing Barnesiella species cures vancomycin-resistant Enterococcus faecium colonization . Infect. Immun. 81 , 965 –973 10.1128/IAI.01197-12 23319552 59. Socransky S.S. and Haffajee A.D. (2005 ) Periodontal microbial ecology . Periodontol 2000 38 , 135 –187 10.1111/j.1600-0757.2005.00107.x 15853940 60. Golubeva A.V. , Joyce S.A. , Moloney G. , Burokas A. , Sherwin E. , Arboleya S. et al. (2017 ) Microbiota-related Changes in Bile Acid & Tryptophan Metabolism are Associated with Gastrointestinal Dysfunction in a Mouse Model of Autism . EBioMedicine 24 , 166 –178 28965876 61. Kverka M. , Zakostelska Z. , Klimesova K. , Sokol D. , Hudcovic T. , Hrncir T. et al. (2011 ) Oral administration of Parabacteroides distasonis antigens attenuates experimental murine colitis through modulation of immunity and microbiota composition . Clin. Exp. Immunol. 163 , 250 –259 10.1111/j.1365-2249.2010.04286.x 21087444 62. Cekanaviciute E. , Yoo B.B. , Runia T.F. , Debelius J.W. , Singh S. , Nelson C.A. et al. (2017 ) Gut bacteria from multiple sclerosis patients modulate human T cells and exacerbate symptoms in mouse models . Proc. Natl. Acad. Sci. U.S.A. 114 , 10713 –10718 10.1073/pnas.1711235114 28893978 63. Del Chierico F. , Nobili V. , Vernocchi P. , Russo A. , De Stefanis C. , Gnani D. et al. (2017 ) Gut microbiota profiling of pediatric nonalcoholic fatty liver disease and obese patients unveiled by an integrated meta-omics-based approach . Hepatology 65 , 451 –464 10.1002/hep.28572 27028797 64. Wang K. , Liao M. , Zhou N. , Bao L. , Ma K. , Zheng Z. et al. (2019 ) Parabacteroides distasonis Alleviates Obesity and Metabolic Dysfunctions via Production of Succinate and Secondary Bile Acids . Cell Rep. 26 , 222e5 –235e5 10.1016/j.celrep.2018.12.028 30605678 65. O'Mahony S.M. , McVey Neufeld K.A. , Waworuntu R.V. , Pusceddu M.M. , Manurung S. , Murphy K. et al. (2019 ) The enduring effects of early-life stress on the microbiota-gut-brain axis are buffered by dietary supplementation with milk fat globule membrane and a prebiotic blend . Eur. J. Neurosci. 51 , 1042 –1058 10.1111/ejn.14514 31339598 66. Sokol H. , Pigneur B. , Watterlot L. , Lakhdari O. , Bermudez-Humaran L.G. , Gratadoux J.J. et al. (2008 ) Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients . Proc. Natl. Acad. Sci. U.S.A. 105 , 16731 –16736 10.1073/pnas.0804812105 18936492 67. Martin R. , Miquel S. , Chain F. , Natividad J.M. , Jury J. , Lu J. et al. (2015 ) Faecalibacterium prausnitzii prevents physiological damages in a chronic low-grade inflammation murine model . BMC Microbiol. 15 , 67 10.1186/s12866-015-0400-1 25888448 68. Vital M. , Howe A.C. and Tiedje J.M. (2014 ) Revealing the bacterial butyrate synthesis pathways by analyzing (meta)genomic data . MBio 5 , e00889 10.1128/mBio.00889-14 24757212