# Dietary supplementation with fermented Mao-tai lees beneficially affects gut microbiota structure and function in pigs

Huan Li, Huawei Li, Peifeng Xie, Zhihua Li, Yulong Yin, Francois F. Blachier, Xiangfeng X. Kong

► To cite this version:

Huan Li, Huawei Li, Peifeng Xie, Zhihua Li, Yulong Yin, et al.. Dietary supplementation with fermented Mao-tai lees beneficially affects gut microbiota structure and function in pigs. AMB Express, 2019, 9, 10.1186/s13568-019-0747-z. http://hal.inrae.fr/hal-02618953

HAL Id: hal-02618953

https://hal.inrae.fr/hal-02618953

Submitted on 25 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

ORIGINAL ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

Check for
updates

# Dietary supplementation with fermented Mao-tai lees beneficially affects gut microbiota structure and function in pigs

Huan Li<fnref n="1" />,²†, Huawei Li<fnref n="1" />†, Peifeng Xie<fnref n="1" />, Zhihua Li<fnref n="1" />, Yulong Yin<fnref n="1" />, Francois Blachier³ and Xiangfeng Kong<fnref n="1" />*

## Abstract

Gut microbiota positively contribute to livestock nutrition and metabolism. The manipulation of these microbes may improve animal health. Some feed additives improve livestock health and metabolism by regulating gut microbiota composition and activity. We fed hybrid pigs diets supplemented with 0% (control), 5% (treat 1), 10% (treat 2), or 15% (treat 3) fermented Mao-tai lees (FML) for 90 days. Short-chain fatty acids (SCFAs), bioamines, and microbial communities found in colonic contents were analyzed to investigate microbiota composition and metabolic profiles. Concentrations of straight-chain fatty acids (e.g., acetate, propionate, and butyrate) and tyramine increased with FML supplementation content. Contrary to the minor effects of 5% and 10% FML on gut microbiota, 15% FML influenced beta diversity (Jaccard or Bray-Curtis dissimilarity) but not alpha diversity (number of operational taxonomic units and Shannon diversity) of pig gut microbial communities compared to the control group. Notably, 15% FML animals were characterized by a higher abundance of potentially beneficial bacteria (Lactobacillus and Akkermansia) but lower abundances of potential pathogens (Escherichia). Numerous genes associated with metabolism (e.g., starch, sucrose, and sulfur-compounds metabolism) showed a higher relative abundance in the 15% FML than in the control group. Additionally, most Phascolarctobacterium, Treponema, Prevotella, and Faecalibacterium bacterial markers in the 15% FML group were positively correlated with straight-chain fatty acid concentrations, suggesting that these bacteria are likely associated with SCFA production. Taken together, our findings demonstrate the beneficial effects of 15% FML on fermentation of undigested compounds and gut microbiota composition in the colon. Thus, 15% FML supplementation in pig feed may possibly represent a way to optimize pig colon health for livestock farming.

**Keywords:** Gut microbiota, Fermented Mao-tai lees, Short-chain fatty acids, Beta diversity, Beneficial bacteria, Livestock farming

## Introduction

Healthy livestock farming technologies and management practices are a central concern in animal farming industry, which must provide high-quality meat products to safely meet human consumption requirements. In particular, pigs are one of the most important economic

livestock species in numerous countries. Pigs represent the largest livestock product in China, with great commercial and economic values worldwide (Bai et al. 2014; Sun and Jia 2015; Yu and Abler 2014). However, intensive or large-scale livestock farming causes various problems such as serious feed resource shortages and grain competition between humans and livestock (Ilea 2009). In addition, the utilization of antibiotics in livestock husbandry potentially poses a threat to livestock health, and it could lead to environmental pollution as well as a decline in animal immunity, gut microbiota disorders, and an increased risk of the spread of antibiotic-resistant bacteria genes in the environment (Baquero et al. 2008; Wu et al. 2010; Zhao et al. 2018). Residues of antibiotics

<figure><img src="image_4.png" /></figure>

<footnote>*Correspondence: nnkxf@isa.ac.cn</footnote>

<footnote>†Huan Li and Huawei Li contributed equally to this work</footnote>

<footnote n="1">¹ Hunan Provincial Key Laboratory of Animal Nutritional Physiology and Metabolic Process, Key Laboratory of Agro-ecological Processes in Subtropical Region, National Engineering Laboratory for Pollution Control and Waste Utilization in Livestock and Poultry Production, Institute of Subtropical Agriculture, Chinese Academy of Sciences, Changsha, Hunan 410125, China</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
in livestock products are likely harmful to human health (Marshall and Levy 2011). However, the restricted use of antibiotics in feed used in livestock farming may increase disease rates in animals and decrease amounts of global products. Therefore, it is imperative to develop safe and effective feed additives that could replace anti-biotics, thus improving livestock health and products. Recently, studies focused on replacing antibiotics in livestock husbandry have examined the effects of probiotics, prebiotics, synbiotics, and other dietary additives on livestock (Abudabos et al. 2017; Markowiak and Slizewska 2018; Marshall and Levy 2011). Moreover, these studies have aroused the interests of scientists globally, since these newly developed feed supplements were found to improve livestock health and growth (Markowiak and Slizewska 2018).

Gut microbiota has several beneficial functions in hosts, including food digestion, energy harvesting, immune regulation, and resistance against pathogens (Rooks and Garrett 2016; Stanley et al. 2016; Tremaroli and Bäckhed 2012). Evidence suggests that probiotics and prebiotics could beneficially modify gut microbiota composition and activity in humans and animals (Barba-Vidal et al. 2018; Moura et al. 2007), thereby improving host metabolism and health, notably at the intestinal level. For example, the probiotics *Lactobacillus (Lactobacillus reuteri avibro)* and *Bacillus (Bacillus subtilis and Bacillus licheniformis)* exhibited the ability to increase nutrient digestibility and animal performance, and they reduced the abundance of pathogens (*Salmonella* and *Escherichia coli*) (Ahmed et al. 2014). In pigs, some prebiotic products, such as xylo-oligosaccharides (XOS), β-mannanase, mannan-oligosaccharides, and yeast cultures have been used to modify gut microbiota composition, improve host immune response, stimulate the growth of more beneficial bacteria, and inhibit the colonization or abundance of pathogens by producing antimicrobial substances in humans and livestock (Barros et al. 2015; De Maesschalck et al. 2015; Liu et al. 2018; Rastall and Gibson 2015). For instance, in broiler chicken, XOS feed additives modulate the gut microbiota and increase the abundance of the potential beneficial bacteria *Lactobacillus* and *Bifidobacterium* (Pourabedin et al. 2015). Moreover, feed containing yeast cultures increased the abundance of *Clostridium* and *Methylobacterium*, which are important cellulose-degrading bacteria that may help herbivorous grass carp (*Ctenopharyngodon idellus*) to degrade ingested plant-based food (Liu et al. 2018).

Recently, liquor lees were been used as promising feed additive candidates, and they were found to be able to shape animal gut microbiota. For example, the content of dried distiller grains with solubles (DDGS) in feed is positively correlated with gut microbial diversity in birds

(Abudabos et al. 2017), and increased microbial diversity is associated with metabolic activity and health in human and animal hosts (Li et al. 2017b, 2018; Tap et al. 2015). However, studies of the effects of liquor lees on livestock gut microbiota are still in their infancy and need further development.

Guizhou Maotai liquor is a fragrant and tasty wine (Wu et al. 2013), and it is renowned as the best Chinese sauce fragrance liquor. The production of Maotai wine was about 38,700 t, and the solid by-products (including lees) were about 110,000 t (Li et al. 2015). DDGS is a suitable substrate for solid-state fermentation (SSF), and only moderate changes were found in its nutritional profile after SSF (Yang et al. 2012). Fermented Maotai lees (FML) also have a high nutritional value (e.g., high protein, cellulose, amino acids, and organic acids), so they may be of interest for livestock farming. Here, we used MiSeq sequencing of 16S rRNA genes to evaluate the effects of FML on the gut microbiota diversity and also measured the gut metabolic products (short-chain fatty acids (SCFAs) and bioamines) in pigs. We solved three questions: (1) Do FML influence the SCFA and bioamine pofiles in the pig gut? (2) Whether FML modify gut microbiota composition and function? (3) Do FML increase the abundance of potential beneficial bacteria and decrease the abundance of specific pathogens?

## Materials and methods

### Experimental design and ethical standards

A total of 24 Duroc × Large White × Landrace hybrid barrows (male), with an initial body weight (BW) of 42.28±1.23 kg (mean±SE), were fed a corn and soy-bean meal-based diet (Additional file 1: Table S1), which met the National Research Council (NRC 2012) requirements for growing-finishing pigs. After these pigs were fed the same diet for 1 week, the animals were randomly arranged to one of the four treatments. Each treatment group had 6 replicates. Pigs in the control group were fed a basic diet, whereas those in the experimental groups were fed the basic diet containing 5%, 10%, or 15% fermented Maotai lees (FML) [0% as Control, 5% (treat 1), 10% (treat 2) and 15% (treat 3)] during the experimental period. Notably, the FML used in the present study were provided by the Road Biological Environmental Co., Ltd., Sichuan, China. The FML products were fermented by saccharomyces yeasts with acid-resisting, overproducing enzyme, and high vigor. After fermentation, the mixture was dried at 60 °C using laboratory incubator, and then smashed. The FML were then stored in -20 °C for experimental use. The determined nutrient levels (%) of the FML based on dry matter content (92.97%) were as follows: ash, 9.28; gross energy (GE), 18.29; crude protein (CP), 23.96; ether extract (EE), 5.39; crude fiber (CF),
17.67; acid detergent fiber (ADF), 38.06; neutral detergent fiber, 47.28; Ca, 0.53; P, 0.55.

Experimental pigs were housed in cages (3.5 m × 5.0 m) equipped with feed intake recording equipment (Beijing Hamoer Automation Equipment Co., Ltd., Beijing, China). The space provided by the equipment allows one pig at a time to have ad libitum access to the diet. Animals had 24 h access to feed and water during the whole course of the experiment. Pigs were labeled with an individual electronic ear marker. The feeding trial lasted for 90 days. The final body weight of each pig was also measured to evaluate the growth status.

The experimental design and procedures in this study were reviewed and approved by the Animal Care and Use Committee of the Institute of Subtropical Agriculture, Chinese Academy of Science. Processing of animal experiments and sample collection strictly followed the relevant guidelines.

## Sample collection

At the end of the feeding trial, the diet was removed 12 h before slaughter. All pigs were transported from the farm to the processing facility (approximately 40 km) at 7:00 h, and sacrificed at 19:00 h under commercial conditions using electrical stunning (120 V, 200 Hz). Colonic contents were collected into 50 ml sterile tubes. After complete mixing, all the gut contents were stored at -20 °C for the following microbiota and metabolic analysis.

## Measurements of short-chain fatty acids (SCFAs) and bioamines

Gut SCFAs (mg/g), including straight-chain fatty acids (acetate, propionate, butyrate, and valerate) and BCFAs (branched-chain fatty acid, including isobutyrate and isovalerate) were analyzed using gas chromatography as previously described (Ji et al. 2018). Bioamines (µg/g), including putrescine, cadaverine, tyramine, spermidine, and spermine were measured using high-performance liquid chromatography as previously described (Ji et al. 2018).

## DNA extraction and MiSeq sequencing of microbial 16S rRNA gene

Genomic DNA of colonic contents was extracted using a QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The DNA concentration of each sample was measured with a NanoDrop® ND-1000 instrument (NanoDrop Technologies Inc., Dover, USA). The extracted DNA was diluted to 10 ng/µl for the following polymerase chain reaction (PCR) amplification. The protocols of PCR amplification, gel extraction, and sequencing library construction were described previously (Li et al. 2016a, b, c). The universal

primers 338F (5'-ACTCCTACGGGAGGCAGCAG-3') and 806R (5'-GGACTACHVGGGTWTCTAAT-3') with 12 nt unique barcode at 5'-end of 338F were used to amplify the V3-V4 region of bacterial 16S rRNA gene (Jiang et al. 2018). After PCR amplification, amplicons were extracted from 1.2% agarose gels and purified using SanPrep DNA Gel Extraction Kit (Sangon Bio-tech, Shanghai, China) and quantified with a NanoDrop® ND-1000 instrument (NanoDrop Technologies Inc., Dover, USA). Purified amplicons were mixed in equal molar together and underwent paired-end sequencing using an Illumina MiSeq sequencer (MiSeq Reagent Kit V2, California, USA, 500 cycles).

## Bioinformatics analysis

Bioinformatics analysis followed the analytical methods of Li et al. (2016b). Briefly, the raw sequences were analyzed using QUIME Pipeline-Version 1.9.0 (http://qiime.org/tutorials/tutorial.html) (Caporaso et al. 2010). All sequences were trimmed and assigned to each sample based on their unique barcodes only if they completely matched their unique barcodes. The FLASH-1.2.8 software (Caporaso et al. 2012) was used to merge the overlapping paired-end reads. The merged sequences with high quality (reads length more than 300 bp, without ambiguous base "N", and average base quality score more than 30) were kept for further analysis. Then the aligned 16S rRNA gene sequences underwent a chimera check using the Uchime algorithm (Edgar et al. 2011). Due to possible contamination of chloroplast sequences in PCR amplification, the Metaxa2 software tool was used to remove chloroplast sequences from our large sequencing datasets (Bengtsson-Palme et al. 2015). After filtering chimeras and chloroplasts, the remaining sequences were clustered into operational taxonomic units (OTUs) at a 97% identity threshold with an open-reference OTU picking method using the Uclust algorithm (Edgar 2010). Those sequences not classifying to bacteria (Eukaryota and Archaea lineages) were filtered out. Singleton OTUs were also removed. The most abundant sequences within each OTU were defined as "representative sequences". Taxonomic classification of representative sequences was implemented using the Ribosomal Database Project classifier in the QUIME platform (Wang et al. 2007).

To minimize the impact of uneven sequencing depth, each sample was rarefied to 32,805 sequences. To evaluate alpha diversity indices, Goods coverage, Chao1, observed OTUs, Shannon diversity and evenness were calculated. The rarefaction curves were generated from the observed OTUs at the OTU level. To assess beta diversity, Jaccard and Bray-Curtis distance metrics were produced through the QUIME pipeline. Jaccard distance is based on the presence/absence of OTUs/species (Jaccard
1912), whereas Bray–Curtis distance is based on the both OTU abundance and presence/absence (Bray and Curtis 1957). Differences in overall bacterial community structure among groups were visualized using the non-metric multidimensional scaling (NMDS) plots of the two similarity matrices.

## Statistical analysis

Permutational multivariate analysis of variance (PERMANOVA) was used to reveal whether the structures of gut microbiota were significantly different among groups based on the Jaccard and Bray–Curtis distance matrices using “adonis” in the R ‘vegan’ package (Li et al. 2017a). The model also comprised other variables, including animal body weight. Taxonomic profiles were evaluated at the phylum and genus levels. Differences in relative abundances of the genera were performed through the *group significance.py* script in the QIIME platform with one-way analysis of variance (one-way ANOVA) with Tukey’s post hoc test. *P*-values were corrected using false discovery rate (FDR). Comparison of the alpha diversity indices among groups was tested using Mann–Whitney *U* tests. Statistical significance was accepted at a *P*<0.05. Linear discriminant analysis effect size (LEfSe) (Nicola et al. 2011) was used to identify bacterial OTU biomarkers at OTU level among groups based on a *P*<0.05 and LDA score > 2.0. This analysis was performed online in the Galaxy workflow framework (http://huttenhower.sph.harvard.edu/galaxy/).

The Spearman correlation analysis among SCFA levels, bioamines, bacterial biomarkers (based on LefSe results) was achieved using function “cor.test” of the package “Stats” in R. One heatmap plot of the correlation values were produced using package “gplots”. One-way ANOVA with Tukey’s post hoc test was used to evaluate the difference of SCFA levels and bioamines among groups. Absolute SCFA or bioamine concentrations (µmol/g) were converted to molar percentages within each individuals (Li et al. 2018), and then differences in overall SCFA or bioamine profile were visualized using the NMDS plots of Bray–Curtis dissimilarity indices. PERMANOVA was also used to reveal whether the structures of SCFA and bioamine profiles were significantly different among groups based on the Bray–Curtis distance matrices using “adonis” in the R ‘vegan’ package.

## Predicted metagenomes

PICRUSTv1.0.0 (Langille et al. 2013) was used to predict abundances of functional gene from OTU abundances based on the 16S rRNA gene sequences. We only focused on the gene functions associated with metabolism between the control and treat 3 groups. Two-tailed *t* tests

(Bonferroni-corrected) were performed to test the differences of gene functions between these two groups.

## Nucleotide sequence accession numbers

The original 16S rRNA data were available at the European Nucleotide Archive by accession number PRJEB27667 (http://www.ebi.ac.uk/ena/data/view/PRJEB27667).

## Results

### The effects of FML on the weight gain, gut SCFAs, and bioamines of pigs

Dietary FML supplementation had no significant effects on the weight gain of pigs (Table 1). However, the concentrations of acetate, propionate, butyrate, valerate, total straight-chain fatty acids, and total SCFAs in the gut contents significantly increased with FML content in the feed (Spearman correlation analysis, all *P* < 0.05; Table 1), whereas there were no significant differences in BCFA (including isobutyrate, isovalerate, and total BCFA) concentration among groups. Over SCFA profile structure showed no significant separation between control and treat groups (Fig. 1a; Additional file 1: Table S2; PERMANOVA, all *P* > 0.05). However, final body weight had significant influences in shaping the SCFA profile (PERMANOVA, $R^2$ = 0.149, *P* = 0.027).

In addition, the concentrations of most bioamines showed no differences between groups, whereas the spermidine concentration of 10% FML group was significantly higher than that of control group (Table 1). Over bioamine profile structure also showed no significant separation between control and treat groups (Fig. 1b; Additional file 1: Table S2; PERMANOVA, all *P* > 0.05). Final body weight also had no significant impacts in shaping the bioamine profile (PERMANOVA, *P* > 0.05).

### The influences of FML on alpha and beta diversity of pig gut microbiota

A total of 2,394,088 sequences were generated from the 24 pig samples. After filtering out low-quality sequences, chimeras, chloroplasts, singletons, and those sequences not classifying to bacteria, we obtained 922,668 valid sequences. To compare samples with different sequencing depth, each sample was rarefied to 32,805 sequences. At a threshold of 97% sequence similarity, 21,920 unique OTUs were identified using Uclust clustering. The rarefaction curve of observed OTUs across all samples reach stable values (Additional file 1: Fig. S1), indicating that our sequencing had captured most species although additional sequencing may detect some rare OTUs. In addition, the Goods coverage (mean ± SE) of gut microbiota at OTU level across all samples was 95.01% ± 0.26%, confirming that our sequencing depth is enough to detect a
<table><caption>Table 1 The comparison of body weight, fecal short chain fatty acids (SCFAs), and bioamine concentrations in pigs with fermented Mao-tai lees (FML) contents</caption><thead><tr><th>Items</th><th>Control</th><th>Treat 1</th><th>Treat 2</th><th>Treat 3</th></tr></thead><tbody><tr><td><strong>Body weight (kg)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Initial Body weight</td><td>46.250 ± 3.080a</td><td>38.620 ± 1.990a</td><td>42.750 ± 2.620a</td><td>41.480 ± 1.310a</td></tr><tr><td>Final body weight</td><td>111.280 ± 5.660a</td><td>108.120 ± 4.330a</td><td>112.950 ± 4.100a</td><td>106.750 ± 3.980a</td></tr><tr><td>Weight gain</td><td>65.030 ± 3.450a</td><td>69.500 ± 3.360a</td><td>70.200 ± 3.240a</td><td>65.270 ± 3.620a</td></tr><tr><td><strong>SCFAs (mg/g)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Acetate</td><td>4.097 ± 0.565ab</td><td>4.256 ± 0.271b</td><td>5.813 ± 0.071a</td><td>6.199 ± 0.871ab</td></tr><tr><td>Propionate</td><td>2.012 ± 0.289a</td><td>2.066 ± 0.260a</td><td>2.833 ± 0.107a</td><td>2.670 ± 0.319a</td></tr><tr><td>Isobutyrate</td><td>0.115 ± 0.022a</td><td>0.163 ± 0.011a</td><td>0.151 ± 0.023a</td><td>0.171 ± 0.029a</td></tr><tr><td>Butyrate</td><td>1.514 ± 0.289a</td><td>1.436 ± 0.221a</td><td>2.256 ± 0.154a</td><td>1.890 ± 0.184a</td></tr><tr><td>Isovalerate</td><td>0.202 ± 0.038a</td><td>0.302 ± 0.023a</td><td>0.274 ± 0.041a</td><td>0.313 ± 0.052a</td></tr><tr><td>Valerate</td><td>0.172 ± 0.022a</td><td>0.237 ± 0.026ab</td><td>0.307 ± 0.018b</td><td>0.262 ± 0.021b</td></tr><tr><td>Total straight-chain fatty acids</td><td>7.622 ± 1.120a</td><td>7.758 ± 0.727a</td><td>10.903 ± 0.315a</td><td>10.760 ± 1.307a</td></tr><tr><td>Total BCFA</td><td>0.317 ± 0.060a</td><td>0.465 ± 0.034a</td><td>0.425 ± 0.064a</td><td>0.484 ± 0.081a</td></tr><tr><td>Total SCFA</td><td>7.940 ± 1.176a</td><td>8.223 ± 0.709a</td><td>11.327 ± 0.366a</td><td>11.244 ± 1.340a</td></tr><tr><td><strong>Bioamines (µg/g)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Putrescine</td><td>5.259 ± 1.065a</td><td>4.322 ± 0.499a</td><td>4.855 ± 0.470a</td><td>3.362 ± 0.691a</td></tr><tr><td>Cadaverine</td><td>11.140 ± 5.438a</td><td>4.188 ± 1.811a</td><td>4.079 ± 0.991a</td><td>2.667 ± 1.342a</td></tr><tr><td>Tyramine</td><td>1.013 ± 0.259a</td><td>1.400 ± 0.257a</td><td>1.294 ± 0.341a</td><td>2.082 ± 0.347a</td></tr><tr><td>Spermidine</td><td>5.584 ± 0.826a</td><td>8.746 ± 0.817ab</td><td>9.240 ± 0.709b</td><td>7.014 ± 1.244ab</td></tr><tr><td>Spermine</td><td>0.547 ± 0.090a</td><td>0.731 ± 0.095a</td><td>0.869 ± 0.161a</td><td>0.454 ± 0.091a</td></tr></tbody></table>

Data was expressed as mean ± SE obtained from 6 individual pigs each. One-way ANOVA with Tukey's post hoc test or Mann-Whitney U-tests was used to evaluate the difference of fecal metabolites in different groups. Significant difference is indicated by different letters. The group control, treat 1, treat 2, and treat 3 signify that feed is supplemented with 0, 5%, 10%, or 15% FML, respectively

majority of gut bacterial species. The alpha diversity values (including Chao1, observed OTUs, Shannon diversity, and evenness) had no significant differences between the control and treat groups (Mann-Whitney U tests, all *P* > 0.05; Table 2).

However, the beta diversity values of gut microbiota had significant differences across groups based on Jaccard (Fig. 1c; Table 3; PERMANOVA, *R*²=0.057, *P*=0.003) and Bray-Curtis distance matrices (Fig. 1d; Table 3; PERMANOVA, *R*²=0.212, *P*=0.008). When we compared the pairwise dissimilarity matrices between any two groups, our results showed that only the community structure of 15% FML group had significant differences with that of the control group (Table 3) based on Jaccard (PERMANOVA, *R*²=0.11, *P*=0.015) or Bray-Curtis distance (PERMANOVA, *R*²=0.171, *P*=0.017), indicating that only 15% FML influences the community structure of pig gut microbiota. In addition, the group 15% FML showed different community structures with 5% FML based on Jaccard (PERMANOVA, *R*²=0.114, *P*=0.023) or Bray-Curtis distance (PERMANOVA, *R*²=0.228, *P*=0.008). However, final body weight had no significant impacts in shaping the beta diversity values of pig gut microbiota (Jaccard, PERMANOVA, *P* > 0.05; Bray-Curtis, *P* > 0.05).

## FML treatments affect the taxonomic composition of gut microbiota

Across all samples, approximately 99% of the total sequences were assigned into 23 phyla. Among these phyla, Firmicutes (62.1%), Bacteroidetes (27.0%), Spirochaetes (3.2%), Proteobacteria (3.0%), Tenericutes (2.2%), and Cyanobacteria (1.2%) were the six dominant bacterial taxa (mean relative abundance > 1%) in the pig gut. Some rare phyla included Actinobacteria, TM7, and WPS-2. Those phyla with mean relative abundance > 0.1% are shown in Fig. 2a. At genus level, S24-7, *Prevotella*, Lachnospiraceae, Ruminococcaceae, *Lactobacillus*, *SMB53*, Clostridiales, *Oscillospira*, Bacteroidales, Clostridiaceae, *Treponema*, *Escherichia*, *Clostridium*, *Ruminococcus*, RF39, *Roseburia*, *p-75-a5*, *Lachnospira*, *Streptococcus*, YS2, and *Coprococcus* were the dominant bacterial genera (Fig. 2b).

We further compared the significant difference of bacterial genera among groups. A total of 14 genera showed significant difference among groups (Table 4). Among these genera, the relative abundance of *Roseburia*, *Desulfovibrio*, *Phascolarctobacterium*, *Oxalobacter*, *Defluviitalea*, and *Prevotella* in the group 15% FML was significantly higher than that of the control group (one-way ANOVA, all *P* < 0.05). A total of 8 genera, such as
<figure><img src="image_5.png" /><figcaption>**Fig. 1** NMDS plots of dissimilarity metrics comparing the profiles of gut microbiota, SCFAs and bioamines among groups. **a** SCFA profile based on Bray-Curtis dissimilarities. **b** Bioamine profile based on Bray-Curtis dissimilarities. **c** Gut microbiota structure based on Jaccard dissimilarities. **d** Gut microbiota structure based on Bray-Curtis dissimilarities</figcaption></figure>

<table><caption>**Table 2** Alpha diversity values of the gut microbiota were evaluated using OTUs defined at 97% similarity threshold</caption><thead><tr><th>Items</th><th>Control</th><th>Treat 1</th><th>Treat 2</th><th>Treat 3</th></tr></thead><tbody><tr><td>Goods coverage</td><td>0.948 ± 0.004a</td><td>0.953 ± 0.005a</td><td>0.954 ± 0.005a</td><td>0.945 ± 0.007a</td></tr><tr><td>Chao1</td><td>6049.249 ± 493.238a</td><td>5584.114 ± 600.391a</td><td>5510.697 ± 578.835a</td><td>6693.918 ± 904.801a</td></tr><tr><td>Observed OTUs</td><td>2974.500 ± 223.060a</td><td>2667.667 ± 285.628a</td><td>2648.667 ± 249.738a</td><td>3184.167 ± 289.860a</td></tr><tr><td>Shannon diversity</td><td>7.574 ± 0.376a</td><td>7.132 ± 0.408a</td><td>7.755 ± 0.216a</td><td>8.216 ± 0.214a</td></tr><tr><td>Evenness</td><td>0.657 ± 0.027a</td><td>0.627 ± 0.028a</td><td>0.684 ± 0.012a</td><td>0.707 ± 0.011a</td></tr></tbody></table>

Data was expressed as mean ± SE obtained from 6 individual pigs each. Mann-Whitney U-tests were used to evaluate the difference in different groups. The group control, treat 1, treat 2, and treat 3 signify that feed is supplemented with 0, 5%, 10%, or 15% fermented Mao-tai lees, respectively. Significant difference is indicated by different letters

*Butyricicoccus*, *Dorea*, *Clostridium*, *SMB53*, and *Peptococcus*, showed a lower abundance than that of the control group (all *P* < 0.05).

To further compare the taxonomic difference between the control and treat groups, we used LEfSe analysis to differential abundance of bacterial taxa at OTU level among different treatments (Fig. 3). The results showed

that a total of 25 bacterial biomarkers/OTUs were differentially abundant among the four groups. OTU54 (belonging to *Lactobacillus reuteri*), OTU57 (*Turicibacter*), OTU73 (Clostridiaceae), OTU72 (Clostridiaceae), OTU74 (*Clostridium butyricum*), OTU75 (*Clostridium celatum*), and OTU76 (*SMB53*) were the dominant microbes in the 5% FML group, while OTU
<table><caption>**Table 3** PERMANOVA showing different community compositions and structures among different groups</caption><thead><tr><th rowspan="2">Items</th><th colspan="2">Jaccard</th><th colspan="2">Bray-Curtis</th></tr><tr><th>R<sup>2</sup></th><th>P</th><th>R<sup>2</sup></th><th>P</th></tr></thead><tbody><tr><td>All</td><td>0.057</td><td><strong>0.003*</strong></td><td>0.212</td><td><strong>0.008*</strong></td></tr><tr><td>Control vs treat 1</td><td>0.090</td><td>0.475</td><td>0.079</td><td>0.554</td></tr><tr><td>Control vs treat 2</td><td>0.103</td><td>0.052</td><td>0.135</td><td>0.113</td></tr><tr><td>Control vs treat 3</td><td>0.110</td><td><strong>0.015*</strong></td><td>0.171</td><td><strong>0.017*</strong></td></tr><tr><td>Treat 1 vs treat 2</td><td>0.099</td><td>0.15</td><td>0.154</td><td>0.053</td></tr><tr><td>Treat 1 vs treat 3</td><td>0.114</td><td><strong>0.023*</strong></td><td>0.228</td><td><strong>0.008*</strong></td></tr><tr><td>Treat 2 vs treat 3</td><td>0.099</td><td>0.104</td><td>0.132</td><td>0.066</td></tr></tbody></table>

Significant differences (*P* < 0.05) are indicated by asterisk (*). The group control, treat 1, treat 2, and treat 3 signify that feed is supplemented with 0, 5%, 10%, or 15% fermented Mao-tai lees, respectively

29 (*Prevotella*) and OTU115 (*Oscillospira*) were mostly in the control group. A total of 13 OTUs, including OTU117 (*Ruminococcus bromii*), OTU121 (*Phascolarcto-bacterium*), OTU132 (*Treponema*), OTU 25, 26, and 31 (*Prevotella*), OTU106 (*Faecalibacterium*), OTU94 (*Rose-buria*), and OTU36 (*Prevotella copri*), were significantly enriched in the 15% FML group.

### Comparison of the potential beneficial bacteria and pathogens abundance

In order to detect the effects of FML on the potential probiotics (*Bifidobacterium*, *Bacillus*, *Lactobacil-lus*, *Akkermansia*, and *Faecalibacterium*) and pathogens (*Escherichia*) in pigs (Barba-Vidal et al. 2018; Rastall and Gibson 2015), we compared the abundance of these gen-era among treatments (Fig. 4). We found that the relative abundance of *Bacillus*, *Lactobacillus*, *Akkermansia*, and *Faecalibacterium* in 15% FML group is higher than that of the control group. However, the 15% FML showed a lower abundance of *Escherichia* than that of the control group.

### The differences of predicted gene functions across groups

Based on the functional prediction of 16S rRNA gene sequences, the difference of gene functions among groups were visualized based on NMDS plot of Bray-Curtis dissimilarity at level 3. Overall functional profiles showed no significant differences among groups (PERMANOVA, R<sup>2</sup> = 0.197, *P* > 0.05, Additional file 1: Fig. S2). When we compared the differences between the 15% FML and control groups, a total of 29 gene functions at level 3 showed significant differences between the two groups. Compared with the control group, the 15% FML was more abundant in 15 gene functions associated with metabolism, such as flavone and flavonol biosynthesis, phenylpropanoid

biosynthesis, starch and sucrose metabolism, carbon fixation in photosynthetic organisms, sulfur metabolism, nicotinate and nicotinamide metabolism, retinol metabolism, and cyanoamino acid metabolism (Additional file 1: Table S3).

### Correlations between bacterial biomarkers and gut metabolites

To understand the relationship between bacterial biomarkers and gut metabolites, we calculated the spearman correlations between those bacterial OTUs (based on LefSe analysis) and SCFAs or bioamines (Fig. 5). Our results showed that acetate was negatively correlated with OTU57, OTU72, OTU75, or OTU76, whereas was positively associated with OTU17, OTU25, OTU26, OTU36, OTU47, OTU82, OTU94, OTU106, OTU118, and OTU121. Propionate only showed positive correlations with OTU17, OTU47, OTU118, and OTU121. Butyrate was related to OTU17, OTU36, OTU47, OTU106, OTU118, and OTU121. These results indicated that most OTUs enriched in the 15% FML group showed correlations with SCFAs. In addition, OTU22 was positively correlated with Putrescine. Tyramine showed positive associations with OTU40, OTU49, and OTU106.

### Discussion

In the field of animal nutrition, the development of new potential prebiotics in animal husbandry industries is of current interest. Recently, high-throughput sequencing of 16S rRNA gene amplicons has been used to investigate potential beneficial effects of various feed additives (e.g., xylo-oligosaccharides and yeast cultures) on the gut microbial communities of farmed animals (Liu et al. 2018; Pourabedin et al. 2015). However, relatively few studies focused on the effects of liquor lees on the composition and function of livestock gut microbiota. This study is the first to report the influence of FML on pig gut microbiota and metabolic profiles. Our results indicated that dietary supplementation with 15% FML significantly influenced gut microbial community structure and bacterial metabolite concentrations in pigs, and those abundant OTUs in the 15% FML group were positively associated with the fermentation of dietary fiber. Moreover, 15% FML improved the relative abundance of potential beneficial bacteria and decreased the abundance of pathogens. These results greatly enhanced our understanding of the potential beneficial effects of FML on pigs.

### FML improve the fermentation ability of gut microbiota for dietary fiber in pigs

Animals acquire nutrients and energy via bacterial fermentation of ingested feed. SCFAs that are major end
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Taxonomic compositions of gut bacterial communities with different dietary treatments with fermented Mao-tai lees in pigs. **a** The relative abundance of bacterial phyla. **b** Their relative abundance of top 21 bacterial taxa within a group at genus level</figcaption></figure>

metabolic products result from the bacterial fermentation of dietary fiber (including cellulose, lignin, and other polysaccharides) in mammalian guts (den Besten et al. 2013). Our results showed that higher FML content in feed were associated with increased concentrations of straight-chain fatty acids, including acetate, propionate, butyrate, and valerate, indicating that FML improve the amount of these SCFAs in the colonic contents. These SCFAs play pivotal roles in several host physiological functions, including nutrient utilization, energy expenditure, gut immunity, and macromolecular synthesis (Koh et al. 2016). For example, butyrate provides 60–70%

energy source for colonic epithelial cells and also protects against colorectal cancer and inflammation (Flint et al. 2012; Scheppach 1994). Butyrate-enriched high-fat diet in mice increased thermogenesis and energy expenditure, and the hosts are resistant to obesity (Gao et al. 2009). Acetate entering peripheral circulation can be metabolized by peripheral tissues and then transported to the liver for cholesterol synthesis (Wolever et al. 1989). Acetate dietary supplements in obese and diabetic rats improve glucose tolerance and reduce weight gain of animals (Yamashita et al. 2014). Most propionate is absorbed by the liver and is a good precursor for gluconeogenesis,
<table><caption>**Table 4** Comparison of mean relative abundance of bacterial genera among groups</caption><thead><tr><th>Genus (%)</th><th>Control</th><th>Treat 1</th><th>Treat 2</th><th>Treat 3</th></tr></thead><tbody><tr><td>Roseburia</td><td>1.058 ± 0.135a</td><td>1.308 ± 0.198a</td><td>2.602 ± 0.494ab</td><td>3.061 ± 0.621b</td></tr><tr><td>Butyricicoccus</td><td>0.166 ± 0.080a</td><td>0.148 ± 0.044a</td><td>0.479 ± 0.113b</td><td>0.145 ± 0.036a</td></tr><tr><td>Peptostreptococcaceae_unclassified genus</td><td>0.186 ± 0.069ab</td><td>0.233 ± 0.065a</td><td>0.061 ± 0.013ab</td><td>0.016 ± 0.054b</td></tr><tr><td>Dorea</td><td>0.080 ± 0.020ab</td><td>0.055 ± 0.017a</td><td>0.175 ± 0.049b</td><td>0.036 ± 0.016a</td></tr><tr><td>Clostridium</td><td>0.169 ± 0.012a</td><td>0.185 ± 0.021a</td><td>0.064 ± 0.031b</td><td>0.024 ± 0.017b</td></tr><tr><td>Desulfovibrio</td><td>0.023 ± 0.019a</td><td>0.007 ± 0.004a</td><td>0.02 ± 0.002a</td><td>0.245 ± 0.110b</td></tr><tr><td>Phascolarctobacterium</td><td>0.180 ± 0.11a</td><td>0.148 ± 0.079a</td><td>0.469 ± 0.080ab</td><td>0.732 ± 0.218b</td></tr><tr><td>SMB53</td><td>8.737 ± 3.998ab</td><td>12.069 ± 3.383a</td><td>2.858 ± 1.065ab</td><td>0.296 ± 0.051b</td></tr><tr><td>Oxalobacter</td><td>0.002 ± 0.002a</td><td>0.006 ± 0.006ab</td><td>0.001 ± 0.001a</td><td>0.042 ± 0.019b</td></tr><tr><td>Defluviitalea</td><td>0.047 ± 0.018a</td><td>0.056 ± 0.030ab</td><td>0.143 ± 0.019b</td><td>0.067 ± 0.020ab</td></tr><tr><td>Clostridiaceae_unclassified genus</td><td>5.549 ± 2.137a</td><td>6.798 ± 2.280a</td><td>2.318 ± 0.359b</td><td>0.303 ± 0.068b</td></tr><tr><td>Peptostreptococcaceae_unclassified genus</td><td>0.146 ± 0.068a</td><td>0.159 ± 0.041a</td><td>0.039 ± 0.013a</td><td>0.007 ± 0.010a</td></tr><tr><td>Prevotella</td><td>4.264 ± 2.183ab</td><td>2.973 ± 0.940a</td><td>9.349 ± 1.695ab</td><td>16.127 ± 5.704b</td></tr><tr><td>Peptococcus</td><td>0.083 ± 0.069a</td><td>0.053 ± 0.019ab</td><td>0.023 ± 0.007ab</td><td>0.012 ± 0.004b</td></tr></tbody></table>

Only those genera that were significant different genera (P < 0.05) among all groups are shown

Data was expressed as mean ± SE obtained from 6 individual pigs each. Significant difference is indicated by different letters between groups. One-way analysis of variance (one-way ANOVA) with Tukey's post hoc test was used to test the differences. P-values were corrected using false discovery rate (FDR). The group control, treat 1, treat 2, and treat 3 signify that feed is supplemented with 0, 5%, 10%, or 15% F fermented Mao-tai lees, respectively

<figure><img src="image_7.png" /><figcaption>**Fig. 3** Differences in the pig gut microbiota among the different fermented Mao-tai lees (FML) treatments. Linear discriminant analysis (LDA) effect size (LEFSe) results show that bacterial OTUs/markers were significantly different in abundance between control and FML-treated groups</figcaption></figure>

liponeogenesis, and protein synthesis (Wolever et al. 1991). Food containing propionate in healthy women increased insulin release and reduced the fasting glucose

level (Venter et al. 1990). These straight-chain fatty acids are indicators of the fermentation of dietary fiber, thus
<figure><img src="image_8.png" /><figcaption>**Fig. 4** The distribution of potential probiotics and pathogens of pigs among treatments. The mean relative abundance of these genera was normalized using Z-score transformation</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 5** The heatmap plot of spearman correlations between bacterial markers and SCFAs or bioamines. Only those correlations with *P* < 0.05 are shown</figcaption></figure>

FML possibly improve the fermentation ability of gut microbiota for dietary fiber.

BCFAs (e.g., isobutyrate and isovalerate) are metabolites that result from protein fermentation via gut bacteria (Birkett et al. 1996), and isobutyrate and isovalerate

originate from L-leucine and L-valine, respectively (Le et al. 2005). BCFA amounts in the gut contents could be regarded as indicators of protein catabolism in the colon (Blachier et al. 2007). Our results showed that FML additive amounts were not associated with BCFAs, indicating
that FML supplements likely failed to improve protein catabolism. Thus, our results indicate that the increased SCFA net concentration in the colon is likely the result of increased fermentation of specific dietary fibers by intestinal microbiota.

Although FML improve the concentrations of straight-chain fatty acids, higher gut SCFA concentrations in 15% FML could result from increased SCFA production or decreased SCFA absorption in the colon, or the fewer bacterial species in the gut utilize the SCFAs as an energy source. SCFA production and profiles in humans or animals were regulated by a lot of different host traits, environmental, dietary or microbial factors (Li et al. 2018; Macfarlane and Macfarlane 2003). The FML groups had higher contents of crude fiber and protein than the control group (Additional file 1: Table S1), thus one possible factor that controls specific SCFA production in our study is possibly nutrient components of available feed. We found that there are several OTUs that were significant different between the control and 15% FML groups and most of these OTUs were positively with the straight-chain fatty acids, thus gut microbiota composition may be another important factor that influences specific gut SCFA. Overall SCFA profile structure was significant influenced by body weight rather than by FML (Fig. 1), indicating that body weight may regulate the relative percentage of these SCFAs. In addition, those unmeasured factors, such as gut transit time, can also cause the difference of SCFA concentration (Cummings 1978) among different groups. Therefore, this may result in a limited understanding for the relationship between food fermentation time and SCFAs.

Bioamines are microbial fermentation products and are derived from amino acids (AA) metabolism. Bioamines can be produced based on AA decarboxylation pathways by specific bacteria (Tuberoso et al. 2015). Bioamines play important roles in various bacterial functions, such as maintaining growth and reproduction of normal cells and reducing bacterial susceptibility to host-derived antimicrobials (Shah and Swiatlo 2008). In addition, bioamines are correlated with numerous host diseases (e.g., psychiatric and neurologic disorders) (Stahl 1977), thus may be an indicator of host immunity and health. However, our results that most of bioamines were not correlated with FML content in feed, indicating that FML may have no impacts on the bioamines. Notably, tyramine was positively correlated with FML content. The possible reason is that FML improve the relative abundance of some tyramine-producing bacteria, such as *Lactobacillus* (De Las et al. 2006) which enriched in those FML groups.

## FML influence the beta diversity and metabolism-associated gene functions of pig gut microbiota

Our results showed that FML had no significant impacts in the alpha diversity, indicating that microbial species diversity is relatively stable in response to FML. However, FML influenced the beta diversity (Jaccard and Bray-Curtis dissimilarity) of pig gut microbiota compared with the control group, indicating that FML may influence species replacement (changes in species taxa) and species sorting (changes in abundance). Notably, only 15% FML was able to shape the community structure compared with the control group. These results suggested that the effect of FML on gut microbiota structure depends on additive amount. One possible reason is that a high-content dietary fiber in the 15% FML modified the gut microbiota structure, as demonstrated by previous study (Tap et al. 2015). However, our data showed that overall gene function profiles showed no difference among groups, indicating that overall functions of gut microbiota were more stable possibly for maintaining the gut performance. Compared with the control group, the 15% FML was enriched in some gene functions associated with metabolism (e.g., starch and sucrose metabolism and sulfur metabolism), thus 15% FML not only influence the gut microbiota structure but also function. In addition to dietary fiber, other nutrients, such as protein, may also impact the structure and function of gut microbiota (Scott et al. 2013).

## FML increase the abundance of several potential beneficial bacteria but decrease the abundance of specific pathogens

Our results showed that the phyla Firmicutes and Bacteroidetes accounted for 89.1% of total sequences, similar to previous findings in the microbiota of large intestines and feces in pigs (Bian et al. 2016; Kelly et al. 2017; Mach et al. 2015). At genus level, we found that *Prevotella* was the most dominant bacterial genera in the pig gut. *Prevotella* is known to produce various enzymes, such as glycoside hydrolases and polysaccharide lyase enzymes (Kaoutari et al. 2013), and the abundance changes of this genus is associated with dietary fiber (Liu et al. 2012). Thus, a higher abundance of *Prevotella* in the 15% FML may result from higher fiber content for 15% FML, so that the fiber ingredients improve the enrichment of this genus. *Prevotella* was also abundant in the gut microbiota from Burkina Faso children, whose diet ingredients contain a large amount of plant fiber (De Filippo et al. 2010). In fact, it has been demonstrated that *Prevotella* was capable of degrading plant cell wall dietary fiber, and thus may produce a great deal of SCFAs that can be utilized by hosts and microbes (De Filippo et al. 2010). In
addition to *Prevotella*, *Roseburia* was also enriched in the 15% FML group. *Roseburia* is butyrate-producing bacteria (Machiels et al. 2014), thus a higher abundance of this genus may lead to higher concentration of butyrate in the 15% FML group.

LEfSe analysis offered high-resolution discrimination at OTU level, and the results showed that those OTUs belonging to *Prevotella* (OTU25, OTU26, OTU31, and OTU36) and *Roseburia* (OTU94), *Ruminococcus bromii* (OTU117), *Treponema* (OTU132), *Phascolarctobacterium* (OTU121), and *Faecalibacterium* (OTU106), were significantly enriched in the 15% FML group. Some members of *Ruminococcus* and *Treponema* are involved in cellulose, lignin, and resistant starch degradation (Niu et al. 2015; Ze et al. 2012). *Faecalibacterium* is one of the most important commensal bacteria in the human gut microbiota. Several members of *Faecalibacterium*, as well as *Phascolarctobacterium* may produce various SCFAs (including acetate and butyrate) through the fermentation of dietary fiber (Lukovac et al. 2014), as demonstrated by our data, which found that most of these abundant OTUs in the 15% FML group were positively correlated with acetate and butyrate. However, it is still difficult to determine which bacteria taxa were responsible for the specific SCFA due to the complex interactions among bacteria, such as cross-feeding (Rey et al. 2010) and resource competition (Mahowald et al. 2009).

Accumulative evidence has demonstrated that some potential prebiotics may stimulate the growth of beneficial bacteria and be resistant against pathogenic bacteria. In this study, we compared the relative abundance of five potential probiotis and one pathogen between the control and treat groups. Our results showed that the potentially beneficial bacterial genera *Bacillus*, *Lactobacillus*, *Akkermansia*, and *Faecalibacterium* were more abundant in the 15% FML, while the potentially pathogenic genera *Escherichia* had a lower abundance. Some members of *Bacillus*, such as *Bacillus subtilis*, may reduce diarrhea scores and modulate microbial diversity in pigs (Bhandari et al. 2008). *Lactobacillus* provided many beneficial effects for pigs, such as improving ileal histomorphology, reducing systemic inflammatory cytokines, and increasing fermentation ability (Guerra-Ordaz et al. 2014). *Akkermansia*, especially *Akkermansia muciniphila*, is a mucin-degrading bacterium derived from mucus layer and shows a positive correlation with metabolic ability and health in mice (Everard et al. 2013). Some reports found that an increase in abundance on *Faecalibacterium* and *Akkermansia* was associated with a healthy gut status. In contrast, *Escherichia* is one known gut pathogens, and can lead to various diseases, such as diarrhea (Sack

1975). These results demonstrated that 15% FML had beneficial effects on the balance of pig gut microbiota.

In conclusion, this study is the first time to investigate the beneficial effects of FML on gut microbiota and metabolic profiles in pigs. Our results clearly demonstrated that 15% FML modulated the gut microbiota structure and function. Importantly, 15% FML increased the abundance of putative beneficial bacteria and decreased the abundance of potential pathogen *Escherichia* spp. in the pig gut, and also contributed to improve the fermentation of dietary fiber or indigestible polysaccharides in feed. Consequently, our results indicated that FML may be developed as a new feed additive for healthy and ecological livestock farming.

## Additional file

Additional file 1. Additional tables and figures.

### Authors' contributions

XFK, HWL and YLY designed the experiments. HL, FB and XFK wrote the manuscript text. HL performed the data analysis. PFX, ZHL, and HWL fed the animals. All authors contributed to the revision of the manuscript. All authors read and approved the final manuscript.

### Author details

¹Hunan Provincial Key Laboratory of Animal Nutritional Physiology and Metabolic Process, Key Laboratory of Agro-ecological Processes in Subtropical Region, National Engineering Laboratory for Pollution Control and Waste Utilization in Livestock and Poultry Production, Institute of Subtropical Agriculture, Chinese Academy of Sciences, Changsha, Hunan 410125, China. ²Institute of Occupational Health and Environmental Health, School of Public Health, Lanzhou University, Lanzhou, Gansu 730000, China. ³UMR PNCA, AgroParis-Tech, INRA, Université Paris-Saclay, Paris, France.

### Acknowledgements

We would like to thank Editage (http://online.editage.cn/) for English language editing.

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and materials

The original 16S rRNA data were available at the European Nucleotide Archive by Accession Number PRJEB27667 (http://www.ebi.ac.uk/ena/data/view/PRJEB27667).

### Ethics approval and consent to participate

Animal ethics approval for the present project was obtained from the Animal Ethics Committee of Institute of Subtropical Agriculture, Chinese Academy of Sciences.

### Funding

The present work was jointly supported by Grants from the National Key Project of Research and Development Plan (2017YFD0500503), the National Nature Science Foundation of China (31772613 and 31572421), the Fundamental Research Funds for the Central Universities (Izujbky-2018-68), the Changsha Key Project of Science and Technology Plan (kq1703007), and Talent Projects of Guangxi Science and Technology Department (AD17195043). All funding contributed to the design of study, collection, analysis, data interpretation and writing the manuscript.
**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 28 January 2019 Accepted: 4 February 2019

Published online: 18 February 2019

**References**

Abudabos AM, Al-Atiyat RM, Albatshan HA, Aljassim R, Aljumaah MR, Alkhulaifi MM, Stanley DM (2017) Effects of concentration of corn distillers dried grains with solubles and enzyme supplementation on cecal microbiota and performance in broiler chickens. *Appl Microbiol Biotechnol* 101(18):7017–7026

Ahmed S, Hoon J, Hong-Seok M, Chul-Ju Y (2014) Evaluation of *Lactobacillus* and *Bacillus*-based probiotics as alternatives to antibiotics in enteric microbial challenged weaned piglets. *Jr Afr Microbiol Res* 8(1):96–104

Bai Z, Ma L, Qin W, Chen C, Oenema O, Zhang F (2014) Changes in pig production in China and their effects on nitrogen and phosphorus use and losses. *Environ Sci Technol* 48:12742–12749

Baquero F, Martinez JL, Canton R (2008) Antibiotics and antibiotic resistance in water environments. *Curr Opin Biotechnol* 19(3):260–265

Barba-Vidal E, Martin-Orue SM, Castillejos L (2018) Review: are we using probiotics correctly in post-weaning piglets? *Animal* 12:1–10

Barros V, Lana G, Lana S, Lana Á, Cunha F, Neto J (2015) β-Mannanase and mannan oligosaccharides in broiler chicken feed. *Ciência Rural* 45:111–117

Bengtsson-Palme J, Hartmann M, Eriksson KM, Pal C, Thorell K, Larsson DG, Nilsson RH (2015) Metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data. *Mol Ecol Resour* 15:1403–1414

Bhandari SK, Xu B, Nyachoti CM, Giesting DW, Krause DO (2008) Evaluation of alternatives to antibiotics using an *Escherichia coli* K88+ model of piglet diarrhea: effects on gut microbial ecology. *J Anim Sci* 86(4):836–847

Bian G, Ma S, Zhu Z, Su Y, Zoetendal EG, Mackie R, Liu J, Mu C, Huang R, Smidt H, Zhu W (2016) Age, introduction of solid feed and weaning are more important determinants of gut bacterial succession in piglets than breed and nursing mother as revealed by a reciprocal cross-fostering model. *Environ Microbiol* 18(5):1566–1577

Birkett A, Muir J, Phillips J, Jones G, ODea K (1996) Resistant starch lowers fecal concentrations of ammonia and phenols in humans. *Am J Clin Nutr* 63:766–772

Blachier F, Mariotti F, Huneau JF, Tomé D (2007) Effects of amino acid-derived luminal metabolites on the colonic epithelium and physiopathological consequences. *Amino Acids* 33:547–562

Bray J, Curtis J (1957) An ordination of the upland forest communities of southern wisconsin. *Ecol Monogr* 27:325–349

Caporaso J, Kuczynski J, Stombaugh J (2010) QIIME allows analysis of high-throughput community sequencing data. *Nat Methods* 7:335–336

Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, Gormley N, Gilbert JA, Smith G, Knight R (2012) Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. *ISME J* 6(8):1621–1624

Cummings J (1978) Diet and transit through the gut. *J Plant Foods* 3:83–95

De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, Collini S, Pieraccini G, Lionetti P (2010) Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. *Proc Natl Acad Sci USA* 107(33):14691–14696

De Las RB, Marcobal A, Carrascosa AV, Muñoz R (2006) PCR detection of foodborne bacteria producing the biogenic amines histamine, tyramine, putrescine, and cadaverine. *J Food Prot* 69:2509–2514

De Maesschalck C, Eeckhaut V, Maertens L, De Lange L, Marchal L, Nezer C, De Baere S, Croubels S, Daube G, Dewulf J, Haesebrouck F, Ducatelle R, Taminau B, Van Immerseel F (2015) Effects of xylo-oligosaccharides on broiler chicken performance and microbiota. *Appl Environ Microbiol* 81(17):5880–5888

den Besten G, van Eunen K, Groen AK, Venema K, Reijingoud DJ, Bakker BM (2013) The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. *J Lipid Res* 54(9):2325–2340

Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. *Bioinformatics* 26:2460–2461

Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* 27(16):2194–2200

Everard A, Belzer C, Geurts L, Ouwerkerk J, Druart C, Bindels L, Guiot Y, Derrien M, Muccioli G, Delzenne N, de Vos W, Cani P (2013) Cross-talk between *Akkermansia muciniphila* and intestinal epithelium controls diet-induced obesity. *Proc Natl Acad Sci USA* 110:9066–9071

FlinT HJ, Scott KP, Louis P, Duncan SH (2012) The role of the gut microbiota in nutrition and health. *Nat Rev Gastroenterol Hepatol* 9(10):577–589

Gao Z, Yin J, Zhang J, Ward RE, Martin RJ, Lefevre M, Cefalu WT, Ye J (2009) Butyrate improves insulin sensitivity and increases energy expenditure in mice. *Diabetes* 58:1509–1517

Guerra-Ordaz AA, González-Ortiz G, La Ragione RM (2014) *Lactulose* and *Lactobacillus plantarum*: a potential complementary symbiotic to control post-weaning colibacillosis in piglets. *Appl Environ Microbiol* 80(16):4879–4886

Ilea RC (2009) Intensive livestock farming: global trends, increased environmental concerns, and ethical solutions. *J Agric Environ Ethics* 22:153–167

Jaccard P (1912) The distribution of the flora in the alpine zone. *New Phytol* 11:37–50

Ji Y, Guo Q, Yin Y, Blachier F, Kong X (2018) Dietary proline supplementation alters colonic luminal microbiota and bacterial metabolite composition between days 45 and 70 of pregnancy in Huanjiang mini-pigs. *J Anim Sci Biotechnol* 9:18

Jiang HY, Zhou YY, Zhou GL, Li YC, Yuan J, Li XH, Ruan B (2018) Gut microbiota profiles in treatment-naïve children with attention deficit hyperactivity disorder. *Behav Brain Res* 347:408–413

Kaoutari AE, Armougom F, Gordon JI, Raoult D, Henrissat B (2013) The abundance and variety of carbohydrate-active enzymes in the human gut microbiota. *Nat Rev Micro* 11(7):497–504

Kelly J, Daly K, Moran AW, Ryan S, Bravo D, Shirazi-Beechey SP (2017) Composition and diversity of mucosa-associated microbiota along the entire length of the pig gastrointestinal tract; dietary influences. *Environ Microbiol* 19(4):1425–1438

Koh A, De Vadder F, Kovatcheva-Datchary P, Backhed F (2016) From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. *Cell* 165(6):1332–1345

Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Thurber RLV, Knight R (2013) Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat Biotechnol* 31(9):814–821

Le PD, Aarnink AJA, Ogink NWM, Becker PM, Verstegen MW (2005) Odour from animal production facilities: its relationship to diet. *Nutr Res Rev* 18:3–30

Li F, Duan Y, Li Y, Tang Y, Geng M, Oladele OA, Kim SW, Yin Y (2015) Effects of dietary n-6:n-3 PUFA ratio on fatty acid composition, free amino acid profile and gene expression of transporters in finishing pigs. *Br J Nutr* 113(5):739–748

Li H, Li T, Beasley DE, Hedenec P, Xiao Z, Zhang S, Li J, Lin Q, Li X (2016a) Diet diversity is associated with beta but not alpha diversity of pika gut microbiota. *Front Microbiol* 7:1169

Li H, Li T, Yao M, Li J, Zhang S, Wirth S, Cao W, Lin Q, Li X (2016b) Pika gut may select for rare but diverse environmental bacteria. *Front Microbiol* 7:1269

Li H, Qu J, Li T, Li J, Lin Q, Li X (2016c) Pika population density is associated with composition and diversity of gut microbiota. *Front Microbiol* 7:758

Li H, Li T, Berasategui A, Rui J, Zhang X, Li C, Xiao Z, Li X (2017a) Gut region influences the diversity and interactions of bacterial communities in pikas (*Ochotona curzoniae* and *Ochotona daurica*). *FEMS Microbiol Ecol* 93(12):fix149

Li T, Li H, Gatesoupe FJ, She R, Lin Q, Yan X, Li J, Li X (2017b) Bacterial signatures of "Red-Operculum" disease in the gut of crucian carp (*Carassius auratus*). *Microb Ecol* 74(3):510–521

Li H, Qu J, Li T, Wirth S, Zhang Y, Zhao X, Li X (2018) Diet simplification selects for high gut microbial diversity and strong fermenting ability in high-altitude pikas. *Appl Microbiol Biotechnol* 102:6739–6751

Liu H, Ivarsson E, Dicksved J, Lundh T, Lindberg J (2012) Dietary inclusion of chicory (*cichorium intybus* L.) affects pig intestinal micro-environment and gut microbiota. *Appl Environ Microbiol* 78:4102–4109
Liu H, Li J, Guo X, Liang Y, Wang W (2018) Yeast culture dietary supplementation modulates gut microbiota, growth and biochemical parameters of grass carp. *Microb Biotechnol* 11(3):551–565

Lukovac S, Belzer C, Pellis L, Keijser BJ, de Vos WM, Montijn RC, Roeselers G (2014) Differential modulation by *Akkermansia muciniphila* and *Faecalibacterium prausnitzii* of host peripheral lipid metabolism and histone acetylation in mouse gut organoids. *mBio* 5(4):e01438

Macfarlane S, Macfarlane G (2003) Regulation of short-chain fatty acid production. *Proc Nutr Soc* 62:67–72

Mach N, Berri M, Estelle J, Levenez F, Lemonnier G, Denis C, Leplat JJ, Chevaleyre C, Billon Y, Dore J, Rogel-Gaillard C, Lepage P (2015) Early-life establishment of the swine gut microbiome and impact on host phenotypes. *Environ Microbiol Rep* 7(3):554–569

Machiels K, Joossens M, Sabino J, De Preter V, Arijs I, Eeckhaut V, Ballet V, Claes K, Van Immerseel F, Verbeke K, Ferrante M, Verhaegen J, Rutgeerts P, Vermeire S (2014) A decrease of the butyrate-producing species *Roseburia hominis* and *Faecalibacterium prausnitzii* defines dysbiosis in patients with ulcerative colitis. *Gut* 63(8):1275–1283

Mahowald MA, Rey FE, Seedorf H, Turnbaugh PJ, Fulton RS, Wollam A, Shah N, Wang C, Magrini V, Wilson RK, Cantarel BL, Coutinho PM, Henrissat B, Crock LW, Russell A, Verberkmoes NC, Hettich RL, Gordon JI (2009) Characterizing a model human gut microbiota composed of members of its two dominant bacterial phyla. *Proc Natl Acad Sci USA* 106(14):5859–5864

Markowiak P, Slizewska K (2018) The role of probiotics, prebiotics and synbiotics in animal nutrition. *Gut Pathog* 10:21

Marshall BM, Levy SB (2011) Food animals and antimicrobials: impacts on human health. *Clin Microbiol Rev* 24(4):718–733

Moura P, Marques S, Alves L, Freire JPB, Cunha LF, Esteves MP (2007) Effect of xylo-oligosaccharides from corn cobs autohydrolysis on the intestinal microbiota of piglets after weaning. *Livest Sci* 108(1–3):244–248

Nicola S, Jacques I, Levi W, Dirk G, Larisa M, Wendy SG, Curtis H (2011) Metagenomic biomarker discovery and explanation. *Genome Biol* 12:R60

Niu Q, Li P, Hao S, Zhang Y, Kim SW, Li H, Ma X, Gao S, He L, Wu W, Huang X, Hua J, Zhou B, Huang R (2015) Dynamic distribution of the gut microbiota and the relationship with apparent crude fiber digestibility and growth stages in pigs. *Sci Rep* 5:9938

NRC (2012) *Nutrient requirements of swine*. The National Academies Press, Washington, DC

Pourabedin M, Guan L, Zhao X (2015) Xylo-oligosaccharides and virginiamycin differentially modulate gut microbial composition in chickens. *Microbiome* 3:15

Rastall RA, Gibson GR (2015) Recent developments in prebiotics to selectively impact beneficial microbes and promote intestinal health. *Curr Opin Biotechnol* 32:42–46

Rey FE, Faith JJ, Bain J, Muehlbauer MJ, Stevens RD, Newgard CB, Gordon JI (2010) Dissecting the in vivo metabolic potential of two human gut acetogens. *J Biol Chem* 285(29):22082–22090

Rooks MG, Garrett WS (2016) Gut microbiota, metabolites and host immunity. *Nat Rev Immunol* 16(6):341–352

Sack RB (1975) Human diarrheal disease caused by enterotoxigenic *Escherichia coli*. *Ann Rev Microbiol* 29:333–354

Scheppach W (1994) Effects of short chain fatty acids on gut morphology and function. *Gut* 35(1):S35–S38

Scott KP, Gratz SW, Sheridan PO, Flint HJ, Duncan SH (2013) The influence of diet on the gut microbiota. *Pharmacol Res* 69(1):52–60

Shah P, Swiatlo E (2008) A multifaceted role for polyamines in bacterial pathogens. *Mol Microbiol* 68(1):4–16

Stahl SM (1977) The human platelet: a diagnostic and research tool for the study of biogenic amines in psychiatric and neurologic disorders. *Arch Gen Psychiatry* 34:509–516

Stanley D, Hughes RJ, Geier M, Moore RJ (2016) Bacteria within the gastrointestinal tract microbiota correlated with improved growth and feed conversion: challenges presented for the identification of performance enhancing probiotic bacteria. *Front Microbiol* 7:187

Sun D, Jia H (2015) China's pork imports in 2014 and the main pork trading market profile in the world Swine. *Ind Sci* 32(2):46–47

Tap J, Furet JP, Bensaada M, Philippe C, Roth H, Rabot S, Lakhdari O, Lombard V, Henrissat B, Corthier G, Fontaine E, Dore J, Leclerc M (2015) Gut microbiota richness promotes its stability upon increased dietary fibre intake in healthy adults. *Environ Microbiol* 17(12):4954–4964

Tremaroli V, Bäckhed F (2012) Functional interactions between the gut microbiota and host metabolism. *Nature* 489(7415):242–249

Tuberoso CI, Congiu F, Serreli G, Mameli S (2015) Determination of dansylated amino acids and biogenic amines in cannonau and vermentino wines by HPLC-FLD. *Food Chem* 175:29–35

Venter CS, Vorster HH, Cummings JH (1990) Effects of dietary propionate on carbohydrate and lipid metabolism in healthy volunteers. *Am J Gastroenterol* 85:549–553

Wang Q, Garrity GM, Tiedje JM, Cole JR (2007) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. *Appl Environ Microbiol* 73(16):5261–5267

Wolever T, Brighenti F, Royall D, Jenkins A, Jenkins D (1989) Effect of rectal infusion of short chain fatty acids in human subjects. *Am J Gastroenterol* 84:1027–1033

Wolever T, Spadafora P, Eshuis H (1991) Interaction between colonic acetate and propionate in humans. *Am J Gastroenterol* 84:1027–1033

Wu X, Ran Z, Kong X, Tang Z, Dong X (2010) Functions of microbiota in monogastric gastrointestinal tract and new practices in animal production. *J Food Agric Environ* 8:695–702

Wu Q, Chen L, Xu Y (2013) Yeast community associated with the solid state fermentation of traditional Chinese Maotai-flavor liquor. *Int J Food Microbiol* 166:323–330

Yamashita H, Fujisawa K, Ito E, Idei S, Kawaguchi N, Kimoto M, Hiemori M, Tsuji H (2014) Improvement of obesity and glucose tolerance by acetate in type 2 diabetic Otsuka long-Evans Tokushima fatty (OLETF) rats. *Biosci Biotechnol Biochem* 71(5):1236–1243

Yang S, Lio J, Wang T (2012) Evaluation of enzyme activity and fiber content of soybean cotyledon fiber and distiller's dried grains with solubles by solid state fermentation. *Appl Biochem Biotechnol* 167(1):109–121

Yu X, Abler D (2014) Where have all the pigs gone? Inconsistencies in pork statistics in China. *China Econ Rev* 30:469–484

Ze X, Duncan SH, Louis P, Flint HJ (2012) *Ruminococcus bromii* is a keystone species for the degradation of resistant starch in the human colon. *ISME J* 6(8):1535–1543

Zhao Y, Su JQ, An XL, Huang FY, Rensing C, Brandt KK, Zhu YG (2018) Feed additives shift gut microbiota and enrich antibiotic resistance in swine gut. *Sci Total Environ* 621:1224–1232