Abstract

The main purpose of the present study was to assess the effect of soluble and insoluble fiber on colonic bacteria and intestinal barrier function in a piglet model. A total of 24 piglets (25 ± 1 d old; 7.50 ± 0.31 kg) were randomly allotted to 4 treatments: basal diet (control, CON), 1% insoluble dietary fiber (IDF) diet, 1% soluble dietary fiber (SDF) diet, and 0.5% insoluble fiber + 0.5% soluble dietary fiber (MDF) diet. The trial lasted 28 days. SDF-fed piglets showed a higher bacterial a-diversity (observed_species, chao1, and ACE) and a higher relative abundance of Proteobacteria and Actinobacteria, Solobacterium, Succinivibrio, Blautia, and Atopobium in colonic digesta than CON, IDF, and MDF groups . At the same time, Bacteroidetes, Euryarchaeota, Phascolarctobacterium, Coprococcus_1, and Prevotella_1 were significantly increased in the IDF group when compared with CON, SDF, and MDF groups . Furthermore, Bacteroidetes and Enterobacteriaceae, Selenomonas, Phascolarctobacterium, and Alloprevotella were significantly higher in the MDF group than those in the other three groups . SDF diet increased the concentrations of short-chain fatty acid (SCFA) in colonic digesta when compared with the CON group and enhanced weight index of the colon than the CON and IDF groups. Furthermore, compared with the CON group, SDF, IDF, and MDF diets all upregulated the mRNA expressions of claudin-1 (CLDN-1) in colonic mucosa , SDF and IDF diets upregulated the mRNA expressions of mucin 2 (MUC2) , SDF diet increased mRNA expressions of zonula occludens 1 (ZO-1) and occludin (OCLN), while the IDF group enhanced the secretory immunoglobulin A (sIgA) concentrations , respectively. IDF and MDF diets decreased expressions of TNF-α . We concluded that the influence of soluble fiber on colonic microbiota was more extensive than that of insoluble fiber. Moreover, soluble fiber could more effectively improve colonic barrier function by upregulating gene expressions of the gut barrier.

1. Introduction

In western countries, colonic disorders are serious health issues [1]. Finding methods to maintain colon health in individuals are of current interest [2]. A healthy colon involves a symbiotic balance among the gut microbiota, the integrity of the intestinal barrier, and minimizing inappropriate inflammatory responses [3, 4]. The controversies exist between dietary fiber intake and colonic disorders [57]. Due to the fermentation in the hindgut, the effect of dietary fiber on hindgut health involved the mutual bacteria and the formation of short-chain fatty acid (SCFA) [8, 9].

Dietary fiber was reported to stimulate health-promoting bacteria (Bifidobacterium, Lactobacillus, Prevotella, etc.) and suppress pathogenic bacterial species (Escherichia coli, etc.) [8, 10]. Short-chain fructooligosaccharide intake increased the abundance of cecal Akkermansia and Blautia [11]. Moreover, fructan and cellulose led to a difference in UniFrac distances of fecal microbiota and a-diversity [12]. Different composition and physicochemical properties of fiber may lead to variations in gut microbiota and SCFA production [13, 14] Dietary fiber has been classified into soluble fiber and insoluble fiber according to the solubility [15] and the influence of two fibers on colonic microbial composition and metabolism, and colonic barrier still needs to be further investigated.

We hypothesized that soluble fiber and insoluble fiber caused different effects on colonic microbiota and regulated colonic barrier. Inulin is a soluble fiber made of fructans, which can be rapidly fermented in the colon. ARBOCEL (a crude fiber concentrate), an insoluble fiber with high water-holding capacity, is made of lignocellulose. These two fibers were selected as supplementary dietary fibers.

All piglets were fed diets for 28 days, and then microbial composition, the content of SCFAs, and the barrier function of the colon were determined. The results could potentially provide some new leads toward understanding the effects of soluble fiber and insoluble fiber on microbial groups and barrier functions in the colon of nonruminant animals and humans.

2. Materials and Methods

2.1. Animals and Experimental Design

All experimental procedures were approved by the Animal Welfare Committee of Sichuan Agricultural University and performed in accordance with the National Research Council’s Guide for the Care and Use of Laboratory Animals.

Twenty-four 25-day-old piglets (Duroc × Landrace × Yorkshire) with body weight (BW) of 7.50 ± 0.31 kg were randomly allotted to four groups with 6 replicates (3 males and 3 females). The four dietary treatments included CON diet (basal diet, maize-soybean meal), and fiber groups including IDF diet (1% insoluble fiber diet), SDF diet (1% soluble fiber diet), and MDF diet (0.5% insoluble fiber + 0.5% soluble fiber diet).

2.2. Animal Feeding and Sample Collection

Experimental diets were formulated to meet or exceed the nutrient requirement recommended by NRC (2012). Ingredient composition of the basal diet is presented in Table 1. 1% maize in the basal diet was replaced by 1% ARBOCEL, 1% inulin, and 0.5% ARBOCEL +0.5% inulin, respectively, in the three fiber groups. ARBOCEL was provided by German J.RETTENMAIER & Söhne Group (Shanghai, China), and inulin was supplied by Ci Yuan Biotech Company Limited (Shanxi, China). The concentrations of inulin and ARBOCEL were 92.15% and 90%, respectively (provided by the manufacturer). The experiment was conducted at the Research Base of the Institute of Animal Nutrition of Sichuan Agricultural University. All pigs were housed individually in metabolism cages (1.5 m × 0.7 m × 1.0 m). The room lighting was natural with the temperature maintained at 26–28°C and relative humidity controlled at 60%–70%. Piglets were fed 4 times daily at 08 : 00, 12 : 00, 16 : 00, and 20 : 00 and had free access to water. The experiment lasted 28 days.

On 29 d, piglets were sacrificed. Digesta from colon was collected to keep in a sterile tube and then frozen at −80°C for the analysis of SCFAs and bacterial community. The colon section was weighed after digesta removed. Mucosal scrapings from the colon were prepared and stored at −80°C to detect gene expression of barrier function.

2.3. Microbial DNA Extraction and Sequencing

Total bacterial DNA of colonic digesta was extracted from each sample by using CTAB/SDS method. 1% agarose gels were used for monitoring the concentration and purity of DNA. Following monitoring, DNA was diluted to 1 ng/μl using sterile water according to its concentration. The bacterial 16S rRNA gene amplification (V3–V4 fragments) was conducted using the barcoded primer pair 343F/806R set (343F: TACGGRAGGCAGCAG, 806R: GGACTACHVGGGTWTCTAAT). All PCRs (30 μL) were carried out with 15 μL of Phusion® High-Fidelity PCR Master Mix (New England Biolabs), 0.2 μM of forward and reverse primers, and about 10 ng template DNA. Thermal cycling reactions were performed by the following cycle program: initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, with a final elongation at 72°C for 5 min. All PCR products were purified using the electrophoresis in agarose gels and SanPrep DNA Gel Extraction Kit (Sangon Biotech, Shanghai, China). Samples with a bright main strip between 400 and 450 bp were chosen for further experiments. All PCR products were mixed in equal density ratios. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. At last, the library was sequenced on an Illumina HiSeq platform. All reads were deposited in the National Center for Biotechnology Information (NCBI) and can be accessed in the Short Read Archive (SRA) under accession number PRJNA493943.

2.4. Bioinformatics Analysis

Paired-end reads from the original DNA fragments are merged by using Fast Length Adjustment of SHort reads (FLASH)—a very fast and accurate analysis tool which is designed to merge paired-end reads when there are overlaps between reads 1 and reads 2 [16]. Paired-end reads were assigned to each sample according to the unique barcodes. Chimeric sequences were removed using USEARCH based on the UCHIME algorithm [17]. The microbial diversity was analyzed using Quantitative Insights Into Microbial Ecology (QIIME) software package [18]. Sequences were assigned to the same Operational Taxonomic Unit (OTU) with a 97% similarity threshold. Alpha diversity was determined based on the number of observed species, Shannon index, Simpson index, Chao1, ACE, goods_coverage, and PD_whole_tree. QIIME calculates both weighted and unweighted UniFrac, which are phylogenetic measures of beta diversity. We used weighted UniFrac for principal coordinate analysis (PCoA). PCoA helps to get principal coordinates and visualize them from complex, multidimensional data.

2.5. qPCR Analysis of Bacterial Groups

Quantitative detection of the relative abundance of Bacteroidetes, Firmicutes, Roseburia, Prevotella, and Ruminococcus [19] in all samples was performed by real-time PCR using SYBR Premix Ex Taq reagents (TaKaRa Biotechnology (Dalian), China). The primers (Table 2) and amplification program were followed as the methods mentioned by Bergström et al. [19]. A reaction was run in a volume of 11 μl with 5.5 μl 2 × SYBR Green PCR Master Mix, 0.4 μL of each primer (100 nmol/L), 2.7 μL nuclease-free water, and 2 μL template DNA. The universal bacterial reference primer set was selected for calculating the abundance of target bacterial groups.

2.6. Detection of SCFA Concentrations

The SCFA concentrations from colonic digesta were evaluated using gas chromatography. Each sample (1 g) was thawed and suspended in 2 ml distilled water in a screw-capped tube. After 30 min at 4°C, the tubes were centrifuged at 5,000 × g for 10 min. 1 mL supernatant was transferred to a new sterile tube and mixed with 0.2 mL 25% metaphosphoric acid and 23.3 μL 210 mmol/L crotonic acid. After stood for 30 min, the tubes were centrifuged at 10,000 × g for 10 min, and 300 μL supernatant was transferred to another sterile tube, and then 900 μL methanol was added. After centrifuged at 10,000 × g for 10 min, 100 μL supernatant was transferred to a sterile tube. The SCFAs (acetate, propionate, and butyrate) were quantified in a gas chromatographic system (VARIAN CP-3800, America).

2.7. ELISA Analysis of Secretory IgA Concentration

According to a 1 : 9 ratio (g/mL), mucosa from the colon was homogenized using physiological saline. The homogenate was centrifuged at 3,500 × g and 4°C for 10 min. Then, the supernatant was taken for the determination of sIgA using commercially available ELISA kit (Chenglin, Beijing).

2.8. RT-qPCR Analysis for Gene Expression

According to the manufacturer’s instructions, total RNA from the colonic mucosa was isolated using TRIzol reagent (TaKaRa Biotechnology (Dalian), China). The yield and purity of total mRNA were measured using a spectrophotometer (Beckman Coulter DU 800, Beckman Coulter Inc, Brea, USA), and an optical density 260: 280 (OD260: OD280) ranging from 1.8 to 2.0 was considered a very low degree of contamination [20]. The integrity of RNA was analyzed by agarose gel electrophoresis. The RNA samples were reversely transcribed into complementary DNA using RT Reagents (TaKaRa Biotechnology (Dalian), China) according to the manufacturer’s instructions. After reverse transcription, gene expressions of zonula occludens 1 (ZO-1), occludin (OCLN), claudin 1 (CLDN-1), mucin 1 (MUC1), mucin 2 (MUC2), interleukin-10 (IL-10), interleukin-1β (IL-1β), and tumor necrosis factor-α (TNF-α) in colon were analyzed by real-time quantitative PCR using SYBR Premix Ex Taq reagents (TaKaRa Biotechnology (Dalian), China) and CFX-96 Real-Time PCR Detection System (Bio-Rad Laboratories, Richmond, CA) as described by Mao [21]. The primers (Table 3) were purchased from Invitrogen (Shanghai, China). Cycling conditions were performed as previously described [21]. Each sample was determined in triplicate on the same PCR plate, and the mean values were used for the statistical analysis. Relative gene expression to the reference gene (β-actin) was used for normalization, and the relative mRNA expression levels of the target gene in comparison with the reference gene were calculated by the 2–ΔΔCT method [22].

2.9. Statistical Analysis

Each piglet was considered as the statistical unit. All data were analyzed as a randomized complete block design using the MIX of SAS (SAS Inst. Inc., Cary, NC). Differences in the weight index of the colon, the relative abundance of certain bacterial phyla, families, or genera, alpha diversity, the concentrations of SCFAs, and gene expressions among treatments were analyzed with one-way ANOVA if the data were in line with a normal distribution (sig > 0.05). Once the results were not in line with normal distribution, they were analyzed by a nonparametric test of significance. The microbiome field results were analyzed using false discovery rate (FDR) to correction with a . The results were presented as means ± SE. was considered a significant difference.

3. Results

3.1. Sequence Analysis

A total of 1,935,292 raw reads with an average of 80,637 ± 10,162 sequences per sample were generated in 24 samples. After removing the low-quality sequences, 1,372,259 valid sequences with an average length of 418 bp were obtained. A total of 1815 OTUs were obtained, which could be identified to 26 bacterial phyla and 312 bacterial genera.

3.2. Microbial Diversity in Different Groups

Alpha diversity, be expressed as observed_species, Shannon’s diversity index, Simpson index, Chao1, ACE, goods_coverage, and PD_whole_tree, was calculated in microbial diversity analysis within the community (Table 4). Chao1 and ACE in the SDF group were the highest among all groups . Observed_species was higher in the SDF group than in the IDF and MDF groups .

Venn diagram showed the shared and unique OTUs among different groups (Figure 1(a)). There were 1099, 921, 1577, and 1126 OTUs in the CON, IDF, SDF, and MDF groups, respectively. The CON group had 26 unique OTUs, the IDF group had 39 unique OTUs, the SDF group had 489 unique OTUs, and the MDF group had 60 unique OTUs. According to PCoA, the colonic bacterial community structures of the SDF group were obviously separated from other groups (Figure 1(b)).

3.3. Composition and Abundance of Bacterial Communities in Different Groups

At the phylum level, a total of 26 phyla were detected in all samples. The most predominant bacterial phyla were Firmicutes and Bacteroidetes, which accounted for 47.9% and 40% of sequences, respectively (Figure 2(a)). The relative abundance of Bacteroidetes in the IDF and MDF groups was significantly higher than that in the SDF group . The relative abundance of Proteobacteria in the SDF group and Euryarchaeota in the IDF group was the highest among all groups, respectively . SDF diet increased the relative abundance of Actinobacteria compared with the CON group .

The most relatively abundant bacterial families were Prevotellaceae and Veillonellaceae. SDF-fed piglets had a lower relative abundance of Prevotellaceae than the CON, IDF, and MDF groups . Piglets from the MDF group had a lower relative abundance of Lachnospiraceae than those from the CON group and had a higher relative abundance of Enterobacteriaceae than those from the IDF group (Table 5 and Figure 2(a)).

A total of 312 different genera were detected. The most relatively abundant genera included Prevotella_9, Succinivibrio, Selenomonas, Alloprevotella, Megasphaera, Dialister, Prevotellaceae_NK3B31_ group, Prevotella_7, Olsenella, and Streptococcus (Figure 2(b)). The relative abundance of Dialister, Pseudobutyrivibrio, Ruminococcus_2, Ruminiclostridium, and Clostridium_sensu_stricto_6 was significantly higher in the CON group than in the IDF, SDF, and MDF groups. Phascolarctobacterium, Coprococcus_1, Prevotella_1, Ruminococcaceae_UCG-008, Leeia, and Treponema_2 were significantly abundant in the IDF group than in the CON, SDF, and MDF groups . Feeding SDF diets increased the relative abundance of Blautia, Solobacterium, Syntrophococcus, Olsenella, Atopobium, Succinivibrio, and Weissella when compared with the CON, IDF, and MDF groups. Selenomonas, Sharpea, Alloprevotella, unidentified_Veillonellaceae, and Phascolarctobacterium were significantly increased in the MDF group when compared with the CON, IDF, and SDF groups (Table 6 and Figure 2(b)).

3.4. Relative Abundance of Specific Bacteria in Different Groups Using qPCR

As shown in Table 7, the relative abundance of Bacteroidetes in the IDF group was higher than that in the SDF group . The CON group had the highest relative abundance of Prevotella .

3.5. Short-Chain Fatty Acid Concentrations in Different Groups

The SCFA concentrations in the colonic digesta are presented in Figure 3. Compared with the CON group, SDF diet increased the concentrations of total SCFAs, acetate, propionate, and butyrate , while IDF diet increased acetate concentration .

3.6. Weight Index and Gene Expressions of Colon in Different Groups

The weight index of colon in the SDF group was higher than that in the CON and IDF groups (Figure 4). The effects of dietary fibers on colonic gene expressions and sIgA concentration are presented in Figure 5. Compared with the control group, SDF diet increased mRNA expressions of ZO-1, CLDN-1, OCLN, and MUC2 while decreased mRNA expression of TNF-α . An increase in CLDN-1 and MUC2 mRNA levels and sIgA concentration was observed in the IDF group when compared with the CON group . Compared with the CON group, MDF diet increased mRNA expression of CLDN-1 and decreased mRNA expression of TNF-α .

4. Discussion

Dietary fiber has been shown to be degraded in the hindgut of animals and influenced the abundance and diversity of intestinal microbiota [23, 24]. Changes of fiber components led to a change in the composition of the microbiota [25, 26]. A remarkably higher relative abundance of the phylum Bacteroidetes and Euryarchaeota and the genus Prevotella, Phascolarctobacterium, Ruminococcaceae, Coprococcus, Leeia, and Treponema was found in the IDF group. SDF diet increased the relative abundance of the phylum Proteobacteria and Actinobacteria and the genus Blautia, Solobacterium, Syntrophococcus, Weissella, Olsenella, Atopobium, and Succinivibrio. MDF diet increased the relative abundance of the phylum Bacteroidetes and the genus Selenomonas, Phascolarctobacterium, Sharpea, and Alloprevotella. All these results indicated that different types of fiber could selectively regulate intestinal bacteria. It was reported that a significant increase in Actinobacteria mostly resulted in an increase in Bifidobacterium [27]. However, there was no significant difference and we did not find an increase in Bifidobacterium although the abundance of Actinobacteria was increased in the SDF group, which may result from the low relative abundance of Bifidobacterium in the colon of pigs.

The α-diversity index (observed_species, chao1, and ACE) of colonic bacteria was significantly increased in the SDF group when compared to other groups. The results were similar with a study in mice showing an increase in bacterial diversity after β-glucan supplementation [28]. These results indicated that the supplementation of soluble fiber but not insoluble fiber may increase the diversity of colonic microbes. Firmicutes, Bacteroidetes, Proteobacteria, and Actinobacteria were the most predominant phyla in all piglets, which were consistent with previous studies in pigs and humans [2931]. Inulin was reported to increase the abundance of Actinobacteria and decrease Bacteroidetes in vitro fermentation, while cellulose increased the abundance of Bacteroidetes and decreased Firmicutes [27]. In the present study, SDF-fed pigs showed a higher abundance of Actinobacteria and a lower abundance of Bacteroidetes; meanwhile, IDF diet increased the abundance of Bacteroidetes and decreased the abundance of Firmicutes. Although the change of Firmicutes was not significant, these results showed that the two fibers have consistent effects on microorganisms in vitro and in vivo.

In our results, IDF diet increased the concentration of acetate, while SDF diet increased the concentration of total SCFAs (acetate, propionate, and butyrate). The amounts of SCFAs in the colon digesta depended on several factors such as the composition of microbiota and types of fiber. The amounts of acetate and propionate correlate positively with Bacteroidetes and genus Blautia within Firmicutes [32, 33]. Higher relative abundances of Bacteroidetes were found in the IDF and MDF groups, while animals from the SDF group showed a higher abundance of Blautia. However, a higher concentration of acetate was only found in piglets fed IDF and SDF diets. The reason may be that feeding MDF diet promoted the growth of non-acetate-producing bacteria. SDF-fed piglets had a higher concentration of butyrate and the abundance of Actinobacteria, which confirmed the study that Actinobacteria produced high amounts of colonic butyrate production [34]. Meanwhile, the production of SCFAs was related to the composition of dietary fiber available for bacteria [35]. All saccharide composition of dietary fiber can be utilized for acetate formation and thus increases the acetate concentration [36]; in addition, butyrate can be produced by the fermentation of fructans [37]. The fermentation of soluble carbohydrates leads to a large amount of propionate production [38, 39]. The concentration of SCFAs (acetate, propionate, and butyrate) in the IDF and SDF groups was consistent with previous studies, which suggested that fermentation patterns of dietary fiber were closely related to the diversity of bacterial community and composition of fiber.

The intestinal barrier is consisted of tight junction proteins (ZO-1, CLDN1, and OCLN), the mucus, and immunological components like sIgA [32, 40]. Previous studies showed that dietary fiber proves intestinal barrier function in humans and animals [41]. High fiber diet increased the capacity of mucin secretion in the gastrointestinal tract [42]. MUC2 is the main secretory mucin in colon [43], and thus, IDF and SDF diets upregulated MUC2 mRNA level, suggesting an increased mucin secretion. Arabinoxylan in wheat was reported to increase intestinal sIgA concentrations in weaned piglets, in line with our results in animals fed IDF, which may be ascribed to a low concentration of toxic products [44]. Tight junction proteins (ZO-1, CLDN1, and OCLN) are highly organized structures that maintain an effective barrier against the invasion of harmful substances [45]. In the present study, different fiber groups did not have the same qualitative or quantitative effects on colonic barrier functions, while SDF diet had a better effect than other diets. Furthermore, SDF and MDF diets might promote barrier function by downregulating TNF-α gene expression since TNF-α was reported to reduce the tight junction protein expression [46]. As fiber is the main substrate for bacterial fermentation, it might regulate mucosal barrier function by supporting more diversified bacterial communities and increasing concentrations of SCFAs [47, 48]. SDF diet increased the microbial diversity and SCFA concentrations in colon, which might be the reason that SDF diet improved the barrier functions more effectively than other groups.

Research studies have increasingly suggested that probiotic bacteria play an important role in regulating gut barrier integrity [49]. Coprococcus was reported to regulate immune responses presumably through the production of IgG [50]. A higher concentration of sIgA found in the IDF group was probably also associated with the increase of Coprococcus. Prevotellaceae was increased in colorectal cancer patients while Blautia and Phascolarctobacterium were reduced [51]. Atopobium was considered to be useful for human health since an inverse correlation between its number and inflammatory bowel disease [52]. In the current study, SDF diet significantly increased the alpha diversity of colonic bacteria and the relative abundance of Blautia and Atopobium and decreased the relative abundance of Prevotellaceae. IDF and MDF diets increased the relative abundance of Phascolarctobacterium. Meanwhile, all fiber groups upregulated gene expressions of ZO-1, OCLN, CLDN-1, MUC1, and MUC2 and downregulated gene expressions of IL-1β and TNF-α in the colon. The results suggested that different fiber supplementation regulated gut barrier function by stimulating the growth of different bacterial species. Probiotic mixture protected the epithelial barrier and increased the OCLN and ZO-1 expression by activating the p38 and ERK signaling pathways, while reversed the effects of TNF-α [53], which suggested that dietary fiber might regulate barrier function by p38 and ERK signaling pathways. L. acidophilus could activate a pathogen-associated molecular pattern receptor, Toll-like receptor 2 (TLR2) in intestinal epithelial cell lines, and enhance the phosphorylation of NF-κB p65 and p38 mitogen-activated protein kinase (MAPK), which indicated another possible mechanism [54]. However, the underlying mechanism needs to be verified by further studies.

The SCFAs (acetate, propionate, and butyrate) produced as end metabolites by the microbiota were reported to improve gut barrier function [55]. An increase in acetate concentration in the IDF and SDF groups might be the partial reason for the increase of MUC2 mRNA level, since acetate has been shown to prevent inflammatory bowel diseases by inducing mucin secretion in mucin-deficient mice [56]. Significant positive correlations between colonic propionate concentrations and TFF expression were observed in rats [43]. Diet supplementation with butyrate inhibited the disruption of the intestinal epithelial barrier induced by high-fat diet via upregulating the gene expression of CLDN-1 [57]. In the present study, the highest concentration of SCFAs might be one of the reasons that SDF diet was more efficient in regulating colonic barrier function than other groups.

5. Conclusions

In summary, different types of fibers had different effects on the colonic barrier function by selectively modulating bacteria and SCFAs. Insoluble fiber like cellulose increased the relative abundance of Bacteroidetes, Euryarchaeota, Phascolarctobacterium, and Coprococcus, while soluble fiber like inulin stimulated Actinobacteria, Proteobacteria, Blautia, and Atopobium. Furthermore, feeding soluble fiber led to a higher concentration of SCFAs, microbial diversity, and community richness than insoluble fiber and then helped to improve the intestinal barrier function.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

The authors thank the following technical personnel for their kind assistance: Huifen Wang, Junning Pu, and Runqi Fu, who assisted in animal care-taking, sampling, and lab analysis. The authors also thank Dr. Hui Yan for revising the manuscript. This study was financially funded by the National Key Research and Development Program of China (2018YFD0500605).