microorganisms Article Encapsulated Mixture of Methyl Salicylate and Tributyrin Modulates Intestinal Microbiota and Improves Growth Performance of Weaned Piglets Yusen Wei 1 , Jiangdi Mao 1 , Jingliang Liu 1 , Yu Zhang 1 , Zhaoxi Deng 1 , Jiaqi Lv 2 , Maolong He 2 , Jianxin Liu 1 and Haifeng Wang 1, * 1 The Key Laboratory of Molecular Animal Nutrition, Ministry of Education, College of Animal Science, Zhejiang University, Hangzhou 310058, China; yusenwei@zju.edu.cn (Y.W.); jiangdimao@126.com (J.M.); jingliangliu@zju.edu.cn (J.L.); zhangyuu@zju.edu.cn (Y.Z.); zhaoxideng@zju.edu.cn (Z.D.); liujx@zju.edu.cn (J.L.) 2 Lucta SA., UAB Research Park, Edifici Eureka, 08193 Bellaterra, Spain; louis.lv@lucta.com (J.L.); maolong.he@lucta.com (M.H.) * Correspondence: haifengwang@zju.edu.cn; Fax: +86-571-88982251 Abstract: Tributyrin and essential oils have been used as alternatives to antimicrobials to improve gut health and growth performance in piglets. This study was to evaluate the effects of a dietary supplement with two encapsulated products containing different combinations of tributyrin with oregano or with methyl salicylate on growth performance, serum biochemical parameters related  to the physiological status, intestinal microbiota and metabolites of piglets. A total of 108 weaned  crossbred piglets (Yorkshire × Landrace, 21 ± 1 d, 8.21 ± 0.04 kg) were randomly divided into three Citation: Wei, Y.; Mao, J.; Liu, J.; groups. Piglets were fed with one of the following diets for 5 weeks: a basal diet as the control Zhang, Y.; Deng, Z.; Lv, J.; He, M.; Liu, (CON); the control diet supplemented with an encapsulated mixture containing 30% of methyl J.; Wang, H. Encapsulated Mixture of salicylate and tributyrin at a dosage of 3 kg/t (CMT); and the control diet supplemented with an Methyl Salicylate and Tributyrin encapsulated mixture containing 30% of oregano oil and tributyrin at a dosage of 3 kg/t (COT). At Modulates Intestinal Microbiota and the end of the feeding trial, six piglets from each group were slaughtered to collect blood and gut Improves Growth Performance of samples for physiological status and gut microbiological analysis. The study found that the CMT Weaned Piglets. Microorganisms 2021, group was larger in feed intake (FI) (p < 0.05), average daily gain (ADG) (p = 0.09), total protein 9, 1342. https://doi.org/10.3390/ (TP), albumin (ALB), glutathione peroxidase (GSH-PX) (p < 0.05), blood total antioxidant capacity microorganisms9061342 (T-AOC) (p < 0.05), and crypt depth in the ileum (p < 0.05) compared with the CON group. The genus Academic Editor: Michael H. Kogut abundance of Tissierella and Campylobacter in the CMT group was significantly decreased compared with the CON group. The CMT group also resulted in significantly higher activity in amino acid Received: 31 May 2021 metabolism and arginine biosynthesis compared with the CON group. The COT group was larger in Accepted: 17 June 2021 T-AOC, and the genus abundance of Streptophyta and Chlamydia was significantly increased in the Published: 21 June 2021 ileum compared with the CON group. Data analysis found a significantly high correlation between the genus abundance of Chlamydia and that of Campylobacter in the ileum. The genus abundance of Publisher’s Note: MDPI stays neutral Campylobacter was also positively correlated with the sorbitol level. In general, the results indicated with regard to jurisdictional claims in that the supplementation of both encapsulated mixtures in diet of weaned piglets could improve published maps and institutional affil- the animal blood antioxidant capacity. Additionally, the encapsulated mixture of methyl salicylate iations. plus tributyrin improved the growth performance and resulted in certain corresponding changes in nutrient metabolism and in the genus abundance of ileum microbial community. Keywords: tributyrin; oregano essential oil; methyl salicylate; growth performance; microbiota; Copyright: © 2021 by the authors. metabolomics; piglets Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons 1. Introduction Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ There were growing evidence that the stress response caused by weaning would 4.0/). destroy the anti-oxidative balance [1], caused intestinal inflammation and damaged the Microorganisms 2021, 9, 1342. https://doi.org/10.3390/microorganisms9061342 https://www.mdpi.com/journal/microorganisms Microorganisms 2021, 9, 1342 2 of 17 intestinal barrier, furthermore worsened the absorptive capacity of the small intestine [2], which would decrease growth performance and eventually cause economic losses [3]. The compatibility of organic acids and plant essential oils to mitigate the negative effects of weaning stress at a lower cost had attracted much attention [4]. Butyric acid as a short-chain fatty acid had a certain therapeutic effect on intestinal diseases by improving intestinal morphology and reducing the rate of apoptosis [5]. Tributyrin, as a source of bu- tyrate, that one molecule releases three molecules of butyrate directly in the small intestine, could be used as an effective feed additive for weaned piglets to improve performance and intestinal health of piglets [6]. Compared with other butyrates, tributyrin has acceptable organoleptic characteristics and provides a higher and slowly released source of butyric acid in the small intestine [7]. Oregano oil was derived from the natural plant oregano (Origanumvuelgare L.) [8]. The main components of the oregano oil extracted from leaves and flowers were pheno- lic compounds, in which carvacrol and thymol had great antibacterial and antioxidant properties [9]. Dietary inclusion of oregano oil has many benefits for the weaned piglets, including strengthening immune response, improving antioxidant capacity and relieving stress [10]. Methyl salicylate was also a kind of essential oil extracted from holly (I. purpurea Hassk.), which has been used for its bactericidal and anti-inflammatory effects [11]. There are few reports on the supplementation of methyl salicylate in animal diets. Additionally, considerable attention has been paid to the potential use of organic acids and plant essential oils to look for benefits in the animal. Based on our previous study, tributyrin at dosage of 1.4 g per kilogram combined with oregano essential oil (OT) or methyl salicylate (MT) respectively 0.6 g per kilogram in powered form may have beneficial effects at the intestinal level of piglets, including the morphological structure and microbiota communities [12]. However, both treatments did not show significant differences in performance parameters. Here, we introduced the encapsulated mixture of essential oil and tributyrin to avoid the lipophilicity and volatility of plant essential oil. This study was aimed to investigate the role of dietary inclusion of encapsulated mixture of tributyrin and essential oil or methyl salicylate in a granular form to modulate the gut microbiota, metabolomics and intestinal morphology, and improving serum antioxidant capacity and growth performance of weaned piglets. 2. Materials and Methods 2.1. Animals and Experimental Diets A total of 108 weaned piglets (Yorkshire × Landrace, 21 ± 1 d, body weight 8.21 ± 0.04 kg), were allowed to three groups according to sex and weight, with 6 pens for each group and 6 piglets for each pen. The piglets were fed with a base diet as the control (CON), an antibiotic-free basal diet composed based on NRC weaned piglet nutritional standards, or the control supplemented with an encapsulated mixture containing 30% of methyl salicylate and tributyrin at a dosage of 3 kg/t (CMT); and the control diet supplemented with an encapsulated mixture containing 30% oregano oil and tributyrin at a dosage of 3 kg/t (COT). Encapsulation was made mainly with calcium palmitate and the sample products were provided by Lucta (Guangzhou) Flavours Co., Ltd. (Guangdong, China). All animal procedures were performed in full accordance with the Regulation for the Use of Experimental Animals in Zhejiang Province, China. This work was specifically approved by the Animal Care and Use Committee of Zhejiang University (ethics code permit no. ZJU20170529). The basic diet composition and nutritional level are shown in Table 1. The experiment lasted for 35 days. The body weight was individually measured on day 0 and day 35. Feed consumption was recorded every 3 days until the end of the experimental period. At day 35, six piglets were selected from each group and euthanized by intra- abdominal injection of 200 mg/kg pentobarbital sodium to collect biological samples. Microorganisms 2021, 9, 1342 3 of 17 Table 1. Basal diet ingredients and chemical composition. Items, g/kg Ingredients Corn 597 Soybean meal 105 Fermented soybean meal 55 Extruded soybean 100 Fish meal 40 Whey powder 30 Glucose 20 Wheat flour 4.4 Soybean oil 10 L-Lysine hydrochloride 3 L-Threonine 1.5 DL-Methionine 0.6 Choline chloride 1 Sodium chloride 3.5 Calcium hydrophosphate 9 Limestone 10 Vitamin-mineral premix 1 10 Nutrition composition Digestible energy † , MJ/kg (calculated) 143.96 Crude protein 190.05 Calcium 10.15 Lysine 13.37 Methionine 4.11 1 Provided per kilogram of diet: retinyl acetate, 1.2 MIU; cholecalciferol, 2700 IU; rac-α-tocopheryl acetate, 75 mg; menadione, 1.25 mg; thiamin, 1.5 mg; riboflavin, 2.5 mg; pantothenic acid, 40 mg; niacin, 34.4 mg; pyridoxol, 2.5 mg; biotin, 0.3 mg; folic acid, 3 mg; cobyrinic acid, 0.04 mg; Zn, 250 mg; Fe, 110 mg; Cu, 100 mg; Mn, 50 mg; I, 0.5 mg; Se, 0.5 mg. † Digestible energy was calculated. 2.2. Collection of Blood and Intestinal Samples Blood samples were collected using a coagulation tube, centrifuged at 3000× g and 4 ◦ C for 15 min. Serum was collected and stored at −20 ◦ C. Ileum and colonic tissues and contents were collected, and immediately frozen in liquid nitrogen and stored at −80 ◦ C. Ileum and colon tissue were also fixed with phosphate-buffered paraformaldehyde (4%, pH 7.6) for further histological measurements. The pH value of the ileum and colon contents was determined immediately. 2.3. Biochemical Analysis and Anti-Oxidatant Capacity Related Parameters of Serum Concentration of serum albumin (ALB), globulin (GLB), albumin/globulin(A/G), glutamic-pyruvic transaminase (ALT), glutamic-oxalacetic transaminase (AST), glucose (GLU), lactate dehydrogenase (LDH), total cholesterol (TC), triglycerides (TG), total protein (TP), uric acid (UA) and blood urea nitrogen (BUN) were determined by automatic blood cell analyzer (XN-2000-A1, SYSMEX, Kobe, Japan). Anti-oxidant capacity related blood parameters included methylene dioxyamphetamine (MDA), superoxide dismutase (SOD), glutathione peroxidase (GSH-PX) and blood total antioxidant capacity (T-AOC), which were measured using commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer’s instructions). 2.4. Morphological Analysis of Ileum The ileum tissue of the piglet was immersed in 10% neutral formalin and covered with wax. The waxed tissue block was manually cut into 3 µm-thick sections, and subjected to deparaffinization and dehydration. Then these sections were treated with different concentrations of alcohol (100%, 95% and 75%) for 15 min, and stained with H&E [13]. An optical microscope system (Olympus Corporation, Tokyo, Japan) was used to obtain a micrograph at a combined magnification of 100× using ImageJ software (National Institute Microorganisms 2021, 9, 1342 4 of 17 of Health, Bethesda, MD, USA). Villus height and crypt depth were measured according to a described method [3]. The ratio of villus height to crypt depth (VH/CD) was calculated. 2.5. Analysis of Short-Chain Fatty Acid in Colonic Content A slightly modified method was used to detect short-chain fatty acids (SCFAs). Colonic contents (0.4 g) were mixed and vortexed with 1.5 mL of phosphate-buffered saline (PBS). The sample was centrifuged at 15,000× g for 15 min at 4 ◦ C. The supernatant was added with 25% metaphosphoric acid at a ratio of 9:1 (v/v), and then the supernatant was passed through a 0.22 µm filter membrane to prepare a colon sample for short-chain fatty acid (SCFA) analysis. A gas chromatograph (model: GC-2010; Shimadzu Corp, Kyoto, Japan) equipped with a chromatographic column (HP-INNOWAX (19091N-133), 30 m × 0.25 mm × 0.25 µm) was used to determine the VFA concentration in the sample [14]. 2.6. 16S rRNA Sequencing Analysis Microbial DNA was extracted from ileum contents using the QIAamp Fast DNA Stool Minikit (Qiagen, Hilden, Germany) according to manufacturer’s protocols. The V3-V4 region of the bacteria 16S rRNA genes were amplified by PCR (95 ◦ C for 3 min, followed by 30 cycles at 98 ◦ C for 20 s, 58 ◦ C for 15 s, and 72 ◦ C for 20 s and a final extension at 72 ◦ C for 5 min) using primers 341F 5’-CCTACGGGRSGCAGCAG)-3’ and 806R 5’- GGACTACVVGGGTATCTAATC-30 . PCR reactions were performed in 30 µL mixture containing 15 µL of 2 × KAPA Library Amplification ReadyMix, 1 µL of each primer (10 µM), 50 ng of template DNA and ddH2 O. Amplicons were extracted from 2% agarose gels and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) according to the manufacturer’s instructions and quantified using Qubit® 2.0 (Invitrogen, Carlsbad, CA, USA). After preparation of library, these tags were sequenced on MiSeq/HiSeq platform (Illumina, Inc., San Diego, CA, USA) for paired-end reads of 500/250 bp. DNA extraction, library construction and sequencing were conducted at Realbio Genomics Institute (Shanghai, China). Tags, trimmed of barcodes and primers, were further checked on their rest lengths and average base quality. 16S tags were restricted between 220 bp and 500 bp such that the average Phred score of bases was no worse than 20 (Q20) and no more than 3 ambiguous N. The copy number of tags was enumerated and redundancy of repeated tags was removed. Only the tags with a frequency of more than 1, which tend to be more reliable, were clustered into OTUs, each of which had a representative tag. Operational taxonomic units (OTUs) were clustered with 97% similarity using UPARSE (http://drive5.com/ uparse/; accessed on 11 November 2020) [15] and chimeric sequences were identified and removed using Userach (version 7.0). Each representative tags was assigned to a taxa by RDP Classifer (http://rdp.cme.msu.edu/; accessed on 11 November 2020) [16,17] against the RDP database (http://rdp.cme.msu.edu/; accessed on 11 November 2020) using a confidence threshold of 0.8. OTU profiling table and alpha/beta diversity analyses were also achieved by python scripts of QIIME (V1.9.1). Anosim was analyzed for the microbial community structure, and Unifrac algorithm was used to calculate the distance between the two samples. Heatmap was created using R (V3.5.1) “g plots” package. Spearman correlation heat map was drawn by the R (V3.5.1) “corrplot” package. PICRUSt (phylogenetic investigation of communities by reconstruction of unobserved states) is based on 16S rRNA and reference sequence database to predict the functional composition of metagenomics. 2.7. GC-TOF-MS and Metabolomics Data Processing After extracting metabolites [18], all ileum contents were analyzed by GC-TOF-MS analysis performed using an Agilent 7890 gas chromatograph coupled with a time-of-flight mass spectrometer (J&W Scientific, Folsom, CA, USA). The system was equipped with a DB-5MS capillary column. One microliter aliquot of sample was injected in splitless mode. Helium was used as the carrier gas, the front inlet purge flow was 3mL min−1 , and the Microorganisms 2021, 9, 1342 5 of 17 gas flow rate through the column was 1 mL min−1 . The initial temperature was kept at 50 ◦ C for 1min, then raised to 310 ◦ C at a rate of 10 ◦ C min−1 , then kept for 8 min at 310 ◦ C. The injection, transfer line, and ion source temperatures were 280 ◦ C, 280 ◦ C and 250 ◦ C, respectively. The energy was −70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with the m/z range of 50–500 at a rate of 12.5 spectra per second after a solvent delay of 6.35 min. Raw data analysis, including peak extraction, baseline adjustment, deconvolution, alignment and integration, was finished with Chroma TOF (V 4.3x, LECO) software and LECO-Fiehn Rtx5 database was used for metabolite identification by matching the mass spectrum and retention index [18]. Finally, the peaks detected in less than half of QC samples or RSD > 30% in QC samples were removed. Then, the missing values were filled up by half of the minimum value. Internal standard normalization method was employed in this data analysis. Data were scaled and logarithmically transformed to minimize the impact of both noise and high variance of the variables. The final dataset containing the information of peak number, sample name and normalized peak area was imported to the MetaboAnalyst 4.0 (http://www.metaboanalyst.ca/; accessed on 10 January 2021) for multivariate analysis. Orthogonal partial least squares discrimination analysis (OPLS-DA), and t-test were performed between the two groups, with the FDR adjusted p value < 0.05 and the VIP > 1.5 being considered as significantly different metabolites. In addition, commercial databases including KEGG (http://www.genome.jp/kegg/; accessed on 10 January 2021) were used for pathway enrichment analysis. 2.8. Statistical Analysis Sample sizes were chosen empirically based on previous experiments. The piglet was used as the experimental unit for blood, intestinal and microbial parameters, and data from the pen for performance parameters. Data was analyzed by using SPSS 20 (SPSS Inc., Chicago, IL, USA). When the data satisfies the normal distribution and homogeneity of variance, one-way of variance (ANOVA) was conducted and means were further compared by using Tukey’s multiple comparison tests. The non-parametric Kruskal–Wallis test with Dunn’s multiple comparisons test was used for a date that was not normally distributed. Results are reported as means ± standard error of the mean (SEM). p < 0.05 was considered as significant indicated by “*” and 0.05 < p < 0.1 was considered as a tendency. 3. Results 3.1. Growth Performance The average body weight of piglets increased from 8.21, 8.10 and 8.31 kg at initial for the CON, COT and CMT groups (Figure 1A), respectively, to 17.64, 17.74 and 19.24 kg in the end (Figure 1B). The CMT group was larger in the FI compared with the CON group (p < 0.05) and the COT group (p < 0.05), which increased by 13.9% and 11.6% respectively (Figure 1C). These correspondingly tended to have larger ADG compared with the CON group (p = 0.09), which increased by 15.9% (Figure 1D). Compared with the CON group, the COT group only slightly increased the FI and the ADG by 2.1% which were not significant. The difference on feed to gain ratio among the groups was not significant, although the CMT group slightly reduced it by 1.5% whereas the COT group increased it by 1.7% compared with the CON group (Figure 1E). Microorganisms 2021, 9, x FOR PEER REVIEW 6 of 18 nificant. The difference on feed to gain ratio among the groups was not significant, alt- Microorganisms 2021, 9, 1342 hough the CMT group slightly reduced it by 1.5% whereas the COT group increased6itofby 17 1.7% compared with the CON group (Figure 1E). Figure 1. Growth performance of piglets in different treatments. (A) Initial body weight. (B) Final body weight. (C) Feed intake. (D) (D) Average Average daily gain. (E) (E) Feed Feed to to gain gain ratio. ratio. Star Star (*) (*) indicates indicates that that there there are are significantly significantly different different (at (at the the 5% 5% threshold of the ANOVA) between groups. threshold of the ANOVA) between groups. 3.2. 3.2. Serum Serum Biochemical Biochemical Indices Indices and and Anti-Oxidative Anti-Oxidative Capacity Capacity Parameters Parameters There There were no significant differences in A/G,LDH, were no significant differences in A/G, LDH,TC, TC,GLU, GLU,ALT, ALT, AST, AST, BUN, BUN, UA UA among the three groups (Figure S1A–H). The content of TG tended to increase among the three groups (Figure S1A–H). The content of TG tended to increase in the CMT in the CMT group comparedwith group compared withthetheCON CON(p(p = 0.09, = 0.09, Figure Figure 2A).2A). ALB ALBandand TP content TP content werewere signifi- significantly cantly higher in the CMT group than the CON group (p < 0.05, Figure higher in the CMT group than the CON group (p < 0.05, Figure 2B,C). MDA level had 2B,C). MDA levela had decreased tendency in the CMT group compared with the CON group (p = 0.08, Figure0.08, a decreased tendency in the CMT group compared with the CON group (p = 2D). Figure 2D).no There was There was no difference significant significantin difference SOD among in SOD among different different treatments treatments (Figure S1I). (Figure Com- Microorganisms 2021, 9, x FOR PEERS1I). pared Compared REVIEW with the CONwith the CONthe group, group, CMTthe CMThad group group had significantly significantly higher GSH-PX higher GSH-PX of(p18 (p < 70.05,
0, Figure 5B,C). The5D). The relativeofabundance visualization beta diversityof bacterial showedtaxa thatatthe theexperi- phylum (Figure mental 5E) andisgenus treatment level (Figure meaningful 5F) in piglets (R > 0, Figure 5D). The from different relative groupsof abundance were shown. bacterial The taxa phyla Firmicutes showed a tendency to increase in abundance in the at the phylum (Figure 5E) and genus level (Figure 5F) in piglets from different groups experimental groups compared were shown. The withphyla the CON groupshowed Firmicutes (Figurea5E). The genus tendency of Clostridium to increase sensuinstricto in abundance also the ex- showed agroups perimental tendency to increase compared withinthe abundance CON group in the experimental (Figure 5E). The groups genus of compared Clostridium with the stricto sensu CON group (Figurea5F). also showed The ratio tendency of Firmicutes/Bacteroidetes to increase in abundance in the (F/B) tended to increase experimental groups in the CMTwith compared groupthecompared CON group with the CON5F). (Figure groupThe (Figure ratio 5G). The heatmap summarized of Firmicutes/Bacteroidetes (F/B)the differences of three groups of bacteria at the genus level namely tended to increase in the CMT group compared with the CON group (Figure 5G). The Streptophyta, Chlamydia, Tissierella heatmap and Campylobacter summarized (Figureof the differences 5H). Spearman three groups of correlation bacteria attest theshowed that namely genus level there was a positive correlation between Campylobacter and Chlamydia (Figure Streptophyta, Chlamydia, Tissierella and Campylobacter (Figure 5H). Spearman correlation5I). Microorganisms 2021, 9, x FOR PEER REVIEW 10 of 18 Microorganisms 2021, 9, 1342 9 of 17 test showed that there was a positive correlation between Campylobacter and Chlamydia (Figure 5I). Figure 5. The high-throughput sequencing of 16S rRNA in ileum contents of piglets. (A) Venn analysis Figure of bacterial 5. The OTUs composition high-throughput sequencingamong of 16StherRNA three in groups. ileum(B) Chao1ofindex. contents (C)(A) piglets. Simpson index. Venn anal- ysis (D) of bacterialanalyses ANOSIM, OTUs composition among of similarities amongthethethree threegroups. groups of(B)samples. Chao1 index. (C) Simpson (E) Proportional index. taxonomic (D) ANOSIM, at assignments analyses of similarities the phylum level fromamong the groups. the three three groups of samples.taxonomic (F) Proportional (E) Proportional taxo- assignments nomic assignments at the phylum level from the three groups. (F) Proportional taxonomic at the genus level from three groups. (G) The ratio of F/B among the three groups. (H) Heatmap assign- ments at the genus level from three groups. (G) The ratio of F/B among the three groups. (H) cluster analysis of differential bacteria at genus level which has a contribution to group differences. Heatmap cluster analysis of differential bacteria at genus level which has a contribution to group (I) Spearman correlation analysis of the differential bacteria. differences. (I) Spearman correlation analysis of the differential bacteria. A total of 190 Kos (Figure S2A) and 50 MetaCyc pathways (Figure S2B) were identified withAsignificant total of 190differences Kos (Figureamong S2A) and the 50 MetaCyc three pathways groups. (Figure Compared withS2B) the were CONidenti- group, fied the expressions of KO0059, KO0626, KO1692, KO0249 and other genes related group, with significant differences among the three groups. Compared with the CON to fatty the expressions acid metabolism of KO0059, KO0626, in the CMT KO1692, and the KO0249were COT groups and other genes related up-regulated to fatty (Figure S2C).acid The expression of quorum sensing and arginine biosynthesis-related genes such as KO2035, KO2034, KO2032, KO2033 and KO5597 were up-regulated in the CMT group compared with the CON group (Figure S2C). We found that the expression of toxin-related genes Microorganisms 2021, 9, 1342 10 of 17 such as K18839, K10953 and K10938 were up-regulated in the COT group compared with the CMT group (Figure S2C). 3.6. Overview of the Intestinal Metabolome In this study, 686 peaks were detected and 482 metabolites were left after relative standard deviation de-noising. Supervised orthogonal projections to latent structures- discriminate analysis (OPLS-DA) were applied (Figure 6A–C) to visualize group separation and find significantly changed metabolites. The OPLS-DA model among different groups were as follows: R2 X = 0.653, R2 Y = 0.824, and Q2 = 0.279 (Figure 6A). The OPLS-DA model between the CMT group and the CON group were as follows: R2 X = 0.67, R2 Y = 0.973 and Q2 = 0.632 (Figure 6B). The OPLS-DA model between the CMT group and the COT group were as follows: R2 X = 0.336, R2 Y = 0.918, and Q2 = 0.568 (Figure 6C). The Q2 values all exceeded 0.4, indicating that these models were more reliable and that consistent modeling and predictability were achieved. However, the OPLS-DA model between the COT group and the CON group were as follows: R2 X = 0.626, R2 Y = 0.646, and Q2 = −0.0417 (Figure S3A). Subsequently, we detected 25 differential metabolites among the three groups. There were 17, 23 and 11 differential metabolites in the CON, CMT and COT groups, respectively (Figure 6D,E, Figure S3B). Exact mass data (m/z) from the KEGG and HMDB databases were used to anno- tate the differential metabolites for pathway enrichment analysis (Figure 6F). Differential metabolic pathways between each two group were shown (Figure 6G,H). The significant differences in the expression of molecules were associated with galactose metabolism, arginine biosynthesis, aminoacyl-tRNA biosynthesis, alanine, aspartate and glutamate metabolism during the trial period. The level of methionine, serine, 2-ketobutyric acid related to cysterine and methionmine metabolism in the CMT group were significantly lower than in the CON group. Regarding arginine biosynthesis, the concentrations of aspartic acid and N-alpha-Acety-L-ornithnine in the CMT group were significantly lower than that in the CON group. Concerning aminoacyl-tRNA biosynthesis, aspartic acid in the CMT group was also less than that in the CON group. Compared with the CMT group, the COT group had a significantly higher concentration of glucose related to starch and sucrose metabolism and neomycin, kanamycin and gentamicin biosynthesis, and also had a significantly higher level of aspartic acid, maleic acid and 4-aminobutyric related to alanine, aspartate and glutamate metabolism. Microorganisms Microorganisms 2021, 2021, 9, 9, x1342 FOR PEER REVIEW 1211of of 18 17 Figure 6. Non-targeted metabonomics derived from the GC/MS metabolite profiles of ileum contents. (A) OPLS-DA score map Figure of separate comparisons 6. Non-targeted metabonomics among thefrom derived threethegroups. GC/MS (B) OPLS-DAprofiles metabolite score map of separate of ileum comparisons contents. (A) OPLS-DA between score map of separate the CMT groupcomparisons and the CON among the(C) group. three groups. (B) OPLS-DA OPLS-DA score map of score mapcomparisons separate of separate comparisons between thebetween the CMT CMT group and group the COTandgroup. the CON (D)group. (C) OPLS-DA Significantly scoremetabolites differential map of separate comparisons between the CON between group andthethe CMTCMTgroup and(E) group. theSignificantly COT group. (D) Significantly differential differential metabolites metabolites between the COT between group and thetheCON CMT group and group. (F)the CMT group. Differential (E) Significantly metabolites differential were classified me- by KEGG tabolites between the COT group and the CMT group. (F) Differential metabolites were classified by KEGG pathway pathway enrichment and significance analysis. (G) Pathway enrichment analysis performed using the significantly different enrichment and significance analysis. (G) Pathway enrichment analysis performed using the significantly different intes- intestinal metabolites between the CON group and the CMT group. (H) Pathway enrichment analysis performed using the tinal metabolites between the CON group and the CMT group. (H) Pathway enrichment analysis performed using the significantly different significantly different intestinal intestinal metabolites metabolites between between thethe COT COT group group and and the the CMT CMT group. group. Exact massofdata 3.7. Correlation (m/z) from Characteristic theMicrobiota Gut KEGG and andHMDB databases were used to annotate Metabolites the differential metabolites for pathway enrichment analysis (Figure 6F). Differential met- The correlation between differential gut bacteria and metabolites among the three abolic pathways between each two group were shown (Figure 6G,H). The significant dif- groups was expounded by Spearman correlation analysis (Figure 7). The genus Strepto- nificantly higher level of aspartic acid, maleic acid and 4-aminobutyric related to alanine, aspartate and glutamate metabolism. 3.7. Correlation of Characteristic Gut Microbiota and Metabolites Microorganisms 2021, 9, 1342 The correlation between differential gut bacteria and metabolites among the12three of 17 groups was expounded by Spearman correlation analysis (Figure 7). The genus Strepto- phyta was strongly positively correlated with 2-deoxytetronic acid, aspartic acid, melibi- ose, maleic acid, sucrose, neohesperidin, phenyl β-D-glucopyranoside hydrate, cytidine phyta was strongly monophosphate, positively correlated 1,5-anhydroglucitol, with 2-deoxytetronic acid, 3,6-anhydro-D-galactose, aspartic tagatose acid, meli- and biphenyl (p biose, maleic acid, sucrose, neohesperidin, phenyl β-D-glucopyranoside hydrate, cytidine < 0.05), and highly positively correlated with 24,25-dihydrolanosterol, ergosterol, lactu- monophosphate, lose and creatine (p1,5-anhydroglucitol, < 0.05). The genus 3,6-anhydro-D-galactose, Campylobacter was stronglytagatose andcorrelated positively biphenyl (p < 0.05), and highly positively correlated with 24,25-dihydrolanosterol, ergosterol, with Na, Na-dimethylhistamine and sorbitol (p < 0.05), and highly positively correlated lactu- lose and creatine (p < 0.05). The genus Campylobacter with 2-deoxytetronic acid and aspartic acid (p < 0.05). was strongly positively correlated with Na, Na-dimethylhistamine and sorbitol (p < 0.05), and highly positively correlated with 2-deoxytetronic acid and aspartic acid (p < 0.05). Figure 7. Figure Heatmap of 7. Heatmap of Spearman Spearman correlation correlation analysis analysis between between the the altered altered microbiota microbiota and and metabolites. metabolites. 4. Discussion 4. Discussion Plant essential Plant essential oil oil which which hadhad many many beneficial beneficial effects effects on on animal animal health health has has been been exca- exca- vated [19]. Tributyrin was a new way to provide butyric acid, which was vated [19]. Tributyrin was a new way to provide butyric acid, which was more efficient in more efficient in intestinal intestinal utilization utilization thanthan butyric butyric acidacid [20]. [20]. UsingUsing the the application application of coating of coating technol- technology ogy could could solvesolve the disadvantage the disadvantage that essential that essential oil could oil could not reach not reach the intestinal the intestinal tract tract effec- effectively [21]. However, little is known either on the potential antibacterial tively [21]. However, little is known either on the potential antibacterial mechanism be- mechanism between tween plant-derived plant-derived compounds compounds andand organic organic acids. acids. Methyl salicylate which has the special flavorofofholly Methyl salicylate which has the special flavor hollyleaves was leaves wasanan additive generally additive gener- ally recognized as safe (GRAS) [22]. We speculated that the additive might have aa similar recognized as safe (GRAS) [22]. We speculated that the additive might have similar effect as effect asaaflavoring flavoringagent agentoror sweetener sweetener to improve to improve the the palatability palatability and and inductivity inductivity of of feed feed or have an antioxidant capacity to help to maintain the quality of feed to promote or have an antioxidant capacity to help to maintain the quality of feed to promote FI. Con- FI. Consequently, the great growth performance in the CMT group was presented by sequently, the great growth performance in the CMT group was presented by alleviating alleviating the decrease of FI caused by weaning stress [23]. However, the COT group had no significant improvement on growth performance, which was consistent with an experiment of three diets with an oregano supplementation at 2 g, 4 g and 8 g per kg feed [24]. The beneficial effect of dietary phytobiotics could be influenced by its dosage and diet composition [25], so more tests are needed to verify the role of different dosage. In addition, there is a limitation that the solo encapsulation vehicle was not determined on the performance of piglets, which will be worth continuing investigation. Serum biochemical indices can reflect growth performance, nutrition and health status to a certain extent. TP is composed of ALB and GLB. ALB revealed a significant correlation with body nutritional status, including ADG and feed conversion [26]. In this study, TP and ALB in the CMT group with better growth performance were significantly higher than the CON group. Consisted with the previous conclusion, in pigs, nutrition, body growth and ALB synthesis were interdependent and ALB was considered one of the most important predictors of performance, especially ADG [6]. The activities of GSH-PX and T-AOC increased significantly in the CMT group than in the CON group. The COT group had significantly higher T-AOC content than the CON group, which might result from the anti-oxidant capacity of the phenolic hydroxyl compounds in plant essential oils. The result was consistent with the previous report of inclusion of plant polyphenols to the Microorganisms 2021, 9, 1342 13 of 17 diet [27]. T-AOC directly reflects the antioxidant capacity of the body, while enhanced SOD and GSH-PX would reduce the concentration of superoxide anion and peroxide in the body, as well as, lower concentration of serum MDA demonstrates the lower degree of oxidative damage took place in the body. Overall, diet supplemented with an encapsulated mixture of methyl salicylate and tributyrin could enhance the antioxidant capacity and accelerate the protein synthesis. Weaning stress would damage gastrointestinal structure [28]. Villus length and crypt depth were correlated with the absorption surface and mucosal secretion [29]. Our results confirmed that the structure of ileum in both the CMT and COT groups had been improved and there was a positive relationship between serum antioxidant parameters and intestinal morphology, which was consistent with the previous study that intestinal injury was significantly reduced by SOD pretreatment [30]. The reason why the CMT group had better intestinal morphology might be attributed to phenolic hydroxyl groups, which have been reported to be beneficial to the maintenance of the intestinal integrity of piglets [31]. Another study has pointed out that tributyrin could strengthen the intestinal barrier via promoting cell differentiation [32]. A previous study showed that SCFAs could affect appetite and energy balance and regulate intestinal morphology [33]. SCFAs molecules could be absorbed by colonic epithelial cells and converted into energy to maintain the intestinal health of animals [34]. We examined the colonic SCFAs and found that valeric acid and isovaleric acid in the CMT group increased compared with the CON group, which was consistent with the conclusion of the previous literature data on the base diet supplemented with 1.5 g/kg organic acid and 30 mg/kg essential oils [35]. However, it was somewhat surprising since total SCFAs in the COT group appeared a tendency to decrease, which is worth investigating. In terms of α-diversity, compared with the CON group, the Simpson index and Chao1 index showed opposite variation in the CMT group. The significance of this finding rests on the fact that there is still a great deal of controversy on the relationship between weight or obesity and intestinal microbial diversity [36]. However, the F/B ratio in the CMT group was increased, which was consistent with characteristics of obesity [37]. It is worth noting that the phyla Proteobacteria containing a wide variety of pathogens, including Escherichia coli, Salmonella and Helicobacter pylori had a relatively high abundance in the CON group. As reported by several studies, oregano essential oil could inhibit the activity of Escherichia colis [38]. In line with this, the COT group showed a decrease in the abundance of Proteobacteria and Escherichia/Shigella. However, the antibacterial mechanism responsible for methyl salicylate in the animals’ intestine requires further discussion. Interestingly, compared with the CON group at the genus level, we found that the abundance of Tissierella decreased in both the COT and CMT groups, but the abundance of Campylobacter decreased only in the CMT group. Compared with positive bacteria, the structure of Campylobacter had stronger resistance to plant essential oils and was more sensitive to organic acids [39]. These different responses might allow for the discrepancy of the COT group and the CMT group on growth performance of piglets in this experiment. Moreover, as described by Qin et al. the drug resistance of Campylobacter in the intestinal tract of livestock and poultry in Asia was generally strong [40]. And through the correlation analysis, the increase of Campylobacter in the COT group might be the cause of the increase in the number of Chlamydia. Thus we tentatively linked the reason to the increasing abundance of Chlamydia, which were associated with different pathologies [41]. One hypothesis was that Veillonellaceae may be used as an index to measure the bactericidal ability of antimicrobials [42]. In this study, the lower abundance of Veillonellaceae indicated that the CMT group had a higher antimicrobial capacity. Additionally, a study found that microbial composition may not be a major factor determining differences in feed efficiency [43]. The CMT group had the higher FI but similar to feed to gain ratio compared with the CON group, which might cause by improving intestinal integrity [44]. Intestinal microorganisms could produce bioactive molecules which affect host metabo- lism [45]. PICRUST analysis showed that the genes related to fatty acid metabolism were Microorganisms 2021, 9, 1342 14 of 17 significantly up-regulated in both the COT and CMT groups. Actually, butyric acid has been proved to improve the energy metabolism of mitochondrial cells through activating Adenosine 5‘-monophosphate (AMP)-activated protein kinase (AMPK), especially fatty acid oxidation [46]. These reports indicated that the supply of butyric acid in the form of tributyrin could alleviate the effect of weaning stress. Notably, in the CMT group, we also found that the quorum sensing of bacteria was interfered, which may reduce the drug resistance of bacteria by changing the ratio of Opp/Dpp about biofilm formation. Consistently, evidence confirmed that plant essential oils could destroy the quorum sensing of bacteria [47]. Whereas the information reflected by PICRUST was limited. Consequently, we applied GC-MS non-targeted metabolomics to explain the functional difference. Via the OPLS-DA model, we found that the metabolic level of the CMT group had significant difference compared with the CON group, while there was no significant difference in the metabolic level between the COT group and the CON group. The data showed that the concentrations of L-serine, L-aspartic acid and L-methionine in the CMT group changed significantly, which indicated bioavailability of amino acids and protein synthesis was improved. Pathway analysis showed that differential metabolites were enriched in cysteine and methionine metabolism, alanine, aspartic acid and glutamic acid metabolism, arginine biosynthesis and so on. Additionally, several reports pointed out that the synthesis of endogenous arginine related with growth performance and antioxidation of piglets [48,49]. Finally, we analyzed the correlation between differential bacteria and metabolites to reveal the interaction between microorganisms and hosts. Although there is too little information about Streptophyta, which was positively correlated with a variety of ma- jor metabolites. We also found that there was a high correlation between sorbitol and Campylobacter. The influence caused by a relative higher abundance of Streptophyta and Campylobacter on the host is worth being investigated more. 5. Conclusions In conclusion, dietary inclusion of an encapsulated mixture of methyl salicylate and tributyrin improved the growth performance, antioxidant capacity, and intestinal villus morphology, and modulated the microbiota and their metabolites. The results may hope- fully serve as useful information for the use of encapsulated mixtures of methyl salicylate and tributyrin as feed additives in piglets. Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/microorganisms9061342/s1, Figure S1: Serum indices of different groups. (A) A/G. (B) LDH. (C) TC. (D) GLU. (E) ALT. (F) AST. (G) BUN. (H) UA. (I) SOD. Figure S2: The differences of gene expression among different groups analyzed by LEfSe. (A) LEfSe analysis based on Genomes (KEGG) ortholog (KO) database. (B) LEfSe analysis based on the MetaCyc database. (C) The expression of the top 20 KO genes respectively in different groups. Figure S3: Non-targeted metabonomics derived from the GC/MS metabolite profiles of ileum contents. (A) OPLS-DA score map of separate comparisons between the COT group and the CON group. (B) Significantly differential metabolites between the COT group and the CON group. Author Contributions: H.W. and M.H. designed the study, analyzed the data and revised the article; Y.W., J.M., Y.Z., J.L. (Jingliang Liu), and Z.D. performed the experiments; J.L. (Jiaqi Lv) Supervision; Y.W. wrote the manuscript; J.L. (Jianxin Liu) revised the article. All authors have read and agreed to the published version of the manuscript. Funding: This study was supported by grants from the Natural Science Foundation of Zhejiang Province (Z19C170001), the Funds of Ten Thousand People Plan of China, the National Key Research and Development Program of China (2017YFD0500502) and the Innovation Research Fund from Lucta Group. Acknowledgments: We would thank Weidong Hua for help in the animal feeding trial and Gemma Tedó for help on the article revision. Microorganisms 2021, 9, 1342 15 of 17 Conflicts of Interest: The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results. Abbreviations A/G: albumin/globulin; ALT, alanine aminotransferase; AST, aspartate aminotransferase; ALB, albu- min; GLB, globulin; GLU, glucose; LDH, lactate dehydrogenase; TC, cholesterol; TG, triglyceride; TP, total protein; UA, uric acid; BUN, blood urea nitrogen; MDA, methylene dioxyamphetamine; SOD, superoxide dismutase; GSH-PX, glutathione peroxidase; T-AOC, blood total antioxidant capacity; VH/CD, villus height/crypt depth; SCFAs, short-chain fatty acids; ADG, average daily gain; FI, feed intake; OTUs, operational taxonomic units; F/B, Firmicutes/Bacteroidetes; g, genus. References 1. Cao, S.T.; Wang, C.C.; Wu, H.; Zhang, Q.H.; Jiao, L.F.; Hu, C.H. Weaning disrupts intestinal antioxidant status, impairs intestinal barrier and mitochondrial function, and triggers mitophagy in piglets. J. Anim. Sci. 2018, 96, 1073–1083. [CrossRef] 2. Moeser, A.J.; Pohl, C.S.; Rajput, M. Weaning stress and gastrointestinal barrier development: Implications for lifelong gut health in pigs. Anim. Nutr. 2017, 3, 313–321. [CrossRef] 3. Hou, Y.Q.; Liu, Y.L.; Hu, J.; Shen, W.H. Effects of lactitol and tributyrin on growth performance, small intestinal morphology and enzyme activity in weaned pigs. Asian Australas. J. Anim. Sci. 2006, 19, 1470–1477. [CrossRef] 4. Zhang, W.; Zhu, Y.H.; Yang, G.Y.; Liu, X.; Xia, B.; Hu, X.; Su, J.H.; Wang, J.F. Lactobacillus rhamnosus GG affects microbiota and suppresses autophagy in the intestines of pigs challenged with salmonella infantis. Front. Microbiol. 2017, 8, 2705. [CrossRef] 5. Diao, H.; Jiao, A.R.; Yu, B.; Mao, X.B.; Chen, D.W. Gastric infusion of short-chain fatty acids can improve intestinal barrier function in weaned piglets. Genes Nutr. 2019, 14, 4. [CrossRef] [PubMed] 6. Sotira, S.; Dell’Anno, M.; Caprarulo, V.; Hejna, M.; Pirrone, F.; Callegari, M.L.; Tucci, T.V.; Rossi, L. Effects of tributyrin supplementation on growth performance, insulin, blood metabolites and gut microbiota in weaned piglets. Animals 2020, 10, 726. [CrossRef] 7. Cresci, G.; Nagy, L.E.; Ganapathy, V. Lactobacillus GG and tributyrin supplementation reduce antibiotic-induced intestinal injury. J. Parenter. Enter. Nutr. 2013, 37, 763–774. [CrossRef] [PubMed] 8. Calo, J.R.; Crandall, P.G.; O’Bryan, C.A.; Ricke, S.C. Essential oils as antimicrobials in food systems—A review. Food Control. 2015, 54, 111–119. [CrossRef] 9. Michiels, J.; Missotten, J.; Dierick, N.; Fremaut, D.; Maene, P.; Smet, S.D. In vitro degradation and in vivo passage kinetics of carvacrol, thymol, eugenol and trans-cinnamaldehyde along the gastrointestinal tract of piglets. J. Sci. Food Agric. 2010, 88, 2371–2381. [CrossRef] 10. Zhang, T.; Zhou, Y.F.; Zou, Y.; Hu, X.M.; Zheng, L.F.; Wei, H.K.; Giannenas, I.; Jin, L.Z.; Peng, J.; Jiang, S.W. Effects of dietary oregano essential oil supplementation on the stress response, antioxidative capacity, and HSPs mRNA expression of transported pigs. Livest. Sci. 2015, 180, 143–149. [CrossRef] 11. Wilson, C.L.; Franklin, J.D.; Otto, B.E. Fruit volatiles inhibitory to monilinia-fructicola and botrytis-cinerea. Plant Disease. 1987, 71, 316–319. [CrossRef] 12. Zhang, W.X.; Zhang, Y.; Zhang, X.W.; Deng, Z.X.; Liu, J.X.; He, M.L.; Wang, H.F. Effects of dietary supplementation with combination of tributyrin and essential oil on gut health and microbiota of weaned piglets. Animals 2020, 10, 180. [CrossRef] [PubMed] 13. Nabuurs, M.J.A.; Hoogendoorn, A.; Vandermolen, E.J.; Vanosta, A.L.M. Villus height and crypt depth in weaned and unweaned pigs, reared under various circumstances in the Netherlands. Res. Vet. Sci. 1993, 55, 78–84. [CrossRef] 14. Zhang, C.; Yu, M.; Yang, Y.; Mu, C.; Su, Y.; Zhu, W. Differential effect of early antibiotic intervention on bacterial fermentation patterns and mucosal gene expression in the colon of pigs under diets with different protein levels. Appl. Microbiol. Biotechnol. 2017, 101, 2493–2505. [CrossRef] [PubMed] 15. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods. 2013, 10, 996–998. [CrossRef] 16. Cole, J.R.; Wang, Q.; Cardenas, E.; Fish, J.; Chai, B.; Farris, R.J.; Kulam-Syed-Mohideen, A.S.; McGarrell, D.M.; Marsh, T.; Garrity, G.M.; et al. The Ribosomal Database Project: Improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37, 141–145. [CrossRef] 17. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [CrossRef] 18. Kind, T.; Wohlgemuth, G.; Lee, D.Y.; Lu, Y.; Palazoglu, M.; Shahbaz, S.; Fiehn, O. FiehnLib: Mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal.Chem. 2009, 81, 10038–10048. [CrossRef] 19. Omonijo, F.A.; Ni, L.; Gong, J.; Wang, Q.; Lahaye, L.; Yang, C. Essential oils as alternatives to antibiotics in swine production. Anim. Nutr. 2018, 4, 126–136. [CrossRef] Microorganisms 2021, 9, 1342 16 of 17 20. Wang, C.C.; Shen, Z.J.; Cao, S.T.; Zhang, Q.H.; Peng, Y.; Hong, Q.H.; Feng, J.; Hu, C.H. Effects of tributyrin on growth performance, intestinal microflora and barrier function of weaned pigs. Anim. Feed. Sci. Technol. 2019, 258, 114311. [CrossRef] 21. Chouhan, S.; Sharma, K.; Guleria, S. Antimicrobial activity of some essential oils—present status and future perspectives. Medicines 2017, 4, 58. [CrossRef] [PubMed] 22. Adams, T.B.; Cohen, S.M.; Doull, J.; Feron, V.J.; Goodman, J.I.; Marnett, L.J.; Munro, I.C.; Portoghese, P.S.; Smith, R.L.; Waddell, W.J.; et al. The FEMA GRAS assessment of hydroxy- and alkoxy-substituted benzyl derivatives used as flavor ingredi- ents. Food Chem. Toxicol. 2005, 43, 1241–1271. [CrossRef] [PubMed] 23. Rattaz, C.; Goubet, N.; Bullinger, A. The calming effect of a familiar odor on full-term newborns. J. Dev. Behav. Pediatr. 2005, 26, 86–92. [CrossRef] [PubMed] 24. Stelter, K.; Frahm, J.; Paulsen, J.; Berk, A.; Kleinwachter, M.; Selmar, D.; Danicke, S. Effects of oregano on performance and immunmodulating factors in weaned piglets. Arch. Anim. Nutr. 2013, 67, 461–476. [CrossRef] [PubMed] 25. Giannenas, I.; Florou-Paneri, P.; Papazahariadou, M.; Christaki, E.; Botsoglou, N.A.; Spais, A.B. Effect of dietary supplementation with oregano essential oil on performance of broilers after experimental infection with Eimeria tenella. Arch. Tierernahr. 2003, 57, 99–106. [PubMed] 26. Fuhrman, M.P.; Charney, P.; Mueller, C.M. Hepatic proteins and nutrition assessment. J. Am. Diet. Assoc. 2004, 104, 1258–1264. [CrossRef] 27. Konstantinos, G.; Dimitrios, S.; Konstantinos, P.; Stylianos, K.; Dimitrios, K.; Panagiotis, G.; Dimitrios, K. Feed supplemented with polyphenolic by product from olive mill wastewater processing improves the redox status in blood and tissues of piglets. Food Chem. Toxicol. 2015, 86, 319–327. 28. Ji, F.J.; Wang, L.X.; Yang, H.S.; Hu, A.; Yin, Y.L. Review: The roles and functions of glutamine on intestinal health and performance of weaning pigs. Animal. 2019, 13, 2727–2735. [CrossRef] 29. Kotunia, A.; Wolinski, J.; Laubitz, D.; Jurkowska, M.; Rome, V.; Guilloteau, P.; Zabielski, R. Effect of sodium butyrate on the small intestine development in neonatal piglets fed [correction of feed] by artificial sow. J. Physiol. Pharmacol. 2004, 55 Suppl 2, 59–68. 30. Droy-Lefaix, M.T.; Drouet, Y.; Geraud, G.; Hosford, D.; Braquet, P. Superoxide dismutase (SOD) and the paf-Antagonist (BN 52021) reduce small intestinal damage induced by ischemia-reperfusion. Taylor & Francis 2009, 13, 725–735. 31. Omonijo, F.A.; Liu, S.; Hui, Q.; Zhang, H.; Lahaye, L.; Bodin, J.-C.; Gong, J.; Nyachoti, M.; Yang, C. Thymol improves barrier function and attenuates inflammatory responses in porcine intestinal epithelial cells during lipopolysaccharide (LPS)-induced inflammation. J. Agric. Food Chem. 2018, 67, 615–624. [CrossRef] 32. Kang, S.N.; Lee, E.; Lee, M.K.; Lim, S.J. Preparation and evaluation of tributyrin emulsion as a potent anti-cancer agent against melanoma. Drug Deliv. 2011, 18, 143–149. [CrossRef] 33. Byrne, C.S.; Chambers, E.S.; Morrison, D.J.; Frost, G. The role of short chain fatty acids in appetite regulation and energy homeostasis. Int. J. Obes. 2015, 39, 1331–1338. [CrossRef] [PubMed] 34. Howard, M.D.; Gordon, D.T.; Pace, L.W.; Garleb, K.A.; Kerley, M.S. Effects of dietary supplementation with fructooligosaccharides on colonic microbiota populations and epithelial cell proliferation in neonatal pigs. J. Pediatr. Gastroenterol. Nutr. 1995, 21, 297–303. [CrossRef] [PubMed] 35. Xu, Y.T.; Liu, L.; Long, S.F.; Pan, L.; Piao, X.S. Effect of organic acids and essential oils on performance, intestinal health and digestive enzyme activities of weaned pigs. Anim. Feed. Sci. Technol. 2018, 235, 110–119. [CrossRef] 36. Le Chatelier, E.; Nielsen, T.; Qin, J.; Prifti, E.; Hildebrand, F.; Falony, G.; Almeida, M.; Arumugam, M.; Batto, J.M.; Kennedy, S.; et al. Richness of human gut microbiome correlates with metabolic markers. Nature 2013, 500, 541–546. [CrossRef] [PubMed] 37. Hua, Y.; Yun, X.; Kelsy, R.; Junjun, W.; Guolong, Z.; Jiangchao, Z.; Yingping, X. Gut microbiota is a major contributor to adiposity in pigs. Front. Microbiol. 2018, 9, 3045. 38. Yang, C.; Zhang, L.; Cao, G.; Feng, J.; Yue, M.; Xu, Y.; Dai, B.; Han, Q.; Guo, X. Effects of dietary supplementation with essential oils and organic acids on the growth performance, immune system, fecal volatile fatty acids, and microflora community in weaned piglets. J. Anim. Sci. 2018, 97, 133–143. [CrossRef] [PubMed] 39. Yang, C.B.; Chowdhury, M.A.K.; Hou, Y.Q.; Gong, J. Phytogenic compounds as alternatives to in-feed antibiotics: Potentials and challenges in application. Pathogens 2015, 4, 137–156. [CrossRef] 40. Qin, S.S.; Wu, C.M.; Wang, Y.; Jeon, B.; Shen, Z.Q.; Wang, Y.; Zhang, Q.; Shen, J.Z. Antimicrobial resistance in Campylobacter coli isolated from pigs in two provinces of China. Int. J. Food Microbiol. 2011, 146, 94–98. [CrossRef] [PubMed] 41. Schautteet, K.; Vanrompay, D. Chlamydiaceae infections in pig. Vet. Res. 2011, 42, 29. [CrossRef] [PubMed] 42. Marchandin, H.; Jumas-Bilak, E. The Family Veillonellaceae. In The Prokaryotes; Rosenberg, E., Ed.; Springer Berlin Heidelberg: Berlin/Heidelberg Germany, 2014; pp. 433–453. 43. McCormack, U.M.; Curiao, T.; Buzoianu, S.G.; Prieto, M.L.; Ryan, T.; Varley, P.; Crispie, F.; Magowan, E.; Metzler-Zebeli, B.U.; Berry, D.; et al. Exploring a possible link between the intestinal microbiota and feed efficiency in pigs. Appl. Environ. Microbiol. 2017, 83, e00380-00317. [CrossRef] 44. Aar, P.; Molist, F.; Klis, J. The central role of intestinal health on the effect of feed additives on feed intake in swine and poultry. Anim. Feed Sci. Technol. 2017, 233, 64–75. 45. Fiorucci, S.; Distrutti, E. Bile acid-activated receptors, intestinal microbiota, and the treatment of metabolic disorders. Trends Mol. Med. 2015, 21, 702–714. [CrossRef] Microorganisms 2021, 9, 1342 17 of 17 46. Pina, M.M.; Giuseppina, M.R.; Gina, C.; Giovanna, T.; Chiara, D.F.; Serena, A.; Marina, P.; Claudio, P.; Francesca, D.G.; Adriano, L.; et al. Butyrate regulates liver mitochondrial function, efficiency, and dynamics in insulin-resistant obese mice. Diabetes 2017, 66, 1405–1418. 47. Zhou, L.; Zheng, H.; Tang, Y.; Yu, W.; Gong, Q. Eugenol inhibits quorum sensing at sub-inhibitory concentrations. Biotechnol. Lett. 2012, 35, 631–637. [CrossRef] [PubMed] 48. Geng, M.M.; Li, T.J.; Kong, X.F.; Song, X.Y.; Chu, W.Y.; Huang, R.L.; Yin, Y.L.; Wu, G.Y. Reduced expression of intestinal N-acetylglutamate synthase in suckling piglets: A novel molecular mechanism for arginine as a nutritionally essential amino acid for neonates. Amino Acids. 2011, 40, 1513–1522. [CrossRef] [PubMed] 49. Dasgupta, T.; Hebbel, R.P.; Kaul, D.K. Protective effect of arginine on oxidative stress in transgenic sickle mouse models. Free Radic. Biol. Med. 2006, 41, 1771–1780. [CrossRef]