If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Koumiss is a type of famous fermented mare milk and considered an important nutritious beverage in central Asian countries. However, the production of koumiss cannot meet public demand in the market due to availability of mare milk. In the present study, 52 lactic acid bacteria and 20 yeast strains from traditional homemade Kazakhstan koumiss were isolated and identified. The isolates were used in a trial that included fermented cow milk, and the flavor profiles, color, and taste to determine their contribution in the co-fermentation of cow milk. Based on the sensory evaluation, KZLAB13 and KZY10 strains were selected as the best cofermentation combinations. The optimal fermentation conditions were confirmed as the ratio of the starter culture 2.4:1.6 % (vol/vol) KZLAB13 strain to KZY10 strain and a temperature of 36°C for 16 h using response surface methodology. After evaluating the quality of the optimized cow-milk koumiss compared with the Kazakhstan koumiss, results suggested that cow milk fermented by these 2 strains possessed a promising taste, flavor, and physicochemical and rheological properties. Altogether, our results showed that cow milk fermented with a combination of KZLAB13 and KZY10 strains can simulate the taste, flavor, and quality of traditional koumiss. Our study provided a novel alternative to mare-milk koumiss and could be used in dairy programs to fulfill the needs of people.
Koumiss (kumis, kymyz, or qymyz) is a fermented mare milk that is famous in central Asian countries. In Kazakhstan, koumiss is a national beverage commonly produced in wooden casks and containers by adding fresh mare milk to a wooden cask containing a small amount of koumiss as a starter culture at room temperature (
). The taste of koumiss is sour; it contains alcohol due to the presence of a high proportion of sugar in mare milk that ferments into lactic acid, alcohol, and other small molecules by microbial fermentation. It was reported that koumiss is fermented by the symbiosis of at least 2 main distinct microorganisms including lactic acid bacteria (LAB) and yeast (
), causing the unique flavor, texture, and acidity level, as well as health benefits of koumiss. Therefore, microorganisms contained in the koumiss are important for the quality of its fermentation (
). Beneficial effects, such as regulation of the alimentary canal activity, improved kidney function, and reduced blood cholesterol levels, also make koumiss a famous beverage in central Asia (
). In this respect, mare milk is very well known in ethnomedicine and is an important component of the traditional diet of Asian steppe people, among whom consumption of mare milk is culturally and environmentally rooted (
). Owing to some beneficial characteristics including low cholesterol levels and a rich amount of polyunsaturated fatty acids, mare milk in the human diet has received more attention in Europe, France, and, in particular, Germany (
The cost of mare milk has been a major factor for its use, and the availability falls short of the demand even in the traditional home countries of koumiss production. Moreover, it is hard to control the quality of koumiss made by traditional craftsmanship. Therefore, it is of extreme importance to exploit novel koumiss alternatives to meet daily life requirements of koumiss consumers (
). Cow milk is commonly used around the world for making dairy products, which is not only considered to be a low-cost dairy milk, but also rich in nutrient content that stimulates our metabolic activity and improves our health. On the other hand, cow milk is a good source of protein, which makes it highly efficient in the fermentation process. This suggests that cow milk may be the most promising mare-milk alternative to produce koumiss, particularly in developing countries.
Because of the growing demand for koumiss consumption numerous strains of LAB and yeast were isolated and identified from Kazakhstan koumiss in the present research. Additionally, 1 LAB and 1 yeast strain were selected for fermenting the cow milk to influence the taste of traditional koumiss by using response surface methodology (RSM). Further, volatile flavor substances formed during cow-milk fermentation were determined by GC-MS methodology and compared with Kazakhstan koumiss. The results of the present study will facilitate the production of koumiss by using cow milk as an alternative to mare milk.
MATERIALS AND METHODS
Collection of Samples
In this study, koumiss was collected from Kazakhstan that was produced by the local artisanal producers. The samples were collected in sterile bottles and transferred to the laboratory for microbiological analysis and stored aseptically at 4°C for further analysis.
Isolation of LAB and Yeasts and Microbiological Analysis
One milliliter of koumiss sample was added to a flask and mixed with 99 mL of sterile distilled water. Appropriate dilutions of 10−1 to 10−7 were carried out in the sterile water. The LAB strains were isolated on the deMan, Rogosa, and Sharpe (MRS) agar medium, and yeasts on the yeast extract peptone dextrose (YPD) agar media. Plates were incubated at 37°C for LAB and 28°C for yeast screening. The predominant, morphologically distinct, and well-isolated colonies of LAB were transferred to MRS agar plates and yeasts to the YPD agar plates for further enrichment.
All isolated strains from koumiss were morphologically characterized by cell morphology and Gram staining (
in: Starr M.P. Stolp H. Truper H.G. The Prokaryotes: A Handbook on Habitats, Isolation, and Identification of Bacteria. vol 2. Springer-Verlag,
New York, NY1981: 1711-1742
). Purified LAB strains were kept in MRS broth and yeasts kept in YPD broth supplemented with 20% (vol/vol) glycerol and stored at −40°C.
Identification of the Strains
The LAB strains were cultured in the MRS medium at 37°C for 24 h, after which cells of the strains were harvested by centrifugation and the genomic DNA of strains were extracted using a genome DNA isolation kit (Sangon Biotech, Shanghai, China). Further, the 16S rRNA gene was amplified using the universal primer 27F (5′-AGAGTTTGATCMTGGCTCAG-3′) and primer 1492R (5′-TACGGYTACCTTGTTACGACTT-3′). We performed PCR by applying the 2 × Taq Master Mixture (Shanghai, China). For amplification reaction, 1 mixture was prepared in a 200-µL tube that contained 6.5 µL of double-distilled (dd) H2O, 2 µL of primer 27F, 2 µL of 1492R (Table 1), 2 µL of genomic DNA, and 12.5 µL of 2 × Taq Master Mixture. The PCR amplification program was performed as described (denaturation at 94°C for 5 min, 35 cycles of denaturation 94°C for 1 min, 72°C for 3 min, and a final extension at 72°C for 10 min) by
and then run as follows. The amplified gene fragment was separated by gel electrophoresis using 1.5% (wt/vol) agarose gel and recovered using a gel extraction kit (Sangon Biotech). Amplified 16S rRNA genes were sequenced by Aoke Company (Yangling, China), and the 16S rRNA gene sequence was identified using the BLAST program developed by the National Center for Biotechnology Information (NCBI; https://blast.ncbi.nlm.nih.gov/Blast.cgi). Phylogenetic analysis was performed by using MEGA X with 100 bootstrap replicates (
For genomic DNA extraction, yeast cells were grown in YPD broth at 28°C on a turning shaker and collected by centrifugation. A genome DNA isolation kit (Sangon Biotech) was used to extract the genomic DNA of yeasts by following the method of
. Amplification of the 26S rRNA gene through PCR was done by using primers ITS1 (5′ TCCGTAGGTGAACCTGCGG 3′) and primer ITS4 (5′ TCCTCCGCTTATTGATATGC 3′). For the amplification reaction, 1 mixture was prepared in a 200-µL tube that contained 6.5 µL of dd H2O, 2 µL of primer ITS1, 2 µL of ITS4, 2 µL of genomic DNA isolated from koumiss, and 12.5 µL 2 × Taq Master Mixture. Amplification was performed in a PTC- 200 Peltier thermal cycler programmed as follows: 94°C for 3 min, then 35 cycles of 94°C for 30 s, 52°C for 30 s, and 74°C for 2 min, followed by 74°C for 10 min. The PCR product was separated by 1.5% (wt/vol) agarose gel electrophoresis, and the DNA was recovered by a DNA gel extraction kit (Sangon). Amplified 26S rRNA genes were sequenced by Aoke Company, and the 26S rRNA gene sequence was identified using the BLAST program. The sequences were used for the phylogenetic analysis using MEGA X with 100 bootstrap replicates (
First, all identified LAB strains (Supplemental Table S1, https://doi.org/10.3168/jds.2020-18527) from traditional koumiss samples were individually inoculated with the MRS medium incubated at 37°C for 24 h. The process was repeated several times umtil the viability of bacteria became stable. For milk inoculation, the tube containing the bacterial suspension was centrifuged (4,000 × g, 4°C, 3 min). Thereafter, cells were harvested and washed twice with 0.1M NaCl solution for milk fermentation. The milk fermentation was processed using low-fat (1.5%) cow milk. All samples were incubated at 37°C until the curd occurred (
Use of lactobacilli strains with probiotic potential in traditional fermented milk and their impact on quality and safety related to Listeria monocytogenes.
). At the end of incubation, the fermented cow milk was evaluated for sensory characteristics such as taste, color, and flavor. The sensory attributes of the samples were evaluated by the 10 panelists, and the strain that contributed to the highest sensory value was selected for further investigation. Second, to find the best combination of the LAB with yeast strain to ferment the cow milk, the previously selected LAB strain was inoculated with different isolated yeasts. All yeast strains were individually inoculated in YPD medium and incubated at 28°C for 24 h. For the selection of the best yeast strain, the same method has done as described above (Supplemental Table S2, https://doi.org/10.3168/jds.2020-18527). The total inoculum size was 6% (vol/vol), and the ratio between LAB and yeast was 1:1 (vol/vol).
Optimization of Cow-Milk Koumiss
One Factor at a Time Experimental Design.
The effects of different fermentation time (7, 10, 13, 16, 19, and 21 h), fermentation temperature (24, 28, 32, 36, 40, and 44°C), and total inoculation size [4, 1, 1, 4, 7, and 10% (vol/vol)] on the sensory value of fermented milk production were investigated by single-factor experiments. The ratios of KZLAB13 strain to KZY10 strain tested were 1.6:2.4, 0.4:0.6, 0.6:0.4, 2.4:1.6, 4.2:2.8, and 6:4% (vol/vol). Strains fermented at a constant temperature at 36°C for 16 h. The sensory evaluation was applied at the end of each process (
). The results of the preliminary screening could provide the required information for the RSM optimization experiment.
Response Surface Methodology: Box-Behnken Design.
To further obtain the optimal level of the main fermentation parameters and study their interactions on the quality of fermented products, the RSM method was applied. Based on the results of the one factor at a time experiment, the 3 most effective factors (X1, fermentation temperature; X2, inoculation size, and X3, fermentation time) were studied in 3-level codes: −1, 0, and +1. In the present study, we used the Box-Behnken design for the optimization of fermentation conditions of the studied fermented product. The design consisted of 17 randomly performed experiments, and the results related to the effect of the different factors on the development of koumiss from cow milk are shown in Table 2. The experiments were conducted in a 250-mL glass jar containing 100 mL of milk. The sensory quality of the fermented cow milk was evaluated as the responses for the studied factors. The quality of the fermented cow-milk koumiss was rated based on the sensory evaluation (Y1) as responded by the panelists. The optimal fermentation conditions were calculated by testing the empirical second-order polynomial regression model. The following second-order polynomial equation describes the relationship between the dependent and independent variables as follows:
where Y is the predicted response; β0 is the intercept; β1, β2, and β3 are the linear coefficients; β11, β22, and β33 are the squared coefficients; and β12, β13, and β23 are the interaction coefficients. Design Expert software (version 8.0.5, Stat-Ease Inc., Minneapolis, MN) was used for the experimental design and statistical analysis, while A, B, and C are coefficients.
Table 2The matrix of Box-Behnken design, analyzed using Design Expert 8 (version 8.0.5, Stat-Ease Inc., Minneapolis, MN)
). Initially, the Bradford solution was prepared and kept at 4°C. We mixed 100 µL of the sample with 5 mL of Bradford solution and incubated it at room temperature for 5 min. After, the absorbance value was measured at 595 nm. A standard solution of BSA (0, 0.0625, 0.125, 0.25, 0.5 and 1 g/L) was used to prepare a calibration curve.
Plate Count Method to Determine LAB and Yeast Population in Traditional Kazakhstan Koumiss and Cow-Milk Koumiss.
The serial dilutions of the samples were made and the suspension was spread onto MRS and YPD agar plates. The plates were then incubated at 37°C and 28°C, respectively. Growth was observed by counting the colonies, and the results were expressed as colony-forming unit per milliliter.
Measurement of pH and Total Acidity.
The pHSJ- 3F meter was used to determine the pH of the samples at room temperature (Leici, Shanghai, China). The total acidity of the fermented samples was determined by the titration method as suggested by
The method described by the Criterion for Sensory Evaluation of Fermented Milk (RHB 04-2020, issued by China's dairy industry standards, http://news.foodmate.net/2020/04/556069.html) was used to evaluate the sensory evaluations and the organoleptic characteristics of the samples. Ten trained panelists were selected for the evaluation of the sensory characteristics of the cow-milk koumiss. The samples were served at 4°C in cups and labeled with numbers. The panel was asked to provide individual scores on a 10 to 30 scale for color (white, milky white, and slightly yellow), a 10 to 40 scale for flavor (light sour smell, sour smell, and bad smell), and a 10 to 30 scale for taste (millet taste, slightly millet taste, and bitter taste). Three replicates were used for each sample in each session, and the mean for each sample was calculated (
The rheometer (AR1000, TA instruments, New Castle, DE) with a 20-mm parallel plate was used for the identification of rheological properties of Kazakhstan koumiss, cow-milk koumiss, and milk samples. The replicate measurements were taken independently, and the RheoWin Pro software package (version 2.94, Thermo Haake, Karlsruhe, Germany) was used for data analysis. Strain sweep (0.01–100% at 1 Hz) was applied to test the linear viscoelastic region of the samples. Further, the frequency dependence of storage modulus (G′) and loss modulus (G″) was determined by a frequency sweep (0.1–10 Hz at 1% strain;
Gas Chromatographic Analysis for Determining the Ethanol Content
About 3 g of Kazakhstan koumiss and cow-milk koumiss samples were weighed and homogenized with a blender, and centrifuged at 4,000 × g for 5 min. Thereafter, the supernatant was filtered through a 0.22-µL disposable syringe membrane filter (Xian, China) and subsequently injected directly into the gas chromatographic apparatus. The reagent ethanol was used as a standard. A diluent-standard solution was prepared by diluting ethanol to 0.05% (vol/vol) with distilled water. The diluent solution was stoppered and stored at −20°C. The ethanol content in the samples was measured by gas chromatography (Perkin Elmer Headspace sampler HS 40XL, Perkin Elmer Auto system Gas Chromatograph equipped with Ionization Detector FID, Akron, OH) as previously used by
Volatile organic compounds (VOC) from Kazakhstan koumiss, cow-milk koumiss, and milk were analyzed using the headspace solid phase microextraction technique (
). Five milliliters of each sample was placed into a 20-mL headspace vial containing 2 g of NaCl and 20 µL of n-hexadecane as the internal standard; after that, vials sealed rapidly. A 50/30 µm DVB/CAR/PDMS fiber (Supelco, Bellefonte, PA) was used to carry the solid phase microextraction procedure. The GC-MS procedure was accomplished by using 7890A-5975C (Agilent Technologies, Santa Clara, CA). A capillary column (HP-5 ms, 30 m × 250 μm × 0.25 μm, Agilent Technologies) was used with helium as a carrier gas at 1.0 mL/min flow rate. The oven temperature was set at 40°C for 3 min, increased to 120°C at a rate of 4°C/min, and then increased further to 240°C at a rate of 6°C/min for 12 min. A final 3-min extension was performed at 250°C. The ion source and the transfer line temperatures were set at 230°C and 250°C, respectively. The National Institute of Standards and Technology (NIST) mass spectral search program (version 2.0; https://www.nist.gov/srd/nist-standard-reference-database-1) was applied to evaluate the results by comparing their mass spectra. The quantitative relative concentration of each volatile compound was calculated by using the internal standard for all samples. To calculate the concentration of the analyte in the sample, the peak area of the analyte in the sample was divided by the peak area of n-hexadecane and multiplied by the concentration of n-hexadecane.
Statistical Analysis
All experiments were carried out in triplicate, and all analyses were carried out once. For optimization of fermentation parameters, the statistical software package Design Expert (Stat-Ease, Inc.) was used, and 3-dimensional response surface plots were produced to elucidate the relationships between the responses and the experimental levels of each independent variable. The statistical competence of the model was determined through ANOVA, presented in Table 3. The F- and P-values were used to test the significance of each factor and regression model. The excellence of the polynomial model equation was determined statistically through the coefficient of determination (R2) and the adjusted R2. Significant differences between the mean values in physicochemical characteristics of Kazakhstan koumiss and cow-milk koumiss were determined by Student's t-test using the Statistix software version 10 (https://www.statistix.com/). For VOC of Kazakhstan koumiss and cow-milk koumiss, multiple comparisons were performed using one-way ANOVA, followed by posthoc testing using Tukey's multiple comparisons test. A value of P < 0.05 was regarded as statistically significant. Data are presented in Table 4, Table 5 as the mean ± standard deviation of at least 3 independent experiments.
Table 3ANOVA of the quadratic model for the sensory evaluation, analyzed using Design Expert 8 (version 8.0.5, Stat-Ease Inc., Minneapolis, MN)
All data were expressed as the mean ± SD of 3 independent experiments. Statistical analyses were carried out by Tukey's multiple comparison test (P < 0.05). The author analyzed experiment data using Statistix version 10 (https://www.statistix.com/).
Different superscript letters in the same row indicate significant differences.
a–c Different superscript letters in the same row indicate significant differences.
1 All data were expressed as the mean ± SD of 3 independent experiments. Statistical analyses were carried out by Tukey's multiple comparison test (P < 0.05). The author analyzed experiment data using Statistix version 10 (https://www.statistix.com/).
). Therefore, MRS and YPD agar plates were used to isolate LAB and yeast strains from Kazakhstan koumiss. All 52 LAB isolates were gram-positive, mesophilic, non-spore–forming, catalase- and oxidase-negative, and rod-shaped, corresponding well with the classical characteristics of Lactobacillus given in Burgey's manual of Systematic Bacteriology (
). The 16S rRNA sequences were used as query sequences to BLAST in the GenBank DNA database (https://blast.ncbi.nlm.nih.gov/Blast.cgi), and the sequences with 98 to 100% identities were selected for the phylogenetic analysis, and all the identified isolates are presented in Figure 1. These results were relatively consistent with
A total of 20 yeast clones were isolated on the YPD agar plate. We sequenced the 26S rRNA genes from 20 isolates, and the similarity was checked by comparing the nucleotide sequence database NCBI in GenBank (https://blast.ncbi.nlm.nih.gov/Blast.cgi). We noted that the obtained sequences showed high similarity, ranging from 98 to 100% to reference sequences. The identified isolates are presented in Figure 2. It has been stated that yeasts are involved in the production of some fermented milk beverages in many countries around the world (
). Yeast is one of the major components of the microbiota of various fermented milk and some other traditional milk products such as fermented milk of camel, goat, cow, and yak (
). It has been investigated that the qualities of fermented milk may be mainly affected by the starter cultures because of its beneficial metabolites (such as some organic acids, exopolysaccharides, and flavor substances). For example, exopolysaccharides produced by some LAB could provide a better texture of yogurt and lower whey separation (
). The LAB strains produce several compounds during lactic acid fermentations such as hydrogen peroxide, ethanol, organic compounds, exopolysaccharides, diacetyl, bacteriocins or bactericidal proteins, and other several enzymes. These compounds are important for various characteristics, such as color, flavor, and aroma, whereas LAB could control the functional characteristics that include fibrinolytic activity, probiotic effects, and antioxidant activity of fermented food (
). In addition, it was further postulated that LAB could affect the growth of filamentous fungi and yeasts that may alter the formation of wine flavor due to the esterification of ethanol and lactic acid or other acids produced by yeasts and filamentous fungi (
Combinations of the Selected LAB and Yeast Strains
To further confirm which LAB strains are the most important bacteria species in Kazakhstan koumiss fermentation, sensory analysis of a single-bacteria fermented milk was conducted. Among all Lactobacillus strains, KZLAB13 strain was selected as the best and used for further investigation. The KZLAB13 strain is in the group with Lactobacillus delbrueckii ssp. delbrueckii DSM 20074 JCM 1012 reference strain and has 93% similarity. According to the results, the milk samples fermented by KZLAB13 strain showed good appearance and a pleasant odor after evaluating by well-trained panelists (Supplemental Table S1, https://doi.org/10.3168/jds.2020-18527). Similarly, L. delbrueckii ssp. bulgaricus is the key organism used as a starter for the preparation of yogurt and may play a major role in the coagulation of milk by rapidly increasing the acidity (
). Also, it has been reported that proteolytic strains of L. delbrueckii ssp. bulgaricus could hydrolyze proteins that may alter the viscosity and trigger the firmness of the yogurt samples (
Influence of addition of proteolytic strains of Lactobacillus delbrueckii ssp. bulgaricus to commercial ABT starter cultures on texture of yoghurt, exopolysaccharide production and survival of bacteria.
). Therefore, the KZLAB13 strain was selected for cofermentation with yeast strains isolated from koumiss. As for the selection of the best yeast strain, sensory evaluation of the sample was conducted after fermenting by each yeast strain combined with the KZLAB13 strain. The sensory score was remarkably higher for the KZLAB13 strain together with the KZY10 strain combination compared with other combinations (Supplementary Table S2, https://doi.org/10.3168/jds.2020-18527). The KZY10 strain is in the same group with Kluveromyces marxianus E4 reference strain and has 88% similarity. Interestingly, K. marxianus was also found to be the most prevalent species in Asian fermented milk, and it was a kind of thermotolerant yeast present in the kefir (
). Kluveromyces marxianus has been considered as one of the most important starter cultures of fermentation due to its attractive organoleptic profile and rapid growth. Furthermore, K. marxianus is an acid tolerant species and can produce a lower level of ethanol. Therefore, K. marxianus species play an important role in food fermentation (
). Thus, KZLAB13 and KZY10 strains were selected as starter cultures for further analysis.
One Factor at a Time Experimental Design
The effect of fermentation parameters such as time, temperature, and inoculum size on the sensory evaluation of the fermented cow milk by the mixture of koumiss microorganisms such as KZLAB13 and KZY10 strains was investigated using a single-factor experiment (Figure 3). It indicated that fermentation time, fermentation temperature, and inoculum size of a starter culture greatly influenced the quality of cow-milk koumiss. The previous study found that the volatile compounds involved in the sensory characteristics were influenced by the inoculum size, which also significantly affects the overall fermentation results (
). We observed that the sensory evaluation scores were significantly increased with the prolonging of the fermentation time from 7 to 16 h for the fermented product; however, a reduction was noted at 19 and 21 h time points. Moreover, the highest sensory evaluation scores were detected at the 16 h time point, whereas the lowest sensory evaluation scores were found at 7 and 21 h, respectively (Figure 3A). The highest sensory evaluation scores were observed at 16 h, suggesting the optimum fermentation time for the fermented product. Furthermore, we studied the influence of fermentation temperature (24, 28, 32, 36, 40, and 44°C) on sensory evaluation of the cow-milk koumiss (Figure 3B). The results showed that fermentation temperatures significantly affected the sensory evaluation of the cow-milk koumiss. The gradual increase was noted for sensory evaluation of the cow-milk koumiss at 24, 28, and 32 °C, and at 36°C cow-milk koumiss had the highest sensory evaluation score compared with other temperatures, whereas gradual decline was detected after 36°C. At 36°C, cow-milk koumiss showed the highest sensory score, indicating the optimum temperature for further optimization. Previously, it has been reported that fermentation temperature also plays an important role in the walnut milk beverage fermentation (
). Therefore, the optimum inoculum size of the starter culture was 4% (vol/vol) for every 100 mL of milk because of its highest sensory evaluation score discovered (Figure 3C). Recent studies postulated that sensory evaluation is directly proportional to the fermentation time. The fermentation time triggered the acidity of the fermented food (
Figure 3Effect of fermentation time (A), fermentation temperature (B), and inoculation size (C) on the sensory evaluation of the cow-milk koumiss. Data represent means ± SD (n = 10 for sensory evaluation scores).
Response Surface Methodology to Optimize the Fermentation Parameters
In the present study, an overall second polynomial equation by multiple regression analysis was developed for the sensory evaluation as represented below in terms of code values:
[2]
where Y is the sensory evaluation, X1 is inoculum size, X2 is fermentation temperature, and X3 is fermentation time, respectively.
In Table 3, the F-value (30.30) of the model indicates that the regression model was statistically significant (P < 0.01). The R2 (0.9750) and adjusted R2 (0.9428) suggested that the regression equation developed has a goodness of fit and could successfully predict the response and explain more than 95% of the variability in sensory evaluation. The P-value for lack of fit (0.1033) implied the lack of fit was not significant compared with the pure error. Furthermore, the coefficient of variation (6.92%) indicated a high degree of precision and better reliability of the experimental results. These results suggested that the model is reliable. The smaller the magnitude of the P-value, the more significant the corresponding term was. The most significant factors of this model were X1,
, and
, with P-values less than 0.001 each; X1, X2, X3, and
are also significant. Information visualization of the relationship between variables and response were carried out through 3-dimensional response surfaces and count plots based on the result of the second-order polynomial function model (Figure 4). As shown in Figure 4A and 4D, inoculum size (X1) and temperature (X2) displayed a significant effect on the sensory evaluation of the cow-milk koumiss. This result suggests that inoculum size and fermentation temperature might be the crucial factors to improve the sensory evaluation of cow-milk koumiss. The mutual interaction between temperature and time (X3) was significant, as the shape of the contour looks like an oval, which is consistent with the results of ANOVA analysis (Figure 4C). On the contrary, circular shapes in Figure 4A indicate that the synergistic effects of inoculum size and temperature were not apparent. Based on the regression equation, the optimal fermentation parameters for the maximum score of the sensory evaluation are the fermentation temperature at 36°C, the inoculum size of starter culture at 4%, and the fermentation time of 16 h.
Figure 4Normal plots and 3-dimensional surface plots indicating the interaction effects of independent variables on sensory evaluation.
Physicochemical Characteristics of Kazakhstan Koumiss and Cow-Milk Koumiss
The average physicochemical and microbiological properties of Kazakhstan koumiss and cow-milk koumiss are shown in Table 4. All data compared with Kazakhstan koumiss and values of acidity, pH, ethanol, protein, yeasts, and LAB contents of cow-milk koumiss significantly varied from Kazakhstan koumiss. The pH and acidity of cow-milk koumiss were 4.3 and 71.6 degrees of titratable acidity, respectively. The cow-milk koumiss had a pH of ~4.3 due to lactic acid production, and another study found similar results for koumiss after inoculation with starter cultures at 25°C from modified bovine milk (
). The koumiss was classified in groups varying on their pH. According to this classification of koumiss as strong (pH 3.3–3.6), moderate (3.9–4.5), and light (4.5–5.0), cow-milk koumiss can be fitted to the moderate koumiss group. The titratable acidity of the Kazakhstan koumiss sample is slightly higher compared with those in the cow-milk koumiss. Our result is in line with the findings reported by
. Based on the previous literature, we observed that several factors affected the acidification parameters. This may be because of the differences in raw material and fermentation procedure. Cow-milk koumiss had higher protein content than the Kazakhstan koumiss. It was noted that the protein contents were 1.9 and 2.7% for Kazakhstan koumiss and cow milk koumiss, respectively, similar to the previous findings by
. Another study also revealed that protein content depends on the milk source; for example, the protein level in cow milk is higher than mare milk, whereas the ratio of casein to whey protein (1:1) is quite similar to human milk (
). The number of microorganisms in Kazakhstan koumiss was 2.24 × 108 (cfu/mL) of LAB and 2.41 × 106 (cfu/mL) of yeast, and in the cow-milk koumiss, 9.9 × 106 (cfu/mL) of LAB and 1.14 × 106 (cfu/mL) of yeast, showing that in Kazakhstan koumiss LAB and yeasts accounted for a high ratio.
Quantification of the VOC
Traditional koumiss is known to contain a large number of LAB and yeast strains. Previous researchers identified more than 50 LAB strains from koumiss (
Isolation and identification of lactic acid bacteria from koumiss in Eastern Inner Mongolia of China. In Proc. International Conference on Materials Science, Resource and Environmental Engineering, Xi'an City, Shaanxi, China. AIP Conference Proceedings 1794:050005.
). Such microbiological biodiversity due to spontaneous fermentation could explain the high number of potential volatile compounds in typical beverages compared with our samples obtained after the optimized fermentation of cow milk by the koumiss starter culture. Table 5 shows the VOC, whose quantity and presence significantly differ from the reference Kazakhstan koumiss sample. Their mean abundance and concentration were quantified, and a total of 41 VOC were detected in the present study. One of the most important compounds for koumiss is ethanol. Compared with Kazakhstan koumiss, which contained 22.98 µg/L of ethanol, our samples of the cow-milk koumiss had a low quantity of alcohol. Indeed, the high level of ethanol usually found in traditional koumiss is usually due to yeast (
Volatile organic compounds in naturally fermented milk and milk fermented using yeasts, lactic acid bacteria and their combinations as starter cultures.
). However, our cow-milk koumiss contained ethanol at a concentration of 20.83 µg/L. This could be due to certain strains of bacteria and because, as mentioned before, K. marxianus produces less alcohol. Additionally, ethanol production could be possible in relatively high quantities with Lactobacillus strains (
Volatile organic compounds profile during milk fermentation by Lactobacillus pentosus and correlations between volatiles flavor and carbohydrate metabolism.
). It is also explained by the fact that mare milk has more sugar than cow milk, and during the fermentation process, sugar converts to alcohol by the starter culture. Some reports investigated that 3-methylbutanoic acid and its corresponding alcohols and aldehydes, namely 3-methylbutanol and 3-methylbutanal, are AA degradation products. It has been identified that the oxidation of 3-methylbutanal can produce 3-methylbutanoic acid (
). The value for 3-methylbutanoic acid was 0.32 µg/L, suggesting that the oxidation of 3-methylbutanal may lead the production of 3-methylbutanoic acid in koumiss samples. The results of our study found a maximum value of 4.7 and 1.25 µg/L for 3-methylbutanol in Kazakhstan koumiss and cow-milk koumiss, respectively. The use of cow milk as a matrix for inoculation could also explain that only 19 compounds from Kazakhstan koumiss were common with the VOC in cow milk koumiss. Certainly, there is an important difference in lactose content between mare milk (60–80 g/L) and cow milk (30–50 g/L); consequently, LAB strains could produce a higher compound variability. The most discriminating molecule oxime methoxy-phenyl (OMP) is a chemical product belonging to the imine group, which has sweetener taste. Generally, OMP is present in large quantities in UHT cow milk (
). In the present study, OMP was significantly higher (P < 0.01) in the milk sample than compared with Kazakhstan koumiss and cow-milk koumiss. It was observed that acetaldehyde is critical for the aroma of yogurt type milk, and thus acetaldehyde has received more attention among all secondary metabolites (
). Herein, we detected that the acetaldehyde contents were higher in the cow-milk koumiss. This higher concentration may be due to a component of the starter isolated from Kazakhstan koumiss, known as KZLAB13 strain. The previous study observed that compared with other LAB, L. delbrueckii strain shared markedly high acetaldehyde production activity in the kefir starter (
). We identified a total of 11 acid compounds among the studied samples. The most abundant acidic compound was acetic acid, followed by octanic acid and butanic acid, respectively. It was reported that acetic acid is the major subproduct of LAB fermentation (
). The range of acetic acid concentration was 0.33 to 2.87 µg/L, where the highest value was detected for Kazakhstan koumiss samples. Furthermore, the butanoic acid and octanoic acid concentrations were only observed for Kazakhstan koumiss samples. Previous studies found that these important acid compounds are found abundantly in dairy foods as well as in fermented milk (
Volatile organic compounds profile during milk fermentation by Lactobacillus pentosus and correlations between volatiles flavor and carbohydrate metabolism.
). We noted that the value of hexanal concentration reached up to 0.024 μg/L in the Kazakhstan koumiss, whereas this value was 0.22 μg/L in cow-milk koumiss, which is significantly lower compared with Kazakhstan koumiss samples. In line with above, recent studies showed that hexanal is an important compound for flavor that can be produced by the oxidation of unsaturated fatty acids in dairy products such as fermented milk (
In situ real-time monitoring of volatile metabolites of fermented milk by dynamic headspace sampling-atmospheric pressure ionization mass spectrometry.
). A total of 6 ketone compounds were detected in the Kazakhstan koumiss, cow-milk koumiss, and milk samples. Among them, some are essential flavor VOC such as 2-heptanone and 9-octadecanone, which were present in considerably high amounts in the studied samples. The concentration of these compounds was 0.3 µg/L, which was higher than compared with cow-milk koumiss.
Rheological Properties
According to the previous studies, the apparent viscosity is considered to be the function of the shear rate, and both (apparent viscosity and shear rate) are inversely proportional to each other. The increase in the shear rate decreases the apparent viscosity, suggesting the pseudoplastic pattern for all the fermented products (
). We observed that the apparent viscosity of cow-milk koumiss was higher than the Kazakhstan koumiss at the beginning of shearing. The casein to whey protein ratio and casein concentration were lower in the mare milk as compared with the cow milk. Furthermore, we noticed clear differences in the casein fractions and the size of casein micelles between mare milk and cow milk. According to earlier studies, the casein of mare milk possesses an equal amount of β- and αS-casein, while αS-casein concentration is slightly higher than β-casein concentration in cow milk (
). These distinguishing properties of cow and mare-milk casein may play an important role in the differentiation of rheological properties of the mare and cow milk. A similar pattern for casein between cow and mare milk was also reported in the present study. The frequency dependence of storage modulus (G′) and loss modulus (G″) was determined by using frequency sweeps at 1% strain (Figure 5). We observed that the apparent viscosity of cow milk koumiss was higher than the Kazakhstan koumiss at the beginning of shearing (Figure 6). The dynamic testing provides useful information on the viscoelastic properties of yogurts (i.e., the G′ and G″), which denote the degree of elastic and viscous behavior, respectively. For all samples, the response was typical for yogurt gels with G′ higher than G″.
Figure 5Frequency sweeps of Kazakhstan koumiss, cow-milk koumiss, and milk at 1% strain. G′ = storage modulus; G″ = loss modulus.
The present study indicates that cow-milk koumiss fermented by a combination of KZLAB13 and KZY10 strains isolated from the traditional homemade Kazakhstan koumiss possessed a satisfactory taste, flavor, and physicochemical and rheological properties. Further, it suggests that this cofermented cow-milk koumiss can be considered as a good alternative product for koumiss by breaking the resource limitation and reducing the material cost. Nevertheless, scale-up fermentation should also be carried out in future studies to support their industrial production. In addition, evaluation of the beneficial health effects of this fermented cow-milk koumiss is of significance to enhance its acceptability among different groups of people.
ACKNOWLEDGMENTS
No funding was received for the work. We thank Zaid Ashiq Khan (Northwest A & F University, Yangling, China) for his kind assistance concerning manuscript language improvement. We also acknowledge the suggestions and constructive comments from the editor and the anonymous reviewers. All authors declare that there is no conflict of interest. Laboratory experiments and methodology were performed by A. R., T. W., W. X., and A. R., and G.X. did formal analysis and software. A. R. wrote the original draft. Review and editing were performed by L. M., Y. H., Y. L., L. X., and Q. Z.. Supervision was performed by X. L..
Isolation and identification of lactic acid bacteria from koumiss in Eastern Inner Mongolia of China. In Proc. International Conference on Materials Science, Resource and Environmental Engineering, Xi'an City, Shaanxi, China. AIP Conference Proceedings 1794:050005.
Volatile organic compounds in naturally fermented milk and milk fermented using yeasts, lactic acid bacteria and their combinations as starter cultures.
In situ real-time monitoring of volatile metabolites of fermented milk by dynamic headspace sampling-atmospheric pressure ionization mass spectrometry.
in: Starr M.P. Stolp H. Truper H.G. The Prokaryotes: A Handbook on Habitats, Isolation, and Identification of Bacteria. vol 2. Springer-Verlag,
New York, NY1981: 1711-1742
Volatile organic compounds profile during milk fermentation by Lactobacillus pentosus and correlations between volatiles flavor and carbohydrate metabolism.
Use of lactobacilli strains with probiotic potential in traditional fermented milk and their impact on quality and safety related to Listeria monocytogenes.
Influence of addition of proteolytic strains of Lactobacillus delbrueckii ssp. bulgaricus to commercial ABT starter cultures on texture of yoghurt, exopolysaccharide production and survival of bacteria.