Hostname: page-component-78c5997874-mlc7c Total loading time: 0 Render date: 2024-11-05T02:11:27.796Z Has data issue: false hasContentIssue false
  • English
  • Français

Investigation of fermentation profiles, bacterial community structure and bacterial β-carotene synthesis of alfalfa silage treated with propionic acid or its combination with squalene

Published online by Cambridge University Press:  17 April 2024

Cheng Zong Open the ORCID record for Cheng Zong [Opens in a new window] ,
Wanqi Jiang ,
Tao Shao  and
Qinhua Liu
Cheng Zong
Affiliation:
Institute of Ensiling and Processing of Grass, College of Agro-grassland Science, Nanjing Agricultural University, Nanjing, China
Wanqi Jiang
Affiliation:
Institute of Ensiling and Processing of Grass, College of Agro-grassland Science, Nanjing Agricultural University, Nanjing, China
Tao Shao
Affiliation:
Institute of Ensiling and Processing of Grass, College of Agro-grassland Science, Nanjing Agricultural University, Nanjing, China
Qinhua Liu*
Affiliation:
Institute of Ensiling and Processing of Grass, College of Agro-grassland Science, Nanjing Agricultural University, Nanjing, China
*
Corresponding author: Qinhua Liu; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

This study examined the fermentation dynamics, bacterial community composition and bacterial β-carotene synthesis in alfalfa that was ensiled for 3, 15, 45 and 90 days without additives (CON), or with a chemical agent (propionic acid, PA); or a combination of PA and squalene, SQPA). The results showed that silage treated with PA had a lower (P < 0.01) pH value than the CON silage in the early ensiling phase (3–15 days). Meanwhile, silage treated with PA had the highest contents of lactic acid, acetic acid and PA after 90 days of ensiling (P < 0.01). The β-carotene in alfalfa was lost seriously in the initial ensiling phase (3 days) and epiphytic Pantoea agglomerans with the ability to produce β-carotene became extinct. With the extension of ensiling time, the loss of β-carotene was alleviated in all silages. PA and SQPA not only lowered bacterial diversity and simplified bacterial networks but also facilitated the emergence of new β-carotene-producing bacteria. The metabolic function prediction indicated that β-carotene synthesis tended to decrease initially and subsequently increase during ensiling. Furthermore, the variance of enzymes involved in the bacterial synthesis of β-carotene in silages was influenced by PA, SQPA and ensiling time. In summary, the impact of solely adding PA demonstrated superior effects on the fermentation quality of alfalfa silage compared to the effects observed with SQPA. Throughout the ensiling of alfalfa, the succession of different β-carotene-producing bacteria resulted in fluctuations in the levels of β-carotene.

Keywords

bacterial community compositionβ-carotene-producing bacteriaensiling timefermentation dynamics
Type
Crops and Soils Research Paper
Information
The Journal of Agricultural Science , Volume 162 , Issue 2 , April 2024 , pp. 118 - 129
Check for updates
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press

Introduction

β-Carotene (pro-vitamin A), the most active biological form among the carotenoids, is beneficial for both human and livestock health (Lindqvist et al., Reference Lindqvist, Nadeau, Jensen and Søegaard2014). Animal products (e.g. meat and milk) play a crucial role in diets worldwide, providing many highly bioavailable essential nutrients (van Hooijdonk and Hettinga, Reference van Hooijdonk and Hettinga2015). The β-carotene concentrations in animal products depend on their diets (Lindqvist et al., Reference Lindqvist, Nadeau, Jensen and Søegaard2014). It makes sense that high β-carotene concentration in feed offers significant advantages not only for livestock health but also for animal products as a β-carotene source for human consumption.

Silage is a primary roughage for ruminants, especially dairy cows, accounting for more than 50% of their rations (Li et al., Reference Li, Zhang, Zhang, Zhang, Usman, Ding, Hao and Guo2022). Silage making is an anaerobic fermentation process that relied mainly on lactic acid bacteria (LAB). Enormous works have confirmed the loss of β-carotene during silage making (Lv et al., Reference Lv, Elsabagh, Obitsu, Sugino, Kurokawa and Kawamura2020; Zong et al., Reference Zong, Wu, Wu, Chen, Dong, Zhao, Shao and Liu2021; Zhang et al., Reference Zhang, Guo, Li, Usman, Zhang and Ding2022). Some reports have shown contradictory outcomes, indicating an increase in β-carotene concentration after ensiling relative to fresh forage material (Lv et al., Reference Lv, Elsabagh, Obitsu, Sugino, Kurokawa and Kawamura2017; Zhao et al., Reference Zhao, Wei, Liu, Kim and Kim2021). Unfortunately, the reason for the diverse findings is unclear. Analogously, Sangija et al. (Reference Sangija, Martin and Matemu2022) explained that lactic acid fermentation increased β-carotene in African nightshade by disrupting the plant matrix and cell cluster during fermentation, leading to carotenoid liberation in food. Additionally, Garrido-Fernandez et al. (Reference Garrido-Fernandez, Maldonado-Barragan, Caballero-Guerrero, Hornero-Mendez and Ruiz-Barba2010) demonstrated that Lactiplantibacillus plantarum could produce carotenoids since it carried genes for carotenoid biosynthesis. However, the above explanations have yet to be tested in silage.

Alfalfa (Medicago sativa L.) is recognized as a preferred legume forage for silage making because of its high nutritional quality (Na et al., Reference Na, Qili, Wu, Sun, Xu, Zhao, Wei, Xue and Tao2022). Nevertheless, alfalfa is difficult to ensile well owing to its high buffering capacity (BC) and low water-soluble carbohydrates (WSC). Spontaneous fermentation of alfalfa silage had poor quality due to undesired microbes, such as Clostridia (Li et al., Reference Li, Jiang, Zheng, Tian, Zheng and Xu2020). To improve alfalfa silage quality, propionic acid (PA) is considered a high-efficiency additive used in silage due to its antimicrobial activity (Dong et al., Reference Dong, Meng and Xiong2006), while the antioxidant squalene can be a carbon source for anaerobic or aerobic microbes (Ghimire et al., Reference Ghimire, Thuan, Koirala and Sohng2016), having the potential to alter bacterial community structure and discover functional bacteria. In our previous work, squalene was detected from fermentation products of alfalfa silage. Moreover, while the addition of squalene prevented the degradation of β-carotene caused by mixed organic acids simulated in a silage environment, thereby increasing the fluorescence intensity of β-carotene, it did not, unfortunately, improve the fermentation quality of alfalfa silage (unpublished data). Considering the advantages of PA and squalene, their combination (SQPA) may have beneficial impacts on both the fermentation quality and β-carotene content of alfalfa silage.

Hence, using single-molecule real-time (SMRT) sequencing technology, the purpose of this study was to: (1) investigate the fermentation profiles, bacterial community structure, and β-carotene variation in alfalfa silage treated with PA or SQPA during ensiling of 3, 15, 45, and 90 days; (2) clarify whether there were functional bacteria that produced β-carotene in alfalfa silage.

Materials and methods

Experimental design and silage preparation

The experiment was conducted using completely randomized blocks, with a 3 × 4 factorial arrangement of 3 additives (alfalfa forage was ensiled without additives as control, CON; 0.4% PA; 0.4% PA + 1% SQ, SQPA), and four ensiling time points (3, 15, 45 and 90 days), with six replicates. The PA and SQPA were diluted in distilled water and applied on a fresh matter (FM) basis. The CON was sprayed with the same volume of distilled water.

Alfalfa was cultivated in the experimental field of Nanjing Agricultural University (32°86′N, 120°02′E) and was harvested from three fields during its full bloom period in May 2020. Fresh alfalfa material was set up with six repetitions for the determination of biochemical properties. The SQ and PA (analytically pure) were both purchased from Shanghai McLean Biochemical Technology Co., Ltd. The harvested alfalfa was delivered right away to the laboratory and chopped into 1–2 cm lengths by a forage chopper. About 720 g of alfalfa for each treatment was filled in a plastic silo (Volume: 1 litre; Density: 206 kg dry matter (DM)/m3). According to the aforementioned experimental design, 2.9 g PA, 10.1 g SQPA, and 10.0 g distilled water were uniformly sprayed on fresh alfalfa. To ensure consistent application rates, an additional 7.1 g of distilled water was added to the PA-treated group. The plastic silos were sealed with a cover and tape, and they were covered in dark plastic bags to block out light. All plastic silos were kept in the laboratory cabinet at room temperature (20–38°C). In the corresponding ensiling time, the silo was opened and sampled. Three random samples of each treatment were used to extract microbial DNA for analysing 16S rRNA genes.

Microbial count and chemical analyses

The plate counting method was used for the determination of the microbial population. Briefly, 10 g fresh alfalfa was blended with 90 ml sterilized saline solution (8.5 g/l NaCl), and serially diluted. The LAB (anaerobiosis, 37°C, 3 days), aerobic bacteria (aerobiosis, 37°C, 3 days), and yeasts (aerobiosis, 37°C, 3 days) were cultured on deMan Rogosa and Sharp agar, nutrient agar, and potato dextrose agar medium, respectively. All microbial data were transformed to log10. The 25 g sample (fresh material or silage) was immersed in 75 ml sterile water and suspended at 4°C overnight. After being filtered through two layers of cheesecloth and a Whatman filter paper, the filtrate was used for determining pH, ammonia-N and organic acids (lactic, acetic, propionic and butyric acid). The pH was measured by a glass electrode pH meter (Hanna Instruments, Inc., Woonsocket, RI, Italy). The ammonia-N was quantified by the phenol-hypochlorite reaction (Broderick and Kang, Reference Broderick and Kang1980). The organic acids were analysed using an HPLC system (Agilent 1260, Agilent Technologies, Santa Clara, CA, USA) equipped with a refractive index detector (Carbomix® H-NP5 column, 2.5 mM H2SO4, 0.5 ml/min). The BC of fresh alfalfa was measured by the titration method (Playne and McDonald, Reference Playne and McDonald1966). Another sample (above 50 g) was dried in a forced-draft oven and DM content was determined by method 930.15 of the Association of Official Analytical Chemists (AOAC, 1990). The dried sample was ground to pass a 1-mm screen with laboratory knife mills (FW100, Taisite Instrument Co., Ltd., Tianjin, China), and the powder sample was used for the analysis of carbohydrates components and crude protein. The total nitrogen (TN) content was measured by method 984.13 (AOAC, 1990), and crude protein was calculated via the TN content multiplied by 6.25. The contents of neutral detergent fibre (aNDFom, aNDFom was assayed with a heat-stable amylase and expressed exclusive of residual ash), acid detergent fibre (aADFom, aADFom was expressed exclusive of residual ash), and acid detergent lignin (ADL) were determined by a ANKOM fibre analyser (ANKOM-200, ANKOM Technologies, Macedon, NY, USA) following the method described by Van Soest et al. (Reference Van Soest, Robertson and Lewis1991). The WSC content was quantified by the anthrone method (Thomas, Reference Thomas1977).

β-carotene analysis

After pre-processing described in our previous study (Liu et al., Reference Liu, Shao and Bai2016), β-carotene content was determined by Agilent HPLC 1260 (Agilent Technologies, Inc., Germany; column: Inertsil ODS-4, GL Science B.V., Japan; detector: Diode Array Detector, DAD VL, Agilent Technologies, Inc., Germany, 450 nm for detecting β-carotene; eluent: 2 ml/min, methanol-acetonitrile 9:1 (v/v); temperature: 45°C).

Measure of enzymes associated with bacterial β-carotene synthesis pathway

The silage samples underwent the following pretreatment: ① The silage samples (freeze-dried and ground, 0.2 g, stem-leaf ratio = 1:3) were added to 1 ml phosphate-buffered saline (pH = 7.4) and mixed thoroughly. ② The mixed solution was centrifuged at 4°C and 3000 × g for 20 min. According to the product manuals provided by Jiangsu Meibiao Biotechnology Co., Ltd (Yancheng, China), enzymes activities, including geranylgeranyl pyrophosphate synthase (EC:2.5.1.29), phytoene synthase (EC:2.5.1.32), phytoene desaturase (EC:1.3.5.5), and lycopene β-cyclase (EC:5.5.1.19), were determined by the corresponding ELISA Kits, respectively.

Microbial sequencing and bioinformatic analysis

The PacBio sequencing process was divided into six steps. Briefly, the TIANamp Bacteria DNA isolation kit (DP302-02, Tiangen, Beijing, China) was used to extract microbial community genomic DNA. Then, the bacterial 16S rRNA genes were amplified using the universal bacterial primers (27F and 1492R). The Polymerase Chain Reactions (PCR) were performed in three samples of each treatment. After electrophoresis, PCR products were purified using the AMPure® PB beads (Pacifc Biosciences, CA, USA) and quantified with Quantus™ Fluorometer (Promega, Madison, WI, USA). Purified products were pooled in equimolar amounts and DNA library was constructed using the SMRTbell® Express Template Prep Kit 2.0 (Pacifc Biosciences, Menlo Park, CA, USA). Purified SMRTbell libraries were sequenced on the Pacbio Sequel II System (Pacifc Biosciences, CA, USA) by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Circular consensus sequence (CCS) was extracted and filtered with the SMRTLink 8.0 (minimum full pass = 3, minimum predicted accuracy = 0.99), and sequences outside of a range of < 1400 bp and > 1600 bp were removed. Optimized-CCS were used for further downstream analyses. Operational taxonomic units (OTUs) with 99% similarity cutoff were clustered via USEARCH (version 7.0), and chimeric sequences were identified and removed (Stackebrandt and Goebel, Reference Stackebrandt and Goebel1994; Edgar, Reference Edgar2013). Each OTU representative sequence was analysed to obtain classification information by RDP Classifier (version 2.13) against the 16S rRNA database (Silva v138) with a confidence threshold of 75% (Wang et al., Reference Wang, Garrity, Tiedje and Cole2007). The raw sequencing data have been deposited in the NCBI database (accession: PRJNA941960).

The α-diversity analysis (Shannon and Chao 1) was carried out by MOTHUR (version 1.30.2). The principal coordinate analysis (PCoA) was graphically constructed by R software (version 3.3.1). Bacterial Community structures were analysed from genus to species levels by R software. To illuminate the correlations between fermentation indices and bacteria (on species level), spearman's correlation heatmaps of different silages were established by R software (pheatmap package). The co-occurrence networks were statistically calculated by Networkx software (version 2.8) based on Spearman's correlation and were visualized by Gephi software (version 0.9.2). Microbial functions based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were measured by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2; https://github.com/picrust/picrust2; Version 2.2.0), which predicted the function of bacterial community based on the proportion of marker gene sequences in samples (Douglas et al., Reference Douglas, Maffei, Zaneveld, Yurgel, Brown, Taylor, Huttenhower and Langille2020).

Statistical analyses

The statistical analysis was performed with IBM SPSS 22.0 for Windows statistical software. Initially, all data relating to fermentation quality, chemical composition, and loss rate of β-carotene of alfalfa silage was subjected to Shapiro–Wilk test for normality and Levene test for homogeneity of variance, respectively, to establish the basic assumptions for conducting analysis of variance (ANOVA, GLM). The effects of additives, ensiling time, and their interaction on the fermentation quality, chemical composition, and loss rate of β-carotene of alfalfa silage were analysed using a two-way ANOVA, according to the statistical model below: Y ij = μ + α i + β j + (αβ × ij) + e ij, where Y ij = dependent variable; μ = overall mean; α i = fixed effect of additives; β j = fixed effect of ensiling time; (αβ × ij) = interaction between additives and ensiling time; and e ij = residual error term. In cases where interaction effects were significant, multiple comparisons for means were conducted using the Tukey's test for data with a single fixed effect.

The indices of Shannon and Chao 1 was analysed using a one-way ANOVA for the following model: Y ij = μ + α i + e ij, where Y ij = dependent variable; μ = overall mean; α i = fixed effect of treatments; and e ij = residual error term. The difference between groups was compared with Tukey's test for determining the significant (P < 0.05).

Beta diversity was assessed using the dispersion level of sample communities based on Bray-Curtis distance, while the dissimilarity among treatments was assessed using the permutational multivariate analysis of variance (PERMANOVA).

The difference of two modules related to β-carotene synthesis (M000364 and M0097) predicted by PICRUSt2 was analysed using a one-way ANOVA for the following model: Y ij = μ + α i + e ij, where Y ij = dependent variable; μ = overall mean; α i = fixed effect of additives; and e ij = residual error term. The multiple comparisons of the means were conducted by Tukey's test at the significance level (P < 0.05).

Results

Biochemical properties of fresh alfalfa

The chemical and microbiological characteristics of fresh alfalfa are shown in Table 1.

Table 1. Chemical and biological composition of fresh alfalfa

a DM, dry matter; FM, fresh matter; mEq, milligram equivalent; cfu, colony-forming units.

Chemical composition, fermentation dynamics and microbial counts of alfalfa silage

The interaction of additives and ensiling time had significant (P < 0.01) effects on DM, crude protein, pH, organic acids (lactic, acetic, propionic and butyric acid), ammonia-N, LAB and aerobic bacteria, while had non-significant (P > 0.05) effects on fiber fractions (aNDFom, aADFom and ADL), WSC and yeasts (Tables 2–4). Throughout the ensiling process, PA- and SQPA-treated silages had a statistically or numerically higher DM content than CON silage. After 3 days of ensiling, PA-treated silage had lower crude protein than CON silage, while its content declined in CON silage from 15 to 90 days. The PA-treated silage had a lower (P < 0.01) pH value than that in CON silage in the early ensiling phase (3–15 days). Meanwhile, PA-treated silage had the highest (P < 0.01) contents of lactic acid, acetic acid and PA after 90 days of ensiling. Ammonia-N accumulated continuously during the ensiling process. After 90 days of ensiling, the ammonia-N in PA-treated silage was higher (P < 0.05) than that in CON and SQPA-treated silages. The PA-treated silages had the lowest (P < 0.01) butyric acid content after 3 days of ensiling. From the perspective of microbe, with the prolonged ensiling time, The LAB counts were on the rise, while aerobic bacteria were exhibiting a downward trend. The LAB counts in all treated silages peaked after 45 days of ensiling, and PA-treated silage had more (P < 0.01) LAB counts than that in CON and SQPA-treated silages. The PA-treated silage had lower (P < 0.01) aerobic bacteria counts than that in CON silage after 45 and 90 days of ensiling.

Table 2. Chemical composition of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

a FM, fresh matter; DM, dry matter; aNDFom, neutral detergent fibre; aADFom, acid detergent fibre; ADL, acid detergent lignin.

b CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene.

c SEM, standard error of the means.

d A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time.

a-bValues with different small letters show significant differences among additives in the same ensiling days (P < 0.05).

A-CValues with different capital letters show significant differences among ensiling days in the same additives (P < 0.05).

Table 3. Fermentation features of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

a DM, dry matter; TN, total nitrogen; WSC, water-soluble carbohydrates.

b CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene.

c SEM, standard error of the means.

d A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time.

a-cValues with different small letters show significant differences among additives in the same ensiling days (P < 0.05).

A-CValues with different capital letters show significant differences among ensiling days in the same additives (P < 0.05).

Table 4. Microbial counts of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

a LAB, lactic acid bacteria; FM, fresh matter; cfu, colony-forming units.

b CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene.

c SEM, standard error of the means.

d A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time.

a-bValues with different small letters show significant differences among additives in the same ensiling days (P < 0.05).

A-DValues with different capital letters show significant differences among ensiling days in the same additives (P < 0.05).

Dynamic changes of β-carotene and lipoxygenase in alfalfa silage

The β-carotene loss occurred in alfalfa silage after ensiling, and its dynamic change is presented in Fig. 1. The interaction of additives and ensiling time had significant (P < 0.01) effects on β-carotene loss. The lower (P < 0.01) β-carotene loss occurred in PA-treated silage for 3 to 45 days of ensiling and in SQPA-treated silage for 3 and 45 days than that in CON silage. The higher (P < 0.01) β-carotene loss occurred in PA- and SQPA-treated silages for 90 days than that in CON silage.

Figure 1. The dynamic of loss rate of β-carotene in alfalfa silage treated with propionic acid or a combination of propionic acid and squalene. Loss rate of β-carotene (%) = (initial content − final content)/initial content × 100%; initial content refers to the β-carotene in fresh alfalfa; final content refers to the β-carotene after the corresponding ensiling time. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time; Values with different small letters show significant differences among additives in the same ensiling days (P < 0.05); Values with different capital letters show significant differences among ensiling days in the same additives (P < 0.05).

Diversity of the bacterial community in alfalfa silage

The bacterial community structure of fresh alfalfa and silages is presented in Fig. 2. The α-diversity analysis (Fig. 2a) showed that Shannon and Chao 1 indices were significantly (P < 0.01) higher in fresh alfalfa than that in silages. The abundance heatmap displayed that the top 15 of bacteria belonged to Weissella, Lactobacillus, Enterococcus, Pantoea, Lactococcus and Enterobacter in fresh alfalfa and silages. On 3 days, predominant LAB were Lactobacillus, Weissella, Enterococcus and Lactococcus in CON silage, whereas predominant LAB were only Weissella in PA- and SQPA-treated silages. With the extension of ensiling time, Weissella was succeeded by Lactobacillus in PA- and SQPA-treated silages. Concretely, on species level, as for PA- and SQPA-treated silages, specie Weissella cibaria dominated after 3 days of ensiling. Species L. plantarum, Lacticaseibacillus paracasei, Apilactobacillus kunkeei, Lentilactobacillus parabuchneri and Fructilactobacillus kullabergensis dominated from 15 to 90 days of ensiling.

Figure 2. The α-diversity (a) and abundance heatmap (top 15, b) of fresh alfalfa and after ensiling. FA, fresh alfalfa; CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45, and 90 days; *, P < 0.05; **, P < 0.01; ***, P < 0.001.

The β-diversity analysis showed the similarities and differences between bacterial communities (Fig. 3). In every treatment cohort throughout the ensiling process, silage was completely (P < 0.01) separated from fresh alfalfa (Figs 3a–d). The PA- and SQPA-treated silages were obviously (P < 0.01) separated from CON silage after ensiling for 3 and 15 days (Figs 3a and b). Subsequently, all silages were gathered together (Figs 3c and d).

Figure 3. The β-diversity based on principal coordinates analysis of bacterial communities in fresh alfalfa and alfalfa ensiled with different treatments for 3 days (a), 15 days (b), 45 days (c), and 90 days (d), respectively. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45, and 90 days.

Correlation between bacterial flora and fermentation products as well as their internal network in alfalfa silage

The correlation analysis between bacterial community and fermentation products in different silages is shown in Fig. 4. In CON silage, notable (P < 0.05, 1 > |r| > 0.6) correlations were observed between LAB species (A. kunkeei and F. kullabergensis) and pH, WSC, butyric acid and acetic acid. Significant (P < 0.05, 1 > |r| > 0.6) correlations were found between L. plantarum and butyric acid and lactic acid. In PA-treated silage, Levilactobacillus brevis and L. plantarum showed a significant (P < 0.01, 1 > |r| > 0.6) correlation with pH. Significant (P < 0.05, 1 > |r| > 0.6) correlations were found between LAB species (A. kunkeei and F. kullabergensis) and WSC and acetic acid. In SQPA-treated silage, Lactobacillus helsingborgensis, Lactobacillus apis, F. kullabergensis, A. kunkeei and Bombilactobacillus mellis had significant (P < 0.05, 1 > |r| > 0.6) correlations with pH, lactic acid, WSC and acetic acid. Besides, Bifidobacterium asteroids had significant (P < 0.05, 1 > |r| > 0.6) correlations with WSC and acetic acid. Lacticaseibacillus paracasei and L. parabuchneri had significant (P < 0.05, 1 > |r| > 0.6) correlations with lactic acid.

Figure 4. Correlation analysis between bacterial communities (top 15 species) and fermentation characteristics in different alfalfa silages. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; AA, acetic acid; LA, lactic acid; PA, propionic acid; BA, butyric acid; WSC, water-soluble carbohydrates; *, P < 0.05; **, P < 0.01; ***, P < 0.001.

As exhibited in Fig. 5, distinct co-occurrence patterns were observed in the bacterial communities of various treated silages. The node and edge numbers of the network reflected the network complexity. The node and edge numbers underwent irregular changes with ensiling time. Overall, the edges of PA- and SQPA-treated silages were numerically less than those of CON silage. Additionally, the ratio of negative and positive (N/P) represented the network stability. The value of N/P in PA-treated silage was numerically lower than that in CON silage.

Figure 5. The bacterial co-occurrence networks (on the species level) in different alfalfa silages. Nodes represent individual species; edges represent significant Spearman correlations (|ρ| > 0.5, P < 0.05). The size of each node is proportional to the number of connected edges, and the nodes are labelled by species and coloured by genus. The colour of the edges corresponds to a positive (red) or negative (blue) relationship; CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene.

Bacterial metabolic function prediction and analysis of enzymes associated with bacterial β-carotene synthesis pathway in alfalfa silage

Utilizing the KEGG pathways database with PICRUSt2, predictive bacterial function profiles are shown in Fig. 6a. According to the findings, the metabolism of cofactors and vitamins was among the top 10 of all metabolic pathways in alfalfa silage. Further results showed that (Fig. 6b), regarding bacteria, two modules related to β-carotene synthesis, C10-C20 isoprenoid biosynthesis (M00364) and β-carotene biosynthesis (M00097), were altered during ensiling. Ensiling time had significant (P < 0.01) effects on both change of M00364 and M00097. The modules of M00364 and M00097 drastically reduced in the initial stage of ensiling, followed by a descent remission as the ensiling time extended. After ensiling for 90 days, M00364 and M00097 showed a positive increase.

Figure 6. The bacterial metabolic function prediction in fresh alfalfa and silages were analysed by PICRUSt2 (a), the change of C10-C20 isoprenoid biosynthesis and β-carotene biosynthesis before and after ensiling (b), and principal coordinates analysis of enzymes (geranylgeranyl pyrophosphate synthase, phytoene synthase, phytoene desaturase, lycopene β-cyclase) associated with bacterial β-carotene synthesis pathway in different alfalfa silages (c). Change (C10-C20 isoprenoid biosynthesis or β-carotene biosynthesis) = relative abundance (after ensiling) – relative abundance (before ensiling). CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45 and 90 days; A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time.

Using PCoA (Fig. 6c), the variances of enzymes (EC:2.5.1.29, EC:2.5.1.32, EC:1.3.5.5, EC:5.5.1.19) involved in the bacterial synthesis of β-carotene in different silages were characterized. The results exhibited significant differences in enzyme activities in different treated silages (PERMANOVA, R 2 = 0.54, P < 0.01). Additionally, ensiling time resulted in varying degrees of alterations (P < 0.01) in enzyme activities.

Discussion

Fermentability of pre-ensiled alfalfa

In this study, fresh alfalfa possessed appropriate BC, which was far below the value from some existing reports (Turan and Onenc, Reference Turan and Onenc2018; Yang et al., Reference Yang, Zhao, Wang, Fan, Wang and Feng2020), but was close to the result stated by Wang et al. (Reference Wang, Li, Zhao, Dong, Dong and Shao2021a) (317 mEq/kg DM). Low BC reduces the resistance to pH decrease, which was conducive to fermentation. Moreover, the WSC (>50 g/kg DM) content met the requirement for acceptable silage (Li et al., Reference Li, Zhang, Gou, Cheng, Bai and Cai2019). In addition to chemical characteristics, the microbial composition is crucial for fermentability. Epiphytic microflora is responsible for silage fermentation, of which LAB are the most concerned becoming a significant factor in predicting the adequacy of fermentation (Lin et al., Reference Lin, Bolsen, Brent, Hart, Dickerson, Feyerherm and Aimutis1992). According to Cai (Reference Cai1999), the minimum amount of epiphytic LAB to obtain well-preserved silage was 5 log cfu/g FM. In addition to LAB, epiphytic microflora includes many undesirable microbes, which can bring serious fermentation loss (Wróbel et al., Reference Wróbel, Nowak, Fabiszewska, Paszkiewicz-Jasińska and Przystupa2023). Therefore, insufficient epiphytic LAB and a high number of aerobic bacteria were unfavourable for spontaneous fermentation (Cai et al., Reference Cai, Benno, Ogawa, Ohmomo, Kumai and Nakase1998).

Chemical composition and fermentation profiles of alfalfa silage

With regards to non-carbohydrate components, the WSC was fermented well by LAB during ensiling, leading to no significant difference in residual WSC among all treated silages (Huo et al., Reference Huo, Zhang, Zhang, Shen, Chen, Liu, Zhang, Wang and Guo2022). As for structural carbohydrates, the contents of aNDFom, aADFom and ADL were not significantly affected by PA and SQPA, which was in accordance with the study conducted by Zhang et al. (Reference Zhang, Wang, Li, Lin, Yang and Ni2019). Fibre degradation depends on enzymatic or strong acid hydrolysis (Huisden et al., Reference Huisden, Adesogan, Kim and Ososanya2009). It could be speculated that fibre degradation was limited by insufficient enzyme activity and acidity in this study. The high production of lactic acid and rapid decline in pH value in all treated silages after 3 days of ensiling indicated strong lactic acid fermentation (Wang et al., Reference Wang, Li, Zhao, Dong, Dong and Shao2021a). The addition of PA directly formed an acidic environment, which was conducive to the growth of LAB in the early stage (Kim et al., Reference Kim, Amanullah, Lee, Joo and Kim2015). This was consistent with the outcome of this study, statistically or numerically more LAB counts in PA- and SQPA-treated silages than that in CON silage after 3 and 15 days of ensiling. Additionally, the high accumulation of acetic acid during ensiling was related to low DM content (244–272 g/kg DM). Excessive acetic acid (>4–6%) is often detected in extremely wet (>70% moisture) silages (Kung et al., Reference Kung, Shaver, Grant and Schmidt2018). Interestingly, compared with the CON, the accumulation rate of acetic acid was slower in PA-treated silage during the early to mid-stage of fermentation (3–45 days), but faster during the later stage. Adding organic acid as an additive created conditions for subsequent growth of heterofermentative LAB in the later stage of fermentation (Tyrolová et al., Reference Tyrolová, Bartoň and Loučka2017). Due to multiple uncontrollable factors, the change of PA in different treated silages was hard to explain with the current evidence available. As a volatile substance, PA loss is inevitable during the fermentation process. Meanwhile, extra PA was generated during the fermentation process. Theoretically, PA should be higher in silages where it was added. Nevertheless, these uncertainties can impact the final outcome. Kim et al. (Reference Kim, Amanullah, Lee, Joo and Kim2015) found that PA was higher in silages treated with PA at early stages of fermentation (2–7 days), but did not persist to the end (100 days). Excessive proteolysis frequently occurs in poor-quality alfalfa silage. It is caused by the multiplication of clostridia, leading to the formation of ammonia-N and butyric acid (Kung et al., Reference Kung, Shaver, Grant and Schmidt2018). In this study, the concentrations of ammonia-N and butyric acid detected in all silages were less than the threshold of 150 g/kg TN and 5 g/kg DM, respectively, indicating well-conserved legume silage (Kung et al., Reference Kung, Shaver, Grant and Schmidt2018). After 90 days of ensiling, PA- and SQPA-treated silages reduced the formation of ammonia-N. This was attributed to the rapid reduction of pH value due to the great production of lactic acid, suppressing the decarboxylation or deamination of amino acids by undesirable clostridia (Oliveira et al., Reference Oliveira, Weinberg, Ogunade, Cervantes, Arriola, Jiang, Kim, Li, Goncalves, Vyas and Adesogan2017). Adding PA increased the lactic acid content, indicating that silage quality was improved, which is consistent with the findings of Jia et al. (Reference Jia, Yun and Yu2021). On the other hand, the combination of squalene and PA weakened the effect of PA in improving the fermentation quality. The reason was that SQPA promoted the growth of undesirable Enterobacter sp. (Fig. 2b).

The PA, SQPA, and ensiling time influenced β-carotene loss in alfalfa silage

In this study, with the extension of fermentation, the β-carotene loss was alleviated, especially PA-treated silage for 45 days and CON silage for 90 days. There are two explanations for this phenomenon. The first explanation is that the extraction of β-carotene derives from stable lipoprotein complexes due to fermentation. Another explanation is that β-carotene is synthesized by LAB. Fermentation disrupts the plant matrix and cell cluster, contributing to carotenoid liberation and bioaccessibility (Sangija et al., Reference Sangija, Martin and Matemu2022). Furthermore, some LAB are known to synthesize carotenoids, including L. plantarum, Enterococcus gilvus and Lactobacillus gasseri (Garrido-Fernandez et al., Reference Garrido-Fernandez, Maldonado-Barragan, Caballero-Guerrero, Hornero-Mendez and Ruiz-Barba2010; Hagi et al., Reference Hagi, Kobayashi and Nomura2015; Xu et al., Reference Xu, Hlaing, Glagovskaia, Augustin and Terefe2020). Therefore, silage could be a candidate resource bank for β-carotene-synthesizing LAB.

The PA and SQPA simplified the microflora and altered the keystone taxon

Anaerobic fermentation resulted in the stability of microflora in alfalfa silage. This was because dominant species constrained species diversity (Polley et al., Reference Polley, Wilsey and Derner2007). In the initial ensiling phase (3 days), PA and SQPA simplified the microbial community composition by promoting the growth of Weissella (W. cibaria), supported by lower diversity indices (Chao 1 and Shannon, Fig. 2a). This deduced that resistance to low pH of Weissella was stronger than other cocci LAB, such as Enterococcus and Lactococcus. Weissella gender are thought to be early colonizers and then they are outcompeted by acid-tolerant Lactobacillus gender due to the pH drop as fermentation progresses (Ogunade et al., Reference Ogunade, Jiang, Pech Cervantes, Kim, Oliveira, Vyas, Weinberg, Jeong and Adesogan2018). Enterococcus and Weissella species constitute a significant portion of the microbial flora typically found on the surfaces of various forage plants (Cai, Reference Cai1999). Due to the lower acid tolerance, the high abundance of Enterococcus in CON in initial ensiling phase indicated that spontaneous fermentation did not achieve a rapid decrease in pH. Enterobacter species are usually undesirable bacteria, which can compete with LAB for the available sugars and can also break down proteins (Wróbel et al., Reference Wróbel, Nowak, Fabiszewska, Paszkiewicz-Jasińska and Przystupa2023). Hence, with the extension of fermentation, the succession from W. cibaria to L. plantarum, L. brevis, L. paracasei, A. kunkeei, L. parabuchneri and F. kullabergensis occurred in PA- and SQPA-treated silages (Fig. 2b). The β-diversity analysis revealed that PA and SQPA influenced the differences of bacterial flora in the early ensiling phase (3 and 15 days). The PA could directly acidify the anaerobic environment and inhibit the growth of undesirable microbes in the early stage of ensiling and accelerated the succession of the bacterial community. As mentioned above, W. cibaria swiftly took over as the dominating bacterial species after 3 days of ensiling. It could not be ignored that microflora separation occurred in the early ensiling stages (3 and 15 days) while aggregation happened in the mid-to-late ensiling stages (45 and 90 days) among all silages. This implied that distinct treatment groups possessed varying microflora in the early ensiling stages, while convergence in microflora is observed in the mid-to-late ensiling stages. Silage additives started acting early in the fermentation process, yet well-fermented silage harboured homogeneous microflora in the stable stage of late fermentation. Likely, Eliasson et al. (Reference Eliasson, Sun, Lundh, Hojer, Saeden, Hetta and Gonda2023) stated that inoculation resulted in well fermented silages comprising a homogenous microbiota.

The performance of bacteria impacted the fermentation quality. In this study, Lactobacillus was the key bacteria for encouraging fermentation. Notably, in PA- and SQPA-treated silages, heatmaps showed that different species belonging to the Lactobacillus played a pivotal role in improving fermentation quality. Xu et al. (Reference Xu, Wang, Rinne, Ke, Weinberg, Da, Bai, Zhang, Li and Guo2021) pointed out that keystone taxon in different silages were entirely distinct. The construction of bacterial networks is commonly employed to illustrate the micro-ecology in silage, wherein the fermentation quality of silage can be highlighted by the complexity of the network. Bai et al. (Reference Bai, Ding, Ke, Xu, Wang, Huang, Zhang, Liu and Guo2021) stated that bacterial correlation structures in silage with a high fermentation quality were more straightforward with lower network complexity. Thus, from the perspective of the microbes, adding PA and SQPA simplified the bacterial network, improving fermentation quality. Likewise, compared with SQPA, PA further simplified the bacterial network. Namely, PA-treated silage had better fermentation quality. Besides, Bai et al. (Reference Bai, Ding, Su, Wang, Cheng, Xie and Guo2022) indicated that fermentation quality was negatively correlated with bacterial community stability. This was in line with this study, substantiated by a lower N/P in PA-treated silage.

Prolonging ensiling time promoted bacterial β-carotene synthesis

The functional prediction could be used to evaluate the effect of bacterial communities on metabolic pathways that underlie silage production (Bai et al., Reference Bai, Ding, Ke, Xu, Wang, Huang, Zhang, Liu and Guo2021). The metabolism of cofactors and vitamins was one of the top 10 microbial metabolic processes of silage. The modules (M00364 and M00096) involving β-carotene synthesis declined in the initial ensiling phase (3 days), indicating the limitation of bacterial β-carotene synthesis. This was in line with the fact that the anaerobic fermentation process weakened the β-carotene related microbes on the surface of the grass. Among these bacteria, Pantoea agglomerans carried the gene crtEBIY for β-carotene synthesis and was a β-carotene producer (Wang et al., Reference Wang, Liu, Jiang and Mao2021b). With the extension of ensiling time, bacterial community succession occurred and novel β-carotene producers might emerge. This was confirmed by the recovery of bacterial β-carotene synthesis. This could explain the phenomenon reported by Zhao et al. (Reference Zhao, Wei, Liu, Kim and Kim2021), who observed an increase in β-carotene in rye silage after 60 days of ensiling. Although the PA and SQPA did not have a significant effect on M00364 and M00096, based on PCoA, it was demonstrated that PA and SQPA could affect the bacterial synthesis of enzymes involved in the β-carotene pathway. Similarly, Lee et al. (Reference Lee, Lee, Cho, Kim and Kim2018) found that red pepper fermented with L. parabuchneri can change carotenoids markedly during the fermentation process. Moreover, many parameters (carbon source, temperature and etc) can affect the biotechnological synthesis of β-carotene (Gupta et al., Reference Gupta, Adin, Panda and Mujeeb2022). The use of additives altered the conditions for bacterial synthesis of β-carotene. These findings implied that further investigation into the microbial viewpoint on β-carotene regulation in silage was promising.

Conclusion

Adding PA alone had a better effect than SQPA on improving fermentation quality. The bacterial β-carotene synthesis occurred during alfalfa ensiling and prolonging the ensiling time benefited its synthesis. The PA and SQPA mainly exerted roles in bacterial β-carotene synthesis in the initial and middle stage. Future studies will concentrate on screening and purifying β-carotene-producing bacteria, thereby contributing to animal nutrition and productivity.

Author contributions

Qinhua Liu and Tao Shao: Conceptualization, Methodology, Writing – review & editing, Funding acquisition, Project administration. Cheng Zong: Investigation, Data curation, Visualization, Writing-original draft. Wanqi Jiang: Investigation, Data curation, Visualization.

Funding statement

This work was supported by the National Natural Science Foundation of China (31971765), the Key research and development projects of Hainan province (ZDYF2022XDNY153), and the Postgraduate Research and Practice Innovation Programme of Jiangsu Province (KYCX22_0782).

Competing interests

The authors declare that they have no financial interests or personal relationships that could influence the work reported in this paper.

Ethical standards

This article does not contain any studies with human participants or animals performed by any of the authors.

References

AOAC (1990) Official Methods of Analysis, 15th Edn. Arlington, VA, USA: Association of Official Analytical Chemists.Google Scholar
Bai, J, Ding, Z, Ke, W, Xu, D, Wang, M, Huang, W, Zhang, Y, Liu, F and Guo, X (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. Microbial Biotechnology 14, 1171–1182.10.1111/1751-7915.13785CrossRefGoogle ScholarPubMed
Bai, J, Ding, Z, Su, R, Wang, M, Cheng, M, Xie, D and Guo, X (2022) Storage temperature is more effective than lactic acid bacteria inoculations in manipulating fermentation and bacterial community diversity, co-occurrence and functionality of the whole-plant corn silage. Microbiology Spectrum 10, e0010122.10.1128/spectrum.00101-22CrossRefGoogle ScholarPubMed
Broderick, GA and Kang, JH (1980) Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media. Journal of Dairy Science 63, 64–75.10.3168/jds.S0022-0302(80)82888-8CrossRefGoogle ScholarPubMed
Cai, Y (1999) Identification and characterization of Enterococcus species isolated from forage crops and their influence on silage fermentation. Journal of Dairy Science 82, 2466–2471.10.3168/jds.S0022-0302(99)75498-6CrossRefGoogle ScholarPubMed
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. Applied and Environmental Microbiology 64, 2982–2987.10.1128/AEM.64.8.2982-2987.1998CrossRefGoogle ScholarPubMed
Dong, XL, Meng, QX and Xiong, YQ (2006) A comparison of buffered propionic acid and Propionibacterium acidipropionici as additives for high oil maize stover silage. Journal of Animal and Feed Sciences 15, 669–683.10.22358/jafs/66939/2006CrossRefGoogle Scholar
Douglas, GM, Maffei, VJ, Zaneveld, J, Yurgel, SN, Brown, JR, Taylor, CM, Huttenhower, C and Langille, MGI (2020) PICRUSt2: an improved and customizable approach for metagenome inference. BioRxiv 672295, 1–16.Google Scholar
Edgar, RC (2013) UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature Methods 10, 996–998.10.1038/nmeth.2604CrossRefGoogle ScholarPubMed
Eliasson, T, Sun, L, Lundh, A, Hojer, A, Saeden, KH, Hetta, M and Gonda, H (2023) Epiphytic microbiota in Swedish grass-clover herbage and the effect of silage additives on fermentation profiles and bacterial community compositions of the resulting silages. Journal of Applied Microbiology 134, 1–17.10.1093/jambio/lxad196CrossRefGoogle ScholarPubMed
Garrido-Fernandez, J, Maldonado-Barragan, A, Caballero-Guerrero, B, Hornero-Mendez, D and Ruiz-Barba, JL (2010) Carotenoid production in Lactobacillus plantarum. International Journal of Food Microbiology 140, 34–39.10.1016/j.ijfoodmicro.2010.02.015CrossRefGoogle ScholarPubMed
Ghimire, GP, Thuan, NH, Koirala, N and Sohng, JK (2016) Advances in biochemistry and microbial production of squalene and its derivatives. Journal of Microbiology and Biotechnology 26, 441–451.10.4014/jmb.1510.10039CrossRefGoogle Scholar
Gupta, I, Adin, SN, Panda, BP and Mujeeb, M (2022) beta-Carotene-production methods, biosynthesis from Phaffia rhodozyma, factors affecting its production during fermentation, pharmacological properties: a review. Biotechnology and Applied Biochemistry 69, 2517–2529.10.1002/bab.2301CrossRefGoogle ScholarPubMed
Hagi, T, Kobayashi, M and Nomura, M (2015) Aerobic conditions increase isoprenoid biosynthesis pathway gene expression levels for carotenoid production in Enterococcus gilvus. FEMS Microbiology Letters 362, fnv075.10.1093/femsle/fnv075CrossRefGoogle ScholarPubMed
Huisden, CM, Adesogan, AT, Kim, SC and Ososanya, T (2009) Effect of applying molasses or inoculants containing homofermentative or heterofermentative bacteria at two rates on the fermentation and aerobic stability of corn silage. Journal of Dairy Science 92, 690–697.10.3168/jds.2008-1546CrossRefGoogle ScholarPubMed
Huo, W, Zhang, Y, Zhang, L, Shen, C, Chen, L, Liu, Q, Zhang, S, Wang, C and Guo, G (2022) Effect of lactobacilli inoculation on protein and carbohydrate fractions, ensiling characteristics and bacterial community of alfalfa silage. Frontiers in Microbiology 13, 1070175.10.3389/fmicb.2022.1070175CrossRefGoogle ScholarPubMed
Jia, T, Yun, Y and Yu, Z (2021) Propionic acid and sodium benzoate affected biogenic amine formation, microbial community, and quality of oat silage. Frontiers in Microbiology 12, 750920.10.3389/fmicb.2021.750920CrossRefGoogle ScholarPubMed
Kim, DH, Amanullah, SM, Lee, HJ, Joo, YH and Kim, SC (2015) Effect of microbial and chemical combo additives on nutritive value and fermentation characteristic of whole crop barley silage. Asian-Australasian Journal of Animal Science 28, 1274–1280.10.5713/ajas.15.0106CrossRefGoogle ScholarPubMed
Kung, L, Shaver, RD, Grant, RJ and Schmidt, RJ (2018) Silage review: interpretation of chemical, microbial, and organoleptic components of silages. Journal of Dairy Science 101, 4020–4033.10.3168/jds.2017-13909CrossRefGoogle ScholarPubMed
Lee, SM, Lee, JY, Cho, YJ, Kim, MS and Kim, YS (2018) Determination of volatiles and carotenoid degradation compounds in red pepper fermented by Lactobacillus parabuchneri. Journal of Food Science 83, 2083–2091.10.1111/1750-3841.14221CrossRefGoogle ScholarPubMed
Li, P, Zhang, Y, Gou, W, Cheng, Q, Bai, S and Cai, Y (2019) Silage fermentation and bacterial community of bur clover, annual ryegrass and their mixtures prepared with microbial inoculant and chemical additive. Animal Feed Science and Technology 247, 285–293.10.1016/j.anifeedsci.2018.11.009CrossRefGoogle 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. Science Reports 10, 17782.10.1038/s41598-020-74958-1CrossRefGoogle ScholarPubMed
Li, F, Zhang, B, Zhang, Y, Zhang, X, Usman, S, Ding, Z, Hao, L and Guo, X (2022) Probiotic effect of feruloyl esterase-producing Lactobacillus plantarum inoculated alfalfa silage on digestion, antioxidant, and immunity status of lactating dairy goats. Animal Nutrition 11, 38–47.10.1016/j.aninu.2022.06.010CrossRefGoogle Scholar
Lin, C, Bolsen, KK, Brent, BE, Hart, RA, Dickerson, JT, Feyerherm, AM and Aimutis, WR (1992) Epiphytic microflora on alfalfa and whole-plant corn. Journal of Dairy Science 75, 2484–2493.10.3168/jds.S0022-0302(92)78010-2CrossRefGoogle ScholarPubMed
Lindqvist, H, Nadeau, E, Jensen, SK and Søegaard, K (2014) α-Tocopherol and β-carotene contents of forage species in a four-cut system. Grass and Forage Science 69, 356–364.10.1111/gfs.12058CrossRefGoogle Scholar
Liu, QH, Shao, T and Bai, YF (2016) The effect of fibrolytic enzyme, Lactobacillus plantarum and two food antioxidants on the fermentation quality, alpha-tocopherol and beta-carotene of high moisture napier grass silage ensiled at different temperatures. Animal Feed Science and Technology 221, 1–11.10.1016/j.anifeedsci.2016.08.020CrossRefGoogle Scholar
Lv, R, Elsabagh, M, Obitsu, T, Sugino, T, Kurokawa, Y and Kawamura, K (2017) Effects of nitrogen fertilizer and harvesting stage on photosynthetic pigments and phytol contents of Italian ryegrass silage. Animal Science Journal 88, 1513–1522.10.1111/asj.12810CrossRefGoogle ScholarPubMed
Lv, R, Elsabagh, M, Obitsu, T, Sugino, T, Kurokawa, Y and Kawamura, K (2020) Effect of varying fermentation conditions with ensiling period and inoculum on photosynthetic pigments and phytol content in Italian ryegrass (Lolium multiflorum Lam.) silage. Animal Science Journal 91, e13309.10.1111/asj.13309CrossRefGoogle ScholarPubMed
Na, N, Qili, M, Wu, N, Sun, L, Xu, H, Zhao, Y, Wei, X, Xue, Y and Tao, Y (2022) Bacterial community and fermentation quality of ensiling alfalfa with commercial lactic acid bacterial additives. Frontiers in Microbiology 13, 836899.10.3389/fmicb.2022.836899CrossRefGoogle ScholarPubMed
Ogunade, IM, Jiang, Y, Pech Cervantes, AA, Kim, DH, Oliveira, AS, Vyas, D, Weinberg, ZG, Jeong, KC and Adesogan, AT (2018) Bacterial diversity and composition of alfalfa silage as analyzed by Illumina MiSeq sequencing: effects of Escherichia coli O157:H7 and silage additives. Journal of Dairy Science 101, 2048–2059.10.3168/jds.2017-12876CrossRefGoogle ScholarPubMed
Oliveira, AS, Weinberg, ZG, Ogunade, IM, Cervantes, AAP, Arriola, KG, Jiang, Y, Kim, D, Li, X, Goncalves, MCM, Vyas, D and Adesogan, AT (2017) Meta-analysis of effects of inoculation with homofermentative and facultative heterofermentative lactic acid bacteria on silage fermentation, aerobic stability, and the performance of dairy cows. Journal of Dairy Science 100, 4587–4603.10.3168/jds.2016-11815CrossRefGoogle ScholarPubMed
Playne, MJ and McDonald, P (1966) The buffering constituents of herbage and of silage. Journal of the Science of Food and Agriculture 17, 264–268.10.1002/jsfa.2740170609CrossRefGoogle Scholar
Polley, WH, Wilsey, BJ and Derner, JD (2007) Dominant species constrain effects of species diversity on temporal variability in biomass production of tallgrass prairie. Oikos 116, 2044–2052.10.1111/j.2007.0030-1299.16080.xCrossRefGoogle Scholar
Sangija, F, Martin, H and Matemu, A (2022) Effect of lactic acid fermentation on the nutritional quality and consumer acceptability of African nightshade. Food Science & Nutrition 00, 1–15.Google Scholar
Stackebrandt, E and Goebel, BM (1994) Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. International Journal of Systematic Bacteriology 44, 846–849.Google Scholar
Thomas, A (1977) An automated procedure for the determination of soluble carbohydrates in herbage. Journal of the Science of Food and Agriculture 28, 639–642.10.1002/jsfa.2740280711CrossRefGoogle Scholar
Turan, A and Onenc, SS (2018) Effect of cumin essential oil usage on fermentation quality, aerobic stability and in vitro digetibility of alfalfa silage. Asian-Australasian Journal of Animal Science 31, 1252–1258.10.5713/ajas.17.0834CrossRefGoogle ScholarPubMed
Tyrolová, Y, Bartoň, L and Loučka, R (2017) Effects of biological and chemical additives on fermentation progress in maize silage. Czech Journal of Animal Science 62, 306–312.10.17221/67/2016-CJASCrossRefGoogle Scholar
van Hooijdonk, T and Hettinga, K (2015) Dairy in a sustainable diet: a question of balance. Nutrition Reviews 73(Suppl 1), 48–54.10.1093/nutrit/nuv040CrossRefGoogle Scholar
Van Soest, PJ, Robertson, JB and Lewis, BA (1991) Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. Journal of Dairy Science 74, 3583–3597.10.3168/jds.S0022-0302(91)78551-2CrossRefGoogle ScholarPubMed
Wang, Q, Garrity, GM, Tiedje, JM and Cole, JR (2007) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology 73, 5261–5267.10.1128/AEM.00062-07CrossRefGoogle ScholarPubMed
Wang, S, Li, J, Zhao, J, Dong, Z, Dong, D and Shao, T (2021 a) Dynamics of the bacterial communities and predicted functional profiles in wilted alfalfa silage. Journal of Applied Microbiology 132, 2613–2624.10.1111/jam.15417CrossRefGoogle ScholarPubMed
Wang, L, Liu, Z, Jiang, H and Mao, X (2021 b) Biotechnology advances in β-carotene production by microorganisms. Trends in Food Science & Technology 111, 322–332.10.1016/j.tifs.2021.02.077CrossRefGoogle Scholar
Wróbel, B, Nowak, J, Fabiszewska, A, Paszkiewicz-Jasińska, A and Przystupa, W (2023) Dry matter losses in silages resulting from epiphytic microbiota activity – a comprehensive study. Agronomy 13, 450.10.3390/agronomy13020450CrossRefGoogle Scholar
Xu, Y, Hlaing, MM, Glagovskaia, O, Augustin, MA and Terefe, NS (2020) Fermentation by probiotic Lactobacillus gasseri strains enhances the carotenoid and fibre contents of carrot juice. Foods 9, 1803.10.3390/foods9121803CrossRefGoogle ScholarPubMed
Xu, D, Wang, N, Rinne, M, Ke, W, Weinberg, ZG, Da, M, Bai, J, Zhang, Y, Li, F and Guo, X (2021) The bacterial community and metabolome dynamics and their interactions modulate fermentation process of whole crop corn silage prepared with or without inoculants. Microbial Biotechnology 14, 561–576.10.1111/1751-7915.13623CrossRefGoogle ScholarPubMed
Yang, F, Zhao, S, Wang, Y, Fan, X, Wang, Y and Feng, C (2020) Assessment of bacterial community composition and dynamics in alfalfa silages with and without Lactobacillus plantarum inoculation using absolute quantification 16S rRNA sequencing. Frontiers in Microbiology 11, 629894.10.3389/fmicb.2020.629894CrossRefGoogle ScholarPubMed
Zhang, YC, Wang, XK, Li, DX, Lin, YL, Yang, FY and Ni, KK (2019) Impact of wilting and additives on fermentation quality and carbohydrate composition of mulberry silage. Asian-Australasian Journal of Animal Science 33, 254–263.10.5713/ajas.18.0925CrossRefGoogle ScholarPubMed
Zhang, X, Guo, XS, Li, FH, Usman, S, Zhang, YX and Ding, ZT (2022) Antioxidant, flavonoid, α-tocopherol, β-carotene, fatty acids, and fermentation profiles of alfalfa silage inoculated with novel Lactiplantibacillus plantarum and Pediococcus acidilactici strains with high-antioxidant activity. Animal Feed Science and Technology 288, 115301.10.1016/j.anifeedsci.2022.115301CrossRefGoogle Scholar
Zhao, GQ, Wei, SN, Liu, C, Kim, HJ and Kim, JG (2021) Effect of harvest dates on beta-carotene content and forage quality of rye (Secale cereale L.) silage and hay. Journal of Animal Science and Technology 63, 354–366.10.5187/jast.2021.e28CrossRefGoogle ScholarPubMed
Zong, C, Wu, Q, Wu, A, Chen, S, Dong, D, Zhao, J, Shao, T and Liu, Q (2021) Exploring the diversity mechanism of fatty acids and the loss mechanisms of polyunsaturated fatty acids and fat-soluble vitamins in alfalfa silage using different additives. Animal Feed Science and Technology 280, 115044.10.1016/j.anifeedsci.2021.115044CrossRefGoogle Scholar
Figure 0

Table 1. Chemical and biological composition of fresh alfalfa

Figure 1

Table 2. Chemical composition of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

Figure 2

Table 3. Fermentation features of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

Figure 3

Table 4. Microbial counts of alfalfa silage treated with propionic acid or a combination of propionic acid and squalene

Figure 4

Figure 1. The dynamic of loss rate of β-carotene in alfalfa silage treated with propionic acid or a combination of propionic acid and squalene. Loss rate of β-carotene (%) = (initial content − final content)/initial content × 100%; initial content refers to the β-carotene in fresh alfalfa; final content refers to the β-carotene after the corresponding ensiling time. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time; Values with different small letters show significant differences among additives in the same ensiling days (P < 0.05); Values with different capital letters show significant differences among ensiling days in the same additives (P < 0.05).

Figure 5

Figure 2. The α-diversity (a) and abundance heatmap (top 15, b) of fresh alfalfa and after ensiling. FA, fresh alfalfa; CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45, and 90 days; *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 6

Figure 3. The β-diversity based on principal coordinates analysis of bacterial communities in fresh alfalfa and alfalfa ensiled with different treatments for 3 days (a), 15 days (b), 45 days (c), and 90 days (d), respectively. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45, and 90 days.

Figure 7

Figure 4. Correlation analysis between bacterial communities (top 15 species) and fermentation characteristics in different alfalfa silages. CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; AA, acetic acid; LA, lactic acid; PA, propionic acid; BA, butyric acid; WSC, water-soluble carbohydrates; *, P < 0.05; **, P < 0.01; ***, P < 0.001.

Figure 8

Figure 5. The bacterial co-occurrence networks (on the species level) in different alfalfa silages. Nodes represent individual species; edges represent significant Spearman correlations (|ρ| > 0.5, P < 0.05). The size of each node is proportional to the number of connected edges, and the nodes are labelled by species and coloured by genus. The colour of the edges corresponds to a positive (red) or negative (blue) relationship; CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene.

Figure 9

Figure 6. The bacterial metabolic function prediction in fresh alfalfa and silages were analysed by PICRUSt2 (a), the change of C10-C20 isoprenoid biosynthesis and β-carotene biosynthesis before and after ensiling (b), and principal coordinates analysis of enzymes (geranylgeranyl pyrophosphate synthase, phytoene synthase, phytoene desaturase, lycopene β-cyclase) associated with bacterial β-carotene synthesis pathway in different alfalfa silages (c). Change (C10-C20 isoprenoid biosynthesis or β-carotene biosynthesis) = relative abundance (after ensiling) – relative abundance (before ensiling). CON, alfalfa ensiled without additive; PA, alfalfa ensiled with propionic acid; SQPA, alfalfa ensiled with a combination of propionic acid and squalene; 3, 15, 45 and 90, alfalfa was ensiled for 3, 15, 45 and 90 days; A, effect of additive; D, effect of ensiling time; A × D, interaction of additive and ensiling time.