- 1College of Animal Science and Technology, Northeast Agricultural University, Harbin, China
- 2College of Food Science, Northeast Agricultural University, Harbin, China
- 3Lethbridge Research and Development Centre, Lethbridge, AB, Canada
- 4School of Nursing and School of Public Health, Yangzhou University, Yangzhou, China
The aim of this study was to investigate the gas production (GP), dry matter disappearance (DMD), fermentation parameters, and rumen microbiota in response to Clostridium butyricum (CB) supplementation in batch culture using a high forage substrate. The doses of CB were supplemented at 0 (Control), 0.5 × 106, 1 × 106, and 2 × 106 CFU/bottle, respectively, at either media pH 6.0 or pH 6.6. The 16S rRNA gene sequencing was used to detect the microbiota of fermentation culture in control and 1 × 106 CFU/bottle after 24 h of incubation. The results showed that the GP (p < 0.001), DMD (p = 0.008), total volatile fatty acid (VFA) concentration (p < 0.001), acetate to propionate ratio (p < 0.001), and NH3-N concentration (p < 0.001) were greater at media pH 6.6 than pH 6.0. Furthermore, the linearly increased DMD (pH 6.0, p = 0.002; pH 6.6, p < 0.001) and quadratically increased butyrate proportion (pH 6.0, p = 0.076; pH 6.6, p < 0.053) and NH3-N concentration (pH 6.0, p = 0.003; pH 6.6, p = 0.014) were observed with increasing doses of CB. The Alpha diversity indexes of OTU number and Chao1 were higher (p = 0.045) at media pH 6.6 than pH 6.0, but they were not affected by CB supplementation. The PCoA analysis (unweighted uniFrac) demonstrated that the clustering of the bacterial microbiota of control and CB were distinctly separated from each other at media pH 6.0. At the phylum level, the abundance of Bacteroidota (p < 0.001) decreased, whereas that of Firmicutes (p = 0.026) increased when the media pH was elevated from 6.0 to 6.6. Supplementation of CB increased relative abundances of Rikenellaceae_RC9_gut_group (p = 0.002), Christensenellaceae_R-7_group (p < 0.001), and NK4A214_group (p = 0.002) at genus level. Interactions between media pH and CB addition were observed for bacteria at both phylum and genus levels. These results indicated that increasing the media pH level and CB supplementation increased in vitro rumen digestibility, and altered the ruminal fermentation pattern (by media pH) and microbiota.
Introduction
The use of antibiotics has contributed to the improvement of animal production and animal health for several decades. Effective antibiotic growth promoters have been commonly fed to livestock to improve feed efficiency and animal health (Lana et al., 1997; Ítavo et al., 2011). However, there are various concerns involving the use of antibiotics including environmental pollution, food safety issues, and risks of developing antibiotic resistance in microbes associated with animal or human diseases. Consequently, discussion around the potential ban of antibiotics has become a trend, but this would impose a considerable challenge for livestock production. Hence, it is crucial to develop economically competitive alternatives (e.g., probiotics, prebiotics, or plant extracts) to antibiotics to meet the requirements of improving animal health and products (Jiao et al., 2017, 2021).
Probiotics, defined as “live microorganisms,” have been considered as potential alternatives to antibiotics to improve growth performance and immunity of pigs (Roselli et al., 2017), sheep (Mani et al., 2021), beef cattle (Ran et al., 2018), or dairy cows (Chiquett et al., 2012). In ruminants, probiotics are considered oxygen consumers or nutrient providers to microorganisms in the rumen for the modulation of ruminal fermentation, and inhabitation of detrimental microorganisms by producing some antimicrobial substances like bacteriocins (Chaucheyras-Durand et al., 2008; McAllister et al., 2011), thereby improving growth performance. Our previous studies confirmed the beneficial effect of probiotics with yeast or lactic acid bacteria on stimulating ruminal fermentation and increasing the copy numbers of cellulolytic bacteria. However, the effects of probiotics varied considerably among different experimental studies, and these differences are likely due to a number of factors that can impact the final outcomes, such as type of probiotics, supplementation dose, media pH level, and type of substrate used (Jiao et al., 2017, 2018, 2019).
As a gram-positive anaerobe, Clostridium butyricum (CB) is known as a producer of butyric acid that forms spores, and is widely used in monogastric feeding to improve growth performance, feed efficiency, and animal health (Wang et al., 2020; Zhang et al., 2020). Because CB can survive in low pH and high bile concentrations, previous studies mainly focused on its role in the modulation of gut microbial composition and improvement of intestinal health (Cao et al., 2019; Li et al., 2019). In ruminants, a recent study conducted by Li et al. (2021) showed that dietary supplementation of CB improved growth performance and enhanced ruminal fermentation by increasing the abundance of rumen microbiota and molar proportion of propionate and butyrate in Holstein heifers. The common diet of lactating dairy cows usually involves highly fermatable concentrate and finely chopped silage (i.e., low physical effectiveness) to meet the energy requirement for high milk production. Therefore, low rumen pH (<5.8) could be observed, particularly following feed ingestion. Similarly, Cai et al. (2021a) reported that the CB supplementation increased total VFA concentration and improved growth performance of heat-stressed goats. Although those studies (Cai et al., 2021a; Li et al., 2021) have investigated the effects of CB supplementation on rumen fermentation characteristics, ruminal pH was at a high level (≥6.5), and limited information on rumen microbiota was obtained. We have identified a research gap when it comes to the impact of CB supplementation on gas production (GP) and rumen microbiota, particularly with varying supplementation doses and ruminal pH level. We hypothesized that supplementing C. butyricum would increase dry matter disappearance (DMD) and rumen fermentation characteristics, and improve rumen microbiota at both phylum and genus levels, depending on media pH level and supplemental dose. Hence, the objective of this study was to evaluate the effects of CB supplementation varying with media pH levels (6.0 or 6.6) and supplemental doses on gas production (GP), fermentation parameters, and microbiota in batch culture.
Materials and methods
Clostridium butyricum resource
The CB strain evaluated in this study was obtained from Greensnow Biological Biotechnology Co., Ltd. (Wuhan, China). The concentration of the CB sample was 2.0 × 108 CFU/g.
Experimental design, substrate, and inoculum
The experiment was conducted using a completely randomized design with two media pH levels (6.0 and 6.6) × 4 supplemental dosages of CB at control (0 CFU/bottle), low (0.5 × 106 CFU/bottle), medium (1 × 106 CFU/bottle), and high (2 × 106 CFU/bottle). The media pH of 6.0 and 6.6 were chosen based on the averaged rumen pH of dairy cow diet (Palmonari et al., 2010). The supplemental doses of CB were outlined according to the manufacturer’s recommendations. A typical dairy cow diet was used as substrate (Table 1). Rumen inoculum was provided by two ruminally fistulated nonlactating dairy cows fed a diet based on 55% forage and 45% concentrate (DM basis). A total mixed ration was prepared daily and offered twice at 06:00 h in the morning and 18:00 h in the afternoon. All experimental procedures involving animals were performed according to the experimental license (protocol number: NEAU-[2011]-9) of Northeast Agricultural University (Harbin, China).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-t001.jpg)
Table 1. Ingredients and chemical composition of experimental diets.a
Batch culture procedures
The batch culture was conducted using 125 ml glass bottles, and rubber stoppers and aluminum caps were used to seal the bottles. Approximately 0.5 g (DM basis) of ground substrate was passed through a 1-mm screen and weighed into the acetone-washed and preweighed filter bags (F57; Ankom Technology, Macedon, NY, United States). The CB samples were added into the bottles in three replicates at desired dosages. Rumen contents were collected 2 h before the morning feeding from four different locations within the rumen, composited and squeezed through a nylon mesh. The strained rumen fluid was kept in an insulated air tight container. The rumen fluid was re-strained through four layers of cheesecloth to remove particles, and warmed at 39°C in a water bath in the laboratory. The pH of rumen fluid was 6.74, 6.81, and 6.69 for the first, second, and third run, respectively. Forty-five milliliters of freshly prepared warm buffer and 15 ml of strained rumen fluid were added into each bottle with continuously flushing of carbon dioxide into the bottle to remove the air from the headspace. The bottles were sealed with rubber stoppers and aluminum caps, and then transferred to a shaking incubator (SPH-2102C, Shanghai Huyueming Scientific Instrument Co., Ltd., Shanghai, China) at 39°C for 24 h. The two media pH levels (6.0 and 6.6) were achieved according to the method described by Jiao et al. (2019) that was adjusted to the volume of sodium bicarbonate in the buffer solution. At each media pH level, three additional bottles containing the same volume of buffer and rumen fluid without substrate or CB additive were used as blanks for correcting GP due to fermentation media. The same batch culture was repeated weekly for another two runs and a run was used as an experimental replicate.
Gas pressure was recorded at 3, 6, 9, 12, and 24 h of incubation using a pressure transducer (model HT-935, Hongcheng Technology, Shenzhen, China) with a 23-gauge needle (0.6 mm) passing through rubber stoppers. Pressure values, corrected for the amount of substrate DM incubated and the gas released from the blanks, were used to generate gas volume (GV) using equation of Romero-Perez and Beauchemin (2018): GV = 4.7047 × (gas pressure) + 0.0512 × (gas pressure2). After 24 h of incubation, the bottles were placed in cold water to stop fermentation. The pH of the fermentation fluid was measured using a portable pH meter (PHS-3C; Nanjing Nanda Analytical Instrument Application Research Institute, Nanjing, China) after removing the aluminum caps and rubber stoppers. Samples of 5 ml of fermentation fluid were collected and preserved with 1 ml of metaphosphoric acid (25%, w/v) and sulfuric acid (1%, v/v), and stored at −20°C for the late analysis of VFA profile and NH3-N concentration. The fermentation liquid samples from three bottles of each treatment combination were pooled by run, and only the control and medium dose of CB were used for DNA extract. Two samples of 5 ml of fermentation liquid were collected and immediately frozen in liquid nitrogen and preserved at −80°C until DNA extraction.
Microbial DNA isolation and 16S rRNA gene sequencing
The rumen fluid samples for determining bacterial community were delivered to LC-Sciences (Hangzhou, China) for 16S rRNA gene sequencing. In brief, DNA from fermentation fluid was extracted using the E.Z.N.A. ®Stool DNA Kit (D4015, Omega, Inc., United States) according to the manufacturer’s instructions. The V4–V5 region of the bacteria 16S rRNA gene was amplified by PCR. The amplification was performed using the forward primer 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and the reverse primer 907R (5′-CCGTCAATTCMTTTRAGTTT-3′). The PCR reactions were performed in triplicate in a 25 μl mixture containing 25 ng of template DNA, 12.5 μl PCR Premix, 2.5 μl of each primer, and PCR-grade water to 25 μl. The PCR reactions were conducted using the following program: (i) initial denaturation at 98°C for 30 s; (ii) 32 cycles of denaturation at 98°C for 10 s, (iii) annealing at 54°C for 30 s, and (iv) extension at 72°C for 45 s; and (v) a final extension step at 72°C for 10 min. The amplified products were detected using 2% agarose gel electrophoresis. The PCR products were purified with AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, United States) and quantified using Qubit (Invitrogen, Carlsbad, United States). The pooled purified libraries were measured for size and quantity with Agilent 2100 Bioanalyzer (Agilent, United States) and further quantified using the Library Quantification Kit for Illumina (Kapa Biosciences, Woburn, MA, United States). Sequencing was performed on NovaSeq PE250 platform.
Paired-end reads were assigned to samples based on their unique barcode and truncated by cutting off the barcode and primer sequence. Paired-end reads were merged using FLASH (version 1.2.8).1 Quality filtering on the raw reads was performed under specific filtering conditions to obtain the high-quality clean tags according to the fqtrim (version 0.94).1 Chimeric sequences were filtered using Vsearch software (version 2.3.4).1 Alpha diversity and beta diversity were calculated by being normalized to the same sequences randomly. Then, according to SILVA (Release 138)1 classifier, feature abundance was normalized using the relative abundance of each sample. Alpha diversity is applied in analyzing the complexity of species diversity for a sample through five indices, including Chao1, OTU, Evenness, Shannon, Simpson, and all these indices in the samples were calculated with QIIME2.1 Beta diversity was calculated by QIIME2. The blast was used for sequence alignment, and the feature sequences were annotated with the SILVA database for each representative sequence.
Chemical analyses
Substrate was analyzed for DM (method 930.15), ash and organic matter (method 942.05), and crude protein (method 990.03; N × 6.25) according to AOAC (2005). Neutral detergent fibre was determined as described by Van Soest et al. (1991) with heat-stable α-amylase and sodium sulfite used. Acid detergent fibre was determined according to AOAC (2005), method 973.18. The VFA concentration was determined using a gas chromatograph (GC-2010, Shimadzu) equipped with a capillary column (30 m × 0.32 mm i.d., 1-μm phase thickness, Zebron ZB-FAAP, Phenomenex, Torrance, CA, United States) according to Beauchemin et al. (2003). The NH3-N concentrations were determined using the method described by Broderick and Kang (1980).
Statistical analysis
Data for GP, DMD, and fermentation characteristics were subjected to statistical analysis using the MIXED procedure of SAS (SAS Inst. Inc. Cary, NC, United States) with the model including fixed effects of media pH, a supplemental dose of CB, the interactions, and the random effects of the run. Data for microbiota were analyzed using the MIXED procedure of SAS (SAS Inst. Inc. Cary, NC, United States) with the model including fixed effects of media pH, CB supplementation (control vs. medium dose), the interactions, and the random effects of the run. Tukey’s multiple comparison test was used to examine the significance among treatments within a media pH level for the data of GP, DMD, and fermentation characteristics as well as for the microbiota data. The CONTRAST statement of SAS with linear and quadratic orthogonal contrasts was used to determine the effect of increasing the dose of CB. Differences were declared significant at p ≤ 0.05 and trends were discussed at 0.05 < p ≤ 0.10 unless otherwise stated.
Results
GP, DMD, and rumen fermentation characteristics
There was no interaction between media pH level and dose of CB on GP at 3, 6, 9, and 24 h, but the interaction was significant for the GP at 12 h (p = 0.029; Table 2); the GP increased with increasing CB supplementation at 9 h (linear; p = 0.011), 12 h (linear; p = 0.002), and 24 h (quadratic; p < 0.001) at media pH 6.0; whereas it was not affected by increasing CB supplementation except at 24 h (quadratic; p = 0.044) at media pH 6.6. Moreover, at media pH 6.6, the GP was greater with supplementation of CB at low and high doses compared with the control and medium doses at 3 h (p = 0.022), 12 h (p = 0.045), and 24 h (p = 0.009) of incubation. Linearly (p < 0.001) increasing CB addition also increased DMD without interaction of pH level with the CB supplemental dose. Overall, the GP (p < 0.001) and DMD (p = 0.008) were greater at media pH 6.6 vs. pH 6.0.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-t002.jpg)
Table 2. Effects of supplemental doses of Clostridium butyricum (CB) on gas production (GP), DM disappearance (DMD), and rumen fermentation characteristics at low and high media pH after 24 h of in vitro fermentation.
In comparison with media pH 6.0, the greater (p < 0.001) total VFA concentration (averaged 60.3 vs. 53.0 mM), molar proportions of acetate (averaged 54.2 vs. 53.2%), A:P ratio (averaged 2.49 vs. 2.31) and NH3-N concentration (averaged 14.5 vs. 13.6 mM), and lower (p < 0.001) molar proportion of propionate (21.8 vs. 23.1%) were observed at media pH 6.6 (Table 2). Furthermore, there was no interaction between pH level and CB supplementation on total VFA concentration, individual VFA profiles, and NH3-N concentration. Increasing the supplementation of CB did not affect the total VFA concentrations and the molar proportion of acetate at either media pH 6.0 or pH 6.6. However, increasing the supplemental dose of CB quadratically changed the proportion of propionate (p = 0.055 at pH 6.0; p = 0.039 at pH 6.6), and linearly increased (p = 0.027) the molar proportion of butyrate at pH 6.0. The fermentation pattern, expressed as A:P ratio was not influenced by supplementing CB. It was also observed that the NH3-N concentrations linearly increased (pH 6.0, p < 0.001; pH 6.6, p = 0.008) with increasing supplementation of CB.
Microbiome diversity analysis
A total of 2,028,448 raw tags with an average of 84,519 raw tags for each sample were generated by 16S rRNA gene sequencing. In addition, a total of 1,854,300 valid tags with an average of 77,263 valid tags for each sample were obtained, and the ratio of valid data to raw data was 89.01–93.00%. The Alpha diversity results showed that the indexes of OTU number (p = 0.045) and Chao1 (p = 0.044) were lower, whereas Simpson (p < 0.001) and Evenness (p = 0.005) were greater at media pH 6.0 than pH 6.6 (Table 3). Supplementation of CB increased the indexes of OTU number (p = 0.047) and Chao1 (p = 0.048), and decreased Simpson (p = 0.008) and Evenness (p = 0.028) at media pH 6.6 than pH 6.0. All indexes of Alpha diversity were unaffected with supplementation of CB at both media pH levels. The interaction between pH and CB tended to be significant for OTU number (p = 0.052) and Chao 1 (p = 0.056). A multiple comparison did not clearly show the synergism effect with the combination of high media pH and CB supplementation compared with other treatment combinations.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-t003.jpg)
Table 3. Bacterial alpha diversity indexes of fermentation fluid supplemented with CB at low or high media pH after 24 h of in vitro incubation.
As for the Beta diversity of the bacterial communities in the fermentation culture, the PcoA analysis (unweighted uniFrac) illustrated that the clustering of bacterial microbiota from control and CB overlapped and no clear distinction was observed at media pH 6.6. Differently, clustering of bacterial microbiota of control and CB were distinctly separated from each other at media pH 6.0 (Figure 1). The Venn diagram in the fermentation culture samples showed that the control and CB groups shared 2,971 OTUs, and had 2,164 and 1,963 exclusive OTUs (respectively) media pH 6.0 (Figure 2). Moreover, there were 3,131 shared OTUs, and 1,937 and 2,355 exclusive OTUs for the control and CB groups at media pH 6.6.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-g001.jpg)
Figure 1. Principal coordinates analysis (PCoA) based on unweighted uniFrac distance matrices of microbiota in ruminal fermentation fluid supplemented with CB at low or high media pH after 24 h of in vitro incubation. CONH, control at media pH 6.6; CONL, control at media pH 6.0; CBH, CB supplemented at 1 × 106 CFU/bottle at media pH 6.6; and CBL, CB supplemented at 1 × 106 CFU/bottle at media pH 6.0.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-g002.jpg)
Figure 2. Analysis of OTUs exclusiveness of microbiota in ruminal fermentation fluid supplemented with or without CB at low or high media pH after 24 h of in vitro incubation. CONH, control at media pH 6.6; CONL, control at media pH 6.0; CBH, CB supplemented at 1 × 106 CFU/bottle at media pH 6.6; and CBL, CB supplemented at 1 × 106 CFU/bottle at media pH 6.0.
Bacterial community
The relative abundances of the 10 most abundant bacterial phyla were detected in the fermentation culture of control and CB (medium dose) groups (Table 4). The relative abundances of Bacteroidota (p < 0.001), Bacteroidetes (p = 0.004), and Desulfobacterota (p = 0.001) were lower, but Firmicutes (p = 0.026), Proteobacteria (p = 0.035), Verrucomicrobiota (p < 0.001), and Cyanobacteria (p < 0.001) were greater at media pH 6.6 than pH 6.0. At media pH 6.0, the relative abundances of Bacteroidota (p < 0.001), Bacteroidetes (p = 0.031), and Desulfobacterota (p = 0.001) were greater, whereas Verrucomicrobiota (p < 0.001), Spirochaetota (p < 0.001), and Cyanobacteria (p < 0.001) were lower with supplementation of CB compared with media pH 6.6. Supplementation of CB decreased the relative abundances of Proteobacteria (p < 0.001) and Verrucomicrobiota (p = 0.004), and increased Desulfobacterota (p = 0.023) and Cyanobacteria (p = 0.003). Interaction (p < 0.001) between pH and CB were observed for Spirochaetota. A multiple comparison among four means did not find a consistent optimal pH and CB combination.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-t004.jpg)
Table 4. The relative abundance (%) of bacterial phylum in fermentation fluid supplemented with CB and without CB (Control) at low or high media pH after 24 h of in vitro incubation.
As shown in Table 5, the relative abundances of 27 predominant bacterial genera (>1% at least in one group) were detected in the fermentation culture. There were interactions for Christensenellaceae_R-7_group (p = 0.001), NK4A214_group (p = 0.003), Succinivibrionaceae_UCG-002 (p = 0.030), and Muribaculaceae_unclassified (p < 0.001). Furthermore, adding CB increased Christensenellaceae_R-7_group (p = 0.001), NK4A214_group (p = 0.003), and Muribaculaceae_unclassified (p < 0.001), and decreased Succinivibrionaceae_UCG-002 (p = 0.002) at media pH 6.0, but they were not affected at media pH 6.6. In comparison with media pH 6.6, the relative abundances of Prevotella (p < 0.001), Rikenellaceae_RC9_gut_group (p < 0.001), Christensenellaceae_R-7_group (p < 0.001), NK4A214_group (p < 0.001), Ruminococcus (p = 0.040), Muribaculaceae_unclassified (p < 0.001), Prevotellaceae_UCG-001 (p < 0.001), Clostridia_UCG-014_unclassified (p < 0.001), Bacteroidales_unclassified (p < 0.001), Bacteroidetes_unclassified (p < 0.001) and Desulfovibrio (p = 0.001) were greater, whereas CAG-352 (p < 0.001), WCHB1-41_unclassified (p < 0.001), Bacteroidales_RF16_group_unclassified (p = 0.003), Firmicutes_unclassified (p = 0.001), Ruminococcaceae_unclassified (p < 0.001), Treponema (p < 0.001), and UCG-010_unclassified (p < 0.001) were lower at media pH 6.0. The relative abundances of Prevotella (p < 0.001), Rikenellaceae_RC9_gut_group (p < 0.001), Christensenellaceae_R-7_group (p < 0.001), NK4A214_group (p < 0.001), Muribaculaceae_unclassified (p < 0.001), Prevotellaceae_UCG-001 (p = 0.002), Clostridia_UCG-014_unclassified (p < 0.001), Bacteroidales_unclassified (p < 0.001), Acinetobacter (p = 0.008), Bacteroidetes_unclassified (p < 0.001), and Desulfovibrio (p = 0.001) were greater, whereas CAG-352 (p < 0.001), WCHB1-41_unclassified (p < 0.001), Bacteroidales_RF16_group_unclassified (p = 0.005), Firmicutes_unclassified (p = 0.005), Ruminococcaceae_unclassified (p < 0.001), Treponema (p < 0.001), and UCG-010_unclassified (p < 0.001) were lower with supplementation at media pH 6.0 than pH 6.6. The CB supplementation increased the relative abundances of Rikenellaceae_RC9_gut_group (p = 0.002), Clostridia_UCG-014_unclassified (p = 0.042), Bacteroidales_unclassified (p = 0.047) and Bacteroidetes_unclassified (p = 0.028), and Desulfovibrio (p = 0.022), and decreased WCHB1-41_unclassified (p = 0.007), and Bacteroidales_RF16_group_unclassified (p = 0.030) and Prevotellaceae_UCG-003 (p = 0.036). A multiple comparison among four means was unable to conclude a consistent optimal pH and CB combination on the relative abundance of bacterial genus.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/960623/fmicb-13-960623-HTML/image_m/fmicb-13-960623-t005.jpg)
Table 5. The relative abundance (%) of bacterial genus (>1% at least in one group) in ruminal fermentation fluid supplemented with CB and without CB (control) at low or high media pH after 24 h of vitro incubation.
Discussion
Gas production, DMD, and fermentation characteristics
The ruminal pH dynamics and pH maximum occurred just before feeding, and the pH minimum was detected approximately 6–7 h after feeding. The mean pH ranged from 5.89 to 6.63 in the rumen when the dairy cows were fed a corn-silage based total mixed ration (Palmonari et al., 2010). Russell and Wilson (1996) suggested that at ruminal pH below 6.0 could be considered a threshold value for subacute rumen acidosis. Therefore, the target media pH was selected at 6.0 and 6.6 for low and high media pH, respectively, also in accordance with the guidelines of our previous in vitro studies with pH 5.8 and 6.5 for low and high media pH (Jiao et al., 2019, 2021). In these studies, the low media pH was set at 5.8 because the high-grain substrate (>85%, DM basis) was used and lower pH was expected. In the present study, the final media pH was 6.17 and 6.69, which were close to our target media low pH of 6.0 and high pH of 6.6, indicating high buffering of fermentation media despite fermentation acids produced from the substrate. The overall greater GP, DMD, and total VFA concentration with pH 6.6 than pH 6.0 in the present study confirmed the increased microbial activity and feed degradability with high rumen pH. In fact, the higher proportion of acetate and lower proportion of propionate that resulted in a higher A:P ratio with media pH 6.6 vs. 6.0 suggested an improvement in fibre digestion. The results were consistent with the experiment conducted by Jiao et al. (2019), who observed an increase in the GP, DMD, total VFA concentration, and A:P ratio from low (5.8) to high media pH (6.5) in batch culture using a diet containing 60% forage. In addition, our outcomes revealed that phylum Proteobacteria involved in fibre digestion (Qiu et al., 2020; Wei et al., 2022) and its relative abundance was decreased at media pH 6.0 compared with pH 6.6, indicating an inhibition effect of low media pH on fibre digestion. Thus, the varied results of GP, DMD and fermentation characteristics could be attributed to the inhibition of growth of phylum Proteobacteria at low media pH vs. high media pH.
Cai et al. (2021a) reported that supplementation of CB at 2 × 104 and 4 × 104 CFU/bottle increased the in vitro digestibilities of DM, NDF, and ADF using rumen inoculum from goats fed a diet containing 56% forage. Our study presents similar results and detected a linear increase in the accumulative GP and DMD at 24 h after increasing the supplemental dose of CB at either media pH 6.0 or pH 6.6. The improved digestibility by adding CB may be explained by its ability to provide short-chain fatty acids, nutritional factors such as enzymes (exo-pectate lyase, pectin methylesterase, and endopectatelyase) and vitamins B and E that could improve the activities of cellulolytic enzymes (e.g., avicelase, CMCase, cellobiase, and xylanase) in the rumen (Cai et al., 2021a,b). Furthermore, the increased butyrate proportion due to CB addition at media pH 6.0 (linear) and pH 6.6 (quadratic) suggested in this study could be partially explained by the butyrate produced by CB. In addition, although the proportion of propionate quadratically changed after increasing the supplemented CB at both media pH levels, its differences among treatments were minimal from a biological perspective. Li et al. (2021) reported increased proportions of propionate and butyrate without changing the proportion of acetate in Holstein heifers supplemented with CB at 2 × 108 CFU/kg dietary DM. Cai et al. (2021b) demonstrated an increase in acetate concentration with dietary supplementation of CB (1 × 108 CFU/g) at 0.05% of the basal diet in heat-stressed goats. The inconsistent effects of CB addition on the VFA profiles found in different studies may be attributed to the different supplemental doses, animal species, animal health status, or experimental approaches (e.g., in vitro vs. in vivo) used in each study. Heat-stressed animals often decrease feed intake, which may positively affect fibre digestibility (i.e., an increase of acetate concentration) due to potentially longer rumen retention. Moreover, the increased NH3-N concentration by supplementation of CB regardless of media pH levels is consistent with the results presented by Cai et al. (2021b) who found an increase of NH3-N concentration with supplementation of CB in the diet of goats. The results suggested promotion of protein degradation as increasing DMD was observed after CB supplementation, but a better microbial protein synthesis was unlikely to happen because the GP and VFA concentration were not lower, and NH3-N concentration kept surging (Blümmel and Ørskov, 1993).
Ruminal bacterial microbiota
Compared with the moderate high forage diet, feeding ruminant with a high concentrate diet would affect fermentation characteristics and the bacterial community, leading to decreased bacterial alpha diversity indexes and rumen pH (Hook et al., 2011; Ma et al., 2021). Zhang et al. (2017) found that the rumen pH linearly decreased when the concentrate level was gradually increased from 20 to 80% in the diet of Holstein heifers, resulting in the decreased bacterial alpha diversity indexes of observed species, PD whole tree, and Pielou in the rumen at lower pH. Consistent with previous studies, our results showed that the bacterial alpha diversity indexes measured as OTU number and Chao 1 increased at media pH 6.6 compared with pH 6.0. Moreover, the PCoA unweighted uniFrac analysis showed that the control groups significantly differed between media pH 6.0 and pH 6.6. The current results of bacterial diversity could be explained by the inhibition of some gram-negative bacteria which were sensitive to the low rumen pH (Khafipour et al., 2009). However, it is important to acknowledge that data on CB supplementation on rumen bacterial diversity is limited. It has been reported that the bacterial diversity was not affected by probiotic supplementation, such as lactic acid bacteria in the diet of sheep (Mani et al., 2021) or active dried yeast in the diet of dairy cows (AlZahal et al., 2017).
The predominant bacterial phyla in the rumen were Bacteroidota and Firmicutes regardless of media pH levels or CB supplementation involved in the present study. Our result are aligned with the findings presented by Li et al. (2022) and Liu et al. (2022) who reported that Bacteroidota and Firmicutes were the dominant bacteria in the rumen of dairy cows. We also identified that 8 in 10 abundances of bacterial phyla and 18 in 27 abundances of bacterial genus in the culture media were affected by media pH levels, which may be partly explained by the death and lysis of Gram-negative bacteria caused by low media pH (Nagaraja and Titgemeyer, 2007). Similar studies have reported that the increased Firmicutes:Bacteroidota ratio in the rumen would favor the capacity for harvesting energy from the diet and has been linked to milk fat yield in dairy cows (Lima et al., 2014; Tong et al., 2022). The increased relative abundance of Firmicutes and decreased Bacteroidota at media pH 6.6 increased Firmicutes:Bacteroidota ratio compared with media pH 6.0, which may contribute to the greater DMD and total VFA concentration at media pH 6.6 than pH 6.0 in the present study.
Several studies demonstrated that supplementation of CB in the diets had an effect on the modulation of bacterial communities in the feces of weaned piglets (Liang et al., 2021), intestine of chickens (Huang et al., 2019), or rumen and rectum of goats (Zhang et al., 2022). In the present study, the relative abundances of five phyla (10 in total) and 12 genera (27 in total) were changed with the supplementation of CB, which may be due to the nutritional factors that change with supplementation (Li et al., 2021). Although the relative abundance of Bacteroidota did not alter, the Rikenellaceae_RC9_gut_group, which belongs to this phylum, increased with CB supplementation. Previous studies reported that Rikenellaceae_RC9_gut_group contains the genes for secreting cellulase and hemicellulase, resulting in an improvement of energy gain from the fibre degradation (Escobar et al., 2014; Fan et al., 2020). In addition, the surging relative abundance of Christensenellaceae_R-7_group with increased the CB supplementation may explain the elevated DMD in the current study. In fact, Huang et al. (2021) reported that the increased abundance of Christensenellaceae_R-7_group favored rumen development by increasing nipple width, epithelial thickness, and stratum corneum thickness and contributed to the absorption and digestion of nutrients in the rumen. Moreover, as cellulose-degrading bacteria, the increased abundance of NK4A214_group by supplementation of CB suggested an improvement of cellulose degradation that explained the increased DMD (Guo et al., 2020; Yu et al., 2020). It is worth mentioning that there were interactions between media pH and CB supplementation for some species of bacteria either in phylum (2/10) and genus (9/27), indicating that the effects of CB could be different when rumen pH differs, particularly when the response to CB has different direction depending on the media pH levels. Therefore, it can be speculated that the mode of action of CB at different pH levels may be different for some bacteria which are sensitive to media pH levels.
Conclusion
This study showed that the GP, DMD, total VFA concentrations, and A:P ratio increased by elevating media pH from 6.0 to 6.6, indicating stimulation of rumen fermentation, which may be partially due to changes in their bacterial diversity indexes and relative abundance associated with fibre degradation at phylum or genus level. Linearly increasing supplementation of CB increased DMD and quadratically changed the proportions of propionate and butyrate and NH3-N concentration at both media pH levels, revealing the dose-dependent traits of the CB effect. The supplementation of CB affected the rumen microbiota at either phylum or genus level, and the increased abundances of Rikenellaceae_RC9_gut_group, Christensenellaceae_R-7_group, and NK4A214_group were considered favorable to cellulose degradation, therefore increasing the DMD of the diet containing moderate high forage. Furthermore, the interaction of the effects of CB with media pH level on the abundance of microbiota is warranted to search for the optimum rumen pH to maximize CB activity.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA844474.
Author contributions
MZ and PJ conceived and designed the experiments. MZ, GL, XZ, XL, SL, and XW performed the experiments. MZ, WY, and YY analyzed the data. MZ and PJ contributed to the writing of the manuscript. All authors reviewed and approved the manuscript.
Funding
This research was supported by the Natural Science Foundation of Heilongjiang Province (YQ2021C018), Postdoctoral Foundation of Heilongjiang Province (LBH-Z21100 and LBH-Z21042), Young Talents Project of Northeast Agricultural University (19QC42), and the University Nursing Program for Young Scholars with Creative Talents in Heilongjiang Province (UNPYSCT-2020106).
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.
Footnotes
1. ^http://ccb.jhu.edu/software/FLASH/
1. ^http://ccb.jhu.edu/software/fqtrim/
1. ^https://github.com/torognes/vsearch
References
AlZahal, O., Li, F., Guan, L. L., Walker, N. D., and McBride, B. W. (2017). Factors influencing ruminal bacterial community diversity and composition and microbial fibrolytic enzyme abundance in lactating dairy cows with a focus on the role of active dry yeast. J. Dairy Sci. 100, 4377–4393. doi: 10.3168/jds.2016-11473
AOAC. (2005). Official Methods of Analysis 16th Edn.. Washington, DC: Association of Analytical Chemists.
Beauchemin, K. A., Yang, W. Z., Morgavi, D. P., Ghorbani, G. R., Kautz, W., and Leedle, J. A. Z. (2003). Effects of bacterial direct-fed microbials and yeast on site and extent of digestion, blood chemistry, and subclinical ruminal acidosis in feedlot cattle. J. Anim. Sci. 81, 1628–1640. doi: 10.2527/2003.8161628x
Blümmel, M., and Ørskov, E. R. (1993). Comparison of in vitro gas production and nylon bag degradability of roughages in predicting feed intake in cattle. Anim. Feed Sci. Technol. 40, 109–119. doi: 10.1016/0377-8401(93)90150-I
Broderick, G. A., and Kang, J. H. (1980). Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media. J. Dairy Sci. 63, 64–75. doi: 10.3168/jds.S0022-0302(80)82888-8
Cai, L., Hartanto, R., Zhang, J., and Qi, D. (2021a). Clostridium butyricum improves rumen fermentation and growth performance of heat-stressed goats in vitro and in vivo. Animals 11:3261. doi: 10.3390/ani11113261
Cai, L., Yu, J., Hartanto, R., and Qi, D. (2021b). Dietary supplementation with Saccharomyces cerevisiae, clostridium butyricum and their combination ameliorate rumen fermentation and growth performance of heat-stressed goats. Animals 11:2116. doi: 10.3390/ani11072116
Cao, M., Li, Y., Wu, Q. J., Zhang, P., Li, W. T., Mao, Z. Y., et al. (2019). Effects of dietary clostridium butyricum addition to sows in late gestation and lactation on reproductive performance and intestinal microbiota. J. Anim. Sci. 97, 3426–3439. doi: 10.1093/jas/skz186
Chaucheyras-Durand, F., Walker, N. D., and Bach, A. (2008). Effects of active dry yeasts on the rumen microbial ecosystem: past. Present and future. Anim. Feed Sci. Technol. 145, 5–26. doi: 10.1016/j.anifeedsci.2007.04.019
Chiquett, J., Allison, M. J., and Rasmussen, M. (2012). Use of Prevotella bryantii 25A and a commercial probiotic during subacute acidosis challenge in midlactation dairy cows. J. Dairy Sci. 95, 5985–5995. doi: 10.3168/jds.2012-5511
Escobar, J. S., Klotz, B., Valdes, B. E., and Agudelo, G. M. (2014). The gut microbiota of Colombians differs from that of Americans. Eur. Asian. BMC Microbiol. 14, 1–14. doi: 10.1186/s12866-014-0311-6
Fan, Q., Wanapat, M., Yan, T., and Hou, F. (2020). Altitude influences microbial diversity and herbage fermentation in the rumen of yaks. BMC Microbiol. 20, 370–313. doi: 10.1186/s12866-020-02054-5
Guo, J., Li, P., Liu, S., Miao, B., Zeng, B., Jiang, Y., et al. (2020). Characterization of the rumen microbiota and volatile fatty acid profiles of weaned goat kids under shrub-grassland grazing and indoor feeding. Animals 10:176. doi: 10.3390/ani10020176
Hook, S. E., Steele, M. A., Northwood, K. S., Dijkstra, J., France, J., and Wright, A. D. (2011). Impact of subacute ruminal acidosis (SARA) adaptation and recovery on the density and diversity of bacteria in the rumen of dairy cows. FEMS Microbiol. Ecol. 78, 275–284. doi: 10.1111/j.1574-6941.2011.01154.x
Huang, C., Ge, F., Yao, X., Guo, X., Bao, P., Ma, X., et al. (2021). Microbiome and metabolomics reveal the effects of different feeding systems on the growth and ruminal development of yaks. Front. Microbiol. 12:682989. doi: 10.3389/fmicb.2021.682989
Huang, T., Peng, X. Y., Gao, B., Wei, Q. L., Xiang, R., Yuan, M. G., et al. (2019). The effect of clostridium butyricum on gut microbiota, immune response and intestinal barrier function during the development of necrotic enteritis in chickens. Front. Microbiol. 10:2309. doi: 10.3389/fmicb.2019.02309
Ítavo, C. C. B. F., Morais, M. G., Costa, C., Ítavo, L. C. V., Franco, G. L., da Silva, J. A., et al. (2011). Addition of propolis or monensin in the diet: behavior and productivity of lambs in feedlot. Anim. Feed Sci. Technol. 165, 161–166. doi: 10.1016/j.anifeedsci.2011.02.020
Jiao, P. X., He, Z. X., Ding, S., Walker, N. D., Cong, Y. Y., Liu, F. Z., et al. (2018). Impact of strain and dose of live yeast and yeast derivatives on in vitro ruminal fermentation of a high-grain diet at two pH levels. Can. J. Anim. Sci. 98, 477–487. doi: 10.1139/cjas-2017-0079
Jiao, P. X., Liu, F. Z., Beauchemin, K. A., and Yang, W. Z. (2017). Impact of strain and dose of lactic acid bacteria on in vitro ruminal fermentation with varying media pH levels and feedsubstrates. Anim. Feed Sci. Technol. 224, 1–13. doi: 10.1016/j.anifeedsci.2016.11.005
Jiao, P. X., Ma, F. C., Beauchemin, K. A., AlZaha, O., Xie, X. L., and Yang, W. Z. (2021). Effect of mixed live yeast and lactic acid bacteria on in vitro fermentation with varying media pH using a high-grain or high-forage diet. Can. J. Anim. Sci. 101, 370–380. doi: 10.1139/cjas-2020-0138
Jiao, P., Wei, C., Sun, Y., Xie, X., Zhang, Y., Wang, S., et al. (2019). Screening of live yeast and yeast derivatives for their impact of strain and dose on in vitro ruminal fermentation and microbial profiles with varying media pH levels in high-forage beef cattle diet. J. Sci. Food Agric. 99, 6751–6760. doi: 10.1002/jsfa.9957
Khafipour, E., Krause, D. O., and Plaizier, J. C. (2009). Alfalfa pellet-induced subacute ruminal acidosis in dairy cows increases bacterial endotoxin in the rumen without causing inflammation. J. Dairy Sci. 92, 1712–1724. doi: 10.3168/jds.2008-1656
Lana, R. P., Fox, D. G., Russell, J. B., and Perry, T. C. (1997). Influence of monensin on Holstein steers fed high-concentrate diets containing soybean meal or urea. J. Anim. Sci. 75, 2571–2579. doi: 10.2527/1997.75102571x
Li, Y., Wang, Y., Lv, J., Dou, X., and Zhang, Y. (2021). Effects of dietary supplementation with clostridium butyricum on the amelioration of growth performance, rumen fermentation, and rumen microbiota of Holstein heifers. Front. Nutr. 8:763700. doi: 10.3389/fnut.2021.763700
Li, Y., Wang, J., Mei, J., Huang, L., and Liu, H. (2022). Effects of mulberry branch and leaves silage on microbial community, rumen fermentation characteristics, and Milk yield in lactating dairy cows. Fermentation 8:86. doi: 10.3390/fermentation8020086
Li, H. Q., Zhou, Y., Ling, H. Y., Luo, L., Qi, D. S., and Feng, L. (2019). The effect of dietary supplementation with clostridium butyricum on the growth performance, immunity, intestinal microbiota and disease resistance of tilapia (Oreochromis niloticus). PLoS One 14:e0223428. doi: 10.1371/journal.pone.0223428
Liang, J., Kou, S., Chen, C., Raza, S. H. A., Wang, S., Ma, X., et al. (2021). Effects of clostridium butyricum on growth performance, metabonomics and intestinal microbial differences of weaned piglets. BMC Microbiol. 21, 85–16. doi: 10.1186/s12866-021-02143-z
Lima, F. S., Oikonomou, G., Lima, S. F., Bicalho, M. L. S., and Bicalho, R. C. (2014). Prepartum and postpartum rumen fluid microbiomes: characterization and correlation with production traits in dairy cows. Appl. Environ. Microbiol. 81, 1327–1337. doi: 10.1128/AEM.03138-14
Liu, Z., Wang, K., Nan, X., Cai, M., Yang, L., Xiong, B., et al. (2022). Synergistic effects of 3-nitrooxypropanol with fumarate in the regulation of propionate formation and methanogenesis in dairy cows in vitro. Appl. Environ. Microbiol. 88:e01908-2. doi: 10.1128/aem.01908-21
Ma, Y., Wang, C., Zhang, H., Yu, L., Dong, L., Gong, D., et al. (2021). Illumina sequencing and metabolomics analysis reveal thiamine modulation of ruminal microbiota and metabolome characteristics in goats fed a high-concentrate diet. Front. Microbiol. 12:653283. doi: 10.3389/fmicb.2021.653283
Mani, S., Aiyegoro, O. A., and Adeleke, M. A. (2021). Characterization of rumen microbiota of two sheep breeds supplemented with direct-fed lactic acid bacteria. Front. Vet. Sci. 7:570074. doi: 10.3389/fvets.2020.570074
McAllister, T. A., Beauchemin, K. A., Alazzeh, A. Y., Baah, J., Teather, R. M., and Stanford, K. (2011). Review: the use of direct fed microbials to mitigate pathogensand enhance production in cattle. Can. J. Anim. Sci. 91, 193–211. doi: 10.4141/cjas10047
Nagaraja, T. G., and Titgemeyer, E. C. (2007). Ruminal acidosis in beef cattle: the current microbiological and nutritional outlook. J. Dairy Sci. 90, E17–E38. doi: 10.3168/jds.2006-478
National Research Council (NRC) (2001). Nutrient Requirements of Dairy Cattle. 7th Edn. Washington, DC: Academies Press.
Palmonari, A., Stevenson, D. M., Mertens, D. R., Cruywagen, C. W., and Weimer, P. J. (2010). pH dynamics and bacterial community composition in the rumen of lactating dairy cows. J. Dairy Sci. 93, 279–287. doi: 10.3168/jds.2009-2207
Qiu, Q., Gao, C., Aziz ur Rahman, M., Cao, B., and Su, H. (2020). Digestive ability, physiological characteristics, and rumen bacterial community of Holstein finishing steers in response to three nutrient density diets as fattening phases advanced. Microorganisms 8:335. doi: 10.3390/microorganisms8030335
Ran, T., Shen, Y., Saleem, A. M., AlZahal, O., Beauchemin, K. A., and Yang, W. (2018). Using ruminally protected and nonprotected active dried yeast as alternatives to antibiotics in finishing beef steers: growth performance, carcass traits, blood metabolites, and fecal Escherichia coli. J. Anim. Sci. 96, 4385–4397. doi: 10.1093/jas/sky272
Romero-Perez, A., and Beauchemin, K. (2018). Estimating gas volume from headspace pressure in a batch culture system. Can. J. Anim. Sci. 98, 593–596. doi: 10.1139/cjas-2017-0100
Roselli, M., Pieper, R., Rogel-Gaillard, C., de Vries, H., Bailey, M., Smidt, H., et al. (2017). Immunomodulating effects of probiotics for microbiota modulation, gut health and disease resistance in pigs. Anim. Feed Sci. Technol. 233, 104–119. doi: 10.1016/j.anifeedsci.2017.07.011
Russell, J. B., and Wilson, D. B. (1996). Why are ruminal cellulolytic bacteria unable to digest cellulose at low pH? J. Dairy Sci. 79, 1503–1509. doi: 10.3168/jds.S0022-0302(96)76510-4
Tong, F., Wang, T., Gao, N. L., Liu, Z., Cui, K., Duan, Y., et al. (2022). The microbiome of the buffalo digestive tract. Nat. Commun. 13, 1–16. doi: 10.1038/s41467-022-28402-9
Van Soest, P. J., Robertson, J. B., and Lewis, B. A. (1991). Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 74, 3583–3597. doi: 10.3168/jds.S0022-0302(91)78551-2
Wang, W. W., Wang, J., Zhang, H. J., Wu, S. G., and Qi, G. H. (2020). Effects of clostridium butyricum on production performance and intestinal absorption function of laying hens in the late phase of production. Anim. Feed Sci. Technol. 264:114476. doi: 10.1016/j.anifeedsci.2020.114476
Wei, X., Ouyang, K., Long, T., Liu, Z., Li, Y., and Qiu, Q. (2022). Dynamic variations in rumen fermentation characteristics and bacterial community composition during in vitro fermentation. Fermentation 8:276. doi: 10.3390/fermentation8060276
Yu, S., Zhang, G., Liu, Z., Wu, P., Yu, Z., and Wang, J. (2020). Repeated inoculation with fresh rumen fluid before or during weaning modulates the microbiota composition and cooccurrence of the rumen and colon of lambs. BMC Microbiol. 20, 1–15. doi: 10.1186/s12866-020-1716-z
Zhang, J., Shi, H., Yang, Y., Li, S., Cao, Z., Ji, S., et al. (2017). Effect of dietary forage to concentrate ratios on dynamic profile changes and interactions of ruminal microbiota and metabolites in Holstein heifers. Front. Microbiol. 8:2206. doi: 10.3389/fmicb.2017.02206
Zhang, C., Yu, Q., Wang, J., Yu, Y., Zhang, T., and Sun, Y. (2022). Effects of dietary supplementation with clostridium butyricum on growth performance, apparent digestibility, blood metabolites, ruminal fermentation and bacterial communities of fattening goats. Front. Nutr. 9:888191. doi: 10.3389/fnut.2022.888191
Keywords: Clostridium butyricum, batch culture, media pH level, microbiota, rumen fermentation
Citation: Zhang M, Liang G, Zhang X, Lu X, Li S, Wang X, Yang W, Yuan Y and Jiao P (2022) The gas production, ruminal fermentation parameters, and microbiota in response to Clostridium butyricum supplementation on in vitro varying with media pH levels. Front. Microbiol. 13:960623. doi: 10.3389/fmicb.2022.960623
Edited by:
Junjun Wang, China Agricultural University, ChinaReviewed by:
Zhaohai Wu, Institute of Animal Sciences (CAAS), ChinaYanfen Cheng, Nanjing Agricultural University, China
Copyright © 2022 Zhang, Liang, Zhang, Lu, Li, Wang, Yang, Yuan and Jiao. 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: Peixin Jiao, peixin.jiao@neau.edu.cn; Yuan Yuan, 006707@yzu.edu.cn
†These authors have contributed equally to this work