Next Article in Journal
Development of Simultaneous Analytical Method for Imidazolinone Herbicides from Livestock Products by UHPLC-MSMS
Next Article in Special Issue
Production of Indigo by Recombinant Escherichia coli with Expression of Monooxygenase, Tryptophanase, and Molecular Chaperone
Previous Article in Journal
Extrusion Simulation for the Design of Cereal and Legume Foods
Previous Article in Special Issue
Characterization of Key Odor-Active Compounds in Sun-Dried Black Tea by Sensory and Instrumental-Directed Flavor Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr

by
Thanyaporn Srimahaeak
1,
Mikael Agerlin Petersen
1,
Søren K. Lillevang
2,
Lene Jespersen
1 and
Nadja Larsen
1,*
1
Department of Food Science, University of Copenhagen, Rolighedsvej 26, 1958 Frederiksberg, Denmark
2
Arla Foods Innovation Centre, Agro Food Park 19, 8200 Aarhus, Denmark
*
Author to whom correspondence should be addressed.
Foods 2022, 11(12), 1776; https://doi.org/10.3390/foods11121776
Submission received: 18 May 2022 / Revised: 9 June 2022 / Accepted: 14 June 2022 / Published: 16 June 2022
(This article belongs to the Special Issue Microbiological Safety and Quality of Fermented Products)

Abstract

:
This study investigated the spoilage potential of yeast strains Kluyveromyces marxianus (Km1, Km2 and Km3), Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1 grown in skyr in cold storage. Yeast strains were isolated from skyr and identified by sequencing of the 26S rRNA gene. K. marxianus yeasts were grown in skyr to high numbers, generating large amounts of volatile organic compounds (VOC) associated with off-flavours, among them were alcohols (3-methyl-1-butanol, 2-methyl-1-propanol and 1-hexanol), esters (ethyl acetate and 3-methylbutyl acetate) and aldehydes (hexanal, methylbutanal and methylpropanal). Growth of P. kudriavzevii Pk1 led to moderate increases in several alcohols and esters (mostly, 3-methyl-1-butanol and ethyl acetate), whereas only minor shifts in VOCs were associated with T. delbrueckii Td2. The levels of the key aroma compounds, diacetyl and acetoin, were significantly decreased by all K. marxianus strains and P. kudriavzevii Pk1. In contrast to the other yeast species, K. marxianus was able to utilize lactose, producing ethanol and carbon dioxide. Based on the overall results, K. marxianus was characterised by the highest spoilage potential. The study revealed the differences between the yeast species in fermentative and spoilage activities, and clarified the role of yeast metabolites for off-flavour formation and quality defects in skyr during cold storage.

1. Introduction

Skyr is a traditional Icelandic dairy product, which nowadays has become popular in many countries of Europe and North America. In the UK and Denmark, the current share of skyr in the total yoghurt market is as high as 24% [1]. Skyr is characterised by its high protein content (10–15%), creamy smooth texture and typical sour dairy flavour [2]. Though skyr is classified in the market as a yoghurt-type product, it is technically a fresh cheese, traditionally made with addition of rennet, an enzyme used for milk coagulation and curd formation [2,3]. Today, most of the modernized dairies solely use starter cultures instead of rennet in skyr production. In contrast to yoghurt, which is typically made from full fat milk, skyr is produced from skimmed milk, and whey is removed after milk fermentation by either centrifugation or ultrafiltration [3]. The major starter cultures in skyr and yoghurt comprise strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. bulgaricus. The role of starter cultures for yoghurt quality has been widely studied, while for skyr this knowledge is insufficient. In yoghurt, a desired aroma depends on the balanced amounts of organic acids and the key aroma compounds such as acetaldehyde, diacetyl, acetoin, acetone and 2-butanone [4,5]. Tian and co-workers reported that the appealing flavour of yoghurt depended on the appropriate concentration of diacetyl, acetaldehyde and acetoin, with the optimal ratio of 4:16:32 mg/L, respectively [6]. Other VOCs, including aldehydes and ketones produced by lipid oxidation, as well as fatty acids and the products of the Maillard reaction (e.g., furan derivatives), have been associated with off-flavours in yoghurts [5,6,7].
In former times, the yeasts Candida famata (teleomorph Debaryomyces hansenii, formally Torulopsis candida), Zygosaccharomyces bailii (formerly Saccharomyces bailii), Saccharomyces delbrueckii var. mongolicus (syn. Torulaspora delbrueckii) and Saccharomyces cerevisiae were used in skyr manufacture to create a product with a unique texture and flavour, lower acidity and a tiny amount of alcohol [3]. In the modern dairy industry, yeasts are rather considered as contaminants and potential spoilage microorganisms. The common yeast spoilers in yoghurt and skyr belong to the genera Kluyveromyces, Debaryomyces, Pichia, Candida and Saccharomyces. Whether a given yeast species can be regarded as a contaminant depends on its ability to propagate at refrigerated temperatures, acidic pH and low water activity during the storage of the dairy products [8,9]. When grown in high numbers (typically more than 105 CFU/g), yeasts can cause quality deterioration and shorten the product’s shelf life. The most common quality defects, resulting from the metabolic activity of yeasts, include the production of gas, off-flavours, discolouration and changes in texture [10,11].
Kluyveromyces marxianus is well known for its biotechnological applications, such as the production of various enzymes, bioethanol, phenylethanol and fructose [12]. Due to the respiratory-fermentative type of metabolism, similar to Saccharomyces cerevisiae, K. marxianus has been successfully used as a baker’s yeast for the improvement of bread sensorial qualities [13,14,15]. In yoghurt-type products, K. marxianus is considered as a major spoilage yeast [10,11,12,13,14,15,16]. Though K. marxianus is designated as GRAS (generally recognized as safe) and rarely associated with health concerns, the quality and preservation of the dairy products contaminated with this species would still be of concern. It is well known that K. marxianus has the ability to metabolize lactose, citrate, protein and fat in milk, and produce a variety of secondary metabolites (e.g., acetaldehydes, carbon dioxide and ethanol), which might impart aroma development in dairy products [17,18]. Total sensory evaluation revealed an extensive negative impact of K. marxianus on the sensory qualities of yoghurts, which was even greater than the non-lactose-fermenting yeast D. hansenii [19]. Fermentation of residual lactose within the yoghurt curd by K. marxianus and the release of carbon dioxide formed cracks, loosening the curd texture [19]. Zhang and co-workers reported that K. marxianus promoted proteolysis in fermented milk, generating high amounts of free amino acids and small peptides [20]. Other yeast species used in this study, P. kudriavzevii (formerly Issatchenkia orientalis and anamorph Candida krusei) and Torulaspora delbrueckii, have been described as non-lactose-fermenting yeasts, sporadically found in kefir grains, yoghurts, cheeses and African fermented foods [21,22,23]. Due to a strong proteolytic activity, P. kudriavzevii can be part of the ripening cultures in some cheeses (e.g., Kazak artisanal cheese in China and Armanda cheese in Spain), producing a range of desired aroma compounds [24]. In other dairy products, specifically at a low pH, P. kudriavzevii might form surface biofilms and present a major cause of fermentative spoilage [17].
Most of the published studies have been primarily focused on yeast biocontrol, growth characteristics and identification in dairy products. Only a few studies with fermented milks, particularly cheeses, examined the production of flavour compounds by K. marxianus. To our knowledge, data on the impact of yeasts on the quality and aroma development in skyr during storage have not been published so far. The aim of the present study was to evaluate the spoilage potential of the yeasts K. marxianus, P. kudriavzevii and T. delbrueckii, grown in skyr at a cold storage temperature, with a focus on the production of volatile aroma compounds, organic acids and other metabolites related to off-flavours and other quality defects.

2. Materials and Methods

2.1. Yeast Strains and Reagents

The strains of Kluyveromyces marxianus Km1, Km2 and Km3 were isolated from contaminated batches of skyr provided by a Danish dairy. Before isolation, skyr was incubated at 25 °C for 2 weeks to promote the growth of yeast contaminants. Afterwards, 20 g of skyr was mixed with 180 mL peptone–saline solution (SPO; 5.0 g NaCl; 0.3 g Na2HPO4·2H2O; 1.0 g Bacto Peptone per litre of distilled water; pH 5.6), homogenized in a stomacher (Struers Chem A/S, Hovedstaden, Denmark) for 120 s at high speed and plated on malt yeast glucose peptone agar (MYGP; 3.0 g malt extract; 3.0 g yeast extract; 10 g glucose; 5.0 g Bacto peptone; 15 g agar, pH 5.6). Then, 100 mg·L−1 chloramphenicol and 50 mg·L−1 chlortetracycline was added to inhibit bacterial growth. MYGP agar plates were incubated at 25 °C for 48 h. Yeast strains Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1 were previously isolated from skyr production and provided by Arla Foods (Copenhagen, Denmark). Yeast isolates were stored at −80 °C in MYGP broth containing 20% (w/v) glycerol. All reagents for the experiments were purchased from Sigma-Aldrich (Copenhagen, Denmark) or Thermo Fisher Scientific (Copenhagen, Denmark) unless otherwise stated.

2.2. Sequencing of 26S rRNA Gene and 5.8S-ITS Region

Total yeast DNA was extracted from the colony material using an InstaGene Matrix DNA extraction kit (Bio-Rad, Hercules, CA, USA). The PCR was performed in a 50 μL mixture, containing 25 μL of Taq DNA Polymerase 2× Master Mix RED (Ampliqon, Odense, Denmark), 5 μL primer mix and 3 μL of the sample DNA. The D1/D2 region of the 26S rRNA gene was amplified, using the primers NL-1 (5′-GCATATCAATAAGCGGAGGAAAAG-3′) and NL-4 (5′-GGTCCGTGTTTCAAGACGG-3′) [25]. The PCR temperature cycling conditions included initial denaturation at 95 °C for 5 min, 35 cycles of 95 °C for 1 min, annealing at 52 °C for 45 s and extension at 72 °C for 1 min, followed by final extension at 72 °C for 7 min. Additionally, sequencing of the internal transcribed spacer region (ITS1-5.8S rRNA-ITS2) was applied to differentiate the species of K. marxianus and K. lactis. Primers used for amplification of the 5.8S-ITS region were ITS-1 (5′-TCC GTA GGT GAA CCT GCG G-3′) and ITS-4 (5′- TCC TCC GCT TAT TGA TAT GC-3′) [26]. The PCR conditions included initial denaturation at 95 °C for 3 min, 35 cycles of 95 °C for 30 s, 60 °C for 45 s and 72 °C for 45 s, followed by final extension at 72 °C for 10 min. The PCR products were sequenced by Macrogen Europe (Amsterdam, The Netherlands).
To identify yeast isolates, sequences were aligned with 26S rRNA gene and ITS1-5.8S rRNA-ITS2 region sequences in the NCBI Genbank database, using the basic local alignment search tool (BLAST) (http://www.ncbi.nlm.nih.gov/BLAST/, accessed on 20 September 2021).

2.3. Pulsed-Field Gel Electrophoresis

Pulsed-field gel electrophoresis (PFGE) was performed to characterise the chromosome length polymorphism among the dairy yeast isolates, compared to the type strains of K. marxianus CBS1553 and K. lactis CBS845 [27]. In brief, yeast chromosomal DNA was loaded into agarose plugs using the CHEF Yeast Genomic DNA plug kit (BIO-RAD, Coppenhagen, Denmark). PFGE was carried out on a PFGE DR-III unit (BIO-RAD, Hercules, CA, USA) at 10 °C in TBE buffer (45 mM Tris-base, 44 mM boric acid and1 mM EDTA). The running conditions were 150 V at a 200 s switch interval for 24 h, followed by 100 V at a 700 s switch interval for 48 h. Saccharomyces cerevisiae 345S (BIO-RAD, Hercules, CA, USA) and Hansenula wingei (BIO-RAD, Hercules, CA, USA) were used as chromosomal markers. The gel was stained in ethidium bromide and visualized with UV transillumination (Alpha-InnoTec GmbH, Kasendorf, Germany). BioNumerics software (version 2.50, Applied Maths, Kortrijk, Belgium) was used for cluster analysis of the DNA profiles.

2.4. Yeast Propagation in Skyr

Skyr (natural skyr, 0.2% fat and 11% protein) was received from a Danish dairy. Skyr was produced with the use of the starter cultures Lactobacillus delbrueckii subsp. bulgaricus and Streptococcus thermophilus, and a bioprotective culture Lacticaseibacillus paracasei. Before the experiments, yeasts were propagated in MYGP broth at 25 °C overnight, washed and resuspended in SPO. Skyr was aseptically distributed in stomacher bags (200 g in each) and inoculated with yeast cell suspensions (2 mL) to a final concentration of 104 CFU·g−1 or with SPO (2 mL) for the control experiments. The mixtures were homogenized in a stomacher (Struers Chem A/S, Coppenhagen, Denmark) for 120 s at high speed, distributed into 50 mL sterile polypropylene boxes (Berry Superfoss 618 PP plast, Dansk Transport Emballage A/S, Vojens, Denmark) and incubated at 8 °C. Samples were withdrawn at the start (Day 0) and after 7, 14 and 21 days of incubation. Samples were stored at −20 °C before the analyses.

2.5. Measurements of Viable Counts, pH and Carbon Dioxide

Viable counts of yeasts in skyr were determined by plating serial dilutions in SPO onto MYGP agar plates with added antibiotics (100 mg·L−1 chloramphenicol and 50 mg·L−1 chlortetracycline), which were incubated at 25 °C for 48 h. Lactic acid bacteria (LAB) were enumerated using de Man, Rogosa Sharpe (MRS, pH 6.2), containing an antifungal agent (0.01% cycloheximide). The MRS agar plates were incubated anaerobically using AnaeroGen sachet (Thermo Fisher Scientific, Roskilde, Denmark) at 37 °C for 48 h. The pH in skyr samples was recorded using a digital pH meter (Mettler Toledo, Glostrup, Denmark). Production of CO2 was calculated from the total weight loss of skyr during incubation and expressed in nmol·g−1 of CO2 [28].

2.6. GC-MS Analysis of the Volatile Compounds and Data Processing

Volatile compounds in skyr were analysed using headspace gas chromatography coupled to mass spectrometry (GC-MS), as previously described [29]. Briefly, a mixture of 30 g of skyr, 30 mL of tap water and 1 mL of internal standard (4-methyl-1-pentanol, 5 mg·L−1) was homogenized using an Ultra-Turrax T25 homogenizer (Janke & Kunkel, Staufen, Germany) at 12,000 rpm for 90 s. VOCs were collected onto Tenax-TA traps (Markes International, Llantrisant, UK) by purging with nitrogen (100 mL·min−1) for 45 min and a dry purge of the traps for 10 min. The volatiles were desorbed by an automatic thermal desorption unit (TurboMatrix 350, Perkin Elmer, Shelton, CT, USA) into a GC-MS (7890A GC-system interfaced with a 5975C VL MSD with triple-Axis detector, Agilent Technologies, Palo Alto, California) and separated on a DB-Wax capillary column. The mass spectra were recorded at 70 eV within a mass/charge ratio of 15–300 m/z. Data were extracted from the chromatograms using the PARADISe software [30]. Volatile compounds were identified from the retention indices (RI), and calculated using a series of n-alkanes (C6–C22) and the NIST05 database. Identifications were confirmed by comparison with the RI of authentic reference compounds or published RI values. Semi-quantitative data for volatiles were obtained by multiplying the internal standard concentration with the ratio of the peak area of the volatile compound to the internal standard peak area [31].

2.7. High Performance Liquid Chromatography

Organic acids (lactic, citric, acetic, orotic and uric acids), carbohydrates (lactose, glucose and galactose) and ethanol were analysed using high performance liquid chromatography (HPLC) [32,33]. Briefly, 5 g of the skyr test samples and controls were diluted with 0.01 M H2SO4 to 25 mL and centrifuged at 4000× g, 20 min at 4 °C (Eppendorf, Hamburg, Germany). The supernatant was collected by filtering through a 0.2 μm membrane (Sarstedt, Germany) and stored at −20 °C before the analysis. The chromatographic separation was conducted with 10 mM H2SO4 as a mobile phase at a flow rate of 0.6 mL·min−1 on an Aminex HPX-87H column (300 mm × 7.8 mm). The HPLC equipment consisted of a quaternary pump, an online degasser, a UV-visible detector and a refractive index detector (RID; Agilent 1100 series, Waldbronn, Germany). The UV detector was set at 210 nm (citric, lactic and acetic acids) and 290 nm (orotic and uric acid) for the detection of organic acids. The RID detector was used for the analyses of carbohydrates. Data were collected and processed using Agilent ChemStation software (Agilent Technologies, Waldbronn, Germany). Quantification of organic acids was carried out using external standards of the relevant compounds. The calibration curves were constructed for the individual organic acids by plotting the peak heights after duplicate injections of standard solutions, at a broad range of concentrations. Identifications were made from the retention times of the individual standards.

2.8. Phenotypic Tests

Yeasts were tested for the ability to ferment and assimilate glucose, lactose and galactose [34]. Thereafter, yeast colony material was resuspended in the media for fermentation tests (0.45% w/v yeast extract, 0.75% w/v peptone, 2.66 mg·100 mL−1 bromothymol blue, 2% w/v carbohydrate added, pH 6.0 ± 0.2) and for assimilation tests (0.067% w/v Difco Yeast Nitrogen Base, 0.05% w/v carbohydrate added, pH 5.6 ± 0.1). Durham tubes were used to detect gas production in fermentation tests. Gas production and colour change of the indicator (bromothymol blue) was recorded weekly. Yeast growth was scored based on the turbidity of the cell suspensions using Wickerham cards.

2.9. Statistical Analysis

All experiments were conducted in triplicate. Significant differences between the results (confidence level of 95%) were tested using one-way analysis of variance (ANOVA) followed by a Tukey’s post-hoc test (JMP software, SAS Institute Inc., Cary, NC, USA) and the unpaired t-test (GraphPad Prism 9, GraphPad Software, San Diego, CA, USA).

2.10. Nucleotide Sequence Accession Number

The sequences of the D1/D2 region of the 26S rRNA gene were deposited in the DDBJ, EMBL and GenBank nucleotide sequence databases under accession numbers OM854769-OM854773 (Table 1). Sequences of the ITS1-5.8S rRNA gene-ITS2 fragments for K. marxianus Km1, Km2 and Km3 were deposited under accession numbers OM406350, OM406351 and OM406352, respectively.

3. Results

3.1. Yeasts Identification and Growth

Yeasts isolates from skyr were identified as K. marxianus, P. kudriavzevii and T. delbrueckii based on the 26S rRNA gene sequencing (Table 1). In addition, sequencing of the 5.8S-ITS region confirmed the strain annotation to K. marxianus (100% identity). Only minor differences in the chromosomal length polymorphism were observed for the strains Km1 and Km2, which could be easily separated from the strain Km3 (Supplementary Figure S1). Based on the PFGE profiles, all three dairy strains were, as expected, more closely related to the type strain of K. marxianus CBS1553 than to the type strain of K. lactis CBS845. All yeast strains were able to grow in skyr at 8 °C, reaching the stationary phase after 7 days of incubation (Figure 1). The highest counts, after 21 days of incubation, were observed for the K. marxianus strains (6.95–7.17 Log CFU·g−1), while P. kudriavzevii Pk1 and T. delbrueckii Td1 were grown to lower numbers (5.41 ± 0.15 and 5.65 ± 0.11 Log CFU·g−1, respectively). Viable counts of LAB in skyr were stable within the incubation period (9.55–9.74 Log CFU·g−1) and were not affected by yeasts. Similarly, pH was basically maintained during the storage period, showing a slight tendency to decrease from 4.16 ± 0.04 (Day 0) to 4.02 ± 0.03 (Day 21).

3.2. Production of Volatile Compounds

In total, 40 volatile aroma compounds of the main VOC types, including alcohols (12 compounds), aldehydes (eight compounds), esters (14 compounds) and ketones (six compounds), were identified in skyr during the growth of yeasts. The PCA plots (Figure 2) describe variation between the yeasts and incubation time (0, 7, 14 and 21 days) in VOC production. The first two components (PC1 and PC2) explained 65.4% of the variation in the data. Yeasts strains of K. marxianus (Km1, Km2 and Km3) were grouped together according to the incubation time and clearly discriminated from the controls. A separate group comprised P. kudriavzevii Pk1 samples after 14 and 21 days of incubation. Samples of T. delbrueckii Td1 were poorly distinguished from skyr controls, indicating similar VOC patterns.
Figure 3 presents the relative distribution of the VOC types, calculated by summing the averages of the relevant compounds at 0, 7, 14 and 21 days of incubation. The proportions of alcohols (10% at baseline) were gradually increased during the growth of K. marxianus strains (up to 76–78%) and P. kudriavzevii Pk1 (up to 45%), along with the amounts of esters (from 2% at baseline to 13–14% and 25%, respectively). Ketones, comprising highly abundant group at the start of incubation (80%), declined to 7–8% for K. marxianus strains and to 28% for P. kudriavzevii Pk1 at the end of incubation. The initial levels of aldehydes (7–9%) were increased after 7 days incubation of K. marxianus strains to 22–26%, followed by a decrease to 2% (Day 21), while for P. kudriavzevii Pk1, a gradual decrease in aldehydes was observed (to 2%). Incubation of T. delbrueckii Td1 led, primarily, to a slight increase in ketones (by 3–8%) on the account of aldehydes.
Concentration of the individual aroma compounds at the start and end of incubation (Day 0 and Day 21), retention indices and odour thresholds in water are shown in Table 2. The results of VOC quantification at Day 7 and 14 are listed in Supplementary Materials (Tables S1 and S2). Figure 4 shows the heatmap and hierarchical cluster analysis of the total VOC data. Alcohols and esters comprised the major VOC cluster in the heatmap (Cluster I, Figure 4). Most of the alcohols were highly increased (more than 10-fold at Day 21) for K. marxianus strains compared to controls, including 2-butanol (from 2.1 to 17–21 µg·kg−1), 1-hexanol (from 4.5 to 46–56 µg·kg−1), 3-methyl-1-butanol or isoamylalcohol (from 21 to 1780–1918 µg·kg−1), 2-methyl-1-propanol (from 0.6 to 446–591 µg·kg−1), 1-propanol (from 0.4 to 14–15 µg·kg−1) and 2-phenylethanol (from 0.4 to 3.2–4.7 µg·kg−1) (Table 2). Other alcohols (1-pentanol, 1-penten-3-ol, 1-octanol) were moderately increased (four- to nine-fold) for K. marxianus strains. The levels of 3-methyl-1-butanol and 1-pentanol were close to the threshold values. Along with the alcohols, ethyl esters of formic, lactic, propionic and fatty acids, as well as 3-methylbutyl acetate (isoamylacetate) and other acetic acid esters (ethyl-, hexyl, pentenyl etc.), were produced in significantly higher amounts by K. marxianus (Table 2, Tables S1 and S2). The major increases above the thresholds were recorded for ethyl acetate (from 0 to 278–291 µg·kg−1), ethyl butyrate (from 0.51 to 34–81 µg·kg−1), ethyl hexanoate (from 0.1 to 19–32 µg·kg−1), ethyl octanoate (from 0.6 to 10–23 µg·kg−1), and 3-methylbutyl acetate (from 0.2 to 29–57 µg·kg−1) (Table 2). Cluster II (Figure 4) was primarily composed of ketones (2-heptanone, 2-butanone and 2-propanone) and straight-chain aldehydes (butanal, heptanal and hexanal), which were present at high concentration at Day 0 and generally decreased during the incubation of K. marxianus and controls. At the end of the incubation period, 2-butanone and 2-propanone (acetone) were found in significantly lower amounts for the K. marxianus strains (between 40–62 µg·kg−1 and 8–19 µg·kg−1, respectively) compared to controls (136 µg·kg−1 and 83 µg·kg−1, respectively), while 2-nonanone was preserved in K. marxianus samples. Cluster III combined branched-chain aldehydes (2-methylbutanal, 3-methylbutanal and 2-methylpropanal), which were dramatically increased during the exponential growth of K. marxianus strains (Figure 4 and Table S1). The highest increase was recorded for 2-methylbutanal and 2-methylpropanal (from 5–7 µg·kg−1 at Day 0 to 237–283 µg·kg−1 and 156–231 µg·kg−1, respectively, at Day 7). The levels of these compounds at Day 21 were still two- to four-fold higher than in the controls and close to the threshold values (Table 2). Cluster IV ketones, 2,3-butanedione (diacetyl) and 3-hydroxy-2-butanone (acetoin) were decreased during the incubation of K. marxianus strains (from 178 to 48–49 µg·kg−1 and from 152 to 55–79 µg·kg−1, respectively), whereas they were accumulated in skyr controls.
Several differences in VOC production were observed between the K. marxianus strains. Thus, Km1 differed from other strains by the high production of ethyl butyrate, pentanol, hexanol and hexanal. Strain Km2 was associated with increased levels of several esters (ethyl hexanoate, ethyl octanoate and 2-phenylethyl acetate), while Km3 was distinguished by the amounts of 2-methyl propanal (highest) and ethyl propionate (lowest). Growth of P. kudriavzevii Pk1 was associated with the increased production of 2-butanol, 3-methyl-1-butanol, 2-propanol, ethyl propionate and ethyl acetate, together with decreases in ketones (2,3-butanedione, 2-butanone and 2-propanone). Only a few significant shifts in VOC production were observed for T. delbrueckii Td1, including increases in 1-penten-3-ol and 2-nonanone, as well as decreases in 1-butanol, 3-methyl-2-butenal and 2,3-butanedione, compared to controls (Table 2).

3.3. Utilization of Lactose and Galactose

The amount of lactose in skyr (initially 3.10 ± 0.16 g·100 g−1) was gradually decreased during incubation of K. marxianus strains (Figure 5A). The lowest concentration of lactose was found for Km1 and Km2 (0.84 ± 0.03 and 0.95 ± 0.03 g·100 g−1, respectively, at Day 21). Lactose concentration for P. kudriavzevii Pk1 and T. delbrueckii Td1 was reduced at Day 7 and remained stable afterwards (2.26 ± 0.16 and 2.16 ± 0.28 g·100 g−1 at Day 21, respectively; Figure 5A). Concurrently, the levels of galactose (initially 0.37 ± 0.02 g·100 g−1) were reduced in yeast incubations and controls, largest for K. marxianus Km1 (0.09 ± 0.01 g/100 g−1 at Day 21; Figure 5B). The concentration of glucose was negligible and not detected by the applied HPLC. Phenotypic tests showed that K. marxianus strains had an ability to ferment and assimilate lactose and galactose, in contrast to P. kudriavzevii Pk1 and T. delbrueckii Td1 (Table S3).

3.4. Production of Organic Acids, Ethanol and Carbon Dioxide

Table 3 shows the amounts of organic acids in skyr during yeast incubation. The major acids at Day 0 were lactic (224 ± 8.5 mg·100 g−1), citric (120 ± 9.1 mg·100 g−1) and acetic (54 ± 1.4 mg·100 g−1) acids. Additionally, orotic and uric acids were detected in lower amounts (5.1 ± 0.16 mg·100 g−1 and 3.7 ± 0.36 mg/100 g−1, respectively). The levels of lactic acid were generally increased by yeasts and in controls (257–285 mg·100 g−1 at Day 21). Citric acid was significantly reduced by yeasts already after 7 days incubation (73–86 mg·100 g−1) and remained stable afterwards, except for K. marxianus Km1 (further decrease to 48 ± 8.1 mg·100 g−1 at Day 21). Similarly, orotic acid was mostly decreased by yeasts, largest by Km1 (3.3 ± 0.14 mg·100 g−1 at day 21), while the levels of uric acid were generally stable over time. The reduction in citric and orotic acids, though to a lesser extent, was also observed in skyr controls. Differences between the yeast species were found for acetic acid at Day 21, which was decreased for K. marxianus Km1 and Km2 (48 ± 2.4 and 50.4 ± 2.5 mg·100 g−1, respectively), opposite to P. kudriavzevii Pk1 and T. delbrueckii Td1 (increased to 60.2 ± 2.9 and 58 ± 4.1 mg·100 g−1, respectively).
Ethanol was detected for K. marxianus at Day 21 in a range of 0.11–0.13%; other samples tested negative for ethanol. Carbon dioxide was continuously produced by K. marxianus Km1, Km2 and Km3, reaching concentrations of 7.20 ± 1.26, 7.49 ± 1.01 and 8.30 ± 0.31 nmol·g−1, respectively, at Day 21 (Figure 6). Production of CO2 by P. kudriavzevii Pk1 and T. delbrueckii Td1 was comparatively low (0.60 ± 0.25 and 0.69 ± 0.09 nmol·g−1 at Day 21, respectively).

4. Discussion

This study investigated the spoilage activity of the yeasts K. marxianus (Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1, proliferating in skyr in cold storage. We were particularly focused on the most common spoilage species K. marxianus, of which strains were isolated from skyr and identified in this study. As expected, all strains of K. marxianus were able to grow in skyr to the levels close to those reported for yoghurts (6–8 Log CFU·g−1) at refrigerated temperatures [22,35,36]. High counts of K. marxianus were due to the ability of this species to ferment lactose and galactose in skyr, with the production of CO2 and ethanol, promoted by oxygen limitation conditions. These results are consistent with earlier studies on yoghurts [36,37,38]. In contrast to K. marxianus, P. kudriavzevii Pk1 and T. delbrueckii Td1 were unable to utilize lactose and galactose, leading to poorer growth and metabolite release.
The flavour of skyr, likewise with other fermented milk products, is created by a balanced mixture of various essential VOCs. It is well recognized that the shifts in VOC amounts, rather than the absence or presence of a particular compound, give rise to off-flavours in dairy products during storage [39,40]. In this study, VOC profiles in skyr highly depended on the yeast species and incubation time. The absolutely largest diversity and changes in VOCs were observed for the strains of K. marxianus, in accordance with their capacity to grow in skyr to high numbers. The impact of P. kudriavzevii Pk1 on individual VOCs was less pronounced, while only minor significant differences from controls were recorded for T. delbrueckii Td1, indicating that the spoilage potential is a combination of both cell numbers and the off-flavour profile. The following discussion will focus mostly on VOC production by K. marxianus in relation to the undesired aroma and spoilage potential.
Higher (fusel) alcohols were accumulated during K. marxianus propagation in skyr, comprising a major VOC fraction. Higher alcohols are produced by yeasts via the Ehrlich pathways through conversion of α-keto acids into aldehydes, which can be either further reduced to higher alcohols or dehydrogenated to organic acids [41]. The activity of the Ehrlich pathways is known to be growth-phase dependent [41,42]. Accordingly, shifts in the levels of aldehydes were recorded at the exponential growth phase of K. marxianus (0–7 days). A reduction in the straight chain aldehydes and temporary increase in the branched chain aldehydes was accompanied by an increase in respective alcohols (e.g., 1-hexanol, 3-methyl-1-butanol and 2-methyl-1-propanol). The increase in the alcohols continued through the stationary growth (14–21 days), indicating that the conversion of aldehydes into alcohols was still active.
The high production of fusel alcohols indicated fermentative dissimilation of lactose by K. marxianus, taking place at oxygen limitation conditions. Generally, alcohols have lower odour activity compared to aldehydes and are considered to be of less importance for yoghurt aroma [5,43]. The individual alcohols in this study were detected in amounts lower than the thresholds (100–40,000 µg·kg−1). Despite this, the impact of alcohols on the overall skyr aroma should not be undervalued, considering their considerable accumulative levels in K. marxianus samples. Specifically, alcohols with relatively low thresholds, such as highly abundant isoamylalcohol (whiskey, malt aroma) and 2-methyl-1-propanol (alcoholic aroma), as well as the less abundant 1-pentanol (fusel) and 1-octanol (orange), might contribute to the off-flavours in skyr. Supporting this assumption, Güler reported a correlation between increases in isoamylalcohol and 2-methyl-1-propanol with a flavour decline in salted yogurt during storage [44]. One of the functions of the Ehrlich pathways is to form quorum-sensing compounds that can induce differentiation and participate in the adaptation of yeast cells to environmental changes [41]. Alcohol 2-phenylethanol, produced by K. marxianus in moderate amounts, is the most well-studied quorum-sensing molecule in dimorphic yeasts (e.g., K. marxianus and S. cerevisiae) [45,46]. Balbino et al. showed that exposure of K. marxianus to 2-phenylethanol provoked changes in cell morphology and fermentative metabolism, and impaired yeast viability [47]. Isoamylalcohol is another signalling compound, which has been reported to induce pseudohyphal growth and filamentation in S. cerevisiae [48]. Regarding the major release of isoamylalcohol by K. marxianus in this study, the significance of this compound and 2-phenylethanol for yeast growth and differentiation during skyr spoilage needs to be further examined.
Aldehydes, even at low levels, may affect a product’s flavour due to their low threshold values. The branched chain aldehydes, methylbutanals (apple-like, malt aroma) and 2-methylpropanal (sharp, pungent aroma), along with hexanal (sharp, fruity aroma) and heptanal (fatty, pungent aroma) would probably account for the initial off-flavours in skyr contaminated with K. marxianus. The same compounds have been referred as primary sources of undesired malty, nutty flavours in cheese [49,50,51]. Carrillo-Carrión and co-workers proposed that the total concentration of volatile aldehydes (C5-C9) can be used as a marker of yogurt degradation at inadequate storage conditions [52]. Among them, hexanal was positively associated with rancidity in foods and negatively associated with the “degree of likeness” in yoghurts [53,54].
Esters composed the second largest group of aroma compounds produced by K. marxianus. Esters are formed by a condensation reaction between alcohols (ethanol or higher alcohols) and acetyl-CoA (for acetate esters) or acyl-CoA (for ethyl esters) [55,56]. It is well known that ester production is positively affected by precursor availability, i.e., the alcohols released by K. marxianus in this study, and the presence of lactose in skyr required for acetyl-CoA formation [43,57]. Esters are characterized by relatively low threshold values and typically impart a fruity aroma to the dairy products [43,44,58]. In this study, ethyl acetate (solvent aroma), ethyl butyrate (pineapple), ethyl hexanoate (aniseed, apple), ethyl octanoate (sour apple) and isoamylacetate (banana, pineapple) were produced in large quantities, close to or higher than the thresholds. The ability of K. marxianus strains to produce esters, in particular ethyl acetate, has been frequently reported [38,56,59,60]. Previously, Leclercq-Perlat et al. identified similar aroma compounds produced by K. marxianus in a cheese medium (e.g., hexanal, 2 (3)-methylbutanal, butanol, 2 (3)-methylbutanol, 2-methylpropanol) and reported a link between the levels of isoamylacetate and ethyl acetate, and ester-like fruity odours [61].
The vicinal ketones diacetyl and acetoin, together with methyl ketones 2-butanone and 2-propanone, were mostly decreased by K. marxianus strains and P. kudriavzevii Pk1, in contrast to the skyr controls (accumulated or maintained). These ketones have been referred to as the key aroma compounds, associated with butter-like, creamy and fruity flavours in yoghurts [7,53,62,63]. The diacetyl flavour is particularly important due to a low threshold value [5,64]. Diacetyl can be produced by lactic acid bacteria and yeasts from lactose by oxidative decarboxylation of α-acetolactate, and further converted to acetoin by intracellular reductases [5,65]. Regarding the low reductase activity reported for the starter cultures S. thermophilus and L. bulgaricus, accumulation of diacetyl and acetoin in skyr was probably enhanced by the metabolic activity of L. paracasei (bioprotective culture in skyr) [66,67]. Lower levels of diacetyl and acetoin for the K. marxianus strains and P. kudriavzevii Pk1 can be explained by the redirection of pyruvate metabolism towards the production of carbon dioxide and ethanol, and/or towards the Ehrlich pathways [65]. Reduction in methyl ketones in yeast fermentations might be due to their conversion to the corresponding alcohols (2-propanol and 2-butanol). Furthermore, the presence of lactose, an easily assimilated carbon source, might result in repression of lipid oxidation and, consequently, the production of methyl ketones [61].
Different strains of K. marxianus exhibited, more or less, similar VOC patterns with a few variations in specific compounds. Regarding these variations, K. marxianus Km1 is expected to be associated with more fusel alcohol and a pungent flavour (higher release of 1-hexanol, 1-pentanol, and ethyl butyrate), while K. marxianus Km2 will probably impart a hint of fruity apple-like odours (larger amounts of ethyl hexanoate and ethyl octanoate). Off-flavours related to P. kudriavzevii Pk1 may occur from increases in isoamyl acetate and ethyl acetate (both higher than their threshold values), along with a decline in ketones and ethyl esters. Although changes in VOCs for T. delbrueckii Td1 were mostly insignificant, the combined impact, e.g., as shifts in essential ketones, diacetyl, 2-heptanone and 2-nonanone, should not be excluded. It should be noted, however, that thresholds in this study has been reported in water and would therefore differ from the actual values, since the main components in skyr (carbohydrates, proteins and fat) would counteract VOC release [49,68,69]. Thus, further sensory evaluation of aroma perception in skyr is needed to draw conclusions of yeast-related off-flavours.
Organic acids contribute to the sensory characteristics and act as preservatives in fermented milk products [70,71,72]. Citric acid or citrate is normally present in raw milk in high amounts and can be degraded by LAB to acetate, lactate, ethanol and acetoin as the major end products [4,73]. The largest reduction in citric acid was observed in yeast incubations, especially for K. marxianus Km1, indicating strain-specific differences in citrate utilization. These results are in accordance with previous findings, showing the ability of some K. marxianus strains to metabolize citric acid [74,75]. Acetic acid is produced from pyruvate by LAB and yeasts via heterolactic fermentation pathways [73]. In accordance with VOC analysis, lower amounts of acetic acid in K. marxianus incubations can be related to acetate conversion into esters, and/or its utilization as a carbon source for the production of ethanol and carbon dioxide. Orotic acid is an intermediate compound in pyrimidine biosynthetic pathways that can be used as a growth factor by the starter cultures [70,76]. Thus, a reduction in orotic acid by yeasts in the exponential growth phase indicated the activation of nucleotide biosynthesis essential for cell growth. Similar to this study, a minor decrease in orotic acid has been reported in yoghurts in cold storage [76,77,78].

5. Conclusions

Spoilage yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 were able to proliferate in skyr during refrigerated storage, exhibiting species-specific growth characteristics and metabolite production. Strains of K. marxianus were characterized by the largest viable counts and major release of various aroma compounds. Many of the VOCs, produced in amounts higher than the threshold values, were previously related to off-flavours and product deterioration. In contrast to the other species, the growth of K. marxianus in skyr resulted in the production of carbon dioxide and ethanol. Extensive growth capacity, and the formation of undesired aroma and carbon dioxide indicated a high spoilage potential of K. marxianus contaminants. The impact of other yeast species on skyr aroma was less pronounced. In particular, T. delbrueckii Td1 had low spoilage activity, showing only minor changes in metabolite profiles, as compared to non-contaminated skyr. The current study clarifies metabolic activities in yeast species, associated with the spoilage of skyr, and highlights the importance of distinguishing between contaminating yeast species, and even strains, when evaluating their impact on product quality and shelf-life. A paradigm shift in yeast quality control of fermented milk products is recommended, moving from solely determining colony forming units towards more detailed studies, focusing on yeast species identification, growth patterns and metabolite profiles.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods11121776/s1, Figure S1: Dendrograms, showing clustering of Kluyveromyces marxianus Km1, Km2, Km3 and the type strains of Kluyveromyces marxianus CBS1553 and Kluyveromyces lactis CBS845; Table S1: Concentration of volatile compounds (VOC) in skyr after 7 days incubation of yeasts Kluyveromyces marxianus Km1, Km2 and Km3, Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1, and in skyr controls (without added yeasts); Table S2 Concentration of volatile compounds (VOC) in skyr after 14 days incubation of yeasts Kluyveromyces marxianus Km1, Km2 and Km3, Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1, and in skyr controls (without added yeasts), and Table S3: Fermentation and assimilation of lactose, glucose and galactose by yeasts Kluyveromyces marxianus Km1, Km2 and Km3, Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1.

Author Contributions

Experimental work, data analysis and original draft preparation, T.S. and N.L. editing, supervision and project administration, L.J.; methodology and data analysis, M.A.P., conceptualization and review, S.K.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Danish Dairy Research Foundation (project “Improve dairy life”), and the Royal Thai Government Scholarship.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arla Foods Ingredients. Skyr is predicted a promising future in Europe. Int. Dairy Mag. 2020, 15, 1–36. [Google Scholar]
  2. Gudmundsson, G.; Kristbergsson, K. Modernization of skyr processing: Icelandic acid-curd soft cheese. In Modernization of Traditional Food Processes and Products, 1st ed.; McElhatton, A., el Idrissi, M.M., Eds.; Springer: New York, NY, USA, 2016; pp. 45–53. [Google Scholar] [CrossRef]
  3. Valsdóttir, Þ.; Sveinsson, Þ.E. Uniqueness of Traditional Skyr; Matís: Reykjavík, Iceland, 2011. [Google Scholar]
  4. Chen, C.; Zhao, S.; Hao, G.; Yu, H.; Tian, H.; Zhao, G. Role of lactic acid bacteria on the yogurt flavour: A review. Int. J. Food Prop. 2017, 20, S316–S330. [Google Scholar] [CrossRef] [Green Version]
  5. Cheng, H. Volatile flavor compounds in yogurt: A review. Crit. Rev. Food Sci. Nutr. 2010, 50, 938–950. [Google Scholar] [CrossRef] [PubMed]
  6. Tian, H.; Yu, B.; Yu, H.; Chen, C. Evaluation of the synergistic olfactory effects of diacetyl, acetaldehyde, and acetoin in a yogurt matrix using odor threshold, aroma intensity, and electronic nose analyses. J. Dairy Sci. 2020, 103, 7957–7967. [Google Scholar] [CrossRef]
  7. Miyaji, K.; Kuwano, Y.; Murakami, Y.; Hirata, S.; Imayoshi, Y.; Maruyama, H.; Koizumi, R.; Inoue, H.; Azuma, N. Off-flavors generated during long-term ambient storage of pasteurized drinking yogurt from skim milk. Biosci. Biotechnol. Biochem. 2020, 85, 391–400. [Google Scholar] [CrossRef]
  8. Garnier, L.; Valence, F.; Pawtowski, A.; Auhustsinava-Galerne, L.; Frotté, N.; Baroncelli, R.; Deniel, F.; Coton, E.; Mounier, J. Diversity of spoilage fungi associated with various French dairy products. Int. J. Food Microbiol. 2017, 241, 191–197. [Google Scholar] [CrossRef]
  9. Awasti, N.; Anand, S. The role of yeast and molds in dairy industry: An update. In Dairy Processing: Advanced Research to Applications, 1st ed.; Minj, J., Sudhakaran, A.V., Kumari, A., Eds.; Springer: Singapore, 2020; pp. 243–262. [Google Scholar] [CrossRef]
  10. Garnier, L.; Valence, F.; Mounier, J. Diversity and Control of Spoilage Fungi in Dairy Products: An Update. Microorganisms 2017, 5, 42. [Google Scholar] [CrossRef] [Green Version]
  11. Geronikou, A.; Srimahaeak, T.; Rantsiou, K.; Triantafillidis, G.; Larsen, N.; Jespersen, L. Occurrence of Yeasts in White-Brined Cheeses: Methodologies for Identification, Spoilage Potential and Good Manufacturing Practices. Front. Microbiol. 2020, 11, 582778. [Google Scholar] [CrossRef]
  12. Karim, A.; Gerliani, N.; Aïder, M. Kluyveromyces marxianus: An emerging yeast cell factory for applications in food and biotechnology. Int. J. Food Microbiol. 2020, 333, 108818. [Google Scholar] [CrossRef]
  13. Caballero, R.; Olguín, P.; Cruz-Guerrero, A.; Gallardo, F.; García-Garibay, M.; Gómez-Ruiz, L. Evaluation of Kluyveromyces marxianus as baker’s yeast. Food Res. Int. 1995, 28, 37–41. [Google Scholar] [CrossRef]
  14. Struyf, N.; Laurent, J.; Verspreet, J.; Verstrepen, K.J.; Courtin, C.M. Saccharomyces cerevisiae and Kluyveromyces marxianus Cocultures Allow Reduction of Fermentable Oligo-, Di-, and Monosaccharides and Polyols Levels in Whole Wheat Bread. J. Agric. Food Chem. 2017, 65, 8704–8713. [Google Scholar] [CrossRef] [Green Version]
  15. Sá, A.G.A.; Meneses, A.; de Araújo, P.H.H.; Oliveira, D. A review on enzymatic synthesis of aromatic esters used as flavor ingredients for food, cosmetics and pharmaceuticals industries. Trends Food Sci. Technol. 2017, 69, 95–105. [Google Scholar] [CrossRef]
  16. Mayoral, M.B.; Martín, R.; Sanz, A.; Hernández, P.E.; González, I.; García, T. Detection of Kluyveromyces marxianus and other spoilage yeasts in yoghurt using a PCR-culture technique. Int. J. Food Microbiol. 2005, 105, 27–34. [Google Scholar] [CrossRef]
  17. Fleet, G.H. Yeast Spoilage of Foods and Beverages. In The Yeasts, 5th ed.; Kurtzman, C.P., Fell, J.W., Boekhout, T., Eds.; Elsevier B.V.: Amsterdam, The Netherlands, 2011; pp. 53–63. [Google Scholar] [CrossRef]
  18. Fonseca, G.G.; Heinzle, E.; Wittmann, C.; Gombert, A.K. The yeast Kluyveromyces marxianus and its biotechnological potential. Appl. Microbiol. Biotechnol. 2008, 79, 339–354. [Google Scholar] [CrossRef]
  19. Šalomskienė, J.; Mačionienė, I. The influence of contamination yoghurt, quark and semi-hard cheese by yeasts on their sensory properties. Vet. Med. Zoot 2009, 48, 72–76. [Google Scholar]
  20. Zhang, D.-D.; Liu, J.-L.; Jiang, T.-M.; Li, L.; Fang, G.-Z.; Liu, Y.-P.; Chen, L.-J. Influence of Kluyveromyces marxianus on proteins, peptides, and amino acids in Lactobacillus-fermented milk. Food Sci. Biotechnol. 2017, 26, 739–748. [Google Scholar] [CrossRef]
  21. Maïworé, J.; Ngoune, L.T.; Piro-Metayer, I.; Montet, D. Identification of yeasts present in artisanal yoghurt and traditionally fermented milks consumed in the northern part of Cameroon. Sci. Afr. 2019, 6, e00159. [Google Scholar] [CrossRef]
  22. Lourens-Hattingh, A.; Viljoen, B.C. Survival of dairy-associated yeasts in yoghurt and yoghurt-related products. Food Microbiol. 2002, 19, 597–604. [Google Scholar] [CrossRef]
  23. Johansen, P.G.; Owusu-Kwarteng, J.; Parkouda, C.; Padonou, S.W.; Jespersen, L. Occurrence and Importance of Yeasts in Indigenous Fermented Food and Beverages Produced in Sub-Saharan Africa. Front. Microbiol. 2019, 10, 1789. [Google Scholar] [CrossRef] [Green Version]
  24. Zheng, X.; Li, K.; Shi, X.; Ni, Y.; Li, B.; Zhuge, B. Potential characterization of yeasts isolated from Kazak artisanal cheese to produce flavoring compounds. MicrobiologyOpen 2017, 7, e00533. [Google Scholar] [CrossRef] [Green Version]
  25. Gori, K.; Ryssel, M.; Arneborg, N.; Jespersen, L. Isolation and Identification of the Microbiota of Danish Farmhouse and Industrially Produced Surface-Ripened Cheeses. Microb. Ecol. 2012, 65, 602–615. [Google Scholar] [CrossRef] [Green Version]
  26. Kühle, A.V.D.A.; Jespersen, L. The Taxonomic Position of Saccharomyces boulardii as Evaluated by Sequence Analysis of the D1/D2 Domain of 26S rDNA, the ITS1-5.8S rDNA-ITS2 Region and the Mitochondrial Cytochrome-c Oxidase II Gene. Syst. Appl. Microbiol. 2003, 26, 564–571. [Google Scholar] [CrossRef]
  27. Jespersen, L.; Nielsen, D.S.; Hønholt, S.; Jakobsen, M. Occurrence and diversity of yeasts involved in fermentation of West African cocoa beans. FEMS Yeast Res. 2005, 5, 441–453. [Google Scholar] [CrossRef] [Green Version]
  28. Renault, P.E.; Albertin, W.; Bely, M. An innovative tool reveals interaction mechanisms among yeast populations under oenological conditions. Appl. Microbiol. Biotechnol. 2013, 97, 4105–4119. [Google Scholar] [CrossRef]
  29. Zhang, L.; Huang, C.; Johansen, P.G.; Petersen, M.A.; Poojary, M.M.; Lund, M.N.; Jespersen, L.; Arneborg, N. The utilisation of amino acids by Debaryomyces hansenii and Yamadazyma triangularis associated with cheese. Int. Dairy J. 2021, 121, 105135. [Google Scholar] [CrossRef]
  30. Johnsen, L.G.; Skou, P.B.; Khakimov, B.; Bro, R. Gas chromatography—Mass spectrometry data processing made easy. J. Chromatogr. A 2017, 1503, 57–64. [Google Scholar] [CrossRef]
  31. Xiao, J.; Chen, Y.; Li, J.; Shi, X.; Deng, L.; Wang, B. Evaluation of the Effect of Auxiliary Starter Yeasts With Enzyme Activities on Kazak Cheese Quality and Flavor. Front. Microbiol. 2020, 11, 3154. [Google Scholar] [CrossRef]
  32. Leite, A.; Leite, D.C.D.A.; Del Aguila, E.; Alvares, T.; Peixoto, R.; Miguel, M.; Silva, J.; Paschoalin, V. Microbiological and chemical characteristics of Brazilian kefir during fermentation and storage processes. J. Dairy Sci. 2013, 96, 4149–4159. [Google Scholar] [CrossRef]
  33. Zeppa, G.; Conterno, L.; Gerbi, V. Determination of Organic Acids, Sugars, Diacetyl, and Acetoin in Cheese by High-Performance Liquid Chromatography. J. Agric. Food Chem. 2001, 49, 2722–2726. [Google Scholar] [CrossRef] [PubMed]
  34. Kurtzman, C.P.; Fell, J.W.; Boekhout, T.; Robert, V. Methods for isolation, phenotypic characterization and maintenance of yeasts. In The Yeasts: A Taxonomic Study, 5th ed.; Kurtzman, C.P., Fell, J.W., Boekhout, T., Eds.; Elsevier: Amsterdam, The Netherlands, 2011; pp. 87–110. [Google Scholar] [CrossRef]
  35. Canganella, F.; Ovidi, M.; Paganini, S.; Vettraino, A.; Bevilacqua, L.; Trovatelli, L. Survival of undesirable micro-organisms in fruit yoghurts during storage at different temperatures. Food Microbiol. 1998, 15, 71–77. [Google Scholar] [CrossRef]
  36. Moineau-Jean, A.; Guévremont, E.; Champagne, C.P.; Roy, D.; Raymond, Y.; LaPointe, G. Fate of Escherichia coli and Kluyveromyces marxianus contaminants during storage of Greek-style yogurt produced by centrifugation or ultrafiltration. Int. Dairy J. 2017, 72, 36–43. [Google Scholar] [CrossRef]
  37. Belloch, C.; Querol, A.; Barrio, E. Yeast and Moulds|Kluyveromyces spp. In Encyclopedia Dairy Science, 2nd ed.; Fuquay, J.W., Ed.; Academic Press: Cambridge, MA, USA, 2011; pp. 754–764. [Google Scholar] [CrossRef]
  38. Fasoli, G.; Tofalo, R.; Lanciotti, R.; Schirone, M.; Patrignani, F.; Perpetuini, G.; Grazia, L.; Corsetti, A.; Suzzi, G. Chromosome arrangement, differentiation of growth kinetics and volatile molecule profiles in Kluyveromyces marxianus strains from Italian cheeses. Int. J. Food Microbiol. 2015, 214, 151–158. [Google Scholar] [CrossRef] [PubMed]
  39. Bosset, J.; Gauch, R. Comparison of the volatile flavour compounds of six european ‘AOC’ cheeses by using a new dynamic headspace GC-MS method. Int. Dairy J. 1993, 3, 359–377. [Google Scholar] [CrossRef]
  40. Mounchili, A.; Wichtel, J.; Bosset, J.; Dohoo, I.; Imhof, M.; Altieri, D.; Mallia, S.; Stryhn, H. HS-SPME gas chromatographic characterization of volatile compounds in milk tainted with off-flavour. Int. Dairy J. 2005, 15, 1203–1215. [Google Scholar] [CrossRef]
  41. Hazelwood, L.A.; Daran, J.-M.; van Maris, A.J.; Pronk, J.T.; Dickinson, J.R. The Ehrlich Pathway for Fusel Alcohol Production: A Century of Research on Saccharomyces cerevisiae Metabolism. Appl. Environ. Microbiol. 2008, 74, 2259–2266. [Google Scholar] [CrossRef] [Green Version]
  42. Wittmann, C.; Hans, M.; Bluemke, W. Metabolic physiology of aroma-producingKluyveromyces marxianus. Yeast 2002, 19, 1351–1363. [Google Scholar] [CrossRef]
  43. Cadwallader, K.R.; Singh, T.K. Flavours and off-flavours in milk and dairy products. In Advanced Dairy Chemistry, 3rd ed; McSweeney, P., Fox, P.F., Eds.; Springer: New York, NY, USA, 2009; pp. 631–690. [Google Scholar] [CrossRef]
  44. Güler, Z. Changes in salted yoghurt during storage. Int. J. Food Sci. Technol. 2007, 42, 235–245. [Google Scholar] [CrossRef]
  45. Chen, H.; Fink, G.R. Feedback control of morphogenesis in fungi by aromatic alcohols. Genes Dev. 2006, 20, 1150–1161. [Google Scholar] [CrossRef] [Green Version]
  46. Gori, K.; Knudsen, P.B.; Nielsen, K.F.; Arneborg, N.; Jespersen, L. Alcohol-based quorum sensing plays a role in adhesion and sliding motility of the yeast Debaryomyces hansenii. FEMS Yeast Res. 2011, 11, 643–652. [Google Scholar] [CrossRef] [Green Version]
  47. Balbino, T.R.; da Silveira, F.A.; Ventorim, R.Z.; Nascimento, A.G.D.; de Oliveira, L.L.; da Silveira, W.B. Adaptive responses of Kluyveromyces marxianus CCT 7735 to 2-phenylethanol stress: Alterations in membrane fatty-acid composition, ergosterol content, exopolysaccharide production and reduction in reactive oxygen species. Fungal Genet. Biol. 2021, 151, 103561. [Google Scholar] [CrossRef]
  48. Kern, K.; Nunn, C.D.; Pichová, A.; Dickinson, J.R. Isoamyl alcohol-induced morphological change in Saccharomyces cerevisiae involves increases in mitochondria and cell wall chitin content. FEMS Yeast Res. 2004, 5, 43–49. [Google Scholar] [CrossRef] [Green Version]
  49. Chen, C.; Zhou, W.; Yu, H.; Yuan, J.; Tian, H. Evaluation of the Perceptual Interactions among Aldehydes in a Cheddar Cheese Matrix According to Odor Threshold and Aroma Intensity. Molecules 2020, 25, 4308. [Google Scholar] [CrossRef]
  50. Hickey, C.; O’Sullivan, M.; Davis, J.; Scholz, D.; Kilcawley, K.; Wilkinson, M.; Sheehan, J. The effect of buttermilk or buttermilk powder addition on functionality, textural, sensory and volatile characteristics of Cheddar-style cheese. Food Res. Int. 2018, 103, 468–477. [Google Scholar] [CrossRef] [Green Version]
  51. McSweeney, P.L.H.; Sousa, M.J. Biochemical pathways for the production of flavour compounds in cheeses during ripening: A review. Lait 2000, 80, 293–324. [Google Scholar] [CrossRef]
  52. Carrillo-Carrión, C.; Cárdenas, S.; Valcárcel, M. Vanguard/rearguard strategy for the evaluation of the degradation of yoghurt samples based on the direct analysis of the volatiles profile through headspace-gas chromatography–mass spectrometry. J. Chromatogr. A 2007, 1141, 98–105. [Google Scholar] [CrossRef]
  53. Sfakianakis, P.; Tzia, C. Flavour profiling by gas chromatography–mass spectrometry and sensory analysis of yoghurt derived from ultrasonicated and homogenised milk. Int. Dairy J. 2017, 75, 120–128. [Google Scholar] [CrossRef]
  54. Hanková, M.; Čížková, H. Sensory active substances causing off-odour in liquid whey during storage. Potravinarstvo Slovak J. Food Sci. 2020, 14, 729–734. [Google Scholar] [CrossRef]
  55. Morrissey, J.P.; Etschmann, M.M.W.; 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 2014, 32, 3–16. [Google Scholar] [CrossRef]
  56. Hoffmann, A.; Kupsch, C.; Walther, T.; Löser, C. Synthesis of ethyl acetate from glucose by Kluyveromyces marxianus, Cyberlindnera jadinii and Wickerhamomyces anomalus depending on the induction mode. Eng. Life Sci. 2020, 21, 154–168. [Google Scholar] [CrossRef]
  57. Saerens, S.M.G.; Delvaux, F.R.; Verstrepen, K.J.; Thevelein, J.M. Production and biological function of volatile esters in Saccharomyces cerevisiae. Microb. Biotechnol. 2009, 3, 165–177. [Google Scholar] [CrossRef] [Green Version]
  58. Moineau-Jean, A.; Raymond, Y.; Sabik, H.; Graveline, N.; Champagne, C.P.; Roy, D.; LaPointe, G. Effect of manufacturing processes and storage on aroma compounds and sensory properties of yoghurt. Int. Dairy J. 2020, 105, 104662. [Google Scholar] [CrossRef]
  59. Padilla, B.; Manzanares, P.; Belloch, C. Yeast species and genetic heterogeneity within Debaryomyces hansenii along the ripening process of traditional ewes’ and goats’ cheeses. Food Microbiol. 2014, 38, 160–166. [Google Scholar] [CrossRef] [Green Version]
  60. Huang, Z.; Huang, L.; Xing, G.; Xu, X.; Tu, C.; Dong, M. Effect of Co-Fermentation with Lactic Acid Bacteria and K. marxianus on Physicochemical and Sensory Properties of Goat Milk. Foods 2020, 9, 299. [Google Scholar] [CrossRef] [Green Version]
  61. Leclercq-Perlat, M.-N.; Corrieu, G.; Spinnler, H.-E. Comparison of Volatile Compounds Produced in Model Cheese Medium Deacidified by Debaryomyces hansenii or Kluyveromyces marxianus. J. Dairy Sci. 2004, 87, 1545–1550. [Google Scholar] [CrossRef] [Green Version]
  62. Benozzi, E.; Romano, A.; Capozzi, V.; Makhoul, S.; Cappellin, L.; Khomenko, I.; Aprea, E.; Scampicchio, M.; Spano, G.; Märk, T.D.; et al. Monitoring of lactic fermentation driven by different starter cultures via direct injection mass spectrometric analysis of flavour-related volatile compounds. Food Res. Int. 2015, 76, 682–688. [Google Scholar] [CrossRef]
  63. Zhao, L.; Feng, R.; Ren, F.; Mao, X. Addition of buttermilk improves the flavor and volatile compound profiles of low-fat yogurt. LWT 2018, 98, 9–17. [Google Scholar] [CrossRef]
  64. Cárcoba, R.; Delgado, T.; Rodríguez, A. Comparative performance of a mixed strain starter in cow’s milk, ewe’s milk and mixtures of these milks. Eur. Food Res. Technol. 2000, 211, 141–146. [Google Scholar] [CrossRef]
  65. Dzialo, M.C.; Park, R.; Steensels, J.; Lievens, B.; Verstrepen, K.J. Physiology, ecology and industrial applications of aroma formation in yeast. FEMS Microbiol. Rev. 2017, 41, S95–S128. [Google Scholar] [CrossRef] [Green Version]
  66. Tamime, A.Y.; Robinson, R.K. Biochemistry of fermentation. In Tamime and Robinson’s Yoghurt, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2007; pp. 535–607. [Google Scholar] [CrossRef]
  67. Aunsbjerg, S.D.; Honoré, A.; Marcussen, J.; Ebrahimi, P.; Vogensen, F.; Benfeldt, C.; Skov, T.; Knøchel, S. Contribution of volatiles to the antifungal effect of Lactobacillus paracasei in defined medium and yogurt. Int. J. Food Microbiol. 2015, 194, 46–53. [Google Scholar] [CrossRef] [PubMed]
  68. Mende, S.; Rohm, H.; Jaros, D. Influence of exopolysaccharides on the structure, texture, stability and sensory properties of yoghurt and related products. Int. Dairy J. 2016, 52, 57–71. [Google Scholar] [CrossRef]
  69. Tian, H.; Xu, X.; Sun, X.; Chen, C.; Yu, H. Evaluation of the perceptual interaction among key aroma compounds in milk fan by gas chromatography−olfactometry, odor threshold, and sensory analyses. J. Dairy Sci. 2020, 103, 5863–5873. [Google Scholar] [CrossRef] [PubMed]
  70. Serra, M.; Trujillo, A.J.; Guamis, B.; Ferragut, V. Flavour profiles and survival of starter cultures of yoghurt produced from high-pressure homogenized milk. Int. Dairy J. 2009, 19, 100–106. [Google Scholar] [CrossRef]
  71. Adhikari, K.; Grün, I.; Mustapha, A.; Fernando, L. Changes in the Profile of Organic Acids in Plain Set and Stirred Yogurts during Manufacture and Refrigerated Storage. J. Food Qual. 2002, 25, 435–451. [Google Scholar] [CrossRef]
  72. Nuryana, I.; Andriani, A.; Lisdiyanti, P. Yopi Analysis of organic acids produced by lactic acid bacteria. IOP Conf. Series: Earth Environ. Sci. 2019, 251, 012054. [Google Scholar] [CrossRef]
  73. Gänzle, M. Lactic metabolism revisited: Metabolism of lactic acid bacteria in food fermentations and food spoilage. Curr. Opin. Food Sci. 2015, 2, 106–117. [Google Scholar] [CrossRef]
  74. Lachance, M.A. Kluyveromyces van der Walt (1971). In The Yeasts, 5th ed.; Kurtzman, C.P., Fell, J.W., Boekhout, T., Eds.; Elsevier B.V.: Amsterdam, The Netherlands, 2011; pp. 471–481. [Google Scholar] [CrossRef]
  75. Tofalo, R.; Fasoli, G.; Schirone, M.; Perpetuini, G.; Pepe, A.; Corsetti, A.; Suzzi, G. The predominance, biodiversity and biotechnological properties of Kluyveromyces marxianus in the production of Pecorino di Farindola cheese. Int. J. Food Microbiol. 2014, 187, 41–49. [Google Scholar] [CrossRef]
  76. Fernández-García, E.; McGregor, J.U. Determination of Organic Acids During the Fermentation and Cold Storage of Yogurt. J. Dairy Sci. 1994, 77, 2934–2939. [Google Scholar] [CrossRef]
  77. Vénica, C.I.; Perotti, M.C.; Bergamini, C.V. Organic acids profiles in lactose-hydrolyzed yogurt with different matrix composition. Dairy Sci. Technol. 2014, 94, 561–580. [Google Scholar] [CrossRef] [Green Version]
  78. La Torre, L.; Tamime, A.Y.; Muir, D.D. Rheology and sensory profiling of set-type fermented milks made with different commercial probiotic and yoghurt starter cultures. Int. J. Dairy Technol. 2003, 56, 163–170. [Google Scholar] [CrossRef]
Figure 1. Viable counts (log CFU·g−1) of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 grown in skyr at 8 °C for 21 days. Mean values and SD from triplicate experiments are presented.
Figure 1. Viable counts (log CFU·g−1) of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 grown in skyr at 8 °C for 21 days. Mean values and SD from triplicate experiments are presented.
Foods 11 01776 g001
Figure 2. Principal component analysis (PCA) of volatile organic compound (VOC) dataset presented by a score plot of the first two PCs. VOCs were quantified during growth of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr at 8 °C for 21 days, and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Mean values from triplicate experiments were used as input in PCA.
Figure 2. Principal component analysis (PCA) of volatile organic compound (VOC) dataset presented by a score plot of the first two PCs. VOCs were quantified during growth of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr at 8 °C for 21 days, and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Mean values from triplicate experiments were used as input in PCA.
Foods 11 01776 g002
Figure 3. Distribution of the major types of volatile compounds (alcohols, aldehydes, esters and ketones) produced by yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 grown in skyr at 8 °C for 21 days and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Mean values from triplicate experiments are presented.
Figure 3. Distribution of the major types of volatile compounds (alcohols, aldehydes, esters and ketones) produced by yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 grown in skyr at 8 °C for 21 days and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Mean values from triplicate experiments are presented.
Foods 11 01776 g003
Figure 4. Hierarchical cluster analysis (Pearson’s correlation) and a heat map of VOC profiles shown in relation to yeast strains and incubation time. VOCs were determined during growth of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr at 8 °C for 21 days and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Heat map was constructed using mean values from triplicates.
Figure 4. Hierarchical cluster analysis (Pearson’s correlation) and a heat map of VOC profiles shown in relation to yeast strains and incubation time. VOCs were determined during growth of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr at 8 °C for 21 days and in controls (skyr without yeasts). Sample codes denote the strain ID and the incubation period (0, 7, 14 and 21 days). Heat map was constructed using mean values from triplicates.
Foods 11 01776 g004
Figure 5. Concentration of lactose (A) and galactose (B) in skyr after 7 days (grey columns), 14 days (white columns) and 21 days (black columns) incubation of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 at 8 °C, and in controls (skyr without yeasts). The dashed lines show initial concentration of lactose and galactose in skyr (mean values at Day 0). Mean values and SD (bars) from triplicate experiments are presented. Different subscripts denote statistical differences (p < 0.05) between the yeast strains within the same period of incubation determined by the one-way ANOVA (Tukey’s post-hoc analysis). Asterisks indicate values significantly different from Day 0 (Student’s t-test, p < 0.05).
Figure 5. Concentration of lactose (A) and galactose (B) in skyr after 7 days (grey columns), 14 days (white columns) and 21 days (black columns) incubation of yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 at 8 °C, and in controls (skyr without yeasts). The dashed lines show initial concentration of lactose and galactose in skyr (mean values at Day 0). Mean values and SD (bars) from triplicate experiments are presented. Different subscripts denote statistical differences (p < 0.05) between the yeast strains within the same period of incubation determined by the one-way ANOVA (Tukey’s post-hoc analysis). Asterisks indicate values significantly different from Day 0 (Student’s t-test, p < 0.05).
Foods 11 01776 g005
Figure 6. Production of carbon dioxide (CO2, nmol·g−1) by yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr, and in controls (without addition of yeasts) after 7 days (grey columns), 14 days (white columns) and 21 days (black columns) incubation at 8 °C. Mean values and SD (bars) from triplicate experiments are presented. Hashes denote significant differences between yeasts and controls within the same day of incubation, and asterisks indicate significant differences across the time points determined by the Student’s t-test (p < 0.05).
Figure 6. Production of carbon dioxide (CO2, nmol·g−1) by yeasts K. marxianus (strains Km1, Km2 and Km3), P. kudriavzevii Pk1 and T. delbrueckii Td1 in skyr, and in controls (without addition of yeasts) after 7 days (grey columns), 14 days (white columns) and 21 days (black columns) incubation at 8 °C. Mean values and SD (bars) from triplicate experiments are presented. Hashes denote significant differences between yeasts and controls within the same day of incubation, and asterisks indicate significant differences across the time points determined by the Student’s t-test (p < 0.05).
Foods 11 01776 g006
Table 1. Yeasts in this study identified by 26S rRNA gene sequencing.
Table 1. Yeasts in this study identified by 26S rRNA gene sequencing.
SpeciesStrain IDIdentity, %GeneBank Accession nr.Origin
Kluyveromyces marxianusKm1100OM854769This study
Kluyveromyces marxianusKm2100OM854770This study
Kluyveromyces marxianusKm399.8OM854771This study
Pichia kudriavzeviiPk1100OM854772Arla Foods, Denmark
Torulaspora delbrueckiiTd1100OM854773Arla Foods, Denmark
Table 2. Concentration of volatile compounds in skyr at start (VOCt0) and after 21 days (VOCt21) incubation of yeasts Kluyveromyces marxianus Km1 Km2 and Km3, Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1, and in skyr controls (without added yeasts), VOC retention indices (RI) and threshold values.
Table 2. Concentration of volatile compounds in skyr at start (VOCt0) and after 21 days (VOCt21) incubation of yeasts Kluyveromyces marxianus Km1 Km2 and Km3, Pichia kudriavzevii Pk1 and Torulaspora delbrueckii Td1, and in skyr controls (without added yeasts), VOC retention indices (RI) and threshold values.
CompoundRI
Calculated
RI 1
Reference
VOCt0, µg·kg−1 3VOCt21, µg·kg−1 2Odour Threshold
µg·kg−1 3
Odour Description 4
ControlK. marxianus
Km1
K. marxianus
Km2
K. marxianus
Km3
P. kudriavzevii Pk1T. delbrueckii Td1
Alcohols
1-Butanol1166116637 ± 6.939 ± 1.4 A22 ± 1.9 B23 ± 1.9 B23 ± 5.1 B14 ± 4.3 B19 ± 2.7 B500-2500rancid, sweet
2-Butanol1031988–10532.1 ± 0.51.6 ± 0.05 B19 ± 3.2 A21 ± 0.62 A17 ± 3.4 A28 ± 7.9 A1.6 ± 0.17 B3300vinous, fruity
1-Hexanol135513723.8 ± 1.04.5 ± 1.4 C56 ± 3.1 A46 ± 2.4 B48 ± 2.1 B4.8 ± 0.53 C5.9 ± 1.3 C2500sweet alcohol
3-Methyl-1-butanol12271174–125521 ± 4.934 ± 13 C1780 ± 48 A1918 ± 100 A1781 ± 163 A412 ± 63 B30 ± 5 C3060whiskey, malt, burnt
3-Methyl-3-buten-1-ol12621226–127710.1 ± 2.432 ± 1.826 ± 1.227 ± 1.626 ± 3.425 ± 2.428 ± 2.6250sweet fruity
2-Methyl-1-propanol111011100.61 ± 0.251.4 ± 0.62 B446 ± 58 A562 ± 17 A591 ± 102 A38 ± 8.9 AB0.90 ± 0.15 B40,000sweet musty, fusel
1-Pentanol126712736.5 ± 1.68.1 ± 0.68 C30 ± 0.76 A24 ± 2.1 B23 ± 2.4 B9.2 ± 0.85 C9.9 ± 1.7 C120fusel
1-Penten-3-ol117911770.28 ± 0.120.40 ± 0.15 B2.2 ± 0.15 A1.4 ± 0.19 AB1.9 ± 0.36 A0.99 ± 0.25 AB2.3 ± 0.83 A500pungent, fruity
2-Phenylethanol192319300.36 ± 0.161.1 ± 1.04 BC3.2 ± 0.20 AB4.7 ± 0.41 A4.5 ± 0.80 A1.04 ± 0.71 BC0.36 ± 0.41 C750–1100mild rose
1-Propanol104310410.35 ± 0.080.35 ± 0.05 B14 ± 3.4 A15 ± 2.3 A14 ± 2.02 A0.69 ± 0.01 B0.34 ± 0.06 B9000alcoholic
2-Propanol936884–9631.1 ± 0.200.65 ± 0.41 B2.7 ± 0.81 AB3.3 ± 1.3 AB3.5 ± 0.53 AB6.9 ± 4.01 A0.68 ± 0.06 B40,000-78,000rubbing alcohol
1-Octanol155915705.4 ± 1.75.8 ± 2.4 C15 ± 0.80 B20 ± 1.7 AB21 ± 1.4 A4.8 ± 0.44 C7.0 ± 1.9 C110–130aromatic, orange
Aldehydes
Butanal875830–9112.9 ± 0.920.82 ± 0.170.59 ± 0.130.45 ± 0.020.60 ± 0.150.33 ± 0.150.70 ± 0.177–37pungent
Heptanal1191118912 ± 3.31.1 ± 0.14 B3.9 ± 1.1 A2.5 ± 0.22 AB2.3 ± 0.63 AB0.76 ± 0.10 B1.1 ± 0.34 B3–60fatty, pungent, fruity
Hexanal1079108221 ± 5.71.3 ± 1.1 B8.6 ± 0.55 A4.2 ± 0.22 AB3.8 ± 0.79 B0.64 ± 0.07 B2.7 ± 1.9 B5–9sharp, fruity, grass
2-Methylbutanal9129124.9 ± 1.13.3 ± 0.07 B10.5 ± 0.72 A11 ± 1.7 A14 ± 3.5 A0.35 ± 0.33 B3.6 ± 0.34 B1–6apple, malt
3-Methylbutanal9169167.6 ± 1.92.9 ± 1.0 B9.3 ± 0.45 A10.2 ± 0.75 A9.8 ± 1.8 A12 ± 2.3 A4.03 ± 0.24 B1–3apple, malt
3-Methyl-2-butenal11991189–12366.8 ± 1.52.5 ± 0.45 A0.46 ± 0.01 C0.66 ± 0.07 C0.67 ± 0.15 C2.1 ± 0.30 AB1.3 ± 0.42 BCnrmetallic, herbaceous
2-Methylpropanal8118127.1 ± 1.15.3 ± 0.53 D10.2 ± 1.6 BC12 ± 0.32 B19 ± 2.2 A5.4 ± 0.90 D6.2 ± 0.70 CD1.5sharp, pungent
Trans-2-Decenal166016624.9 ± 1.72.9 ± 2.3 AB5.3 ± 0.33 AB6.9 ± 0.55 A6.8 ± 1.3 A1.9 ± 0.28 B3.2 ± 1.3 AB0.4orange, tallow
Esters
Butyl acetate107010820.019 ± 0.0360.51 ± 0.10 AB1.8 ± 0.71 AB2.8 ± 1.6 A0.92 ± 0.15 AB1.01 ± 0.60 AB0.094 ± 0.023 B10–100sweet, fruity
Ethyl acetate891891ndnd B289 ± 4.7 A283 ± 15 A291 ± 0.8 A278 ± 15 AndB5–5000pineapple, solvent
Ethyl butyrate103310360.51 ± 0.150.87 ± 0.23 D81 ± 4.2 A44 ± 2.2 B34 ± 2.3 C1.63 ± 0.48 D0.57 ± 0.06 D1pineapple
Ethyl hexanoate124412510.093 ± 0.0390.072 ± 0.077 C23 ± 2.1 B32 ± 1.5 A19 ± 1.2 Bnd C0.18 ± 0.06 C1fruity, aniseed, apple
Ethyl octanoate141914450.55 ± 0.240.26 ± 0.37 C10.1 ± 1.2 B23 ± 5.7 A14 ± 2.1 Bnd C0.11 ± 0.15 C15fruity, sour apple
Ethyl formate8238230.27 ± 0.170.038 ± 0.054 B5.7 ± 2.9 A7.5 ± 0.36 A8.5 ± 1.5 A0.39 ± 0.08 B0.27 ± 0.04 B150,000fruity, pungent
Ethyl lactate134413560.20 ± 0.071.3 ± 0.89 C21 ± 1.3 A17 ± 1.7 AB14 ± 2.4 B0.58 ± 0.26 C0.82 ± 0.68 C14,000buttery, fruity
Ethyl propionate956956nd0.076 ± 0.069 B6.2 ± 2.6 A7.6 ± 4.9 AB1.3 ± 0.15 B2.1 ± 1.1 ABnd B10fruity, rum, pineapple
Hexyl acetate128112910.13 ± 0.040.26 ± 0.04 B1.3 ± 0.23 A1.7 ± 0.30 A1.2 ± 0.12 A0.21 ± 0.01 B0.24 ± 0.04 B2sweet, fruity, pear
2-Methylpropyl acetate101310180.011 ± 0.0071.78 ± 0.333.6 ± 1.35.1 ± 3.52.1 ± 0.120.87 ± 1.010.33 ± 0.1466fruity, floral
3-Methylbutyl acetate113611420.17 ± 0.0910.30 ± 0.31 C33 ± 8.1 AB57 ± 21 A29 ± 0.49 AB3.8 ± 2.02 BC0.31 ± 0.31 C5sweet, fruity, banana
4-Pentenyl acetate119912040.078 ± 0.0360.08 ± 0.02 B2.08 ± 1.04 AB3.6 ± 2.4 A0.92 ± 0.36 AB0.13 ± 0.11 B0.01 ± 0.01 Bnrgreen, plastic, weedy
Pentyl acetate118311850.22 ± 0.110.58 ± 0.18 AB0.85 ± 0.24 AB0.94 ± 0.31 A0.59 ± 0.01 AB0.36 ± 0.12 AB0.23 ± 0.04 B80banana
2-Phenylethyl acetate182118350.074 ± 0.0410.10 ± 0.13 C2.3 ± 0.44 B5.03 ± 0.28 A2.9 ± 0.47 B0.04 ± 0.05 C0.10 ± 0.07 C650sweet, floral, fruity
Ketones
2,3-Butanedione984985178 ± 36526 ± 11 A48 ± 4.3 C49 ± 10 C48 ± 5.4 C72 ± 9.3 C251 ± 47 B1.1–6.5chlorine, buttery
2-Butanone905906191 ± 25136 ± 5.6 A40 ± 7.9 B44 ± 2.9 B62 ± 11 B47 ± 11 B142 ± 10.5 A8400ethereal, fruity
2-Heptanone1189118960 ± 946 ± 4.2 B47 ± 1.4 B41 ± 1.9 B40.2 ± 2.7 B46 ± 0.22 B60.1 ± 2.6 A140–3000fruity, banana, spicy
3-Hydroxy-2-butanone12901302152 ± 28201 ± 10 A55 ± 0.92 B56 ± 15 B79 ± 15 B128 ± 30.3 AB173 ± 48 A800buttery, woody, yogurt
2-Nonanone1380139831 ± 100.41 ± 0.58 B12 ± 0.80 AB20.5 ± 1.5 A18 ± 2.5 A2.5 ± 0.33 B24 ± 10.6 A5-200hot milk, soap, green
2-Propanone815815109 ± 1283 ± 8.2 A8.3 ± 2.9 B10.4 ± 1.4 B19 ± 1.5 B38 ± 28 B86 ± 8.7 A20,000fruity
1 Reference RI are from analysis of authentic reference compounds or reported intervals of RI (https://pubchem.ncbi.nlm.nih.gov). 2 Relative abundance of each compound was calculated from ratio of the peak area of the compound to the internal standard (4-Methyl-1-pentanol). Different superscript within a raw denote statistically different values (p < 0.05) determined by the one-way ANOVA (Tukey’s post-hoc analysis). 3 Odour thresholds from https://pubchem.ncbi.nlm.nih.gov, or and http://www.leffingwell.com/odorthre.htm; nr–not recorded. 4 Odour description from https://pubchem.ncbi.nlm.nih.gov, or http://www.thegoodscentscompany.com/.
Table 3. Concentration of organic acids in skyr during yeast incubation.
Table 3. Concentration of organic acids in skyr during yeast incubation.
Organic AcidsInitial Concentration, mg·100 g−1Yeast SpeciesConcentration, mg·100 g−1 1
Day 7Day 14Day 21
Acetic acid54 ± 1.4Control57 ± 2.957 ± 2.267 ± 4.2 A *
K. marxianus Km154 ± 3.252 ± 1.148 ± 2.4 C *
K. marxianus Km261 ± 4.756 ± 0.67 *50.4 ± 2.5 BC
K. marxianus Km356 ± 4.854 ± 1.0958 ± 7.2 ABC
P. kudriavzevii Pk158 ± 2.454 ± 3.460.2 ± 2.9 AB
T. delbrueckii Td155 ± 6.355 ± 0.9458 ± 4.1 ABC
Citric acid120 ± 9.1Control94 ± 3.4 A *98 ± 13 A96 ± 5.7 A *
K. marxianus Km186 ± 1.8 AB *74 ± 1.1 B *48 ± 8.1 B *
K. marxianus Km281 ± 6.9 AB *74 ± 1.6 B *80 ± 8.4 A *
K. marxianus Km382 ± 7.7 AB *74 ± 2.7 B *83 ± 7.9 A *
P. kudriavzevii Pk181 ± 7.3 AB *81 ± 7.3 AB *75 ± 14 A *
T. delbrueckii Td173 ± 3.9 B *71 ± 3.1 B *79 ± 0.64 A *
Lactic acid224 ± 8.5Control244 ± 13 A267 ± 1.4 A *287 ± 31 *
K. marxianus Km1223 ± 4.8 AB233 ± 3.1 B257 ± 5.3 *
K. marxianus Km2218 ± 5.2 B245 ± 18 AB261 ± 6.4 *
K. marxianus Km3219 ± 7.9 B236 ± 0.57 B285 ± 19 *
P. kudriavzevii Pk1223 ± 6.1 AB241 ± 9.7 B283 ± 12 *
T. delbrueckii Td1221 ± 9.2 B231 ± 1.8 B271 ± 16 *
Orotic acid5.1 ± 0.16Control4.7 ± 0.384.6 ± 0.09 A *4.8 ± 0.57 A
K. marxianus Km14.3 ± 0.09 *4.1 ± 0.03 B *3.3 ± 0.14 B *
K. marxianus Km24.2 ± 0.17 *5.1 ± 0.29 A4.4 ± 0.17 A *
K. marxianus Km34.1 ± 0.19 *4.1 ± 0.08 B *4.8 ± 0.41 A
P. kudriavzevii Pk14.2 ± 0.17 *4.2 ± 0.23 B *5.01 ± 0.29 A
T. delbrueckii Td14.02 ± 0.39 *4.1 ± 0.06 B *4.6 ± 0.34 A
Uric acid3.7 ± 0.36Control3.1 ± 0.913.5 ± 0.20 A3.7 ± 0.73
K. marxianus Km12.7 ± 0.382.5 ± 0.20 B *2.6 ± 0.19 *
K. marxianus Km23.1 ± 0.243.4 ± 0.20 A2.8 ± 0.58
K. marxianus Km33.03 ± 0.293.2 ± 0.40 A3.2 ± 0.74
P. kudriavzevii Pk13.02 ± 0.403.2 ± 0.18 A *3.2 ± 0.60
T. delbrueckii Td12.9 ± 0.45 *3.01 ± 0.21 AB3.5 ± 0.40
1 Different subscripts within a column denote values statistically different (p < 0.05) between the treatments, as determined by one-way ANOVA and Tukey’s post-hoc analysis. Asterisks denote statistical differences (p < 0.05) compared to the initial values (Student’s t-test).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Srimahaeak, T.; Petersen, M.A.; Lillevang, S.K.; Jespersen, L.; Larsen, N. Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr. Foods 2022, 11, 1776. https://doi.org/10.3390/foods11121776

AMA Style

Srimahaeak T, Petersen MA, Lillevang SK, Jespersen L, Larsen N. Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr. Foods. 2022; 11(12):1776. https://doi.org/10.3390/foods11121776

Chicago/Turabian Style

Srimahaeak, Thanyaporn, Mikael Agerlin Petersen, Søren K. Lillevang, Lene Jespersen, and Nadja Larsen. 2022. "Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr" Foods 11, no. 12: 1776. https://doi.org/10.3390/foods11121776

APA Style

Srimahaeak, T., Petersen, M. A., Lillevang, S. K., Jespersen, L., & Larsen, N. (2022). Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr. Foods, 11(12), 1776. https://doi.org/10.3390/foods11121776

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

Article Metrics

Back to TopTop