Advertisement
Research| Volume 104, ISSUE 1, P151-166, January 2021

Isolation and identification of microorganisms in Kazakhstan koumiss and their application in preparing cow-milk koumiss

Open ArchivePublished:November 05, 2020DOI:https://doi.org/10.3168/jds.2020-18527

      ABSTRACT

      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.

      Key words

      INTRODUCTION

      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 (
      • Hui Y.H.
      • Evranuz E.Ö.
      Handbook of animal-based fermented food and beverage technology.
      ). 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 (
      • Montanari G.
      • Zambonelli C.
      • Grazia L.
      • Kamesheva G.K.
      • Shigaeva M.K.
      Saccharomyces unisporus as the principal alcoholic fermentation microorganism of traditional koumiss.
      ;
      • Cagno R.D.
      • Tamborrino A.
      • Gallo G.
      • Leone C.
      • Angelis M.D.
      • Faccia M.
      • Amirante P.
      • Gobbetti M.
      Uses of mares' milk in manufacture of fermented milks.
      ), 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 (
      • Dönmez N.
      • Kısadere İ.
      • Balaban C.
      • Kadiralieva N.
      Effects of traditional homemade koumiss on some hematological and biochemical characteristics in sedentary men exposed to exercise.
      ). It has been found that protein, fat, lactose, minerals, enzymes, vitamins, and pigments are the main components of koumiss (
      • Malacarne M.
      • Martuzzi F.
      • Summer A.
      • Mariani P.
      Protein and fat composition of mare's milk: Some nutritional remarks with reference to human and cow's milk.
      ;
      • Danova S.
      • Petrov K.
      • Pavlov P.
      • Petrova P.
      Isolation and characterization of Lactobacillus strains involved in koumiss fermentation.
      ). Compared with other fermented milk, koumiss possesses almost all essential AA and high vitamin C levels that are important for humans (
      • Park Y.W.
      • Haenlein G.F.
      • Wendorff W.L.
      ;
      • Bornaz S.
      • Guizani N.
      • Sammari J.
      • Allouch W.
      • Sahli A.
      • Attia H.
      Physicochemical properties of fermented Arabian mares' milk.
      ). 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 (
      • Cagno R.D.
      • Tamborrino A.
      • Gallo G.
      • Leone C.
      • Angelis M.D.
      • Faccia M.
      • Amirante P.
      • Gobbetti M.
      Uses of mares' milk in manufacture of fermented milks.
      ;
      • Mu Z.
      • Yang X.
      • Yuan H.
      Detection and identification of wild yeast in Koumiss.
      ;
      • Dhewa T.
      • Mishra V.
      • Kumar N.
      • Sangu K.
      Koumiss: The nutritional and therapeutic values.
      ). 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 (
      • Langlois B.
      The history, ethnology and social importance of mare's milk consumption in central Asia.
      ;
      • Ishii S.
      • Hosino B.
      • Komiyama H.
      • Uehara A.
      • Nurtazin S.
      Study on production and properties of kumiss of herders in Mongolian dry steppe.
      ). 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 (
      • Park Y.W.
      • Zhang H.
      • Zhang B.
      • Zhang L.
      Mare milk.
      ;
      • Doreau M.
      • Martin-Rosset W.
      Animals that produce dairy foods Horse.
      ; Salimei and Park, 2017).
      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 (
      • Kücükcetin A.
      • Yaygin H.
      • Hinrichs J.
      • Kulozik U.
      Adaptation of bovine milk towards mares' milk composition by means of membrane technology for koumiss manufacture.
      ). 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 (
      • Norris J.
      • Berkeley R.
      • Logan N.
      • O'donnell A.
      The genera Bacillus and Sporolactobacillus.
      ) as well as biochemically characterized for catalase and oxidase testing before stock preparation (
      • Nair P.S.
      • Surendran P.K.
      Biochemical characterization of lactic acid bacteria isolated from fish and prawn.
      ). 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
      • Dertli E.
      • Mercan E.
      • Arıcı M.
      • Yılmaz M.T.
      • Sağdıç O.
      Characterisation of lactic acid bacteria from Turkish sourdough and determination of their exopolysaccharide (EPS) production characteristics.
      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 (
      • Saitou N.
      • Nei M.
      The neighbor-joining method: A new method for reconstructing phylogenetic trees.
      ;
      • Tamura K.
      • Peterson D.
      • Peterson N.
      • Stecher G.
      • Nei M.
      • Kumar S.
      MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.
      ).
      Table 1Primers used for 16S rRNA and 26S rRNA sequencing
      PrimerSequences (5′–3′)
      27F5′-AGAGTTTGATCMTGGCTCAG-3′
      1492R5′-TACGGYTACCTTGTTACGACTT-3′
      ITS15′-TCCGTAGGTGAACCTGCGG-3′
      ITS45′-TCCTCCGCTTATTGATATGC-3′
      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
      • Makimura K.
      • Murayama S.Y.
      • Yamaguchi H.
      Detection of a wide range of medically important fungi by the polymerase chain reaction.
      . 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 (
      • Saitou N.
      • Nei M.
      The neighbor-joining method: A new method for reconstructing phylogenetic trees.
      ).

      Selection of Best LAB and Yeasts Combinations

      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 (
      • Papadopoulou O.S.
      • Argyri A.A.
      • Varzakis E.
      • Sidira M.
      • Kourkoutas Y.
      • Galanis A.
      • Tassou C.
      • Chorianopoulos N.G.
      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 (
      • Wang Q.
      • Liu C.
      • Jing Y.
      • Fan S.
      • Cai J.
      Evaluation of fermentation conditions to improve the sensory quality of broomcorn millet sour porridge.
      ). 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:
      Y = β0 + β1A + β2B + β3C + β11A2 + β22B2 + β33C2 + β12AB + β13AC + β23BC,
      [1]


      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)
      StandardRunInoculum size (X1)Temperature (X2)Time (X3)Actual (sensory evaluation)Predicted value
      611.00 (7 g/L)0.00 (36°C)−1.00(12 h)6241.62
      32−1.00 (1 g/L)1.00 (40°C)0.00(16 h)6774.32
      1030.00 (4 g/L)1.00−1.007763.63
      1440.000.000.009080.38
      1250.001.001.00(20 h)5644.25
      56−1.000.00−1.003960.50
      970.00−1.00 (32°C)−1.004536.50
      1680.000.000.009769.75
      891.000.001.007543.13
      4101.001.000.007775.13
      17110.000.000.009161.88
      2121.00−1.000.007157.88
      713−1.000.001.003591.80
      11140.00−1.001.006091.80
      115−1.00−1.000.004591.80
      13160.000.000.009291.80
      15170.000.000.008991.80

      The Bradford Method to Determine Protein Content.

      The protein contents were determined using the Bradford assay (
      • Cheng Y.
      • Wei H.
      • Sun R.
      • Tian Z.
      • Zheng X.
      Rapid method for protein quantitation by Bradford assay after elimination of the interference of polysorbate 80.
      ). 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
      • Dan T.
      • Chen H.
      • Li T.
      • Tian J.
      • Ren W.
      • Zhang H.
      • Sun T.
      Influence of Lactobacillus plantarum P-8 on fermented milk flavor and storage stability.
      .

      Sensory Evaluation

      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 (
      • Wang Q.
      • Liu C.
      • Jing Y.
      • Fan S.
      • Cai J.
      Evaluation of fermentation conditions to improve the sensory quality of broomcorn millet sour porridge.
      ).

      Rheological Measurements

      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;
      • Wang X.
      • Chen Q.
      • Lü X.
      Pectin extracted from apple pomace and citrus peel by subcritical water.
      ).

      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
      • Macciola V.
      • Candela G.
      • De Leonardis A.
      Rapid gas-chromatographic method for the determination of diacetyl in milk, fermented milk and butter.
      .

      Solid Phase Microextraction GC–MS Analysis

      Volatile organic compounds (VOC) from Kazakhstan koumiss, cow-milk koumiss, and milk were analyzed using the headspace solid phase microextraction technique (
      • Dan T.
      • Wang D.
      • Jin R.
      • Zhang H.
      • Zhou T.
      • Sun T.
      Characterization of volatile compounds in fermented milk using solid-phase microextraction methods coupled with gas chromatography-mass spectrometry.
      ). 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)
      SourceSum of squaresdfMean squareF-valueP-valueSignificant
      Model6,157.489684.1630.300.0001
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X1-Inoculum1,225.1311,225.1354.260.0002
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X2-Temperature392.001392.0017.360.0042
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X3-Time1.1211.120.0500.8297
      X1 X264.00164.002.830.1361
      X1 X372.25172.253.200.1168
      X2 X3324.001324.0014.350.0068
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X1
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      1,184.8411,184.8452.480.0002
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X2
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      423.161423.1618.740.0034
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      X3
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      2,089.1612,089.1692.530.0001
      **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      Residual158.05722.58
      Lack of fit119.25339.754.100.1033
      Pure error38.8049.70
      Cor total
      Cor = corrected total sum of squares.
      6,315.5316
      R20.9750
      Adjusted R20.9428
      CV (%)6.92
      1 Cor = corrected total sum of squares.
      * **, and *** represent P < 0.05, P < 0.01, and P < 0.001, respectively.
      Table 4Comparison of physicochemical and microbiological properties of Kazakhstan koumiss and cow-milk koumiss
      Data are presented as the mean ± SD of 3 independent experiments. Statistical analyses were carried out by the Student's t-test.
      SampleAcidity (°T)pHEthanol (%)Protein (%)Lactic acid bacteria (cfu/mL)Yeast (cfu/mL)
      Kazakhstan koumiss97.40 ± 8.83
      Superscript letters in the same column indicate significant differences (P < 0.05).
      3.41 ± 0.08
      Superscript letters in the same column indicate significant differences (P < 0.05).
      2.12 ± 0.24
      Superscript letters in the same column indicate significant differences (P < 0.05).
      1.95 ± 0.13
      Superscript letters in the same column indicate significant differences (P < 0.05).
      2.24 ± (0.26) × 10
      Superscript letters in the same column indicate significant differences (P < 0.05).
      2.41 ± (0.40) × 10
      Superscript letters in the same column indicate significant differences (P < 0.05).
      Cow-milk koumiss71.67 ± 8.02
      Superscript letters in the same column indicate significant differences (P < 0.05).
      4.33 ± 0.01
      Superscript letters in the same column indicate significant differences (P < 0.05).
      0.50 ± 0.08
      Superscript letters in the same column indicate significant differences (P < 0.05).
      2.77 ± 0.24
      Superscript letters in the same column indicate significant differences (P < 0.05).
      9.90 ± (0.16) × 10
      Superscript letters in the same column indicate significant differences (P < 0.05).
      1.14 ± (0.19) × 10
      Superscript letters in the same column indicate significant differences (P < 0.05).
      a,b Superscript letters in the same column indicate significant differences (P < 0.05).
      1 Data are presented as the mean ± SD of 3 independent experiments. Statistical analyses were carried out by the Student's t-test.
      Table 5Volatile organic compounds (μg/L) quantified from Kazakhstan koumiss, cow-milk koumiss, and milk
      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/).
      CompoundKoumissCow-milk koumissMilk
      Alcohol
       Ethanol22.98 ± 0.14
      Different superscript letters in the same row indicate significant differences.
      20.83 ± 1.35
      Different superscript letters in the same row indicate significant differences.
      0.88 ± 0.22
      Different superscript letters in the same row indicate significant differences.
       1-Butanol, 3-methyl-4.70 ± 0.54
      Different superscript letters in the same row indicate significant differences.
      12.58 ± 1.31
      Different superscript letters in the same row indicate significant differences.
       ND
      ND = not detected.
       1-Butanol, 2-methyl-0.23 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      3.56 ± 0.40
      Different superscript letters in the same row indicate significant differences.
      ND
       6-Tridecanol, 3,9-diethyl-0.85 ± 0.02NDND
       2-Cyclopropyl-2-nitro-1-phenyl-ethanol0.60 ± 0.07NDND
       1-PropanolND0.16 ± 0.03ND
       1-Propanol, 2-methyl-ND1.52 ± 0.80ND
       Phenylethyl alcoholND0.23 ± 0.05ND
       n-Tridecan-1-olNDND0.04 ± 0.03
      Aldehyde
       Hexanal0.24 ± 0.18
      Different superscript letters in the same row indicate significant differences.
      0.21 ± 0.03
      Different superscript letters in the same row indicate significant differences.
      0.07 ± 0.23
      Different superscript letters in the same row indicate significant differences.
       Heptanal0.77 ± 0.22
      Different superscript letters in the same row indicate significant differences.
      0.68 ± 0.07
      Different superscript letters in the same row indicate significant differences.
      0.13 ± 0.32
      Different superscript letters in the same row indicate significant differences.
       Octanal0.84 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.77 ± 0.03
      Different superscript letters in the same row indicate significant differences.
      0.01 ± 0.20
      Different superscript letters in the same row indicate significant differences.
       Nonanal0.53 ± 0.04
      Different superscript letters in the same row indicate significant differences.
      0.50 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.05 ± 0.07
      Different superscript letters in the same row indicate significant differences.
       2,5-Dihydroxybenzaldehyde, 2TMS derivative0.18 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.15 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.20 ± 0.03
      Different superscript letters in the same row indicate significant differences.
       Decanal0.30 ± 0.03
      Different superscript letters in the same row indicate significant differences.
      0.22 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      ND
       Benzaldehyde, 2,4-dimethyl-0.22 ± 0.006NDND
       Undecanal0.13 ± 0.005NDND
       Acetaldehyde0.24 ± 0.17
      Different superscript letters in the same row indicate significant differences.
      1.16 ± 0.29
      Different superscript letters in the same row indicate significant differences.
      ND
       Propanal, 2-methyl-ND0.13 ± 0.02ND
      Ketone
       2-NonadecanoneND0.77 ± 0.04ND
       2-Heptanone0.24 ± 0.058
      Different superscript letters in the same row indicate significant differences.
      0.17 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.33 ± 0.10
      Different superscript letters in the same row indicate significant differences.
       9-Octadecanone0.11 ± 0.08
      Different superscript letters in the same row indicate significant differences.
      0.69 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.59 ± 0.12
      Different superscript letters in the same row indicate significant differences.
       2-OctanoneND0.2 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.32 ± 0.07
      Different superscript letters in the same row indicate significant differences.
       2-NonanoneND0.16 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.20 ± 0.03
      Different superscript letters in the same row indicate significant differences.
       4-HeptadecanoneND0.44 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.40 ± 0.07
      Different superscript letters in the same row indicate significant differences.
      Ester
       Butanoic acid, ethyl ester0.31 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.16 ± 0.002
      Different superscript letters in the same row indicate significant differences.
      ND
       Propanoic acid, 2-hydroxy-, ethyl ester4.67 ± 0.17NDND
       Ethyl acetate0.69 ± 0.43
      Different superscript letters in the same row indicate significant differences.
      9.11 ± 0.74
      Different superscript letters in the same row indicate significant differences.
      0.23 ± 0.04
      Different superscript letters in the same row indicate significant differences.
       Octanoic acid, ethyl ester1.53 ± 0.04NDND
       Decanoic acid, ethyl ester0.24 ± 0.05NDND
       Diethyl phthalate1.22 ± 0.08NDND
      Organic acid
       L-Lactic acid0.85 ± 0.05NDND
       Octanoic acid4.36 ± 0.39NDND
       n-Decanoic acid0.23 ± 0.01NDND
       Acetic acid2.87 ± 0.05
      Different superscript letters in the same row indicate significant differences.
      0.33
      Different superscript letters in the same row indicate significant differences.
      ± 0.01
      ND
       Butanoic acid, 3-methyl-0.32 ± 0.00NDND
      Imine
       Oxime-, methoxy-phenyl-0.51 ± 0.25
      Different superscript letters in the same row indicate significant differences.
      0.09 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.14 ± 0.03
      Different superscript letters in the same row indicate significant differences.
      Phenol
       Phenol, 3,5-bis(1,1-dimethylethyl)-0.13 ± 0.02
      Different superscript letters in the same row indicate significant differences.
      0.10 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      0.14 ± 0.01
      Different superscript letters in the same row indicate significant differences.
      Amine
       Benzeneethanamine0.17 ± 0.01NDND
       l-Alanine ethylamide, (S)-0.57 ± 0.08
      Different superscript letters in the same row indicate significant differences.
      0.47 ± 0.03
      Different superscript letters in the same row indicate significant differences.
      0.58 ± 0.14
      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/).
      2 ND = not detected.

      RESULTS AND DISCUSSION

      Screening of LAB and Yeast Strains from Kazakhstan Koumiss

      The previous literature revealed that LAB and yeast strains are 2 major microorganisms in traditional fermented foods (
      • Tamang J.P.
      • Fleet G.H.
      Yeasts diversity in fermented foods and beverages.
      ;
      • Lv X.-C.
      • Weng X.
      • Zhang W.
      • Rao P.-F.
      • Ni L.
      Microbial diversity of traditional fermentation starters for Hong Qu glutinous rice wine as determined by PCR-mediated DGGE.
      ;
      • Tamang J.P.
      • Watanabe K.
      • Holzapfel W.H.
      Diversity of microorganisms in global fermented foods and beverages.
      ;
      • García C.
      • Rendueles M.
      • Díaz M.
      Liquid-phase food fermentations with microbial consortia involving lactic acid bacteria.
      ). 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 (
      • Forouhandeh H.
      • Vahed A.
      • Hejazi M.
      • Nahaei M.
      • Dibavar M.A.
      Isolation and phenotypic characterization of Lactobacillus species from various dairy products.
      ). 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
      • Hao Y.
      • Zhao L.
      • Zhang H.
      • Zhai Z.
      • Huang Y.
      • Liu X.
      • Zhang L.
      Identification of the bacterial biodiversity in koumiss by denaturing gradient gel electrophoresis and species-specific polymerase chain reaction.
      . Moreover, previous studies have observed the presence of these strains in traditional dairy products, including koumiss (
      • Bilige M.
      • Liu W.
      • Rina W.
      • Wang L.
      • Sun T.
      • Wang J.
      • Li H.
      • Zhang H.
      Evaluation of potential probiotics properties of the screened Lactobacilli isolated from home-made koumiss in Mongolia.
      ). The microbial communities of koumiss are similar to those in kefir, cheese, and yogurt (
      • Buttriss J.
      Nutritional properties of fermented milk products.
      ;
      • Pan D.D.
      • Zeng X.Q.
      • Yan Y.T.
      Characterisation of Lactobacillus fermentum SM-7 isolated from koumiss, a potential probiotic bacterium with cholesterol-lowering effects.
      ).
      Figure thumbnail gr1
      Figure 1Phylogenetic tree based on 16S rRNA sequences constructed by the maximum likelihood method.
      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 (
      • Wang S.-Y.
      • Chen H.-C.
      • Liu J.-R.
      • Lin Y.-C.
      • Chen M.-J.
      Identification of yeasts and evaluation of their distribution in Taiwanese kefir and viili starters.
      ). 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 (
      • Lore T.A.
      • Mbugua S.K.
      • Wangoh J.
      Enumeration and identification of microflora in suusac, a Kenyan traditional fermented camel milk product.
      ;
      • Watanabe K.
      • Fujimoto J.
      • Sasamoto M.
      • Dugersuren J.
      • Tumursuh T.
      • Demberel S.
      Diversity of lactic acid bacteria and yeasts in Airag and Tarag, traditional fermented milk products of Mongolia.
      ). Another study shows that yeast plays an important role in defining the qualitative traits of koumiss (
      • Mu Z.
      • Yang X.
      • Yuan H.
      Detection and identification of wild yeast in Koumiss.
      ). 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 (
      • Han X.
      • Yang Z.
      • Jing X.
      • Yu P.
      • Zhang Y.
      • Yi H.
      • Zhang L.
      Improvement of the texture of yogurt by use of exopolysaccharide producing lactic acid bacteria.
      ). 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 (
      • Kavitake D.
      • Kandasamy S.
      • Devi P.B.
      • Shetty P.H.
      Recent developments on encapsulation of lactic acid bacteria as potential starter culture in fermented foods–A review.
      ;
      • Rakhmanova A.
      • Khan Z.
      • Shah K.
      A mini review fermentation and preservation: Role of Lactic Acid Bacteria.
      ). 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 (
      • Liu S.-Q.
      • Holland R.
      • Crow V.
      Synthesis of ethyl butanoate by a commercial lipase in aqueous media under conditions relevant to cheese ripening.
      ;
      • Wang C.
      • Shi D.
      • Gong G.
      Microorganisms in Daqu: A starter culture of Chinese Maotai-flavor liquor.
      ).
      Figure thumbnail gr2
      Figure 2Phylogenetic tree based on 26S rRNA sequences constructed by the maximum likelihood method.

      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 (
      • Yu J.
      • Wang W.
      • Menghe B.
      • Jiri M.
      • Wang H.
      • Liu W.
      • Bao Q.
      • Lu Q.
      • Zhang J.
      • Wang F.
      • Xu H.Y.
      • Sun T.S.
      • Zhang H.P.
      Diversity of lactic acid bacteria associated with traditional fermented dairy products in Mongolia.
      ). 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 (
      • Shihata A.
      • Shah N.
      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 (
      • Gethins L.
      • Rea M.C.
      • Stanton C.
      • Ross R.P.
      • Kilcawley K.
      • O'Sullivan M.
      • Crotty S.
      • Morrissey J.P.
      Acquisition of the yeast Kluyveromyces marxianus from unpasteurised milk by a kefir grain enhances kefir quality.
      ). Some studies found that the cooperation of LAB and yeast is crucial for the production of acidified milk during fungal-lactic fermentations (
      • Tamang J.P.
      • Watanabe K.
      • Holzapfel W.H.
      Diversity of microorganisms in global fermented foods and beverages.
      ). 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 (
      • Lane M.M.
      • Morrissey J.P.
      Kluyveromyces marxianus: A yeast emerging from its sister's shadow.
      ;
      • Morrissey J.P.
      • Etschmann M.M.
      • Schrader J.
      • de Billerbeck G.M.
      Cell factory applications of the yeast Kluyveromyces marxianus for the biotechnological production of natural flavour and fragrance molecules.
      ). 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 (
      • Helland M.H.
      • Wicklund T.
      • Narvhus J.A.
      Growth and metabolism of selected strains of probiotic bacteria, in maize porridge with added malted barley.
      ). 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 (
      • Cui X.-H.
      • Chen S.-J.
      • Wang Y.
      • Han J.-R.
      Fermentation conditions of walnut milk beverage inoculated with kefir grains.
      ). 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 (
      • Wang Q.
      • Liu C.
      • Jing Y.
      • Fan S.
      • Cai J.
      Evaluation of fermentation conditions to improve the sensory quality of broomcorn millet sour porridge.
      ). Moreover, it has been observed that the fermentation time significantly affects the rate of fermentation (
      • Pereira A.L.F.
      • Maciel T.C.
      • Rodrigues S.
      Probiotic beverage from cashew apple juice fermented with Lactobacillus casei.
      ;
      • Cui X.-H.
      • Chen S.-J.
      • Wang Y.
      • Han J.-R.
      Fermentation conditions of walnut milk beverage inoculated with kefir grains.
      ).
      Figure thumbnail gr3
      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:
      Y=91.80+12.38X1+7X2+0.37X3-4X1X2+4.25X1X3-9X2X3-16.78X12-10.03X22-22.27X32,
      [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, X12, and X32, with P-values less than 0.001 each; X1, X2, X3, and X22 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 thumbnail gr4
      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 (
      • Kücükcetin A.
      • Yaygin H.
      • Hinrichs J.
      • Kulozik U.
      Adaptation of bovine milk towards mares' milk composition by means of membrane technology for koumiss manufacture.
      ). 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
      • Chen Y.
      • Wang Z.
      • Chen X.
      • Liu Y.
      • Zhang H.
      • Sun T.
      Identification of angiotensin I-converting enzyme inhibitory peptides from koumiss, a traditional fermented mare's milk.
      and
      • Ishii S.
      • Hosino B.
      • Komiyama H.
      • Uehara A.
      • Nurtazin S.
      Study on production and properties of kumiss of herders in Mongolian dry steppe.
      . 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
      • Ha S.
      • Leng A.
      • Mang L.
      Koumiss and its medicinal values.
      . 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 (
      • Ha S.
      • Leng A.
      • Mang L.
      Koumiss and its medicinal values.
      ). 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 (
      • Bai L.
      • Ji S.
      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 (
      • Gadaga T.H.
      • Viljoen B.C.
      • Narvhus J.A.
      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 (
      • Pan D.D.
      • Wu Z.
      • Peng T.
      • Zeng X.
      • Li H.
      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 (
      • Smit B.A.
      • Engels W.J.
      • Smit G.
      Branched chain aldehydes: Production and breakdown pathways and relevance for flavour in foods.
      ). In dairy foods, the quantity of 3-methylbutanol is considerably high and imparts an “alcoholic and floral” flavor to fermented milk (
      • Molimard P.
      • Spinnler H.
      Compounds involved in the flavor of surface mold-ripened cheeses: Origins and properties.
      ). 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 (
      • Dursun A.
      • Güler Z.
      • Şekerli Y.E.
      Characterization of volatile compounds and organic acids in ultra-high-temperature milk packaged in tetra brik cartons.
      ). 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 (
      • Beshkova D.
      • Simova E.
      • Frengova G.
      • Simov Z.
      • Dimitrov Z.P.
      Production of volatile aroma compounds by kefir starter cultures.
      ). 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 (
      • Simova E.
      • Beshkova D.
      • Angelov A.
      • Hristozova T.
      • Frengova G.
      • Spasov Z.
      Lactic acid bacteria and yeasts in kefir grains and kefir made from them.
      ). 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 (
      • Smit B.A.
      • Engels W.J.
      • Smit G.
      Branched chain aldehydes: Production and breakdown pathways and relevance for flavour in foods.
      ). 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 (
      • Pan D.D.
      • Wu Z.
      • Peng T.
      • Zeng X.
      • Li H.
      Volatile organic compounds profile during milk fermentation by Lactobacillus pentosus and correlations between volatiles flavor and carbohydrate metabolism.
      ). It was also hypothesized that volatile compounds may affect the sensory characteristics including flavor and taste in fermented milk products (
      • Cheng H.
      Volatile flavor compounds in yogurt: A review.
      ). 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 (
      • Liang H.
      • Yang S.
      • Liu W.
      • Yuan L.
      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 (
      • Campanella O.
      • Dorward N.
      • Singh H.
      A study of the rheological properties of concentrated food emulsions.
      ;
      • Abu-Jdayil B.
      • Mohameed H.
      Experimental and modelling studies of the flow properties of concentrated yogurt as affected by the storage time.
      ). 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 (
      • Ochirkhuyag B.
      • Chobert J.-M.
      • Dalgalarrondo M.
      • Haertlé T.
      Characterization of mare caseins. Identification of αS1- and αS2-caseins.
      ;
      • Malacarne M.
      • Martuzzi F.
      • Summer A.
      • Mariani P.
      Protein and fat composition of mare's milk: Some nutritional remarks with reference to human and cow's milk.
      ). It was further investigated that the concentration of κ-casein is higher in mare milk as compared with cow milk (
      • Egito A.
      • Girardet J.-M.
      • Miclo L.
      • Mollé D.
      • Humbert G.
      • Gaillard J.-L.
      Susceptibility of equine κ-and β-caseins to hydrolysis by chymosin.
      ). The casein micelles of mare milk are larger and less porous than cow milk (
      • Buchheim W.
      • Lund S.
      • Scholtissek J.
      Vergleichende Untersuchungen zur Struktur und Grösse von Caseinmicellen in der Milch verschiedener Species.
      ). 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 thumbnail gr5
      Figure 5Frequency sweeps of Kazakhstan koumiss, cow-milk koumiss, and milk at 1% strain. G′ = storage modulus; G″ = loss modulus.

      CONCLUSIONS

      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..

      Supplementary Material

      REFERENCES

        • Abu-Jdayil B.
        • Mohameed H.
        Experimental and modelling studies of the flow properties of concentrated yogurt as affected by the storage time.
        J. Food Eng. 2002; 52: 359-365
        • Bai L.
        • Ji S.
        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.
        • Beshkova D.
        • Simova E.
        • Frengova G.
        • Simov Z.
        • Dimitrov Z.P.
        Production of volatile aroma compounds by kefir starter cultures.
        Int. Dairy J. 2003; 13: 529-535
        • Bilige M.
        • Liu W.
        • Rina W.
        • Wang L.
        • Sun T.
        • Wang J.
        • Li H.
        • Zhang H.
        Evaluation of potential probiotics properties of the screened Lactobacilli isolated from home-made koumiss in Mongolia.
        Ann. Microbiol. 2009; 59: 493-498
        • Bornaz S.
        • Guizani N.
        • Sammari J.
        • Allouch W.
        • Sahli A.
        • Attia H.
        Physicochemical properties of fermented Arabian mares' milk.
        Int. Dairy J. 2010; 20: 500-505
        • Buchheim W.
        • Lund S.
        • Scholtissek J.
        Vergleichende Untersuchungen zur Struktur und Grösse von Caseinmicellen in der Milch verschiedener Species.
        Kieler Milchwirtschaftliche Forschungsberichte. 1989; 41: 253-265
        • Buttriss J.
        Nutritional properties of fermented milk products.
        Int. J. Dairy Technol. 1997; 50: 21-27
        • Cagno R.D.
        • Tamborrino A.
        • Gallo G.
        • Leone C.
        • Angelis M.D.
        • Faccia M.
        • Amirante P.
        • Gobbetti M.
        Uses of mares' milk in manufacture of fermented milks.
        Int. Dairy J. 2004; 14: 767-775
        • Campanella O.
        • Dorward N.
        • Singh H.
        A study of the rheological properties of concentrated food emulsions.
        J. Food Eng. 1995; 25: 427-440
        • Chen Y.
        • Wang Z.
        • Chen X.
        • Liu Y.
        • Zhang H.
        • Sun T.
        Identification of angiotensin I-converting enzyme inhibitory peptides from koumiss, a traditional fermented mare's milk.
        J. Dairy Sci. 2010; 93 (20172208): 884-892
        • Cheng H.
        Volatile flavor compounds in yogurt: A review.
        Crit. Rev. Food Sci. Nutr. 2010; 50 (21108074): 938-950
        • Cheng Y.
        • Wei H.
        • Sun R.
        • Tian Z.
        • Zheng X.
        Rapid method for protein quantitation by Bradford assay after elimination of the interference of polysorbate 80.
        Anal. Biochem. 2016; 494 (26545323): 37-39
        • Cui X.-H.
        • Chen S.-J.
        • Wang Y.
        • Han J.-R.
        Fermentation conditions of walnut milk beverage inoculated with kefir grains.
        Lebensm. Wiss. Technol. 2013; 50: 349-352
        • Dan T.
        • Chen H.
        • Li T.
        • Tian J.
        • Ren W.
        • Zhang H.
        • Sun T.
        Influence of Lactobacillus plantarum P-8 on fermented milk flavor and storage stability.
        Front. Microbiol. 2019; 9 (30687239): 3133-3146
        • Dan T.
        • Wang D.
        • Jin R.
        • Zhang H.
        • Zhou T.
        • Sun T.
        Characterization of volatile compounds in fermented milk using solid-phase microextraction methods coupled with gas chromatography-mass spectrometry.
        J. Dairy Sci. 2017; 100 (28161169): 2488-2500
        • Danova S.
        • Petrov K.
        • Pavlov P.
        • Petrova P.
        Isolation and characterization of Lactobacillus strains involved in koumiss fermentation.
        Int. J. Dairy Technol. 2005; 58: 100-105
        • Dertli E.
        • Mercan E.
        • Arıcı M.
        • Yılmaz M.T.
        • Sağdıç O.
        Characterisation of lactic acid bacteria from Turkish sourdough and determination of their exopolysaccharide (EPS) production characteristics.
        Lebensm. Wiss. Technol. 2016; 71: 116-124
        • Dhewa T.
        • Mishra V.
        • Kumar N.
        • Sangu K.
        Koumiss: The nutritional and therapeutic values.
        in: Fermented Fermented Milk and Dairy Products. CRC Press, Boca Raton, FL2015: 483-494
        • Dönmez N.
        • Kısadere İ.
        • Balaban C.
        • Kadiralieva N.
        Effects of traditional homemade koumiss on some hematological and biochemical characteristics in sedentary men exposed to exercise.
        Biotech. Histochem. 2014; 89 (24893332): 558-563
        • Doreau M.
        • Martin-Rosset W.
        Animals that produce dairy foods Horse.
        in: Encyclopedia of Dairy Sciences. second ed. Academic Press, Cambridge, MA2011: 358-364
        • Dursun A.
        • Güler Z.
        • Şekerli Y.E.
        Characterization of volatile compounds and organic acids in ultra-high-temperature milk packaged in tetra brik cartons.
        Int. J. Food Prop. 2017; 20: 1511-1521
        • Egito A.
        • Girardet J.-M.
        • Miclo L.
        • Mollé D.
        • Humbert G.
        • Gaillard J.-L.
        Susceptibility of equine κ-and β-caseins to hydrolysis by chymosin.
        Int. Dairy J. 2001; 11: 885-893
        • Forouhandeh H.
        • Vahed A.
        • Hejazi M.
        • Nahaei M.
        • Dibavar M.A.
        Isolation and phenotypic characterization of Lactobacillus species from various dairy products.
        Curr. Res. Bacteriol. 2010; 3: 84-88
        • Gadaga T.H.
        • Viljoen B.C.
        • Narvhus J.A.
        Volatile organic compounds in naturally fermented milk and milk fermented using yeasts, lactic acid bacteria and their combinations as starter cultures.
        Food Technol. Biotechnol. 2007; 45: 195-200
        • García C.
        • Rendueles M.
        • Díaz M.
        Liquid-phase food fermentations with microbial consortia involving lactic acid bacteria.
        Food Res. Int. 2019; 119 (30884650): 207-220
        • Gethins L.
        • Rea M.C.
        • Stanton C.
        • Ross R.P.
        • Kilcawley K.
        • O'Sullivan M.
        • Crotty S.
        • Morrissey J.P.
        Acquisition of the yeast Kluyveromyces marxianus from unpasteurised milk by a kefir grain enhances kefir quality.
        FEMS Microbiol. Lett. 2016; 363 (27369085)fnw165
        • Ha S.
        • Leng A.
        • Mang L.
        Koumiss and its medicinal values.
        China Journal of Chinese Materia. Medica. 2003; 28: 11-14
        • Han X.
        • Yang Z.
        • Jing X.
        • Yu P.
        • Zhang Y.
        • Yi H.
        • Zhang L.
        Improvement of the texture of yogurt by use of exopolysaccharide producing lactic acid bacteria.
        BioMed Res. Int. 2016; 2016 (27294135)7945675
        • Hao Y.
        • Zhao L.
        • Zhang H.
        • Zhai Z.
        • Huang Y.
        • Liu X.
        • Zhang L.
        Identification of the bacterial biodiversity in koumiss by denaturing gradient gel electrophoresis and species-specific polymerase chain reaction.
        J. Dairy Sci. 2010; 93 (20412906): 1926-1933
        • Helland M.H.
        • Wicklund T.
        • Narvhus J.A.
        Growth and metabolism of selected strains of probiotic bacteria, in maize porridge with added malted barley.
        Int. J. Food Microbiol. 2004; 91 (14984778): 305-313
        • Hui Y.H.
        • Evranuz E.Ö.
        Handbook of animal-based fermented food and beverage technology.
        Vol. 1. CRC press, Boca Raton, FL2012
        • Ishii S.
        • Hosino B.
        • Komiyama H.
        • Uehara A.
        • Nurtazin S.
        Study on production and properties of kumiss of herders in Mongolian dry steppe.
        J. of Arid Land Stud. 2014; 24: 196-197
        • Kavitake D.
        • Kandasamy S.
        • Devi P.B.
        • Shetty P.H.
        Recent developments on encapsulation of lactic acid bacteria as potential starter culture in fermented foods–A review.
        Food Biosci. 2018; 21: 34-44
        • Kücükcetin A.
        • Yaygin H.
        • Hinrichs J.
        • Kulozik U.
        Adaptation of bovine milk towards mares' milk composition by means of membrane technology for koumiss manufacture.
        Int. Dairy J. 2003; 13: 945-951
        • Lane M.M.
        • Morrissey J.P.
        Kluyveromyces marxianus: A yeast emerging from its sister's shadow.
        Fungal Biol. Rev. 2010; 24: 17-26
        • Langlois B.
        The history, ethnology and social importance of mare's milk consumption in central Asia.
        J. Life Sci. 2011; 5: 863-872
        • Liang H.
        • Yang S.
        • Liu W.
        • Yuan L.
        In situ real-time monitoring of volatile metabolites of fermented milk by dynamic headspace sampling-atmospheric pressure ionization mass spectrometry.
        Shipin Kexue. 2012; 33: 307-310
        • Liu S.-Q.
        • Holland R.
        • Crow V.
        Synthesis of ethyl butanoate by a commercial lipase in aqueous media under conditions relevant to cheese ripening.
        J. Dairy Res. 2003; 70 (12916832): 359-363
        • Lore T.A.
        • Mbugua S.K.
        • Wangoh J.
        Enumeration and identification of microflora in suusac, a Kenyan traditional fermented camel milk product.
        Lebensm. Wiss. Technol. 2005; 38: 125-130
        • Lv X.-C.
        • Weng X.
        • Zhang W.
        • Rao P.-F.
        • Ni L.
        Microbial diversity of traditional fermentation starters for Hong Qu glutinous rice wine as determined by PCR-mediated DGGE.
        Food Control. 2012; 28: 426-434
        • Macciola V.
        • Candela G.
        • De Leonardis A.
        Rapid gas-chromatographic method for the determination of diacetyl in milk, fermented milk and butter.
        Food Control. 2008; 19: 873-878
        • Makimura K.
        • Murayama S.Y.
        • Yamaguchi H.
        Detection of a wide range of medically important fungi by the polymerase chain reaction.
        J. Med. Microbiol. 1994; 40 (8176723): 358-364
        • Malacarne M.
        • Martuzzi F.
        • Summer A.
        • Mariani P.
        Protein and fat composition of mare's milk: Some nutritional remarks with reference to human and cow's milk.
        Int. Dairy J. 2002; 12: 869-877
        • Molimard P.
        • Spinnler H.
        Compounds involved in the flavor of surface mold-ripened cheeses: Origins and properties.
        J. Dairy Sci. 1996; 79: 169-184
        • Montanari G.
        • Zambonelli C.
        • Grazia L.
        • Kamesheva G.K.
        • Shigaeva M.K.
        Saccharomyces unisporus as the principal alcoholic fermentation microorganism of traditional koumiss.
        J. Dairy Res. 1996; 63: 327-331
        • Morrissey J.P.
        • Etschmann M.M.
        • Schrader J.
        • de Billerbeck G.M.
        Cell factory applications of the yeast Kluyveromyces marxianus for the biotechnological production of natural flavour and fragrance molecules.
        Yeast. 2015; 32 (25393382): 3-16
        • Mu Z.
        • Yang X.
        • Yuan H.
        Detection and identification of wild yeast in Koumiss.
        Food Microbiol. 2012; 31 (22608237): 301-308
        • Nair P.S.
        • Surendran P.K.
        Biochemical characterization of lactic acid bacteria isolated from fish and prawn.
        J. Culture. Collection. 2005; 4: 48-52
        • Norris J.
        • Berkeley R.
        • Logan N.
        • O'donnell A.
        The genera Bacillus and Sporolactobacillus.
        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
        • Ochirkhuyag B.
        • Chobert J.-M.
        • Dalgalarrondo M.
        • Haertlé T.
        Characterization of mare caseins. Identification of αS1- and αS2-caseins.
        Lait. 2000; 80: 223-235
        • Pan D.D.
        • Wu Z.
        • Peng T.
        • Zeng X.
        • Li H.
        Volatile organic compounds profile during milk fermentation by Lactobacillus pentosus and correlations between volatiles flavor and carbohydrate metabolism.
        J. Dairy Sci. 2014; 97 (24359834): 624-631
        • Pan D.D.
        • Zeng X.Q.
        • Yan Y.T.
        Characterisation of Lactobacillus fermentum SM-7 isolated from koumiss, a potential probiotic bacterium with cholesterol-lowering effects.
        J. Sci. Food Agric. 2011; 91 (21218486): 512-518
        • Papadopoulou O.S.
        • Argyri A.A.
        • Varzakis E.
        • Sidira M.
        • Kourkoutas Y.
        • Galanis A.
        • Tassou C.
        • Chorianopoulos N.G.
        Use of lactobacilli strains with probiotic potential in traditional fermented milk and their impact on quality and safety related to Listeria monocytogenes.
        Int. Dairy J. 2019; 98: 44-53
        • Park Y.W.
        • Haenlein G.F.
        • Wendorff W.L.
        Park W.Y. Haenlein G.F.W. Handbook of Milk of Non-Bovine Mammals. Blackwell Publishing, Oxford, UK2006: 512
        • Park Y.W.
        • Zhang H.
        • Zhang B.
        • Zhang L.
        Mare milk.
        in: Park W.Y. Haenlein G.F.W. Handbook of Milk of Nonbovine Mammals. Blackwell Publishing, Oxford, UK2006: 275-296
        • Pereira A.L.F.
        • Maciel T.C.
        • Rodrigues S.
        Probiotic beverage from cashew apple juice fermented with Lactobacillus casei.
        Food Res. Int. 2011; 44: 1276-1283
        • Rakhmanova A.
        • Khan Z.
        • Shah K.
        A mini review fermentation and preservation: Role of Lactic Acid Bacteria.
        Food Process Technol. 2018; 6: 414-417
        • Saitou N.
        • Nei M.
        The neighbor-joining method: A new method for reconstructing phylogenetic trees.
        Mol. Biol. Evol. 1987; 4 (3447015): 406-425
        • Shihata A.
        • Shah N.
        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.
        Int. Dairy J. 2002; 12: 765-772
        • Simova E.
        • Beshkova D.
        • Angelov A.
        • Hristozova T.
        • Frengova G.
        • Spasov Z.
        Lactic acid bacteria and yeasts in kefir grains and kefir made from them.
        J. Ind. Microbiol. Biotechnol. 2002; 28 (11938463): 1-6
        • Smit B.A.
        • Engels W.J.
        • Smit G.
        Branched chain aldehydes: Production and breakdown pathways and relevance for flavour in foods.
        Appl. Microbiol. Biotechnol. 2009; 81 (19015847): 987-999
        • Tamang J.P.
        • Fleet G.H.
        Yeasts diversity in fermented foods and beverages.
        in: Yeast Biotechnology: Diversity and Applications. Springer, New York, NY2009: 169-198
        • Tamang J.P.
        • Watanabe K.
        • Holzapfel W.H.
        Diversity of microorganisms in global fermented foods and beverages.
        Front. Microbiol. 2016; 7 (27047484): 377
        • Tamura K.
        • Peterson D.
        • Peterson N.
        • Stecher G.
        • Nei M.
        • Kumar S.
        MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.
        Mol. Biol. Evol. 2011; 28 (21546353): 2731-2739
        • Wang C.
        • Shi D.
        • Gong G.
        Microorganisms in Daqu: A starter culture of Chinese Maotai-flavor liquor.
        World J. Microbiol. Biotechnol. 2008; 24: 2183-2190
        • Wang Q.
        • Liu C.
        • Jing Y.
        • Fan S.
        • Cai J.
        Evaluation of fermentation conditions to improve the sensory quality of broomcorn millet sour porridge.
        Lebensm. Wiss. Technol. 2019; 104: 165-172
        • Wang S.-Y.
        • Chen H.-C.
        • Liu J.-R.
        • Lin Y.-C.
        • Chen M.-J.
        Identification of yeasts and evaluation of their distribution in Taiwanese kefir and viili starters.
        J. Dairy Sci. 2008; 91 (18832201): 3798-3805
        • Wang X.
        • Chen Q.
        • Lü X.
        Pectin extracted from apple pomace and citrus peel by subcritical water.
        Food Hydrocoll. 2014; 38: 129-137
        • Watanabe K.
        • Fujimoto J.
        • Sasamoto M.
        • Dugersuren J.
        • Tumursuh T.
        • Demberel S.
        Diversity of lactic acid bacteria and yeasts in Airag and Tarag, traditional fermented milk products of Mongolia.
        World J. Microbiol. Biotechnol. 2008; 24: 1313-1325
        • Yu J.
        • Wang W.
        • Menghe B.
        • Jiri M.
        • Wang H.
        • Liu W.
        • Bao Q.
        • Lu Q.
        • Zhang J.
        • Wang F.
        • Xu H.Y.
        • Sun T.S.
        • Zhang H.P.
        Diversity of lactic acid bacteria associated with traditional fermented dairy products in Mongolia.
        J. Dairy Sci. 2011; 94 (21700007): 3229-3241