Zhihui Fu
Lin Sun
ZhiJun Wang
Yichao Liu1
Junfeng Hao
Gentu Ge- 1College of Grassland, Resources and Environment, Key Laboratory of Forage Cultivation, Processing and High Efficient Utilization of Ministry of Agriculture, and Key Laboratory of Grassland Resources, Inner Mongolia Agricultural University, Ministry of Education, Hohhot, China
- 2Inner Mongolia Academy of Agricultural and Animal Husbandry Sciences, Hohhot, China
Introduction: Italian ryegrass is less studied in northern China due to high-quality forage grass has not been fully utilized. Full utilization of high-quality forage grass helps to alleviate the shortage of forage grass in winter and spring season and guarantee stable development of livestock production. Consequently, this study was aimed to evaluate the effects of different regions in northern China on the fermentative products, bacterial community compositions, and metabolic pathways and metabolites of Italian ryegrass silage.
Methods: The Italian ryegrass was harvested from three regions (Ordos-WK; Hohhot-AK; Ulanqab-SYK) and ensiled for 60 days. Single molecule real-time (SMRT) sequencing and ultra-high performance liquid chromatography-mass spectrometry (UHPLC–MS/MS) were used to analyze bacterial communities and metabolites, respectively.
Results: After 60 d of fermentation, the SYK group had the lowest pH (4.67), the highest lactic acid contents (95.02 g/kg DM) and largest lactic acid bacteria populations (6.66 log10 cfu/g FM) among the treatment groups. In addition, the SYK group had the highest abundance of Lactiplantibacillus plantarum (63.98%). In SYK group, isoquinoline alkaloid biosynthesis was the significantly enriched (p < 0.05) and high-impact value (0.0225) metabolic pathway. In AK group, tryptophan metabolism the was the significantly enriched (p < 0.001) and high-impact value (0.1387) metabolic pathway. In WK group, citrate cycle (TCA cycle) was the significantly enriched (p < 0.001) and high-impact value (0.1174) metabolic pathway. Further, Lactiplantibacillus plantarum was positively correlated with cinnamic acid, tetranor 12-HETE, D-Mannitol, (2S)-2-amino-4-methylpentanoic acid L-Leucine, guanine, isoleucyl-aspartate and 3,4-Dihydroxyphenyl propanoate, but negatively correlated with isocitrate and D-mannose.
Discussion: In conclusion, this study can improve our understanding of the ensiling microbiology and metabolomics in different regions to further regulate the fermentation products and promote livestock production.
Introduction
Ensiling has become a common and effective method for long-term conservation of livestock forages. At present, ensiling of forages is widely carried out around the world to continuously provide feeds for ruminants (Wang et al., 2022a). Silage quality largely depends on the characteristics of forage at harvest (Boever et al., 2013), because during the ensiling process, anaerobic microorganisms mainly composed of LAB ferment the available chemical component of the forage and use WSC to produce LA, reduce pH, acidize silage, and inhibit harmful microorganisms to ensure long-term preservation of the silage (Na et al., 2022). The chemical composition of forage is important for silage fermentation (Yin et al., 2022), but the chemical components of forage varies with the environment in which it is grown (Johnson et al., 1999). In addition to the chemical composition, the epiphytic microbiota of the aboveground parts also change with the change of the growing environment (Handley, 2019). However, with the interaction between forage and climate, epiphytic microorganisms gradually become the specific microorganisms of their hosts (Yin et al., 2022). Epiphytic microbial communities are important in ensiling forage during fermentation; composed of a variety of microorganisms, they create a variety of metabolites (Kung et al., 2018).
Italian ryegrass (Lolium multiflorum Lam.) is a globally significant fodder crop. It is currently extensively dispersed throughout temperate parts of the world, and cultivated in Europe, America, and Asia (Parvin et al., 2010). Italian ryegrass is one of the most significant and frequent forages for dairy cows in temperate climates (Lv et al., 2021). As the common forage for grazing ruminants, it is characterized by high forage production, nutritional value, digestibility, and good ensiling qualities—it has especially high levels of soluble and degradable nitrogen and carbohydrates (Stergiadis et al., 2015). Making silage from Italian ryegrass can help with bridging the gap between year-round livestock output and the seasonal imbalance of available feed (Wright et al., 2000). However, at present, as a high-quality forage with good ensiling qualities providing high-quality feed for ruminants, Italian ryegrass has been understudied for Italian ryegrass silage in northern China because it is difficult to overwinter in northern areas. Moreover, studies on microbial fermentation and metabolic mechanisms of Italian ryegrass silage have not yet been reported. Therefore, we selected Ordos City, Hohhot City, Ulanqab City, the three representative forage producing areas and milk source gold belt to test the cultivation of high-quality Italian ryegrass to promote the development of high-quality forage industry and dairy industry.
In recent years, SMRT technology has been used to track changing communities of microbial composition and to determine the dominant species in silage, as well as the correlation of dominant species with fermentation products (Bai et al., 2021; Du et al., 2022a). Previous studies have integrated using 16S rRNA sequencing and metabolomics to investigate silage microbiomes and metabolome patterns in order to better understand the biological mechanisms underpinning silage (Guo et al., 2018; Xu et al., 2019). Moreover, many hitherto unknown compounds can be identified by metabolomics (Guo et al., 2018; Xu et al., 2019). Metabolomic research tools help to reveal the biochemical network mechanisms of the fermentation process and can also be used to guide the regulation and prediction of component changes in the fermentation process. During the fermentation process, a significant variety of metabolites are formed, including diverse amino acids, fatty acids, oligosaccharides, vitamins, peptides, and aromatic compounds, and the kinds of these metabolites are researched (Wang et al., 2021). It is important to study the types, quantities, and influencing factors of these metabolites for the exploitation and scientific evaluation of lactic acid bacteria fermentation (Wang et al., 2021). In addition, the use of correlation-based analysis of coupled microbiome-metabolome data sets is becoming increasingly common in research, with the goal of identifying microbial drivers of metabolic variance (Noecker et al., 2019).
To the best of our knowledge, few studies have used multi-omics methods to analyze differences in Italian ryegrass silage fermentation quality, microbial community and metabolome characteristics in different regional settings, and few previous studies have used LC–MS to identify the metabolome of Italian ryegrass silage. We hypothesized that Italian ryegrass grown in different regional environments has different microbial community changes and metabolomics characteristics during fermentation due to differences in chemical composition and microbial composition, resulting in different fermentation quality. Therefore, this study aimed to use SMRT sequencing and metabolomics techniques to study the effects of regions on silage quality from a microbiological and metabolomics perspective. Further, we sought to investigate the microbial community structure and metabolomics characteristics of Italian ryegrass silage fermentation in different regions, so as to better use the regulation of microbial community structure and metabolomics to regulate the fermentation quality to produce high-quality silage. In addition, understanding the relationship between metabolites and fermentation bacteria provides new ideas for the screening and utilization of functional lactic acid bacteria, and provides a theoretical basis for improving fermentation quality.
Materials and methods
Experimental design and silage preparation
The planting time, sampling time, and the basic general situation of Italian ryegrass are shown in Table 1. The cultivation and management measures in the three regions are consistent, including irrigation and fertilization. Italian ryegrasses were harvested at booting stages from three plots randomly selected from three regions (WK, AK, SYK) as three replicates, and the collected samples were placed on clean plastic sheets and left to sun-dry. Sun-drying continued until dry matter content was about 35%, whereupon the Italian ryegrass was chopped into 2-cm-long segments with a forage cutter (Mode-8200, Minghong Business, Shandong, China). From each repetition, 300-g samples were collected in sterilized bags (75% alcohol for sterilization) and plated in an ice-boxes, and then immediately sent to the laboratory for determination of microbial community and chemical characteristics (Table 2). The rest of the prepared ryegrass (500 g) was packed into polyethylene plastic bags (size: 300 mm × 400 mm; Embossed Food Saver Bag Co., Ltd., Chengdu, China) and vacuum sealed with a vacuum sealer (DZ-400, Shandong Zhucheng Yizhong Machinery Co., Ltd., Zhucheng, China). There were three replicates for each treatment, and all samples (three treatments × three replicates) were stored at ambient temperature (24–26°C) under sheltered conditions. After 60 days of ensiling, the fermentation quality, microbial community and metabolites of samples from three regions were determined.
Table 1. Sampling time and the basic general situation of sample plots.
Table 2. Chemical and microbial compositions of fresh materials.
Chemical analysis of raw materials and Italian ryegrass silages
Measurement of each key variable (i.e., chemical composition, fermentation characteristics, and microbial counts) of each sample of raw materials and silage was repeated with three replicates. The determination of metabolomics was performed with six replicates. After ensiling, each sample was divided into two equal parts (three replicates × 2 = six replicates). From each sample, 20 g was taken and placed in a sterilized cryopreservation tube. The cryopreservation tube was immediately placed in liquid nitrogen for quick freezing for 15 min, and then stored at −80°C. These samples (three treatments × 6 replicates = 27) were sent to Majorbio Bio-Pharm Technology (Majorbio Bio-Pharm Technology Co. Ltd., Shanghai, China) for metabolomics determination.
Following the method detailed by Du et al. (2021a), the raw materials and Italian ryegrass silages were dried in an oven for 48 h at 65°C until a constant mass was attained, and then ground to pass a 1-mm screen (FM100, Taisite Instrument Co. Ltd., Tianjin, China) for chemical analysis. TN was determined using a Kjeldahl apparatus (Vapodest 50 s, Gerhart, Bonn, Germany) following Patrica (1997), CP content was determined as the TN × 6.25. NDF and ADF were measured according to the method described by Van Soest et al. (1991), using an Ankom A2000i fiber analyzer (A2000i, Ankom Technology, Macedon, NY, United States). The WSC content was analyzed by referring to the anthrone reagent colorimetry Thomas (1977).
In order to analyze the fermentation characteristics of forage, 10 g samples of silage were combined with 90 g of deionized water and kept in a 4°C fridge for 24 h (Du et al., 2021a). The liquid extract was filtered using four layers of gauze and filter paper and used for the following analysis. The pH of the sample was measured with a glass electrode pH meter (pH S-3C, LEICI, Shanghai, China). The organic acid, including LA, AA, PA, and BA content of silage was determined using high performance liquid chromatography (Model: e2695, Waters Corporation, Massachusetts, United States) as described by Cheng et al. (2021). The AN was determined by the phenol-hypochlorite reaction (Broderick and Kang, 1980).
Following Li et al. (2014), 10 g of fresh materials or Italian ryegrass silage was sampled and blended with 90 ml of sterilized water and serially diluted in sterilized distilled water from 10−1 to 10−5, the numbers of LAB, coliform bacteria, yeast, and molds were then measured. The colonies of LAB were anaerobically cultured on De Man Rogosa Sharpe agar medium at 37°C for 48 h, and the colonies were counted. Coliform bacteria were aerobically incubated at 37°C for 48 h on violet-red bile agar culture media. Yeasts and molds were aerobically incubated at 30°C for 48 h on potato dextrose agar culture media. All culture media were from the same manufacturer (Guangzhou Huankai Microbial Science and Technology Co. Ltd., Guangzhou, China).
Microbial diversity analysis
The Fast DNA SPIN for Soil kit (MP Biomedicals, Solon, OH, United States) was used to extract genomic DNA from the microbial communities of the fresh and silage ryegrass samples. The concentration and purity of the DNA were evaluated using a Nanodrop 2000 UV–Vis spectrophotometer (Thermo Scientific, Wilmington, DE, United States). The primer pairs 27F and 1492R across the full-length bacterial 16S rRNA gene were amplified using an ABI GeneAmp® 9700 PCR thermocycler (ABI, CA, United States).
The PCR amplification procedure was performed following the methods of Li et al. (2022). The resultant PCR amplicons were extracted from a 2% (w/v) agarose geland further purified using AMPure® PB beads (Pacifc Biosciences, CA, United States). The purified product was mixed in equal molars and a DNA library was constructed using SMRTbell® Express Template Prep Kit 2.0 (Pacifc Biosciences, CA, United States) according to the manufacturer’s instructions. The purified SMRTbell library was sequenced using SMRT sequencing technology on the Pacbio Sequel II System (Pacifc Biosciences, CA, United States). The original 16S rRNA gene sequencing data were spliced using FLASH. UPARSE was used to cluster OTUs with a similarity cutoff of 97%. Each OTU representative sequence was analyzed using the RDP Classifier (Version 2.2) classifier algorithm and database (nt_v20210917), and the confidence threshold was 70%.
Bioinformatic analysis of the silage microbial diversity was carried out using the Majorbio Cloud platform.1 Following the methods detailed by Schloss et al. (2009), alpha diversity indices including observed OTUs, Chao1 richness, Shannon index, and Good’s coverage were calculated with Mothur (Version 1.30.1) based on the OTUs information. Following OTU clustering, the Venn diagram package (Version 1.2) in R statistical tools was used to draw a Venn diagram (Du et al., 2021b). The similarity among the microbial communities in different samples was determined by PCoA based on Bray–Curtis dissimilarity using Vegan package (Version 2.5–3). The heat map hierarchical cluster analysis and analysis of similarities (ANOSIM) were shown by R statistical pvclust (Version 3.0.2) (Du et al., 2021b). The distance-based redundancy analysis by Vegan package (Version 2.5–3) was used to examine the relationship between environmental factors and nutrient composition, microbial community structure of the raw materials in three regions (Yang X. et al., 2022).
Sequencing and analysis of metabolites
As described by Zhang et al. (2019), 50 mg of solid silage ryegrass samples were put in an EP tube and the metabolites were extracted with 400 μl methanol: water (4: 1, v/v) solution. According to the method described by Fu et al. (2022), a mixed quality control sample (QC) was prepared. The processing and testing of QC samples were the same as those of analytical samples. The instrument platform of this LC–MS study used a Thermo Fisher ultra-high performance liquid chromatography tandem Fourier transform mass spectrometry UHPLC-quadrupole-electrostatic field orbital trap high resolution mass spectrometer HF-X system. UHPLC–MS/MS analysis were performed according to the method described by Yang Z. et al. (2022).
Multivariate statistical analysis of metabolites and identification of differential metabolites
The ropls R package (Version1.6.2) from Bioconductor on Majorbio Cloud Platform (see Footnote 1) was used to perform multivariate statistical analysis, including PCA, PLSDA, and OPLS-DA. The differential metabolites were detected based on OPLS-DA analysis to obtain the first principal component (VIP) of variable importance in the projection. Combined with Student t test, the screening conditions for the main differential metabolites were: p < 0.05 and VIP > 1. Using one-way analysis of variance (ANOVA) (p < 0.05 for significant difference) and Tukey–Kramer’s test method yielded 471 differential metabolites to create metabolic sets for silage samples from three regions.
Metabolic enrichment and pathway analyses were built up and connected to biochemical pathways using database searches (KEGG).2 The KEGG is a knowledge base for systematic analysis of gene function, linking genomic information and functional information. Using the KEGG database, the metabolites in the metabolic concentration can be classified according to the pathways they participate in or the functions they perform. Enrichment analysis is frequently used to analyze a group of metabolites in a function node whether it appears or not. The principle was that the annotation analysis of a single metabolite develops into an annotation analysis of a group of metabolites. KEGG topology analysis was used as the method for relative-betweenness centrality. Scipy.stats (Python packages)3 was used to identify statistically significantly enriched pathways using a Fisher’s exact test.
Statistical analysis
The variance analysis (ANOVA) of fermentation, nutritional characteristics and microbial populations of fresh materials and silage samples were performed using the general linear model program (GLM) to test the differences between samples in SAS software (ver. 9.3; SAS Institute Inc., Cary, NC, United States). One-way analysis of variance (ANOVA) and a Tukey’s honest significant difference (HSD) test were used to examine differences among samples, and the effect was considered significance when p < 0.05. Microbiota and metabolome data analysis were performed using the Majorbio I-Sanger Cloud Platform.4
Results
Chemical composition and microbial population of raw materials
The chemical composition and microbial populations of fresh materials before ensiling are shown in Table 2. The different growth regions had a significant effect (p < 0.05) on contents of DM, CP, NDF, ADF, and WSC, as well as coliform bacteria, yeasts and molds numbers. The DM content of fresh materials ranged from 33.4 to 35.47%, and the CP content ranged from 14.65% DM to 23.54% DM. The CP content (23.54% DM) of SYK was significantly higher than the others (p < 0.05). The WSC content (9.23% DM) of SYK was the highest and significantly higher than that of WK (p < 0.05). There were no significant differences with LAB among the regions, but the LAB numbers of AK were the highest. The counts of undesirable microorganism (e.g., yeast and molds) in WK were the lowest among three groups.
To analyze the effect of environmental factors on nutrient composition, microbial community structure of raw materials, redundancy analysis was performed for the microbial communities of fresh materials. Mean relative humidity, WSC, CP, altitude, mean temperature, and mean sunshine duration had a significant effect on the variability of epiphytic bacterial communities. Among them, mean temperature, altitude and mean relative humidity were positively correlated with CP and WSC, and explained up to 84.55% (57.99 and 26.56% on RDA axes 1 and 2, respectively) of total variance of bacterial community (Figure 1).
Figure 1. Redundancy analysis the relationship between environmental factors and nutrient composition, microbial community structure of raw materials. The axes are labeled with the percentage of total variance explained (%). Red arrow lengths indicate the variance explained by environmental factors and nutrient composition, and the blue arrow lengths indicate the variance explained with the microbial community structure. Different graphics represent raw material samples of three regions. MT, mean temperature; A, altitude; MRH, mean relative humidity; MSD, mean sunshine duration; CP, crude protein; WSC, water-soluble carbohydrate; NDF, neutral detergent fiber; ADF, acid detergent fiber; WK, China, Inner Mongolia, Ordos-wushenqi; AK, China, Inner Mongolia, Hohhot; SYK, China, Inner Mongolia, Ulanqab-siziwangqi.
Fermentation quality of Italian ryegrass silage in different regions
As shown in Table 3, the regions had a significant effect (p < 0.05) on the AN/TN, and the contents of LA, AA, PA, and BA. But we found that there were no significant differences between the numbers of LAB and WSC content after ensiling. The pH of the silages ranged from 4.67 to 5.11, and the SYK achieved lower pH (4.67) as compared to others. Moreover, SYK also had the highest LA content (95.02 g/kg DM). The WK had the highest pH (5.11), while it had the lowest LA content (21.27 g/kg DM) and was significantly lower than the other groups (p < 0.05). The SYK had the highest AA (7.77 g/kg DM) and PA (13.35 g/kg DM) content and was significantly higher than the other groups (p < 0.05). The LAB population of all three groups increased (4.84 vs. 6.50; 5.06 vs. 5.76; 3.67 vs. 6.66) after 60 days of ensiling, whereas the WSC contents (7.50 vs. 5.25% DM; 9.01 vs. 5.63% DM; 9.23 vs. 5.12% DM) decreased.
Table 3. Fermentation quality of ryegrass silage in different regions.
Microbial community of raw materials and Italian ryegrass silages
The diversity and richness indices of bacterial communities in raw materials and Italian ryegrass silages are illustrated in Table 4. All of the samples had coverage values above 0.99. The Sobs, Shannon, Ace, and Chao 1 indices in all silage samples were lower than that of all fresh materials. After ensiling, the Sobs, Shannon, Ace, and Chao 1 indices ranged from 31.33 to 47.00, 1.83 to 2.53, 46.01 to 121.75, 44.75 to 85.44, respectively. In addition, the highest values for these indices of bacterial diversity and community were observed in fresh materials of group AK.
Table 4. Alpha diversity of the bacterial community in fresh materials and ryegrass silage.
A Venn diagram of the OTUs in raw materials and Italian ryegrass silages is shown in Figure 2. The dominant microbiome of the W, A, and SY (Figure 2A) contained 99 shared OTUs, as well as 7, 16, and 6 unique OTUs, respectively. After ensiling, the WK, AK and SYK silages shared 28 OTUs (Figure 2B), as well as 8, 26, and 13 unique OTUs, respectively. Compared to raw materials, the dominant microbiome of Italian ryegrass silages shared OTUs decreased, and this trend was the same as the Sobs, Shannon, Ace, and Chao 1 indices, indicating that the microbial diversity was decreased after ensiling, while the different regions significantly affected microbial diversity.
Figure 2. Venn analysis of shared and unique microbial OTUs of raw materials and ryegrass silages. (A) Raw materials; (B) Ryegrass silages. W, Raw materials of Ordos-wushenqi; A, Raw materials of Hohhot; SY, Raw materials of Ulanqab-siziwangqi; WK, Ryegrass silages of Ordos-wushenqi; AK, Ryegrass silages of Hohhot; SYK, Ryegrass silages of Ulanqab-siziwangqi.
We used β-diversity analysis utilizing PCoA to confirm differences in the bacterial communities of the fresh materials and Italian ryegrass silage (Figure 3). The PCoA was based on Bray–Curtis distances at the OTU level and an ANOSIM test with 999 permutations between the different groups. As shown in Figure 3, we found that there was no significant separation between bacterial communities of raw materials, but after ensiling we observed a clearly separation among the bacterial communities in different regions, showing that different regions significantly affected the β diversity of microbial communities.
Figure 3. The principal coordinate analysis (PCoA) plot showing variation in bacterial community structure of three regions. W, Raw materials of Ordos-wushenqi; A, Raw materials of Hohhot; SY, Raw materials of Ulanqab-siziwangqi; WK, Ryegrass silages of Ordos-wushenqi; AK, Ryegrass silages of Hohhot; SYK, Ryegrass silages of Ulanqab-siziwangqi.
Relative bacterial abundances at the phylum (Figure 4A), genus (Figure 4B) and species (Figure 4C) level in raw materials and Italian ryegrass silages are shown in Figure 4. Proteobacteria was the main phylum in the fresh materials, while Firmicutes were the main phylum in Italian ryegrass silages. The relative abundance of Proteobacteria in W, A and SY were 44.57, 50.8, and 46.96%, respectively. After ensiling, the abundance of Proteobacteria decreased while Firmicutes abundance increased. The relative abundance of Firmicutes in WK, AK, and SYK were 73.42, 95.91, and 93.44%, respectively. Unclassified bacteria, Pantoea, Microbacterium, and Paracoccus were the main genera in raw materials. After fermentation, Enterococcus, Lactiplantibacillus, and Levilactobacillus were the most abundant bacteria. The Enterococcus relative abundance of WK, AK, and SYK were 43.17, 40.88, and 11.24%, respectively. The Lactiplantibacillus relative abundance of WK, AK, and SYK were 4.17, 26.24, and 63.98%, respectively. The Levilactobacillus relative abundance of WK, AK, and SYK were 19.28, 6%, and 14.28%, respectively. The highest abundance of Lactiplantibacillus occurred in SYK (63.98%) while and the highest abundance of Enterococcus occurred in WK and AK (43.17 and 40.88%, respectively). After ensiling, uncultured Enterococcus sp. (34.37%) and uncultured bacterium (23.68%) were the predominant species in WK, while the relative abundance of Lactiplantibacillus plantarum was only 4.17%. The relative abundance of Lactiplantibacillus plantarum and uncultured Enterococcus sp. in AK were 26.23 and 27.76%, respectively. The relative abundance of Lactiplantibacillus plantarum and uncultured Enterococcus sp. in SYK were 63.98 and 10.52%, respectively.
Figure 4. Bacterial community composition on phylum (A), genus (B), and species (C) level in fresh ryegrass and ryegrass silage by SMRT. W, Raw materials of Ordos-wushenqi; A, Raw materials of Hohhot; SY, Raw materials of Ulanqab-siziwangqi; WK, Ryegrass silages of Ordos-wushenqi; AK, Ryegrass silages of Hohhot; SYK, Ryegrass silages of Ulanqab-siziwangqi.
Correlation analysis of fermentation products and microbial community
As shown in Figure 5, the pH was negatively correlated with Lactiplantibacillus plantarum, Latilactobacillus sakei, and Weissella paramesenteroides, but positively correlated with Enterococcus faecalis and uncultured Enterococcus sp. The AN/TN was positively correlated with Lactiplantibacillus plantarum, Latilactobacillus sakei, and Weissella paramesenteroides, but negatively correlated with Enterococcus faecalis, Paracoccus sp., and Loigolactobacillus coryniformis. LA, AA and PA were positively correlated with Lactiplantibacillus plantarum, Latilactobacillus sakei, and Weissella paramesenteroides, but negatively correlated with uncultured Enterococcus sp., Enterococcus faecalis, and Paracoccus sp. BA was positively correlated with Lactococcus sp., Weissella minor, and Fructilactobacillus lindneri, but negatively correlated with uncultured bacterium. LAB was negatively correlated with Agrobacterium rubi, Alkalihalobacillus gibsonii, Fructilactobacillus lindneri, and Enterococcus hermanniensis.
Figure 5. Correlation analyses between bacterial community and terminal fermentation products at species level. Fermentation characteristics are displayed horizontally and the bacterial community information is displayed vertically. The corresponding value of the middle heat map is the Spearman correlation coefficient r, which ranges between −1 and 1; r < 0 indicates a negative correlation (orange), r > 0 indicates a positive correction (green), and “*,” “**,” and “***” represent p < 0.05, p < 0.01, and p < 0.001, respectively. AN/TN, ammonia nitrogen/total nitrogen; LAB, lactic acid bacteria; WSC, water-soluble carbohydrate.
KEGG enrichment analysis and topology analysis of metabolite pathway identified in Italian ryegrass silage
KEGG enrichment analysis of differential metabolites after ensiling from three regions in Figure 6. The ordinate in the figure represents the pathway name, and the abscissa represents the enrichment rate, which indicates the ratio of the numbers of enriched metabolites in the pathway to the background number annotated to the pathway. The larger the value, the greater the degree of enrichment. In SYK group, only isoquinoline alkaloid biosynthesis metabolic pathway was significantly enriched (p < 0.05), and the enrichment rate was 0.0164. In AK group, there were 4 significantly enriched pathways. They were betalain biosynthesis, tryptophan metabolism, tyrosine metabolism and isoquinoline alkaloid biosynthesis, among which tyrosine metabolism were significantly enriched (p < 0.001), and betalain biosynthesis had the highest enrichment rate (0.0417). In WK group, there were 7 significantly enriched pathways. There were nitrogen metabolism, amino sugar and nucleotide sugar metabolism, pyrimidine metabolism, ABC transporters, alanine, aspartate and glutamate metabolism, citrate cycle (TCA cycle) and glyoxylate and dicarboxylate metabolism, among which citrate cycle (TCA cycle) were significantly enriched (p < 0.001), and it had the highest enrichment rate (0.1). According to the classification of KEGG pathway database, except for ABC transporters, which are related to environmental information processing, the remaining metabolites belong to the metabolic class.
Figure 6. KEGG pathway enrichment analysis of the main differential metabolites of ryegrass silage from three different regions in positive and negative ion modes (WK, Ordos; AK, Hohhot; SYK, Ulanqab). M, and EIP are the class names of metabolic pathways in KEGG annotation. M, metabolism; EIP, environmental information processing. Value of p-uncorrected <0.05 and column chart showing value of p values for the top 20 pathways; *p < 0.05; **p < 0.01, ***p < 0.001.
Topology analysis of differential metabolites after ensiling from three regions as shown in Figures 7A–C, and KEGG topology statistics table was shown in Supplementary Table 1. Each bubble in the figure represents a KEGG pathway; the horizontal axis represents the relative importance of the metabolites in the pathway, the size of impact value; the vertical axis represents the enrichment significance of metabolites involved in the pathway-log10 (p-value); the size of the bubble represents the impact value; the larger the bubble, the greater the importance of the pathway. Bubbles labeled with the letter A in the figure represent significantly enriched (p < 0.05) and high-impact value metabolic pathways. In SYK group (Figure 7A), the significantly enriched and high-impact value (0.0225) metabolic pathway was isoquinoline alkaloid biosynthesis. In AK group (Figure 7B), the significantly enriched and high-impact value (0.1387) metabolic pathway was tryptophan metabolism. In WK group (Figure 7C), the significantly enriched and high-impact value (0.1174) metabolic pathway was citrate cycle (TCA cycle). This was consistent with the results of KEGG enrichment analysis.
Figure 7. (A) SYK, Ryegrass silages of Ulanqab-siziwangqi, (B) AK, Ryegrass silages of Hohhot and (C) WK, Ryegrass silages of Ordos-wushenqi; Bubbles labeled with the letter A in the figure represent significantly enriched (p < 0.05) and high-impact value metabolic pathways.
Correlation analysis of high abundance bacteria and fermentation metabolites
Spearman correlation between bacteria and differential fermentation metabolites at the level of species is shown in Figure 8, and those bacteria had the highly relative abundance in three region silages. Lactiplantibacillus plantarum was positively correlated with cinnamic acid, tetranor 12-HETE, D-Mannitol, (2S)-2-amino-4-methylpentanoic acid L-Leucine, guanine, isoleucyl-aspartate and 3,4-Dihydroxyphenyl propanoate, but negatively correlated with isocitrate and D-mannose. Levilactobacillus brevis was negatively correlated with L-phenylalanyl-L-proline. The Enterococcus hermanniensis was positively correlated with chlorogenic acid, 2-hydroxycinnamic acid, ganoderic acid F, methylmalonic acid, and isorhamnetin 3-glucoside.
Figure 8. Correlation analysis of high abundance of species-level bacteria and metabolites in silage from three regions. Bacterial community information is displayed horizontally and the fermentation metabolites are displayed vertically. The corresponding value of the middle heat map is the Spearman correlation coefficient r, which ranges between −1 and 1; r < 0 indicates a negative correlation (orange), r > 0 indicates a positive correction (green), and “*,” “**,” and “***” represent p < 0.05, p < 0.01, and p < 0.001, respectively.
Discussion
Ensiling involves complex interactions between microorganisms and chemical components of raw materials in an anaerobic environment (Wang et al., 2022b). Therefore, the involved species of microorganisms attached to the surface of the raw materials and chemical components of the raw materials are important factors affecting silage quality (Yin et al., 2022). To investigate the contribution of microorganisms and chemical components attached to the raw materials of Italian ryegrass to the fermentation products of Italian ryegrass silage in three different regions, the microbial community characteristics, metabolic pathways and metabolites of Italian ryegrass silages after 60 days ensiling were analyzed.
Chemical composition and microbial population of raw materials
As a high-quality forage, Italian ryegrass has a good nutritional profile, including high levels of CP and WSC (Yan et al., 2019). In this study, the CP content was higher than 14% DM in all regions, where the CP content of the group SYK was the highest, which was higher than that of perennial ryegrass (10.46% DM, Dong et al., 2020), alfalfa (18.3% DM, Wang et al., 2019), and dried corn stover (6.08% DM, Yan et al., 2019). Guan et al. (2018) reported that differences in temperature, precipitation, and soil environment may influence the CP content. The WSC substance of raw materials is basic for lactic acid fermentation, and the ideal WSC substance ought to surpass 5% DM for a satisfactory silage fermentation (Li et al., 2019). The WSC substance (7.5% DM, 9.01% DM, and 9.23% DM) of fresh Italian ryegrass in the three different regions was higher than 5% DM, thereby demonstrating that adequate substrate may be given for microorganisms. Among these, the CP content (23.54% DM) and WSC content (9.23% DM) in the SYK group were the highest, which may because the SYK region had the highest altitude of the studied regions (Ding et al., 2020; Yang X. et al., 2022). Ding et al. (2020) reported that with an increase in altitude, WSC and CP content also increased in fresh grass. This may because low temperatures at the high altitudes restrain plant respiration, which is conducive to the accumulation of WSC and CP in the cell protoplasm. Wang et al. (2017) suggested that the minimum requirement for LAB number on raw material should be more than 5.00 log10 cfu/g FM. In our study, the low number of LAB attached to the raw material may lead to poor silage fermentation.
Fermentation quality of Italian ryegrass silage in different regions
After 60 days of ensiling, the LAB population of all three groups increased, whereas the WSC contents decreased, suggesting that fermentation of LAB consumes WSC during ensiling. According to Muck et al. (2018), LAB used the WSC as a substrate during fermentation to produced LA and preserve forage during silage. Therefore, the WSC content in each treatment group decreased after ensiling, but all were higher than 5% DM. According to Wang et al. (2022b), higher WSC content in silage may be more beneficial for ruminant nutrition, which also proves that Italian ryegrass silage is a high-quality feed source for ruminants. According to a previous study, a pH of 4.2 or less indicates well-fermented silage (Xiong et al., 2022). In this study, among all treatment groups, the SYK group had the lowest pH at 4.67, which is higher than 4.2. According to previous study (Yan et al., 2019), the pH value of natural fermented Italian ryegrass silage without additives is higher than 4.2, which agrees with our findings. The final pH of silage is affected by several factors (e.g., very high protein and ash content) which may account for silages that present with higher than the normal pH values (Kung et al., 2018). In the SYK group, the CP content of the raw material and the AN/TN ratio after ensiling were higher than others, which indicated that the CP was degraded more during the fermentation process, and the amino acid decarboxylation produced basic amines, resulting in a higher pH value (Pessione et al., 2010). After 60 days of ensiling, the SYK group had the highest LA content, which was possibility due to the relatively high abundance of Lactiplantibacillus plantarum (63.98%) and some Levilactobacillus brevis (14.28%) in the SYK group. It is reported that homofermentative LAB (e.g., Lactiplantibacillus plantarum) produce LA to reduce pH during fermentation, while heterofermentative LAB (e.g., Levilactobacillus brevis) produce LA and AA during fermentation to improve aerobic stability (Xu et al., 2021). Silage samples containing heterofermentative LAB had a higher pH and AA concentrations (Xu et al., 2019). However, there was no significant difference in the numbers of LAB in the three different regions, indicating that the fermentation efficiency of LAB in the SYK silage samples may be greater in later fermentation stage (Ding et al., 2020). In our future studies, we can further study the acid-producing ability by screening and purifying the dominant LAB in samples from SYK region. We found that AA and PA content in the SYK group was the highest compared to the other groups, which may be related to the presence of more bacteria producing AA and PA in these raw materials (Guan et al., 2018). Propionibacteria that convert glucose and LA into PA and AA may also exist in SYK silage samples (Kung et al., 2018). The WK group had the highest pH and the lowest LA content, which was due to the lowest relative abundance of Lactiplantibacillus plantarum (4.17%) in the WK group (Xu et al., 2021). In addition, compared with the WSC content of the raw materials, the SYK group consumed the most WSC (9.23% DM vs. 5.32% DM) during the fermentation process, and the WK group (7.50% DM vs. 5.25% DM) consumed the least, indicating that the growth and reproduction of Lactiplantibacillus plantarum consumed more WSC content (Muck et al., 2018). The ammonia nitrogen ratio is an important indicator of silage fermentation quality; the greater the ratio, the more amino acids and proteins are decomposed (Denek et al., 2016). In the present study, the ammonia nitrogen ratios in all treatments were less than 10%, indicating that the protein degradation of silage was acceptable in all regions (He et al., 2022).
Microbial community of raw materials and Italian ryegrass silages
Since the ensiling process depends on the interactions of multiple bacteria, the bacterial community structure directly affects the silage quality (Ni et al., 2017). In the present study, the bacteria in all samples were sequenced by SMRT sequencing technology to accurately assess the microbial community and diversity (Table 4). The coverage values were greater than 0.99 in all samples, indicating that the reliability of species detection was very high and the sequencing depth was sufficient to detect the maximum bacterial communities. The alpha diversity indices (ACE, Chao1, and Shannon) were used to reflect the microbial richness and species diversity of samples (Du et al., 2021b). Shannon’s index is used to measure species diversity with a higher value indicating increased species diversity. Ace and Chao 1 indices are used to measure species richness with lower values indicating lower species richness. In our study, the alpha diversity indices (ACE, Chao1, and Shannon) of the fresh material was higher than that of the Italian ryegrass silage samples, and these indices showed trends similar to the Sobs index. This suggests that the microbial richness and species diversity of Italian ryegrass decreased after silage (Cai et al., 2020). This may be because fresh material exists in an aerobic and neutral environment that suits the proliferation of epiphytic aerobic microorganisms (Xu et al., 2021). However, after 60 days of silage, an acidic anaerobic environment was formed in the silage samples, resulting in a decrease in bacterial diversity, which was also confirmed by the lowest pH values and the Shannon index in SYK. This also suggested that the acidic and anaerobic environment affect the succession of microorganisms in the silage (Wang et al., 2022a). Du et al. (2021b) also made the same observations and found that the lowest alpha diversity was observed in high-quality silage due to the large proportion of beneficial microorganisms. Therefore, the alpha diversity index can be used as a reference for evaluating high-quality silage.
Venn diagrams (Figure 2) showed common microbes in different environments for all samples, and the number of unique and shared microbes in different samples based on OTU taxa (Guan et al., 2018). Therefore, the amount of raw material OTUs were different in the three studied regions (Figure 2A). After 60 day of fermentation, the OTUs shared by regions of silage decreased from 99 to 28, and this was due to the formation of an acidic environment by microbial fermentation, which had been the shared OTU of the dominant microbiota, leading to a decrease in the microbial community (Du et al., 2021b).
The PCoA map clearly illustrated the changes of the microbial community before and after ensiling. The microbial community of the raw materials from the three regions was not clearly separated. The distance of SY and A was closer than W, likely because of the similar geographic location and altitude between group AK and SYK, while the geographic location and altitude of the WK group and these two groups were further apart, indicating that the geographic environment, altitude, and climate differences in three different regions affected the microbial composition (Yang X. et al., 2022). The ensiling samples were significantly separated from the raw material samples, indicating that the microbial community was changed during the silage process, which affected the silage fermentation products and metabolic differences (Ni et al., 2017).
The microbial community composition, structure, and function are key to the process of ensiling fermentation (Cai et al., 2020; Du et al., 2022a). In our study, Proteobacteria was the main phyla in Italian ryegrass raw materials in three regions, and this result was same as Wang et al. (2022b). Proteobacteria play an important role in degrading organic matter and accelerating carbon and nitrogen cycles (Ma et al., 2018). After 60 day of fermentation, the relative abundance of Proteobacteria decreased, while Firmicutes increased to become the most dominant phylum. This may be due to the fact that the acidic and anaerobic environment during silage fermentation inhibited the growth and reproduction of Proteobacteria, but promoted the growth and reproduction of Firmicutes (Keshri et al., 2018). According to Wang et al. (2020), the acid hydrolysis function of Firmicutes plays an important role in anaerobic environment, it can produce a variety of enzymes (e.g., proteases).
In the raw material, in addition to unidentified bacteria, we also found several bacteria, such as Agrobacterium, Pantoea, and Paracoccus. According to Gao et al. (2021), imidacloprid biodegradation in soils was aided by Paracoccus bacteria and it was positively correlated with high degradation activity. According to Walterson and Stavrinides (2015), Pantoea was isolated from Enterobacter genus and fermented lactose while forming slime colonies, and some Pantoea sp. are pathogenic to vegetables (Du et al., 2022a). However, little research has been done on the role of these bacteria in silage. Ogunade et al. (2018) reported Pantoea had an ability to reduce ammonia-nitrogen concentrations in silage, but more studies are required to elucidate their role during ensiling. Based on previous study (Gao et al., 2021), we speculate that these bacteria may originate from the soil in which the Italian ryegrass raw material lives and are carried from the soil during the growth or harvest, but the role of these bacteria in silage fermentation needs further study.
After 60 day of ensiling, Lactiplantibacillus and Enterococcus were the most dominant genera, but their relative abundances varied among treatments. Previous study (Wang et al., 2019) have found that Lactiplantibacillus was the most common bacteria after successful fermentation of silage, because Lactiplantibacillus inhibits the growth of harmful microorganisms by producing LA and lowering the pH of the silage under anaerobic conditions (Cai et al., 2020). In the current study, the SYK group had the highest relative abundance of Lactiplantibacillus plantarum, followed by the AK group, and WK had the lowest relative abundance of Lactiplantibacillus plantarum. The SYK group had the lowest abundance of uncultured Enterococcus sp., and the WK group had the highest relative abundance of uncultured Enterococcus sp. and uncultured bacterium, which indicated that during the fermentation process of SYK, the attached microorganisms that were conducive to fermentation interacted with the chemical components in the raw materials to promote fermentation, thereby increasing the abundance of LAB. However, LAB in the WK group were not fully fermented, resulting in a lower abundance of Lactiplantibacillus plantarum and a higher abundance of uncultured Enterococcus sp. and uncultured bacterium (Du et al., 2021c). Moreover, the WK group had the highest pH value, followed by the AK group, and the SYK group had the lowest pH value. This was because Lactiplantibacillus plantarum is a facultatively homofermentative LAB that grows well in acidic environments and promotes LA fermentation during silage (Li and Nishino, 2013). Santos et al. (2016) reported that most Enterococcus sp. detected in silage were non-pathogenic bacteria that competed with LAB for nutrients to utilize the WSC as a substrate to survive at low pH environment, resulting in nutrient loss.
Correlation analysis of fermentation products and microbial community
Microorganisms influence silage fermentation through a series of metabolites, and these metabolites also affect the community structure of microorganisms. They interact during the fermentation process to jointly regulate the process of silage fermentation and affect silage quality (Du et al., 2022b). In the current study, we performed correlation analyses between the bacterial community and terminal fermentation products at the species level (Figure 5) by using SMRT. After 60 day of ensiling, Lactiplantibacillus plantarum, Latilactobacillus sakei, and Weissella paramesenteroides were negatively correlated with pH, but positively correlated with AN/TN, LA, AA, and PA. Enterococcus faecalis and uncultured Enterococcus sp. was positively correlated with pH, but negatively correlated with LA, AA, and PA. The results are similar to those of Du et al. (2022b), indicating that acid tolerance of Lactiplantibacillus plantarum, Latilactobacillus sakei, and Weissella paramesenteroides, and their critical importance for LA production during fermentation. Lactiplantibacillus plantarum rapidly acidifies silage to inhibit the growth of harmful microorganisms in the late stage of ensiling (Du et al., 2021b). Latilactobacillus sakei was able to use glucose, fructose and different hexose sugars as primary energy sources to produce LA through the glycolytic pathway (Belleggia et al., 2022). In addition, Latilactobacillus sakei was degraded by proteases and aminopeptidases, using free amino acids and nucleotides as energy sources (Belleggia et al., 2020). Du et al. (2021b) reported that Lactiplantibacillus plantarum can reduce the pH of the silage and limit the loss of protein and carbohydrates, thereby reducing the NH3-N concentration and the fermentation loss in the silages. Generally, LAB fermentation can effectively inhibit protein decomposition, whereas the presence of LAB had the opposite effect in this study. This may be because the initial acidification caused by lactic acid fermentation fails to effectively prevent the proliferation and fermentation of Clostridia in silage, resulting in poor quality silage due to the production of propionic acid and butyric acid or the accumulation of ammonia and amine (Li et al., 2020). Tsigkrimani et al. (2022) reported that Enterococcus faecalis, Weissella paramesenteroides, and Enterococcus sp. were resistant to antibiotics (vancomycin). These microorganisms affect animal health by affecting silage quality, which would be an interesting research topic.
KEGG metabolomic pathways and properties of Italian ryegrass silage
Metabolomics technology can accurately identify metabolites in silage. In recent years, many studies (Guo et al., 2018; Xu et al., 2019, 2021) have applied metabolomics technology to identify metabolites during silage fermentation and the correlation between fermenting bacteria and metabolites. At the same time, the KEGG database was used to analyze the biological function and utility of metabolites from a high-level and genomic perspective (Wang et al., 2022a). This can be achieved by displaying the metabolites in the metabolic concentration on the KEGG pathway map, and using KEGG enrichment analysis to explore the most relevant pathways and underlying mechanisms (Yang Z. et al., 2022).
The ensiling process is mediated by microbial metabolic pathways, which have transformed metabolites or degraded substrates (Bai et al., 2022). In this study, we constructed a metabolic set of the identified differential metabolites, and used KEGG enrichment analysis to clarify that interaction between the microorganisms on the surface of Italian ryegrass and the chemical components of Italian ryegrass in three regions had extremely significant effects on biosynthesis of other secondary metabolites (isoquinoline alkaloid biosynthesis, betalain biosynthesis), amino acid metabolism (alanine, aspartate and glutamate metabolism, tyrosine metabolism, tryptophan metabolism), carbohydrate metabolism (glyoxylate and dicarboxylate metabolism, citrate cycle, amino sugar and nucleotide sugar metabolism), membrane transport (ABC transporters), nucleotide metabolism (pyrimidine metabolism) and energy metabolism (nitrogen metabolism). In addition, the biosynthesis of other secondary metabolites (isoquinoline alkaloid biosynthesis), amino acid metabolism (tyrosine metabolism) and carbohydrate metabolism (citrate cycle) were the most significant and important metabolic pathways in the SYK, AK, and WK groups, respectively.
Du et al. (2022a) reported that glycolysis, protein hydrolysis, carbohydrate metabolism, and amino acid metabolism are the main microbial metabolic pathways affecting the flavor and quality of silage. Amino acids are essential substances for forage grasses and its metabolism can be oxidized to carbon dioxide and water through the tricarboxylic acid cycle, releasing energy, and playing an important role in protein synthesis and primary metabolism of forage grasses (Du et al., 2021b, 2022b). Amino acid metabolism plays a crucial role in the formation of metabolites, most of which are essential metabolites required by LAB for growth and protein synthesis during fermentation (Zhou et al., 2022). Phenylalanine is an essential and aromatic amino acid that is oxidized to tyrosine in the body phenylalanine hydroxylase and is involved in sugar metabolism (Pan et al., 2020). In the citric acid cycle, a condensation reaction between oxaloacetic acid and the acetyl group of acetyl-CoA produces citric acid, which is often used as a flavoring agent in foods, and Citric acid-mediated pH adjustment can improve oxidant performance and enzyme activity, thus extending food shelf life (Ke et al., 2017). Understanding the metabolic pathways and biological functions of metabolites can help us regulate the fermentation process to obtain high-quality silage or silage with certain biological functions.
Correlation analysis of high abundance bacteria and fermentation metabolites
Previous research (Guo et al., 2018; Xu et al., 2019; Guan et al., 2020) found that metabolomics technology can more accurately mirror the metabolites composition of the environment, and that it can also be used in silage. Those studies also highlighted links between amino acid metabolism, coenzyme factors and vitamin metabolism, lipid metabolism, carbohydrate metabolism, and terpenoid and ketone compound metabolism in Napier grass silage inoculated by screened LAB. The above studies have greatly enriched our understanding of metabolites in silage. In this study, we also used SMRT sequencing and metabolomics technology to analyze the correlation of differential metabolites in each treatment group with bacteria with higher relative abundance. The contents of cinnamic acid, tetranor 12-HETE, D-Mannitol, (2S)-2-amino-4-methylpentanoic acid L-Leucine, guanine, isoleucyl-aspartate and 3,4-Dihydroxyphenyl propanoate were positively correlated with the abundance of Lactiplantibacillus plantarum, while the contents of isocitrate and D-mannose were negatively correlated with the abundance of Lactiplantibacillus plantarum. Xu et al. (2021) reported that Lactiplantibacillus plantarum was positively correlated with 3-hydroxy fatty acids and could be considered as a species for screening inoculants with potential antifungal activity. Cinnamic acid and its derivatives are widely found in plants, grains, and vegetables, and in recent years, they have attracted much attention due to the positive health effects of various derivatives of cinnamic acid (such as caffeic acid, ferulic acid), including antioxidant, antibacterial and hepatoprotective effects (Yamaga et al., 2022). Ferulic acid, also a hydroxycinnamic acid, has antioxidant activity and is negatively correlated with Lactiplantibacillus plantarum (Xu et al., 2021). The content of L-phenylalanyl-L-proline was negatively correlated with the abundance of Levilactobacillus brevis. L-phenylalanyl-L-proline is a derivative of L-proline, a major amino acid that maintains cell structure and function, as well as a key regulator of cell metabolism and physiology (Liu et al., 2019). The contents of chlorogenic acid, 2-hydroxycinnamic acid, ganoderic acid F, methylmalonic acid and isorhamnetin 3-glucoside were positively correlated with the abundance of Enterococcus hermanniensis. Chlorogenic acid is a phenylpropionic acid compound with antioxidant, anticancer, antibacterial, antihistamine and other health promoting properties (Zhou et al., 2022). Given that Enterococcus hermanniensis was positively correlated with chlorogenic acid, we may consider using it as a candidate strain to produce LAB with antibacterial properties in future study. An in-depth understanding of the relationship between metabolites and fermenting bacteria will help us to understand the fermentation mechanism of silage from the perspective of bioinformatics, so as to provide a theoretical basis for screening functional lactic acid bacteria and producing functional compound lactic acid bacteria silage additives.
Conclusion
The different regional environments can change the microbial community attached to the surface of Italian ryegrass and the chemical composition of Italian ryegrass, thereby altering the fermentation quality and metabolomic characteristics of silage by affecting the interaction between microorganisms and chemical components. After silage, the SYK samples had the lowest pH value and best fermentation quality. The biosynthesis of other secondary metabolites (isoquinoline alkaloid biosynthesis), amino acid metabolism (tyrosine metabolism) and carbohydrate metabolism (citrate cycle) were the most significant and important metabolic pathways in the SYK, AK, and WK groups, respectively. Cinnamic acid content, which has a positive effect on health, increased with the abundance of Lactiplantibacillus plantarum, and chlorogenic acid and 2-hydroxycinnamic acid, which have antibacterial and antioxidant properties, increased with the abundance of Enterococcus hermanniensis. Therefore, using SMRT and metabolomics technology to conduct multi-omics combined analysis of silage will help us to deeply understand the fermentation mechanism of silage. Furthermore, our findings provide new ideas and information for the screening of functional lactic acid bacteria strains and the production of functional compound additives and high-quality functional silage.
Data availability statement
The datasets generated for this study can be found in the NCBI and an accession number SRP373166.
Author contributions
ZF: methodology, visualization, validation, data curation, and wrote the original draft. LS and CG: interpreted the data and edit the language. ZW: conceptualization, acquisition, and writing – review and editing. YL and JH: software. GG: conceptualization and funding acquisition. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the Key Laboratory of Forage Cultivation and the Processing and Highly Efficient Utilization of the Ministry of Agriculture, the Key Laboratory of Grassland Resources of the Ministry of Education, and funded by the National Technical System of Forage Industry for Dry Grass Storage (CARS-34), China.
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1076499/full#supplementary-material
Abbreviations
LAB, lactic acid bacteria; WSC, water-soluble carbohydrates; LA, lactic acid; SMRT, Single Molecule Real Time; UHPLC–MS/MS, ultra-high performance liquid chromatography-mass spectrometry; FW, fresh matter; DM, dry matter; TN, total nitrogen; AN, ammonia nitrogen; CP, crude protein; NDF, neutral detergent fiber; ADF, acid detergent fiber; AA, acetic acid; PA, propionic acid; BA, butyric acid; OTUs, operational taxonomic units; PCoA, principal coordinate analysis; PCA, principal component analysis; PLSDA, partial least-square discriminant analysis; OPLS-DA, orthogonal projections to latent structure discriminant analysis; KEGG, Kyoto Encyclopedia of Genes and Genomes; WK, Italian ryegrass silage from Ordos (China, Inner Mongolia, Ordos-wushenqi); AK, Italian ryegrass silage from Hohhot (China, Inner Mongolia, Hohhot); SYK, Italian ryegrass silage from Ulanqab (China, Inner Mongolia, Ulanqab-siziwangqi).
Footnotes
References
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
Bai, J., Franco, M., Ding, Z., Hao, L., Ke, W., Wang, M., et al. (2022). Effect of Bacillus amyloliquefaciens and Bacillus subtilis on fermentation, dynamics of bacterial community and their functional shifts of whole-plant corn silage. J. Anim. Sci. Biotechnol. 13:7. doi: 10.1186/s40104-021-00649-0
Belleggia, L., Ferrocino, I., Rita Corvaglia, M., Cesaro, C., Milanovic, V., Cardinali, F., et al. (2022). Profiling of autochthonous microbiota and characterization of the dominant lactic acid bacteria occurring in fermented fish sausages. Food Res. Int. 154:110990. doi: 10.1016/j.foodres.2022.110990
Belleggia, L., Milanović, V., Ferrocino, I., Cocolin, L., Haouet, M. N., Scuota, S., et al. (2020). Is there any still undisclosed biodiversity in Ciauscolo salami? A new glance into the microbiota of an artisan production as revealed by high-throughput sequencing. Meat Sci. 165:108128. doi: 10.1016/j.meatsci.2020.108128
Boever, J. D., Dupon, E., Wambacq, E., and Latre, J. (2013). The effect of a mixture of Lactobacillus strains on silage quality and nutritive value of grass harvested at four groeth stages and ensiled for two periods. Agric. Food Sci. 22, 115–126. doi: 10.23986/afsci.6709
Broderick, G. A., and Kang, J. H. (1980). Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media. J. Dairy Sci. 63, 64–75. doi: 10.3168/jds.S0022-0302(80)82888-8
Cai, Y., Du, Z., Yamasaki, S., Nguluve, D., Tinga, B., Macome, F., et al. (2020). Community of natural lactic acid bacteria and silage fermentation of corn stover and sugarcane tops in Africa. Asian Australas J. Anim. Sci. 33, 1252–1264. doi: 10.5713/ajas.19.0348
Cheng, Q. M., Li, P., Xiao, B. X., Yang, F. Y., Li, D. X., Ge, G. T., et al. (2021). Effects of LAB inoculant and cellulase on the fermentation quality and chemical composition of forage soybean silage prepared with corn stover. Grassl. Sci. 67, 83–90. doi: 10.1111/grs.12289
Denek, N., Aydin, S. S., and Can, A. (2016). The effects of dried pistachio (Pistachio veraL.) by-product addition on corn silage fermentation andin vitromethane production. J. Appl. Anim. Res. 45, 185–189. doi: 10.1080/09712119.2016.1141778
Ding, Z., Bai, J., Xu, D., Li, F., Zhang, Y., and Guo, X. (2020). Microbial community dynamics and natural fermentation profiles of ensiled alpine grass Elymus nutans prepared from different regions of the Qinghai-Tibetan plateau. Front. Microbiol. 11:855. doi: 10.3389/fmicb.2020.00855
Dong, L. F., Zhang, H. S., Gao, Y. H., and Diao, Q. Y. (2020). Dynamic profiles of fermentation characteristics and bacterial community composition of Broussonetia papyrifera ensiled with perennial ryegrass. Bioresour. Technol. 310:123396. doi: 10.1016/j.biortech.2020.123396
Du, Z., Lin, Y., Sun, L., Yang, F., and Cai, Y. (2021c). Microbial community structure, co-occurrence network and fermentation characteristics of woody plant silage. J. Sci. Food Agric. 102, 1193–1204. doi: 10.1002/jsfa.11457
Du, Z., Sun, L., Chen, C., Lin, J., Yang, F., and Cai, Y. (2021a). Exploring microbial community structure and metabolic gene clusters during silage fermentation of paper mulberry, a high-protein woody plant. Anim. Feed Sci. Tech. 275:114766. doi: 10.1016/j.anifeedsci.2020.114766
Du, Z., Sun, L., Lin, Y., Chen, C., Yang, F., and Cai, Y. (2022a). Use of Napier grass and rice straw hay as exogenous additive improves microbial community and fermentation quality of paper mulberry silage. Anim. Feed Sci. Tech. 285:115219. doi: 10.1016/j.anifeedsci.2022.115219
Du, Z., Sun, L., Lin, Y., Yang, F., and Cai, Y. (2021b). The use of PacBio SMRT technology to explore the microbial network and fermentation characteristics of woody silage prepared with exogenous carbohydrate additives. J. Appl. Microbiol. 131, 2193–2211. doi: 10.1111/jam.15124
Du, Z., Yamasaki, S., Oya, T., Nguluve, D., Euridse, D., Tinga, B., et al. (2022b). Microbial co-occurrence network and fermentation information of natural Woody-Plant silage prepared with grass and crop by-product in Southern Africa. Front. Microbiol. 13:756209. doi: 10.3389/fmicb.2022.756209
Fu, Z., Sun, L., Hou, M., Hao, J., Lu, Q., Liu, T., et al. (2022). Effects of different harvest frequencies on microbial community and metabolomic properties of annual ryegrass silage. Front. Microbiol. 13:971449. doi: 10.3389/fmicb.2022.971449
Gao, Y., Liu, M., Zhao, X., Zhang, X., and Zhou, F. (2021). Paracoccus and Achromobacter bacteria contribute to rapid biodegradation of imidacloprid in soils. Ecotoxicol. Environ. Saf. 225:112785. doi: 10.1016/j.ecoenv.2021.112785
Guan, H., Shuai, Y., Ran, Q., Yan, Y., Wang, X., Li, D., et al. (2020). The microbiome and metabolome of Napier grass silages prepared with screened lactic acid bacteria during ensiling and aerobic exposure. Anim. Feed Sci. Tech. 269:114673. doi: 10.1016/j.anifeedsci.2020.114673
Guan, H., Yan, Y. H., Li, X. L., Li, X. M., Shuai, Y., Feng, G. Y., et al. (2018). Microbial communities and natural fermentation of corn silages prepared with farm bunker-silo in Southwest China. Bioresour. Technol. 265, 282–290. doi: 10.1016/j.biortech.2018.06.018
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
Handley, K. M. (2019). Determining microbial roles in ecosystem function: redefining microbial food webs and transcending kingdom barriers. mSystems 4:4. doi: 10.1128/mSystems.00153-19
He, L., Wang, Y., Guo, X., Chen, X., and Zhang, Q. (2022). Evaluating the effectiveness of screened lactic acid bacteria in improving crop residues silage: fermentation parameter, nitrogen fraction, and bacterial community. Front. Microbiol. 13:680988. doi: 10.3389/fmicb.2022.680988
Johnson, L., Harrison, J. H., Hunt, C., Shinners, K., Doggett, C. G., and Sapienza, D. (1999). Nutritive value of corn silage as affected by maturity and mechanical processing: a contemporary review. J. Dairy Sci. 82, 2813–2825. doi: 10.3168/jds.S0022-0302(99)75540-2
Ke, W. C., Ding, W. R., Xu, D. M., Ding, L. M., Zhang, P., Li, F. D., et al. (2017). Effects of addition of malic or citric acids on fermentation quality and chemical characteristics of alfalfa silage. J. Dairy Sci. 100, 8958–8966. doi: 10.3168/jds.2017-12875
Keshri, J., Chen, Y. R., Pinto, R., Kroupitski, Y., Weinberg, Z. G., and Sela, S. (2018). Microbiome dynamics during ensiling of corn with and without Lactobacillus plantarum inoculant. Appl. Microbiol. Biot. 102, 4025–4037. doi: 10.1007/s00253-018-8903-y
Kung, L. M., 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
Li, Y., Du, S., Sun, L., Cheng, Q., Hao, J., Lu, Q., et al. (2022). Effects of lactic acid bacteria and molasses additives on dynamic fermentation quality and microbial community of native grass silage. Front. Microbiol. 13:830121. doi: 10.3389/fmicb.2022.830121
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
Li, Y., and Nishino, N. (2013). Changes in the bacterial community and composition of fermentation products during ensiling of wilted Italian ryegrass and wilted guinea grass silages. Anim. Sci. J. 84, 607–612. doi: 10.1111/asj.12046
Li, P., Zhang, Y., Gou, W. L., Cheng, Q. M., Bai, S. Q., and Cai, Y. M. (2019). Silage fermentation and bacterial community of bur clover, annual ryegrass and their mixtures prepared with microbial inoculant and chemical additive. Anim. Feed Sci. Tech. 247, 285–293. doi: 10.1016/j.anifeedsci.2018.11.009
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. Tech. 197, 206–212. doi: 10.1016/j.anifeedsci.2014.06.016
Liu, C., Wu, H., Liu, S., Chai, S., Meng, Q., and Zhou, Z. (2019). Dynamic alterations in yak rumen bacteria community and metabolome characteristics in response to feed type. Front. Microbiol. 10:1116. doi: 10.3389/fmicb.2019.01116
Lv, R. L., Elsabagh, M., Obitsu, T., Sugino, T., and Kurokawa, Y. (2021). Effect of phytol in forage on phytanic acid content in cow's milk. Anim. Biosci. 34, 1616–1622. doi: 10.5713/ab.21.0086
Ma, S., Fang, C., Sun, X., Han, L., He, X., and Huang, G. (2018). Bacterial community succession during pig manure and wheat straw aerobic composting covered with a semi-permeable membrane under slight positive pressure. Bioresour. Technol. 259, 221–227. doi: 10.1016/j.biortech.2018.03.054
Muck, R. E., Muck, R. E., Nadeau, E. M. G., McAllister, T. A., Contreras-Govea, F. E., Santos, M. C., et al. (2018). Silage review: recent advances and future uses of silage additives. J. Dairy Sci. 101, 3980–4000. doi: 10.3168/jds.2017-13839
Na, N., Qili, M., Wu, N., Sun, L., Xu, H., Zhao, Y., et al. (2022). Bacterial community and fermentation quality of ensiling alfalfa with commercial lactic acid bacterial additives. Front. Microbiol. 13:836899. doi: 10.3389/fmicb.2022.836899
Ni, K. K., Wang, F. F., Zhu, B. G., Yang, J. X., Zhou, G. A., Pan, Y., et al. (2017). Effects of lactic acid bacteria and molasses additives on the microbial community and fermentation quality of soybean silage. Bioresour. Technol. 238, 706–715. doi: 10.1016/j.biortech.2017.04.055
Noecker, C., Chiu, H. C., McNally, C. P., and Borenstein, E. (2019). Defining and evaluating microbial contributions to metabolite variation in microbiome-metabolome association studies. mSystems 4:e00579-19. doi: 10.1128/mSystems.00579-19
Ogunade, I. M., Martinez-Tuppia, C., Queiroz, O. C. M., Jiang, Y., Drouin, P., Wu, F., et al. (2018). Silage review: mycotoxins in silage: occurrence, effects, prevention, and mitigation. J. Dairy Sci. 101, 4034–4059. doi: 10.3168/jds.2017-13788
Pan, L., Li, Z., Wang, Y., Zhang, B., Liu, G., and Liu, J. (2020). Network pharmacology and metabolomics study on the intervention of traditional Chinese medicine Huanglian decoction in rats with type 2 diabetes mellitus. J. Ethnopharmacol. 258:112842. doi: 10.1016/j.jep.2020.112842
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. Tech. 160, 160–166. doi: 10.1016/j.anifeedsci.2010.07.010
Patrica, C. (1997). Official Method of Analysis of AOAC International. Washington, DC: Cunniff Publications.
Pessione, A., Lamberti, C., and Pessione, E. (2010). Proteomics as a tool for studying energy metabolism in lactic acid bacteria. Mol. BioSyst. 6, 1419–1430. doi: 10.1039/C001948H
Santos, A. O., Avila, C. L., 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
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Stergiadis, S., Allen, M., Chen, X. J., Wills, D., and Yan, T. (2015). Prediction of nutrient digestibility and energy concentrations in fresh grass using nutrient composition. J. Dairy Sci. 98, 3257–3273. doi: 10.3168/jds.2014-8587
Thomas, T. A. (1977). An automated procedure for the determination of soluble carbohydrates in herbage. J. Sci. Food Agric. 28, 639–642. doi: 10.1002/jsfa.2740280711
Tsigkrimani, M., Panagiotarea, K., Paramithiotis, S., Bosnea, L., Pappa, E., Drosinos, E. H., et al. (2022). Microbial ecology of sheep milk, artisanal feta, and kefalograviera cheeses. Part II: technological, safety, and probiotic attributes of lactic acid bacteria isolates. Foods 11:459. doi: 10.3390/foods11030459
Van Soest, P. V., 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
Walterson, A. M., and Stavrinides, J. (2015). Pantoea: insights into a highly versatile and diverse genus within the Enterobacteriaceae. FEMS Microbiol. Rev. 39, 968–984. doi: 10.1093/femsre/fuv027
Wang, C., He, L. W., Xing, Y. Q., Zhou, W., Yang, F. Y., Chen, X. Y., 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
Wang, S., Li, J., Zhao, J., Dong, Z., and Shao, T. (2022b). Exploring the ensiling characteristics and bacterial community of red clover inoculated with the epiphytic bacteria from temperate gramineous grasses. J. Appl. Microbiol. 132, 177–188. doi: 10.1111/jam.15234
Wang, S., Shao, T., Li, J., Zhao, J., and Dong, Z. (2022a). Fermentation profiles, bacterial community compositions, and their predicted functional characteristics of grass silage in response to epiphytic microbiota on legume forages. Front. Microbiol. 13:830888. doi: 10.3389/fmicb.2022.830888
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
Wang, Y., Wu, J., Lv, M., Shao, Z., Hungwe, M., Wang, J., et al. (2021). Metabolism characteristics of lactic acid bacteria and the expanding applications in food industry. Front. Bioeng. Biotechnol. 9:612285. doi: 10.3389/fbioe.2021.612285
Wang, S., Yuan, X., Dong, Z., Li, J., and Shao, T. (2017). Effects of ensiling corn stover with legume herbags in different proportions on fermentation characteristics, nutritive quality and in vitro digestibility on the Tibetan plateau. Grassl. Sci. 63, 236–244. doi: 10.1111/grs.12173
Wright, D. A., Gordon, F. J., Steen, R. W. J., and Patterson, D. C. (2000). Factors influencing the response in intake of silage and animal performance after wilting of grass before ensiling: a review. Grass Forage Sci. 55, 1–3. doi: 10.1046/j.1365-2494.2000.00198.x
Xiong, Y., Xu, J., Guo, L., Chen, F., Jiang, D., Lin, Y., et al. (2022). Exploring the effects of different bacteria additives on fermentation quality, microbial community and in vitro gas production of forage oat silage. Animals 12:1122. doi: 10.3390/ani12091122
Xu, D. M., Ding, W. R., Ke, W. C., Li, F. H., Zhang, P., and Guo, X. S. (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
Xu, D., Wang, N., Rinne, M., Ke, W., Weinberg, Z. G., Da, M., et al. (2021). The bacterial community and metabolome dynamics and their interactions modulate fermentation process of whole crop corn silage prepared with or without inoculants. Microb. Biotechnol. 14, 561–576. doi: 10.1111/1751-7915.13623
Yamaga, M., Tani, H., and Murota, K. (2022). Metabolic pathways of cinnamic acid derivatives in Brazilian green propolis in rats. Biosci. Biotechnol. Biochem. 86, 1075–1084. doi: 10.1093/bbb/zbac076
Yan, Y. H., Li, X. M., Guan, H., Huang, L. K., Ma, X., Peng, Y., et al. (2019). Microbial community and fermentation characteristic of Italian ryegrass silage prepared with corn stover and lactic acid bacteria. Bioresour. Technol. 279, 166–173. doi: 10.1016/j.biortech.2019.01.107
Yang, X., Bao, Y., Shao, T., Wang, W., Ma, P., Wang, W., et al. (2022). Altitudinal distribution patterns of phyllosphere microbial communities and their contribution to silage fermentation of kobresia pygmaea along the elevation gradient on the Tibetan plateau. Front. Microbiol. 13:874582. doi: 10.3389/fmicb.2022.874582
Yang, Z., Wang, D., Li, Y., Zhou, X., Liu, T., Shi, C., et al. (2022). Untargeted metabolomics analysis of the anti-diabetic effect of red ginseng extract in type 2 diabetes mellitus rats based on UHPLC-MS/MS. Biomed. Pharmacother. 146:112495. doi: 10.1016/j.biopha.2021.112495
Yin, X., Zhao, J., Wang, S., Dong, Z., Li, J., and Shao, T. (2022). Separating the chemical and microbial factors of oat harvested at two growth stages to determine the main factor on silage fermentation. J. Appl. Microbiol. 132, 4266–4276. doi: 10.1111/jam.15566
Zhang, W., Tan, B., Ye, G., Wang, J., Dong, X., Yang, Q., et al. (2019). Identification of potential biomarkers for soybean meal-induced enteritis in juvenile pearl gentian grouper, Epinephelus lanceolatus♂× Epinephelus fuscoguttatus♀. Aquaculture 512:734337. doi: 10.1016/j.aquaculture.2019.734337
Keywords: Italian ryegrass, bacterial community, metabolomic profiles, fermentation, different regions
Citation: Fu Z, Sun L, Wang Z, Liu Y, Hao J, Gao C and Ge G (2023) Effect of different regions on fermentation profiles, microbial communities, and their metabolomic pathways and properties in Italian ryegrass silage. Front. Microbiol. 13:1076499. doi: 10.3389/fmicb.2022.1076499
Edited by:
Jingyang Luo, Hohai University, ChinaReviewed by:
Yang Wu, Tongji University, ChinaBo Wang, Beijing University of Technology, China
Yinglong Su, East China Normal University, China
Copyright © 2023 Fu, Sun, Wang, Liu, Hao, Gao and Ge. 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: Gentu Ge, 
gegentu@163.com