Zhongxia Li1,2†Huibin Zhang1†Li He1Yaqin Hou1Yingjuan Che1Tian Liu1Shaobai Xiong1Xuguang Zhang2Shunjing Luo1Chengmei Liu1*
Tingting Chen1*- 1State Key Laboratory of Food Science and Technology, School of Food Science and Technology, Nanchang University, Nanchang, Jiangxi, China
- 2BYHEALTH Institute of Nutrition and Health, Guangzhou, China
Introduction: Arabinoxylan (AX) is a versatile polysaccharide that shows various effects in modulating gut microbiota and health. The influence of arabinoxylan carbohydrate structural feature and feruloylation on fermentability and the effect of modulation of gut microbiota of AX was not clear.
Methods: Arabinoxylans from rice bran and corn bran (RAX and CAX), and their deferulyolated counterpart dRAX and dCAX were fermented using an in vitro fermentation model. Structural information was determined based on monosaccharide composition. Gas production of fermentation products, SCFAs production, pH change, and microbiota change were measured.
Results: RAX and dRAX posessed lower A/X ratio compared with CAX and dCAX. The gas and total SCFAs production were lower in RAX and dRAX, and the butyrate production were higher in RAX and dRAX compared with CAX and dCAX. Butyrate production was lower at dRAX compared to RAX. On the other hand, butyrate production was higher in dCAX than in CAX. The microbiota shift were different for the four fibers.
Discussion: The AXs from rice have a higher A/X ratio than the AXs from maize, suggesting more branching and a more complex side chain. The structural difference was crucial for the difference in fermentation pattern. Different Bacteroides species are responsible for the utilization of rice AXs and corn AXs. Although feruloylation had a minor effect on the overall fermentation pattern, it significantly affected butyrate production and alpha diversity. dRAX promoted less butyrate than RAX, which is associated with a significantly lower amount of Faecalibacterium prausnitzi. dCAX promoted more butyrate than CAX, which may be associated with a lower amount of Bacteroides ovatus and a higher amount of Blautia in dCAX compared to CAX. The effects of feruloylation on the fermentation pattern and the resulted microbiota shift of AX varied depending on the carbohydrate structure.
1. Introduction
Arabinoxylans (AXs) are an important source of dietary fiber found in various cereals such as wheat, corn, rice, sorghum and others, which possess various biological activities such as boosting immunity and antioxidant effects, strengthening the intestinal epithelial barrier, relieving constipation, and improving lipid and glucose metabolism (Comino et al., 2016; Kim et al., 2017; Wang et al., 2020). Arabinoxylans consist of a-1,4-xylan polymers that are predominantly substituted with arabinose side chains (Kaur et al., 2021) and cannot be degraded by mammalian enzymes (Paesani et al., 2020). The structures of AXs derived from different raw materials or extracted by different methods are not very consistent, each may have different unique structural features (Chen et al., 2019). The carbohydrate structure [different sugar composition, molecular weight (MW), type of branching, etc.] and feruloylation dominate the structures of AXs (Snelders et al., 2013), play a significant role in the ability of bacteria to degrade the polysaccharide (Wang et al., 2020).
The carbohydrate structure is the important factor of properties of AXs, which likely affects the fermentability of the gut fecal microbiota after AXs fermentation, including the fermentation rate of the gut microbiota, short-chain fatty acids (SCFAs) production, and the abundance and composition of the gut microbiota (Damen et al., 2011; Kaur et al., 2011; Rumpagaporn et al., 2015; Xu et al., 2022). Arabinoxylans with less branch or lower ratio of arabinose to xylose (A/X) are more easily fermented than highly branched or substituted AXs (Pollet et al., 2012; Mendis et al., 2016; Tiwari et al., 2019). Meanwhile, compared AXs with lower average degree of polymerization, the higher average degree of polymerization of AXs promoted the growth of probiotics such as Roseburia and Eubacterium rectale (Agathobacter rectalis), and increase in butyrate levels in the caeca of the rats (Damen et al., 2011).
Feruloylation is the binding of ferulic acid (FA) to arabinoxylan. Ferulic acid is the major bound polyphenol of AXs known for its health-promoting biological effects, such as antioxidant activity, anti-cancer activity, and cardiovascular disease, to name a few (Ghosh et al., 2017; Johnson et al., 2020). The bound FA content of AXs is related to the extraction method applied. AXs can be extracted usually using alkali and enzymatic treatments. Enzymatic treatment produced the hydrolytic degradation products of AXs (AXOS) with lower molecule weight and higher bound FA content than the AXs with alkali treated (Chen et al., 2019). The bound FA could delay the overall degradation rate of AXOS during fermentation, promote the production of SCFAs and the growth of beneficial bacteria (Snelders et al., 2014; Gong et al., 2019), while the effect of bound FA on the ability of AXs to affect the gut microbiota still unclear.
Here we assume that carbohydrate structure and feruloylation make a difference to the intestinal probiotic properties of AXs, and no relevant reports have been found. Therefore, we obtain two kinds of AXs (from rice bran and corn bran) with different carbohydrate structures but similar content of bound FA and then degrade the bound FA with feruloyl esterase to obtain two less feruloylated AXs. Investigate the carbohydrate structure and feruloylation of AXs on the intestinal probiotic properties of AX by in vitro.
2. Materials and methods
2.1. Materials
Defatted rice bran and corn bran powder were purchased from the Jiangxi Tianyu Oil Co., Ltd. (Jiangxi Province, China) and the native market, respectively. Human salivary α-amylase (A1031), pepsin (P7012), and pancreatin (P1750) were purchased from Sigma-Aldrich Chemical Co., Ltd. (Shanghai, China). Thermostable α-amylase (S10004), neutral protease (S10013), and 3 kDa dialysis membranes (MD3044) were purchased from Yuanye biological technology Co. (Shanghai, China). Feruloyl esterase (EC 3.1.1.73) from Clostridium thermocellum was purchased from Megazyme. All other chemicals and reagents were of analytical purity.
2.2. Extraction of arabinoxylan with different carbohydrate structures and similar content of bound ferulic acid
AXs extraction followed the procedure described previously with some modifications (Doner and Hicks, 1997). The bran was crushed to pass a 60-mesd sieve. Deionized water was added to the bran (10:1, v/w), the pH was adjusted to 7 with 2 mol/l NaOH, and then stirred for 1 h at 100°C. The starch and protein in the bran were successively removed with heat-stable α-amylase (2 g/100 g bran, hydrolysis for 1 h at 90–95°C) and neutral protease (1 g/100 g bran, hydrolysis for 6 h at 50°C and pH = 6), then the solution was boiled for 15 min to inactivate the enzymes. The reaction solution was centrifuged to obtain supernatant and residue. For the supernatant, add absolute ethanol and dry the precipitate to obtain the water-extractable arabinoxylan. For the residue, add 0.25 mol/l NaOH and stir for 1 h at 60°C, then slowly add 30% H2O2 and hold for 4 h. Centrifuge the mixture, wash the residue once with 1 mol/l NaOH solution, and centrifuge again. Adjust the pH of the aqueous solution to 7 with 6 mol/l HCl, then add 4 volumes of absolute ethanol and dry the precipitate to obtain alkali-extractable arabinoxylan. RAX is the water-extractable arabinoxylan from rice bran and CAX is the alkali-extractable arabinoxylan from corn bran.
2.3. Bound ferulic acid depletion of arabinoxylan
Depletion of bound FA from AXs was performed as previously described with slight modifications (Snelders et al., 2013). FA was enzymatically hydrolyzed from AXs (5 g/100 ml deionized water) with a feruloyl esterase (230 μl/5 g AX) treatment at pH = 7 and 50°C for 24 h, boiled for 15 min to inactivate the feruloyl esterase. Then remove the free FA with Amberlite XAD-4 column, the solution was brought over the column with peristaltic pump and the run was collected and loaded on the column for another two times. Then the column was washed twice with deionized water, the wash fraction and the fraction of the reloaded solution fraction as AX solution, then absolute ethanol was added, and the precipitate was dried to obtain the ferulic acid depleted AX. RAX and CAX depleted the bound FA to obtain dRAX and dCAX, respectively.
2.4. Bound ferulic acid content of arabinoxylan
The bound FA content of AXs was measured as previously described with slight changes (Snelders et al., 2013). For determination of total FA, samples were hydrolyzed by addition of 2 mol/l NaOH and incubated for 2 h to release the ferulate group from the conjugate under N2 atmosphere. Then 6 mol/l HCl was added to adjust the pH to 2.0 and acidify the released ferulate group. The free FA is then extracted from the aqueous solution with three volumes of ethyl acetate three times. The pooled ethyl acetate extracts are then evaporated to dryness and resuspended in a known volume of methanol. For the determination of free FA, no saponification step occurred, and samples were analyzed from the acidification step. The amount of total and free ferulic acid is then quantified using high performance liquid chromatography (Agilent 1260, United States). The system combined with a UV detector at 320 nm and a Select C18 column (Silgreen GHAAL12S05, Greenherbs, China) at 20°C. The mobile phase contained solvent A (1.0% glacial acetic acid-water solution) and solvent B (acetonitrile) at a flow rate of 0.8 ml/min. 80% solvent A and 20% solvent B were eluted isocratically.
2.5. Monosaccharide and uronic acid compositions of arabinoxylan
The monosaccharide and uronic acid composition of AXs was analyzed by high-performance anion-exchange chromatography (ICS-6000, Thermo Scientific, United States) with pulse amperometric detection (Xie et al., 2013). The AX samples (20 mg) were hydrolyzed with 0.5 ml of 12 mol/l H2SO4 in a sealed glass tube and stirred for 30 min with a magnetic stirrer in an ice bath. Then, the H2SO4 was diluted to 2 mol/l for another 4 h of hydrolysis at 105°C. After the solution cooled to room temperature, the concentration of the mixture was diluted to about 100 μg/ml. The monosaccharide standards (including rhamnose, fructose, arabinose, xylose, galactose, mannose, and glucose) and the uronic acid standards (galacturonic acid and glucuronic acid) were dissolved in ultrapure water. The solution was filtered through a 0.22 μm syringe filter before injection into the system. For analysis, a CarboPac™ PA20 guard column (3 mm × 30 mm) and a CarboPacTM PA20 analytical column (3 mm × 150 mm) were connected in series. After injection of 25 μl of the sample, a quaternary gradient system consisting of ultrapure water (eluent A), 10 mmol/l NaOH solution (eluent B), 250 mmol/l NaOH (eluent C), and 500 mmol/l CH3COONa (eluent D) at 30°C was used to separate monosaccharide and uronic acid. UV detection was performed at 280 nm and the flow rate was 0.5 ml/min.
2.6. In vitro digestive pretreatment
Samples were subjected to simulated upper gastrointestinal tract digestion using a previous method with slight modifications (Lebet et al., 1998). The sample (2 g dry weight) was resuspended in 14 ml phosphate buffer saline (20 mM. pH = 6.9). The a-amylase solution (224 μl. 28.4 U/ml) was added and the mixture was then continuously stirred on a magnetic stirrer at 37°C for 10–20 min. After the pH was adjusted to 2.0, pepsin (34 μl, 51 U/ml) was added and incubated at 37°C for 30 min with continuous stirring. The pH was adjusted to 6.9, then trypsin (334 μl, 5 mg/ml) was added and incubated at 37°C for 90 min with constant stirring. After digestion, the suspension was boiled for 5 min to inactivate the enzymes. The cooled suspension was dialyzed against purified water using 3 kDa dialysis membranes for 3 days, changing the water every 6 h, and then freeze-dried for in vitro fermentation.
2.7. In vitro fecal fermentation
The in vitro fermentation was conducted under strictly anaerobic conditions and performed as per the method (Long et al., 2015). Fresh stool was collected from three healthy donors between the ages of 22 and 28 years who had no history of intestinal disease and had not been treated with antibiotics for at least 3 months before stool collection. This study was approved by the Ethics Committee of Nanchang University. All experiments were performed under guidelines and regulations of the IRB protocol. Informed consent was obtained from all donors. Fresh stool samples were collected from participants and immediately transferred to an anaerobic workstation (Bactron300-2, Shellab Bactron Company, United States) containing 5% H2, 5% CO2, and 90% N2. Three volumes of PBS medium (PBS medium was sterilized and deoxygenated, 1 l of PBS medium contains 8 g NaCl, 0.2 g KCL, 1.15 g Na2HPO4, 0.2 g KH2PO4, 0.5 g L-cysteine, 1 ml 1 mg/ml resazurin, and the rest is distilled water) were added to the samples and shaken until dispersion. The slurries were filtered through four layers of gauze cloth to remove visible particles. The fecal inoculum solution (1 ml) and 4 ml of PBS medium (containing 50 mg each of RAX, dRAX, CAX, and dCAX) were mixed. FOS was used as a positive control. The mixtures were incubated at 37°C. After 0, 6, 12, 24, and 48 h, gas production was measured with a syringe, and culture samples were collected and centrifuged at 13000 rpm for 5 min. The supernatants were stored at –80°C until analysis of SCFAs. The precipitates obtained after 24 h from all samples and at the beginning (0 h of Blank) were used for DNA extraction.
2.8. Short-chain fatty acids analysis
SCFAs (containing acetate, propionate, and butyrate) were determined as previously described with slight modifications (Zhao et al., 2006). 400 μl of the supernatant was added to a 2 ml centrifuge tube. 1 mol/l HCl (80 μl) and 10 mmol/l 2-ethylbutanoic acid (65 μl) were added to the supernatant and vortex shaken. Then, filtered through a 0.22 μm syringe filter, 1 μl was injected for gas chromatographic measurements (7080B, Agilent, United States) using a capillary column (HP-FFAP, 30 m × 0.25 mm × 0.25 μm, United States). The flow rates of N2, H2, and air were 1.5, 40, and 400 ml/min, respectively. The program initial temperature was 60°C, rose to 160°C at 10°C/min, maintained at 160°C for 1 min, increased to 200°C at 10°C/min, and held at 220°C for 5 min. The injector and detector temperatures were set at 200 and 240°C, respectively.
2.9. 16S rRNA sequencing and bioinformatics
Precipitates from the ferments were used for DNA extraction as previously described (Yu et al., 2021). DNA was extracted using the OMEGA Soil DNA Kit (D15625163 01; Omega Bio-Tek, Norcross, GA, United States). The concentration of extracted DNA was determined using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States), and purity was determined by agarose gel electrophoresis. Forward primer 338F (5′-ACTCCTACGGGAGGCAGCA-3′) and reverse primer 806R (5′-GGACTACHVGGGTWICTAAT-3′) were amplified at the V3–V4 regions of the bacterial 16S sRNA gene by PCR. Sample-specific 7-bp barcodes were inserted into the primers for multiplex sequencing. PCR amplicons were purified using VAHTS DNA Clean Beads (Vazyme, Nanjing, China) and quantified using the Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, United States). After individual quantification, all PCR amplicons were combined in equal amounts, and pair-end 2 × 250 bp sequencing was performed using the Illumina NovaSeq6000 system according to standard protocols at Shanghai Personal Biotechnology (Shanghai, China).
The obtained fastq sequences were further processed using the DADA2 pipeline with QIIME2 2019.4 (Callahan et al., 2016). Briefly, raw sequence data were demultiplexed using the demux plugin and primers were cut using the cutadapt plugin (Martin, 2011). Sequences were then quality filtered, denoised, merged, and chimeras removed using the DADA2 plugin (Callahan et al., 2016). Non-singleton amplicon sequence variants (ASVs) were aligned using mafft (Katoh et al., 2002) and used to generate a phylogeny using fasttree2 (Price et al., 2009). Alpha diversity metrics, including observed species and Shannon index, and beta diversity metrics were estimated using the diversity plugin, diluting samples to 35,582 sequences per sample. Taxonomy was assigned to ASVs using the classification Klearn naive Bayes taxonomy classifier in the Feature Classifier plugin (Bokulich et al., 2018) against the Greengenes database (v.13.8; DeSantis et al., 2006), which was set as the reference for taxonomic classification of each ASV. The original papers presented in this study are in the public domain. These data can be found at: https://doi.org/10.6084/m9.figshare.21695396.v1.
2.10. Statistical analysis
All experiments were performed in triplicate. Statistical differences between groups in the parameters of physicochemical properties, gas production, SCFAs production and α-diversity were analyzed using one-way analysis of variance (ANOVA) and Tukey’s post hoc test. The LEfSe method was performed at the ASVs level (using one-against-all comparisons, less strict) with LDA score large than 3 and assessed by Wilcoxon text. Data were analyzed using SPSS Statistics software (version 26.0, Ine. Chicago. IL, United States) and R 4.1.3. A significance level of p < 0.05 was set for all tests.
3. Results and discussion
3.1. Structural properties of the arabinoxylan
Two arabinoxylans were obtained from rice bran (RAX) and corn bran (CAX). The monomeric units of RAX and CAX were both predominantly xylose and arabinose (Table 1). The A/X ratio of RAX was 1.11, which was greater than that of CAX (A/X ratio 0.82). A higher A/X ratio was also found in arabinoxylans from rice bran compared to arabinoxylans from corn bran, suggesting that the structure of RAX is more branched and has complex side chains (Rose et al., 2010). RAX and CAX were then incubated with ferulic acid esterase to obtain the deferuloylated counterparts. The content of bound FA was reduced from 0.98 and 1.05 mg/g in RAX and CAX to 0.17 and 0.15 mg/g in dRAX and dCAX, respectively (Table 1).
Table 1. Bound ferulic acid content (mg/g AX) and monosaccharide compositions (%) of RAX, CAX, dRAX, and dCAX.
3.2. Effects of arabinoxylan on the gas production
Gas production is an important indicator of the fiber fermentation profile. The gas accumulation curve showed a two-phase fermentation pattern for all fibers (Figure 1). The slope of the first phase was largest for FOS and smallest for CAX. A flatter slope indicates slower gas accumulation and possibly slower fermentation. AXs were fermented more slowly than FOS. CAX and dCAX were fermented more slowly than RAX and dRAX. Removal of feruloylation did not affect the first phase of fermentation. RAX and dRAX both fermented rapidly within the first 6 h, then slowly until 48 h with little gas formation. After debranching in the first 6 h, the remaining AXs were poorly fermented and gas accumulation reached a plateau (Gu et al., 2021). The total gas accumulation of RAX and dRAX was lower than that of CAX and dCAX, which is related to the structural difference. Rice AXs which has more branches and complex side chains are more difficult to ferment than Corn AXs (Rumpagaporn et al., 2015). The decrease in gas accumulation after 12 h and 24 h in dRAX suggests that removal of feruloylation from RAX altered the fermentation pattern. Feruloylation was found to inhibit fermentation of AXOS. Specifically, removal of FA residues from feruloylated xylooligosaccharides makes them more amenable to further hydrolysis by xylanase and other degradative enzymes (Long et al., 2018). However, the effect of feruloylation on long chain AXs is still not clear. There is no difference between CAX and dCAX in gas accumulation suggesting that the inhibitory effect of FA depends on the arabinoxylan structure.
Figure 1. Gas production of RAX, CAX, dRAX, and dCAX compared to the Blank and FOS (positive control) by in vitro human fecal fermentation. FOS was used as a positive control. Blank did not contain any substrate. RAX and dRAX stand for high and low bound ferulic acid rice bran arabinoxylan extracted with water, respectively. CAX and dCAX stand for high and low bound ferulic acid corn bran arabinoxylan extracted with 0.25 mol/l NaOH, respectively.
3.3. Effects of arabinoxylan on the SCFAs production
To better understand the effects of carbohydrate structure and feruloylation on the fermentation characteristics of AXs, the production of SCFAs during the 48 h fermentation period was evaluated. The overall production of SCFAs largely followed the trend of gas fermentation (Figure 2A). The major SCFAs in all AXs were acetate and propionate, followed by butyrate. Over the course of 48 h, acetate production showed a similar trend to the overall SCFA production (Figure 2B). Propionate production of AXs from rice bran increased significantly within 24 h and then reached a plateau, whereas propionate production of AXs from corn bran increased steadily within 48 h (Figure 2C). This is consistent with previous reports that AX induces an increase in propionate output (Chen et al., 2019; Nguyen et al., 2020). In the first 6 h, butyrate production exhibited a pattern comparable to that of total SCFA production. AXs from rice bran produced more butyrate than AXs from corn bran at all time points (Figure 2D). On the other hand, butyrate production from RAX was higher than that from dRAX and butyrate production from CAX was lower than that from dCAX. FA in the free from has been reported to accumulate butyrate by promoting the growth of butyrate-producing bacteria (Song et al., 2020). Similar to gas accumulation, the effect of bound FA on butyrate production was also not clear. Bound FA was found to inhibit fermentation on AXOS and did not affect butyrate concentration (Snelders et al., 2014). CAX with FA bound at 0.65 mg/g and 1.15 mg/g were similarly fermented and produced similar amounts of butyrate (Zhang et al., 2019). Here, the different effects of FA on butyrate were possibly due to the different structure of arabinoxylans. In summary, structural differences in carbohydrate structure were related to differences in gas and total SCFA production, and feruloylation was related to butyrate production.
Figure 2. (A) Total SCFAs, (B) acetate, (C) propionate, and (D) butyrate production of each group by in vitro human fecal fermentation. Abbreviations are explained in Figure 1.
3.4. Effects of arabinoxylan on the overall structure of gut microbiota
Summary of taxa at the phylum level showed that the relative abundance of Bacteroidetes was increased during fermentation (Figure 3A). All AXs significantly increased the relative abundance of Bacteroidetes, but not the relative abundance of Proteobacteria and Fusobacterium. Another reason for the increase in relative abundance of Bacteroidetes may be the result of a decrease in Firmicutes. The genus Faecalibacterium, which belongs to the phylum Firmicutes, was predominant in the initial fecal sample and decreased significantly after 24 h of fermentation (Figure 3B). Compared with Blank and FOS, AXs promoted Bacteroides. The genera Bacteroides have been shown to possess the genetic and functional traits necessary for arabinoxylan utilization (Nguyen et al., 2020). The xylanolytic bacteria possessed polysaccharide utilization loci (PUL) that encode for the xylan utilization system (Zhang et al., 2014). Prevotella was another xylanolytic bacterium that was not presented due to interindividual differences (Fehlner-Peach et al., 2019). Unlike the selfish Prevotella, Bacteroides supported a complex metabolic web of cross-feeding interactions. Partial breakdown products were released by the action of xylan utilization system and shared with neighboring bacteria, which promoted the increase of alpha diversity (Grondin et al., 2017).
Figure 3. The relative abundance of the microbiota at the (A) phylum levels and the (B) genus levels of each group at 24 h in vitro fermentation. Abbreviations are explained in Figure 1.
The alpha diversity of the microbiota after 24 h fermentation was indicated by Shannon diversity, which is an indicator of both richness and evenness, and observed features, which is a direct indicator of ASV richness (Figure 4). After 24 h of fermentation, the Shannon diversity was significantly increased in the blank group (Figure 4A). This is due to the increase in observed features. It is unlikely that bacterial species increase in a closed ecosystem such as in vitro fermentation. The observed features were indicative of the coverage of the sampling depth (Figure 4B). The low abundance species that were initially not detected grew and became detectable in the Blank group. In comparison, these low abundance species remained present and undetectable in the fiber treatment groups. All AXs slightly increased the Shannon index and observed features compared to Initial. Long-chain AXs have been shown to restore the alpha diversity, which is decreased in high-fat diet fed mice (Neyrinck et al., 2011). Shannon index and observed features were lower in RAX than in CAX, possibly related to the difficulty of fermentation of RAX. Higher alpha diversity was observed in dRAX and dCAX than in RAX and CAX, respectively. Removal of FA is thought to slightly increase alpha diversity during fermentation. The in vivo study also showed that the binding of FA increased the Shannon diversity of the gut microbiota of XOS (Ou et al., 2016).
Figure 4. Alpha-diversity of the bacterial community as measured though (A) Shannon index and (B) observed features at 24 h in vitro fermentation. Abbreviations are explained in Figure 1.
The PCA plot showed that the composition of the gut microbiota had shifted significantly by FOS and AXs after 24 h of fermentation (Figure 5A). On the PC1 axis, separation was observed between AXs from rice bran and corn bran, and the PC1 score was significantly different, explaining the different effect of AXs from two origins on the microbiota (Figure 5B). Bacteroides ASV_2294 has been associated with the degradation of maize AXs and Bacteroides ASV_10074 with the degradation of rice AXs. Although both belong to Bacteroides, these two Bacteroides may differ greatly in PULs. B. thetaiotaomicron and B. ovatus are the two most studied Bacteroides that have different glycan niches. Both Bacteroides harbor about 100 PULs, but almost few in common (Grondin et al., 2017). In addition, the removal of feruloylation also shifted the microbiota on PC1, suggesting feruloylation interference with the degradation of the AXs by Bacteroides.
Figure 5. (A) The principal component analysis (PCA) of community structures and (B) the scores of the primary principle component (PC1) at 24 h in vitro fermentation with each group. Abbreviations are explained in Figure 1.
3.5. Effects of arabinoxylan on the specific bacterial taxa as biomarkers in the microbiota
The bacteria responsible for the difference between groups were selected by LEfSe analysis with a linear discriminant analysis score greater than 3 (Figure 6). Comparing RAX and CAX, unidentified Bacteroides, Faecalibacterium prausnitzi, unidentified Coprococcus, Roseburia inulinivorans, Bacteroides uniformis, unidentified Dialister, Clostridium clostridioforme, and Mitsuokella multacida were more abundant in RAX, whereas Bacteroides ovatus, unidentified Blautia, and Bacteroides acidifaciens were more abundant in CAX. In agreement with PCA analysis, different species of Bacteroides were responsible for the utilization of RAX and CAX, respectively. In addition, butyrogenic bacteria, including Faecalibacterium prausnitzi and Roseburia inulinivorans, which are more abundant in RAX, explained the higher butyrate production of RAX than CAX (Figure 2D). Roseburia and Blautia are the main genera belong to the Lachnospiracea family. Genomic analysis of Lachnospiraceae revealed a considerable capacity to utilize diet derived polysaccharides including starch, inulin and AX. Roseburia species are able to degrade AXs through the activity of xylanase, α-L-arabinofuranosidase and β-xylosidase (Vacca et al., 2020).
Figure 6. Most differentially enriched taxa between groups after 24 hours of in vitro fermentation, identified by LEfSe analysis (A) RAX and CAX, (B) RAX and dRAX, (C) CAX and dCAX. Each color represents a treatment group. The marker points in the taxonomic cladogram represent the difference at the different taxonomic levels between the compared groups. The corresponding taxa were labeled as taxa level followed by the names. The different species with a linear discriminant analysis (LDA) of more than 3 were listed with their respective LDA scores.
Arabinoxylans are primarily degraded by P. copri, B. longum, and certain Bacteroides species. In addition to P. copri, which behaves “selfishly” and recycles AXs itself, B. longum and certain Bacteroides species released AXOS and xylan through the activity of endoxylanase, β-xylosidase, and α-arabinofuranosidase (Despres et al., 2016; Nguyen et al., 2020). It was shown that the complexity of AXs molecules strongly influenced the fermentation rate (Chen et al., 2017). Bifidobacterium spp. are known to degrade xylan, but which were not identified in this study, possibly due to interindividual differences. Bifidobacterium longum is a primary degrader of AXs, capable of cleaving the complex AX structure by AX-degrading glycosidase including β-xylosidase and α-arabinofuranosidase (Komeno et al., 2019). Other Bifidobacterium species did not possess the AX-degrading genes and were therefore unable to degrade long-chain AX themselves, but can degrade AXOS (Riviere et al., 2014).
Compared to dRAX, Faecalibacterium prausnitzi was promoted in RAX, which explained the higher butyrate production in RAX than in dRAX. It is suggested that the bound FA could promote the butyrate-producing bacterium and increase the butyrate production after fermentation of AX. However, comparing dCAX and CAX, the differences were related to saccharolytic species B. ovatus, unidentified Bacteroides and unidentified Blautia. Bacteroides ovatus possesses feruloyl esterase, which leads to a different degradation pattern of the xylanolytic reaction (Kmezik et al., 2020). With the decrease of propionated producing B. ovatus, Blautia, a butyrate-producing bacterium was increased in dCAX, which explained the increase of butyrate production in dCAX compared to CAX. Other ferulic acid esterase active bacteria such as lactic acid generating bacteria (Limosilactobacillus, Lentilactobacillus, Bifidobacterium, etc.) that were identified in this study possibly due to the microbiota composition of the fecal inoculum. It is suggested that feruloylation may interference with the xylanolytic pattern. And the impact of feruloylation differed depending on the AX structures.
Bound FA reduced the fermentability of AXOS because the feruloylated and diferuloylated arabinose substituents limited the access of xylanolytic enzymes to their target sites (Snelders et al., 2014). Compared to AXOS, AX is more complex and requires a larger repertoire of hydrolases to access and utilize the complex structure. AX usually embedded in the complex cell wall matrix with the presence of monomers and dimers of FA. Xylanase is required to dissolve part of the cell wall structure by forming short feruloylated xylooligosaccharides, and the presence of feruloylation hinders xylanase activity. Feruloyl esterases have been used to release FA and diFA from the complex cell wall and led to further degradation of AX (Oliveira et al., 2019). It was believed that the removal of FA residues from feruloylated xylooligosaccharides makes them again more accessible for further hydrolysis by xylanase and other degrading enzymes (Long et al., 2018). Similarly, when the long-chain AXs were cross-linked by FA, the xylan backbone was more resistant to degradation and the arabinose side chains were utilized first. The gel-like structure sterically hindered enzyme activity (Hopkins et al., 2003; Zhang et al., 2019). Compared to AXOS and AX matrix, the presence of bound FA only showed influence on a couple of bacteria for extracted AX polymers. For AX polymers that are released from the cell wall matrix, steric hinder induced by crosslink no longer existed. Compared to feruloylated AXOS, the amount of bound FA in CAX and RAX are relatively low. This also partially explained the only 6 bacteria species are affected by the bound FA.
4. Conclusion
In conclusion, AXs with different structure and feruloylation fermented differently. Fermentability was mainly determined by the carbohydrate backbone and its branches, including pH, gas production, and SCFAs profile. Although feruloylation had a minor effect on the fermentation pattern, it showed a significant effect on specific bacteria. Compared to RAX, removal of feruloylation significantly reduced the amount of Faecalibacterium prausnitzi, a known beneficial bacterium that promotes butyrate production. Compared to CAX, removal of feruloylation significantly reduced the amount of Bacteroides ovatus, a known propionate-producing bacterium that degrades complex carbohydrates, and increased the amount of Blautia, a butyrate-producing bacterium. The effects of feruloylation on the fermentation pattern and the resulted microbiota shift of AX varied depending on the carbohydrate structure.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Nanchang University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
ZL: conceptualization and methodology. HZ: data curation, writing-original draft preparation, and investigation. LH: validation. YH: formal analysis. YC: visualization. TL: resources. SX: software. XZ: reviewing and editing. SL: supervision and investigation. CL: supervision and salidation. TC: conceptualization and writing-reviewing and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was financially supported by the National Natural Science Foundation of China (Nos. 31901703, 32160530, and 32072257), China Postdoctoral Science Foundation (2020M671974), the Natural Science Foundation of Jiangxi Province (No. 20202BAB215015), Central Government Guide Local Special Fund Project for Scientific and Technological Development of Jiangxi Province (20212ZDD02008), and Jiangxi Research Project Funding for Selected Postdoc (2020KY07).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Bokulich, N. A., Kaehler, B. D., Rideout, J. R., Dillon, M., Bolyen, E., Knight, R., et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2's q2-feature-classifier plugin. Microbiome 6:90. doi: 10.1186/s40168-018-0470-z
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chen, Z. Y., Li, S. S., Fu, Y. F., Li, C., Chen, D. W., and Chen, H. (2019). Arabinoxylan structural characteristics, interaction with gut microbiota and potential health functions. J. Funct. Foods 54, 536–551. doi: 10.1016/j.jff.2019.02.007
Chen, T., Long, W., Zhang, C., Liu, S., Zhao, L., and Hamaker, B. R. (2017). Fiber-utilizing capacity varies in Prevotella-versus Bacteroides-dominated gut microbiota. Sci. Rep. 7:2594. doi: 10.1038/s41598-017-02995-4
Comino, P., Collins, H., Lahnstein, J., and Gidley, M. J. (2016). Effects of diverse food processing conditions on the structure and solubility of wheat, barley and rye endosperm dietary fibre. J. Food Eng. 169, 228–237. doi: 10.1016/j.jfoodeng.2015.08.037
Damen, B., Verspreet, J., Pollet, A., Broekaert, W. F., Delcour, J. A., and Courtin, C. M. (2011). Prebiotic effects and intestinal fermentation of cereal arabinoxylans and arabinoxylan oligosaccharides in rats depend strongly on their structural properties and joint presence. Mol. Nutr. Food Res. 55, 1862–1874. doi: 10.1002/mnfr.201100377
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Despres, J., Forano, E., Lepercq, P., Comtet-Marre, S., Jubelin, G., Chambon, C., et al. (2016). Xylan degradation by the human gut Bacteroides xylanisolvens XB1A(T) involves two distinct gene clusters that are linked at the transcriptional level. BMC Genom. 17:326. doi: 10.1186/s12864-016-2680-8
Doner, L., and Hicks, K. (1997). Isolation of hemicellulose from corn fiber by alkaline hydrogen peroxide extraction. Cereal Chem. 74, 176–181. doi: 10.1094/CCHEM.1997.74.2.176
Fehlner-Peach, H., Magnabosco, C., Raghavan, V., Scher, J. U., Tett, A., Cox, L. M., et al. (2019). Distinct polysaccharide utilization profiles of human intestinal Prevotella copri isolates. Cell Host Microbe 26, 680–690.e5. doi: 10.1016/j.chom.2019.10.013
Ghosh, S., Basak, P., Dutta, S., Chowdhury, S., and Sil, P. C. (2017). New insights into the ameliorative effects of ferulic acid in pathophysiological conditions. Food Chem. Toxicol. 103, 41–55. doi: 10.1016/j.fct.2017.02.028
Gong, L. X., Wang, H. N., Wang, T. X., Liu, Y. L., Wang, J., and Sun, B. G. (2019). Feruloylated oligosaccharides modulate the gut microbiota in vitro via the combined actions of oligosaccharides and ferulic acid. J. Funct. Foods 60:103453. doi: 10.1016/j.jff.2019.103453
Grondin, J. M., Tamura, K., Dejean, G., Abbott, D. W., and Brumer, H. (2017). Polysaccharide utilization loci: fueling microbial communities. J. Bacteriol. 199:e00860. doi: 10.1128/JB.00860-16
Gu, I., Lam, W. S., Marasini, D., Brownmiller, C., Savary, B. J., Lee, J. A., et al. (2021). In vitro fecal fermentation patterns of Arabinoxylan from Rice bran on fecal microbiota from normal-weight and overweight/obese subjects. Nutrients 13:2052. doi: 10.3390/nu13062052
Hopkins, M. J., Englyst, H. N., Macfarlane, S., Furrie, E., Macfarlane, G. T., and McBain, A. J. (2003). Degradation of cross-linked and non-cross-linked arabinoxylans by the intestinal microbiota in children. Appl. Environ. Microbiol. 69, 6354–6360. doi: 10.1128/aem.69.11.6354-6360.2003
Johnson, E. M., Lee, H., Jayabalan, R., and Suh, J. W. (2020). Ferulic acid grafted self-assembled fructo-oligosaccharide micro particle for targeted delivery to colon. Carbohydr. Polym. 247:116550. doi: 10.1016/j.carbpol.2020.116550
Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436
Kaur, A., Rose, D. J., Rumpagaporn, P., Patterson, J. A., and Hamaker, B. R. (2011). In vitro batch fecal fermentation comparison of gas and short-chain fatty acid production using "slowly fermentable" dietary fibers. J. Food Sci. 76, H137–H142. doi: 10.1111/j.1750-3841.2011.02172.x
Kaur, A., Singh, B., Yadav, M. P., Bhinder, S., and Singh, N. (2021). Isolation of arabinoxylan and cellulose-rich arabinoxylan from wheat bran of different varieties and their functionalities. Food Hydrocoll. 112:106287. doi: 10.1016/j.foodhyd.2020.106287
Kim, H., Yu, K. W., Hong, H. D., and Shin, K. S. (2017). Effect of arabinoxylan-and rhamnogalacturonan I-rich polysaccharides isolated from young barley leaf on intestinal immunostimulatory activity. J. Funct. Foods 35, 384–390. doi: 10.1016/j.jff.2017.05.052
Kmezik, C., Bonzom, C., Olsson, L., Mazurkewich, S., and Larsbrink, J. (2020). Multimodular fused acetyl-feruloyl esterases from soil and gut Bacteroidetes improve xylanase depolymerization of recalcitrant biomass. Biotechnol. Biofuels 13:60. doi: 10.1186/s13068-020-01698-9
Komeno, M., Hayamizu, H., Fujita, K., and Ashida, H. (2019). Two novel alpha-l-Arabinofuranosidases from Bifidobacterium longum subsp. longum belonging to glycoside hydrolase family 43 cooperatively degrade Arabinan. Appl. Environ. Microbiol. 85:e02582-02518. doi: 10.1128/AEM.02582-18
Lebet, V., Arrigoni, E., and Amado, R. (1998). Digestion procedure using mammalian enzymes to obtain substrates for in vitro fermentation studies. LWT-Food Sci. Technol. 31, 509–515. doi: 10.1006/fstl.1998.0402
Long, W., Xue, Z., Zhang, Q., Feng, Z., Bridgewater, L., Wang, L., et al. (2015). Differential responses of gut microbiota to the same prebiotic formula in oligotrophic and eutrophic batch fermentation systems. Sci. Rep. 5:13469. doi: 10.1038/srep13469
Long, L., Zhao, H., Ding, D., Xu, M., and Ding, S. (2018). Heterologous expression of two aspergillus Niger feruloyl esterases in Trichoderma reesei for the production of ferulic acid from wheat bran. Bioprocess Biosyst. Eng. 41, 593–601. doi: 10.1007/s00449-018-1894-3
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17, 10–12. doi: 10.14806/ej.17.1.200
Mendis, M., Leclerc, E., and Simsek, S. (2016). Arabinoxylans, gut microbiota and immunity. Carbohydr. Polym. 139, 159–166. doi: 10.1016/j.carbpol.2015.11.068
Neyrinck, A. M., Possemiers, S., Druart, C., van de Wiele, T., De Backer, F., Cani, P. D., et al. (2011). Prebiotic effects of wheat Arabinoxylan related to the increase in Bifidobacteria, Roseburia and Bacteroides/Prevotella in diet-induced obese mice. PLoS One 6:e20944. doi: 10.1371/journal.pone.0020944
Nguyen, N. K., Deehan, E. C., Zhang, Z., Jin, M., Baskota, N., Perez-Munoz, M. E., et al. (2020). Gut microbiota modulation with long-chain corn bran arabinoxylan in adults with overweight and obesity is linked to an individualized temporal increase in fecal propionate. Microbiome 8:118. doi: 10.1186/s40168-020-00887-w
Oliveira, D. M., Mota, T. R., Oliva, B., Segato, F., Marchiosi, R., Ferrarese-Filho, O., et al. (2019). Feruloyl esterases: biocatalysts to overcome biomass recalcitrance and for the production of bioactive compounds. Bioresour. Technol. 278, 408–423. doi: 10.1016/j.biortech.2019.01.064
Ou, J. Y., Huang, J. Q., Song, Y., Yao, S. W., Peng, X. C., Wang, M. F., et al. (2016). Feruloylated oligosaccharides from maize bran modulated the gut microbiota in rats. Plant Foods Hum. Nutr. 71, 123–128. doi: 10.1007/s11130-016-0547-4
Paesani, C., Degano, A. L., Salvucci, E., Zalosnik, M. I., Fabi, J. P., Sciarini, L. S., et al. (2020). Soluble arabinoxylans extracted from soft and hard wheat show a differential prebiotic effect in vitro and in vivo. J. Cereal Sci. 93:102956. doi: 10.1016/j.jcs.2020.102956
Pollet, A., Van Craeyveld, V., Van de Wiele, T., Verstraete, W., Delcour, J. A., and Courtin, C. M. (2012). In vitro fermentation of arabinoxylan oligosaccharides and low molecular mass arabinoxylans with different structural properties from wheat (Triticum aestivum L.) bran and psyllium (Plantago ovata Forsk) seed husk. J. Agric. Food Chem. 60, 946–954. doi: 10.1021/jf203820j
Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). Fast tree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077
Riviere, A., Moens, F., Selak, M., Maes, D., Weckx, S., and De Vuyst, L. (2014). The ability of bifidobacteria to degrade arabinoxylan oligosaccharide constituents and derived oligosaccharides is strain dependent. Appl. Environ. Microbiol. 80, 204–217. doi: 10.1128/AEM.02853-13
Rose, D. J., Patterson, J. A., and Hamaker, B. R. (2010). Structural differences among alkali-soluble arabinoxylans from maize (Zea mays), rice (Oryza sativa), and wheat (Triticum aestivum) brans influence human fecal fermentation profiles. J. Agric. Food Chem. 58, 493–499. doi: 10.1021/jf9020416
Rumpagaporn, P., Reuhs, B. L., Kaur, A., Patterson, J. A., Keshavarzian, A., and Hamaker, B. R. (2015). Structural features of soluble cereal arabinoxylan fibers associated with a slow rate of in vitro fermentation by human fecal microbiota. Carbohydr. Polym. 130, 191–197. doi: 10.1016/j.carbpol.2015.04.041
Snelders, J., Dornez, E., Delcour, J. A., and Courtin, C. M. (2013). Ferulic acid content and appearance determine the antioxidant capacity of arabinoxylanoligosaccharides. J. Agric. Food Chem. 61, 10173–10182. doi: 10.1021/jf403160x
Snelders, J., Olaerts, H., Dornez, E., Van de Wiele, T., Aura, A.-M., Vanhaecke, L., et al. (2014). Structural features and feruloylation modulate the fermentability and evolution of antioxidant properties of arabinoxylanoligosaccharides during in vitro fermentation by human gut derived microbiota. J. Funct. Foods 10, 1–12. doi: 10.1016/j.jff.2014.05.011
Song, Y., Wu, M. S., Tao, G., Lu, M. W., Lin, J., and Huang, J. Q. (2020). Feruloylated oligosaccharides and ferulic acid alter gut microbiome to alleviate diabetic syndrome. Food Res. Int. 137:109410. doi: 10.1016/j.foodres.2020.109410
Tiwari, U. P., Singh, A. K., and Jha, R. (2019). Fermentation characteristics of resistant starch, arabinoxylan, and beta-glucan and their effects on the gut microbial ecology of pigs: a review. Anim. Nutr. 5, 217–226. doi: 10.1016/j.aninu.2019.04.003
Vacca, M., Celano, G., Calabrese, F. M., Portincasa, P., Gobbetti, M., and De Angelis, M. (2020). The controversial role of human gut Lachnospiraceae. Microorganisms 8:573. doi: 10.3390/microorganisms8040573
Wang, J., Bai, J. Y., Fan, M. C., Li, T. T., Li, Y., Qian, H. F., et al. (2020). Cereal-derived arabinoxylans: structural features and structure–activity correlations. Trends Food Sci. Technol. 96, 157–165. doi: 10.1016/j.tifs.2019.12.016
Xie, J. H., Shen, M. Y., Nie, S. P., Liu, X., Zhang, H., and Xie, M. Y. (2013). Analysis of monosaccharide composition of Cyclocarya paliurus polysaccharide with anion exchange chromatography. Carbohydr. Polym. 98, 976–981. doi: 10.1016/j.carbpol.2013.07.011
Xu, H. D., Reuhs, B. L., Cantu-Jungles, T. M., Tuncil, Y. E., Kaur, A., Terekhov, A., et al. (2022). Corn arabinoxylan has a repeating structure of subunits of high branch complexity with slow gut microbiota fermentation. Carbohydr. Polym. 289:119435. doi: 10.1016/j.carbpol.2022.119435
Yu, M., Arioglu-Tuncil, S., Xie, Z., Fu, X., Huang, Q., Chen, T., et al. (2021). In vitro fecal fermentation profiles and microbiota responses of pulse cell wall polysaccharides: enterotype effect. Food Funct. 12, 8376–8385. doi: 10.1039/d1fo01231b
Zhang, M., Chekan, J. R., Dodd, D., Hong, P. Y., Radlinski, L., Revindran, V., et al. (2014). Xylan utilization in human gut commensal bacteria is orchestrated by unique modular organization of polysaccharide-degrading enzymes. Proc. Natl. Acad. Sci. U. S. A. 111, E3708–E3717. doi: 10.1073/pnas.1406156111
Zhang, X., Chen, T., Lim, J., Xie, J., Zhang, B., Yao, T., et al. (2019). Fabrication of a soluble crosslinked corn bran arabinoxylan matrix supports a shift to butyrogenic gut bacteria. Food Funct. 10, 4497–4504. doi: 10.1039/c8fo02575d
Keywords: arabinoxylan, feruloylation, carbohydrate structure, gut microbiota, in vitro fermentation
Citation: Li Z, Zhang H, He L, Hou Y, Che Y, Liu T, Xiong S, Zhang X, Luo S, Liu C and Chen T (2023) Influence of structural features and feruloylation on fermentability and ability to modulate gut microbiota of arabinoxylan in in vitro fermentation. Front. Microbiol. 13:1113601. doi: 10.3389/fmicb.2022.1113601
Edited by:
Xiaowei Zhang, University of Shanghai for Science and Technology, ChinaReviewed by:
Tianming Yao, Purdue University, United StatesCheng Li, University of Shanghai for Science and Technology, China
Copyright © 2023 Li, Zhang, He, Hou, Che, Liu, Xiong, Zhang, Luo, Liu and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tingting Chen, ✉ dGluZ3RpbmdjaGVuQG5jdS5lZHUuY24=; Chengmei Liu, ✉ bGl1Y2hlbmdtZWlAYWxpeXVuLmNvbQ==
†These authors have contributed equally to this work