Next Article in Journal
Effects of Lysine Cell Mass Supplementation as a Substitute for L-Lysine·HCl on Growth Performance, Diarrhea Incidence, and Blood Profiles in Weaning Pigs
Previous Article in Journal
Heart Rate and Heart Rate Variability Change with Sleep Stage in Dairy Cows
Previous Article in Special Issue
Evaluation of Nitrogen-Corrected Apparent Metabolizable Energy and Standardized Ileal Amino Acid Digestibility of Different Sources of Rice and Rice Milling Byproducts in Broilers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Different Kefir Source on Fermentation, Aerobic Stability, and Microbial Community of Alfalfa Silage

1
Department of Animal Science, Tekirdag Namık Kemal University, Tekirdag 59030, Turkey
2
Balikesir Directorate of Provincial Agriculture and Forestry, Republic of Turkey Ministry of Agriculture and Forestry, Balikesir 10470, Turkey
3
Department of Agricultural Biotechnology, Tekirdag Namık Kemal University, Tekirdag 59030, Turkey
*
Author to whom correspondence should be addressed.
Animals 2021, 11(7), 2096; https://doi.org/10.3390/ani11072096
Submission received: 25 June 2021 / Revised: 9 July 2021 / Accepted: 13 July 2021 / Published: 14 July 2021
(This article belongs to the Special Issue Feed Safety and Quality Control)

Abstract

:

Simple Summary

Minimizing silage additives cost while increasing silage quality is important for a sustainable livestock enterprise, especially in undeveloped and developing countries. In this study, therefore, commercially available kefir yeast (CK) and homemade kefir culture (HK), as a low-cost additive, was applied at untreated a common control (CON) and three different application doses (5.0, 5.7, and 6.0 log cfu g−1) on wilted alfalfa and evaluated with the fermentation characteristics and aerobic stability. The addition of HK with an application dose greater than 5.0 log cfu g−1 prevents mold formation and inhibits yeast counts in silages. Indeed, both CK and HK improve the silage quality and aerobic stability of alfalfa even with low water-soluble carbohydrate content.

Abstract

The present study has been one of the first attempts to thoroughly examine the effects of different kefir sources on fermentation characteristics, aerobic stability, and microbial communities of alfalfa silages. The effects of commercial kefir (CK) and homemade kefir culture (HK) applied with untreated a common control (CON) and three different application doses (5.0, 5.7, and 6.0 log cfu g−1) on wilted alfalfa and stored at an ambient temperature of 25–30 °C are studied. After 45 days of ensiling, fermentation characteristics and aerobic stability of silages were measured, and bacterial diversity was investigated by 16S ribosomal RNA gene sequencing using the GenomeLab™ GeXP platform. Both CK and HK accelerate more lactic acid production and reduced ammonia nitrogen concentration. Factor analysis of kefir sources suggests that the addition of kefir improves the aerobic stability of silages, even the initial water-soluble carbohydrate (WSC) content is inadequate via its antimicrobial effect on yeast and mold formation. Enterococcus faecium, Pediococcus pentosaceous and Lactobacillus brevis were dominant bacterial species among the treated groups at silo opening, while Lactobacillus plantarum and Lactobacillus brevis became dominant bacterial species after 7 days of aerobic exposure. In conclusion, the application of kefir on alfalfa silages improves fermentation quality and aerobic stability even with low WSC content.

1. Introduction

Alfalfa (Medicago sativa L.), one of the most common perennial forage legumes, plays a pivotal role in meeting the nutritional requirements of ruminants worldwide. It is mainly associated with high protein content and high nutritional quality, yielding high dry matter per acre and broader adaptation capability [1]. The most well-known application to minimize nutrient losses for future use of fresh forage crops is ensiling. Nevertheless, the main challenge faced by many researchers at ensiling process is the high buffering capacity (Bc) and low concentration of water-soluble carbohydrates (WSC), and high moisture content of alfalfa which results in undesirable secondary clostridial fermentation [2]. In such cases, the use of inhibitors or inoculants becomes mandatory to dominate microbial biota by lactic acid bacteria (LAB), which fermented the lysed plant membranes to lactic acid resulting in a lower pH of the ensiled material [3].
Commercial LAB has been used for many years as microbial additives owing to their inhibitory effect on undesirable microorganisms, e.g., Clostridia, Enterobacter and several other bacteria [4]. It is also well known from previous studies that LAB has a favorable contribution to the quality of the silage in terms of flavor and sensory profile, and preservation time of the final products [5,6]. However, as previously stated by Tao et al. [7], adding commercial LAB is one of the most frequently stated problems due to its high cost, especially in undeveloped and developing countries. Schnürer and Jonsson [8] were obviously right to draw our attention to the ingredients of an excellent starter culture for well-preserved silage, recommending a combination form of yeast and LAB; hence, kefir might be used as an alternative silage additive due to its complex symbiotic diversity of microorganisms, which has shown heterofermentative properties including LAB, yeast, and acetic acid bacteria [9]. Furthermore, kefir grains are recognized as a unique natural dairy starter culture that can be used several times, and kefir-derived products are also known as easily prepared, low-cost, and fascinating functional food with health benefits [10,11]. Kefir has also received considerable scholarly attention in recent years due to its antioxidant, anti-inflammatory, antifungal, antibacterial specialties, and various pharmaceutical attributes [12,13,14]. Previous research evaluating the effect of kefir microorganisms on aerobic spoilage of sorghum and wheat-straw silages also strengthens these assumptions [15,16]. However, there is a notable lack of high-quality research focusing specifically on fermentation characteristics, aerobic stability, and the microbial community of alfalfa silages inoculated with different kefir sources.
The aim of this paper was to evaluate the effectiveness of different kefir sources, which are produced commercially or homemade, and their various application doses onto first-cut alfalfa for improving fermentation characteristics and aerobic stability. The specific objective of this study was to determine the microbial community of silages by high-throughput sequencing methods after ensiling and aerobic exposure (AE) to explain the impact of isolated bacteria from alfalfa silages.

2. Materials and Methods

2.1. Forage and Silage Preparation

The plant material of the current study, alfalfa, was grown in an experimental plot of the Field Crops Department of Namık Kemal University (40.59° N and 27.34° S, Tekirdag, Turkey). The total precipitation, long-term mean precipitation, mean average temperature, and long-term mean average temperature of the experimental year was 299.0 mm and 581.5 mm, 15.7 °C and 14.0 °C, respectively. Since the total precipitation was not sufficient, forage water demands were supplied by irrigation fortnightly. First-cut alfalfa was harvested at the early blooming stage (10–20%) on 18 May 2019, by a forage harvester, wilted for 24 h, and manually chopped to approximately 1.5–2.0 cm in length. Wilted alfalfa had 304.6 g kg−1 dry matter (DM), 202.1 g kg−1 DM of crude protein (CP), 15.45 g kg−1 DM of WSC, 445 meq NAOH/kg DM of Bc, 7.50 of pH, 5.30 log cfu g−1 of LAB, and 8.08 log cfu g−1 of yeast and no mold before ensiling.
The silage was made in a laboratory-scale fermentation system: approximately 500 g wilted alfalfa was weighed and packed into polythene bags and then sealed by a vacuum sealer (CAS CVP−260PD). A commercial kefir yeast (MYStarter KF, contains Lactococcus lactis subsp. lactis biovar diacetylactis, Lactobacillus brevis, Leuconostoc mesenteroides subsp. mesenteroides ve Saccharomyces cerevisiae strains) and homemade kefir culture (contains Enterococcus faecalis, Lactobacillus brevis and Micrococcus luteus according to the 16S rRNA gene sequencing) was used for comparison. The silage treatments (each, 10 replicates) were designed as commercial kefir yeast (CF), and homemade kefir culture (HK) with an untreated common control (CON) and three different application doses (5.0, 5.7, and 6.0 log cfu g−1 of fresh matter) for 45 d of ensiling at an ambient temperature of 25–30 °C.

2.2. Chemical Analysis

After 45 d of ensiling, opened silos subsampled for microbial enumeration and aerobic stability determination. A representative 20 g wet silage or pre-ensiled material was taken and gently mixed in 180 mL of distilled water at room temperature for 1 h and then filtered through 4 layers of cheesecloths to determine ammonia nitrogen (NH3-N) and organic acid content. The pH of each silo was measured using silage extract with a pH meter (WTW-inoLab ph 730). The DM of samples was determined by drying at 60 ± 2 °C in an air-forced oven for 48 h, and DM loss was calculated via the weight differences between wilted alfalfa and opened silage samples. The nitrogen (N) content of wilted alfalfa was measured by the Kjeldahl method and multiplied by 6.25 to get the crude protein (CP) ratio by using AOAC methods. Samples were analyzed for NH3-N and WSC as previously reported by Anonymous [17] by using micro distillation and 0.2% anthrone reagent, respectively. The Bc of pre-ensiled alfalfa was determined by Playne and Mc Donald [18].
The organic acid content of silages (Acetic acid, AA; propionic acid, PA; Butyric acid, BA) was evaluated after deproteinization of silage extract with the metaphosphoric acid-formic acid mixture (3:1, v:v) according to the procedure described by Ulger et al. [19] by using a gas chromatograph (Shimadzu GC−2010+, Kyoto, Japan) with a capillary column (Restek, Bellefonte, PA, USA; 30 m, i. d.: 0.25 mm, f.t.: 0.25 μm), and with flame ionization detector (FID) over a temperature range of 45–230 °C. The lactic acid (LA) content of silages was determined using a spectrophotometric method previously described by Koc and Coskuntuna [20].

2.3. Microbial Populations

The LAB, yeast, and mold count was performed by using subsamples immediately after opening the silos and after the 7 d AE according to method described previously by Seale et al. [21]. While the pour plate method and MRS Agar (Merck, Darmstadt, Germany) were used to determine LAB for incubating anaerobically at 30 °C for 3 d, the spread-plate method and potato dextrose agar (Merck, Darmstadt, Germany) were used for yeast and mold enumeration after incubating aerobically at 28 °C for 5 d.

2.4. Aerobic Stability Analysis

Based on trapping CO2 gases into the KOH solution, the bottle system is one of the most common procedures for determining aerobic stability when reaching thermocouples is difficult. Generally, after 5 to 7 days of AE, pH, CO2, LAB, yeast, and mold content of silage are used to assess aerobic stability. In the current study, both the bottle system, previously described by Asbell et al. [22] and thermocouples (HOBO Pendant Temperature/Light 64K Data Logger, Onset Computer Corporation, Bourne, MA, USA) was used to record temperature with a 2 h interval during the 7 d of AE. Aerobic deterioration was considered when the temperature of the silage samples was 2 °C higher than the ambient temperature.

2.5. Microbial Diversity Analysis

After the proliferation of LAB in nutrient broth for 16 h, 1000 μL of aliquot was taken and centrifuged at 1000 g for 10 min. Microbial DNA from each silage sample was extracted after opening the silos and after the 7 d AE according to Liu et al. [23]. The 16S rDNA regions were amplified using primers F: 5′-AGAGTTTGATCCCTGGCTCAG−3′ and R: 5′-CCGTCAATTCCTTTGAGTTT−3′ [24]. PCR amplification reactions; 50 µL PCR volumes included: 10 ng rDNA, 1 µM of each primer, 1× PCR Buffer ((NH4)2SO4), 200 µM dNTP, 2.0 mM MgCl2 and 0.1U i-Taq™ DNA polymerase (5 U/mL) (iNtRON Biotechnology Inc., Burlington, MA, USA). The cycling protocol was 5 min at 95 °C for initial denaturation, 37 cycles of amplification; 95 °C for 45 s, 60 °C annealings for 60 s, 72 °C for the 40 s, and 10 min at 72 °C for final extension (Applied Biosystems ProFlex PCR System (Applied Biosystems, Foster City, CA, USA)). Afterward, the PCR products were run on 1.0% agarose gel using horizontal electrophoresis, and the gels were stained with SafeView™ Classic (Applied Biological Material Inc., Richmond, BC, Canada). PCR products were visualized under UV light in the Gel Documentation System.
The PCR products were sequenced using GenomeLab™ GeXp Genetic Analysis System (Beckman Coulter, Inc., Fullerton, CA, USA) after the precipitation with 3M NaAc. The chromatogram carefully checked the sequencing of the 16S rDNA region for overlapping nucleotide peaks by using ChromasPro Version 2.1.8 (Technelysium Pty. Ltd., South Brisbane, QLD, Australia). The checked sequences file consisting of MSTN fragments was controlled by the MEGA7 software (Molecular Evolutionary Genetics Analysis, version 7.0) [25]. The sequence data reported in this study were archived in The National Center for Biotechnology Information (NCBI) with the accession numbers MZ014989–MZ015001. The nucleotide sequences of studied 16S rDNA region in different species were performed from The National Center for Biotechnology Information webpage (https://www.ncbi.nlm.nih.gov/, accessed on 12 April 2021). The retrieved partial sequences were aligned by Clustal X, and the phylogenetic tree was generated by the Neighbor Joining (NJ) method (Kimura 2) in MEGA 7 software [25].

2.6. Statistical Analysis

Data were previously adjusted for the fixed effects of additive (CK and HK), dose (a common control, 5.0, 5.7, and 6.0 log cfu−1), and the interaction between these effects. The adjustment was made by analysis of variance using procedure PROC MIXED from SAS [26], considering the following statistical model:
Yijk = μ + ai + bj + (ab)ij + eijk
in which, Yijk is the value of measured characteristics; μ is a constant associated with each observation; ai is the effect of additive i; bj is the effect of dose j; (ab)ij is the interaction effect between the additive i and dose j; eijk is the random error of each observation. The differences among treatment means were tested using Tukey’s multiple range test, and significance was established at p < 0.05. Principal component analysis (PCA) was carried out according to the PRIN method of SAS by using the 19 variables and 1 supplementary variable (treatment group) and only PCAs with eigenvalues higher than 1 were retained and interpreted [19]. Furthermore, obtained PCAs were rotated in orthogonal Varimax rotation by FACTOR procedure of SAS, and only with an absolute loading value higher than 0.50 were considered to load on specific extracted PCAs (SAS, 2004).

3. Results

Kefir sources and their application doses did not change the DM, and WSC content of the silages (Table 1). However, comparing the pH and DM loss data reveals that both kefir source and doses significantly affect pH (p < 0.05), and the DM loss values (p < 0.001) of silages. In addition, the HK with an application dose of 6.0 log cfu g−1 (26.4 g kg−1 DM), CON (26.2 g kg−1 DM), and CK with an application dose of 5.0 log cfu g−1 (26.2 g kg−1 DM) presented higher DM loss. On the other hand, no significant results were obtained for WSC content among the silage groups (p > 0.05).
A comparison of the microbiological composition of alfalfa silages after silo openings reveals that the highest LAB count was observed in the CK group with an application dose of 6.0 log cfu g−1 (p < 0.01, Figure 1). Besides, from the figure above we can see that yeast count was affected significantly from different kefir sources (p < 0.001), while the mold count was not (p > 0.05).
From the data in Figure 2, it is apparent that the LA and BA concentrations of silages treated with different kefir sources and application doses were significantly affected (p < 0.001 and p < 0.05, respectively) while the AA and PA were not (p > 0.05). What stands out in this figure is the highest LA/AA ratio observed in the HK group with an application dose of 6.0 log cfu g−1, the lowest ratio observed in control. As shown in Figure 2, the results also indicate no statistically significant differences between the kefir sources and their application doses in terms of NH3-N concentration (p > 0.05).
The DM, pH, and CO2 values were similar among the kefir sources and their application doses, and no significant differences were observed up to 7 d of AE (Figure 3). On the other hand, the differences between the kefir sources and their application doses were significant for yeast (p < 0.05) and mold (p < 0.001). What is interesting about the data in this figure is that the HK with an application dose of higher than 5.0 log cfu g−1 prevents mold formation. Furthermore, a comparison of the data recorded by the data logger reveals that the highest aerobic stability was observed in CK with an application rate of 5.7 log cfu g−1 as >168 h while the lowest observed in HK with an application 5.0 log cfu g−1 as 32 h (Figure 3).
Particularly revealing is how principal factor analysis explains 86.75% of the total variability of original variables (Table 2). Loading vectors associated with the original variables, with eigenvalues higher than 1, are reported and used to interpret the meanings of 5 retained PCs. Biplot ordering using PCA of alfalfa silage characteristics inoculated with different kefir sources is also presented in Figure 4.
As shown in Figure 4, the first 2 components of PCA explain 55.1% of the total variation. Further analysis of the data reveals the variable loadings on PC1 were related to organic acids, such as LA, AA, BA, PA or fermentation characteristics such as NH3-N, or a microbiological composition such as LAB and yeast count (after AE), or DM and DM loss. PC1 had an eigenvalue of 6.428 and explained 33.83% of the total variability. While the CO2, pH (after AE), WSC, aerobic stability, and yeast count (after AE) loaded on PC2; yeast, DM loss, pH, and mold count (after AE) loaded on PC3 and explained 21.26% and 15.11% of the total variability, respectively. PC4 was characterized by LAB and mold count, and explained 9.92% of the total variability with an eigenvalue of 1.885. The last component, PC5, was characterized by NH3-N and pH explaining 6.64% of the total variability with an eigenvalue of 1.261.
We performed 16S rRNA sequencing to identify the microbial communities in the alfalfa silages systematically. The predominant bacterial species isolated at silo opening and after aerobic exposure were summarized in Table 3. Besides, the phylogenetic relationship used to estimate the relationship among various species based on genetic distances between the alfalfa silage treated with different doses of kefir sources was shown in Figure 5. As shown in Table 3, most of the bacteria were detected in alfalfa silage treated with different doses of kefir source. Furthermore, 16S rRNA sequence analysis indicates that while Enterococcus faecium, Pediococcus pentosaceous and L. brevis were dominant bacterial species among the treated groups at silo opening; Lactobacillus plantarum and L. brevis became dominant bacterial species after 7 d of AE.

4. Discussion

The present study was designed to determine fermentation quality and microbial community composition of alfalfa silage treated with different kefir sources and to find the most suitable dosage of use. Prior studies have noted the importance of an adequate substrate for LAB, DM, and WSC content to produce stable silages [4,7,27]. The DM content of wilted alfalfa (304.6 g kg−1 FW) and LAB concentration (5.30 log cfu g−1) before ensiling were consistent with the previous studies [28,29]. However, the WSC of pre-ensiled material (15.45 g kg−1 DM) was not adequate, considering the recommendation of 50 g kg−1 DM as a minimum required to ensure good fermentation during ensiling [1].
The major limitation of this study is the low WSC content (15.45 g kg−1 DM) of wilted alfalfa, which had been not sufficient to initiate lactic acid fermentation (Table 1). One of the most used parameters to determine the level of proteolysis in silages is NH3-N. Not significantly, but numerically, the NH3-N concentration of alfalfa silages was decreased with the addition of both different kefir sources (Figure 2). Therefore, it seems the Bc capacity of alfalfa plays a more significant role than the proteolysis during fermentation, considering the obtained results of this study. It is a well-known fact that proteolysis could account for more than half of the total N in alfalfa silage during fermentation and may result in efficient N utilization by ruminants [27]. Moreover, it is widely acknowledged that different LAB strains alone or combined with fibrolytic enzymes, sugar source, or organic acid can be used as additives to produce good quality alfalfa silage due to high Bc and low WSC concentration [30,31]. While there is limited information on proteolysis and NH3-N formation of kefir treated forage, in our study, the treatment with different kefir source and their various doses conserved more LA (HK produce more LA with an increasing application rate), and reduced NH3-N concentration.
Wang et al. [32] stated that one major factor that affects the extent of fermentation is silage pH. The current study found that kefir sources cannot be reduced forages’ pH (all above 5.0) to desired level alone when the initial WSC of forage was inadequate (Table 1). A possible explanation for these results may be the lack of adequate LAB growth in alfalfa silage by accelerating LA production to decrease pH during the fermentation. Prior studies have noted the importance of ensiling of legume forages with high WSC sources to produce more organic acids and obtain better fermentation quality, which are confirmed in our results [33,34].
It is expected that a rapid decrease in pH value with an increased dose of CK and HK inoculated silages. The LAB count was higher in treated silages than in the control silage except for HK with an application rate of 5.0 log cfu g−1 (p < 0.01, Figure 1). On the other hand, these results were not confirmed by silage pH. This may be explained by the fact that the conversion of LA to BA by undesirable microorganisms in silage, and the growth and metabolism of yeast which is found naturally in kefir fauna, utilized soluble carbohydrates at the initial fermentation stage and resulted in DM losses [35,36].
As previously stated by Kleinschmit and Kung [37], LA/AA ratio could be used as an indicator of an effective homolactic acid fermentation. In general, the study found a tendency for an increasing LA/AA ratio when the application rate was greater than 5 log cfu g−1 which means HK serves mainly homolactic characteristics whereas CK serves mainly heterolactic (Figure 2).
Previous research has shown that DM loss reflects the nutritive value and fermentation quality of silage [38,39]. There were no differences in DM contents between alfalfa silages treated with different kefir sources (p > 0.05). However, the lower DM loss was observed in silages treated with 5.7 log cfu g−1 of CK and 5.0 log cfu g−1 of HK. Furthermore, the silages treated with log 5.7 cfu g−1 of CK had the highest WSC content compared to other silages. As mentioned in the literature review, the amount of DM loss related to fermentation can vary depending on the dominant microbial species and fermented substrate types [38]. The homolactic LAB found in kefir’s microbial flora played an important role at these application rates to reduce DM loss by limiting plant respiration and undesirable microbial growth. A similar result was also found when chemical additives or homolactic LAB was used as silage additives [40]. Moreover, a tendency to increase DM loss with an increased application rate of HK was detected in this study. This may be related to the increased number of lactate assimilating yeasts in silage, resulting in higher DM losses and aerobic spoilage. A broadly similar point has also recently been made by Bai et al. [36].
Studies to improve feed-out stages of silages show the importance of aerobic stability, due to its potential cause of nutrient and DM loss, leads to health risks in animals and humans in terms of mycotoxins produced from undesirable microorganisms [7,41,42]. Several factors may influence the aerobic stability of silages, such as the size of the LAB population, the composition of inoculants, organic acid concentration, and the WSC content of ensiling material [7]. In this study, the highest aerobic stability was observed in CK with an application rate of 5.7 log cfu g−1 according to the indirect measurement of aerobic deterioration (Figure 3). Inconsistent with the expectation, the highest yeast count was observed in CK with an application rate of 5.7 log cfu g−1, possibly related to uncontrolled environmental conditions at the initial fermentation stage. On the other hand, the addition of HK with an application dose greater than 5.0 log cfu g−1 prevents mold formation and inhibits yeast counts in silages, suggesting that it improve aerobic stability. A possible explanation for this might be that metabolites produced from LAB, such as acetic acid, 1, 2-propanediol, and ethanol, improve aerobic stability [2,42]. Another possible explanation for this is that antimicrobial substances originated from alfalfa silages, e.g., saponin, during the fermentation process possess antimicrobial activity against bacteria and yeast, such as Bacillus subtilis, Candida albicans and S. cerevisiae [43].
The five silage quality components were related to organic acids, and thus, PC1 can be considered an organic acid-type factor. Surprisingly, a positive relationship was found between aerobic stability, pH (after AE), WSC, and yeast count (after AE). This finding was unexpected and suggested that the addition of kefir improves the aerobic stability of silages—even the initial WSC content was inadequate via its antimicrobial effect on yeast and mold formation (PC2, Table 2). The PC3 confirms that DM loss during ensiling is possibly related to the metabolism of yeasts, which utilizes WSC and produces ethanol [32]. The factor analysis may also support the negative relationship between LAB and mold during the fermentation process. Another negative correlation was detected between NH3-N and pH in PC5. In general, therefore, it seems that the pH of silages is closely related to the remaining acid concentration after the neutralization of organic acids by NH3-N in the silo. These results are similar to those reported by Bai et al. [44].
As shown in Table 1 and Figure 2, while comparable results were obtained in silages treated with CK and HK for LA content, only comparable results were obtained in silage treated with HK for pH after 45 d of fermentation. Many scholars hold the view that lactic acid-producing cocci, e.g., Pediococci, Streptococci, Enterococci, Lactococci and Leuconostocs, initiate LA fermentation in the early stages of ensiling and then are replaced by more acid-tolerant Lactobacilli, such as L. plantarum and L. brevis [2,45,46]. From 16S rRNA sequencing, it was indicated that E. faecium, P. pentosaceous and L. brevis were dominant bacterial species among the treated groups at silo opening (Table 3). The most obvious finding to emerge from the 16S rRNA sequence analysis is that synergetic effects of E. faecium and P. pentosaceous increased LA concentration and decreased pH value during the initial stages of ensiling. A possible explanation for these results may be the enzymatic hydrolysis of lignocellulosic biomass by the cellulolytic potential of E. faecium during the ensiling process. A similar LA accumulation and pH decline was also reported by Li et al. [47] for Pennisetum sinese (a kind of tropical perennial grass) silage.
In contrast, the monitored aerobic stability of CK with an application rate of 5.7 log cfu g−1 by indirect and direct methods was not confirmed with 16S rRNA sequence analysis. A possible explanation of these inconsistent results was mainly related to the growth of undesirable microorganisms, which results in aerobic deterioration. Although the predominant bacterial species were the same in such silage groups, the aerobic stability was different due to different residual WSC, LA, and AA concentrations found in silages.

5. Conclusions

This study set out to investigate the impact of alfalfa treated with different kefir sources and their various application doses on fermentation characteristics and aerobic stability. A key finding of the present study was that kefir improved the aerobic stability of silages by inhibiting the proliferation of yeast and mold even with the low WSC content. Factor analysis clearly shows that the antimicrobial activity of kefir has an inhibitory effect on yeast count during aerobic exposure. Considerably more work will need to be done to determine the effects of kefir with an adequate level of carbohydrate sources.

Author Contributions

Conceptualization, F.K. and S.E.; methodology, F.K., S.E., E.Ö.Ü. and R.I.; software, S.E. and R.I.; validation, B.O., E.Ö.Ü. and R.I.; formal analysis, B.O., E.Ö.Ü. and R.I.; investigation, F.K. and S.E.; data curation, E.Ö.Ü. and F.K.; writing—original draft preparation, S.E. and F.K.; writing—review and editing, S.E., F.K., E.Ö.Ü. and R.I.; visualization, S.E. and R.I.; supervision, F.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yuan, X.; Wen, A.; Dong, Z.; Desta, S.T.; Shao, T. Effects of formic acid and potassium diformate on the fermentation quality, chemical composition and aerobic stability of alfalfa silage. Grass Forage Sci. 2017, 72, 833–839. [Google Scholar] [CrossRef]
  2. Jiang, D.; Li, B.; Zheng, M.; Niu, D.; Zuo, S.; Xu, C. Effects of Pediococcus pentosaceus on fermentation, aerobic stability and microbial communities during ensiling and aerobic spoilage of total mixed ration silage containing alfalfa (Medicago sativa L.). Grassl. Sci. 2020, 66, 215–224. [Google Scholar] [CrossRef]
  3. Mc Garvey, J.A.; Franco, R.B.; Palumbo, J.D.; Hnasko, R.; Stanker, L.; Mitloehner, F.M. Bacterial population dynamics during the ensiling of Medicago sativa (alfalfa) and subsequent exposure to air. J. Appl. Microbiol. 2013, 114, 1661–1670. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, Q.; Yu, Z.; Wang, X.; Tian, J. Effects of inoculants and environmental temperature on fermentation quality and bacterial diversity of alfalfa silage. Anim. Sci. J. 2018, 89, 1085–1092. [Google Scholar] [CrossRef] [PubMed]
  5. Guo, X.S.; Bai, J.; Li, F.H.; Xu, D.M.; Zhang, Y.X.; Bu, D.P.; Zhao, L.S. Effects of malate, citrate, succinate and fumarate on fermentation, chemical composition, aerobic stability and digestibility of alfalfa silage. Anim. Feed. Sci. Technol. 2020, 268, 114604. [Google Scholar] [CrossRef]
  6. Zhao, S.S.; Wang, Y.P.; Yang, F.Y.; Wang, Y.; Zhang, H. Screening a Lactobacillus plantarum strain for good adaption in alfalfa ensiling and demonstrating its improvement of alfalfa silage quality. J. Appl. Microbiol. 2020, 129, 233–242. [Google Scholar] [CrossRef] [PubMed]
  7. Tao, L.; Zhou, H.; Zhang, N.; Si, B.; Tu, Y.; Ma, T.; Diao, Q. Effects of different source additives and wilt conditions on the pH value, aerobic stability, and carbohydrate and protein fractions of alfalfa silage. Anim. Sci. J. 2017, 88, 99–106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Schnürer, J.; Jonsson, A. Pichia anomala J121: A 30-year overnight near success biopreservation story. Antonie Van Leeuwenhoek 2011, 99, 5–12. [Google Scholar] [CrossRef]
  9. Chen, R.; Yang, S.; Zhang, L.; Zhou, Y.J. Advanced strategies for production of natural products in yeast. iScience 2020, 23, 100879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Chen, Z.; Liu, T.; Ye, T.; Yang, X.; Xue, Y.; Shen, Y.; Zhang, Q.; Zheng, X. Effect of lactic acid bacteria and yeasts on the structure and fermentation properties of Tibetan kefir grains. Int. Dairy J. 2021, 114, 104943. [Google Scholar] [CrossRef]
  11. da Silva, M.S.B.; Okura, M.H. Produtos à base de kefir desenvolvidos e estudados no Brasil. Res. Soc. Dev. 2021, 10, e19010716491. [Google Scholar] [CrossRef]
  12. John, S.M.; Deeseenthum, S. Properties and benefits of kefir—A review. Songklanakarin J. Sci. Technol. 2015, 37, 275–282. [Google Scholar]
  13. Kim, D.H.; Jeong, D.; Kim, H.; Kang, I.B.; Chon, J.W.; Song, K.Y.; Seo, K.H. Antimicrobial activity of kefir against various food pathogens and spoilage bacteria. Korean J. Food Sci. Anim. Resour. 2016, 36, 787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhang, J.; Zhao, X.; Jiang, Y.; Zhao, W.; Guo, T.; Cao, Y.; Teng, J.; Hao, X.; Zhao, J.; Yang, Z. Antioxidant status and gut microbiota change in an aging mouse model as influenced by exopolysaccharide produced by Lactobacillus plantarum YW11 isolated from Tibetan kefir. J. Dairy Sci. 2017, 100, 6025–6041. [Google Scholar] [CrossRef] [PubMed]
  15. Gonda, M.; Garmendia, G.; Rufo, C.; León Peláez, Á.; Wisniewski, M.; Droby, S.; Vero, S. Biocontrol of Aspergillus flavus in ensiled sorghum by water kefir microorganisms. Microorganisms 2019, 7, 253. [Google Scholar] [CrossRef] [Green Version]
  16. Okuyucu, B.; Esen, S. Is it possible to improve the fermentation and nutritional quality of wheat straw silage by replacing commercial inoculant with kefir? KSU J. Agric. Nat. 2021. [Google Scholar] [CrossRef]
  17. Anonymous. The Analysis of Agricultural Material; Reference Book; HMSO: London, UK, 1986; pp. 427–428. [Google Scholar]
  18. Playne, M.J.; Mc Donald, P. The buffering constituent of herbage and of silage. J. Sci. Food Agric. 1966, 17, 264–268. [Google Scholar] [CrossRef]
  19. Ulger, I.; Buyukkılıc Beyzi, S.; Kaliber, M.; Konca, Y. Chemical, nutritive, fermentation profile and gas production of citrus pulp silages, alone or combined with maize silage. S Afr. J. Anim. Sci. 2020, 50, 161–169. [Google Scholar] [CrossRef]
  20. Koc, F.; Coskuntuna, L. The comparison of the two different methods on the determination of organic acids in silage fodders. J. Anim. Prod. 2003, 44, 37–47. [Google Scholar]
  21. Seale, D.R.; Pahlow, G.; Spoelstra, S.F.; Lindgren, S.; Dellaglio, F.; Lowe, J.F. Methods for the Microbiological Analysis of Silage. In Proceedings of the Eurobac Conference, Uppsala, Sweden, 12–16 August 1986; pp. 147–164. [Google Scholar]
  22. Ashbell, G.; Weinberg, Z.G.; Azrieli, A.; Hen, Y.; Horev, B. A simple system to study the aerobic determination of silages. Can. Agric. Eng. 1991, 34, 171–175. [Google Scholar]
  23. Liu, B.; Huan, H.; Gu, H.; Xu, N.; Shen, Q.; Ding, C. Dynamics of a microbial community during ensiling and upon aerobic exposure in lactic acid bacteria inoculation-treated and untreated barley silages. Bioresour. Technol. 2019, 273, 212–219. [Google Scholar] [CrossRef]
  24. Beasley, S.S.; Saris, P.E.J. Nisin-producing Lactococcus lactis strains from human milk. Appl. Environ. Microbiol. 2004, 70, 5051–5053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. SAS. SAS/STAT® User’s Guide: Statistics (Version 9.1); SAS Institute Inc.: Cary, NC, USA, 2004. [Google Scholar]
  27. Li, P.; Ji, S.; Hou, C.; Tang, H.; Wang, Q.; Shen, Y. Effects of chemical additives on the fermentation quality and N distribution of alfalfa silage in south of China. Anim. Sci. J. 2016, 87, 1472–1479. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, Q.H.; Yang, F.Y.; Zhang, J.G.; Shao, T. Characteristics of Lactobacillus parafarraginis ZH1 and its role in improving the aerobic stability of silages. J. Appl. Microbiol. 2014, 117, 405–416. [Google Scholar] [CrossRef]
  29. Lynch, J.P.; Jin, L.; Lara, E.C.; Baah, J.; Beauchemin, K.A. The effect of exogenous fibrolytic enzymes and a ferulic acid esterase-producing inoculant on the fibre degradability, chemical composition and conservation characteristics of alfalfa silage. Anim. Feed. Sci. Technol. 2014, 193, 21–31. [Google Scholar] [CrossRef]
  30. Muck, R.E.; Nadeau, E.M.G.; McAllister, T.A.; Contreras-Govea, F.E.; Santos, M.C.; Kung, L., Jr. Silage review: Recent advances and future uses of silage additives. J. Dairy Sci. 2018, 101, 3980–4000. [Google Scholar] [CrossRef]
  31. Wang, B.; Yu, Z. Effects of moisture content and additives on the ensiling quality and vitamins changes of alfalfa silage with or without rain damage. Anim. Sci. J. 2020, 91, e13379. [Google Scholar] [CrossRef]
  32. Wang, C.; He, L.; Xing, Y.; Zhou, W.; Yang, F.; Chen, X.; Zhang, Q. Fermentation quality and microbial community of alfalfa and stylo silage mixed with Moringa oleifera leaves. Bioresour. Technol. 2019, 284, 240–247. [Google Scholar] [CrossRef] [PubMed]
  33. Zhang, Q.; Zhao, M.; Wang, X.; Yu, Z.; Na, R. Ensiling alfalfa with whole crop corn improves the silage quality and in vitro digestibility of the silage mixtures. Grassl. Sci. 2017, 63, 211–217. [Google Scholar] [CrossRef]
  34. Ni, K.; Zhao, J.; Zhu, B.; Su, R.; Pan, Y.; Ma, J.; Zhou, G.; Tao, Y.; Liu, X.; Zhong, J. Assessing the fermentation quality and microbial community of the mixed silage of forage soybean with crop corn or sorghum. Bioresour. Technol. 2018, 265, 563–567. [Google Scholar] [CrossRef] [PubMed]
  35. Ávila, C.L.S.; Carvalho, B.F.; Pinto, J.C.; Duarte, W.F.; Schwan, R.F. The use of Lactobacillus species as starter cultures for enhancing the quality of sugar cane silage. J. Dairy Sci. 2014, 97, 940–951. [Google Scholar] [CrossRef] [PubMed]
  36. Bai, J.; Xu, D.; Xie, D.; Wang, M.; Li, Z.; Guo, X. Effects of antibacterial peptide-producing Bacillus subtilis and Lactobacillus buchneri on fermentation, aerobic stability, and microbial community of alfalfa silage. Bioresour. Technol. 2020, 315, 123881. [Google Scholar] [CrossRef] [PubMed]
  37. Kleinschmit, D.H.; Kung, L., Jr. The effects of Lactobacillus buchneri 40788 and Pediococcus pentosaceus R1094 on the fermentation of corn silage. J. Dairy Sci. 2006, 89, 3999–4004. [Google Scholar] [CrossRef] [Green Version]
  38. Borreani, G.; Tabacco, E.; Schmidt, R.J.; Holmes, B.J.; Muck, R.E. Silage review: Factors affecting dry matter and quality losses in silages. J. Dairy Sci. 2018, 101, 3952–3979. [Google Scholar] [CrossRef] [Green Version]
  39. Guo, L.; Yao, D.; Li, D.; Lin, Y.; Bureenok, S.; Ni, K.; Yang, F. Effects of lactic acid bacteria isolated from rumen fluid and feces of dairy cows on fermentation quality, microbial community, and in vitro digestibility of alfalfa silage. Front. Microbiol. 2020, 10, 2998. [Google Scholar] [CrossRef] [PubMed]
  40. Auerbach, H.; Nadeau, E. Effects of additive type on fermentation and aerobic stability and its interaction with air exposure on silage nutritive value. Agronomy 2020, 10, 1229. [Google Scholar] [CrossRef]
  41. Schmidt, R.J.; Kung, L., Jr. The effects of Lactobacillus buchneri with or without a homolactic bacterium on the fermentation and aerobic stability of corn silages made at different locations. J. Dairy Sci. 2010, 93, 1616–1624. [Google Scholar] [CrossRef] [PubMed]
  42. Guo, X.S.; Undersander, D.J.; Combs, D.K. Effect of Lactobacillus inoculants and forage dry matter on the fermentation and aerobic stability of ensiled mixed-crop tall fescue and meadow fescue. J. Dairy Sci. 2013, 96, 1735–1744. [Google Scholar] [CrossRef] [Green Version]
  43. Wu, B.; Nishino, N. Identification and isolation of Lactobacillus fructivorans from wilted alfalfa silage with and without molasses. J. Appl. Microbiol. 2015, 120, 543–551. [Google Scholar] [CrossRef]
  44. Bai, J.; Ding, Z.; Ke, W.; Xu, D.; Wang, M.; Huang, W.; Zhang, Y.; Liu, F.; Guo, X. 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. 2021, 14, 1171–1182. [Google Scholar] [CrossRef] [PubMed]
  45. Ni, K.; Wang, F.; Zhu, B.; Yang, J.; Zhou, G.; Pan, Y.; Tao, Y.; Zhong, J. Effects of lactic acid bacteria and molasses additives on the microbial community and fermentation quality of soybean silage. Bioresour. Technol. 2017, 238, 706–715. [Google Scholar] [CrossRef] [PubMed]
  46. Xu, D.M.; Ke, W.C.; Zhang, P.; Li, F.H.; Guo, X.S. Characteristics of Pediococcus pentosaceus Q6 isolated from Elymus nutans growing on the Tibetan Plateau and its application for silage preparation at low temperature. J. Appl. Microbiol. 2018, 126, 40–48. [Google Scholar] [CrossRef] [PubMed]
  47. Li, J.; Yuan, X.; Desta, S.T.; Dong, Z.; Mugabe, W.; Shao, T. Characterization of Enterococcus faecalis JF85 and Enterococcus faecium Y83 isolated from Tibetan yak (Bos grunniens) for ensiling Pennisetum sinese. Bioresour. Technol. 2018, 257, 76–83. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Changes in lactic acid bacteria (LAB, (A)), yeast (B), and mold (C) count of alfalfa silage after 45 d of ensiling. CK: commercial kefir, HK, homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c, d) in each graph are statistically different (p < 0.05), NS: not significant, ** : p < 0.01, *** : p < 0.001.
Figure 1. Changes in lactic acid bacteria (LAB, (A)), yeast (B), and mold (C) count of alfalfa silage after 45 d of ensiling. CK: commercial kefir, HK, homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c, d) in each graph are statistically different (p < 0.05), NS: not significant, ** : p < 0.01, *** : p < 0.001.
Animals 11 02096 g001
Figure 2. Changes in ammonia nitrogen (NH3-N, (A)), lactic acid (LA, (B)), acetic acid (AA, (C)), butyric acid (BA, (D)), propionic acid (PA, (E)), and LA/AA ratio (F) of alfalfa silages after 45 d of ensiling. CK: commercial kefir, HK, homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c, d) in each graph are statistically different (p < 0.05), NS: not significant, * : p < 0.05, ** : p < 0.01, *** : p < 0.001.
Figure 2. Changes in ammonia nitrogen (NH3-N, (A)), lactic acid (LA, (B)), acetic acid (AA, (C)), butyric acid (BA, (D)), propionic acid (PA, (E)), and LA/AA ratio (F) of alfalfa silages after 45 d of ensiling. CK: commercial kefir, HK, homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c, d) in each graph are statistically different (p < 0.05), NS: not significant, * : p < 0.05, ** : p < 0.01, *** : p < 0.001.
Animals 11 02096 g002
Figure 3. Changes in dry matter (DM, (A)), pH (B), yeast (C), mold (D), carbon dioxide (CO2, (E)), and aerobic stability (F) of alfalfa silages after 7 d of aerobic exposure. CK: commercial kefir, HK: homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c) in each graph are statistically different (p < 0.05), NS: not significant, * : p < 0.05, *** : p < 0.001.
Figure 3. Changes in dry matter (DM, (A)), pH (B), yeast (C), mold (D), carbon dioxide (CO2, (E)), and aerobic stability (F) of alfalfa silages after 7 d of aerobic exposure. CK: commercial kefir, HK: homemade kefir, cfu: colony-forming units. The values with different letters (a, b, c) in each graph are statistically different (p < 0.05), NS: not significant, * : p < 0.05, *** : p < 0.001.
Animals 11 02096 g003
Figure 4. Biplot ordering using principal component analysis of alfalfa silage characteristics inoculated with different kefir sources. CK: commercial kefir, HK, homemade kefir, AA: acetic acid, AE: aerobic exposure, AS: aerobic stability, BA: butyric acid, CO2: carbon dioxide, DM: dry matter, DM loss: dry matter loss, LA: lactic acid, LAB: lactic acid bacteria, M: mold, NH3-N: ammonia nitrogen, PA: propionic acid, WSC: water-soluble carbohydrate, Y: yeast.
Figure 4. Biplot ordering using principal component analysis of alfalfa silage characteristics inoculated with different kefir sources. CK: commercial kefir, HK, homemade kefir, AA: acetic acid, AE: aerobic exposure, AS: aerobic stability, BA: butyric acid, CO2: carbon dioxide, DM: dry matter, DM loss: dry matter loss, LA: lactic acid, LAB: lactic acid bacteria, M: mold, NH3-N: ammonia nitrogen, PA: propionic acid, WSC: water-soluble carbohydrate, Y: yeast.
Animals 11 02096 g004
Figure 5. The phylogenic tree of the bacterial community is estimated using the NJ method (Kimura 2) that retrieved sequences provided by the NCBI GenBank database (MZ014989- MZ015001 from this study) (https://www.ncbi.nlm.nih.gov/, accessed on 12 April 2021). (a) comparison of sequences that retrieved from NCBI GenBank database and commercial kefir, (b) comparison of sequences that retrieved from NCBI GenBank database and homemade kefir.
Figure 5. The phylogenic tree of the bacterial community is estimated using the NJ method (Kimura 2) that retrieved sequences provided by the NCBI GenBank database (MZ014989- MZ015001 from this study) (https://www.ncbi.nlm.nih.gov/, accessed on 12 April 2021). (a) comparison of sequences that retrieved from NCBI GenBank database and commercial kefir, (b) comparison of sequences that retrieved from NCBI GenBank database and homemade kefir.
Animals 11 02096 g005
Table 1. Some chemical properties of alfalfa silage treated with different kefir sources after 45 d of ensiling.
Table 1. Some chemical properties of alfalfa silage treated with different kefir sources after 45 d of ensiling.
ItemCON
(n = 10)
CK (Log cfu g−1)HK (Log cfu g−1)SEMp-Value
5.0
(n = 10)
5.7
(n = 10)
6.0
(n = 10)
5.0
(n = 10)
5.7
(n = 10)
6.0
(n = 10)
DM284.5305.2292.3299.4272.5279.5315.49.540.151
DM loss 26.2 a26.2 a21.0c23.1 b20.9c25.6 a26.4 a0.21<0.001
pH5.60 ab5.70 ab5.85 ab5.50 ab5.90 a5.45 b5.55 ab0.080.024
WSC8.006.6211.687.678.9011.307.921.050.607
CON: control, CK: commercial kefir, HK, homemade kefir, DM: dry matter (g kg−1), WSC: water-soluble carbohydrate (g kg−1 DM), cfu: colony-forming units, SEM: standard error of the mean. The values with different letters (a, b, c) in each graph are statistically different (p < 0.05).
Table 2. Loading vectors of original variables after 45 d of ensiling alfalfa treated with different kefir source 1.
Table 2. Loading vectors of original variables after 45 d of ensiling alfalfa treated with different kefir source 1.
Original VariablePC 1PC 2PC 3PC 4PC 5
Yeast (after aerobic exposure)0.780 *0.514 *−0.0500.231−0.030
Butyric acid0.780 *−0.2940.023−0.1800.276
Propionic acid0.759 *−0.2710.1200.0030.382
Acetic acid 0.665 *−0.4300.4290.2720.077
Ammonia nitrogen0.513 *0.4270.4460.049−0.547 *
pH0.345−0.174−0.626 *0.262−0.538 *
pH (after aerobic exposure)0.4820.851 *−0.016−0.053−0.007
Carbon dioxide0.4400.846 *−0.156−0.0910.051
Aerobic stability−0.2130.750 *0.3400.447−0.061
Water soluble carbohydrate0.3750.733 *−0.288−0.0670.240
Mold0.1130.4160.010−0.630 *0.366
Mold (after aerobic exposure)0.223−0.3180.751 *0.4100.212
Dry matter loss−0.521 *0.0500.693 *−0.423−0.145
Lactic acid bacteria−0.619 *0.366−0.0880.507 *0.281
Dry matter (after aerobic exposure)−0.3940.3690.2280.4470.196
Yeast0.181−0.303−0.826 *0.3350.197
Lactic acid/Acetic acid−0.897 *0.095−0.151−0.218−0.037
Lactic acid −0.907 *0.176−0.206−0.0680.077
Dry matter−0.815 *0.0180.0680.2720.124
Eigenvalue6.4284.0382.8701.8851.261
Proportion (%)33.8321.2615.119.926.64
Cumulative (%)33.8355.0970.1980.1186.75
1 Values estimated by principal factor analysis after Varimax rotation of extracted PCAs * Variables with loading vectors higher than 0.50 were considered to load on specific PCA.
Table 3. 16SrRNA sequences isolated from ensiling alfalfa treated with different kefir sources.
Table 3. 16SrRNA sequences isolated from ensiling alfalfa treated with different kefir sources.
TreatmentsSilo OpeningAfter 7 d of Aerobic Exposure
CONLactobacillus brevisEnterococcus gallinarum,
Enterococcus casseliflavus,
Weissella paramesenteroides
CK1Pediococcus pentosaceus,
Enterococcus faecium
Lactobacillus plantarum,
Enterococcus faecalis
CK2Pediococcus pentosaceus,
Enterococcus faecium
Bacillus sp., Enterococcus faecalis
CK3Enterococcus faeciumLactobacillus plantarum,
Lactobacillus brevis
HK1Lactobacillus brevis,
Enterococcus faecium
Lactobacillus plantarum,
Lactobacillus brevis
HK2Pediococcus pentosaceus,
Enterococcus faecium
Weissella paramesenteroides,
Bacillus sp.
HK3Lactobacillus brevisLactobacillus plantarum,
Lactobacillus brevis
CON: Control, CK: commercial kefir, HK: homemade kefir, HK1: 5.0 log cfu g−1; HK2: 5.7 log cfu g−1, HK3: 6.0 log cfu g−1; CK1: 5.0 log cfu g−1, CK2: 5.7 log cfu g−1, CK3: 6.0 log cfu g−1.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koç, F.; Özkan Ünal, E.; Okuyucu, B.; Esen, S.; Işık, R. Effect of Different Kefir Source on Fermentation, Aerobic Stability, and Microbial Community of Alfalfa Silage. Animals 2021, 11, 2096. https://doi.org/10.3390/ani11072096

AMA Style

Koç F, Özkan Ünal E, Okuyucu B, Esen S, Işık R. Effect of Different Kefir Source on Fermentation, Aerobic Stability, and Microbial Community of Alfalfa Silage. Animals. 2021; 11(7):2096. https://doi.org/10.3390/ani11072096

Chicago/Turabian Style

Koç, Fisun, Emel Özkan Ünal, Berrin Okuyucu, Selim Esen, and Raziye Işık. 2021. "Effect of Different Kefir Source on Fermentation, Aerobic Stability, and Microbial Community of Alfalfa Silage" Animals 11, no. 7: 2096. https://doi.org/10.3390/ani11072096

APA Style

Koç, F., Özkan Ünal, E., Okuyucu, B., Esen, S., & Işık, R. (2021). Effect of Different Kefir Source on Fermentation, Aerobic Stability, and Microbial Community of Alfalfa Silage. Animals, 11(7), 2096. https://doi.org/10.3390/ani11072096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop