Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2017.00653 Microbiology Original Research Rice Bran and Probiotics Alter the Porcine Large Intestine and Serum Metabolomes for Protection against Human Rotavirus Diarrhea Nealon Nora Jean 1 Yuan Lijuan 2 Yang Xingdong 3 Ryan Elizabeth P. 1 * 1Nutrition and Toxicology Laboratory, Department of Environmental and Radiological Health Sciences, College of Veterinary Medicine and Biomedical Sciences, Colorado State University, Fort Collins CO, USA 2Yuan Laboratory, Department of Biomedical Sciences and Pathobiology, Virginia-Maryland Regional College of Veterinary Medicine, Virginia Polytechnic Institute and State University, Blacksburg VA, USA 3Laboratory of Infectious Diseases, Viral Pathogenesis and Evolution Section, National Institute of Allergy and Infectious Diseases, National Institute of Health, Bethesda MD, USA

Edited by: Jack Wong, The Chinese University of Hong Kong, Hong Kong

Reviewed by: Hongxia Wang, University of Alabama at Birmingham, USA; Vishvanath Tiwari, Central University of Rajasthan, India

*Correspondence: Elizabeth P. Ryan, e.p.ryan@colostate.edu

This article was submitted to Antimicrobials, Resistance and Chemotherapy, a section of the journal Frontiers in Microbiology

 21 04 2017 2017 8 653 11 10 2016 30 03 2017 Copyright © 2017 Nealon, Yuan, Yang and Ryan. 2017 Nealon, Yuan, Yang and Ryan

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Human rotavirus (HRV) is a leading cause of severe childhood diarrhea, and there is limited vaccine efficacy in the developing world. Neonatal gnotobiotic pigs consuming a prophylactic synbiotic combination of probiotics and rice bran (Pro+RB) did not exhibit HRV diarrhea after challenge. Multiple immune, gut barrier protective, and anti-diarrheal mechanisms contributed to the prophylactic efficacy of Pro+RB when compared to probiotics (Pro) alone. In order to understand the molecular signature associated with diarrheal protection by Pro+RB, a global non-targeted metabolomics approach was applied to investigate the large intestinal contents and serum of neonatal gnotobiotic pigs. The ultra-high performance liquid chromatography-tandem mass spectrometry platform revealed significantly different metabolites (293 in LIC and 84 in serum) in the pigs fed Pro+RB compared to Pro, and many of these metabolites were lipids and amino acid/peptides. Lipid metabolites included 2-oleoylglycerol (increased 293.40-fold in LIC of Pro+RB, p = 3.04E-10), which can modulate gastric emptying, andhyodeoxycholate (decreased 0.054-fold in the LIC of Pro+RB, p = 0.0040) that can increase colonic mucus production to improve intestinal barrier function. Amino acid metabolites included cysteine (decreased 0.40-fold in LIC, p = 0.033, and 0.62-fold in serum, p = 0.014 of Pro+RB), which has been found to reduce inflammation, lower oxidative stress and modulate mucosal immunity, and histamine (decreased 0.18-fold in LIC, p = 0.00030, of Pro+RB and 1.57-fold in serum, p = 0.043), which modulates local and systemic inflammatory responses as well as influences the enteric nervous system. Alterations to entire LIC and serum metabolic pathways further contributed to the anti-diarrheal and anti-viral activities of Pro+RB such as sphingolipid, mono/diacylglycerol, fatty acid, secondary bile acid, and polyamine metabolism. Sphingolipid and long chain fatty acid profiles influenced the ability of HRV to both infect and replicate within cells, suggesting that Pro+RB created a protective lipid profile that interferes with HRV activity. Polyamines act on enterocyte calcium-sensing receptors to modulate intracellular calcium levels, and may directly interfere with rotavirus replication. These results support that multiple host and probiotic metabolic networks, notably those involving lipid and amino acid/peptide metabolism, are important mechanisms through which Pro+RB protected against HRV diarrhea in neonatal gnotobiotic pigs.

 probiotic prebiotic synbiotic human rotavirus diarrhea metabolomics rice bran pigs OPP1015267 Bill and Melinda Gates Foundation10.13039/100000865 Introduction

Globally, HRV is the most common cause of severe diarrhea in children under 5, and annually it is responsible for approximately 450,000 deaths worldwide (Clarke and Desselberger, 2015). Live attenuated vaccines, including RV5 and Rotarix, have greatly reduced overall mortality and gastroenteritis in many international vaccine programs, yet they have a markedly lower efficacy in developing countries where most HRV outbreaks still occur (Yen et al., 2014). Underlying pathological processes in the host intestinal tract may contribute to this reduced HRV vaccine efficacy and include environmental enteric dysfunction, where abnormal intestinal barrier and mucosal immune functions result primarily or secondarily from dysbiosis of the microbiome (Keusch et al., 2016). The combination of limited vaccine efficacy and compromised gut function warrants the need for more effective preventive strategies against HRV, especially in developing countries.

Probiotic bacteria, including LGG and EcN, represent a safe, alternative therapeutic that can reduce the severity of HRV diarrhea. In humans and animals, probiotics modulate mucosal immunity, reduced HRV binding and infection, and produced antimicrobial peptides with anti-viral activity against HRV (Chenoll et al., 2016; Olaya Galan et al., 2016; Vlasova et al., 2016; Wang et al., 2016). In gnotobiotic pigs, LGG demonstrated a strong adjuvant effect when supplemented with an oral attenuated HRV vaccine and increased mucosal populations of HRV-specific IFN-gamma producing T lymphocytes (Wen et al., 2015). Collectively, these results support that probiotics can function as an effective component of a prophylactic HRV treatment.

The ability of probiotics to combat HRV can be enhanced when they are combined into a synbiotic with appropriate prebiotics, which are exclusive nutrient sources for probiotic bacteria that cannot be digested by the host. RB, the outer covering of the rice grain, represents an affordable, sustainable, and globally produced source of prebiotics. Prior research demonstrated that gnotobiotic neonatal pigs fed a prophylactic synbiotic combination of LGG/EcN and RB (Pro+RB) did not experience diarrhea after oral challenge with a high dose of live, virulent HRV, whereas pigs consuming LGG/EcN alone (Pro) or RB showed only reductions in diarrhea (Yang et al., 2015). In piglets consuming Pro+RB, diarrhea elimination was associated with normal tissue histology and healthy levels of serum gut permeability markers throughout the HRV challenge (Yang et al., 2015). Additional studies with this pig model have demonstrated the ability of Pro+RB to substantially reduce human norovirus-associated diarrhea, while stimulating more potent adaptive immune responses against viral antigens, and preserving colonic tissue architecture during infection (Lei et al., 2016). These findings suggested that the synbiotic combination of a probiotic and a prebiotic functioned through multiple immune, gut barrier protective and anti-diarrheal mechanisms to enhance protection against HRV diarrhea. The increased efficacy may arise in part through RB modulations of probiotic metabolism (Henderson et al., 2012).

Metabolomics, the systematic study of metabolites, metabolic pathways, and their interconnected networks, can be utilized to evaluate how RB alters probiotic function to prevent diarrhea. The objective of this study was to compare the LIC and serum metabolome of pigs consuming a prophylactic combination of Pro+RB to those consuming probiotics alone. This metabolic fingerprint shows how alterations to metabolic networks and pathways function to protect the host against HRV diarrhea. The hypothesis is that the LIC and serum metabolome of neonatal gnotobiotic pigs consuming Pro+RB contains pathways and networks of metabolites associated with enhanced anti-diarrheal and gut mucosal immune-modulatory activities.

 Materials and Methods Experimental Design and Sample Collection

Neonatal piglets (Large White crossbred) were reared in gnotobiotic isolators and maintained on experimental treatments as described previously by Yang et al. (2015). Briefly, piglets were started on a diet of ultra-pasteurized bovine milk, and from PND 3 through 40, the animals were separated into dietary treatments. Piglets in Pro group were inoculated with a sub-therapeutic dosage of 1 × 104 LGG/EcN on PND 3, 5, and 7 to ensure adequate colonization of the intestinal tract. Animals in the Pro+RB treatment group were also fed heat-stabilized Calrose RB daily at 10% of their caloric intake. On PND 40, piglets were euthanized and LIC and serum samples were collected from each animal. The experimental design, including a time course for sample collection, is included as Figure 1. All animal care and use procedures were approved by the Institutional Animal Care and Use Committee at Virginia Polytechnic Institute and State University and all ample collection and experimental procedures were conducted with accordance to the approved guidelines.

Experimental design and sample collection. ∗PND indicates postnatal day. ∗∗LGG/EcN refers to a 1:1 cocktail of Lactobacillus rhamnosus GG/Escherichia coli Nissle, which was provided to the appropriate treatment groups on postnatal days 3, 5, and 7. RB refers to Calrose rice bran.

 Metabolome Sample Preparation

To evaluate differences between Pro+RB and Pro treatments, a metabolomic analysis was performed by Metabolon Inc. © (Durham, NC, USA). Briefly, 0.5 mL of LIC or serum were collected in triplicate from each animal, stored at -80°C in microcentrifuge tubes, and were shipped to Metabolon on dry ice. Samples were stored at -80°C until processing. Prior to metabolite detection, a methanol solution was added to each sample to extract proteins and improve the recovery of small molecules. Samples were vigorously shaken for 2 min and then centrifuged. The resulting extracts were split into five parts for analysis via UPLC-MS/MS: two replicates for reverse phase UPLC-MS/MS using positive ion mode electrospray ionization, a third replicate for reverse phase UPLC-MS/MS with negative ion mode electrospray ionization, a fourth replicate for HPLC-MS/MS negative ion mode electrospray ionization, and a fifth replicate saved as a back-up. To remove any remaining organic solvent from the fractionation, extracts were placed on a TurboVap © (Zymark), and were then stored in liquid nitrogen prior to analysis.

 UPLC-MS/MS Analysis

Metabolite detection was accomplished using a Waters ACQUITY UPLC, a Thermo Scientific Q-exactive heated electrospray ionization source, and an Orbitrap mass analyzer operated at a 35,000 mass resolution. Each extract was dried and reconstituted in solvents appropriate for each detection method. Internal standards were added for quality control purposes. Acidic positive ion mode replicates were optimized for detection of either hydrophobic or hydrophilic compounds by gradient eluting in a mobile phase of water, methanol and 0.05% perfluoropentanoic acid and 0.1% formic acid (hydrophilic compounds), or methanol, acetonitrile, water, 0.05% perfluoropentanoic acid and 0.01% formic acid. All acidic positive ion mode replicates were ran on a C18 column (Waters UPLC BEH C100 18-2.1 mm × 150 mm, 1.7 μm). One basic ion extracts were eluted using an identical, separate C18 column with methanol, water, and 6.5 mM of Ammonium Bicarbonate at pH 8. The fourth extract was eluted using a high-performance liquid chromatography column (Waters UPLC BEH Amide 2.1 mm × 150 mm, 1.7 μm) with water, acetonitrile and 10 mM ammonium formate at pH 10.8. Mass spectrometry was performed using MS and data-dependent MSn scans using dynamic exclusion, and the total scan range covered 70–10,000 m/z.

 Data Extraction and Compound Identification

Raw data was extracted, peak-identified and quality-controlled as described previously (Brown et al., 2016). Briefly, compound identifications were compared to an internal library of over 3300 purified standards or recurrent unknown entities and were matched to library entities based on a narrow retention time/index, a mass to charge ratio within 10 parts per million of the standard, and chromatographic data including MS/MS forward and reverse scores between the experimental data and authentic standards. The raw counts of each sample were quantified using area under the curve and converted into relative abundances which were median-scaled to one. For each metabolite, fold difference was calculated by dividing the scaled relative abundance of Pro+RB by Pro.

 Metabolic Pathway Analysis

Identified metabolites were organized into metabolic pathways based on their biochemical and physiological properties. Pathway enrichment scores were calculated by dividing the number of significant metabolites in a pathway (k) by the total number of detected metabolites in the pathway (m). This value was then divided by the fraction of the total number of significant metabolites in the data set (n) over the total number of detected metabolites in the complete dataset (N):

k/mn/N

Metabolic pathways with enrichment scores greater than one indicated that the pathway contained more metabolites with statistically significant fold differences compared to all other pathways within the matrix. In this investigation, there were 20 LIC and 16 serum pathways with a score greater than one. An additional layer to score interpretation was applied such that pathways containing only metabolites with established functions related to anti-diarrheal activity were selected for significance in results interpretations and further network analysis.

For both LIC and serum lipid metabolites, pathway enrichment scores were calculated for the following pathways: steroid, lysolipid and lysoplasmalogen, acyl choline, glycine and carnitine, phospholipid, dicarboxylate fatty acid, mevalonate, endocannabinoid, primary bile acid, medium chain fatty acid, amide and branched fatty acid, mono/dihydroxy fatty acid, sterols, polyunsaturated (n3 and n6) fatty acid, long chain fatty acid, secondary bile acid, mono/diacylglycerol, sphingolipid, ketone and glycerolipid. For LIC and serum amino acids/peptides, pathway enrichment scores were calculated for the following pathways: leucine isoleucine and valine, tryptophan, dipeptide and dipeptide derivative, phenylalanine and tyrosine, glycine, serine and threonine, urea, arginine and proline, lysine, glutamate, methionine, cysteine, taurine and S-adenosyl methionine, glutathione, polyamine, histidine, creatine, alanine and aspartate, gamma-glutamyl amino acid, and guanidino and acetamido. Pathway enrichment scores were visualized using GraphPad Prism 6.07 (San Diego, CA, USA).

 Metabolic Network Analysis

Pathways containing metabolites with established functions related to anti-diarrheal activity were selected for visualization using Cytoscape 2.8.3 © software as previously described (Brown et al., 2016). Visualization of Pro+RB versus Pro included all metabolites from selected lipid and amino acid/peptide pathways in either the LIC or serum matrix. Node diameter was determined by the magnitude of the fold difference between Pro+RB and Pro, and node color was determined by the direction of fold difference when comparing Pro+RB and Pro. Pathway enrichment scores are presented as a number on each pathway node. The LIC lipid pathways visualized in Cytoscape were: phospholipid, sphingolipid, mono/diacylglycerol, polyunsaturated (n3 and n6) fatty acid, endocannabinoid, sterol, long chain fatty acid, mevalonate, medium chain fatty acid, dihydroxy fatty acid, primary bile acid, and secondary bile acid. A visual of the serum lipid pathways included: mono/dihydroxy fatty acid, monoacylglycerol, endocannabinoid, phospholipid, and dicarboxylate fatty acid. A visual of the LIC amino acid/peptide pathways included: glutamate, histidine, urea, arginine and proline, tryptophan, polyamine, alanine and aspartate, methionine, cysteine, S-adenosyl methionine and taurine, carnitine, dipeptide/derivatives, and gamma-glutamyl amino acid. A visual of the serum amino acid/peptide pathways included: guanidino and acetamido, creatine, histidine, glutamate, urea, arginine and proline, and methionine, cysteine, taurine and S-adenosyl methionine.

 Statistical Analysis

Statistical analysis was performed by Metabolon Inc. © using ArrayStudio (Omicsoft, Cary, NC, USA). Briefly, the raw abundance of each metabolite from each sample of Pro or Pro+RB was median-scaled into a relative abundance. In both the LIC and serum, the contrast of Pro+RB versus Pro was selected for analysis from a larger, multi-treatment data set where a 2-way ANOVA was applied across all treatments (data not shown). Specifically, scaled metabolite abundance values were compared between Pro+RB and Pro using a Welch’s two-sample t-test. Statistical significance was determined at the level of p < 0.05. Due to the multiple comparisons inherent to the 2-way ANOVA analysis, a false discovery rate (q) was calculated for each metabolite

 Results Non-targeted, Global Metabolomics Revealed Differences between Probiotic+Rice Bran and Probiotic Treatments

Across Pro+RB and Pro treatments, 625 total metabolites were detected in the LIC and 532 metabolites were detected in the serum (Table 1). Seven hundred forty-six unique metabolites were detected across both matrices, where 415 were shared between LIC and serum, 214 were exclusive to the LIC, and 117 were exclusive to the serum. Of the 625 LIC metabolites, 170 were classified as amino acids, 31 peptides, 30 carbohydrates, 10 TCA cycle, 238 lipids, 39 nucleotides, 32 cofactors and vitamins, 54 phytochemicals/other. Of the 532 serum metabolites, 153 were classified as amino acids, 22 peptides, 21 carbohydrates, 9 TCA cycle, 229 lipids, 39 nucleotides, 25 cofactors and vitamins, and 32 phytochemicals.

Number of identified metabolites from large intestinal contents and serum for pigs consuming probiotics in the presence and absence of rice bran (p < 0.05).

Metabolite classification Large intestinal contents Serum
Amino acids 170∗ (↑7, ↓78)∗∗ 153 (↑14, ↓15)
Peptides 31 (↑2, ↓14) 24 (↑0, ↓2)
Carbohydrates 32 (↑2, ↓11) 21 (↑0, ↓1)
Energy metabolism 10 (↑0, ↓3) 9 (↑0, ↓0)
Lipids 257 (↑47, ↓58) 229 (↑15, ↓16)
Nucleotides 39 (↑1, ↓14) 39 (↑6, ↓1)
Cofactors and vitamins 32 (↑4,↓16) 25 (↑8, ↓0)
Phytochemical/Other 54 (↑15, ↓21) 32 (↑3, ↓3)
Total number of identified metabolites 625 (293) 532 (84)
∗Indicates the total number of metabolites detected in a given metabolite class. ∗∗Values in parentheses indicate the number of significantly different (p < 0.05) metabolites that increased (↑) or decreased (↓) when comparing the fold difference of its scaled relative abundance in Pro+RB by Pro.

There were 293 LIC metabolites and 84 serum metabolites with statistically significant (p < 0.05) fold differences when comparing Pro+RB to Pro (Table 1). In the LIC, 78 metabolites had a significantly higher scaled relative abundance in Pro+RB compared to Pro and included: 7 amino acids, 2 peptides, 2 carbohydrates, 47 lipids, 1 nucleotide, 4 cofactors/vitamins, and 15 phytochemical/other. Two hundred fifteen LIC metabolites had a significantly lower scaled relative abundance in Pro+RB compared to Pro and included: 78 amino acids, 14 peptides, 11 carbohydrates, 3 TCA cycle, 58 lipids, 14 nucleotides, 16 cofactors/vitamins, and 21 phytochemical/other.

In the serum, 45 metabolites had a significantly higher scaled relative abundance in Pro+RB compared to Pro and included: 14 amino acids, 15 lipids, 6 nucleotides, 8 cofactors/vitamins, and 3 phytochemical/other. Thirty-eight metabolites had a significantly lower scaled relative abundance in Pro+RB compared to Pro and included: 15 amino acids, 1 peptide, 1 carbohydrate, 17 lipids, 1 nucleotide, and 1 cofactor/vitamin.

 Lipid Metabolite Profiles Were Differentially Expressed between Probiotic+Rice Bran and Probiotic Treatment Groups and Were Associated with Anti-diarrheal Activity

Across both the LIC and serum, lipids represented ∼41 and 43% of the identified metabolites, respectively, and thus were examined more closely for their prospective anti-diarrheal-related properties. Supplementary Table 1 lists the significant (p < 0.05) lipids with the fold difference and p-values provided across both LIC and serum matrices, and Supplementary Table 3 presents a comprehensive literature search into the anti-diarrheal functions the lipids with statistically significant fold differences between Pro+RB and Pro. Analysis revealed 26 metabolites that possessed immunomodulatory, gut barrier protective and anti-viral functions related to anti-diarrheal activity. In the LIC, these lipids included long chain fatty acid (12 significant metabolites), polyunsaturated fatty acid (8 significant metabolites), phospholipid (5 significant metabolites), mono/diacylglycerol (14 significant metabolites), endocannabinoid (3 significant metabolites), sterol (4 significant metabolites), mevalonate (1 significant metabolite), primary bile acid (5 significant metabolites), secondary bile acid (7 significant metabolites), dihydroxy fatty acid (2 significant metabolites), medium chain fatty acid (3 significant metabolites), and sphingolipids (16 significant metabolites). In the serum, these lipids included mono/dihydroxy fatty acid (5 significant metabolites), endocannabinoid (1 significant metabolite), phospholipid (8 significant metabolites), and dicarboxylate fatty acid (5 significant metabolites).

 Amino Acid and Peptide Metabolite Profiles Differ between Probiotic+Rice Bran and Probiotic Treatment Groups and Are Associated with Anti-diarrheal Activity

The LIC and serum amino acid/peptide metabolites represented ∼27 and 28% of the total number of identified metabolites, respectively, and like lipids, were examined more closely for compounds with potential anti-diarrheal activity. Supplementary Table 2 lists the amino acid/peptide metabolites with statistically significant (p < 0.05) fold differences and p-values between Pro+RB and Pro across both LIC and serum matrices. Supplementary Table 4 presents a comprehensive literature search of the anti-diarrheal functions of amino acids/peptides with statistically significant fold differences between Pro+RB and Pro. Analysis revealed 22 metabolites that possessed immunomodulatory, gut barrier protective and anti-viral functions related to anti-diarrheal activity (Supplementary Table 4). In the LIC these amino acids/peptides included glutamate (6 significantly different metabolites), histidine (10 significantly different metabolites), urea cycle arginine and proline (10 significantly different metabolites), tryptophan (5 significantly different metabolites), polyamine (6 significantly different metabolites), alanine and aspartate (6 significantly different metabolites), methionine, cysteine, SAM, and taurine (11 significantly different metabolites), dipeptide/derivatives (6 significantly different metabolites), gamma-glutamyl amino acid (10 significantly different metabolites), and carnitine (2 significantly different metabolites). In the serum, these amino acid/peptide metabolites included urea cycle, arginine, and proline (5 significantly different metabolites), phenylalanine and tyrosine (2 significantly different metabolites), methionine, cysteine, SAM, and taurine (4 significantly different metabolites), histidine (6 significantly different metabolites), and glutamate (3 significantly different metabolites).

 Metabolic Pathways in the LIC and Serum Containing Anti-diarrheal Metabolites Were Associated with High Pathway Enrichment Scores

To examine how individual metabolic pathways contributed to Pro+RB versus Pro differences in protection from HRV diarrhea in pigs, pathway enrichment scores were calculated for each individual lipid (Figure 2A) and amino acid/peptide (Figure 2B) pathway in the LIC and serum as described in the methods. In general, many of these pathways exhibited a high pathway enrichment score (>1.00), supporting that these pathways are major contributors to treatment differences between Pro+RB and Pro. Lipid pathways in the LIC that contained anti-diarrheal metabolites included: phospholipid (score 0.41), sphingolipid (score 1.90), mono/diacylglycerol (score 1.20), polyunsaturated fatty acid (score 1.43), endocannabinoid (score 0.92), sterol (score 1.51), long chain fatty acid (score 1.84), mevalonate (score 0.75), medium chain fatty acid (score 1.07), mono/dihydroxy fatty acid (score 1.20), primary bile acid (score 0.97), and secondary bile acid (score 1.84). In the serum, lipid metabolic pathways containing anti-diarrheal metabolites included mono/dihydroxy fatty acid (score 2.11), endocannabinoid (score 3.16), phospholipid (score 1.49), dicarboxylate fatty acid (score 2.43). For lipids, metabolic pathways that contained anti-diarrheal metabolites in both the LIC and serum matrices included dicarboxylate fatty acid, endocannabinoid, phospholipid, and mono/dihydroxy fatty acid.

Lipid and amino acid/peptide pathway enrichment scores in the large intestinal contents and serum when comparing Probiotic + Rice Bran versus Probiotic. Pathway enrichment scores of (A) lipid pathways and (B) amino acid/peptide pathways are plotted for both large intestinal (white bars) and serum matrices (black bars). ∗ Indicates a pathway containing metabolites identified as having anti-diarrheal activity. In (A), FA indicates fatty acid, BA indicates bile acid, MAG/DAG indicates monoacyl/diacylglycerol, PUFA indicates polyunsaturated fatty acid, and LCFA indicates long chain fatty acid. In (B), SAM indicates S-adenosyl methionine, and AA indicates amino acid.

Amino acid/peptide metabolic pathways in the LIC containing anti-diarrheal metabolites included: glutamate (score 1.28), histidine (score 1.43), urea cycle, arginine, and proline (score 1.26), tryptophan (score 0.59), polyamine (score 1.43), alanine and aspartate (score 1.61), methionine, cysteine, SAM, and taurine (score 1.31), dipeptide/derivatives (score 0.64), gamma-glutamyl amino acids (score 1.95), and carnitine (score 1.07). In the serum, amino acid/peptide metabolic pathways containing anti-diarrheal metabolites included: urea cycle, arginine, and proline (score 1.86), phenylalanine and tyrosine (score 0.79), methionine, cysteine, SAM, and taurine (score 1.40), histidine (score 3.80), and glutamate (score 2.11). For amino acids, metabolic pathways that contained anti-diarrheal metabolites in both the LIC and serum matrices included histidine, glutamate, urea, arginine and proline, and methionine (Figure 2B).

 Cytoscape Analysis Reveals Network-wide Alterations in Lipid and Amino Acid Metabolite Profiles Associated with Anti-diarrheal Activity

The Cytoscape visualization tool classified individual metabolites into metabolic networks of their corresponding pathways. Lipid metabolites and pathways containing metabolites with anti-diarrheal activity were visualized in the LIC (Figure 3) and the serum (Figure 4). In the LIC, the lipid metabolic network displayed large increases in multiple lipid metabolites from Pro+RB when compared to Pro. Many of these increases corresponded to metabolites from mono/diacylglycerol, and sphingolipid pathways and included the anti-diarrheal associated metabolites 2-palmitoylglycerol (increased 7.08-fold in Pro+RB, p = 6.33E-05), 1-linoleoyl glycerol (increased 320.72-fold in Pro+RB, p = 1.04E-08), 2-linoleoylglycerol (increased 187.08-fold in Pro+RB, p = 9.27E-10), palmitoyl sphingomyelin (d18:1/16:0) (increased 5.42-fold in Pro+RB, p = 0.017), and sphingosine (decreased 0.29-fold in Pro+RB, p = 0.029) (Figure 3). There were LIC lipid pathways with large decreases in multiple lipid metabolites from Pro+RB when compared to Pro involving medium chain, long chain, and polyunsaturated fatty acids. Associated anti-diarrheal metabolites included stearate (decreased 0.65-fold in Pro+RB, p = 0.015), palmitate (decreased 0.59-fold in Pro+RB, p = 0.0087), oleate/vaccenate (decreased 0.56-fold in Pro+RB, p = 0.030), caprate (decreased 0.24-fold in Pro+RB, p = 0.0012), docaspentaenoate (22:5n3) (decreased 0.26-fold in Pro+RB, p = 0.032), eicosapentaenoate (20:5n3) (decreased 0.26-fold in Pro+RB, p = 0.041), and palmitoleate (16:1n7) (decreased 0.22-fold in Pro+RB, p = 0.00026) (Figure 3). In the serum, the lipid metabolic network displayed both large increases and decreases in lipid metabolites when Pro+RB was compared to Pro. Specifically, mono/dihydroxy and dicarboxylate fatty acid pathways exhibited large increases in the anti-diarrheal metabolites 12,13-dihydroxyoctadecenoic acid (increased 15.03-fold in Pro+RB, p = 0.00031), pimelate (increased 1.58-fold in Pro+RB, p = 0.010), and azelate (increased 1.73-fold in Pro+RB, p = 0.023) in Pro+RB compared to Pro. Phospholipid metabolism resulted in large decreases in 1-palmitoyl-2-palmitoleyol-glycerophosphocholine (decreased 0.71-fold in Pro+RB, p = 0.028), 1-palmitoyl-2-linoleyol glycerophosphocholine (decreased 0.59-fold in Pro+RB, p = 0.015), 1-stearoyl-2-oleoyl glycerophosphocholine (decreased 0.73-fold in Pro+RB, p = 0.045), and 1,2-dioleoyl glycerophosphocholine (decreased 0.65-fold in Pro+RB, p = 0.014) when comparing Pro+RB and Pro (Figure 4).

Cytoscape network analysis of select lipid pathways in the large intestinal contents when comparing Probiotic + Rice Bran and Probiotic. Cytoscape visualization of lipid metabolites in LIC metabolic pathways associated with anti-diarrheal activity across Pro+RB and Pro. For each metabolite, node diameter is proportional to the magnitude of the fold difference in Pro+RB compared to Pro. Node color indicates the direction of a metabolite’s fold difference: red indicates metabolites with a higher scaled abundance in Pro+RB (p < 0.05), blue indicates lower abundance in Pro+RB (p < 0.05), pink indicates trending higher in Pro+RB (0.05 < p < 0.10), and light blue indicates trending lower in Pro+RB (0.05 < p < 0.10). Black nodes indicate metabolites with fold differences that were not significantly altered between treatments. The pathway enrichment score is the number in the circle for each pathway.

Cytoscape network analysis of select lipid pathways in the serum when comparing Probiotic + Rice Bran and Probiotic. Cytoscape visualization of lipid metabolites in serum metabolic pathways associated with anti-diarrheal activity across Pro+RB and Pro. For each metabolite, node diameter is proportional to the magnitude of the fold difference in Pro+RB compared to Pro. Node color indicates the direction of a metabolite’s fold difference: red indicates metabolites with a higher scaled abundance in Pro+RB (p < 0.05), blue indicates lower abundance in Pro+RB (p < 0.05), pink indicates trending higher in Pro+RB (0.05 < p < 0.10), and light blue indicates trending lower in Pro+RB (0.05 < p < 0.10). Black nodes indicate metabolites with fold differences that were not significantly altered between treatments. The pathway enrichment score is the number in the circle for each pathway.

Figures 5, 6 display amino acid/peptide LIC and serum metabolic pathways, respectively, containing anti-diarrheal-associated metabolites. In the LIC, Pro+RB displayed large decreases in multiple metabolites, across numerous pathways, when compared to Pro. These decreases corresponded to metabolites from gamma-glutamyl amino acid, polyamine, alanine and aspartate, histidine, glutamate, and methionine, cysteine, SAM, and taurine metabolism and included the anti-diarrheal associated metabolites gamma-glutamyl valine (decreased 0.15-fold in Pro+RB, p = 0.0011), gamma-glutamyl-epsilon-lysine (decreased 0.061-fold in Pro+RB, p = 2.18E-05), spermidine (decreased 0.038-fold in Pro+RB, p = 1.64E-10), alanine (decreased 0.40-fold in Pro+RB, p = 0.026), histidine (decreased 0.48-fold in Pro+RB, p = 0.012), histamine (decreased 0.18-fold in Pro+RB, p = 0.00030), 3-methylhistidine (decreased 0.26-fold in Pro+RB, p = 2.70E-05), gamma-aminobutyrate (GABA) (decreased 0.46-fold in Pro+RB, p = 0.023), glutamate (decreased 0.36-fold in Pro+RB, p = 2.38E-05), cysteine (decreased 0.40-fold, p = 0.033), homocysteine (decreased 0.41-fold in Pro+RB, p = 0.029), and taurine (decreased 0.21-fold in Pro+RB, p = 0.00034) (Figure 5). In the serum, amino acid/peptide metabolic network displayed both multiple increases in histidine metabolites of Pro+RB compared to Pro, including the anti-diarrheal associated metabolite histamine (increased 1.57-fold in Pro+RB, p = 0.043) (Figure 6).

Cytoscape network analysis of select amino acid/peptide pathways in the large intestinal contents when comparing Probiotic + Rice Bran and Probiotic. Cytoscape visualization of amino acid/peptide metabolites in LIC metabolic pathways associated with anti-diarrheal activity across Pro+RB and Pro. For each metabolite, node diameter is proportional to the magnitude of the fold difference in Pro+RB compared to Pro. Node color indicates the direction of a metabolite’s fold difference: red indicates metabolites with a higher scaled abundance in Pro+RB (p < 0.05), blue indicates lower abundance in Pro+RB (p < 0.05), pink indicates trending higher in Pro+RB (0.05 < p < 0.10), and light blue indicates trending lower in Pro+RB (0.05 < p < 0.10). Black nodes indicate metabolites with fold differences that were not significantly altered between treatments. The pathway enrichment score is the number in the circle for each pathway.

Cytoscape network analysis of select amino acid/peptide pathways in the serum when comparing Probiotic + Rice Bran and Probiotic. Cytoscape visualization of amino acid/peptide metabolites in serum metabolic pathways associated with anti-diarrheal activity across Pro+RB and Pro. For each metabolite, node diameter is proportional to the magnitude of the fold difference in Pro+RB compared to Pro. Node color indicates the direction of a metabolite’s fold difference: red indicates metabolites with a higher scaled abundance in Pro+RB (p < 0.05), blue indicates lower abundance in Pro+RB (p < 0.05), pink indicates trending higher in Pro+RB (0.05 < p < 0.10), and light blue indicates trending lower in Pro+RB (0.05 < p < 0.10). Black nodes indicate metabolites with fold differences that were not significantly altered between treatments. The pathway enrichment score is the number in the circle for each pathway.

 Discussion

Non-targeted global metabolomics revealed substantial alterations to both lipid and amino acid/peptide metabolism in the LIC and serum (Table 1), and multiple metabolite and metabolic pathway contributions to prophylactic diarrhea protection (Figure 2 and Supplementary Tables 2, 3). This manuscript further highlights those compounds associated with protection against HRV diarrhea related to increased production of gut-protective, immunomodulatory, and anti-viral compounds.

Rice bran-mediated alterations to probiotic lipid metabolism represent one route through which a combination of probiotics and prebiotics create a gut-protective environment against HRV infection. RB provides probiotic bacteria with a source of complex lipids that can be catabolized to release additional lipids into the large intestinal lumen or converted into bacterial lipid products (Buccioni et al., 2012). Alterations to LIC mono/diacylglycerol metabolism (pathway enrichment score 1.20) represent one area of lipid metabolism associated with anti-diarrheal activity (Figure 2). Elevations of multiple mono/diacylglycerols in the LIC of Pro+RB included 1-linoleoyl glycerol (increased 320.72 fold in Pro+RB, p = 1.04E-08), 2-linoleoyl glycerol (increased 187.09 fold in Pro+RB, p = 9.27E-10), and 2-palmitoyl glycerol (increased 7.08 fold in Pro+RB, p = 6.33E-05) (Supplementary Tables 1, 3). 1-linoleoyl glycerol has been demonstrated to systemically reduce inflammation in porcine models (Han et al., 2015). 2-linoleoyl glycerol serves as an important precursor in the formation of 2-arachidonyl glycerol that can broadly modulate inflammation through direct binding to endocannabinoid receptors, or via its catabolism into ligands of cyclooxygenases, cytochrome P450 enzymes, and lipoxygenase enzymes (Mechoulam et al., 1995; Turcotte et al., 2015). Similarly, 2-palmitoylglycerol is produced by mammalian and bacterial cells and influences cannabinoid receptor expression in the gastrointestinal tract (Aviello et al., 2008; Murataeva et al., 2016; Yuan et al., 2016). Collectively, probiotic catabolism of complex RB triacylglycerides may release mono/diacylglycerols into the intestinal lumen and increase their bioavailability in the intestinal lumen, where they can be utilized in host enterocyte metabolism to exert their anti-diarrhea effects. Reducing inflammation along the gastrointestinal tract may provide the host with improved barrier function (De Jong et al., 2016; Lei et al., 2016), and modulating cannabinoid receptor function can alter intestinal motility (Taschler et al., 2017) that can help prevent them against HRV-induced diarrhea. Ultimately, these released mono/diacylglycerols could be utilized in host enterocyte metabolism and promote anti-inflammatory and anti-viral activities that protect the host from HRV diarrhea.

In a healthy gut, primary bile acids play key roles in facilitating the digestion and absorption of fats, nutrients, and vitamins (Martinot et al., 2017). A major function of probiotic bacteria is the bioconversion of host-derived primary bile acids into bacterially modified secondary bile acids which function as antimicrobials, regulate host lipid and glucose metabolism, and influence the growth and activity of gut flora (Nie et al., 2015). Analysis of LIC primary and secondary bile acid metabolites in Pro+RB versus Pro revealed anti-diarrheal activities of the primary bile acid chenodeoxycholate (decreased 0.054-fold in Pro+RB, p = 0.0078), and the secondary bile acid hyodeoxycholate (decreased 0.054-fold in Pro+RB, p = 0.0040) (Supplementary Tables 1, 3). Investigations have suggested that lower luminal levels of chenodeoxycholate are associated with decreased gastric emptying times and reduced expression of the pro-inflammatory cytokines TNF-a, IL-1B, IL-6, and IL-8 (Goyal et al., 2015; Peleman et al., 2016). The decrease in chenodeoxycholate may be associated with metabolism of hyodeoxycholate a secondary bile acid resulting from bacterial metabolism of chenodeoxycholate (Sabbatini et al., 2013; Theriot et al., 2016). In the colonic lumen, rodent studies suggest that hyodeoxycholate acts directly on goblet cells to increase mucus production, thereby modulating barrier function (Barcelo et al., 2001). RB supplementation may have increased probiotic conversion of chenodeoxycholate into hyodeoxycholate, which provided, that ultimately was taken-up by enterocytes, out of the lumen, and used to mediate anti-inflammatory processes in the cell, and thus improve barrier function that provided protection during challenge with HRV.

Probiotic production of fatty acids, along with bioconversion or release of RB fatty acids, may influence immune and metabolic processes that contribute to diarrheal protection. The dicarboxylic fatty acid pimelate was elevated in 5.17-fold the LIC (p = 0.00033) and 1.58-fold the serum (p = 0.01) of Pro+RB versus Pro (Supplementary Tables 1, 3) and is produced by a variety of bacterial species. Past investigations have demonstrated that pimelate produced by the probiotic Bifidobacterium spp. in mice can alter gut luminal metabolism by increasing the bioavailability of biotin to enterocytes, modulating the development of T-regulatory cells, and influencing the composition of the gut microbiome (Sugahara et al., 2015). In healthy human and animal enterocytes, biotin is an essential co-factor produced by gut flora and it has multiple roles in lipid and amino acid metabolism, including as a modulator of mucosal immune responses and gut microbial populations (Ghosal et al., 2015). Although not established in this gnotobiotic model of disease, HRV infections are thought to contribute to disruption of normal intestinal microflora (Zhang et al., 2014), suggesting HRV pathogenesis may indirectly alter host biotin metabolism. Pimelate may have enhanced biotin uptake into enterocytes to promote gut metabolic, microbial and immune homeostasis, thereby improving barrier function and reducing inflammation, and may have thus allowed the host to better resist diarrhea.

Rice bran modification to probiotic amino acid/peptide metabolism also provided prophylactic diarrhea protection. Probiotic bacteria can catabolize complex RB and host proteins to release smaller bioactive peptides and amino acids, or bioconvert plant and host amino acids into bacterial-derived amino acids (Neis et al., 2015; Pessione and Cirrincione, 2016) to create a novel profile of amino acids/peptides that participate in anti-diarrheal and anti-viral responses against HRV. Probiotic and host histidine metabolism was differentially modulated between Pro+RB and Pro groups. Histidine, which was decreased 0.12-fold in the LIC of Pro+RB when compared to Pro (p = 0.012) (Supplementary Tables 2, 4), functioned as an anti-secretory agent to reduce Cholera-associated diarrhea when supplemented into a rice-based oral rehydration salt (Rabbani et al., 2005), suggesting that it could function similarly to reduce the onset of HRV diarrhea.

Histamine, a downstream metabolite of histidine that can be produced by both animals and probiotics (Gao et al., 2015), decreased 0.19-fold in the LIC (p = 0.00030) and increased 1.57-fold in the serum (p = 0.043) of Pro+RB compared to Pro (Supplementary Tables 2, 4). Histamine is a well-characterized mediator of inflammation, mucosal immunity, and the enteric nervous system (Smolinska et al., 2014), yet it can suppress these processes when bound to specific receptors in the colon (Gao et al., 2015). Previous investigations with Lactobacillus reuteri demonstrated that probiotics can use a bacterial-derived histidine deactylase to convert dietary histidine into histamine in the colonic lumen, which acts on H2 receptors to suppress colonic inflammation and reduce damage to the colonic epithelium (Gao et al., 2015). Probiotics may have converted RB histidine into histamine that left the LIC to bind H2 receptors, caused local modulation of colonic inflammatory responses, and increased barrier function, all of which ultimately interfered with HRV pathogenesis and provided diarrhea protection.

Other contributors to osmotic balance and thus likely anti-diarrheal mediators include cysteine and cystine metabolites. Cysteine, decreased 0.40-fold in the LIC (p = 0.033) and 0.62-fold in the serum (p = 0.014) of Pro+RB compared to Pro (Supplementary Tables 2, 4), reduces oxidative stress, decreases inflammation, lowers permeability, and modulates the mucosal immune system in the intestinal tract of neonatal pigs (Bauchart-Thevret et al., 2009; Ruth and Field, 2013). Cystine, increased 1.30-fold (p = 0.044) in the serum of Pro+RB compared to Pro (Supplementary Tables 2, 4), serves as an excitatory neurotransmitter in the brainstem of humans and other animals, where it functions peripherally to delay gastric emptying and suppress food intake (Khan et al., 2014; McGavigan et al., 2015). During HRV infections, inflammation can cause enterocyte damage and increase gut permeability which both exacerbate HRV diarrhea (Yang et al., 2015). Prophylactic modulations to cysteine and cysteine induced by Pro+RB may favorably shape the mucosal immune response and enteric nervous system to protect enterocytes and gut flora against oxidative damage and maintain normal gut motility during HRV infections to diarrhea.

Pro+RB alterations to entire pathways, containing multiple metabolites working in synergy, have been theorized to have anti-diarrheal effects on HRV. This was most evident in the LIC lipid pathways including long chain fatty acid (pathway enrichment score 1.84), sphingolipid (score 1.90, LIC), and secondary bile acid (score 1.84, LIC) pathways (Figures 2A, 3), all of which contained multiple anti-diarrheal-associated metabolites and had pathway enrichment scores higher than one. This suggests that these pathways were contributors to treatment differences contributing to the anti-diarrheal enhancement of Pro+RB versus Pro.Systems of long chain fatty acids have been reported to influence HRV infectivity and replication in the intestinal tract (Kim et al., 2012), whereas pathway-wide shifts in sphingolipid metabolism can influence the composition of host membrane lipid rafts (Barenholz, 2004), where HRV must successfully bind to infect enterocytes. Similarly, entire profiles of secondary bile acids have been associated with reductions in HRV diarrhea in Kenyan children, and likely function to disrupt HRV function by altering the pH of the intestinal lumen (Tazume et al., 1990).

Entire amino acid/peptide pathways associated with anti-diarrheal function included gamma-glutamyl amino acid (pathway enrichment score 1.95, LIC) and polyamine metabolism (score 1.43, LIC). HRV diarrhea results in part due to impaired transport of gamma-glutamyl amino acids from the intestinal lumen and into enterocytes (Katyal et al., 1999). Large differences in the LIC gamma-glutamyl amino acid profile between Pro+RB and Pro (Figures 2B, 5) suggest that prophylactic alterations to luminal amino acid transport are associated with anti-diarrheal effects. In addition, HRV de-regulates host calcium metabolism and induces massive intracellular calcium increases to support HRV replication (Hyser et al., 2013). Pathway-wide modulation of LIC polyamine metabolites (Figures 2B, 5) in Pro+RB can work to reduce diarrhea and interfere with HRV replication. Past investigations support that polyamines in the intestinal lumen act on the enterocyte calcium-sensing receptor to inhibit the intracellular release of calcium and reduce chloride secretion into the intestinal lumen (Rogers et al., 2015), simultaneously interfering with rotavirus replication and altering luminal osmolarity to reduce the onset of diarrhea. Collectively, analysis of entire lipid and amino acid/peptide pathways supports that a more comprehensive understanding of anti-diarrheal activity can be achieved when clusters of closely related metabolites are considered as part of an entire, functioning metabolic system.

While many pathways appeared to be differentially regulated exclusively in the LIC or serum, multiple lipid and amino acid/peptide pathways contained anti-diarrheal-associated metabolites in both LIC and serum matrices. These pathways included dicarboxylate fatty acid, endocannabinoid, phospholipid, and mono/dihydroxy fatty acid, histidine, glutamate, urea, arginine and proline, and methionine. When coupled with their high pathway enrichment scores, this suggests that these pathways may be important both locally and systemically for promoting enhanced diarrheal protection in Pro+RB compared to Pro. For example, the glutamate metabolism (pathway enrichment score 1.28, LIC and 2.11 serum) (Figure 2B), works through the central and enteric nervous systems to modulate gut motility, intestinal fluid balance, and inflammation (Kondoh et al., 2009; Li et al., 2012; Auteri et al., 2015). In addition, endocannabinoid metabolism, (score 0.92 LIC and 3.16 serum) works locally in the intestinal lumen to improve barrier function, modulate motility, and improve mucosal integrity and systemically to influence immune-mediated inflammatory responses and central nervous system functioning that could ultimately influence the course of diarrheal diseases in the gastrointestinal tract (Goodgame, 1999; Lee et al., 2016; Gertsch, 2017).

Currently, variable vaccine efficacy and environmental enteric dysfunction exacerbate HRV outbreaks in the developing countries and warrant the development of stronger preventive strategies. The combination of probiotics and RB has the potential to serve as a safe, natural, prophylactic treatment aimed at reducing HRV diarrhea, and accomplishes this in part via multiple modulations to bacterial and host lipid and amino acid/peptide metabolism. Probiotic metabolism of RB can exert direct anti-viral activity against HRV by interfering with HRV infection, replication and pathogenesis, but also interfere with HRV pathogenicity by improving gut barrier function, and modulating mucosal immune responses, especially inflammatory processes that influence the course of HRV diarrhea. Direct modulations to host metabolism and immunity also indicate that while Pro+RB is effective against HRV, it is likely that this combination can provide broad protection against a wide range of pathogenic diarrheal diseases. Although the current investigation focused primarily on modulations to lipid and amino acid/peptide metabolism, there are likely other classes of metabolites that contribute to the anti-viral and anti-diarrheal enhancement associated with combinations of probiotics and RB. An evaluation of these other metabolic pathways can serve as future investigations into the molecular mechanisms of diarrhea protection by this combination. Probiotic and RB prophylactic therapies address a global need for more accessible, natural, antimicrobial agents, and have the potential to protect people and animals against a wide variety of enteric diseases.

 Author Contributions

LY, XY and ER designed animal studies, and LY and XY conducted animal studies. NN and ER performed metabolomic data analysis and wrote the manuscript. All authors approved the final manuscript for submission.

 Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding. This work was funded by a Grand Challenges in Global Health Award from the Bill and Melinda Gates Foundation (OPP1015267).

The authors would like to thank Anna McClung from the USDA Rice Research Unit for supplying rice bran, the research group of Dr. Lijuan Yuan for their assistance with sample collection, and Dustin G. Brown and Colette R. Worcester for assistance with metabolomics data interpretation.

 Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00653/full#supplementary-material

 References Abe M. Iriki T. Kaneshige K. Kuwashima K. Watanabe S. Sato H. (2001). Adverse effects of excess lysine in calves. J. Anim. Sci. 79 1337–1345. 10.2527/2001.7951337x Akula M. K. Shi M. Jiang Z. Foster C. E. Miao D. Li A. S. (2016). Control of the innate immune response by the mevalonate pathway. Nat. Immunol. 17 922–929. 10.1038/ni.3487 Alp Avci G. (2016). Selection of superior bifidobacteria in the presence of rotavirus. Braz. J. Med. Biol. Res. 49:e5562. 10.1590/1414-431X20165562 Auteri M. Zizzo M. G. Serio R. (2015). GABA and GABA receptors in the gastrointestinal tract: from motility to inflammation. Pharmacol. Res. 93 11–21. 10.1016/j.phrs.2014.12.001 Aviello G. Romano B. Izzo A. A. (2008). Cannabinoids and gastrointestinal motility: animal and human studies. Eur. Rev. Med. Pharmacol. Sci. 12(Suppl. 1) 81–93. Baraniuk J. N. El-Amin S. Corey R. Rayhan R. Timbol C. (2013). Carnosine treatment for gulf war illness: a randomized, controlled trial. Glob. J. Health Sci. 5 69–81. 10.5539/gjhs.v5n3p69 Barcelo A. Claustre J. Toumi F. Burlet G. Chayvialle J. A. Cuber J. C. (2001). Effect of bile salts on colonic mucus secretion in isolated vascularly perfused rat colon. Dig. Dis. Sci. 46 1223–1231. 10.1023/A:1010607127822 Barenholz Y. (2004). Sphingomyelin and cholesterol: from membrane biophysics and rafts to potential medical applications. Subcell. Biochem. 37 167–215. 10.1007/978-1-4757-5806-1_5 Bauchart-Thevret C. Stoll B. Chacko S. Burrin D. G. (2009). Sulfur amino acid deficiency upregulates intestinal methionine cycle activity and suppresses epithelial growth in neonatal pigs. Am. J. Physiol. Endocrinol. Metab. 296 E1239–E1250. 10.1152/ajpendo.91021.2008 Brocker C. Thompson D. C. Vasiliou V. (2012). The role of hyperosmotic stress in inflammation and disease. Biomol. Concepts 3 345–364. 10.1515/bmc-2012-0001 Brown D. G. Rao S. Weir T. L. O’malia J. Bazan M. Brown R. J. (2016). Metabolomics and metabolic pathway networks from human colorectal cancers, adjacent mucosa, and stool. Cancer Metab. 4 11. 10.1186/s40170-016-0151-y Brundige D. R. Maga E. A. Klasing K. C. Murray J. D. (2010). Consumption of pasteurized human lysozyme transgenic goats’ milk alters serum metabolite profile in young pigs. Transgenic Res. 19 563–574. 10.1007/s11248-009-9334-4 Buccioni A. Decandia M. Minieri S. Molle G. Cabiddu A. (2012). Lipid metabolism in the rumen: new insights on lipolysis and biohydrogenation with an emphasis on the role of endogenous plant factors. Anim. Feed Sci. Technol. 174 1–25. 10.1016/j.anifeedsci.2012.02.009 Chenoll E. Casinos B. Bataller E. Buesa J. Ramon D. Genoves S. (2016). Identification of a peptide produced by Bifidobacterium longum CECT 7210 with antirotaviral activity. Front. Microbiol. 7:655. 10.3389/fmicb.2016.00655 Clarke E. Desselberger U. (2015). Correlates of protection against human rotavirus disease and the factors influencing protection in low-income settings. Mucosal Immunol. 8 1–17. 10.1038/mi.2014.114 Collino S. Montoliu I. Martin F. P. Scherer M. Mari D. Salvioli S. (2013). Metabolic signatures of extreme longevity in northern Italian centenarians reveal a complex remodeling of lipids, amino acids, and gut microbiota metabolism. PLoS ONE 8:e56564. 10.1371/journal.pone.0056564 Cordaro M. Impellizzeri D. Gugliandolo E. Siracusa R. Crupi R. Esposito E. (2016). Adelmidrol, a palmitoylethanolamide analogue, as a new pharmacological treatment for the management of inflammatory bowel disease. Mol. Pharmacol. 90 549–561. 10.1124/mol.116.105668 D’Argenio G. Calvani M. Della Valle N. Cosenza V. Di Matteo G. Giorgio P. (2005). Differential expression of multiple transglutaminases in human colon: impaired keratinocyte transglutaminase expression in ulcerative colitis. Gut 54 496–502. 10.1136/gut.2004.049411 De Jong P. R. González-Navajas J. M. Jansen N. J. G. (2016). The digestive tract as the origin of systemic inflammation. Crit. Care 20 279. 10.1186/s13054-016-1458-3 Engels F. Kessels G. C. Henricks P. A. Nijkamp F. P. (1996). Preferential formation of 13-hydroxylinoleic acid by human peripheral blood eosinophils. Prostaglandins 52 117–124. 10.1016/0090-6980(96)00057-3 Frigolet M. E. Gutiérrez-Aguilar R. (2017). The role of the novel lipokine palmitoleic acid in health and disease. Adv. Nutr. 8 173S–181S. 10.3945/an.115.011130 Fujimura K. E. Sitarik A. R. Havstad S. Lin D. L. Levan S. Fadrosh D. (2016). Neonatal gut microbiota associates with childhood multisensitized atopy and T cell differentiation. Nat. Med. 22 1187–1191. 10.1038/nm.4176 Gao C. Major A. Rendon D. Lugo M. Jackson V. Shi Z. (2015). Histamine H2 receptor-mediated suppression of intestinal inflammation by probiotic Lactobacillus reuteri. MBio 6 6. 10.1128/mBio.01358-15 Gertsch J. (2017). Cannabimimetic phytochemicals in the diet–an evolutionary link to food selection and metabolic stress adaptation? Br. J. Pharmacol. 10.1111/bph.13676 [Epub ahead of print]. Ghosal A. Jellbauer S. Kapadia R. Raffatellu M. Said H. M. (2015). Salmonella infection inhibits intestinal biotin transport: cellular and molecular mechanisms. Am. J. Physiol. Gastrointest. Liver Physiol. 309 2. 10.1152/ajpgi.00112.2015 Goodgame R. W. (1999). Viral infections of the gastrointestinal tract. Curr. Gastroenterol. Rep. 1 292–300. 10.1007/s11894-999-0112-5 Goyal N. Rana A. Bijjem K. R. Kumar P. (2015). Effect of chenodeoxycholic acid and sodium hydrogen sulfide in dinitro benzene sulfonic acid (DNBS)-Induced ulcerative colitis in rats. Pharmacol. Rep. 67 616–623. 10.1016/j.pharep.2014.12.018 Gutiérrez C. Villa S. Mota F. R. Calva J. J. (2007). Does an L-glutamine-containing, Glucose-free, oral rehydration solution reduce stool output and time to rehydrate in children with acute diarrhoea? a double-blind randomized clinical trial. J. Health Popul. Nutr. 25 278–284. Han X. Wang T. Zhang J. Liu X. Li Z. Wang G. (2015). Apolipoprotein CIII regulates lipoprotein-associated phospholipase A2 expression via the MAPK and NFκB pathways. Biol. Open 4 661–665. 10.1242/bio.201410900 Hasobe M. McKee J. G. Ishii H. Cools M. Borchardt R. T. De Clerq E. (1989). Elucidation of the mechanism by which homocysteine potentiates the anti-vaccinia virus effects of the S-adenosylhomocysteine hydrolase inhibitor 9-(trans-2’,trans-3’-dihydroxycyclopent4’-enyl)-adenine. Mol. Pharmacol. 36 490–496. Henderson A. J. Kumar A. Barnett B. Dow S. W. Ryan E. P. (2012). Consumption of rice bran increases mucosal immunoglobulin A concentrations and numbers of intestinal Lactobacillus spp. J. Med. Food 15 469–475. 10.1089/jmf.2011.0213 Hyser J. M. Utama B. Crawford S. E. Broughman J. R. Estes M. K. (2013). Activation of the endoplasmic reticulum calcium sensor STIM1 and store-operated calcium entry by rotavirus requires NSP4 viroporin activity. J. Virol. 87 13579–13588. 10.1128/JVI.02629-13 Ishida T. Nishiumi S. Tanahashi T. Yamasaki A. Yamazaki A. Akashi T. (2013). Linoleoyl ethanolamide reduces lipopolysaccharide-induced inflammation in macrophages and ameliorates 2,4-dinitrofluorobenzene-induced contact dermatitis in mice. Eur. J. Pharmacol. 699 6–13. 10.1016/j.ejphar.2012.11.030 Jourdan N. Brunet J. P. Sapin C. Blasi A. Cotte-Laffitte J. Forestier F. (1998). Rotavirus infection reduces sucrase-isomaltase expression in human intestinal epithelial cells by perturbing protein targeting and organization of microvillar cytoskeleton. J. Virol. 72 7228–7236. Katyal R. Rana S. V. Vaiphei K. Ohja S. Singh K. Singh V. (1999). Effect of rotavirus infection on small gut pathophysiology in a mouse model. J. Gastroenterol. Hepatol. 14 779–784. 10.1046/j.1440-1746.1999.01948.x Keusch G. T. Fischer Walker C. Das J. K. Horton S. Habte D. (2016). “Diarrheal Diseases,” in Reproductive, Maternal, Newborn and Child Health: Disease Control Priorities eds Black R. Laxminarayan R. Temmerman M. Walker N. (Washington, DC: The World Bank). Khan M. T. Van Dijil J. M. Harmsen H. J. M. (2014). Antioxidants keep the potentially probiotic but highly oxygen-sensitive human gut bacterium Faecalibacterium prausnitzii alive at ambient air. PLoS ONE 9:e96097. 10.1371/journal.pone.0096097 Kim Y. George D. Prior A. M. Prasain K. Hao S. Le D. D. (2012). Novel triacsin C analogs as potential antivirals against rotavirus infections. Eur. J. Med. Chem. 50 311–318. 10.1016/j.ejmech.2012.02.010 Kondoh T. Mallick H. N. Torii K. (2009). Activation of the gut-brain axis by dietary glutamate and physiologic significance in energy homeostasis. Am. J. Clin. Nutr. 90 832S–837S. 10.3945/ajcn.2009.27462V Kong W. Wang J. Ping X. Shen J. Ni X. Liu F. (2015). Biomarkers for assessing mucosal barrier dysfunction induced by chemotherapy: identifying a rapid and simple biomarker. Clin. Lab. 61 371–378. 10.7754/Clin.Lab.2014.140712 Lee Y. Jo J. Chung H. Y. Pothoulakis C. Im E. (2016). Endocannabinoids in the gastrointestinal tract. Am. J. Physiol. Gastrointest. Liver Physiol. 311 G655–G666. 10.1152/ajpgi.00294.2015 Lei S. Ramesh A. Twitchell E. Wen K. Bui T. Weiss M. (2016). High protective efficacy of probiotics and rice bran against human norovirus infection and diarrhea in gnotobiotic pigs. Front. Microbiol. 7:1699. 10.3389/fmicb.2016.01699 Li Y. Xiang Y. Y. Lu W. Y. Liu C. Li J. (2012). A novel role of intestine epithelial GABAergic signaling in regulating intestinal fluid secretion. Am. J. Physiol. Gastrointest. Liver Physiol. 303 G453–G460. 10.1152/ajpgi.00497.2011 Lin H. C. Van Citters G. W. Heimer F. Bonorris G. (2001). Slowing of gastrointestinal transit by oleic acid: a preliminary report of a novel, nutrient-based treatment in humans. Dig. Dis. Sci. 46 223–229. 10.1023/A:1005676811378 Liu G. Fang T. Yan T. Jia G. Zhao H. Huang Z. (2014). Metabolomic strategy for the detection of metabolic effects of spermine supplementation in weaned rats. J. Agric. Food Chem. 62 9035–9042. 10.1021/jf500882t Lustgarten M. S. Fielding R. A. (2016). Metabolites associated with circulating interleukin-6 in older adults. J Gerontol. A Biol. Sci. Med. Sci. 10.1093/gerona/glw039 [Epub ahead of print]. Mandøe J. M. Hansen K. B. Hartmann B. Rehfeld J. F. Holst J. J. Hansen H. S. (2015). The 2-monoacylglycerol moiety of dietary fat appears to be responsible for the fat-induced release of GLP-1 in humans. Am. J. Clin. Nutr. 102 548–555. 10.3945/ajcn.115.106799 Martinot E. Sèdes L. Baptissart M. Lobaccaro J. M. Caira F. Beaudoin C. (2017). Bile acids and their receptors. Mol. Aspects Med. 10.1016/j.mam.2017.01.006 [Epub ahead of print]. Matty A. J. Seshardi B. (1965). Effect of thyroxine on the isolated rat intestine. Gut 6 200–202. 10.1136/gut.6.2.200 McGavigan A. K. O’hara H. C. Amin A. Kinsey-Jones J. Spreckley E. Alamshah A. (2015). L-cysteine suppresses ghrelin and reduces appetite in rodents and humans. Int. J. Obes. 39 447–455. 10.1038/ijo.2014.172 Mechoulam R. Ben-Shabat S. Hanus L. Ligumsky M. Kaminski N. E. Schatz A. R. (1995). Identification of an endogenous 2-monoglyceride, present in canine gut, that binds to cannabinoid receptors. Biochem. Pharmacol. 50 83–90. 10.1016/0006-2952(95)00109-D Mehmood M. H. Anila N. Begum S. Syed S. A. Siddiqui B. S. Gilani A. H. (2014). Pharmacological basis for the medicinal use of Carissa carandas in constipation and diarrhea. J. Ethnopharmacol. 153 359–367. 10.1016/j.jep.2014.02.024 Mounce B. C. Cesaro T. Moratorio G. Hooikas P. J. Yakoleva A. Werneke S. W. (2016). Inhibition of polyamine biosynthesis is a broad-spectrum strategy against RNA viruses. J. Virol. 90 9683–9692. 10.1128/JVI.01347-16 Murataeva N. Dhopeshwarkard A. Yin D. Mitjavila J. Bradshaw H. Straiker A. (2016). Where’s my entourage? The curious case of 2-oleoylglycerol, 2-linolenoylglycerol, and 2-palmitoylglycerol. Pharmacol. Res. 110 173–180. 10.1016/j.phrs.2016.04.015 Naouar M. S. Mekki L. Z. Charfi L. Boubaker J. Filali A. (2016). Preventive and curative effect of Pistacia lentiscus oil in experimental colitis. Biomed. Pharmacother. 83 577–583. 10.1016/j.biopha.2016.07.021 Neis E. P. Dejong C. H. Rensen S. S. (2015). The role of microbial amino acid metabolism in host metabolism. Nutrients 7 2930–2946. 10.3390/nu7042930 Nie Y. F. Hu J. Yan X. H. (2015). Cross-talk between bile acids and intestinal microbiota in host metabolism and health. J. Zhejiang Univ. Sci. B 16 436–446. 10.1631/jzus.B1400327 Nielsch A. S. Ward F. W. Coates M. E. Walker R. Rowland I. R. (1991). Influence of dietary protein and gut microflora on endogenous synthesis of nitrate induced by bacterial endotoxin in the rat. Food Chem. Toxicol. 29 387–390. 10.1016/0278-6915(91)90078-L Nilsson Å (2016). Role of sphingolipids in infant gut health and immunity. J. Pediatr. 173 S53–S59. 10.1016/j.jpeds.2016.02.076 Olaya Galan N. N. Ulloa Rubiano J. C. Velez Reyes F. A. Fernandez Duarte K. P. Salas Cardenas S. P. Gutierrez Fernandez M. F. (2016). In vitro antiviral activity of Lactobacillus casei and Bifidobacterium adolescentis against rotavirus infection monitored by NSP4 protein production. J. Appl. Microbiol. 120 1041–1051. 10.1111/jam.13069 Pandey G. Gupta S. S. Bhatia A. Sidhu O. P. Rawat A. K. Rao C. V. (2016). Grilling enhances antidiarrheal activity of Terminalia bellirica Roxb. fruits. J. Ethnopharmacol. 10.1016/j.jep.2016.12.003 [Epub ahead of print]. Patra F. C. Sack D. A. Islam A. Alam A. N. Mazumber R. N. (1989). Oral rehydration formula containing alanine and glucose for the treatment of diarrhoea: a controlled trial. Br. Med. J. 298 1353–1356. 10.1136/bmj.298.6684.1353 Patten G. S. Abeywardena M. Y. Mcmurchie E. J. Jahangiri A. (2002). Dietary fish oil increases acetylcholine- and eicosanoid-induced contractility of isolated rat ileum. J. Nutr. 132 2506–2513. Pearl D. S. Masoodi M. Eiden M. Brummer J. Gullick D. Mckeever T. M. (2014). Altered colonic mucosal availability of n-3 and n-6 polyunsaturated fatty acids in ulcerative colitis and the relationship to disease activity. J. Crohns Colitis 8 70–79. 10.1016/j.crohns.2013.03.013 Peleman C. Camilleri M. Busciglio I. Burton D. Donato L. Zinsmesiter A. R. (2016). Colonic transit and bile acid synthesis or excretion in patients with irritable bowel syndrome-diarrhea without bile acid malabsorption. Clin. Gastroenterol. Hepatol. 10.1016/j.cgh.2016.11.012 [Epub ahead of print]. Pessione E. Cirrincione S. (2016). Bioactive molecules released in food by lactic acid bacteria: encrypted peptides and biogenic amines. Front. Microbiol. 7:876. 10.3389/fmicb.2016.00876 Plaza-Zamora J. Sabater-Molina M. Rodríguez-Palmero M. Rivero M. Bosch V. Nadal J. M. (2013). Polyamines in human breast milk for preterm and term infants. Br. J. Nutr. 110 524–528. 10.1017/S0007114512005284 Rabbani G. H. Sack D. A. Ahmed S. Peterson J. W. Saha S. K. Marni F. (2005). Antidiarrheal effects of L-histidine-supplemented rice-based oral rehydration solution in the treatment of male adults with severe cholera in Bangladesh: a double-blind, randomized trial. J. Infect. Dis. 191 1507–1514. 10.1086/428449 Ren W. Wu M. Luo W. Huang R. Yin Y. Li Y. (2013). Dietary supplementation with proline confers a positive effect in both porcine circovirus-infected pregnant and non-pregnant mice. Br. J. Nutr. 110 1492–1499. 10.1017/S0007114513000652 Rezaei R. Knabe D. A. Tekwe C. D. Dahanayaka S. Ficken M. D. Fielder S. E. (2013). Dietary supplementation with monosodium glutamate is safe and improves growth performance in postweaning pigs. Amino Acids 44 911–923. 10.1007/s00726-012-1420-x Rogers A. C. Mcdermott F. D. Mohan H. M. O’connell P. R. Winter D. C. Baird A. W. (2015). The effects of polyamines on human colonic mucosal function. Eur. J. Pharmacol. 764 157–163. 10.1016/j.ejphar.2015.07.006 Ruth M. R. Field C. J. (2013). The immune modifying effects of amino acids on gut-associated lymphoid tissue. J. Anim. Sci. Biotechnol. 4:27. 10.1186/2049-1891-4-27 Sabbatini P. Filipponi P. Sardella R. Natalini B. Nuti R. Macchiarulo A. (2013). Synthesis and quantitative structure-property relationships of side chain-modified hyodeoxycholic acid derivatives. Molecules 18:10497. 10.3390/molecules180910497 Schirmer M. Smeekens S. P. Vlamakis H. Jaeger M. Oosting M. Franzosa E. A. (2016). Linking the human gut microbiome to inflammatory cytokine production capacity. Cell 167 1125–1136.e8. 10.1016/j.cell.2016.10.020 Shimizu M. Zhao Z. Ishimoto Y. Satsu H. (2009). Dietary taurine attenuates dextran sulfate sodium (DSS)-induced experimental colitis in mice. Adv. Exp. Med. Biol. 643 265–271. 10.1007/978-0-387-75681-3_27 Sinha R. Ahn J. Sampson J. N. Shi J. Yu G. Xiong X. (2016). Fecal microbiota, fecal metabolome, and colorectal cancer interrelations. PLoS ONE 11:e0152126. 10.1371/journal.pone.0152126 Smolinska S. Jutel M. Crameri R. O’mahony L. (2014). Histamine and gut mucosal immune regulation. Allergy 69 273–281. 10.1111/all.12330 Soderholm J. D. Olaison G. Peterson K. H. Franzen L. E. Lindmark T. Wiren M. (2002). Augmented increase in tight junction permeability by luminal stimuli in the non-inflamed ileum of Crohn’s disease. Gut 50 307–313. 10.1136/gut.50.3.307 Somekawa S. Hayashi N. Nijima A. Uneyama H. Torii K. (2012). Dietary free glutamate prevents diarrhoea during intra-gastric tube feeding in a rat model. Br. J. Nutr. 107 20–23. 10.1017/S0007114511002674 Srinivas S. R. Prasad P. D. Umapathy N. S. Ganapathy V. Shekhawat P. S. (2007). Transport of butyryl- carnitine, a potential prodrug, via the carnitine transporter OCTN2 and the amino acid transporter ATB. Am. J. Physiol. Gastrointest. Liver Physiol. 293 G1046–G1053. Sugahara H. Odamaki T. Fukuda S. Kato T. Xiao J. Z. Abe F. (2015). Probiotic Bifidobacterium longum alters gut luminal metabolism through modification of the gut microbial community. Sci. Rep. 5:13548. 10.1038/srep13548 Swami U. Goel S. Mani S. (2013). Therapeutic targeting of CPT-11 induced diarrhea: a case for prophylaxis. Curr. Drug Targets 14 777–797. 10.2174/1389450111314070007 Taschler U. Hansenoehrl C. Storr M. Schicho R. (2017). Cannabinoid receptors in regulating the GI tract: experimental evidence and therapeutic relevance. Handb. Exp. Pharmacol. 10.1007/164-2016-105 [Epub ahead of print]. Tazume S. Takeshi K. Saidi S. M. Ichoroh C. G. Mutua W. R. Waiyaki P. G. (1990). Ecological studies on intestinal microbial flora of Kenyan children with diarrhoea. J. Trop. Med. Hyg. 93 215–221. Theriot C. M. Bowman A. A. Young V. B. (2016). Antibiotic-induced alterations of the gut microbiota alter secondary bile acid production and allow for clostridium difficile spore germination and outgrowth in the large intestine. mSphere 1:e00045–15. 10.1128/mSphere.00045-15 Timmons J. Chang E. T. Wang J. Y. Rao J. N. (2012). Polyamines and gut mucosal homeostasis. J. Gastrointest. Dig. Syst. 2 001. Trabelsi I. Ktari N. Ben Slima S. Hamden K. Ben Salah R. (2016). Effect of a probiotic Lactobacillus plantarum TN8 strain on trinitrobenzene sulphonic acid-induced colitis in rats. J. Anim. Physiol. Anim. Nutr. (Berl) 101 311–319. 10.1111/jpn.12536 Turcotte C. Chouinard F. Lefebvre J. S. Flamand N. (2015). Regulation of inflammation by cannabinoids, the endocannabinoids 2-arachidonoyl-glycerol and arachidonoyl-ethanolamide, and their metabolites. J. Leukoc. Biol. 97 1049–1070. 10.1189/jlb.3RU0115-021R Vlasova A. N. Kandasamy S. Chattha K. S. Rajashekara G. Saif L. J. (2016). Comparison of probiotic lactobacilli and bifidobacteria effects, immune responses and rotavirus vaccines and infection in different host species. Vet. Immunol. Immunopathol. 172 72–84. 10.1016/j.vetimm.2016.01.003 Wang A. Ling Z. Yang Z. Kiela P. R. Wang T. Wang C. (2015). Gut microbial dysbiosis may predict diarrhea and fatigue in patients undergoing pelvic cancer radiotherapy: a pilot study. PLoS ONE 10:e0126312. 10.1371/journal.pone.0126312 Wang H. Gao K. Wen K. Allen I. C. Li G. Zhang W. (2016). Lactobacillus rhamnosus GG modulates innate signaling pathway and cytokine responses to rotavirus vaccine in intestinal mononuclear cells of gnotobiotic pigs transplanted with human gut microbiota. BMC Microbiol. 16:109. 10.1186/s12866-016-0727-2 Wee L. H. Ueno K. Takeda H. Soma K. Kitagawa H. Satoh T. (1988). Ethyl eicosapentaenoate restored the immunosuppression in mice fed fat-free diet. Res. Commun. Chem. Pathol. Pharmacol. 62 49–66. Wen K. Liu F. Li G. Bai M. Kocher J. Yang X. (2015). Lactobacillus rhamnosus GG dosage affects the adjuvanticity and protection against rotavirus diarrhea in gnotobiotic pigs. J. Pediatr. Gastroenterol. Nutr. 60 834–843. 10.1097/MPG.0000000000000694 Yang X. Twitchell E. Li G. Wen K. Weiss M. Kocher J. (2015). High protective efficacy of rice bran against human rotavirus diarrhea via enhancing probiotic growth, gut barrier function, and innate immunity. Sci. Rep. 5:15004. 10.1038/srep15004 Yeh K. Y. Yeh M. Holt P. R. (1989). Differential effects of thyroxine and cortisone on jejunal sucrase expression in suckling rats. Am. J. Physiol. 256 G604–G612. Yen C. Tate J. E. Hyde T. B. Cortese M. M. Lopman B. A. Jiang B. (2014). Rotavirus vaccines: current status and future considerations. Hum. Vaccin. Immunother. 10 1436–1448. 10.4161/hv.28857 Yuan D. Wu Z. Wang Y. (2016). Evolution of the diacylglycerol lipases. Prog. Lipid Res. 64 85–97. 10.1016/j.plipres.2016.08.004 Zavala-Mendoza D. Alarcon-Aguilar F. J. Perez- Gutierrez S. Escobar-Villanueva M. C. Zavala-Sanchez M. A. (2013). Composition and antidiarrheal activity of bidens odorata cav. Evid. Based Complement. Alternat. Med. 2013:170290. 10.1155/2013/170290 Zhang H. Wang H. Shepherd M. Wen K. Li G. Yang X. (2014). Probiotics and virulent human rotavirus modulate the transplanted human gut microbiota in gnotobiotic pigs. Gut Pathog. 6 39. 10.1186/s13099-014-0039-8 Abbreviations EcN

Escherichia coli Nissle 1917

 HPLC-MS/MS

high-performance liquid chromatography coupled with tandem mass spectrometry

 HRV

human rotavirus

 LGG

Lactobacillus rhamnosus GG

 LIC

large intestinal contents

 Pro

probiotic

 Pro+RB

combination treatment of probiotic + rice bran

 PND

postnatal day

 RB

rice bran

 UPLC-MS/MS

ultra-high performance liquid chromatography coupled with tandem mass spectrometry