Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 07 April 2022
Sec. Microbiotechnology

The Effect of Early and Delayed Harvest on Dynamics of Fermentation Profile, Chemical Composition, and Bacterial Community of King Grass Silage

\r\nShihe Long&#x;Shihe Long1†Xuefeng Li&#x;Xuefeng Li2†Xianjun YuanXianjun Yuan3Rina SuRina Su4Junxin PanJunxin Pan1Ye ChangYe Chang1Mengli ShiMengli Shi1Zhihai CuiZhihai Cui1Naixin HuangNaixin Huang1Jian Wang*\r\nJian Wang1*
  • 1College of Animal Science and Technology, Hainan University, Haikou, China
  • 2College of Forestry, Hainan University, Haikou, China
  • 3Institute of Ensiling and Processing of Grass, Nanjing Agricultural University, Nanjing, China
  • 4State Key Laboratory of Grassland Agro-ecosystems, School of Life Sciences, Lanzhou University, Lanzhou, China

The objective of this study was to assess the effect of harvesting time on the fermentation characteristics, chemical composition, and microbial community of king grass silage. King grass was harvested at three growth periods of 90 days (KN90S), 110 days (KN110S), and 130 days (KN130S); chopped into 2–3-cm particle size; and ensiled in polyethylene bags (20 × 30 cm). The fermentation quality and chemical composition of silages were analyzed after 1, 3, 7, 14, 30, and 60 days of ensiling. Bacterial community of silage ensiled for 60 days was profiled using next generation sequencing (NGS) technology. The KN110S showed the most extensive lactic acid (LA) fermentation during 7 days of fermentation compared to KN90S and KN130S. After 60 days of fermentation, the KN110S showed the lowest pH and the highest lactic acid content among the three treatments. The butyric acid and ammonia nitrogen contents of KN90S and KN130S were significantly greater than those of KN110S (p < 0.05). After a timespan of 60 days of ensiling, the bacterial community of king grass silage was predominantly populated by Proteobacteria in phylum level, whereas unclassified Enterobacteriaceae genus remained dominant in all silages. A higher relative abundance of Clostridium was observed in KN90S and KN130S, but not in KN110S, and greater abundance of Lactococcus appeared in KN110S and KN130S silages than KN90S. It is concluded that harvesting time had an important effect on the fermentation quality and microbial community of king grass silage.

Introduction

King grass (Pennisetum purpureum Rich × Pennisetum americanum) is a major feed resource for ruminants in tropical regions due to its high nutritional values (Li et al., 2014) and can be harvested for five to eight times per year, with high biomass productions as much as 75–180 t of fresh matter (FM)/ha per year, which makes it an important potential resource in the development of animal feeding systems, production of biofuels, composts, or as a substrate for biodigestion (Botero-Londoño et al., 2021). In Hainan, utmost weather events such as typhoons and extreme rainfall frequently occurred in recent years, forcing farmers to adjust the forage harvest time to avoid agricultural losses. Therefore, proper preservation of king grass is one of the critical strategies to provide forage with consistent nutrients for ruminants and ensure a continuous supply of bioenergy production.

Ensiling is a prevailing forage preservation method for animal production due to its long storage duration, good palatability, and easy operation (McDonald et al., 1991). The provision of high-quality silage is the key to sustaining high productivity in intensive livestock pasture-based systems. The growth stage had been characterized as the main factor affecting silage quality because the dry matter (DM) production, ensilability, chemical composition, and nutritive value of forage are varied along with its change (Nazli et al., 2019). Early harvest crops were more suitable for ensiling because of higher water-soluble carbohydrates (WSC) (Stirling et al., 2021), but they showed a lower DM yield than delayed harvest (Xie et al., 2012). Besides, clostridial activity will be encouraged by the high moisture content of early harvested crops, resulting in considerable DM and energy losses (McDonald et al., 1991). Delayed harvest increases neutral detergent fiber (NDF) and acid detergent fiber (ADF) contents of crops, which is inconducive to preservation. At the same time, the WSC content tends to decrease and slow down the fermentation, retarding the decrease of pH necessary for efficient preservation (Kalač, 2011).

Complex microbial interactions directly determine the natural preservation process in terms of ensiling. To get an insight into the microbial community of silages is critical to distinguish valuable microorganisms and understand the underlying scientific basis for optimal silage making (Nishino, 2015; Guo et al., 2018). In this context, multiple molecular techniques have been applied to investigate the microbial ecology of silages (McAllister et al., 2018), among them is the next generation sequencing (NGS) technique, an advanced tool used for identifying microbial community in many crops and forages (Zhao et al., 2017; Keshri et al., 2019; Xu et al., 2019; Bai et al., 2021). However, it is still unclear whether different growth stages affect bacterial community in king grass silage.

It is hypothesized that the growth stage affected the chemical composition and microbial community of king grass, which in turn affected its fermentation quality. In the present study, we estimate the response of chemical composition and microbial community of king grass silage at different growth stages related to harvest time, which could provide a theoretical basis for improving king grass silage via control of harvest time.

Materials and Methods

Raw Materials and Silage Preparation

The king grass (cultivar Reyan No. 4) was planted in the experimental field of Hainan University (Haikou, China, 20°06′N, 110°32′E) on June 20, 2020. After 90, 110, and 130 days of growth, three plots of the forage were harvested at noon to avoid the variation in the diurnal rhythm of nutrients induced by photosynthesis and respiration of herbage. Fresh forage was immediately transported to the laboratory, cut into a theoretic length of 20–30 mm using a forage chopper, and mixed evenly. Approximately 250 g of fresh mix was ensiled in polyethylene bags (20 × 30 cm; Dongguan Bojia Packing Products Co., Ltd., China) and vacuumed tightly using a vacuum sealer (YS-10 L; Fujian Jiamei Machinery Co., Ltd., China). Each treatment uses 18 silage bags (six ensiling durations × three replicates). All bags were stored at an ambient temperature and sampled after 1, 3, 7, 14, 30, and 60 days of ensiling. Bacterial community was profiled in silage ensiled for 60 days. Additionally, a subsample of fresh materials was sampled for the analysis of chemical characteristics and microbial counts.

Fermentation Parameter, Chemical Composition, and Microbial Enumeration

Silage samples of 20 g were homogenized with 180 ml of distilled water and stored at 4°C overnight, followed by filtering through four layers of cheesecloth. The filtrate was used to determine pH and the concentrations of organic acid and ammonia nitrogen (AN). The pH was measured immediately using a glass electrode pH meter (PHS-3C; INESA Scientific Instrument, Shanghai). The organic acids, including lactic acid (LA), acetic acid (AA), propionic acid (PA), and butyric acid (BA), were measured using high-performance liquid chromatography (LC20A, Shimadzu, Ltd., Tokyo, Japan), equipped with a UV detector (SPD-M10AVP) (column, Carbomix H-NP; mobile phase, 2.5 mmol/l of H2SO4, pH 2.6; flow rate, 0.6 ml/min; temperature, 55°C; injection volume, 5 μl; and SPD, 210 nm). The AN concentration was analyzed as described by Weatherburn (1967). To determine the DM content, a subsample of fresh and silage was oven-dried at 65°C for 72 h. The dried samples were then ground through a 1-mm screen by a feed grinder (YBF-2500; Wande, China) for the following chemical composition analyses. Crude protein (CP) concentration was quantified by an automatic Kjeldahl nitrogen analyzer (K9860, Shandong Haineng Scientific Instrument Co., Ltd., China). The concentrations of NDF and ADF were measured as reported by Van Soest et al. (1991), and WSC were measured according to Owens et al. (1999).

Using the pour plate technique for microbial enumeration (Man et al., 1960), 20 g of fresh samples were blended with 180 ml of sterilized saline (8.5 g/l) and serially diluted to 10–fold. Lactic acid bacteria (LAB) were anaerobically incubated at 37°C for 48 h on de Man, Rogosa and Sharpe agar (Solarbio Life Sciences Co., Ltd., Beijing, China). Coliform bacteria were aerobically incubated at 37°C for 48 h on Violet Red Bile agar (Solarbio Life Sciences Co., Ltd., Beijing, China). Yeasts and molds were aerobically cultured at 25°C for 72 h on Rose Bengal agar (Solarbio Life Sciences Ltd. Co., Beijing, China).

Microbial Community Analysis

After 60 days of fermentation, 20 g of silage sample was homogenized with 100 ml of sterile sodium chloride solution at 120 for 15 min, followed by filtering through four layers of cheesecloth. Total DNA was extracted by TIANamp Bacteria DNA isolation kit (DP302–02, Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer’s protocols. The quality of the extracted DNA was checked by a NanoPhotometer N50 UV-Vis spectrophotometry (Implen GmbH, Germany). All extracted DNA samples were frozen at −80°C for further analysis. The V3–V4 region of 16S rRNA was amplified with primers 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATCTAATCC-3′). The reads were merged into a sequence using PEAR (version 0.9.8), removing the bases with a mass value below 20 at the tail of reads through PRINSEQ (version 0.20.4) to obtain clean data. Operational taxonomic unit (OTU) clustering was performed with 97% of similarity through the Usearch program (version 11.0.667).1 The Ribosome Database Project (RDP) algorithm was used to classify the representative sequences of each OTU at different levels. The alpha diversity including Shannon index, Chao richness estimator, Ace, and Good’s coverage index was calculated using Mothur (version 1.43.0).2 Venn diagrams (version 1.6.20) was used to construct the Venn diagram, and the principal coordinates analysis (PCoA) presented in this study was conducted with the VEGEN package in R software (version 3.6.0) (Dixon, 2003).

Statistical Analysis

Microbial populations were estimated as colony-forming units (CFU)/gram of forage, and logarithmic conversion was performed before statistical analysis. The fermentation profiles of the silages were analyzed in a completely randomized block design with 3 × 6 factorial treatment arrangements with three replicates per treatment. The prime effect of the growth period, ensiling duration, and their interactions were analyzed by the general linear model (GLM) procedure of SAS 9.3 software (SAS Institute Inc., 2019). Significance was declared at p < 0.05 and tendencies at 0.05 < p ≤ 0.10.

Results

Chemical and Microbial Composition of King Grass Before Ensiling

The chemical and microbial composition of king grass before ensiling are depicted in Table 1. The growth period had a significant effect (p < 0.01) on contents of DM, WSC, CP, ADF, and NDF, as well as LAB, molds, yeasts, and coliform bacteria numbers. The WSC content of king grass ranged from 73.8 to 94.7 g/kg of DM. The undesirable microorganism (e.g., mold, yeast, and coliform bacteria) counts of KN110S were the lowest among three groups. The readings of LAB (3.8–4.4 log10 CFU/g of FM), ADF (385.8–439.5 g/kg of DM), and NDF (647.3–710.7 g/kg of DM) increased (p < 0.05), whereas CP (92.3–66.7 g/kg of DM) content decreased (p < 0.05) as the growth time prolonged.

TABLE 1
www.frontiersin.org

Table 1. Chemical characteristics and attached microorganisms of king grass before ensiling.

Fermentation Characteristics and Chemical Composition of King Grass Silage

As shown in Tables 2, 3 the ensiling duration and growth period had a significant effect (p < 0.01) on all fermentative parameters except for LA/AA, and there were interactions on the levels of pH, DM, WSC, CP, ADF, LA, AA, LA/AA, PA, BA, and AN. The DM contents of all king grass silage gradually decreased with ensiling time, and the marked decline (p < 0.05) was observed in KN90S. The CP concentrations of KN90S, KN110S, and KN130S (91.7 vs. 73.1 g/kg of DM; 86.4 vs. 80.2 g/kg of DM; 66.0 vs. 52.5 g/kg of DM, respectively) statistically reduced (p < 0.05) after 60 days of ensiling, whereas the ADF (average of 396.0 vs. 440.7 g/kg of DM) and NDF (average of 690.8 vs. 718.4 g/kg of DM) contents increased. For KN110S, the CP concentration decreased slightly to a smaller extent than those in KN90S and KN130S. The WSC concentrations in KN90S (71.4 vs. 28.6 g/kg of DM), KN110S (65.2 vs. 43.9 g/kg of DM), and KN130S (79.3 vs. 58.8 g/kg of DM) were dramatically reduced within the early stage of fermentation.

TABLE 2
www.frontiersin.org

Table 2. Changes in the chemical composition of king grass silage during ensiling.

TABLE 3
www.frontiersin.org

Table 3. Changes in the fermentation quality of king grass silage during ensiling.

Different extents of LA fermentation occurred among the three groups. The pH value of KN110S reached 4.5 after 7 days of ensiling, while the pH value of KN90S and KN130S remained above 5. AA contents of KN110S and KN130S presented an increasing trend during ensiling, and LA/AA values decreased from days 14 to 60, while a significant decrease in LA/AA exhibited in KN90S from days 1 to 3. A significant increase (p < 0.05) of PA, BA, and AN content was observed in king grass silage on day 60 compared with day 1, and the final PA contents of the three treatments were around 1.1 g/kg of DM. KN110S appeared to have a slower accumulation of BA and AN concentrations than KN90S and KN130S during ensiling and the lowest final BA (4.7 vs. 10.2 and 9.9 g/kg of DM) and AN (89.4 vs. 244.7 and 232.9 g/kg of TN) at 60 days.

Bacterial Diversity of King Grass Silage After 60 Days of Ensiling

The bacterial alpha diversity of king grass silage is shown in Table 4. A total of 500,278 sequences of bacteria were obtained using NGS technology. The Good’s coverage values of all groups were 0.99, and a non-significant difference was observed in OTUs, Shannon, Chao1, Ace, and Simpson indexes among KN90S, KN110S, and KN130S. The core OTUs were displayed in the Venn diagram (Figure 1). The overlapping OTUs for the bacterial communities among king grass silages at three different harvest times were 298, accounting for 91.69% of total OTUs. The bacterial community of king grass silages can be clearly demonstrated by PCoA. As shown in Figure 2, principal component 1 (PCoA1) and component 2 (PCoA2) can be used to explain 66.7 and 16.7% of the total variance, respectively, indicating that the bacterial composition of silages from different growth periods could be well separated from each other.

TABLE 4
www.frontiersin.org

Table 4. Alpha diversity of bacterial diversity for king grass silage.

FIGURE 1
www.frontiersin.org

Figure 1. Venn diagram of the core OTUs of king grass silage (KN90S, king grass silage harvested at growth periods of 90 days; KN110S, king grass silage harvested at growth periods of 110 days; KN130S, king grass silage harvested at growth periods of 130 days; OTU, operational taxonomic unit).

FIGURE 2
www.frontiersin.org

Figure 2. Principal coordinates analysis (PCoA) of the bacterial community for king grass silages (KN90S, king grass silage harvested at growth periods of 90 days; KN110S, king grass silage harvested at growth periods of 110 days; KN130S, king grass silage harvested at growth periods of 130 days).

The bacterial communities at the phylum level revealed that Proteobacteria was the most abundant phylum (average of 65.17%) in king grass silage (Figure 3). Conversely, Firmicutes only occupied an average of 34.83% of the population. The bacterial composition of king grass silage at the genus level was also analyzed to better understand the bacterial community of king grass silage (Figure 4). The top three genera of the KN110S silage were unclassified-Enterobacteriaceae (55.0%), Lactobacillus (27.9%), and Lactococcus (6.9%), while unclassified-Enterobacteriaceae (52.8 and 67.0%), Lactobacillus (28.0 and 7.6%), and Clostridium (7.4 and 9.8%) were the first three dominating genera observed only in KN90S and KN130S. A greater abundance of Lactobacillus than Lactococcus and Pediococcus (average of 22.794.15, and 0.80%, respectively) appeared in king grass silage. The relative abundance of Lactococcus in KN90S, KN110S, and KN130S was 0.9, 6.9, and 3.7%, respectively. A higher relative abundance of Clostridium was found in KN90S (7.4%) and KN130S (9.8%) silages, which was consistent with the greater BA and AN content of these groups.

FIGURE 3
www.frontiersin.org

Figure 3. Differences in the bacterial community of king grass silage at the phylum level (KN90S, king grass silage harvested at growth periods of 90 days; KN110S, king grass silage harvested at growth periods of 110 days; KN130S, king grass silage harvested at growth periods of 130 days).

FIGURE 4
www.frontiersin.org

Figure 4. Differences in the bacterial community of king grass silage at the genus level (KN90S, king grass silage harvested at growth periods of 90 days; KN110S, king grass silage harvested at growth periods of 110 days; KN130S, king grass silage harvested at growth periods of 130 days).

Discussion

Chemical and Microbial Composition of King Grass

Sufficient epiphytic LAB (105 CFU/g FM) and WSC contents (5% DM) of fresh material were considered as essential factors for assuring the acceptable fermentation quality (McDonald et al., 1991). The WSC content of king grass in the present study may provide enough fermentable substrate for the extensive fermentation. However, the number of LAB was lower than the theoretical requirement for good preservation. Previous findings reported that the number of microorganisms on the surface of plants varied with the growth periods (Xie et al., 2012; Lounglawan et al., 2014), and the epiphytic LAB increased with delayed harvest dates (Müller, 2009). The chemical composition and nutritional value of forages change along with the growth period (Nazli et al., 2019). The decline in CP and the increase in ADF and NDF contents in king grass with the extension of growth days were in line with that of Botero-Londoño et al. (2021) in king grass with different cutting ages. They found a linear increase in fiber levels with the age of cutting and a significant decrease in crude protein levels. The decrease in CP was due to an increase in the proportion of stems (Ribeiro et al., 2014), while the increase in fiber concentration was due to the lignification process as harvest dates progressed (Darby and Lauer, 2002).

Fermentation Characteristics and Chemical Composition of King Grass Silage

The loss of dry matter during ensiling will reduce the nutritive value of silage. DM showed significant decreases after 60 days of ensiling, which is in agreement with the report by Wang et al. (2020) in Italian ryegrass silages. Savoie et al. (2002) reported that mature corn silage (108 days of growth) produced less effluent (1.41 vs. 2.47% of the original wet mass) than premature corn silage (94 days of growth) on a laboratory scale. The decrease in DM content during ensiling was due to the activity of undesirable bacteria, which metabolized soluble carbohydrates to CO2, causing DM losses (Ávila et al., 2014), while the severe decline in DM in KN90S during 7 days of fermentation was due to effluent production. High-moisture silage often bears a high risk of effluent production, leading to high DM loss (Zhao et al., 2021), and the peak effluent of silage typically occurs within 10 days after ensiling given the time required for plant cell to rupture (Gebrehanna et al., 2014). This is consistent with the view of Savoie et al. (2002), who found that the amount of daily effluent collected from the timothy silage was relatively high within 7 days of ensiling, and the effluent flow declined in subsequent days.

The degradation of protein is inevitable when silage is made, as a result a combination of plant enzymes and microbial metabolism (McDonald et al., 1991). According to the higher BA and AN content, statistically reduced CP concentrations in this study reflect that king grass underwent extensive fermentation of Clostridium butyricum and protein degradation during ensiling (Kung et al., 2018). This might be explained by the fact that silage with a high moisture content bears a great risk of extensive clostridial fermentation, resulting in a large amount of protein hydrolysis and BA occurrence (He et al., 2020a). A relatively low degree of CP degradation in KN110S could be explained by a rapid drop in pH value. A previous study reported that protein breakdown and high concentrations of AN in silages are typically caused by a slow drop in pH (Ogunade et al., 2018). Similarly, Li et al. (2020) found that fast acidification of silage was conducive to preventing clostridial fermentation. Likewise, the protein degradation in alfalfa silage has been effectively inhibited by inoculated L. plantarum on account of the fast drop in pH that inhibits the activity of proteolytic bacteria (Filya et al., 2007). Fiber can hardly be degraded in silage due to its complex fibrous structures (Nair et al., 2020). In present study, we found an increased trend of ADF and NDF during ensiling, which was similar to those reported by Wang C. et al. (2019), but it was contrary to the study of Zi et al. (2021) who found a numerically declined ADF and NDF contents of king grass silage without additives during fermentation. This may be due to the loss of DM during ensiling, increasing the relative levels of ADF and NDF (Ferrero et al., 2019). Furthermore, fermentation of sugars into CO2 by microbial metabolism at the ensiling process can lead to increases in the fibrous fraction of silage (Ávila et al., 2008), and the metabolism of fermented substrates to metabolites (mainly organic acids) by microorganisms anaerobically can also increase the contents of ADF and NDF. The consumption rate of WSC correlates with the extent of fermentation. The rapid decline of WSC during the initial stage of ensiling indicates that all king grass underwent extensive fermentation.

Rapid acidification of silage depends on LA production by LAB, which is imperative to prevent the propagation of early undesirable microorganisms and the loss of nutrients (Liu et al., 2016). The most rapid decrease in pH and increase in LA content were found in KN110S, followed by KN130S and KN90S. This might be attributed to the lowest attached undesirable microorganisms in KN110S, and the lower BA concentration of KN110S also confirms the extensive LA fermentation of KN110S silages (Shao et al., 2007). However, the final pH value of all silages was still above the ideal silage pH of 4.2 at day 60. The finding was consistent with a previous work by Li D. et al. (2019) in several tropical forage silage after 60 days of fermentation including king grass. The rates of diffusion of WSC from intact and ruptured cells into the aqueous phase are rather important than the total amounts in the crop for silage making (McDonald et al., 1991), considering that there was adequate fermentable substrate of fresh king grass in this study. Distinct from temperate crops, the higher incorporation of air in the ensiling and the restriction of cell breakdown and juice release due to the properties of a tropical grass with coarse porosity and stemmy structures delayed the onset of the fermentation phase during the very early stage of the ensiling by extending respiration and aerobic microbial activity (Shao et al., 2005), causing the higher loss of WSC, which is detrimental to subsequent anaerobic lactic acid fermentation to drop pH (Dunière et al., 2013). A previous study has reported that AA-type fermentation was considered the primary pattern in tropical silage (Nishino et al., 2012), and AA was the main fermentation product in various tropical silages including paspalum, stylo, and cassava (Li D. et al., 2019; Li M. et al., 2019). The rapid increase in AA concentration and the marked decrease in LA/AA in KN90S silage are also reported by Yuan et al. (2020) in Napier grass silage from 3 to 5 days of ensiling. The production of AA may be partly related to the metabolism of AA-producing bacteria (e.g., Enterobacter) and the ability of homofermentative LAB to convert LA anoxically into AA and 1,2-propanediol (Oude Elferink et al., 2001). The production of PA in silages was likely caused by several microorganisms such as Propionibacterium sp. and Clostridium propionicum (Santos et al., 2016). A significant increase in the PA content was observed in king grass silage on day 60 compared to day 1, although this may not be of practical significance because of their low concentration. The final PA content of the three treatments were around 1.1 g/kg of DM, which is consistent with previous reports that observed a close PA content in king grass silage without additives (Li D. et al., 2019; Zi et al., 2021).

Bacterial Diversity of King Grass Silage

This study revealed the relative abundance and diversity of bacteria of king grass silage from different growth stages. The analysis of the alpha diversity of the bacterial community revealed the variance of the microbial community. The coverage values of all samples were 0.99, suggesting that the sequencing had adequately reflected the profile of the bacterial community. The highly overlapping OTUs in the present study reflect that many similar microbes were involved in silage. PCoA clearly demonstrates the variance of the microbial community by the clear separation of the three treatments, which indicates that growth periods had a remarkable effect on the bacterial community of king grass silage and that ensiling could be the main factor affecting fermentation. The most dominant phyla in the 60-day king grass silage were Proteobacteria, followed by Firmicutes. Proteobacteria in silage is comprised of the families Enterobacteriaceae, Sphingomonadaceae, Xanthomonadaceae, and Aurantimonadaceae and the genera Sphingomonas, Pantoea, Pseudomonas, and Stenotrophomonas (Ogunade et al., 2018). This phenomenon was inconsistent with previous studies reporting that the dominant phylum was shifted from Proteobacteria to Firmicutes along with the ensiling process (Xu et al., 2017; Yuan et al., 2020), although Wang et al. (2020) reported a dominant Proteobacteria (62.45%) phylum at 60 days in ensiling sorghum forages. To our knowledge, the ensiling conditions favored the selection of species belonging to Firmicutes because these species thrive under low pH and anaerobic conditions (Dong et al., 2019). This may be related to the high pH environment (>4.60), which is not low enough to enhance the shift from Proteobacteria to Firmicutes phylum after fermentation. Enterobacteria is a non-spore forming, facultative anaerobe that can ferment LA to AA and other products (Santos et al., 2016) and was usually outcompeted by Lactobacillus species with a decrease in pH during ensiling (Parvin et al., 2010), but we found a high relative abundance of the unclassified-Enterobacteriaceae that appeared in 60-day silage, which could be attributed to the high moisture of typical tropical grass. Wang Y. et al. (2019) suggested that a lower abundance of Enterobacter spp. might be connected to the high DM content of Moringa oleifera leaf (448 g/kg of FM) silages. Li D. et al. (2019) reported that Enterobacter was the predominant genus in high-moisture tropical forage including king grass, paspalum, and stylo silages, accompanied by high AA production. He et al. (2020a) made a similar report in which Enterobacter (39.3%) was the dominant genus of corn straw silage with high moisture content, and its relative abundance decreased from 39.3 to 28.9% and the AA content decreased significantly when the moisture content was lowered by mixing corn stalk with Bauhinia variegata flower. The high relative abundance of unclassified-Enterobacteriaceae in all silage samples might account for the high AA content (Nishino et al., 2012; Wang et al., 2020). Previous reports suggested that cocci such as Lactococcus, Leuconostoc, Pediococcus, and Enterococcus initiated lactic acid fermentation at the early stage of fermentation (Cai et al., 1998), and they were outcompeted by more acid-tolerant Lactobacillus species with the decline of pH as the ensiling time prolonged (Graf et al., 2016), which might partially explain the higher level of Lactobacillus than Lactococcus and Pediococcus in this study. The portion of Lactococcus remains high after ensiling, especially in the KN110S and KN130S silages. Xu et al. (2017) determined the bacterial community in corn stover silage by high-throughput sequencing and found that the abundance of Lactococcus in natural silage was higher than that of inoculated silage (5.69% vs. 0.37%), and the pH value of the inoculated group decreased more efficiently in the early stages of ensiling. Mu et al. (2020) detected greater relative abundances of Lactococcus in untreated mixed silage of amaranth and rice straw compared to those treated with additives. Likewise, Wang Y. et al. (2019) reported a minor presence of Lactococcus (less than 1%) in Moringa oleifera leaf silage with a pH below 4.2 after 60 days of ensiling. As aforementioned, a relatively slowly drop in pH and higher 60-day silage pH might explain the greater relative abundance of Lactococcus in KN110S and KN130S silages, whereas a low relative abundance of Lactococcus in KN90S may be due to the loss of soluble sugars caused by the effluent, resulting in a deficit of substrates that obstructed competitive colonization of Lactococcus in the early stage of fermentation. Clostridium was often detected in silages, and its presence is undesirable since it can utilize carbohydrates and proteins as a substrate for their growth and ferment LA to BA, resulting in DM and energy losses (McDonald et al., 1991). The higher relative abundance of Clostridium presented in this study is consistent with the greater BA and AN content. This might be because high-moisture grasses always face a greater risk of Clostridium fermentation (He et al., 2020b; Mu et al., 2020). Clostridium inhibition in KN110S (0.5%) supported the speculation that clostridial fermentation was restricted by a faster rate of pH decline compared to KN90S and KN130S.

Conclusion

The KN110S silage displayed a relatively higher LA content, lactic acid bacteria abundance at the genus level, lower pH values, contents of AN and BA, and abundance of Clostridium compared with the KN90S and KN130S silage after ensiling. The results of this study indicated that harvesting date plays a paramount role in the fermentation parameter and microbial community of king grass silage. We recommend that king grass with a growth period of 110 days provides silage with better fermentation quality.

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 below: https://www.ncbi.nlm.nih.gov/, PRJNA800375.

Author Contributions

SL and JW designed the experiments. SL, JP, YC, MS, ZC, and NH conducted the experiments. SL, YC, and MS analyzed the data. SL wrote the manuscript. SL, XY, XL, and RS reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 31860682).

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://drive5.com/uparse/
  2. ^ http://www.mothur.org/wiki/Classify.seqs

References

Ávila, C. L. S., Carvalho, B. F., Pinto, J. C., Duarte, W. F., and Schwan, R. F. (2014). The use of Lactobacillus species as starter cultures for enhancing the quality of sugar cane silage. J. Dairy Sci. 97, 940–951. doi: 10.3168/jds.2013-6987

PubMed Abstract | CrossRef Full Text | Google Scholar

Ávila, C. L. S., Pinto, J. C., Sugawara, M. S., Silva, M. S., and Schwan, R. F. (2008). Qualidade da silagem de cana-de-açúcar inoculada com uma cepa de Lactobacillus buchneri. Acta Sci. 30, 255–261. doi: 10.4025/actascianimsci.v30i3.1564

CrossRef Full Text | Google Scholar

Bai, J., Ding, Z., Ke, W., Xu, D., Wang, M., Huang, W., et al. (2021). Different lactic acid bacteria and their combinations regulated the fermentation process of ensiled alfalfa: ensiling characteristics, dynamics of bacterial community and their functional shifts. Microb. Biotechnol. 14, 1171–1182. doi: 10.1111/1751-7915.13785

PubMed Abstract | CrossRef Full Text | Google Scholar

Botero-Londoño, J. M., Celis-Celis, E. M., and Botero-Londoño, M. A. (2021). Nutritional quality, nutrient uptake and biomass production of Pennisetum purpureum cv. King grass. Sci. Rep. 11:13799. doi: 10.1038/s41598-021-93301-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Y., Benno, Y., Ogawa, M., Ohmomo, S., Kumai, S., and Nakase, T. (1998). Influence of Lactobacillus spp. from an inoculant and of Weissella and Leuconostoc spp. from forage crops on silage fermentation. Appl. Environ. Microbiol. 64, 2982–2987. doi: 10.1128/AEM.64.8.2982-2987.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

Darby, H. M., and Lauer, J. G. (2002). Harvest date and hybrid influence on corn forage yield, quality, and preservation. Agron. J. 94, 559–566. doi: 10.2134/agronj2002.0559

CrossRef Full Text | Google Scholar

Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x

CrossRef Full Text | Google Scholar

Dong, Z., Li, J., Chen, L., Wang, S., and Shao, T. (2019). Effects of freeze-thaw event on microbial community dynamics during red clover ensiling. Front. Microbiol. 10:1559. doi: 10.3389/fmicb.2019.01559

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunière, L., Sindou, J., Chaucheyras-Durand, F., Chevallier, I., and Thévenot-Sergentet, D. (2013). Silage processing and strategies to prevent persistence of undesirable microorganisms. Anim. Feed Sci. Technol. 182, 1–15. doi: 10.1016/j.anifeedsci.2013.04.006

CrossRef Full Text | Google Scholar

Ferrero, F., Prencipe, S., Spadaro, D., Gullino, M. L., Cavallarin, L., Piano, S., et al. (2019). Increase in aflatoxins due to Aspergillus section Flavi multiplication during the aerobic deterioration of corn silage treated with different bacteria inocula. J. Dairy Sci. 102, 1176–1193. doi: 10.3168/jds.2018-15468

PubMed Abstract | CrossRef Full Text | Google Scholar

Filya, I., Muck, R. E., and Contreras-Govea, F. E. (2007). Inoculant effects on alfalfa silage: fermentation products and nutritive value. J. Dairy Sci. 90, 5108–5114. doi: 10.3168/jds.2006-877

PubMed Abstract | CrossRef Full Text | Google Scholar

Gebrehanna, M. M., Gordon, R. J., Madani, A., VanderZaag, A. C., and Wood, J. D. (2014). Silage effluent management: a review. J. Environ. Manage. 143, 113–122. doi: 10.1016/j.jenvman.2014.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Graf, K., Ulrich, A., Idler, C., and Klocke, M. (2016). Bacterial community dynamics during ensiling of perennial ryegrass at two compaction levels monitored by terminal restriction fragment length polymorphism. J. Appl. Microbiol. 120, 1479–1491. doi: 10.1111/jam.13114

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, X. S., Ke, W. C., Ding, W. R., Ding, L. M., Xu, D. M., Wang, W. W., et al. (2018). Profiling of metabolome and bacterial community dynamics in ensiled Medicago sativa inoculated without or with Lactobacillus plantarum or Lactobacillus buchneri. Sci. Rep. 8:357. doi: 10.1038/s41598-017-18348-0

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Wang, C., Xing, Y., Zhou, W., Pian, R., Chen, X., et al. (2020a). Ensiling characteristics, proteolysis and bacterial community of high-moisture corn stalk and stylo silage prepared with bauhinia variegate flower. Bioresour. Technol. 296:122336. doi: 10.1016/j.biortech.2019.122336

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Chen, N., Lv, H., Wang, C., Zhou, W., Chen, X., et al. (2020b). Gallic acid influencing fermentation quality, nitrogen distribution and bacterial community of high-moisture mulberry leaves and stylo silage. Bioresour. Technol. 295:122255. doi: 10.1016/j.biortech.2019.122255

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalač, P. (2011). The required characteristics of ensiled crops used as a feedstock for biogas production: a review. J. Agrobiol. 28, 85–96. doi: 10.2478/v10146-011-0010-y

CrossRef Full Text | Google Scholar

Keshri, J., Chen, Y., Pinto, R., Kroupitski, Y., Weinberg, Z. G., and Sela, S. S. (2019). Bacterial dynamics of wheat silage. Front. Microbiol. 10:1532. doi: 10.3389/fmicb.2019.01532

PubMed Abstract | CrossRef Full Text | Google Scholar

Kung, L. Jr., Shaver, R. D., Grant, R. J., and Schmidt, R. J. (2018). Silage review: interpretation of chemical, microbial, and organoleptic components of silages. J. Dairy Sci. 101, 4020–4033. doi: 10.3168/jds.2017-13909

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Ni, K., Zhang, Y., Lin, Y., and Yang, F. (2019). Fermentation characteristics, chemical composition and microbial community of tropical forage silage under different temperatures. Asian Australas. J. Anim. Sci. 32, 665–674. doi: 10.5713/ajas.18.0085

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Zi, X., Tang, J., Zhou, H., and Cai, Y. (2019). Silage fermentation, chemical composition and ruminal degradation of king grass, cassava foliage and their mixture. Grassl. Sci. 65, 210–215. doi: 10.1111/grs.12235

CrossRef Full Text | Google Scholar

Li, M., Zi, X., Zhou, H., Hou, G., and Cai, Y. (2014). Effects of sucrose, glucose, molasses and cellulase on fermentation quality and in vitro gas production of king grass silage. Anim. Feed Sci. Technol. 197, 206–212. doi: 10.1016/j.anifeedsci.2014.06.016

CrossRef Full Text | Google Scholar

Li, R., Jiang, D., Zheng, M., Tian, P., Zheng, M., and Xu, C. (2020). Microbial community dynamics during alfalfa silage with or without clostridial fermentation. Sci. Rep. 10:17782. doi: 10.1038/s41598-020-74958-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q. H., Li, X. Y., Desta, S. T., Zhang, J. G., and Shao, T. (2016). 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

Lounglawan, P., Lounglawan, W., and Suksombat, W. (2014). Effect of cutting interval and cutting height on yield and chemical composition of king napier grass (Pennisetum Purpureum x Pennisetum Americanum). APCBEE Proc. 8, 27–31. doi: 10.1016/j.apcbee.2014.01.075

CrossRef Full Text | Google Scholar

Man, J. C. D., Rogosa, M., and Sharpe, M. E. (1960). A medium for the cultivation of lactobacilli. J. Appl. Bacteriol. 23, 130–135. doi: 10.1111/j.1365-2672.1960.tb00188.x

CrossRef Full Text | Google Scholar

McAllister, T. A., Dunière, L., Drouin, P., Xu, S., Wang, Y., Munns, K., et al. (2018). Silage review: using molecular approaches to define the microbial ecology of silage. J. Dairy Sci. 101, 4060–4074. doi: 10.3168/jds.2017-13704

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald, P., Henderson, A. R., and Heron, S. J. E. (1991). The Biochemistry of Silage. Marlow: Chalcombe Publication.

Google Scholar

Mu, L., Xie, Z., Hu, L., Chen, G., and Zhang, Z. (2020). Cellulase interacts with Lactobacillus plantarum to affect chemical composition, bacterial communities, and aerobic stability in mixed silage of high-moisture amaranth and rice straw. Bioresour. Technol. 315:123772. doi: 10.1016/j.biortech.2020.123772

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, C. E. (2009). Influence of harvest date of primary growth on microbial flora of grass herbages and haylage, and on fermentation and aerobic stability of haylage conserved in laboratory silos. Grass Forage Sci. 64, 328–338. doi: 10.1111/j.1365-2494.2009.00695.x

CrossRef Full Text | Google Scholar

Nair, J., Huaxin, N., Andrada, E., Yang, H. E., Chevaux, E., Drouin, P., et al. (2020). Effects of inoculation of corn silage with Lactobacillus hilgardii and Lactobacillus buchneri on silage quality, aerobic stability, nutrient digestibility, and growth performance of growing beef cattle. J. Anim. Sci. 98:skaa267. doi: 10.1093/jas/skaa267

PubMed Abstract | CrossRef Full Text | Google Scholar

Nazli, M. H., Halim, R. A., Abdullah, A. M., Hussin, G., and Samsudin, A. A. (2019). Potential of four corn varieties at different harvest stages for tropical silage production in Malaysia. Asian Australas. J. Anim. Sci. 32, 224–232. doi: 10.5713/ajas.18.0175

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishino, N. (2015). “New trends in silage technology,” in Proceedings of the XVII International Silage Conference, (Piracicaba: ESALQ USP). doi: 10.3390/microorganisms9020274

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishino, N., Li, Y., Wang, C., and Parvin, S. (2012). Effects of wilting and molasses addition on fermentation and bacterial community in guinea grass silage. Lett. Appl. Microbiol. 54, 175–181. doi: 10.1111/j.1472-765X.2011.03191.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogunade, I. M., Jiang, Y., Pech Cervantes, A. A., Kim, D. H., Oliveira, A. S., Vyas, D., et al. (2018). Bacterial diversity and composition of alfalfa silage as analyzed by Illumina MiSeq sequencing: Effects of Escherichia coli O157:H7 and silage additives. J. Dairy Sci. 101, 2048–2059. doi: 10.3168/jds.2017-12876

PubMed Abstract | CrossRef Full Text | Google Scholar

Oude Elferink, S. J. W. H., Krooneman, J., Gottschal, J. C., Spoelstra, S. F., Faber, F., and Driehuis, F. (2001). Anaerobic conversion of lactic acid to acetic acid and 1, 2-propanediol by Lactobacillus buchneri. Appl. Environ. Microbiol. 67, 125–132. doi: 10.1128/AEM.67.1.125-132.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Owens, V. N., Albrecht, K. A., Muck, R. E., and Duke, S. H. (1999). Protein degradation and fermentation characteristics of red clover and alfalfa silage harvested with varying levels of total nonstructural carbohydrates. Crop Sci. 39, 1873–1880. doi: 10.2135/cropsci1999.3961873x

PubMed Abstract | CrossRef Full Text | Google Scholar

Parvin, S., Wang, C., Li, Y., and Nishino, N. (2010). Effects of inoculation with lactic acid bacteria on the bacterial communities of Italian ryegrass, whole crop maize, guinea grass and rhodes grass silages. Anim. Feed Sci. Technol. 160, 160–166. doi: 10.1016/j.anifeedsci.2010.07.010

CrossRef Full Text | Google Scholar

Ribeiro, G. O. Jr., Teixeira, A. M., Velasco, F. O., Faria Júnior, W. G., Pereira, L. G. R., Chaves, A. V., et al. (2014). Production, nutritional quality and in vitro methane production from andropogon gayanus grass harvested at different maturities and preserved as hay or silage. Asian Australas. J. Anim. Sci. 27, 330–341. doi: 10.5713/ajas.2013.13161

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, A. O., Ávila, C. L. S., Pinto, J. C., Carvalho, B. F., Dias, D. R., and Schwan, R. F. (2016). Fermentative profile and bacterial diversity of corn silages inoculated with new tropical lactic acid bacteria. J. Appl. Microbiol. 120, 266–279. doi: 10.1111/jam.12980

PubMed Abstract | CrossRef Full Text | Google Scholar

SAS Institute Inc (2019). SAS OnlineDoc§9.3.1, Edn. SAS Institute Inc. 9.3, 1 Edn. Cary, NC: SAS Institute Inc.

Google Scholar

Savoie, P., Amyot, A., and Thériault, R. (2002). Effect of moisture content, chopping, and processing on silage effluent. Trans. ASAE. 45, 907–914. doi: 10.13031/2013.9937

CrossRef Full Text | Google Scholar

Shao, T., Zhang, Z. X., Shimojo, M., Wang, T., and Masuda, Y. (2005). Comparison of fermentation characteristics of Italian ryegrass (Lolium multiflorum Lam.) and guineagrass (Panicum maximum Jacq.) during the early stage of ensiling. Asian Australas J Anim Sci. 18, 1727–1734. doi: 10.5713/ajas.2005.1727

CrossRef Full Text | Google Scholar

Shao, T., Zhanga, L., Shimojo, M., and Masuda, Y. (2007). Fermentation quality of Italian ryegrass (Lolium multiflorum Lam.) silages treated with encapsulated-glucose, glucose, sorbic acid and pre-fermented juices. Asian Australas. J. Anim. Sci. 20, 1699–1704. doi: 10.5713/ajas.2007.1699

CrossRef Full Text | Google Scholar

Stirling, S., Díaz, J. E., Repetto, J. L., Pla, M., Arroyo, J. M., and Cajarville, C. (2021). Growth stage and ensiling: impact on chemical composition, conservation quality and in situ ruminal degradability of whole-crop oat. J. Sci. Food Agric. doi: 10.1002/jsfa.11619

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. doi: 10.3168/jds.S0022-0302(91)78551-2

CrossRef Full Text | Google Scholar

Wang, C., He, L., Xing, Y., Zhou, W., Yang, F., Chen, X., et al. (2019). Fermentation quality and microbial community of alfalfa and stylo silage mixed with moringa oleifera leaves. Bioresour. Technol. 284, 240–247. doi: 10.1016/j.biortech.2019.03.129

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Sun, Y., Zhao, J., Dong, Z., Li, J., Nazar, M., et al. (2020). Assessment of inoculating various epiphytic microbiota on fermentative profile and microbial community dynamics in sterile Italian ryegrass. J. Appl. Microbiol. 129, 509–520. doi: 10.1111/jam.14636

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., He, L., Xing, Y., Zheng, Y., Zhou, W., Pian, R., et al. (2019). Dynamics of bacterial community and fermentation quality during ensiling of wilted and unwilted moringa oleifera leaf silage with or without lactic acid bacterial inoculants. mSphere 4:e00341–19. doi: 10.1128/mSphere.00341-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Weatherburn, M. W. (1967). Phenol-hypochlorite reaction for determination of ammonia. Anal. Chem. 39, 971–974. doi: 10.1021/ac60252a045

CrossRef Full Text | Google Scholar

Xie, Z. L., Zhang, T. F., Chen, X. Z., Li, G. D., and Zhang, J. G. (2012). Effects of maturity stages on the nutritive composition and silage quality of whole crop wheat. Asian Australas. J. Anim. Sci. 25, 1374–1380. doi: 10.5713/ajas.2012.12084

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, D., Ding, W., Ke, W., Li, F., Zhang, P., and Guo, X. (2019). Modulation of metabolome and bacterial community in whole crop corn silage by inoculating homofermentative lactobacillus plantarum and heterofermentative Lactobacillus buchneri. Front. Microbiol. 9:3299. doi: 10.3389/fmicb.2018.03299

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z. S., He, H. Y., Zhang, S. S., and Kong, J. (2017). Effects of inoculants Lactobacillus brevis and Lactobacillus parafarraginis on the fermentation characteristics and microbial communities of corn stover silage. Sci. Rep. 7:13614. doi: 10.1038/s41598-017-14052-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, X. J., Dong, Z. H., Li, J. F., and Shao, T. (2020). Microbial community dynamics and their contributions to organic acid production during the early stage of the ensiling of Napier grass (Pennisetum purpureum). Grass Forage Sci. 75, 37–44. doi: 10.1111/gfs.12455

CrossRef Full Text | Google Scholar

Zhao, S. S., Yang, F. Y., Wang, Y., Fan, X. M., Feng, C. S., and Wang, Y. P. (2021). Dynamics of fermentation parameters and bacterial community in high-moisture alfalfa silage with or without lactic acid bacteria. Microorganisms 9:1225. doi: 10.3390/microorganisms9061225

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X. L., Liu, J. H., Liu, J. J., Yang, F. Y., Zhu, W. B., Yuan, X. F., et al. (2017). Effect of ensiling and silage additives on biogas production and microbial community dynamics during anaerobic digestion of switchgrass. Bioresour. Technol. 241, 349–359. doi: 10.1016/j.biortech.2017.03.183

PubMed Abstract | CrossRef Full Text | Google Scholar

Zi, X., Li, M., Chen, Y., Lv, R., Zhou, H., and Tang, J. (2021). Effects of citric acid and lactobacillus plantarum on silage quality and bacterial diversity of king grass silage. Front. Microbiol. 12:631096. doi: 10.3389/fmicb.2021.631096

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: king grass, silage, growth period, bacterial community, fermentation quality

Citation: Long S, Li X, Yuan X, Su R, Pan J, Chang Y, Shi M, Cui Z, Huang N and Wang J (2022) The Effect of Early and Delayed Harvest on Dynamics of Fermentation Profile, Chemical Composition, and Bacterial Community of King Grass Silage. Front. Microbiol. 13:864649. doi: 10.3389/fmicb.2022.864649

Received: 28 January 2022; Accepted: 23 February 2022;
Published: 07 April 2022.

Edited by:

Mukesh Kumar Awasthi, Northwest A&F University, China

Reviewed by:

Yuli Yanti, Sebelas Maret University, Indonesia
Yong Sun, Northeast Agricultural University, China

Copyright © 2022 Long, Li, Yuan, Su, Pan, Chang, Shi, Cui, Huang and Wang. 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: Jian Wang, 990807@hainanu.edu.cn

These authors have contributed equally to this work

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.