Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 05 November 2021
Sec. Microbiotechnology

Changes in Carbohydrate Composition in Fermented Total Mixed Ration and Its Effects on in vitro Methane Production and Microbiome

\r\nYang LiYang Li1Jingyi LvJingyi Lv1Jihong WangJihong Wang1Shuang ZhouShuang Zhou1Guangning ZhangGuangning Zhang1Bingdong WeiBingdong Wei2Yukun SunYukun Sun1Yaxue LanYaxue Lan1Xiujing Dou*Xiujing Dou1*Yonggen Zhang*Yonggen Zhang1*
  • 1College of Animal Sciences and Technology, Northeast Agriculture University, Harbin, China
  • 2Jilin Academy of Agricultural Sciences, Changchun, China

The purpose of this experiment was to investigate the changes of carbohydrate composition in fermented total mixed diet and its effects on rumen fermentation, methane production, and rumen microbiome in vitro. The concentrate-to-forage ratio of the total mixed ration (TMR) was 4:6, and TMR was ensiled with lactic acid bacteria and fibrolytic enzymes. The results showed that different TMRs had different carbohydrate compositions and subfractions, fermentation characteristics, and bacterial community diversity. After fermentation, the fermented total mixed ration (FTMR) group had lower contents of neutral detergent fiber, acid detergent fiber, starch, non-fibrous carbohydrates, and carbohydrates. In addition, lactic acid content and relative abundance of Lactobacillus in the FTMR group were higher. Compared with the TMR group, the in vitro ammonia nitrogen and total volatile fatty acid concentrations and the molar proportion of propionate and butyrate were increased in the FTMR group. However, the ruminal pH, molar proportion of acetate, and methane production were significantly decreased in the FTMR group. Notably, we found that the relative abundance of ruminal bacteria was higher in FTMR than in TMR samples, including Prevotella, Coprococcus, and Oscillospira. At the same time, we found that the diversity of methanogens in the FTMR group was lower than that in the TMR group. The relative abundance of Methanobrevibacter significantly decreased, while the relative abundances of Methanoplanus and vadinCA11 increased. The relative abundances of Entodinium and Pichia significantly decreased in the FTMR group compared with the TMR group. These results suggest that FTMR can be used as an environmentally cleaner technology in animal farming due to its ability to improve ruminal fermentation, modulate the rumen microbiome, and reduce methane emissions.

Introduction

Methane is a greenhouse gas with an abundance second only to carbon dioxide, and the greenhouse efficiency is 25 times that of the equivalent amount of carbon dioxide (IPCC, 2007), which can enhance radiative forcing and greenhouse effects, and further aggravate disasters such as climate deterioration, global warming, and land desertification. Methanogenesis is a fundamental rumen metabolic process; it is also the main cause of nutrient loss and energy waste in feed (Takahashi et al., 1997; Takahashi, 2001) and contributes to 11–17% of the global greenhouse gas emissions (Beauchemin et al., 2009; Goel and Makkar, 2012). Therefore, whether from the perspective of environmental protection or animal production, it is extremely necessary to reduce methane emissions from ruminants.

Methane reduction is an important research topic not only in the field of ruminant nutrition but also in environmental protection. In animal production, methane emissions can be reduced through dietary manipulation and feed additives. However, the ideal methane inhibitor needs to have many characteristics, such as being extremely specific and long-lasting, safe for use in animals, and not leaving behind residues in edible products (Van Nevel and Demeyer, 1996). In addition, ionophore antibiotics (monensin) have been banned in the production of heifers and dairy cows in some countries. Therefore, it is a safer and more effective method to reduce methane production by using ruminant diet nutrition control methods.

Fermented total mixed ration (FTMR) is a new type of ruminant feeding technology, in which the finished total mixed rations (TMR) are tightly wrapped by a professional baler and a special plastic stretch film, and the TMR is sealed by a specific single or compound strain (Liu et al., 2016a,b) for anaerobic fermentation so that the TMR can be stored for a certain time. Our previous studies proved that FTMR could improve the feed efficiency and lactation performance of dairy cows (Zhang et al., 2020b); we also discovered that the contents of structural carbohydrates (NDF and ADF) in FTMR were significantly reduced with an increased lactic acid content (Zhang et al., 2020a,b). Diet composition can have a major effect on methane production (Benchaar et al., 2001). Highly concentrated diets should reduce methane production (Granja-Salcedo et al., 2016). The type of carbohydrate in rumen fermentation regulates methane production mainly by affecting the rumen pH value and the total amount of rumen microbiome. Reducing the proportion of structural carbohydrates in the diet reduces rumen pH and thus reduces methane production by inhibiting the activity of methanogens (Russell, 1998). Changes in fermentation substrates in TMR also affect rumen fermentation types and volatile fatty acid synthesis pathways (Cao et al., 2010). The inhibition of methane production is normally accompanied by increased in propionate production (Wolin, 1975), which uses hydrogen and lactic acid. Furthermore, the relative abundance of rumen bacteria and protozoa, especially methanogens, is an important reason to explain the variation in ruminal methane production. It is important to fully understand the microbial changes involved in rumen methane production to clarify the mechanism by which FTMR regulates it. Therefore, we hypothesized that FTMR might decrease methane production by altering the rumen fermentation pattern and inhibiting the major methanogens or bacteria and protozoa related to methanogenesis in the rumen.

So the current study aimed to evaluate the effects of fermentation on the carbohydrate composition and fermentation index of TMR and investigate the effects of FTMR on in vitro ruminal fermentation, methane production, and rumen microbiome. The results from this study are expected to improve the carbohydrate composition of TMR concomitant with decreased CH4 production for cleaner environmental animal agriculture.

Materials and Methods

Preparation of Total Mixed Ration

Total mixed ration in this study was prepared using compound feed in the Ruminant Nutrition Laboratory from Northeast Agriculture University. Corn silage, alfalfa hay, steam-flaked corn, soybean meal, rapeseed meal, dry corn gluten feed, distillers dried grains with soluble (DDGS), and premix were procured from the Songhuajiang Dairy Farm (Harbin city, Heilongjiang province, China). The ingredient and chemical composition are shown in Table 1. The TMR was formulated with a forage-to-concentrate ratio of 60:40 (DM basis) and met the animals’ requirements for nutrition based on the Cornell-Penn-Miner dairy model (Version 3.0.10; Cornell University, Ithaca, NY; University of Pennsylvania, Kennett Square, PA; and William H. Miner Agricultural Research Institute, Chazy, NY). Alfalfa hay was chopped to a length of 3–4 cm.

TABLE 1
www.frontiersin.org

Table 1. Ingredient and nutrient composition of the experimental ration (DM basis).

The TMR was ensiled with a combination of lactic acid bacteria [LAB; Lactobacillus plantarum CGMCC10516 and Lactobacillus buchneri BNCC189797, applied at a ratio of 1:1; theoretical final application rate of 106 colony-forming units (cfu)/g of fresh matter (FM)] and fibrolytic enzyme (1 g fibrolytic enzyme per 1 kg FM [EN; 10,000 U/g activity, XS Biotechnology Co., Ltd., Beijing, China)]. Additives were homogenously mixed into TMR using a hand sprayer. After proper mixing, the TMR was packed into a 20 × 30 cm plastic laboratory fermentation bag (Chuangjia Packaging Material Co., Ltd., Wenzhou, China). The air was subsequently removed by a vacuum packing machine (Maige Automation Equipment Co., Ltd., Qingdao, China). Eight bags (1 ensiling day × 1 treatment × 8 replicates) were prepared and incubated indoor at ambient temperature (20–30°C) for 30 days. Three bags (3 replicates) of non-fermented total mixed and FTMR were opened to analyze them for fermentation characteristics, chemical compositions, and in vitro fermentation. Five samples (5 replicates) per two treatments were analyzed to determine the microbial community.

Fermentation Index, Carbohydrate Composition, and Subfractions of Carbohydrates Using Cornell Net Carbohydrate and Protein System Analysis

After 30 days, 50 g of the wet ensiled silage samples of FTMR and the non-fermented TMR were homogenized with 200 mL of sterilized distilled water and stored at 4°C overnight (Cai et al., 1999). The pH was measured using a pH meter (Sartorius basic pH meter, Göttingen, Germany). The concentrations of lactic acid (LA), acetic acid (AA), and butyric acid (BA) were measured by high-performance liquid chromatography (Yuan et al., 2016), while ammonia nitrogen (NH3-N) was determined according to the phenol/hypochlorite method (Broderick and Kang, 1980).

Fermented or non-fermented TMR were dried at 60°C for 48 h and ground with a high-speed universal mill to pass through a 1-mm sieve grind. The dry matter (DM, AOAC 930.15) content of fermented or non-fermented TMR samples was determined by the methods of AOAC (2005). The contents of neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), and residual ash were measured according to the filtration method (Van Soest et al., 1991) using an Ankom 220 fiber analyzer (Ankom Technology Corp., Macedon, NY, United States), and α-amylase and sodium sulfite were used for the NDF procedure. Starch content was determined using the Megazyme Total Starch Assay Procedure (product no: K-TSTA; Megazyme International Ireland Ltd., Wicklow, Ireland). All analyses were conducted in triplicate. The lignin content was calculated as ADL—residual ash, non-fibrous carbohydrate (NFC) as OM—CP—NDF—fat; and carbohydrate (CHO) as 100—EE—CP—ash. The carbohydrate subfractions were partitioned into CA, CB1, CB2, and CC according to the Cornell Net Carbohydrate and Protein System (CNCPS) (Sniffen et al., 1992).

In vitro Incubation

The in vitro incubation procedures were as previously described (Li et al., 2019). The substrates were prepared with fermented or non-fermented TMR by drying and grinding through a 0.45 mm sieve. Rumen fluid was collected via rumen cannula before the morning feeding from three steers of the yellow cattle breed, filtered with four-layer cheesecloth, and mixed with preheated artificial saliva (Menke and Steingass, 1988) at a ratio of 2:1 (buffer: ruminal fluid, v:v). The ruminal fluid (150 mL) that was buffered was dispensed into prewarmed 200-mL incubation flasks. Two grams of each substrate was blended with the buffered ruminal fluid in each incubation flask, which was incubated at 39°C for 24 h in water. The methane production was measured by a real-time in vitro fermentation system (produced by Jilin Academy of Agricultural Sciences, code Qtfxy-6), which was tested for the effluent gas discharged from each incubation flask. Nitrogen (purity 99.99%) was passed into the incubation flask from the bottom at a speed of 200 mL/min. Methane was carried by nitrogen into an AGM10 sensor (Sensors Europe GmbH, Erkrath, FRG), and the concentration of methane was measured and recorded every 6 min (Sun et al., 2017). After determining methane production, the pH of the culture liquor was determined using a pH meter (Sartorius basic pH meter, Göttingen, Germany). For analysis of NH3-N and volatile fatty acids (VFAs), 1 ml of 25% meta-phosphoric acid was added to 5 ml of culture liquor and stored at −20°C until analysis. The VFA concentration was determined by gas chromatography as previously described (Stewart and Duncan, 1985). Ammonia nitrogen was determined according to the phenol/hypochlorite method (Broderick and Kang, 1980). The microbial crude protein (MCP) was determined following the procedures of Hu et al. (2005). All analyses were conducted in triplicate.

Microbial Community Diversity Analysis

The remaining silage subsample (10 g) from fermented or non-fermented TMR was mixed with 40 mL of saline solution (NaCl, 0.90 g/g) and shaken at 120 r/m for 2 h. The filtered liquor through gauze was centrifuged at 10,000 r/m for 10 min at 4°C. The supernatant was discarded, but the deposit was suspended in 3 mL of saline solution (Zhang et al., 2020a). This procedure was repeated for good precipitation. According to the manufacturer’s protocol, genomic DNA from fermented or non-fermented TMR was extracted using Fast DNA SPIN extraction kits (MP Biomedicals, Santa Ana, CA, United States). Metagenomic DNA was extracted from ruminal samples using the cetyltrimethylammonium bromide (CTAB) method but with bead beating (Wang et al., 2019). The quality of the DNA extracts was evaluated using agarose (1%) electrophoresis, while the DNA concentration was determined using the Qubit dsDNA BR Assay Kit (Invitrogen Corporation, United States) on a Qubit 2.0 fluorometer (Invitrogen Corporation, United States). The extracted DNA was subjected to PCR amplification in triplicate using the Accuprime Taq DNA Polymerase System (Invitrogen, Carlsbad, CA, United States). For bacterial analysis, the V3-V4 region of the 16S rRNA gene was amplified using primers 338F (5′-ACTCCTRCGGGAGGCAGCAG-3′) and 806R (5′-GGACTACCVGGGTATCTAAT-3′) (Mao et al., 2015). For ciliate protozoal analysis, V3-V4 and signature regions 1–2 of the 18S rRNA gene were amplified using primers PSSU-316F (5′-GCTTTCGWTGGTAGTGTATT-3′) (Sylvester et al., 2004) and GIC758R (5′-CAACTGTCTCTATKAAYCG-3′) (Ishaq and Wright, 2014). For archaeal analysis, the V3-V4 region of the 16S rRNA gene was amplified using primers F341 (5′-CTACGGGGYGCASCAG-3′) (Wei et al., 2013) and R806 (5′-GGACTACVVGGGTATCTAATC-3′). After purification and quantification of the PCR products, the DNA fragments of the community were sequenced by double-terminal (paired-end) sequencing using an Illumina MiSeq platform (Wuhan Frasergen Bioinformatics Co., Ltd., Wuhan, China). Sequencing data were processed using the Quantitative Insights Into Microbial Ecology (QIIME, v1.8.0)1 pipeline as previously described (Caporaso et al., 2010). Alpha diversity measurements, including Chao1 richness estimates, Shannon and Simpson diversity indices, and observed species, were calculated for each sample. The microbiome was compared as beta diversity using the distance matrices generated from weighted UniFrac and principal coordinates analysis (PCoA). Taxa abundance at the genus level was statistically compared between the samples for fermented or non-fermented TMR and their abundance of rumen microbiome.

Statistical Analyses

Data were analyzed using SAS software (version 9.4, SAS Institute Inc., Cary, NC). Data for rumen microbiome were subjected to the non-parametric Wilcoxon Two-Sample Test to examine the relative abundance of the genus-level bacterial community in TMR or FTMR and the effect of TMR after fermentation. The other data were analyzed using the one-way ANOVA procedure. Significance was declared at P < 0.05 and trends at 0.05 ≤ P < 0.10. Standard errors of the mean and standard deviation (rumen microbiome) are reported.

Results and Discussion

Compared with the control group, the FTMR group had lower DM (P = 0.0009), NDF (P = 0.006), ADF (P = 0.001), starch (P = 0.0002), NFC (P < 0.0001), CHO (P = 0.0006) contents, and CB1 (P = 0.0003) and CB2 (P = 0.03) subfractions. Moreover, the fractions of CA1 (P = 0.002) and CC (P = 0.01) in TMR increased significantly after fermentation (Table 2). This indicates that fermentation improves the availability of carbohydrates in TMR. It can be seen from the results of this experiment that the changes in carbohydrate composition of TMR by cofermentation of lactic acid bacteria and cellulose enzymes are obvious, which also lay a foundation for the change in rumen fermentation type and the reduction in methane production. A high-concentrate diet should decrease methane emissions (Granja-Salcedo et al., 2016); therefore, a diet with a higher proportion of roughage was designed in this study to emphasize the role of FTMR in decreasing methane emissions. The low DM content in fermentation TMR may be due to the decomposition of nutrients into liquids and gases by microorganisms (Kim et al., 2016). Cao et al. (2010) found that the contents of DM, NFC, NDF, and ADF in FTMR were lower than those in TMR, which was consistent with the results of this study. The contents of starch, NDF, and ADF in FTMR were low, which might be caused by lactic acid bacteria hydrolyzing digestible carbohydrates in the feed to produce lactic acid. In addition, the degradation of plant cell walls by cellulase to increase the content of water-soluble carbohydrates could also be a major reason for the decrease in NDF and ADF (Liu et al., 2016a). The digestible carbohydrate components in TMR will be used as substrates for lactic acid production to facilitate the fermentation of the diet (Xing et al., 2009). The carbohydrate components of CNCPS could reflect the degradability of dietary carbohydrates in the rumen. The increase in the CA fraction and decrease in the fractions of CB1 and CB2 in FTMR are related to the conversion of starch and digestible fiber in feed to lactate by lactobacillus. The lower content of NDF and ADF in the FTMR diet may be due to the fact that lactic acid bacteria hydrolyzed more digestible plant cells to produce acid during silage (Liu et al., 2016a). Moreover, inoculated LAB have higher fermentation efficiency than epiphytic LAB to transform sugars into lactic acid (Yuan et al., 2015). In this study, a fibrolytic enzyme was added to TMR to degrade the plant cell wall and produce lactate as a substrate (Xing et al., 2009); however, DM loss due to fermentation increased the CC fraction.

TABLE 2
www.frontiersin.org

Table 2. Chemical profiles and carbohydrate subfractions in fermented or non-fermented TMR.

As shown in Table 3, fermentation reduced the pH value of TMR to below 4.2. Compared with TMR, ammonia nitrogen, lactic, and acetic acid concentration in FTMR was significantly increased. Butyric acid was detected in both TMR groups, indicating that the quality of both TMR groups was good. The rapid decline in pH depends on lactic acid produced by lactic acid bacteria and cellulases degrading carbohydrates, which is important in reducing nutrient loss caused by the growth of undesirable microorganisms during the early stage of fermentation (Liu et al., 2016a). When the pH value of silage is lower than 4.2, the quality of silage can be guaranteed (Wen et al., 2017). Large amounts of starch and fermentable carbohydrates in TMR provide sufficient fermentation substrates for lactic acid bacteria and cellulase, leading to a significant increase in lactic acid content in FTMR. Previous studies have found that LAB containing Lactobacillus plantarum and Lactobacillus buchneri inoculated in TMR produced a large amount of lactic and acetic acid during fermentation (Zhang et al., 2020b). Liu et al. (2016a) also found that the addition of microbial preparations and enzyme preparations in FTMR could increase the fermentation of lactic acid bacteria. With the extension of fermentation time, lactic acid is converted into acetic acid, which is consistent with the results of Alli et al. (2010). Although the pH of TMR was significantly decreased, we did not find any negative effects of FTMR on the performance of dairy cows in previous experiments (Zhang et al., 2020b). The ammonia-nitrogen content in FTMR was significantly increased, which proved that there was more decomposition of true protein in the feed after fermentation (Xing et al., 2009), but the increase in non-protein nitrogen content also provided an abundant nitrogen source for the rumen.

TABLE 3
www.frontiersin.org

Table 3. Fermentation characteristics of fermented or non-fermented TMR.

Table 4 and Figure 1 show that lactic acid bacteria became the absolute dominant bacterial genus in the TMR after fermentation (P < 0.0001). The PCoA plots (Figure 2) clearly indicate the microbial community variability and show that the distance between the FTMR and the TMR was far. Lactobacillus and Acetobacter were the two main bacterial genera that changed after fermentation. In this experiment, the addition of exogenous Lactobacillus increased the relative abundance of lactic acid bacteria in the TMR. Moreover, the addition of cellulase degraded the plant cell wall in the diet and increased the relative abundance of Lactobacillus (Xing et al., 2009). In the non-fermented TMR, it may be that the Acetobacter was attached to the raw materials of the TMR, making Lactobacillus the dominant bacteria in mixing and stirring. In conclusion, Lactobacillus and cellulase catalyze the fermentation process, promote the growth of Lactobacillus, and increased the content of lactic acid and acetic acid, thereby reducing the pH value (pH < 4.2), inhibiting the growth of other bacteria (Acetobacter, Gluconacetobacter, and Paenibacillus), and making Lactobacillus the dominant bacteria.

TABLE 4
www.frontiersin.org

Table 4. Relative abundance (%) of 10 most predominant genera in fermented or non-fermented TMR.

FIGURE 1
www.frontiersin.org

Figure 1. Relative abundance of bacterial under genera level in fermented or non-fermented TMR.

FIGURE 2
www.frontiersin.org

Figure 2. Principal Co-ordinates analysis of the bacterial community in fermented or non-fermented TMR.

In the in vitro gas production test, we further found that FTMR increased ammonia nitrogen (P = 0.003), total volatile fatty acid concentration (P = 0.0048), and molar proportion of propionate (P = 0.03) and butyrate (P = 0.001) and decreased rumen pH (P = 0.03), molar proportion of acetate (P = 0.006), acetate to propionate ratio (P = 0.01) and methane production (P = 0.0007) (Table 5). Figure 3 shows the real-time methane production in vitro with fermented or non-fermented TMR. The higher non-protein nitrogen content in FTMR may lead to an increase in ammonia-nitrogen content in the rumen, and this result was consistent with that of Cao et al. (2010). In this experiment, it was found that the content of NDF and ADF in TMR decreased, and the acetic acid yield decreased after fermentation. More structural carbohydrates in FTMR are converted into rapidly degradable carbohydrates. However, the increase in the CA fraction in FTMR may increase volatile fatty acid concentration and a decrease in pH in vitro. There was a negative correlation between rumen pH and fermentation gas production (Waghorn, 1991). Reducing the proportion of structural carbohydrates in the diet decreases rumen pH and thus reduces methane production by inhibiting the activity of methanogens (Russell, 1998). Bacteria degraded the feed concentrate and CA in TMR to produce propionic acid, which decreased rumen pH and the acetate to propionate ratio. At the same time, lactic acid from FTMR and rumen fermentation may have further increased the amount of propionic acid using the hydrogen produced by the fermentation reaction (Cao et al., 2010). The inhibition of methane production by FTMR may be due to indirect or direct inhibition (or both) of methanogens via a decline in H2 production due to reduced acetate and more propionate production (Cieslak et al., 2013; Kumar et al., 2014). In addition, the high content of lactic acid in FTMR is also key to regulating methane production. Propionic acid was produced when lactic acid in the feed was secondarily fermented in the rumen by lactate-utilizing bacteria (Megasphaera elsdenii, Selenomonas ruminantium, and Veillonella parvula) (Dawson et al., 1997; Russell and Wallace, 1997). Propionic acid fermentation uses electrons to reduce methane production. Moreover, lactic acid in the rumen uses hydrogen for conversion to propionic acid, and the reduction of hydrogen also inhibits the formation of methane from hydrogen and carbon dioxide (Moss et al., 2000). In addition, shifting ruminal fermentation to more propionate would inhibit hydrogen-producing bacteria (Ruminococcus albus, Ruminococcus flavefaciens, and protozoa) (Cobellis et al., 2016), so it is necessary to understand further the effect of FTMR on the rumen microbiome in vitro.

TABLE 5
www.frontiersin.org

Table 5. In vitro pH, NH3-N, MCP, VFA, and methane production of fermented or non-fermented TMR.

FIGURE 3
www.frontiersin.org

Figure 3. Methane production curve in vitro with fermented or non-fermented TMR.

As seen in Table 6, FTMR improved the Chao1 (P = 0.001) and observed species (P = 0.002) of rumen bacteria and tended to reduce the Shannon index of methanogens (P = 0.07), which proved that FTMR had higher ruminal bacterial abundance and lower methanogen diversity. In addition, we found that FTMR significantly increased the relative abundance of Prevotella (P < 0.05), Coprococcus (P < 0.05) and Oscillospira (P < 0.05); and significantly reduced the relative abundance of Entodinium (P < 0.05) and Pichia (P < 0.05). Although it increased the relative abundance of Methanoplanus (P < 0.05) and vadinCA11 (P < 0.05), it had a significant reduction effect on Methanobrevibacter (P < 0.05) (Table 7). In an experiment to regulate methane production, bacterial genera were analyzed to understand the effects of exogenous substances on the rumen microbiome (Martínez-Fernández et al., 2014; Zhou et al., 2017, 2020). In the present study, 16S rRNA sequencing was used to comparatively examine the influence of FTMR on ruminal bacteria, protozoa, and methanogens. Prevotella, a gram-negative genus of Bacteroidetes, increased with FTMR supplementation, possibly due to increased fermentation substrates such as non-protein nitrogen and reduced competition from other bacteria. It has been suggested that FTMR could increase ruminal protein degradation and ammonia concentrations (Zhang et al., 2020b). We also observed that FTMR significantly increased ammonia nitrogen, which is in agreement with a previous study (Cao et al., 2010). Studies have shown that Coprococcus is involved in an important metabolic pathway in the rumen (Martinez-Fernandez et al., 2017). In cows with higher feeding efficiency, the abundance of genes related to the acrylate pathway has increased, which involves the conversion of lactic acid to propionic acid (Martinez-Fernandez et al., 2017). In general, the main pathway for propionic acid synthesis in highly efficient animals is the acrylate pathway rather than the succinate pathway (Shabat et al., 2016). Shabat et al. (2016) found that the relative abundance of Coprococcus in the rumen of cows with low CH4 emissions and high feed efficiency was higher, consistent with the present experimental results. Furthermore, Mackie et al. (2003) found that the abundance of Oscillospira, which is the bacterium involved in the degradation of the plant cell wall (Yanagita et al., 2003), in the rumen was diet-dependent, and that Oscillospira may play an important role in producing or utilizing amino acids (Hua et al., 2017). Therefore, we think that the change in the relative abundance of Oscillospira was related to the change in ammonia nitrogen and structural carbohydrate content in FTMR. The number of rumen protozoa is closely related to the number of methanogens and methane production, and Entodinium had the highest contribution rate to methane production (Ranilla et al., 2007). Eugène et al. (2004) conducted a meta-analysis of data in the published literature and found that defaunation could increase the concentration of propionic acid, and the increase in propionic acid production would reduce methane production through the utilization of hydrogen. Moreover, almost all rumen ciliates have methanogens on the surface. It has also been found that there are endosymbiotic methanogens in the cytosol of rumen ciliates, and the number of methanogens in the cytosol is far greater than the number of methanogens bound to the surface of the ciliates (Finlay et al., 1994). Thirty-seven percent of methane in ruminants is produced by protozoan-related methane bacteria (Finlay et al., 1994). In this experiment, the relative abundance of Entodinium and Pichia decreased, and the rumen methanogens decreased, thereby reducing rumen methane production. Although changes in the abundance of rumen bacteria and protozoa have a positive impact on in vitro methane production, studies have indicated that ruminal methane production may be much more influenced by the relative abundance of archaea rather than by the microbial population structure (Wallace, 2004; Duarte et al., 2017). Methanobacterium, Methanobrevibacter, Methanomicrobium, and Methanosarcine are four methanogens widely existing in the rumen of ruminants, among which Methanobrevibacter is the most dominant genus (Jarvis et al., 2000). In this study, FTMR significantly reduced the relative abundance of Methanobrevibacter in the rumen, which may be due to the reduction in the abundance of protozoa, resulting in the lack of attachments, which in turn inhibited the activity of methanogens and reduced methane production. It is also possible that FTMR decreased the rumen pH and inhibited the activity of methanogens in the rumen, thus reducing methane production. Interestingly, we found that FTMR increased the relative abundance of Methanoplanus and Vadinca11, which may be caused by the decrease in the abundance of other methanogens. However, the causes of the increase in these genera are not well understood, and further studies are needed to elucidate their mechanism.

TABLE 6
www.frontiersin.org

Table 6. Changes in richness and diversity of rumen microbiome in fermenters fed fermented or non-fermented TMR.

TABLE 7
www.frontiersin.org

Table 7. Relative abundance (%) of 10 most predominant genera of rumen microbiome in fermenters fed fermented or non-fermented TMR.

Conclusion

Compared to non-fermented TMR, FTMR has a lower structural carbohydrate content and a higher CA fraction, lactic acid content and relative abundance of Lactobacillus. FTMR increases the molar proportion of propionate and butyrate, decreases ruminal pH and the ratio of acetate to propionate and lowers methane emissions. The effect of FTMR on reducing methane emissions seems to result from the conversion of lactic acid to propionic acid and the decrease in protozoa and methanogens in the rumen. FTMR could be used as an environmentally cleaner technology in animal farming due to its ability to improve ruminal fermentation of diets and reduce methane production.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The animal study was reviewed and approved by all animal studies were conducted according to the animal care and use guidelines of the Animal Care and Use Committee of Animal Science and Technology College, Northeast Agricultural University.

Author Contributions

YLi, XD, and YZ designed the experiment. JL, JW, and SZ conducted the experiment and collected samples. YS and YLa analyzed the samples and data. YLi wrote the manuscript. YLi, JL, GZ, and BW revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This study was financially supported by the China Postdoctoral Science Foundation Grant for the 65th batch (2019M651252), the Young Talents Project of the Northeast Agricultural University (18QC35), and the Shuangcheng Nestle Dairy Farming Institute (NEAU-JKYZYJY-2020-1).

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.

Abbreviations

AA, acetic acid; ADL, acid detergent lignin; ADF, acid detergent fiber; BA, butyric acid; CH4, methane; cfu, colony-forming unit; CP, crude protein; CHO, carbohydrate; CNCPS, Cornell net carbohydrate and protein system; DM, dry matter; DNA, deoxyribonucleic acid; DDGS, dried distiller grains with solubles; EE, ether extract; EN, enzyme; FM, fresh matter; FTMR, fermented total mixed ration; GLM, general linear model; LA, lactic acid; LAB, lactic acid bacteria; NDF, neutral detergent fiber; NFC, non-fibrous carbohydrate; NH3-N, ammonia nitrogen; OM, organic matter; PCoA, principal coordinates analysis; TMR, total mixed ration; VFA, volatile fatty acid.

Footnotes

  1. ^ http://qiime.org/

References

Alli, I., Fairbairn, R., Noroozi, E., and Baker, B. E. (2010). The effects of molasses on the fermentation of chopped whole−plant leucaena. J. Sci. Food Agric. 35, 285–289. doi: 10.1002/jsfa.2740350307

CrossRef Full Text | Google Scholar

AOAC. (2005). Official Methods of Analysis, sixteenth Edn. Arlington, VA: Association of Official Analytical Chemists.

Google Scholar

Beauchemin, K. A., McAllister, T. A., and McGinn, S. M. (2009). Dietary mitigation of enteric methane from cattle. CAB Rev. 35, 111–118. doi: 10.1079/PAVSNNR20094035

CrossRef Full Text | Google Scholar

Benchaar, C., Pomar, C., and Chiquette, J. (2001). Evaluation of dietary strategies to reduce methane production in ruminants: a modelling approach. Can. J. Anim. Sci. 81, 563–574. doi: 10.4141/A00-119

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Cai, Y., Benno, Y., Ogawa, M., and Kumai, S. (1999). Effect of applying lactic acid bacteria isolated from forage crops on fermentation characteristics and aerobic deterioration of silage. J. Dairy Sci. 82, 520–526. doi: 10.3168/jds.S0022-0302(99)75263-X

CrossRef Full Text | Google Scholar

Cao, Y., Takahashi, T., Horiguchi, K. I., Yoshida, N., and Cai, Y. M. (2010). Methane emissions from sheep fed fermented or non-fermented total mixed ration containing whole-crop rice and rice bran. Anim. Feed Sci. Technol. 157, 72–78. doi: 10.1016/j.anifeedsci.2010.02.004

CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Cieslak, A., Szumacher-Strabel, M., Stochmal, A., and Oleszek, W. (2013). Plant components with specific activities against rumen methanogens. Animal 7, 253–265. doi: 10.1017/S1751731113000852

PubMed Abstract | CrossRef Full Text | Google Scholar

Cobellis, G., Trabalza-Marinucci, M., Marcotullio, M. C., and Yu, Z. (2016). Evaluation of different essential oils in modulating methane and ammonia production, rumen fermentation, and rumen bacteria in vitro. Anim. Feed Sci. Technol. 215, 25–36. doi: 10.1016/j.anifeedsci.2016.02.008

CrossRef Full Text | Google Scholar

Dawson, K. A., Rasmussen, M. A., and Allison, M. J. (1997). “Digestive disorders and nutritional toxicity,” in The Rumen Microbial Ecosystem, 2nd Edn, eds P. N. Hobson and C. S. Stewart (London, England: Blackie Academic and Professional), 633–660. doi: 10.1007/978-94-009-1453-7_14

CrossRef Full Text | Google Scholar

Duarte, A. C., Holman, D. B., Alexander, T. W., Durmic, Z., Vercoe, P. E., and Chaves, A. V. (2017). The type of forage substrate preparation included as substrate in a RUSITEC system affects the ruminal microbiota and fermentation characteristics. Front. Microbiol. 8:704. doi: 10.3389/fmicb.2017.00704

PubMed Abstract | CrossRef Full Text | Google Scholar

Eugène, M., Archimède, M., and Sauvant, D. (2004). Quantitative meta-analysis on the effects of defaunation of the rumen on growth, intake and digestion in ruminants. Livest. Prod. Sci. 85, 81–97. doi: 10.1016/S0301-6226(03)00117-9

CrossRef Full Text | Google Scholar

Finlay, B. J., Esteban, G., Clarke, K. J., Williams, A. G., Embley, T. M., and Hirt, R. P. (1994). Some rumen ciliates have endosymbiotic methanogens. FEMS Microbiol. Lett. 117, 157–161. doi: 10.1111/j.1574-6968.1994.tb06758.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Goel, G., and Makkar, H. P. S. (2012). Methane mitigation from ruminants using tannins and saponins. Trop. Anim. Health Pro. 44, 729–739. doi: 10.1007/s11250-011-9966-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Granja-Salcedo, Y. T., Ribeiro Júnior, C. S., De Jesus, R. B., Gomez-Insuasti, A. S., Rivera, A. R., Messana, J. D., et al. (2016). Effect of different levels of concentrate on ruminal microorganisms and rumen fermentation in Nellore steers. Arch. Anim. Nutr. 70, 17–32. doi: 10.1080/1745039X.2015.1117562

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W. L., Liu, J. X., Ye, J. A., Wu, Y. M., and Guo, Y. Q. (2005). Effect of tea saponin on rumen fermentation in vitro. Anim. Feed Sci. Technol. 120, 333–339. doi: 10.1016/j.anifeedsci.2005.02.029

CrossRef Full Text | Google Scholar

Hua, C. F., Tian, J., Ping, T., Cong, R. H., Luo, Y. W., Geng, Y. L., et al. (2017). Feeding a high concentration diet induces unhealthy alterations in the composition and metabolism of ruminal microbiota and host response in a goat model. Front. Microbiol. 8:138. doi: 10.3389/fmicb.2017.00138

PubMed Abstract | CrossRef Full Text | Google Scholar

IPCC. (2007). “Climate Change 2007: Synthesis Report,” in Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, eds R. K. Pachauri and A. Reisinger (Geneva, Switzerland: IPCC), 104.

Google Scholar

Ishaq, S. L., and Wright, A. D. (2014). Design and validation of four new primers for next-generation sequencing to target the 18S rRNA genes of gastrointestinal ciliate protozoa. Appl. Environ. Microbiol. 80, 5515–5521. doi: 10.1128/AEM.01644-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarvis, G. N., Strömpl, C., Burgess, D. M., Skillman, L. C., Moore, E. R. B., and Joblin, K. N. (2000). Isolation and identification of ruminal methanogens from grazing cattle. Curr. Microbiol. 40, 327–332. doi: 10.1007/s002849910065

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. S., Lee, Y. H., Kim, Y. I., Ahmadi, F., Oh, Y. K., Park, J. M., et al. (2016). Effect of microbial inoculant or molasses on fermentative quality and aerobic stability of sawdust-based spent mushroom substrate. Bioresour. Technol. 216, 188–195. doi: 10.1016/j.biortech.2016.05.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Choudhury, P. K., Carro, M. D., Griffith, G. W., Dagar, S. S., Puniya, M., et al. (2014). New aspects and strategies for methane mitigation from ruminants. Appl. Microbiol. Biotechnol. 98, 31–44. doi: 10.1007/s00253-013-5365-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R. H., Teng, Z. W., Lang, C. L., Zhou, H. Z., Zhong, W. G., Ban, Z. B., et al. (2019). Effect of different forage-to-concentrate ratios on ruminal bacterial structure and real-time methane production in sheep. PLoS One 14:e0214777. doi: 10.1371/journal.pone.0214777

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q. H., Li, X. Y., Desta, S. T., Zhang, J. G., and Shao, T. (2016a). Effects of Lactobacillus plantarum and fibrolytic enzyme on the fermentation quality and in vitro digestibility of total mixed rations silage including rape straw. J. Integr. Agric. 15, 2087–2096. doi: 10.1016/S2095-3119(15)61233-3

CrossRef Full Text | Google Scholar

Liu, Q. H., Shao, T., and Bai, Y. F. (2016b). The effect of fibrolytic enzyme, Lactobacillus plantarum and two food antioxidants on the fermentation quality, alpha-tocopherol and beta-carotene of high moisture napier grass silage ensiled at different temperatures. Anim. Feed Sci. Technol. 221, 1–11. doi: 10.1016/j.anifeedsci.2016.08.020

CrossRef Full Text | Google Scholar

Mackie, R. I., Aminov, R. I., Hu, W. P., Klieve, A. V., Ouwerkerk, D., Sundset, M. A., et al. (2003). Ecology of uncultivated oscillospira species in the rumen of cattle, sheep, and reindeer as assessed by microscopy and molecular mpproaches. Appl. Environ. Microbiol. 69, 6808–6815. doi: 10.1128/AEM.69.11.6808-6815.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, S., Zhang, M., Liu, J., and Zhu, W. (2015). Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: membership and potential function. Sci. Rep. 5:16116. doi: 10.1038/srep16116

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Fernández, G., Abecia, L., Ramos-Morales, E., Martin-García, A. I., Molina-Alcaide, E., and Yáñez-Ruiz, D. R. (2014). Effects of propyl propane thiosulfinate on nutrient utilization, ruminal fermentation, microbial population and methane emissions in goats. Anim. Feed Sci. Technol. 191, 16–25. doi: 10.1016/j.anifeedsci.2014.01.008

CrossRef Full Text | Google Scholar

Martinez-Fernandez, G., Denman, S. E., Cheung, J., and Mcsweeney, C. S. (2017). Phloroglucinol degradation in the rumen promotes the capture of excess hydrogen generated from methanogenesis inhibition. Front. Microbiol. 8:1871. doi: 10.3389/fmicb.2017.01871

PubMed Abstract | CrossRef Full Text | Google Scholar

Menke, K. H., and Steingass, H. (1988). Estimation of the energetic feed value obtained from chemical analysis and in vitro gas production using rumen fluid. Anim. Res. Dev. 28, 7–55.

Google Scholar

Ministry of Agriculture of P. R. China [MOA] (2004). Feeding Standard of Dairy Cattle (NY/T 34-2004). Beijing: MOA.

Google Scholar

Moss, A. R., Jouany, J. P., and Newbold, J. (2000). Methane production by ruminants: its contribution to global warming. Ann. Zootech. 49, 231–253. doi: 10.1051/animres:2000119

CrossRef Full Text | Google Scholar

Ranilla, M. J., Jouany, J. P., and Morgavi, D. P. (2007). Methane production and substrate degradation by rumen microbial communities containing single protozoal species in vitro. Lett. Appl. Microbiol. 45, 675–680.

Google Scholar

Russell, J. B. (1998). The Importance of pH in the regulation of ruminal acetate to propionate ratio and methane production in vitro. J. Dairy Sci. 81, 3222–3230.

Google Scholar

Russell, J. B., and Wallace, R. J. (1997). “Energy-yielding and energy-consuming reactions,” in The Rumen Microbial Ecosystem, 2nd Edn, eds P. N. Hobson and C. S. Stewart (London, England: Blackie Academic and Professional), 246–282.

Google Scholar

Shabat, S. K. B., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Miller, M. E. B., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972.

Google Scholar

Sniffen, C. J., O’Connor, J. D., Van Soest, P. J., Fox, D. G., and Russell, J. B. (1992). A net carbohydrate and protein system for evaluating cattle diets: II. Carbohydrate and protein availability. J. Anim. Sci. 70, 3562–3577.

Google Scholar

Stewart, C. S., and Duncan, S. H. (1985). The effect of avoparcin on cellulolytic bacteria of the ovine rumen. Microbiology 131, 427–435.

Google Scholar

Sun, Y. K., Yan, X. G., Ban, Z. B., Yang, H. M., Hegarty, R. S., and Zhao, Y. M. (2017). The effect of cysteamine hydrochloride and nitrate supplementation on in-vitro and in-vivo methane production and productivity of cattle. Anim. Feed Sci. Technol. 232, 49–56. doi: 10.1016/j.anifeedsci.2017.03.016

CrossRef Full Text | Google Scholar

Sylvester, J. T., Karnati, S. K. R., Yu, Z. T., Morrison, M., and Firkins, J. L. (2004). Development of an assay to quantify rumen ciliate protozoal biomass in cows using real-time PCR. J. Nutr. 134, 3378–3384.

Google Scholar

Takahashi, J. (2001). Nutritional manipulation of methanogenesis in ruminants. Asian-Austral. J. Anim. Sci. 14, 131–135.

Google Scholar

Takahashi, J., Chaudhry, A. S., Beneke, R. G., Suhubdy, and Young, B. A. (1997). Modification of methane emission in sheep by cysteine and a microbial preparation. Sci. Total Environ. 204, 117–123.

Google Scholar

Van Nevel, C. J., and Demeyer, D. I. (1996). Control of rumen methanogenesis. Environ. Monit. Assess. 42, 73–97. doi: 10.1007/BF00394043

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Waghorn, G. C. (1991). Bloat in cattle 47. Relationships between intra-ruminal pressure, distension, and the volume of gas used to simulate bloat in cows. N. Z. J. Agric. Res. 34, 213–220.

Google Scholar

Wallace, R. J. (2004). Antimicrobial properties of plant secondary metabolites. Proc. Nutr. Soc. 63, 621–629.

Google Scholar

Wang, L. J., Zhang, G. N., Xu, H. J., Xin, H. S., and Zhang, Y. G. (2019). Metagenomic analyses of microbial and carbohydrate-active enzymes in the rumen of holstein cows fed different forage-to-concentrate ratios. Front. Microbiol. 10:649. doi: 10.3389/fmicb.2019.00649

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, M. M., Zhang, R. B., Wang, Y. G., Ji, H. G., Zheng, Jia, Chen, X. H., et al. (2013). Microbial community structure and diversity in deep-sea hydrothermal vent sediments along the Eastern Lau Spreading Centre. Acta Oceanol. Sin. 32, 42–51. doi: 10.1007/s13131-013-0276-6

CrossRef Full Text | Google Scholar

Wen, A. Y., Yuan, X. J., Wang, J., Desta, S. T., and Shao, T. (2017). Effects of four short-chain fatty acids or salts on dynamics of fermentation and microbial characteristics of alfalfa silage. Anim. Feed Sci. Technol. 223, 141–148.

Google Scholar

Wolin, M. (1975). “Interactions between the bacterial species of the rumen,” in Digestion & Metabolism in the Ruminant, eds I. W. McDonald and A. C. I. Warner (Armidale, Australia: University of New England Publication), 134–148.

Google Scholar

Xing, L., Chen, L. J., and Han, L. J. (2009). The effect of an inoculant and enzymes on fermentation and nutritive value of sorghum straw silages. Bioresour. Technol. 100, 488–491. doi: 10.1016/j.biortech.2008.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Yanagita, K., Manome, A., Meng, X. Y., Hanada, S., Kanagawa, T., Tsuchida, T., et al. (2003). Flow cytometric sorting, phylogenetic analysis and in situ detection of Oscillospira guillermondii, a large, morphologically conspicuous but uncultured ruminal bacterium. Int. J. Syst. Evol. Microbiol. 53, 1609–1614. doi: 10.1099/ijs.0.02541-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, X. J., Guo, G., Wen, A. Y., Desta, S. T., Wang, J., Wang, Y., et al. (2015). The effect of different additives on the fermentation quality, in vitro digestibility and aerobic stability of a total mixed ration silage. Anim. Feed Sci. Technol. 207, 41–50. doi: 10.1016/j.anifeedsci.2015.06.001

CrossRef Full Text | Google Scholar

Yuan, X. J., Wen, A. Y., Desta, S. T., Dong, Z. H., and Shao, T. (2016). Effects of 4 short-chain fatty acids or salts on dynamics of nitrogen transformations and intrinsic protease activity of alfalfa silage: Effect of organic salts on proteolysis. J. Sci. Food Agric. 97, 2759–2766.

Google Scholar

Zhang, G. N., Li, Y., Fang, X. P., Cai, Y. M., and Zhang, Y. G. (2020b). Lactation performance, nitrogen utilization, and profitability in dairy cows fed fermented total mixed ration containing wet corn gluten feed and corn stover in combination replacing a portion of alfalfa hay. Anim. Feed Sci. Technol. 269:114687. doi: 10.1016/j.anifeedsci.2020.114687

CrossRef Full Text | Google Scholar

Zhang, G. N., Fang, X. P., Feng, G. Z., Li, Y., and Zhang, Y. G. (2020a). Silage fermentation, bacterial community, and aerobic stability of total mixed ration containing wet corn gluten feed and corn stover prepared with different additives. Animals 10:1775.

Google Scholar

Zhou, R., Wu, J. P., Lang, X., Liu, L. S., Casper, D. P., Wang, C. L., et al. (2020). Effects of oregano essential oil on in vitro ruminal fermentation, methane production, and ruminal microbial community. J. Dairy Sci. 103, 2303–2314.

Google Scholar

Zhou, Z. M., Fang, L., Meng, Q. X., Li, S. L., Chai, S. T., Liu, S. J., et al. (2017). Assessment of ruminal bacterial and archaeal community structure in Yak (Bos grunniens). Front. Microbiol. 8:179. doi: 10.3389/fmicb.2017.00179

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: fermented total mixed ration, carbohydrate component, methane yield, ruminal fermentation, rumen microbiome

Citation: Li Y, Lv J, Wang J, Zhou S, Zhang G, Wei B, Sun Y, Lan Y, Dou X and Zhang Y (2021) Changes in Carbohydrate Composition in Fermented Total Mixed Ration and Its Effects on in vitro Methane Production and Microbiome. Front. Microbiol. 12:738334. doi: 10.3389/fmicb.2021.738334

Received: 26 July 2021; Accepted: 08 October 2021;
Published: 05 November 2021.

Edited by:

Mukesh Kumar Awasthi, Northwest A&F University, China

Reviewed by:

Mariana Ornaghi, State University of Maringá, Brazil
Yury Tatiana Granja-Salcedo, Colombian Corporation for Agricultural Research (AGROSAVIA), Colombia

Copyright © 2021 Li, Lv, Wang, Zhou, Zhang, Wei, Sun, Lan, Dou and Zhang. 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: Xiujing Dou, douxiujing@neau.edu.cn; Yonggen Zhang, zhangyonggen@sina.com

Disclaimer: 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.