The Good, the Bad, and the Ugly: Tales of Mold-Ripened Cheese

Deciphering microbial communities of three Savoyard raw milk cheeses along ripening and regarding the cheese process

Different Savoyard cheeses are granted with PDO (Protected Designation or Origin) and PGI (Protected Geographical Indication) which guarantees consumers compliance with strict specifications. The use of raw milk is known to be crucial for specific flavor development. To unravel the factors influencing microbial ecosystems across cheese making steps, according to the seasonality (winter and summer) and the mode of production (farmhouse and dairy factory ones), gene targeting on bacteria and fungus was used to have a full picture of 3 cheese making technologies, from the raw milk to the end of the ripening. Our results revealed that Savoyard raw milks are a plenteous source of biodiversity together with the brines used during the process, that may support the development of specific features for each cheese. It was shown that rinds and curds have very contrasted ecosystem diversity, composition, and evolution. Ripening stage was selective for some bacterial species, whereas fungus were mainly ubiquitous in dairy samples. All ripening stages are impacted by the type of cheese technologies, with a higher impact on bacterial communities, except for fungal rind communities, for which the technology is the more discriminant. The specific microorganism's abundance for each technology allow to see a real bar-code, with more or less differences regarding bacterial or fungal communities. Bacterial structuration is shaped mainly by matrices, differently regarding technologies while the influence of technology is higher for fungi. Production types showed 10 differential bacterial species, farmhouses showed more ripening taxa, while dairy factory products showing more lactic acid bacteria. Meanwhile, seasonality looks to be a minor element for the comprehension of both microbial ecosystems, but the uniqueness of each dairy plant is a key explicative feature, more for bacteria than for fungus communities.

Analysis of Microbiota Persistence in Quebec's Terroir Cheese Using a Metabarcoding Approach

Environmental short amplicon sequencing, or metabarcoding, is commonly used to characterize the bacterial and fungal microbiota of cheese. Comparisons between different metabarcoding studies are complicated by the use of different gene markers. Here, we systematically compare different metabarcoding molecular targets using V3–V4 and V6–V8 regions of the bacterial 16S rDNA and fungal ITS1 and ITS2 regions. Taxonomic profiles varied depending on the molecular markers used. Based on data quality and detection capacity of the markers toward microorganisms usually associated with the dairy environment, the ribosomal regions V3–V4 and ITS2 were selected and further used to evaluate variability in the microbial ecosystem of terroir cheeses from the province of Quebec in Canada. Both fungal and bacterial ecosystem profiles were described for 32 different ready-to-eat bloomy-, washed- and natural-rind specialty cheese varieties. Among them, 15 were studied over two different production years. Using the Bray–Curtis dissimilarity index as an indicator of microbial shifts, we found that most variations could be explained by either a voluntary change in starter or ripening culture composition, or by changes in the cheesemaking technology. Overall, our results suggest the persistence of the microbiota between the two years studied—these data aid understanding of cheese microbiota composition and persistence during cheese ripening.

In-Depth Investigation of the Safety of Wooden Shelves Used for Traditional Cheese Ripening

The main goal of this research was to characterize the bacterial diversity of the wooden boards used for aging traditional Sicilian cheeses and to evaluate whether pathogenic bacteria are associated with these surfaces. Eighteen cheese dairy factories producing three traditional cheese typologies (PDO Pecorino Siciliano, PDO Piacentinu Ennese, and Caciocavallo Palermitano) were selected within the region of Sicily. The wooden shelf surfaces were sampled by a destructive method to detach wood splinters as well as by a nondestructive brushing to collect microbial cells. Scanning electron microscopy showed the presence of almost continuous bacterial formations on the majority of the shelves analyzed. Yeasts and fungal hyphae were also visualized, indicating the complexity of the plank communities. The amplicon library of the 16S rRNA gene V3-V4 region was paired-end sequenced using the Illumina MiSeq system, allowing the identification of 14 phyla, 32 classes, 52 orders, 93 families, and 137 genera. Staphylococcus equorum was identified from all wooden surfaces, with a maximum abundance of 64.75%. Among cheese-surface-ripening bacteria, Brevibacterium and Corynebacterium were detected in almost all samples. Several halophilic (Halomonas, Tetragenococcus halophilus, Chromohalobacter, Salimicrobium, Marinococcus, Salegentibacter, Haererehalobacter, Marinobacter, and Idiomarinaceae) and moderately halophilic (Salinicoccus, Psychrobacter, and Salinisphaera) bacteria were frequently identified. Lactic acid bacteria (LAB) were present at low percentages in the genera Leuconostoc, Lactococcus, Lactobacillus, Pediococcus, and Streptococcus. The levels of viable microorganisms on the wooden shelves ranged between 2.4 and 7.8 log CFU/cm2. In some cases, LAB were counted at very high levels (8.2 log CFU/cm2). Members of the Enterobacteriaceae family were detected in a viable state for only six samples. Coagulase-positive staphylococci, Salmonella spp., and Listeria monocytogenes were not detected. Seventy-five strains belonged to the genera Leuconostoc, Lactococcus, Pediococcus, Enterococcus, Lactobacillus, and Weissella. IMPORTANCE This study provides evidence for the lack of pathogenic bacteria on the wooden shelves used to ripen internal bacterially ripened semihard and hard cheeses produced in Sicily. These three cheeses are not inoculated on their surfaces, and surface ripening is not considered to occur or, at least, does not occur at the same extent as surface-inoculated smear cheeses. Several bacterial groups identified from the wooden shelves are typically associated with smear cheeses, strongly suggesting that PDO Pecorino Siciliano, PDO Piacentinu Ennese, and Caciocavallo Palermitano cheese rind contributes to their final organoleptic profiles.

Composition, Succession, and Source Tracking of Microbial Communities throughout the Traditional Production of a Farmstead Cheese

Prior to the advent of milk pasteurization and the use of defined-strain starter cultures, the production and ripening of cheese relied on the introduction and growth of adventitious microbes from the environment. This study characterized microbial community structures throughout a traditional farmstead cheese production continuum and evaluated the role of the environment in microbial transfer. In total, 118 samples (e.g., raw milk, cheese, and environmental surfaces) were collected from milk harvesting through cheese ripening. Microbial communities were characterized based on amplicon sequencing of bacterial 16S rRNA and fungal internal transcribed spacer genes using the Illumina MiSeq platform. Results indicated that the environment in each processing room harbored unique microbial ecosystems and consistently contributed microbes to milk, curd, and cheese. The diverse microbial composition of milk was initially attributed to milker hands and cow teats and then changed substantially following overnight ripening in a wooden vat to one dominated by lactic acid bacteria, including Lactococcus lactis, Lactobacillus, and Leuconostoc, as well as fungi such as Exophiala, Kluyveromyces, and Candida. Additional microbial contributions were attributed to processing tools, but the composition of the cheese paste remained relatively stable over 60 days of ripening. In contrast, rind communities that were largely influenced by direct contact with bamboo aging mats showed a distinct succession pattern compared to the interior sections. Overall, these findings highlight the critical role of traditional tools and practices in shaping the microbial composition of cheese and broaden our understanding of processing environments as important sources of microbes in food. IMPORTANCE Throughout the 20th century, especially in the United States, sanitation practices, pasteurization of milk, and the use of commercial defined-strain starter cultures have enhanced the safety and consistency of cheese. However, these practices can reduce cheese microbial diversity. The rapid growth of the artisanal cheese industry in the United States has renewed interest in recapturing the diversity of dairy products and the microbes involved in their production. Here, we demonstrate the essential role of the environment, including the use of wooden tools and cheesemaking equipment, as sources of dominant microbes that shape the fermentation and ripening processes of a traditional farmstead cheese produced without the addition of starter cultures or direct inoculation of any other bacteria or fungi. These data enrich our understanding of the microbial interactions between products and the environment and identify taxa that contribute to the microbial diversity of cheese and cheese production.

Occurrence of Ochratoxin A in Different Types of Cheese Offered for Sale in Italy

The detection of Ochratoxin A (OTA) in the milk of ruminants occurs infrequently and at low levels, but its occurrence may be higher in dairy products such as cheese. The aim of this study was to investigate the presence of OTA in cheeses purchased in the metropolitan city of Bologna (Italy) and the surrounding area. For the analysis, a LC-MS/MS method with a limit of quantification (LOQ) of 1 µg/kg was used. OTA was detected in seven out of 51 samples of grated hard cheese (concentration range 1.3-22.4 µg/kg), while it was not found in the 33 cheeses of other types which were also analysed. These data show a low risk of OTA contamination for almost all types of cheese analysed. To improve the safety of cheese marketed in grated form, more regulations on cheese rind, which is the part most susceptible to OTA-producing moulds, should be implemented or, alternatively, producers should consider not using the rind as row material for grated cheese. It would be interesting to continue these investigations particularly on grated hard cheeses to have more data to update the risk assessment of OTA in cheese, as also suggested by EFSA in its 2020 scientific opinion on OTA.

Modelling the Radial Growth of Geotrichum candidum: Effects of Temperature and Water Activity

Modelling the growth of microorganisms in relation to environmental factors provides quantitative knowledge that can be used to predict their behaviour in foods. For this reason, the effects of temperature and water activity (aw) adjusted with NaCl on the surface growth of two isolates and one culture strain of Geotrichum candidum were studied. A dataset of growth parameters obtained from almost 600 growth curves was employed for secondary modelling with cardinal models (CMs). The theoretical minimal temperature resulting from the modelling of the mycelium proliferation rate ranged from -5.2 to -0.4 °C. Optimal and maximal temperatures were calculated and found to have narrow ranges of 25.4 to 28.0 °C and 34.2 to 37.6 °C, respectively. Cardinal aw values associated with radial growth (awmin from 0.948-0.960 and awopt from 0.992-0.993) confirmed the salt sensitivity of the species. Model goodness-of-fit was evaluated by the coefficient of determination R2, which ranged from 0.954 to 0.985, and RMSE, which ranged from 0.28 to 0.42. Substantially higher variability accompanied the lag time for growth modelling than the radial growth rate modelling despite the square root transformation of the reciprocal lag phase data (R2 = 0.685 to 0.808). Nevertheless, the findings demonstrate that the outputs of growth modelling can be applied to the quantitative evaluation of the roles of G. candidum in fresh cheese spoilage as well as the ripening of Camembert-type cheeses or various artisanal cheeses. Along with validation, the interactions with lactic acid bacteria can be included to improve the predictions of G. candidum in the future.

Strain-Level Diversity Impacts Cheese Rind Microbiome Assembly and Function

Diversification can generate genomic and phenotypic strain-level diversity within microbial species. This microdiversity is widely recognized in populations, but the community-level consequences of microbial strain-level diversity are poorly characterized. Using the cheese rind model system, we tested whether strain diversity across microbiomes from distinct geographic regions impacts assembly dynamics and functional outputs. We first isolated the same three bacterial species (Staphylococcus equorum, Brevibacterium auranticum, and Brachybacterium alimentarium) from nine cheeses produced in different regions of the United States and Europe to construct nine synthetic microbial communities consisting of distinct strains of the same three bacterial species. Comparative genomics identified distinct phylogenetic clusters and significant variation in genome content across the nine synthetic communities. When we assembled each synthetic community with initially identical compositions, community structure diverged over time, resulting in communities with different dominant taxa. The taxonomically identical communities showed differing responses to abiotic (high salt) and biotic (the fungus Penicillium) perturbations, with some communities showing no response and others substantially shifting in composition. Functional differences were also observed across the nine communities, with significant variation in pigment production (light yellow to orange) and in composition of volatile organic compound profiles emitted from the rinds (nutty to sulfury).IMPORTANCE Our work demonstrated that the specific microbial strains used to construct a microbiome could impact the species composition, perturbation responses, and functional outputs of that system. These findings suggest that 16S rRNA gene taxonomic profiles alone may have limited potential to predict the dynamics of microbial communities because they usually do not capture strain-level diversity. Observations from our synthetic communities also suggest that strain-level diversity has the potential to drive variability in the aesthetics and quality of surface-ripened cheeses.

Fungal networks shape dynamics of bacterial dispersal and community assembly in cheese rind microbiomes

Most studies of bacterial motility have examined small-scale (micrometer–centimeter) cell dispersal in monocultures. However, bacteria live in multispecies communities, where interactions with other microbes may inhibit or facilitate dispersal. Here, we demonstrate that motile bacteria in cheese rind microbiomes use physical networks created by filamentous fungi for dispersal, and that these interactions can shape microbial community structure. Serratia proteamaculans and other motile cheese rind bacteria disperse on fungal networks by swimming in the liquid layers formed on fungal hyphae. RNA-sequencing, transposon mutagenesis, and comparative genomics identify potential genetic mechanisms, including flagella-mediated motility, that control bacterial dispersal on hyphae. By manipulating fungal networks in experimental communities, we demonstrate that fungal-mediated bacterial dispersal can shift cheese rind microbiome composition by promoting the growth of motile over non-motile community members. Our single-cell to whole-community systems approach highlights the interactive dynamics of bacterial motility in multispecies microbiomes.

Interactions with other microbes may inhibit or facilitate the dispersal of bacteria. Here, Zhang et al. use cheese rind microbiomes as a model to show that physical networks created by filamentous fungi can affect the dispersal of motile bacteria and thus shape the diversity of microbial communities.

Biodiversity of the Surface Microbial Consortia from Limburger, Reblochon, Livarot, Tilsit, and Gubbeen Cheeses

Comprehensive collaborative studies from our laboratories reveal the extensive biodiversity of the microflora of the surfaces of smear-ripened cheeses. Two thousand five hundred ninety-seven strains of bacteria and 2,446 strains of yeasts from the surface of the smear-ripened cheeses Limburger, Reblochon, Livarot, Tilsit, and Gubbeen, isolated at three or four times during ripening, were identified; 55 species of bacteria and 30 species of yeast were found. The microfloras of the five cheeses showed many similarities but also many differences and interbatch variation. Very few of the commercial smear microorganisms, deliberately inoculated onto the cheese surface, were reisolated and then mainly from the initial stages of ripening, implying that smear cheese production units must have an adventitious "house" flora. Limburger cheese had the simplest microflora, containing two yeasts, Debaryomyces hansenii and Geotrichum candidum, and two bacteria, Arthrobacter arilaitensis and Brevibacterium aurantiacum. The microflora of Livarot was the most complicated, comprising 10 yeasts and 38 bacteria, including many gram-negative organisms. Reblochon also had a very diverse microflora containing 8 yeasts and 13 bacteria (excluding gram-negative organisms which were not identified), while Gubbeen had 7 yeasts and 18 bacteria and Tilsit had 5 yeasts and 9 bacteria. D. hansenii was by far the dominant yeast, followed in order by G. candidum, Candida catenulata, and Kluyveromyces lactis. B. aurantiacum was the dominant bacterium and was found in every batch of the 5 cheeses. The next most common bacteria, in order, were Staphylococcus saprophyticus, A. arilaitensis, Corynebacterium casei, Corynebacterium variabile, and Microbacterium gubbeenense. S. saprophyticus was mainly found in Gubbeen, and A. arilaitensis was found in all cheeses but not in every batch. C. casei was found in most batches of Reblochon, Livarot, Tilsit, and Gubbeen. C. variabile was found in all batches of Gubbeen and Reblochon but in only one batch of Tilsit and in no batch of Limburger or Livarot. Other bacteria were isolated in low numbers from each of the cheeses, suggesting that each of the 5 cheeses has a unique microflora. In Gubbeen cheese, several different strains of the dominant bacteria were present, as determined by pulsed-field gel electrophoresis, and many of the less common bacteria were present as single clones. The culture-independent method, denaturing gradient gel electrophoresis, resulted in identification of several bacteria which were not found by the culture-dependent (isolation and rep-PCR identification) method. It was thus a useful complementary technique to identify other bacteria in the cheeses. The gross composition, the rate of increase in pH, and the indices of proteolysis were different in most of the cheeses.

Diversity of yeast and mold species from a variety of cheese types

To generate a comprehensive profile of viable fungi (yeasts and molds) on cheese as it is purchased by consumers, 44 types of cheese were obtained from a local grocery store from 1 to 4 times each (depending on availability) and sampled. Pure cultures were obtained and identified by DNA sequence of the ITS region, as well as growth characteristics and colony morphology. The yeast Debaryomyces hansenii was the most abundant fungus, present in 79 % of all cheeses and 63 % of all samples. Penicillium roqueforti was the most common mold, isolated from a variety of cheeses in addition to the blue cheeses. Eighteen other fungal species were isolated, ten from only one sample each. Most fungi isolated have been documented from dairy products; a few raise potential food safety concerns (i.e. Aspergillus flavus, isolated from a single sample and capable of producing aflatoxins; and Candida parapsilosis, an emerging human pathogen isolated from three cheeses). With the exception of D. hansenii (present in most cheese) and P. roqueforti (a necessary component of blue cheese), no strong correlation was observed between cheese type, manufacturer, or sampling time with the yeast or mold species composition.