If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, CanadaLaboratoire de mycologie alimentaire (LMA), Université Laval, Québec, G1V 0A6, Canada
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, CanadaLaboratoire de mycologie alimentaire (LMA), Université Laval, Québec, G1V 0A6, Canada
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, CanadaLaboratoire de mycologie alimentaire (LMA), Université Laval, Québec, G1V 0A6, Canada
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, Canada
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, CanadaAgriculture and Agri-Food Canada, Saint-Hyacinthe Research and Development Centre, Saint-Hyacinthe, J2S 8E3, Canada
Department of Food Sciences, Institute of Nutrition and Functional Foods (INAF), STELA Dairy Research Centre, Université Laval, Québec, G1V 0A6, CanadaLaboratoire de mycologie alimentaire (LMA), Université Laval, Québec, G1V 0A6, Canada
Indigenous microorganisms are important components of the complex ecosystem of many dairy foods including cheeses, and they are potential contributors to the development of a specific cheese's sensory properties. Among these indigenous microorganisms are the yeasts Cyberlindnera jadinii, Pichia kudriavzevii, and Kazachstania servazzii, which were previously detected using traditional microbiological methods in both raw milk and some artisanal specialty cheeses produced in the province of Québec, Canada. However, their levels across different cheese varieties are unknown. A highly specific and sensitive real-time quantitative PCR assay was developed to quantitate these yeast species in a variety of specialty cheeses (bloomy-rind, washed-rind, and natural-rind cheeses from raw, thermized, and pasteurized milks). The specificity of the quantitative PCR assay was validated, and it showed no cross-amplification with 11 other fungal microorganisms usually found in bloomy-rind and washed-rind cheeses. Cyberlindnera jadinii and P. kudriavzevii were found in the majority of the cheeses analyzed (25 of 29 and 24 of 29 cheeses, respectively) in concentrations up to 104 to 108 gene copies/g in the cheese cores, which are considered oxygen-poor environments, and 101 to 104 gene copies/cm2 in the rind. However, their high abundance was not observed in the same samples. Whereas C. jadinii was present and dominant in all core and rind samples, P. kudriavzevii was mostly present in cheese cores. In contrast, K. servazzii was present in the rinds of only 2 cheeses, in concentrations ranging from 101 to 103 gene copies/cm2, and in 1 cheese core at 105 gene copies/g. Thus, in the ecosystems of specialty cheeses, indigenous yeasts are highly frequent but variable, with certain species selectively present in specific varieties. These results shed light on some indigenous yeasts that establish during the ripening of specialty cheeses.
A consistent cheese quality and flavor can be difficult to achieve from batch to batch, potentially due to the complex microbial ecosystem and the presence of indigenous microorganisms (
). Recent studies have demonstrated that the cheese ecosystem is composed of hundreds of adventitious microbial species that emerge during the ripening process (
Yeasts are natural members of the cheese ecosystem because of their capacity to grow in wide ranges of pH and temperature and their tolerance to high salt concentrations and low water activity (
). Because of their common presence, multiple techniques are used to detect, identify, and quantify yeasts in cheese, including traditional microbiology techniques, quantitative (q)PCR, and omics methods (
Metatranscriptome analysis of fungal strains Penicillium camemberti and Geotrichum candidum reveal cheese matrix breakdown and potential development of sensory properties of ripened Camembert-type cheese.
). To date, few studies have addressed the frequency and distribution of indigenous yeasts in cheeses, despite their potential effects on the physicochemical composition of the cheeses and their microbial ecosystem (
The presence of indigenous yeasts Cyberlindnera jadinii (synonyms: Pichia jadinii, Saccharomyces jadinii; anamorph: Candida utiliz), Pichia kudriavzevii (synonym: Issatchenkia orientalis; anamorph: Candida krusei), and Kazachstania servazzii (synonym: Saccharomyces servazzii) has been previously established in various fermented foods. Kazachstania servazzii has been found in pickles and kefir (
Production of traditional sorghum beer “Ikigage” using Saccharomyces cerevisae, Lactobacillus fermentum and Issatckenkia orientalis as starter cultures.
). In dairy foods, K. servazzii is found less frequently than C. jadinii and P. kudriavzevii, and fewer studies have been conducted to evaluate its role in the dairy environment (
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
). In cheese slurries, C. jadinii, present in Münster cheese, imparts a distinct odor in addition to an elevated pH, and P. kudriavzevii influences the chemical environment of cheese by increasing pH and proteolysis (
In raw milk or artisanal specialty cheeses from the province of Québec, P. kudriavzevii was isolated and identified, using ITS1–5.8S–ITS2 rDNA sequencing, in numerous cheese samples, whereas C. jadinii and K. servazzii were found less frequently (
). Because the approach was not quantitative, the abundance of these yeasts across different cheeses and their ability to influence the cheese ripening process remain unclear.
The aim of this study was to develop a qPCR method to determine the concentrations of C. jadinii, P. kudriavzevii, and K. servazzii in different areas (core and rind) of several types of commercial specialty cheeses.
MATERIALS AND METHODS
Biological Material and Genomic DNA Extraction
Reference fungal strains from the Laboratoire de Mycologie Alimentaire (LMA, Université Laval, Québec, Canada) collection were used to develop a qPCR assay (Table 1). Strains were cultured from culture stock stored at −80°C in YEG medium (10 g/L yeast extract, Bio Basic; 10 g/L glucose, EMD Chemicals) containing 15% glycerol, and were plated on YEG agar (15 g/L bacto agar; BD Diagnostics). Individual colonies were cultured in YEG liquid medium and incubated at 25°C for 48 h. Yeasts were pelleted by centrifugation (5 min at 20,000 × g at room temperature), which were used for genomic DNA extraction using a phenol-chloroform extraction protocol (
Metatranscriptome analysis of fungal strains Penicillium camemberti and Geotrichum candidum reveal cheese matrix breakdown and potential development of sensory properties of ripened Camembert-type cheese.
Triplicates of each of 29 ready-to-eat specialty cheese varieties from 22 cheese plants (for a total of 87 individual cheeses sampled) from the province of Québec, Canada, were sampled. They were classified according to their variety: washed rind (n = 15), bloomy rind (n = 12) or natural rind (n = 2), and by heat treatment of the milk: raw (n = 6), thermized (n = 14), or pasteurized (n = 9), and were coded for confidentiality (Table 2). Thermization was conducted by each cheesemaker according to the regulation of the province of Québec; that is, between 15 and 20 s at 57 to 68°C. The cheeses were stored at 4°C for 3 d. A total of 50 cm2 was sampled and weighed for each cheese rind; in brief, 25 cm2 was scraped from each flat surface to recover the microbiota and a minimum of the cheese matrix underneath (<1 mm). In addition, 10 g of each cheese variety's core was sampled for 1 of 3 specimens for each of the 29 cheeses. All rind and core samples were stored at −80°C before grinding using a liquid nitrogen–cooled CryoMill (Retsch). Grinding consisted of a precooling cycle at 5 Hz, followed by a 4-min grinding step at 25 Hz.
Table 2Distribution of specialty cheeses (C1 to C29) selected from the province of Québec
Cheese total DNA extraction was performed in triplicate for each cheese sample using a Purelink viral RNA/DNA kit (Invitrogen/Life Technologies) with the following modifications: a total of 30 to 40 mg of each ground cheese sample was homogenized in 200 µL of 0.9% NaCl, before the addition of proteinase K and lysis buffer from the kit. After DNA elution, RNA was removed using a 0.5 mg/mL RNase A treatment for 1 h at 37°C. Yield and quality of the extracted DNA were measured using a Synergy HI microplate reader and Gen5 software (BioTek).
Target Gene Sequencing and Design of Real-Time qPCR Primers and Probes
Specific gene sequences for the targeted yeast species were selected based on their availability in GenBank. The genes coding for malic enzyme (CME1; accession number DQ173437.1), glutathione S-transferase Y-1 (GSTY-1; accession number AB021655.1), and centromere H3 (CenH3; accession number DQ826421.1) were selected as qPCR gene targets and sequenced for C. jadinii (forward primer 5′-CTGAAGGTGACGCCATTGCT-3′ and reverse primer 5′-GTGGACGACCAGTGGTGACA-3′), P. kudriavzevii (forward primer 5′- CTGAAGGTGACGCCATTGCT-3′ and reverse primer 5′-AACAATGTGCTGTAAGGTTACGTA TAG G-3′), and K. servazzii (forward primer 5′-TTCCCCCAGTGACCTTGATG-3′ and reverse primer 5′-TTCGCATGTTCCAGTAAACCAA-3′), respectively. Based on the target gene sequences obtained, qPCR-specific primers and TaqMan probes were designed using Primer Express v2.0 (ThermoFisher Scientific) and analyzed with OligoAnalyzer v3.1 (Integrated DNA Technologies; Table 3). Primer specificity was verified by PCR on several strains of the targeted yeast species (Table 1), and resulting PCR amplicons were sequenced and verified using BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi). Each qPCR design was assessed for cross-amplification against 10 other fungal species commonly found in dairy products (Table 1).
Table 3Primers and TaqMan probes used for quantitative (q)PCR quantification of indigenous yeasts and efficiency of each qPCR assay
qPCR Absolute Quantification Using Standard Curve Assays
The qPCR target genes (CEM1, GSTY-1, and CenH3) were individually PCR amplified, cloned in the pCR-II TOPO vector, and transformed in Escherichia coli TOP10 using a TOPO TA Cloning kit (Invitrogen/Life Technologies). Recombinant plasmid DNA was purified with a Plasmid DNA Maxiprep kit (Qiagen), and then sequenced. To use the recombinant plasmid DNA as a quantification standard, the concentration was measured with a Synergy HI microplate reader and Gen5 software (BioTek). The number of gene copies was calculated using the following formula:
where Amt. = amount, and 660 g/mol is the average mass of 1 bp of double-stranded DNA.
When used as a quantification standard during qPCR assays, the recombinant plasmid DNA was diluted 6- to 10-fold and deposited in triplicate in 96-well plates. These serial dilutions allowed us to determine a limit of detection of 3 gene copies per PCR reaction. Standard curves were generated to determine the qPCR efficiency and dynamic range of each reaction (dynamic range for target gene CME1: 3 to 3 × 108 gene copies per reaction; for GSTY-1 and CenH3: 3 to 3 × 109 gene copies per reaction).
Quantitative PCR was performed on an ABI7500 Fast system apparatus (Applied Biosystems/ThermoFisher Scientific). The qPCR plates were set up using an epMotion 5075 VAC automated pipetting system (Eppendorf). Optimal dilutions of DNA samples were chosen to fit the dynamic range of the standard curves and to ensure the absence of PCR inhibitors. The qPCR reagents were used as described by
. The qPCR program began with an initial denaturation step at 95°C for 20 s, followed by 40 cycles of denaturation 95°C for 3 s and annealing/extension at 60°C for 30 s. The qPCR data were analyzed using 7500 software v2.0.6 (ThermoFisher Scientific).
Statistical Analysis
A one-way ANOVA was used to analyze the effect of the type of cheese on the detection of each indigenous yeast being tested. Because of the non-normality of the data, P-values were adjusted using a permutation distribution. Variety of cheese was considered a fixed factor, and comparisons were made between specific types of cheese (rind type and heat treatment). SAS/STAT software (2015; SAS Institute Inc.) was used for analysis; specifically, the GLM procedure was followed by the MULTTEST procedure.
RESULTS
qPCR Assay Efficiency
Highly specific and sensitive qPCR assays were developed for the quantification of 3 indigenous yeast species (C. jadinii, P. kudriavzevii, and K. servazzii) in samples originating from both the rind and core of specialty cheeses. The specificity of qPCR amplification was confirmed by sequencing the obtained qPCR amplicons and by negative qPCR detection using 10 other fungal species isolated from dairy products. BLAST searches of the CME1, GSTY-1, and CenH3 target genes revealed that these genes were present in only one copy in the genomes of C. jadinii [whole-genome sequencing (WGS) project BAEL01], P. kudriavzevii (WGS project JQFK01), and K. servazzii (unpublished WGS project LMA-647, Laboratoire de Mycologie Alimentaire), respectively. Resulting average qPCR efficiencies were 98.4, 97.8, and 96.6% for C. jadinii, P. kudriavzevii, and K. servazzii, respectively (R2 values ≥0.990; Table 3).
Distribution of C. jadinii in Quebec's Specialty Cheeses
Cyberlindnera jadinii was detected in 25 of 29 cheese cores at concentrations of 105 to 108 gene copies/g (Table 4), and was detected in the rinds of 24 of 29 cheeses at 10° to 104 gene copies/cm2 in all cheese varieties. It was absent in both the core and rind of only one cheese, C1, a washed-rind cheese made with raw milk. The population of C. jadinii in the cheese core was not dependent on the cheese variety (P > 0.05) but was correlated with the milk heat treatment applied (P < 0.01, Figure 1), especially in bloomy-rind cheeses made from raw and thermized milk (P < 0.001). At the cheese surface, high C. jadinii counts were obtained in bloomy-rind cheeses (12/29 cheeses) compared with washed-rind cheeses (15/29 cheeses; P ≤ 0.001, Table 4), although this difference was not observed in core samples taken from the same cheese varieties. Cheeses made from raw and thermized milk had more C. jadinii on their surface than pasteurized milk cheeses (P < 0.001), but we detected no differences between raw and thermized milk cheeses (Figure 2).
Table 4Quantitative PCR of Cyberlindnera jadinii in specialty cheeses
Cheeses with the same letter (A to H) were produced by the same cheese plant. Cheeses that do not have a letter to identify the cheese plant do not share a production site with any other cheeses analyzed in this study. Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
<LOD
C2
4.7 ± 1.9 × 101
1.5 ± 1.3 × 106
C3 (B)
3.6 ± 6.2
5.4 ± 9.4 × 105
Thermized
C4 (A)
<LOD
1.6 ± 1.6 × 106
C5 (C)
8.2 ± 3.4 × 102
1.6 ± 0.2 × 108
C6 (D)
<LOD
1.6 ± 2.8 × 106
C7
0.6 ± 1.1 × 101
<LOD
C8 (D)
2.4 ± 4.2 × 101
1.3 ± 1.8 × 106
Pasteurized
C9 (E)
0.6 ± 1.0 × 101
<LOD
C10 (E)
1.7 ± 0.8 × 102
3.6 ± 1.3 × 106
C11
1.0 ± 1.8 × 101
1.0 ± 0.1 × 107
C12 (F)
<LOD
1.9 ± 1.1 × 106
C13
<LOD
1.5 ± 2.1 × 106
C14
1.0 ± 4.0 × 101
1.5 ± 1.6 × 106
C15
4.1 ± 7.0
2.8 ± 4.9 × 106
Bloomy
Raw
C16 (G)
8.8 ± 3.0
5.6 ± 0.8 × 106
C17 (B)
1.6 ± 0.4 × 104
3.8 ± 1.2 × 107
C18 (G)
1.7 ± 2.9 × 101
1.4 ± 2.3 × 106
Thermized
C19 (C)
9.2 ± 2.3 × 103
1.0 ± 0.1 × 108
C20 (H)
0.9 ± 1.6 × 101
<LOD
Pasteurized
C21
9.4 ± 8.1
3.4 ± 0.7 × 106
C22
1.0 ± 1.7 × 102
2.9 ± 2.0 × 106
C23
6.6 ± 5.5 × 101
1.6 ± 0.4 × 106
C24
0.9 ± 1.6 × 101
2.7 ± 0.4 × 106
C25
3.6 ± 3.6 × 101
5.1 ± 2.6 × 106
C26
1.0 ± 1.7 × 101
2.8 ± 0.3 × 106
C27 (F)
3.7 ± 2.1 × 101
2.4 ± 0.5 × 106
Natural
Thermized
C28
5.6 ± 7.5
2.7 ± 0.7 × 106
C29 (H)
4.9 ± 8.4
3.0 ± 5.2 × 106
1 Cheeses with the same letter (A to H) were produced by the same cheese plant. Cheeses that do not have a letter to identify the cheese plant do not share a production site with any other cheeses analyzed in this study. Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
2 LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
Figure 2Quantitative PCR enumeration (gene copies/cm2) of Cyberlindnera jadinii on the surface of bloomy-rind specialty cheeses. Different letters (a, b) indicate significant differences (P < 0.001) among surface samples. Error bars represent SE.
The presence of C. jadinii in multiple cheeses manufactured in the same cheese plant was verified by analyzing 2 different cheese varieties from each of 8 artisanal cheese plants (plant A: C1 and C4; plant B: C3 and C17; plant C: C5 and C19; plant D: C6 and C8; plant E: C9 and C10; plant F: C12 and C27; plant G: C16 and C18; plant H: C20 and C29, Table 4). Counts of C. jadinii on the surface of a bloomy-rind and a washed-rind cheese (C17 and C3) produced at the same cheese plant were significantly different (104 gene copies/cm2 and 1 gene copy/cm2, respectively; P < 0.001). The same significant difference in C. jadinii counts was observed in the cores of these cheeses: C. jadinii was more abundant in bloomy-rind cheeses than in washed-rind cheeses (P < 0.001; Table 4).
Distribution of P. kudriavzevii in Quebec's Specialty Cheeses
Pichia kudriavzevii was detected on the surface of 6 varieties of cheese at 0.9 ± 1.5 to 1.2 ± 2.1 × 102 gene copies/cm2, and in the core of 24 cheeses from 2.4 ± 4.2 × 105 to 4.15 ± 1.03 × 107 gene copies/g (Table 5). It was not detected in the rind or core of 4 cheeses (C1, C12, C15, and C20). The type of heat treatment applied to the milk used for each cheese did not affect the presence or levels of P. kudriavzevii in the cores of washed-rind cheeses. In the bloomy-rind cheeses, however, a significantly higher count of P. kudriavzevii was detected in raw milk cheeses than in thermized and pasteurized milk cheeses (P < 0.05; Figure 3). This indigenous yeast was detected in the core of 11 of 12 bloomy-rind cheeses in amounts ranging from 105 to 107 gene copies/g. Finally, different cheeses produced in the same plant showed different levels of P. kudriavzevii. For example, a difference was observed between washed-rind and bloomy-rind cheeses C3 and C17 from cheese plant B (105 and 107 gene copies/g, respectively). Only 1 of the 2 cheese varieties from cheese plant H, which produced the bloomy-rind cheese C20 (< limit of detection) and the natural-rind cheese C29 (107 gene copies/g), contained P. kudriavzevii.
Table 5Quantitative PCR of Pichia kudriavzevii in specialty cheeses
Cheese samples in which Pichia kudriavzevii was not detected (in both the rind and core) are not presented. Cheeses with the same letter (A to H) were produced by the same cheese plant. (Cheeses with letters not repeated in this table indicate that another cheese produced by this cheese plant was analyzed but the microorganism was absent.) Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
C3 (B)
<LOD
2.4 ± 4.2 × 105
Thermized
C4 (A)
<LOD
1.08 ± 0.20 × 107
C5 (C)
<LOD
1.3 ± 1.8 × 107
C6 (D)
<LOD
4.3 ± 7.4 × 105
C7
<LOD
4.27 ± 7.40 × 106
C8 (D)
<LOD
0.6 ± 1.1 × 104
Pasteurized
C9 (E)
1.4 ± 2.4 × 101
9.6 ± 15.3 × 106
C10 (E)
<LOD
3.71 ± 0.54 × 107
C11
<LOD
2.26 ± 1.63 × 107
C13
<LOD
1.6 ± 2.6 × 107
C14
<LOD
5.2 ± 8.9 × 105
Bloomy
Raw
C16 (G)
0.9 ± 1.5
1.72 ± 2.09 × 107
C17 (B)
<LOD
4.15 ± 1.03 × 107
C18 (G)
5.1 ± 8.9
1.33 ± 1.15 × 107
Thermized
C19 (C)
<LOD
8.76 ± 1.29 × 106
Pasteurized
C21
<LOD
6.28 ± 4.72 × 106
C22
<LOD
7.7 ± 7.3 × 106
C23
1.2 ± 2.1 × 102
3.3 ± 4.7 × 105
C24
<LOD
1.59 ± 1.19 × 107
C25
<LOD
3.33 ± 3.25 × 106
C26
<LOD
1.4 ± 1.4 × 107
C27 (F)
<LOD
1.3 ± 2.0 × 107
Natural
Thermized
C28
<LOD
1.08 ± 0.18 × 107
C29 (H)
7.1 ± 6.2 × 101
2.41 ± 1.25 × 107
1 Cheese samples in which Pichia kudriavzevii was not detected (in both the rind and core) are not presented. Cheeses with the same letter (A to H) were produced by the same cheese plant. (Cheeses with letters not repeated in this table indicate that another cheese produced by this cheese plant was analyzed but the microorganism was absent.) Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
2 LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
Figure 3Quantitative PCR enumeration (gene copies/g) of Pichia kudriavzevii in the cores of bloomy-rind cheeses. Values were calculated from means of each variety of cheese (all 3 biological replicates of each individual cheese were included). Different letters (a, b) indicate significant differences (P < 0.05) among core samples. Error bars represent SE.
Distribution of K. servazzii in Quebec's Specialty Cheeses
Kazachstania servazzii was detected in the rind of only 2 cheeses: natural-rind cheeses C28 and C29 (both made from thermized milk), at 1.89 ± 0.51 × 103 gene copies/cm2 and 0.6 ± 1.0 × 101 gene copies/cm2, respectively (Table 6). Kazachstania servazzii was detected in only 1 cheese core: cheese C12 at 0.9 ± 1.6 × 105 gene copies/g.
Table 6Quantitative PCR of Kazachstania servazzii in specialty cheeses
Cheese samples in which Kazachstania servazzii was not detected (in both the rind and core) are not presented. Cheeses with the same letter (A to H) were produced by the same cheese plant. Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
0.9 ± 1.6 × 105
Natural
Thermized
C28
1.89 ± 0.51 × 103
<LOD
C29 (H)
0.6 ± 1.0 × 101
<LOD
1 Cheese samples in which Kazachstania servazzii was not detected (in both the rind and core) are not presented. Cheeses with the same letter (A to H) were produced by the same cheese plant. Means and SE were calculated from means of each of 3 biological replicates (3 technical replicates were performed per biological replicate) for the rinds of each cheese variety, and from the means of values obtained from 3 technical replicates for each core (1 biological replicate was sampled per cheese variety).
2 LOD = limit of detection, determined using serial dilutions, and set to 3 gene copies per PCR reaction.
This study aimed to establish the frequency and distribution of the 3 indigenous yeast species C. jadinii, P. kudriavzevii, and K. servazzii in specialty cheeses. These yeasts have been previously detected in cheeses from Québec (
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
). The use of qPCR in this study provided the first quantitative assessment of indigenous yeast presence in both the rind and the core of cheeses with different rind varieties (bloomy-, washed-, and natural-rind cheeses), and from milk undergoing different heat treatments (raw, thermized, and pasteurized).
Before this study, the presence of C. jadinii was rarely reported in cheese. It had been isolated from, but not quantified in, the rind of a Livarot cheese (a soft, red washed-rind cheese made from raw milk), during an early stage of ripening before salting (
). In an earlier study by our group, using traditional microbiological methods, C. jadinii was isolated in only a single raw milk sample and a single semi-hard washed-rind cheese core made with thermized milk (
). Interestingly, the qPCR method detected and quantified C. jadinii in 83% of surfaces and 86% of cores, a higher frequency than previously reported. The low detection rate in previous studies could be explained by the lack of specific detection methods (
) or by the presence of C. jadinii as a subdominant species of the cheese ecosystem because it is masked by other microorganisms. Given that C. jadinii was detected in relatively large quantities in several cheese varieties, the effect of C. jadinii on the physicochemical characteristics and microbial ecosystem of cheeses during ripening warrants further investigation.
Kazachstania servazzii was detected in only 3 varieties of specialty cheese analyzed in this study: on 2 natural-rind cheeses (C28 and C29) and in 1 washed-rind cheese core (C12). Kazachstania servazzii was found in the rind (surface) of all 3 biological and 27 technical replicates of C28. It was also isolated from the same cheese in 2012 (
), suggesting that this yeast is well established in the cheese plant where C28 is produced. Kazachstania servazzii is rarely reported in cheese, with the exception of a Camembert-type cheese produced using a Tibetan freeze-dried kefir coculture as a starter (
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
). This ecosystem clearly evolved during the ripening process, as D. hansenii was found in higher concentrations at the end of the ripening period (d 10 to 35) of the cheese. These results, taken in conjunction with the current report, suggest that K. servazzii is not well adapted to survive in the cheese matrix or that it is displaced by other yeasts during ripening. This could be explained by the ability of K. servazzii to metabolize galactose but not lactate (
). Depending on techniques used during ripening (e.g., lactic vs. rennet coagulum, soft vs. hard cheese), residual galactose can be limited relative to lactate during the ripening process (
In previous work, P. kudriavzevii was found in raw milk from 8 farms and was the most frequently isolated yeast species from ripened cheese: 7 isolates came from cheese rinds and 12 from cheese cores (
). Pichia kudriavzevii is naturally present in milk, can survive the cheesemaking process, and can be found in cheese. In this study, P. kudriavzevii was found in 26 of 29 cheeses analyzed (6 surface and 24 cores), in all cheese rind varieties (washed, bloomy, and natural rind) and in cheese from milk submitted to all types of heat treatment (raw, thermized, and pasteurized). Our results confirm that P. kudriavzevii is frequently found in cheeses; it has been reported in Austrian fresh and sour curd cheeses (
). Pichia kudriavzevii was found to be one of the dominant species on the surface of Cotija raw milk cheese, as identified by PCR-denaturing gradient gel electrophoresis (PCR-DGGE;
). This evidence, in conjunction with our results, suggests that the presence of P. kudriavzevii in cheese might have been underestimated in previous studies using traditional microbiological methods (
Although P. kudriavzevii appears to be common in cheese microbial ecosystems, its contribution to the cheese ripening process remains unclear. It is generally accepted that yeasts are found in 10- to 1,000-fold-higher concentrations on the surface of cheese than in the core because of their oxygen and nutritional requirements (
). However, our qPCR quantification method showed that the indigenous yeasts C. jadinii and P. kudriavzevii can be found in high concentrations (up to 104–108 gene copies/g) in the cheese core compared with the surface (101–104 gene copies/cm2, equivalent to 102–105 gene copies/g). As the presence and metabolism of indigenous yeasts in the core and surface could affect the flavor development of the cheese (
), further investigations of the role of C. jadinii and P. kudriavzevii in cheese ripening are warranted. Furthermore, the presence of these particular yeasts may reflect intrinsic properties, such as the ability to grow under low oxygen conditions or in the presence of lactate and citrate (
Previous studies have shown that the type of heat treatment applied to milk affects the biodiversity of microorganisms found in cheese made from that milk (
). The quantity of C. jadinii in the cheese core increased with milk thermization compared with cheese from raw and pasteurized milk (P < 0.01; Figure 1). Conversely, P. kudriavzevii was more abundant in bloomy-rind cheese made from raw milk (P < 0.05; Figure 3). However, caution should be exercised in the interpretation of these results because initial concentrations of yeasts in the milk were not established for the cheeses studied. Changes in indigenous yeast counts from the raw milk to the ripened cheese surface and core should be further studied to provide insight into the sources and effects of these microorganisms.
CONCLUSIONS
Using species-specific primers and probes, we developed a highly specific and sensitive real-time quantitative PCR method to analyze the presence and distribution of Cyberlindnera jadinii, Pichia kudriavzevii, and Kazachstania servazzii in specialty cheeses. Results confirmed that these 3 indigenous yeasts are present in different commercial cheeses produced in the province of Québec and revealed that C. jadinii and P. kudriavzevii are found in especially large quantities in the cheese core. This study highlights the importance of research not only on the microbial dynamics of the cheese core but also on the effects of these indigenous yeast species on the physicochemical properties of specialty cheeses.
ACKNOWLEDGMENTS
This work was supported by Novalait Inc., the Fonds de Recherche Nature et Technologies du Québec (FRQNT; Québec, Canada), the Québec Consortium for Industrial Bioprocess Research and Innovation (CRIBIQ; Québec, Canada), and the Natural Sciences and Engineering Research Council of Canada (NSERC; Ottawa, Canada; Grant ID 2016-LG-195927). The authors are grateful to Laval University's Service de Consultation Statistique for their assistance. The authors have not stated any conflicts of interest.
REFERENCES
Al-Samarrai T.H.
Schmid J.
A simple method for extraction of fungal genomic DNA.
Metatranscriptome analysis of fungal strains Penicillium camemberti and Geotrichum candidum reveal cheese matrix breakdown and potential development of sensory properties of ripened Camembert-type cheese.
Production of traditional sorghum beer “Ikigage” using Saccharomyces cerevisae, Lactobacillus fermentum and Issatckenkia orientalis as starter cultures.
Microbial diversity of a Camembert-type cheese using freeze-dried Tibetan kefir coculture as starter culture by culture-dependent and culture-independent methods.
00523-9&id=gr1.jpg){kind=link}
00523-9&id=gr2.jpg){kind=link}
00523-9&id=gr3.jpg){kind=link}