Dietary Inulin Supplementation Modulates Short-Chain Fatty Acid Levels and Cecum Microbiota Composition and Function in Chickens Infected With Salmonella

The current study investigated the effects of inulin on the gut microbiota, microbiome functions, and short-chain fatty acids (SCFAs) levels in specific pathogen-free (SPF) chickens infected with Salmonella enteritidis (SE). SPF Arbor Acres chickens (n = 240, 1-day-old) were divided into four groups: a control group (CON) fed a basal diet without inulin supplementation or SE infection, and three groups fed a basal diet supplemented with inulin 0, 0.5, and 1% (SE, 0.5%InSE, 1%InSE, respectively) up to 28-days-old, followed by SE challenge at 28 days of age. Cecal SCFA contents and microbiome composition and function were analyzed at 1-day post-infection. The results showed that SE infection significantly decreased cecal butyrate concentrations compared with the CON group (p < 0.05), while inulin supplementation reversed these changes compared with the SE group (p < 0.05). Inulin supplementation at 1% significantly increased the abundances of Lactobacillus and Streptococcus, and significantly decreased the abundances of Subdoligranulum and Sellimonas compared with the SE group (p < 0.05). The functional profiles of microbial communities based on metagenomic sequencing analysis showed that SE infection significantly increased the abundances of pathways related to carbohydrate metabolism, amino acid metabolism, energy metabolism, metabolism of cofactors and vitamins, and glycan biosynthesis and metabolism (p < 0.05), and significantly decreased the abundances of pathways related to nucleotide metabolism, translation, and replication and repair compared with the CON group (p < 0.05), and these effects were reversed by inulin supplementation (0.5 and 1%) (p < 0.05). In conclusion, inulin modulated the dysbiosis induced by SE infection via affecting SCFA metabolism and microbial functional profiles.


INTRODUCTION
Chicken meat is one of the most important and economical sources of animal protein for human consumption worldwide. However, the modern poultry industry uses large amounts of antimicrobials to prevent disease and enhance growth and productivity, resulting in antimicrobial resistance and a risk to human health (Shang et al., 2018). Numerous studies have demonstrated that many bacterial pathogens isolated from poultry meat, such as Escherichia coli and Salmonella, possess resistance genes (Abraham et al., 2019). Salmonella Enteritis (SE) infection in poultry leads to reduced growth and lower economic benefits, and poses a threat to food safety and human public health as a result of possible foodborne outbreaks in humans (Nair et al., 2018). Antibiotic-based growth promoters will consequently be totally banned in China in 2020. Alternatives to antibiotic growth promoters in animal feed are therefore receiving increasing attention, with the aim of improving meat quality and safety, while promoting animal productivity and welfare (Markowiak andŚliżewska, 2018).
Inulin-type fructans are one of the best-studied prebiotics and are used to improve poultry growth performance and meat quality (Guaragni et al., 2019), enhance nutrient utilization and immune function (Xia et al., 2019), regulate gut microbiota-related metabolism disorders, and alleviate intestinal inflammation (Singh et al., 2019). On one hand, inulin can be recognized directly by carbohydrate receptors on the surface of the intestinal epithelial cell and immune cells (such as dendritic cells), which consequently modulate the host immune response and inhibit pathogen colonization Teng and Kim, 2018). Inulin has been reported to downregulate the gene expression of proinflammatory cytokines, such as tumor necrosis factor-α (TNF-α), interferon-γ, and interleukin (IL)-6, and inhibit pro-inflammatory immunity (Teng and Kim, 2018), which provides metabolic advantages to gut commensal bacteria to compete with Salmonella (Sun et al., 2018). On the other hand, inulin has indirect regulatory effects that modulate dysbiosis and the host immune response based on microbiota-dependent mechanisms (Sivaprakasam et al., 2016). Generally, inulin stimulates the growth and activity of lactic acid bacteria and can also be fermented by gut commensal bacteria, serving as a nutrient source for the bacteria, which generate short-chain fatty acids (SCFAs), mainly acetate, propionate, and butyrate. These in turn suppress infection by acid-sensitive pathogens through creating an acidic environment (Teng and Kim, 2018). The regulatory mechanisms of SCFAs have been studied extensively, including their functions in nutrition, energy metabolism, microbiota regulation, histone deacetylase inhibition, lysine deacetylase inhibition, G-protein coupled receptor activation, and acylation of bacterial virulence regulators . However, most studies in poultry have focused on microbial metabolic regulation, especially the proliferation of beneficial microbiota such as Bifidobacterium and Lactobacillus, and the inhibition of colonization by pathogenic bacteria, such as Salmonella (Micciche et al., 2018).
The addition of inulin may thus provide a strategy for reducing the abundance of SE via a microbiota-dependent immunomodulatory mechanism. Indeed, SE infection potentially changes gut microbial communities in poultry, resulting in gut dysbiosis. SE challenge decreased the abundance of SCFAproducing Clostridia, Bifidobacterium, and Lactobacillus, leading to decreased SCFA levels, attenuating the intestinal innate immune response, and increasing gut inflammation (Rivera-Chávez et al., 2016). Inulin also increased the concentration of SCFAs, which in turn enhanced the abundance of the major beneficial bacteria Bifidobacteria and Lactobacillus in the cecum, which competed with pathogenic bacteria and reduced the SE populations (Micciche et al., 2018). The composition of the intestinal microbiota may play a fundamental role in the defensive mechanism of inulin against SE infection. However, the cecum is the most densely populated microbial habitat in the chicken gut (up to 10 11 cells/g), but its microbial composition remains largely unknown (Stanley et al., 2013). Therefore, the current study aimed to use metagenomics approaches to determine how SE infection affected SCFA production and functional microbial profiles in the chicken cecum, and to clarify how inulin modulated the dysbiosis induced by SE infection via SCFA metabolism and microbial functional profiles. These results will provide new opportunities for improving the gut microecosystem balance and reducing the presence of pathogens to maintain gut health and food safety.

MATERIALS AND METHODS
The study protocol was approved by the Animal Welfare Committee of the Institute of Animal Sciences (Chinese Academy of Agricultural Sciences, Beijing, China), Ethical approval regarding animal survival was given by the animal ethics committee of IAS-CAAS (approval number: IASCAAS-AE20140615).

Animal Feeding and Management
A total of 240, 1-day-old Arbor Acres specific pathogen-free (SPF) chicks were acquired from a local hatchery (Merial Vital Laboratory Animal Technology Co., Ltd., Beijing, China). All birds were housed in a sterile isolation chamber (IPQ-Type 3 negative pressure isolator; China Agricultural University, Beijing, China) and fed a corn/soybean-based sterile diet on the day of arrival. The birds were randomly assigned to four groups (six replicates of 10 birds each): a control group (CON) fed a basal diet, and three groups fed with basal diet supplemented with different doses of inulin (0, 0.5, and 1%) up to 28 days of age, followed by challenge at day 28 with SE (SE, 0.5%InSE, and 1%InSE groups, respectively). The inulin with DP of >10 was obtained from Cichorium intybus (Sigma-Aldrich, St. Louis, MO, United States). Details of the basal diet ingredients and calculated nutrient contents are provided in Supplementary Table 1. The temperature in the isolator was maintained at 35 • C for the first week and then decreased by 2 • C each week until the end of the experiment (day 29). All chickens had free access to feed and water (sterilized at 121 • C for 15 min).

SE Challenge
Salmonella enteritidis (CMCC50041, China Institute of Veterinary Drugs Control, Beijing) was obtained from the China Center of Industrial Culture Collection. Briefly, chickens in the SE, 0.5%InSE, and 1%InSE groups were infected with 1.0 mL of an actively growing culture of SE of 1 × 10 8 colony-forming units/mL by oral gavage at 28 days of age. Birds in the CON group received an equal volume of phosphate-buffered saline.

Cecal Sample Collection
At 1-day post-infection, one chicken from each replicate was randomly selected and sacrificed. The cecum on both sides was opened using sterile scissors and tweezers, and the cecal contents were squeezed out into two frozen tubes and stored at −80 • C for SCFA and DNA extraction for 16S rDNA and metagenomics analysis.

SCFA Concentration Analysis
Cecal SCFA concentrations were measured by gas chromatography-mass spectrometry. Briefly, cecal contents (50 mg) were homogenized in 200 µL of 24% metaphosphoric acid using a vortex mixer for 15 min and centrifuged at 5,000 × g for 5 min at 4 • C. The supernatants were analyzed using a DB-FFAP column (Agilent Technologies, Santa Clara, CA, United States) and mass-selective flame ionization detector (FID) with the following parameters: injector volume, 2 µL; split ratio, 50:1; injector temperature, 250 • C; FID temperature, 280 • C; and carrier gas, helium. The initial oven temperature was 70 • C for 10 min, increased by 5 • C/min to 210 • C, and then maintained for 12 min. The acetic acid, propionic acid, and butyric acid peaks in a standard solution (Sigma-Aldrich, MO, United States) were analyzed using the same parameters. The peaks of individual SCFAs in each cecal sample were acquired using the same parameters. The molar concentration of each SCFA was calculated using the ratio of the peak area of the individual SCFA and the peak area of the standard solution multiplied by the concentration of the standard solution.
High-Throughput Sequencing of the Microbial 16S rRNA Gene The cecal microbiota was determined with 16S rRNA sequencing analysis as described previously (Tan et al., 2019). Briefly, bacterial DNA was extracted from 22 cecal digested samples (5 from CON and SE group respective and 6 from 0.5%InSE, and 1%InSE groups respective) using an E.Z.N.A. R Stool DNA Kit (Omega Bio-tek, Norcross, GA, United States) according to the manufacturer's protocols. The DNA concentration and purification were determined using TBS-380 and a NanoDrop 2000 UV-vis spectrophotometer (Thermo Scientific, Wilmington, DE, United States), respectively. DNA quality was examined by 1% agarose gel electrophoresis. The V3-V4 hypervariable regions of the bacterial 16S rRNA gene were amplified with primers 338F (5 -ACTCCTACGGGAGGCAGCAG-3 ) and 806R (5 -GGACTACHVGGGTWTCTAAT-3 ) using a thermocycler polymerase chain reaction (PCR) system (GeneAmp 9700, ABI, CA, United States). The resulting PCR products were extracted from a 2% agarose gel and further purified using an AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) and quantified using QuantiFluor TM -ST (Promega, WI, United States) according to the manufacturer's protocol. Purified amplicons were pooled in equimolar amounts and subjected to paired-end sequencing on an Illumina MiSeq platform (Illumina, San Diego, CA, United States) according to standard protocols, by Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China).

Processing of Sequencing Data and Diversity Analysis
Raw fastq files were demultiplexed and quality-filtered to obtain high quality clean tags according to QIIME (V1.7.0; San Diego, CA, United States) with the following criteria (Magoč and Salzberg, 2011): (i) 300 bp reads were truncated at any site with an average quality score <20 over a 50 bp sliding window, discarding truncated reads <50 bp; (ii) exact barcode matching, a two nucleotide mismatch in primer matching and reads containing ambiguous characters were removed; and (iii) only sequences with >10 bp overlap were assembled according to their overlap sequence. Reads that could not be assembled were discarded. Operational taxonomic units (OTUs) were clustered with 97% similarity cut-off using UPARSE (version 7.1 1 ), and chimeric sequences were identified and removed using the UCHIME algorithm. The taxonomy of each 16S rRNA gene sequence was analyzed by RDP Classifier 2 against the Silva (SSU115) 16S rRNA database, with a confidence threshold of 70%. Alpha and beta diversity were analyzed based on these output normalized data. Chao1, Simpson, Shannon, and ACE indexes were calculated to assess the alpha-diversity in each sample using QIIME (V1.7.0; San Diego, CA, United States) and displayed using R software (Version 3.3.1, R Core Team, Vienna, Austria). Beta diversity based on unweighted UniFrac was performed to differentiate among the samples in terms of species complexity using QIIME software (V1.7.0; San Diego, CA, United States). Principal coordinates analysis (PCoA) was applied to reduce the dimension of the original variables and executed with QIIME and R software using the unweighted pair-group method with arithmetic means.

Library Construction and Metagenomics Sequencing
The cecal microbiota was further investigated using a metagenomic sequencing method as described previously . DNA was extracted from 12 samples (3 from each group respective) for metagenomics analysis, using an E.Z.N.A. R Stool DNA Kit (Omega Bio-tek, Inc., Norcross, GA, United States) according to the manufacturer's protocols. DNA from each sample was then fragmented to an average size of approximately 300 bp using Covaris M220 (Gene Company Limited, Beijing, China) for paired-end library construction, using a TruSeq TM DNA Sample Prep Kit (Illumina Inc., San Diego, CA, United States). Adapters containing the full complement of sequencing primer hybridization sites were ligated to the blunt-end fragments. Paired-end sequencing was performed on an Illumina HiSeq4000 platform (Illumina Inc., San Diego, CA, United States) at Majorbio Bio-Pharm Technology Co., Ltd., using a HiSeq 3000/4000 PE Cluster Kit and HiSeq 3000/4000 SBS Kits, according to the manufacturer's instructions 3 .

Sequence Quality Control and Genome Assembly
Reads were aligned to the Gallus gallus genome (version 5.0) by BWA 4 and hits associated with the reads and their paired reads were removed. The 3 and 5 ends were stripped using SeqPrep 5 . Low quality reads (length < 50 bp, quality value < 20, or having N bases) were removed using Sickle 6 . De Bruijngraph-based assembler SOAP de novo ( 7 Version 1.06) was employed to assemble short reads (Li et al., 2008). K-mers, varying from 1/32/3 of read length, were tested for each sample. Scaffolds > 500 bp were retained for statistical tests. The quality and quantity of scaffolds generated by each assembly were evaluated and the K-mer that yielded the minimum scaffold number and the maximum values of N50 and N90 was finally selected. Scaffolds > 500 bp were extracted and broken into contigs without gaps and the contigs were used for further gene prediction and annotation.

Gene Prediction, Taxonomy, and Functional Annotation
Open reading frames (ORFs) from each metagenomic sample were predicted using Metagene 8 . Predicted ORFs ≥ 100 bp were retrieved and translated into amino acid sequences using the NCBI translation table 9 . All sequences from gene sets with a 95% sequence identity (90% coverage) were clustered as the non-redundant gene catalog using CD-HIT 10 . After quality control, reads were mapped to the representative genes with 95% identity using SOAP aligner (see text footnote 7), and the gene abundance in each sample was evaluated. Taxonomic annotations were carried out using BLASTP (Version 2.2.28+ 11 ) by aligning non-redundant gene catalogs against the NCBI NR database, with an e-value cut-off of 1e −5 . Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was conducted using BLASTP search (Version 2.2.28+) against the KEGG database 12 , with an e-value cut-off of 1e −5 . Based on the comparison results, functional annotation was performed using KOBAS 2.0 (KEGG language-based annotation system). The corresponding gene abundance sum of KEGG orthologies (KOs) and KEGG were then used to calculate the abundance of the corresponding functional category.

Statistical Analysis
Data are shown as mean ± standard deviation. Differences in SCFA contents, alpha diversity indices, relative abundance of bacterial taxa among treatment groups, and bacterial functional pathways were analyzed using SAS 9.2 (SAS Institute, Cary, NC, United States), followed by Duncan's multiple comparison tests. Differences between two groups were analyzed using Student's t-test. A significant difference was declared when p < 0.05. Spearman's correlations between SCFA contents and the top 15 significantly different genera of cecal microbiota were assessed using the PROC CORR procedure in SAS 9.2. A correlation matrix was created and visualized in a heatmap format using R package version 3.3.1. Cecal microbiota abundance and SCFA concentration were considered to be correlated with each other when the absolute value of the correlation coefficient (r) was > 0.55 and p < 0.05. PCoA, hierarchical clustering analysis, and heat map analyses were conducted using R package version 3.3.1.

Inulin Affects the Production of SCFA
Cecal digesta SCFA levels are shown in Table 1. SE challenge significantly decreased the butyrate content and increased the acetate and propionate contents (p < 0.05) of the cecal digesta compared with the CON group, while inulin supplementation (0.5 and 1%) significantly increased acetate and butyrate concentrations (p < 0.05), and 1% inulin supplementation decreased the propionate concentration compared with the SE group (p < 0.05).

Analysis of the Microbial 16S rRNA High-Throughput Sequencing
The 16S rRNA gene sequencing was performed to analyze the cecal microbial community structures in the different groups. After assigning OTUs and removing chimeras, there were 35,000-50,000 sequencing reads in each sample (mean  Figure 1). Good's coverage index (Supplementary Table 2) showed that sufficient sequencing coverage was acquired. Alpha diversity richness indexes (Sobs, Chao1, ACE) and diversity indexes (Shannon and Simpson) of the cecal microbiota are shown in Table 2. The Shannon index was significantly higher in the CON and inulin supplementation groups than in the SE group (p < 0.05), and the Sobs, ACE, and Chao1 indexes were significantly higher in inulin supplementation groups than in the CON group (p < 0.05).
The beta diversity analysis was measured by unweighted UniFrac distances and illustrated in a PCoA plot, which was used to estimate the similarity among different groups (Figure 1). Similar to PCoA, clustering analysis revealed short UniFrac distances between the two inulin groups and long distances between the samples and the CON, and SE and inulintreated groups (Supplementary Figure 2). The dissimilarity in cecal microbial communities was greatest between the SE and inulin groups.

Functional Capacity Profiling of Cecum Microbiome Based on Metagenomic Sequencing
The functional capacity of the gut microbiome was investigated using metagenomic sequencing data. Based on the annotation of ORFs predicted from the assembled contigs, a total of 8,608,302 ORFs with an average length of 595.64 bp was found to infer with the functional capacity. The predicted genes were classified by aligning them to the KEGG database. A total of 6725 KOs were identified and assigned to 402 KEGG pathways. A total of 1341 differential KOs between different treatments were analyzed, and the relative abundances of the top 30 KOs showed different enrichments (Figure 5)   The cecal microbiota showed differential abundance at the phylum level among four groups according to one-way ANOVA. (C) The cecal microbiota showed differential abundance at the phylum level between the CON and SE groups according to t-tests. (D) The cecal microbiota showed differential abundance at the phylum level between the SE and 1%InSE groups according to t-test. The cecal microbiota showed differential abundance at the genus level between CON and SE group according to t-test. (C) The cecal microbiota showed differential abundance at the genus level between the SE and 0.5%InSE groups according to t-test. (D) The cecal microbiota showed differential abundance at the genus level between the SE and 1%InSE groups according to t-test.  (Figure 5). These changes were also reversed by inulin supplementation compared with the SE group (p < 0.05). We further identified the top15 KEGG microbial functions at level 2 of the KEGG pathways terms showing distinct enrichments between inulin-supplemented chickens with and without SE infection ( Table 3). The abundances of pathways related to carbohydrate metabolism, energy metabolism, metabolism of cofactors and vitamins, and glycan biosynthesis and metabolism were significantly increased in the SE group compared with the CON group (p < 0.05), while pathways related to nucleotide metabolism, translation, and replication and repair were significantly decreased (p < 0.05). However, the abundance of amino acid metabolism pathways was similar in the CON and SE groups. Inulin supplementation (0.5 and 1%) significantly reversed the changes in abundances of these KEGG pathways (p < 0.05). We also identified the top 30 differently enriched KEGG microbial functions and compared the KEGG pathways at level 3 ( Figure 6 and Table 4). The abundances of 12 KEGG functional pathways, including one pathway in membrane transport ("Bacterial secretion system"), one in carbohydrate metabolism ("Propanoate metabolism"), one in glycan biosynthesis and metabolism ("Peptidoglycan biosynthesis"), five in replication and repair ("Homologous recombination, " "DNA replication, " and "Mismatch/Base excision/Nucleotide excision"), two in translation ("Aminoacyl-tRNA biosynthesis" and "Ribosome"), and two pathways in nucleotide metabolism ("Purine/Pyrimidine metabolism") were significantly reduced in the SE group compared with the CON group, while inulin addition significantly increased the abundances of these KEGG pathways compared with the SE group (p < 0.05). However, a total of 18 KEGG pathways, including seven pathways in carbohydrate metabolism ("Starch and sucrose/Galactose/Fructose and mannose/C5-Branched dibasic acid metabolism, " "Citrate cycle (TCA cycle), " and "Pentose and glucuronate interconversions"), one in glycan biosynthesis and metabolism ("Glycosphingolipid biosynthesis"), six pathways in amino acid metabolism ("Alanine. aspartate and glutamate/Histidine/Arginine and proline/Phenylalanine/Tyrosine metabolism/Valine, leucine and isoleucine degradation"), two in energy metabolism ("Nitrogen/Sulfur metabolism"), and two in metabolism of cofactors and vitamins ("Pantothenate and CoA/Folate biosynthesis") were over-represented in the SE group compared with the CON group (p < 0.05), but were significantly reduced in the 0.5%InSE and 1%InSE groups compared with the SE group (p < 0.05).

DISCUSSION
The cecum is an important region for microbial activity in the chicken gut, partly because the greatest density of microbial communities (up to 10 11 cells/g) occurs here, the digesta remains in the cecum for a relatively long time (10-12h), and microorganisms reside here longer than in the small intestine (Hong et al., 2019). Moreover, the cecal microbiome plays important nutritional and immune roles, such as nitrogen recycling, providing vitamins, amino acids, and SCFA to their host, especially by supplying energy to the host and competitively inhibiting pathogenic infections (Sood et al., 2019). In the current study, we observed that the cecal microbiome was mainly composed of Firmicutes, Bacteroidetes, and Actinobacteria at the phylum level, in accordance with other studies (Ocejo et al., 2019). Although some researchers also showed that Proteobacteria and Tenericutes were common phyla in the chicken cecum, and the abundances of these phyla were low in the current study. Among the six most abundant genera, only Bifidobacterium belonged to the phylum Actinobacteria, while Faecalibacterium, Bacteroides, Lactobacillus, Ruminococcus, and Streptococcus all belonged to the Firmicutes phylum, showing that the predominant cecal bacteria were Firmicutes. The chicken cecal microbiome has been found to be affected by various factors, including diet, sex, genotype, feed additives, housing conditions, and pathogen infection (Clavijo and Flórez, 2018). Notably, pathogenic infections, such as Clostridium perfringens, Eimeria species, and SE, resulted in significant changes in the cecal microbiota (Juricova et al., 2013). In the present study, the abundances of Firmicutes, Actinobacteria, and Tenericutes phyla were reduced by SE infection; however, the abundance of Firmicutes (mainly Lactobacillus and Streptococcus) was increased by inulin supplementation. Bifidobacterium and Lactobacillus are considered as beneficial bacteria in many animal species, because of their ability to produce natural bacteriocines. Lactobacillus has also been reported to produce organic acids (lactic and acetic acids) in poultry, to promote growth performance and combat infection by pathogens, such as Salmonella, Shigella, Clostridium, and Listeria (Yadav and Jha, 2019). Bifidobacterium also acts as an immunostimulant, as well as competing with pathogenic bacteria for cell adhesion sites and producing essential SCFAs for energy production by the host (Yan et al., 2017). Inulin-type fructan, which is an important dietary fiber, can modulate the gut microbial population, especially by increasing the abundances of Bifidobacterium and Lactobacillus and inhibiting pathogen proliferation (E. coli and Salmonella) by promoting a low pH in the cecum of broiler chickens (Mesa et al., 2017). We previously showed that dietary inulin supplementation protected the gut mucosa by improving immune responses, SCFA production, and intestinal morphology and suppressing the colonization of SE in SPF broiler chickens (Song et al., 2020). However, the effects of inulin on the intestinal microbiota in chickens requires further investigation.
In the present study, inulin increased the abundance of Lactobacillus, but 0.5% inulin decreased the abundance of Bifidobacterium compared with SE infection. Xia et al. (2019) also showed that inulin lowered the relative abundances of Lactobacillus and Bifidobacterium species at an early age, but subsequently increased their relative abundances. Moreover, as a commensal organism in chickens, the genus Streptococcus played FIGURE 6 | Hierarchical clustering analysis and heat map analysis on the relative abundance of KEGG microbial function at level-3 pathways in the CON, SE, 0.5%InSE and 1%InSE groups. Red, high KO abundance; blue, low KO abundance. CON (n = 3): basal diet without inulin supplementation; SE (n = 3): Salmonella-infected chicken fed with basal diet; 0.5% InSE (n = 3): Salmonella-infected chicken fed with 0.5% inulin; 1% InSE (n = 3): Salmonella-infected chicken fed with 1% inulin. a vital role in reducing the colonization and incidence of SE (Adhikari et al., 2019). In this study, inulin also increased the abundance of Streptococcus. Inulin may therefore competitively inhibit Salmonella colonization and strengthen the immune system partly because of its ability to increase the proliferation of beneficial microbiota, such as Bifidobacteria, Lactobacilli, and Streptococcus, in the gastrointestinal tract (Micciche et al., 2018).
Short-chain fatty acids, including acetate, propionate, and butyrate, are the end-products of fermentation of inulin by gut microbes and have been extensively investigated for their roles in ameliorating the intestinal mucosal barrier and suppressing intestinal inflammation via several different mechanisms (Yang et al., 2016). SE infection decreased the concentration of butyrate, possibly as a result of a reduction in the abundance of butyrate-producing Clostridia induced by SE (Rivera-Chávez et al., 2016). The present results also showed that acetate and butyrate were the main cecal SCFAs that contributed to the acidic environment in the gut, to suppress pathogen growth (Teng and Kim, 2018). Inulin supplementation increased acetate and butyrate levels in the present study, in accordance with our observations on the changes in SCFAs (Song et al., 2018). Butyrate was reported to be the main energy source for the proliferation of intestinal epithelium cells and to increase villus height (De Vadder et al., 2014). Moreover, an increase in butyrate may increase oxygen consumption by colonocytes, leading to inhibition of the proliferation and expansion of aerobic bacteria, such as SE Litvak et al., 2018). Numerous reports have shown that SCFAs, especially butyrate, stimulate the growth of Bifidobacteria and Lactobacilli in poultry, which in turn promote an increase in SCFAs, lowering the pH and enhancing the mucosal immune function and inhibiting the growth of pathogenic bacteria such as Salmonella and E. coli (Micciche et al., 2018).
In the present study, Clostridiales, including Sellimonas, Peptococcaceae, Lachnospiraceae, and Subdoligranulum, had a negative relationship with SCFA production, while vadinBB60, Bacteroides, and Defluviitaleaceae_UCG-011 had positive relationships with SCFA production. Moreover, vadinBB60 was positively related to propionate production, and some reports showed that it was also associated with butyrate production (Van den Abbeele et al., 2010). In accordance with previous studies, the current study found that inulin increased the abundances of Streptococcus and Lactobacillus, which were shown to be positively related to butyrate production (Yan et al., 2017). Streptococcus, Lactobacillales, and Firmicutes exhibited a positive relationship with acetate production in the present study, playing an important role in reducing the colonization and incidence of SE through reducing the pH and competitive exclusion (Adhikari et al., 2019). The functional profiles of microbial communities based on metagenomic sequencing analysis showed that the most abundant functions were carbohydrate metabolism, amino acid metabolism, nucleotide metabolism, membrane transport, and energy metabolism with significant differences in level 2 KEGG pathways among the four treatments. We also showed that the abundances of functional pathways related to membrane transport (bacterial secretion system) and replication and repair (mismatch repair, DNA replication) were reduced, while amino acid metabolism, carbohydrate metabolism, and energy metabolism were increased in the SE group compared with the CON group, and there were significant differences in the level 3 KEGG pathways between the SE and CON groups. Membrane transport pathways, such as bacterial secretion, are essential for protein secretion, enhancing cell attachment, and for the survival of bacteria in the gut ecosystem (Green and Mecsas, 2016). Pilmis et al. (2020) demonstrated that the bacterial secretion system would be related to defending endogenous bacteria (Bacteroidetes and Proteobacteria) against the overgrowth of certain exogenous pathogenic bacteria species, such as SE. In the current study, the abundance of the bacterial secretion system pathway, including the associated KOs (K16785, K15580), was significantly reduced in the SE group compared with the CON group. Similarly, Mon et al. (2020) showed that the bacterial secretion system was down-regulated in SE-infected chickens, pointing toward a possible downshift in the ability of intestinal microbes for adhesion and attachment of mucosal cells. However, inulin supplementation significantly increased the abundance of the bacterial secretion system and the related KOs. These results suggested that inulin might play an important role in affecting the gut microbiota composition in birds challenged with SE through the microbial functional metabolism.
The present results also showed that amino acid metabolism and energy metabolism at level 3, including the related KOs (K00266, K01652), were enriched in the SE-infection group, indicating that pathways associated with nutrients and energy metabolism were activated to meet the nutrient and energy requirements of the cecal microbiome in SE-infected chickens. However, the gut microbiota shifts toward the utilization of less favorable residual peptides and proteins to gain energy (proteolytic fermentation) when SCFAs are insufficient to meet their energy needs for gut microbial growth and proliferation because of SE infection (Makki et al., 2018;Canfora et al., 2019). SE infection resulted in a decrease in butyrate production in the present study, which contributed to enrichment of amino-acid-metabolism-related pathways involving "Phenylalanine metabolism, " "Tyrosine metabolism, " and "Valine, leucine and isoleucine degradation" in the SE group. Notably, microbial proteolytic fermentation can produce branched-chain amino acids (e.g., valine, leucine, and isoleucine) and aromatic amino acids (e.g., tyrosine, phenylalanine, and tryptophan), which can be further metabolized by microbiota to meet the energy requirement (Canfora et al., 2019). We therefore speculated that microbial fermentation of amino acids may also provide sources for energy metabolism (nitrogen and sulfur metabolism), which could explain why energy metabolismrelated pathways involving "nitrogen metabolism" and "sulfur metabolism" were enriched in the SE-infection group. Studies have therefore suggested that infected chickens should be supplied with the required energy by providing mixtures of dietary fibers (Bortoluzzi et al., 2017;Canfora et al., 2019). Some studies also showed that prebiotics, such as fructooligosaccharide and xylooligosaccharide, had an inhibitory effect on proteolysis (Lecerf et al., 2012). These results therefore suggested that inulin might decreasing the availability of amino acids for proteolytic bacteria, thus reducing ammonia production . Microbial metabolites derived from the fermentation of amino acids or proteins, mainly ammonia, amines, indole compounds, and hydrogen sulfide, are potentially detrimental to gut integrity, which could be responsible for SE-infectionrelated impairment of the mucosal barrier (Song et al., 2020). In the present study, we also found that fermentation of inulin increased butyrate and acetate production, which may decrease the energy sources from the fermentation of amino acids or protein for utilization by bacteria (Canfora et al., 2019). All these results indicated that the main strains responsible for proteolytic fermentation are Bifidobacterium, Clostridium, Enterococcus, Lactobacillus, Pediococcus, and Streptococcus, which may in turn be depleted by metabolizing the amino acids (Pugin et al., 2017), possibly helping to explain the decreases in abundances of Bifidobacterium, Clostridium, Lactobacillus, and Streptococcus in the SE group and the increases in inulin supplementation group in this study.

CONCLUSION
The results of this study demonstrated that dietary supplementation of inulin reversed the decrease in the alpha diversity (Shannon index) and richness (Chao1 index) induced by SE-infection. Moreover, Inulin increased butyrate and acetate production, which may decrease the energy sources from the fermentation of amino acids or protein for utilization by microbe. Accordingly, the enrichments of metabolic pathways associated with amino acid metabolism and energy metabolism were reduced and the abundances of Bifidobacterium, Clostridium, Lactobacillus, and Streptococcus responsible for proteolytic fermentation were increased, which may be the main cause of inulin improving the cecal microbiome of SE-infected chickens.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI PRJNA645511 and PRJNA645563.

ETHICS STATEMENT
The animal study was reviewed and approved by the Animal Welfare Committee of the Institute of Animal Sciences (Chinese Academy of Agricultural Sciences, Beijing, China), Ethical approval regarding animal survival was given by the animal ethics committee of IAS-CAAS (Approval number: IASCAAS-AE20140615).

AUTHOR CONTRIBUTIONS
JS formulated the research question, performed the study, and drafted the manuscript. QL and RL contributed to suggestion of data analysis and interpretation. MZ and GZ contributed to the schedule the study. NE and JW contributed to study design and modify the manuscript. All authors contributed to the article and approved the submitted version.