Pretreatment of Rapeseed Meal Increases Its Recalcitrant Fibre Fermentation and Alters the Microbial Community in an in Vitro Model of Swine Large Intestine


 Background: The aim of current study was to investigate whether degradation of rapeseed meal (RSM) which was modified by a cellulase, two pectinase, or alkaline treatment was improved by the swine gut microbiota compared to untreated RSM, and whether the microbiota composition was changed. Methods: An in vitro study was performed to assess how enzymatic and chemical pretreated rapeseed meal (RSM) influences the fibre fermentation and microbial community in the swine large intestine. RSM was processed enzymatically by a cellulase (CELL), two pectinases (PECT), or chemically by an alkaline (ALK) treatment. 16S rRNA gene sequencing data was performedto evaluate changes in the gut microbiota composition, whereas short chain fatty acid production (ion-chromatography) and non-starch polysaccharides (NSP) composition(using monoclonal antibodies; mAbs) were used to assess fibre degradation.Results: The results showed that ALK, CELL, PECT1, and PECT2changed microbial community composition, increased the abundance of microbial fibre-degrading enzymes and pathways,andincreased acetic acid, propionic acid, butyric acid, and total SCFA production. The increased genera also positively correlated with SCFA production. The cell wall polysaccharide structuresof RSM shiftedafter ALK, CELL, PECT1, and PECT2 treatment. The degradation of NSPduring the fermentation period was dynamic, and not continuous based on the epitope recognition by mAbs.Conclusion: This study provides the first detailed analysis of changes in the swine intestinal microbiota due to RSM modified by ALK, CELL, PECT1 and PECT2. ALK, CELL, PECT1 and PECT2 altered microbial community structure, shifted the predicted functional metagenomic profile and subsequently increased total SCFA production. Our findings that ALK, CELL, PECT1 and PECT2increased fiber degradability in RSM could help guide feed additive strategies to improve efficiency and productivity in swine industry. The current study gave insight into how feed enzyme modulate microbial status, which provides good opportunity to develop novel carbohydrase, particularly in swine feed.


Introduction
Rapeseed meal (RSM), a by-product of rapeseed oil production, is not only a suitable protein source for swine feed but also a potentially energy source. RSM contains 20 to 40% non-starch polysaccharides (NSP) (Simbaya et al. 1995;Slominski and Campbell 1990). The primary cell walls of RSM consist of pectin and xyloglucan and its cellulose micro brils are interlinked with xyloglucan via hydrogen bonds forming a stiff network (Carpita and Gibeaut 1993). Pectins are linked to each other and cross-linked between pectin and hemicellulose, and between pectin and cellulose (Broxterman and Schols 2018). Pectins, consisting of homogalacturonan, rhamnogalacturonan, xylogalacturonan, arabinogalactan and arabinan, are the major polysaccharides present in the dehulled rapeseed meal (Eriksson, Andersson, and Åman 1997). In the secondary cell wall of RSM, the main carbohydrates are 4-O-methyglucuronoxylan, xyloglucan, and cellulose. A drawback of using RSM in animal feed is that the complex cell wall polysaccharides cannot be utilized by endogenous enzymes from monogastric animals (e.g. pigs), and also can only partly be fermented by the microbial community in the gastrointestinal tract (GIT) of the pig.
Therefore, the animal feed sector seeks opportunities to enhance degradability of NSP of feedstuffs, in order to improve its potential as a nutrient source for domestic animals. Previous research showed that physical processing technologies, such as hammer milling, pelleting, wet-milling, extrusion, and mild hydrothermal acid treatment, had limited effect on recalcitrant NSP structures (de Vries et al. 2013; De . As a result, more e cient solutions are needed to modify the cell wall architecture and allow the gut microbiota to utilize the complex carbohydrates. Former research has shown that NSPdegrading enzymes, such as cellulase and pectinases, could open the cell wall structure and improve NSP degradability (Giraldo et al. 2008; Giraldo et al. 2007). Previous studies showed that addition of pectolytic enzymes improved degradability of NSP of RSM in vitro ) and in broilers (Pustjens et al. 2014;De Vries et al. 2014). However, none of the above studies investigated how the gut microbiota was affected by the modi ed RSM or by the feed enzymes. It is important to know this, as NSP can only be fermented by microbes. Bindelle et al. (Bindelle et al. 2011) demonstrated that NSP-degrading enzymes increased abundances of cellulolytic Ruminococcusand xylanolytic Clostridium-like bacteria and altered fermentation patterns of barley cultivars and wheat products. Torok and colleagues (Torok et al. 2008) investigated changes in gut microbial population in response to the supplementation of an NSP-degrading enzyme (containing β-glucanase, xylanase, and protease activities) in a barley-based diet in chickens, and the results showed that microbial composition revealed distinct clusters correlating with un-supplemented and enzyme supplemented birds. Previous research reported that the pretreatment of feed stuffs with carbohydrases can cause the release of reducing sugars and other hydrolysis materials, promoting chemotactic response in speci c bacteria, and stimulating their attachment to feed particles, and thereby growth of these microbes (Giraldo et al. 2007; Beauchemin et al. 2003; Ribeiro et al. 2015; Ribeiro et al. 2018).
In the present study, RSM (predigested with digestive enzymes) was treated independently with two kinds of pectinases (PECT1 and PECT2), one cellulase (CELL), or alkaline (ALK), and afterwards the untreated and treated RSM preparations were fermented in the Swine Large Intestine in vitro Model (SLIM) (Long, de Vries, and Venema 2020). The aim of the current study was to investigate whether fermentation by the swine gut microbiota of treated RSM was improved compared to untreated RSM, and whether the microbiota composition and activity were changed.

Substrate preparation
Rapeseed meal (Brassica napus, Cargill N.V., Antwerp, Belgium; 2011) was obtained from a commercial feed mill (Agri rm B.V., Utrecht, the Netherlands). Preparation method I (predigestion of RSM after carbohydrase or alkaline treatment) [ Figure 1]: to 200 g of RSM 40 mL 10*gastric electrolyte concentrate solution (GES, 310 g sodium chloride, 110 g potassium chloride, 15 g calcium chloride di-hydrate, and 4840 g ultrapure water) and 360 mL ultrapure water were added. The pH was adjusted to 5.5 and then nothing (CON), 10 mL of alkaline (ALK, 6 M NaOH), or the following carbohydrases were added CELL (Accellerase 1000, Sigma-Aldrich, Missouri, United States), PECT1 (Pectinex Ultra SP, Novozymes A/S, Bagsvaerd, Denmark), or PECT2 (Multifect Pectinase, DuPont Industrial Biosciences, Genencor division, Rochester, NY). Enzyme preparations were incubated at 37 ℃ for 2 h, with occasional shaking (every 30 min), while ALK was incubated overnight at 4 ℃. Enzyme preparations were then heated at 100 ℃ for 5 min to inactive enzymes. Afterwards, for all ve samples, 120 mL GES was added and pH adjusted to 3 to continue with the gastric incubation according to the predigestion protocol as described elsewhere (Sáyago-Ayerdi, Zamora-Gasga, and Venema 2017). After predigestion, the slurry was centrifuged (8.000 g, at 4 ºC, for 20 minutes) and dialysis was performed for the supernatants. For dialysis, a dialyzer (Sure ux, Nipro Europe Group Companies, Mechelen, Belgium) was used with a peristaltic pump to remove small digestion products and water. After reduction of the total volume to ~ 450-500 mL, supernatant was mixed with pellet, and freeze-dried. Method II (predigestion of RSM before carbohydrase or alkaline treatment) [ Figure 1]: four batches of 200 g RSM were predigested as described before (Sáyago-Ayerdi, Zamora-Gasga, and Venema 2017) and then dialyzed. Afterwards, 55 mL 10*GES was added, and pH adjusted to 5.5, after which 10 mL of CELL, PECT1, PECT2, or ALK treatment commenced, respectively. Enzyme preparations were incubated at 37 ℃ for 2 hours with occasional shaking (every 30 min), and ALK was incubated overnight at 4 ℃. Afterwards, enzyme preparations were heated at 100 ℃ for 5 min to inactive enzymes, and pH was neutralized to 6.5-7 with HCl or NaOH, and the samples were freeze-dried. Samples are differentiated by the su x _B (for before) or _A (for after) (e.g. PECT1_A) for carbohydrase-or ALK-treatment prior to and after digestion, respectively.

Fermentation in The Swine in vitro Large Intestinal Model (SLIM)
The setup of SLIM was as follows: a fully computer-controlled in vitro model based on TIM-2 (Minekus et al. 1999) was used to mimic the swine large intestine (Long, de Vries, and Venema 2020). The pH (5.9) was controlled by continuous addition of 2 M sodium hydroxide. Standard ileal e ux medium of pigs (SIEMP) was used to simulate the materials entering the colon. The SIEMP, adapted from (Gibson, Cummings, and Macfarlane 1988)  and 0.4 g cysteine•HCl, plus 1 mL of the vitamin mixture. All medium components were purchased at Tritium Microbiology (Eindhoven, The Netherlands). The pig fecal inoculum was a standardized microbiota from growing pigs, freshly collected from the oor (48 pens with 6 pigs/pen, Hypor Libra x Hypor Maxter, Hendrix Genetics, Boxmeer, The Netherlands), but only material from the top was selected (so not toughing the oor). Feces was pooled and mixed with dialysate as described before (Long, de Vries, and Venema 2020).
In order to create a complete anaerobic environment, SLIM with 90 mL dialysate in each of the 4 individual units was ushed with gaseous nitrogen for at least 3 hours before incorporating the standardized microbiota. Thirty mL of the standardized microbiota was added to each SLIM-unit, making the total volume 120 mL. Figure 1 shows the experimental set-up for bre addition to SLIM. The microbiota was adapted to the model with SIEMP for 16 hours. During the adaptation phase, SIEMP was added into each SLIM-unit at a rate of 2.5 mL/h through the feeding syringe. At the end of the adaptation period, a 2-hour starvation period was performed, which was used to allow all the carbohydrates within SIEMP to be fermented. Afterwards, a shot of 5 grams of the different RSMs was given to the system at time point 0 hours, and incubation was continued for 24 h after that (Fig. 1).
Sequencing of the V3-V4 region of the 16S rRNA gene Microbial DNA extraction and sequencing of the V3-V4 region of the 16S rRNA gene were performed by BaseClear B.V. (Leiden, The Netherlands). Brie y, genomic DNA extraction was performed using the Quick-DNA™ Fecal/Soil Microbe Miniprep Kit (Zymo Research, California, USA) according to the manufacturer's instructions. Barcoded amplicons from the V3-V4 region of 16S rRNA genes were generated using a 2step PCR. 10-25 ng genomic DNA was used as template for the rst PCR with a total volume of 50 µL using the 341F (5'-CCTACGGGNGGCWGCAG-3') and the 785R (5'-GACTACHVGGGTATCTAATCC-3') primers (Klindworth et al. 2013) appended with Illumina adaptor sequences. PCR products were puri ed (QIAquick PCR Puri cation Kit, Venlo, The Netherlands) and the size of the PCR products were checked on a Fragment analyzer (Advanced Analytical, Ankeny, US) and quanti ed by uorometric analysis. Puri ed PCR products were used for the 2nd PCR in combination with sample-speci c barcoded primers (Nextera XT index kit, Illumina, city, CA, USA). Subsequently, PCR products were puri ed, checked on a Fragment analyzer and quanti ed, followed by multiplexing, clustering, and sequencing on an Illumina MiSeq with the paired-end (2x) 300 bp protocol and indexing. The sequencing run was analyzed with the Illumina CASAVA pipeline (v1.8.3) and demultiplexed based on sample-speci c barcodes.
The demultiplexed raw sequences obtained from BaseClear were processed using the QIIME2 pipeline (Bolyen et al. 2019). In short, reads were imported and quality ltered and dereplicated with q2-dada2 (Callahan et al. 2016). Next, dada2 was performed with paired-end reads and truncations parameters were as follows: the rst 17 and 14 base pairs were trimmed off in forward and reverse reads, respectively. And at position 280 base pairs the fragment was truncated in forward reads, and at position 230 base pairs for the reverse reads. The processed sequences were used for all the downstream analyses. Alpha-diversity (Shannon index) and β-diversity (weighted and unweighted UniFrac) were analyzed by the q2-phylogeny plugin (https://github.com/qiime2/q2-diversity).  (Pattathil et al. 2010). In brief, samples prepared above were applied (50 µL of 100 µg/mL in deionized water per well) to Costar 3598 96-well plates (Corning Life Sciences, city, country) and were dried to the well surfaces by evaporation overnight at 37℃. Control wells contained deionized water. The plates were blocked with 200 µl of 1% (w/v) bovine serum albumin (BSA) in Tris-buffered saline (50 mM Tris-HCl, pH 7.6, containing 100 mM sodium chloride) for 1 h. Blocking agent was removed by aspiration, and 50 µl of undiluted hybridoma supernatant were added to each well and incubated for 1 h at room temperature.
Supernatant was removed and wells were washed three times with 200 µl of 0.1% (w/v) BSA in Trisbuffered saline (wash buffer). Peroxidase-conjugated goat anti-mouse IgG, anti-mouse IgM, goat anti-rat IgG, or goat anti-rat IgM antibodies (Sigma-Aldrich), depending on the primary antibody used, were diluted 1:5,000 in wash buffer, and 50 µl were added to each well and incubated for 1 h. Wells were then washed ve times with 200 µl of wash buffer. Next, 3,3',5,5'-tetramethylbenzidine (TMB) solution (Sigma-Aldrich, St. Louis, USA) was freshly prepared according to the manufacturer's instructions, and 50 µL were added to each well. After 20 min, the reaction was stopped by adding 50 µl of 0.5 N sulfuric acid to each well. The OD of each well was read at a wavelength of 450 nm using a Multi-mode microplate reader (BioTek Synergy HTX, Abcoude, The Netherlands).
Polysaccharide panel screening. Polysaccharide panel screening of mAbs was carried out by ELISA against all lumen samples immobilized to 96-well plates. Duplicate preparations of each polysaccharide were used for all experiments reported here. Binding data was visualized in R via ComplexHeatmap package (Gu, Eils, and Schlesner 2016).

Statistics
Kruskal-Wallis Rank Sum Test (one-way ANOVA on ranks) was applied to compare α-diversities (Shannon index) among different RSM treatments and time points, and Wilcoxon Rank Sum Test was used for pairwise comparison in R version 3.5.3 (https://www.r-project.org/). Bonferroni was used to correct Pvalues. Permutational multivariate analysis of variance (PERMANOVA; REFERENCE) was performed to test the signi cance of β-diversity distances in QIIME2 (weighted and unweighted UniFrac) between nonprocessed and processed RSM. The results were visualized in R (R version 3.5.3).
The amplicon sequence variant (ASV) table (feature table of QIIME2) was normalized and ltered in R, and statistical analysis and visualizations were performed in R and STAMP (Parks, Tyson et al. 2014). The table was normalized via division by the sum of sequences in a given sample and multiplied by the minimum sum across all samples. Relative abundances were ltered as follows: values below a relative abundance threshold of 0.01% were not taken into account; taxa with a median relative abundance < 1% in all groups were not considered for statistical analysis. White's non-parametric t-test was applied to compare between CON and treatments. P-values were corrected using the Benjamini-Hochberg method. A q-value (corrected P-value) < 0.05 was considered signi cant.
Pearson correlations between continuous meta-variables and taxonomic variables were calculated and visualized in R. Parameters were set as follows: Missing values for meta-variables were handled as NO imputation (replacing missing data with substituted); zeros were kept for the calculation of correlation; a minimum number of 0.1% was considered for calculation; a minimum of 4 paired observations were required for calculation of correlations. P-values were corrected using the Benjamini-Hochberg method. A q-value (corrected P-value) < 0.05 was considered signi cant.
SCFA production between CON and the treated RSM substrates were compared and visualized in R.

Characteristics of non-processed and processed RSM
A comprehensive set of 155 plant cell wall glycan-directed monoclonal antibodies was used to screen untreated RSM by a ELISA-based assay (Pattathil et al. 2010;Pattathil et al. 2012), and 34 antibodies reacted with RSM (data not shown). These were subsequently used in the current study to obtain information on the presence and relative abundance of speci c epitopes that are characteristic of the different types of polymers in untreated RSM and RSM processed by ALK, CELL, PECT1, and PECT2. Figure 2 shows that both increases and decreases in epitope recognition occurred in ALK, CELL, PECT1, and PECT2 compared to CON. Samples from after (_A) and before (_B) predigestion clustered together according to each treatment, which indicated _A and _B from the same treatment had similar epitope accessibility. ALK strongly increased binding of non-fucosylated XG mAbs, while CELL, PECT1, and PECT2 led to disappearance of those compared to CON, regardless of _A and _B treatment. All the treatments increased the binding of "Linseed Mucilage RG-I group" directed mAbs, but had little effect on Xylan-2 and RG-Ic group compared to CON. Binding of MAC204 (AG-1), which is binding to gum tragacanth and to lettuce and green tomato RG-I preparations (arabinogalactan), disappeared with ALK_A, ALK_B, and CELL_B, while increased binding of CCRC-M107 (AG-2), which binds to linear and branched arabinans and RG-I preparations from diverse plants but does not bind to larch arabinogalactan (Pattathil et al. 2010), was observed in ALK_A and ALK_B. CELL_A, CELL_B, PECT1_A, PECT1_B, and PECT2_B led to disappearance of the binding of mAbs of "pectic backbone group" and CCRC-M 133, which also binds to linear and branched arabinans and RG-I preparations from diverse plants but do not bind to larch arabinogalactan. ALK, PECT1, and PECT2 increased binding of mAbs directed against the arabinogalactan side chains of RG-I (RG-I/AG). ALK, CELL, PECT1, and PECT2 signi cantly changed microbiota composition compared to CON To determine the changes in composition of the gut microbiota fed with a shot of 5 g CON, ALK, CELL, PECT1, or PECT2, a comparison of microbiota based on sequencing the V3-V4 region of the 16S rRNA gene was performed. Shannon indeces signi cantly decreased at t4, t6, t8 and t24, compared to that of t0 ( Figure S1A). When data from all of the time points were pooled, there were no signi cant differences among CON, ALK, CELL, PECT1, and PECT2 in Shannon index ( Figure S1B). Phylogeny based UniFrac methodology was then used to compare the β-diversity of the microbial communities between microbiota fed with non-processed and processed RSM. Unweighted UniFrac analysis (Fig. 3) shows that samples from processed RSM (ALK, CELL, PECT1, and PECT2) signi cantly (P = 0.004) separated from nonprocessed RSM (CON), and samples from different processed method clustered together. Samples from CON, ALK, CELL, PECT1, and PECT2 all clustered together (P = 0.125) with respect to weighted UniFrac ( Figure S2).
There were no signi cant differences between microbiota fed with RSM predigested before or after carbohydrase or ALK treatment with respect to both α-diversity ( Figure S3), or β-diversity (data not shown), which indicated predigesting before or after processing RSM had little effect on microbiota composition. Glycome pro ling of RSM (Fig. 2) also shows that samples from after and before processing clustered together according to each treatment, which indicated their polysaccharide compositions were similar to each other. Therefore, they were treated as duplicates in the following (microbial relative abundance) analyses.
Relative abundances of taxa within the pig microbiotas fed with non-processed and processed RSM were compared to identify signi cantly different bacterial taxa. At genus level (Fig. 4), seven genera were signi cant higher in relative abundance after ALK, CELL, PECT1, and PECT2 treatment compared to CON.
PICRUSt2 analyses revealed that microbial functional abundances related to carbohydrate metabolism and SCFA production were signi cantly increased with processed RSM compared to CON PICRUSt2 was performed to the 16S rRNA gene data to predict metagenomic functional pro les. In this study we focused on carbohydrate metabolism related microbial functions (Fig. 5) Six microbial pathways related to fermentation were signi cant higher in relative abundance when microbiotas were fed with ALK, CELL, PECT1, and PECT2 compared to microbiota fed with CON. These pathways were pyruvate dehydrogenase E1 component, short-chain fatty acids transporter, mannose-1phosphate guanylytransferase, superpathway of glucose and xylose degradation, sucrose degradation IV, and L-lysine fermentation to acetate and butanoate. The relative abundance of lactose/L-arabinose transport system permease protein was signi cant higher in ALK, PECT1 and PECT2 compared to CON. Figure 6 shows that the cumulative acetic, propionic, and butyric acid and total SCFA production were higher in ALK, CELL, PECT1, and PECT2 compared to CON. For acetic acid, more than 2 times greater production was observed when the microbiota was fed with ALK, CELL, PECT1, and PECT2 compared to when the microbiota was fed with CON. The production of propionic, butyric acid and total SCFA in ALK, CELL, PECT1, and PECT2 were more than 1.6 times higher than that in CON, except for propionic (1.3 times), and butyric acid (1.4 times) production in ALK.
Glycome pro ling shows that binding of mAbs in lumen digests were dynamic during the in vitro fermentation in SLIM To investigate the dynamic changes of polysaccharides structure in CON, ALK, CELL, PECT1, and PECT2 during in vitro fermentation, a time series of sampling was performed and the set of 34 mAbs was used to screen the lumen digests. Figure 7 shows that no or few binding signals were observed in Nonfucosylated XG, AG-2, and Pectic Backbone mAbs upon feeding CON, ALK, CELL, PECT1, and PECT2 at all time points.
No binding of CCRC-M83 that speci cally bind to Linseed Mucilage RG-I was observed at t0 in all treatments (just prior to addition of the bre shots), and binding signals appeared afterwards. For ALK, binding of CCRC-M83 decreased from t0.5 to t2, increased again at t4, and decreased to the lowest value at t24. For CELL, binding of CCRC-M83 decreased from t0.5 to t2, stabilized from t4 to t8, and decreased to the lowest value at t24. For CON, binding increased from t0.5 to t1, decreased from t2 to t6, increased at t8, and decreased to the lowest value at t24. For PECT1, binding increased from t0.5 to t4, and then decreased progressively until t24. For PECT2, binding increased slightly from t0.5 to t8, and then decreased to the lowest value at t24.
In terms of Xylan-2 recognizing mAb (i.e. CCRC-M105), weak bindings were detected at all time points after t0 with CELL and CON. Increased binding of CCRC-M105 was observed from t0.5 to t1 in ALK and the binding strength decreased from t2 to t4, slightly increased again at t6, and then decreased until t24.
For PECT1, increased binding of CCRC-M105 was observed from t0.5 to t6, which decreased afterward until t24. Binding for PECT2 was dynamic from t0.5 to t24, but the lowest binding was observed at t24.
As for the RG-Ic recognizing mAb (i.e. CCRC-M30), binding of CCRC-M30 was lower in ALK, PECT1, and PECT2 according to each time point compared to CELL and CON, but binding over time was dynamic. With respect to the AG-4 mAbs (recognizing arabinogalactans), weak and dynamic binding was observed at each time point in all treatments.
With respect to RG-I/AG mAbs, more active mAbs were observed compared to other groups of mAbs in all treatments. Within RG-I/AG mAbs, binding strengths of CCRC-M25 and CCRC-M60 were stronger than other RG-I/AG mAbs in all treatments, and their binding strengths were uctuating during the whole fermentation period and still existed at t24 in all treatments.
The correlation between binding of mAbs and relative abundance of microbial genera was also analyzed. Within the mAbs that recognizing non-fucocylated XG, CCRC-M93 had signi cant negative correlation with Bacteroidales S24-7 group.

Discussion
Cell wall polysaccharides composition of differently pretreated RSM Cell wall polysaccharides of RSM (Brassica napus) were comprehensively studied by chemical methods , which showed that they generally consisted of arabinan, homogalacturonan, Rhamnogalacturoan I (RG-I), type II arabinogalactan (AG), and xyloglucan (XG) (Siddiqui and Wood 1977). In a previous study (Long, de Vries and Venema submitted), we determined the monosaccharide constituent composition, which was in line with Pustjens et al. (2013). Our current ndings with the mAbs are also in line with this, since nonfucosylated XG-, RG-I-, pectic backbone-, and AG-recognizing mAbs bound to CON. By similar reasoning, xylan was also detected in CON, which was not reported before, although it is not entirely clear how much cross-reactivity the mAbs show.
ALK treatment intensively increased binding of mAbs that speci cally bind to nonfucosylated XG (Fig. 2). This is consistent with previous nding that alkali could extract XG from hemicellulose ), which led to XG being detected by the mAbs in the current study. Previous research also demonstrated that XG is linked via hydrogen bridges to the surface of cellulose micro brils, which is extractable with alkali, but not accessible by enzyme (Pauly et al. 1999). However, enzymatic treatments (PECT1, PECT2, and CELL) led to a disappearance of binding nonfucosylated XG-speci c mAbs compared to CON, which was unexpected. It could be that XG might become entrapped by other cell wall structures after pectinase or cellulase treatment. It is unknown whether this kind of the hidden phenomenon would reduce the degradability of XG.
Enzymatic treatments increased binding of linseed mucilage RG-I-speci c mAbs and some mAbs directed against the arabinogalactan side chains of RG-I (RG-I/AG), and reduced binding of some mAbs directed to pectic backbone, which indicated that PECT1 and PECT2 broke down pectic backbone and exposed RG-I and its arabinogalactan side chains. These results are supported by a previous study showing that PECT1 and PECT2 degraded pectic backbone and released RG-I and its arabinogalactan side chains ). Moreover, RG-I attached to cellulose micro brils can be released by cellulase (Oechslin, Lutz, and Amadò 2003), while at the same time it might block the accessibility of pectic backbone because of the shifting of the polysaccharides structure.
ALK treatment led to the disappearance of binding of mAbs that speci cally bind to AG-1 while it increased binding of mAbs that speci cally binding to AG-2, which indicated that ALK could cause arabinogalactan to be physically entrapped in other cell wall structures, while linear and branched arabinan become accessible (Pattathil et al. 2010).
They were lower in binding strengths of mAbs in AG-4 from ALK_B, PECT1_B, PECT2_B, and CELL_B, compared to ALK_A, PECT1_A, PECT2_A, and CELL_A, respectively. This indicated that small fragments were produced in AG-4 group after enzymatic and chemical treating RSM, which were dialyzed out in ALK_B, PECT1_B, PECT2_B, and CELL_B during the subsequent predigestion treatment. Overall, both enzymatic and chemical processed RSM can release some polysaccharides but physically entrap or shield some others at the same time.
Community structure of swine microbiota fed with differently pretreated RSM The PCoA analysis showed signi cant difference of both enzymatic and chemical treatments on swine microbial community (P = 0.004) compared to CON, based on unweighted UniFrac metric (Fig. 3); and overall there were no signi cant differences in this metric between treatments (Fig. 3). Dietary bre is known to have a considerable effect on gut microbiota composition (Yatsunenko et  These observations might explain the speci c genus changes in microbiota composition after microbiotas fed with ALK, PECT1, PECT2, and CELL compare to CON (Fig. 4), despite the differential effect of the different treatments as assessed by the glycome pro ling.
In order to understand which microbes were signi cantly in uenced by dietary supplementation of processed RSM, microbial genera were compared among the ve treatments (Fig. 4). The relative abundances of Ruminococcaceae NK4A214 group, Ruminococcaceae UCG-002, Ruminococcaceae UCG-005, Roseburia, Bi dobacterium, and Christensenellaceae R-7 group were signi cantly increased in ALK, PECT1, PECT2, and CELL (Fig. 4). Research has shown that genera of family Ruminococcaceae contain major (hemi)cellulolytic and pectinolytic species (Nyonyo, Shinkai, and Mitsumori 2014;Pettipher and Latham 1979). Thus, microbes from family Ruminococcaceae play an important role in degrading (hemi)cellulose and pectin in their activity against recalcitrant bre (Flint et al. 2008;Shinkai and Kobayashi 2007). Previous research demonstrated that the higher the relative abundance of Ruminococcaceae NK4A214 group, the higher the bre degradability (Hartinger, Gresner, and Südekum 2019; ). This observation was consistent with the current study, which showed that the SCFA production (which is usually used as evaluation of bre degradability in vitro) was also higher in ALK, CELL, PECT1, and PECT2 compared to CON. Furthermore, our study also showed that the relative abundance of Ruminococcaceae NK4A214 group, as well as Ruminococcaceae UCG-002 and Ruminococcaceae UCG-005, signi cantly correlated with acetic acid, propionic acid, butyric acid, and (not surprisingly) total SCFA production. Roseburia is a well-known butyrate-producing bacteria (Duncan et al. 2002), and a primary degrader of β-mannans (La Rosa et al. 2019). The signi cant increase in abundance of this genus was accordance with the increased butyric acid production in the current study ( Fig. 6), and increased predicted abundance of beta-mannosidases (EC.3.2.1.25) in the processed RSM groups (Fig. 5). Bi dobacterium is suggested to be a common cross-feeding bacteria for sugar utilization (Cockburn and Koropatkin 2016). Previous research has been shown that Bi dobacterium bi dum relies on the presence of a primary degrader in order to grow with either resistant starch or xylan both in vitro  and in vivo (Turroni et al. 2016). It has also been reported that numerous Bi dobacterium species grew to higher cell densities accompanied by upregulating their respective saccharolytic pathways when grown in co-culture compared to their growth in monoculture . This nding was consistent with the current study, where the relative abundance of Bi dobacterium was much higher than that of other genera, especially in ALK (data not shown). Alternatively, synergy existed among genera that were targeting the same substrates, possibly by specializing in degrading different motifs within the molecule.
The relative abundance of Christensenellaceae R-7 group increased after enzymatic and chemical treatment, which was consistent with a previous study where rumen microbiota was fed with bolytic enzyme-treated wheat straw (Ribeiro et al. 2020). Christensenellaceae plays an important role in degrading bre (Mao et al. 2015) and producing acetic and butyric acid (Morotomi, Nagai, and Watanabe 2012). All the observations above were supported by former reports that feed processed by carbohydrase enzymes could stimulate growth of speci c microbes (Wang et al. 2001;Giraldo et al. 2007;Giraldo et al. 2008). It can be speculated that the pretreatment of feed stuffs with carbohydrases causes the release of hydrolysis materials (presumably oligosaccharides), which promote the chemotactic response of speci c bacteria, and stimulates their attachment to feed particles, thereby leading to growth of these microbes

Fibre Degradation And Scfa Production
Enzymatic and chemical treatment on RSM increased the amount of acetic acid, propionic acid, butyric acid, and thereby total SCFA production (Fig. 6). This observation was consistent with the predicted functional pro les related to carbohydrate metabolism, where the relative abundance of bre breakdown and fermentation enzymes and pathways increased (Fig. 5). This nding was in accordance with previous studies that addition of pectolytic enzymes improved degradability of non-starch polysaccharides (NSP) of RSM in vitro ) and in broilers (Pustjens et al. 2014;De Vries et al. 2014). Giraldo et al. also reported that supplementation of exogenous cellulase increased SCFA production (Giraldo et al. 2007). Previous research demonstrated that supplementation of carbohydrase on feed before feeding could increase microbial protein production, ruminal cellulolytic bacterial numbers, and ruminal brolytic activity (Wang et al. 2001). Thus, the ndings above indicate that the enzymatic and chemical treatment on RSM could su ciently open cell wall architecture (refer to Fig. 2), to enable effective accessibility of NSP to bacterial degradation enzymes, and subsequently stimulate expression of microbial saccharolytic pathways.
No bindings of mAbs recognizing non-fucolysated XG were observed in lumen samples after the microbiota was fed with ALK or CON (Fig. 7), which was unexpected as binding signals were seen in the substrates themselves prior to addition to SLIM (Fig. 2). The hypothesis could be entertained that XG was immediately utilized by bacteria before our rst sampling time point (after 30 minutes), or that the XG structures were unable to be recognized by the mAbs after supplementing them to lumen, due to entrapment by other molecules. Binding signals still existed at t24 for all mAbs that showed signal at the t0.5 time point, which indicated that these structures cannot be degraded any further or more fermentation time was needed. Glycome pro ling showed that binding of mAbs recognizing each polysaccharide structure/epitope were dynamic during the 24 h fermentation period (Fig. 7). It is not unlikely that certain polysaccharide structure, such as (hemi)cellulose and pectin, were exposed to microbes stage by stage, due to ever increasing degradation of the cell wall structures over time. For instance, bacteria should break down side-chain AG before they can utilize RG-I. However, this hypothesis should be validated in future study.
The increased in abundance of Ruminococcaceae NK4A214 group, Ruminococcaceae UCG-002, and Ruminococcaceae UCG-005 in enzymatic and chemical treated RSM groups showed negative correlations with XG, RG-I, and RG-I/AG (arabinogalactan side chains of RG-I), and positively correlated with AG-1, AG-2, and pectic backbone. These observations suggested that ALK, CELL, PECT1, and PECT2 stimulated these genera to utilize XG, RG-I, and arabinogalactan side chains of RG-I, and exposed AG-1, AG-2, and pectic backbone to other bacteria in the current study. The abundance of Roseburia and Succinivibrionaceae UCG-001 in enzymatic and chemical treated RSM groups had negative correlations with XG and AG-1, which suggested that the treatments on RSM stimulated these genera to degrade XG and AG-1.

Conclusion
The present study clearly demonstrated that both enzymatic and chemical pretreatment on RSM shifted its cell wall polysaccharide structure, subsequently altering microbial community composition and functional pro le compared to untreated RSM, and eventually increased bre degradability as evaluated by SCFA production. Furthermore, glycome pro ling showed that the abundance of cell wall polysaccharides were dynamically changed during fermentation, and did not continuously decrease during the fermentation period as one might expect. Our ndings that ALK, CELL, PECT1 and PECT2 increased ber degradability in RSM could help guide feed additive strategies to improve e ciency and productivity in swine industry. The current study gave insight into how feed enzyme modulate microbial status, which provides good opportunity to develop novel carbohydrase, particularly in swine feed.     Principal coordinate analysis (PCoA) plot generated based on the calculated distances in unweighted UniFrac matrix. Samples were grouped by shape and color in terms of treatment and time point they belonged to, respectively: CON(sphere), ALK(plus), PECT1 (triangle), PECT2(diamond),CELL(square); a red-green-purple scale was used to indicate the fermentation time (red and purpledepict start and end of the fermentation period).

Figure 3
Principal coordinate analysis (PCoA) plot generated based on the calculated distances in unweighted UniFrac matrix. Samples were grouped by shape and color in terms of treatment and time point they belonged to, respectively: CON(sphere), ALK(plus), PECT1 (triangle), PECT2(diamond),CELL(square); a red-green-purple scale was used to indicate the fermentation time (red and purpledepict start and end of the fermentation period).

Figure 5
Signi cantly different metagenomic functions in relative abundance among ALK,CELL, CON, PECT1,and PECT2. * P < 0.05; ** P < 0.01. Signi cantly different metagenomic functions in relative abundance among ALK,CELL, CON, PECT1,and PECT2. * P < 0.05; ** P < 0.01. Cumulative Acetic (A), Propionic (B), and Butyric acid (C), and total short-chain fatty acid (D)production during fermentation of ALK, PECT1, PECT2, and CELL compared to CON. Fold change of short-chain fatty acid in ALK, PECT1, PECT2, and CELL compared to CON (E). Cumulative Acetic (A), Propionic (B), and Butyric acid (C), and total short-chain fatty acid (D)production during fermentation of ALK, PECT1, PECT2, and CELL compared to CON. Fold change of short-chain fatty acid in ALK, PECT1, PECT2, and CELL compared to CON (E).   Correlation between core bacterial genera,and SCFAand binding of mAbs. Statistical signi cance was determined for all pairwise comparisons using Pearson'smethod. * P < 0.05; ** P < 0.01. Size of the circles indicatescorrelation strength, which bigger sizes meaning higher correlations. Blue circles represent positive correlations (correlation coe cients from 0 to 1), whereas red circles represent negative correlations(correlation coe cients from 0 to -1).total, total SCFA production. The mAbs are grouped based on the cell wall glycans they predominantlyrecognize as shown in the panel on the left-hand side of the gure. XG, xyloglucan; RG-I, rhamnogalacturonan I; AG, arabinogalactan; RG-I/AG, arabinogalactanside chains of RG-I.

Figure 8
Correlation between core bacterial genera,and SCFAand binding of mAbs. Statistical signi cance was determined for all pairwise comparisons using Pearson'smethod. * P < 0.05; ** P < 0.01. Size of the circles indicatescorrelation strength, which bigger sizes meaning higher correlations. Blue circles represent positive correlations (correlation coe cients from 0 to 1), whereas red circles represent negative correlations(correlation coe cients from 0 to -1).total, total SCFA production. The mAbs are grouped based on the cell wall glycans they predominantlyrecognize as shown in the panel on the left-hand side of the gure. XG, xyloglucan; RG-I, rhamnogalacturonan I; AG, arabinogalactan; RG-I/AG, arabinogalactanside chains of RG-I.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download.