DNA-Based Methods for Cheese Origin Authentication: Comparison
Please note this is a comparison between Version 2 by Dean Liu and Version 1 by Marco Cardin.

Metabolic activities of cheese microbiota play a crucial role in the development of cheese typicity. For geographical origin authentication, microbiota fingerprinting is therefore of high interest as traditional and artisanal cheeses are produced with a more diversified microbiota associated with the cheese-making process (e.g., use of raw milk, starter, brine, equipment and materials, and ripening rooms). While cheese microbial diversity was traditionally investigated using culture-dependent methods, hence overlooking unculturable or subdominant species, nowadays, culture-independent methods (HTS) have unravelled this diversity and provided further means to connect microbiota composition to cheese quality and typicity, but also origin. This success is due to the availability of new sequencing platforms, bioinformatic pipelines, and a continuous decrease in cost. While complete microbial genome sequencing (i.e. metagenomics) can be used, metabarcoding (also known as metagenetics and corresponding to the selective amplification of taxonomically relevant followed by NGS sequencing) is the most widely reported in the scientific literature. Herein we report recent scientific literature about metabarcoding analysis of cheese microbiota, highlighting what factors contribute to its formation and how it could be used to authenticate cheese origin.

  • cheese
  • geographical origin
  • authentication
  • next-generation sequencing

1. Introduction

To perform DNA metabarcoding, cheese samples are first homogenized, then total DNA is most frequently extracted using commercial kits, ad hoc protocols or a combination thereof [1][2]. Hypervariable regions of taxonomically relevant genes (e.g., 16S rDNA for bacteria and archaea, ITS, 18S rDNA, 26S rDNA for fungi) are amplified by PCR reactions, while a second amplification step tags amplicons with specific DNA fragments –barcodes- and dedicated adapters for the final sequencing step using next-generation technologies (e.g., Illumina, Pacbio, Iontorrent, or Nanopore). Typically, 16S rDNA and ITS (internal transcribed spacer) markers, targeting bacteria and fungi, respectively, are employed to generate compositional data describing microbial taxa and their relative abundance in cheese microbial communities. After sequencing, two complementary but different ways can be used for amplicon clustering from quality-checked data, namely, operational taxonomic units (OTUs) and amplicon sequencing variants (ASVs) [3]. Subsequently, taxonomic assignment is performed using a specific classifier tool (BLAST, RDP, UCLUST, SortMeRNA) against various reference databases, such as Greengenes, SILVA, and UNITE [4][5]. Generally, clustered OTUs/ASVs are analyzed from the phylum to the genus level since they can be less precise at the species level [6].

12. Main Factors Affecting Cheese Microbial Diversity

In a recent study by Kamimura et al. [7], the microbiota of 578 traditional Brazilian cheeses were analyzed by metabarcoding. Bacterial communities were distinctly clustered with PCA by cheese type and regional origins while, at the genus level, hierarchical cluster analysis separated production regions. These authors were thus able to identify specific origin-related microbiota. The core microbiota of Brazilian traditional cheeses displayed different relative abundances and oligotypes (i.e. closely related but distinct bacterial taxa) of LAB belonging to the Enterococcus, Lactococcus, Streptococcus, Leuconostoc, and Lactobacillus sensu lato genera, as well as other taxa belonging to the Enterobacteriaceae family and Staphylococcus genus. Within the same regional area, microbiota analysis differentiated the origin of traditional cheeses, namely Cerrado, Araxà, Canastra, Campos des Vertenes and Serro, produced with a similar natural whey starter and ripening period (17 and 22 days). These findings were in agreement with those of another study that analyzed 97 samples of Minas artisanal cheese from 6 different producers located in the same region [8]. Starter cultures, consisting of Streptococcus, Lactococcus and Lactobacillus sensu lato spp., constituted the core microbiota in all farms. However, significant differences in family- and genus-level bacterial community relative abundances were observed between the studied farms due to environmental factors including their geographical location. Even when dominant genera may be inferred to the natural whey cultures used, meta-analysis from amplicon sequencing data of traditional artisanal cheeses from Italy, Belgium, and Kalmykia indicated differences in bacterial structures between cheeses produced in different geographic areas (unweighted PCoA cluster MANOVA, p < 0.001 [9]). Those produced using natural milk cultures showed improved acidification without an effect on the typical cheese microbiota. Indeed, at the end of the ripening period, cheese origins clustered according to producer facilities (PCoA on Bray–Curtis) [10].
Another study was performed by Zago et al. [11], 118 Grana Padano samples were analyzed after 7-8 months of ripening and a common core microbiota composed of Lactobacillus-, Lactococcus-, Lacticaseibacillus-, Limosilactobacillus- and Streptococcus-dominant genera was observed. More precisely, differences in bacterial abundance, richness, and evenness were found for dominant and sub-dominant groups according to production region, a result also confirmed by PERMANOVA b-diversity analysis. The authors also identified specific species that could be linked to several production areas; however, no species biomarkers were identified regardless of production area and non-metric multidimensional scaling did not show any clear clustering profile.
Some cheeses are produced in very small geographical zones by a limited number of producers. This is the case for Plaisentif and Historic Rebel cheeses from the mountainous regions in Italy that are only produced during specific seasons (violet blooming season and grazing season) by 14 and 12 producers, respectively [12][13]. Both are raw milk cheeses produced without starter adjunction. Bacterial amplicon sequencing analysis (16S rDNA V4 region) for Plaisentif cheese identified dominant genera and, more importantly, differences in bacterial community profiles between producers suggesting fraudulent starter addition in some cheeses [13].
Based on a similar analysis for Historical Rebel cheese, the core microbiota was composed of 5 different genera -Streptococcus, Lactobacillus, Lactococcus, Leuconostoc and Pediococcus- with Streptococcus relative abundances ranging from 60% to 85% [12]. Richness and other a-diversity parameters differed among producers as well as in multivariate analysis (PCoA on unweighted Unifrac) and, based on the observed significant differences, pasture area could be linked to the different Historic Rebel cheese producers.

23. Climatic and Environmental Condition

Other factors that can impact microbial communities of traditional cheeses are the climatic and dairy environment conditions that are directly associated with geographical origins. This was observed for traditional Chinese Rushan cheese produced using Chaenomeles sinensis boiled extract as a clotting agent in three different regions. Even if the same UHT milk and production equipment were used, geographical origins significantly impacted the relative abundance of 12 dominant genera, namely Lactobacillus, Acinetobacter, Acetobacter, Lactococcus, Enterobacter, Moraxella, Enterococcus, Streptococcus, Kocuria, Staphylococcus, Chryseobacterium and Exiguobacterium [14], and is likely related to specific house microbiota and open-air drying. This result was also confirmed by PCoA clusters and Anosim analysis.
As typical cheese microbiota can be in part acquired from the specific raw materials used, traditional tools, environmental and production conditions, cheese-making process and geographical area, a comparison between traditional and industrial cheeses may provide additional information to authenticate cheese origin. Noteworthy, some authors have reported that commercial starters, inoculated at ~106 CFU/mL, prevent resident microbiota from developing, especially during ripening [15][16]. Overall, milk pasteurization, use of similar commercial starters, similar industrial equipment, and standardized recipes for cheese production are crucial factors that decrease cheese microbial complexity and biodiversity, and lead to highly standardized productions and thus lower variance is detected (Figure 1). This hypothesis is in accordance with the study by Kamilari et al. [17], in which a significant decrease in bacterial diversity was observed in industrially produced Haloumi cheeses vs. artisanal products. However, the microbial diversity observed for artisanal Haloumi cheeses could not be linked to the producer geographical origins. In another study aiming to authenticate cheese origin at the producer level, no distinction could again be made [18].
Figure 1. Technological factors affecting cheese microbiota biodiversity. Green and purple lines show combinations of technological factors during cheese-making that increase or decrease this biodiversity.

34. Cheese Ripening

Cheese ripening is another factor that affects microbial community diversity and cheese typicity. Ripening can be considered as a selection process that leads to cheese microbial composition changes. According to Gobetti et al. [19], intentionally added microorganisms used in cheese-making include primary starters (natural milk culture, natural whey culture, or lyophilized commercial starters), secondary or adjunct LAB starters, and milk autochthonous microbiota (Non-starter lactic acid bacteria -NSLAB- and others). These are the main ripening agents in intermediate to long ripening times, which mainly explain the observed diversity and typicity of the produced cheeses. The relationship between primary starters and NSLAB during maturation is well-known, and involves a progressive reduction of the former in favor of the latter. The role of NSLAB is crucial in maturation and the development of the typical characteristics of traditional cheese. In this case, useful insights can be gained by comparing the genomic features of primary starters for the presence of genes involved in the metabolic pathways important for cheese maturation. While primary starters have important genetic features for the utilization of lactose (mainly connected with their acidification ability) NSLAB possess many more genes coding for peptidases, peptide transporters within cells, and amino acid catabolism that can represent an advantage during cheese maturation [20]. Moreover, compared to primary starters, NSLAB tend to adapt better to the hostile conditions of the cheese ripening environment, such as temperature, salt content, pH, and redox potential. In fact, NSLAB can adopt alternative metabolic pathways to produce energy from unconventional sources while resisting acid conditions. Therefore, NSLAB present in raw milk at a sufficient inoculum to colonize ripened raw milk cheeses, or acquired from the house microbiota, could be an indicator of geographical origin at the producer level. Beyond NSLAB, other microorganisms belonging to various groups can influence cheese ripening. This is the case of fungal communities in many traditional cheeses, such as Queijo de Azeitão in Portugal [21], Tomme d’Orchies in France [22], and Robiola di Roccaverano in Italy [23]. Indeed, fungal communities are well-known for their decisive role in flavor and texture of white and blue-veined mold-ripened cheeses due to lipolytic, proteolytic, and glycolytic activities, leading to high production of aromatic ketones and alcohols [24][25][26]. Generally, fungal species such as Penicillium camemberti, Penicillium roqueforti, Debaryomyces hansenii, Kluyveromyces marxianus, Candida catenulata, Galactomyces geotrichum, and Mucor lanceolatus are either deliberately added as technological adjunct cultures or present in the production environment [26][27][28][29]. Nevertheless, in traditional cheeses, fungal communities were reported to be more diverse than the used starter but, at the same time, not connected with geographical origin [1][21]. Table 1 reports advantages and limitations of amplicon sequencing for cheese origin authentication.
Advantages Reference Limitations Reference
Time- and cost-effective processing of large sample numbers [30] Analyses could be biased by sample processing, DNA extraction methods, and equimolar library preparation [31]
Consolidated pipeline for data analysis [32] PCR amplification steps include errors, e.g., PCR specificity and variation of 16S rRNA copy number per genome [33]
Identification of taxonomic groups associated with typical flavor and cheese-making technology [34] Under- or over-estimation of microbial community diversity [32][35]
Allows improvement of cheese-making to ensure safety while preserving typicity [7] Lack of absolute abundance [36][37]
Evaluation of core microbiome describing facility-associated microbial groups [38][39] Limited and uneven taxonomic resolutions [40][41]
May pinpoint new biotypes [17] DNA amplicon sequencing typically does not discriminate between live and dead microorganisms (except if DNA stains such as propidium monoazide are used) [42]
The mentioned studies showed that in traditional cheeses, the combination of artisanal cheese-making, specific raw materials, and characteristic environmental conditions shape microbial community diversity according to geographical origin. Most analyses conducted using 16S rDNA amplicon sequencing discriminated cheese origin, although taxonomic classification was still limited to genus/family-level descriptions and only a few cheese types per study were considered. To further assess the unicity of typical cheeses against food fraud, more in-depth studies, including meta-analyses on all available cheese data and increased depth of microbial population descriptions (e.g., metagenomics), are of interest.

References

  1. Benjamin E. Wolfe; Julie E. Button; Marcela Santarelli; Rachel J. Dutton; Cheese Rind Communities Provide Tractable Systems for In Situ and In Vitro Studies of Microbial Diversity. Cell 2014, 158, 422-433, 10.1016/j.cell.2014.05.041.
  2. Ilhan Cem Duru; Pia Laine; Margarita Andreevskaya; Lars Paulin; Soila Kananen; Soile Tynkkynen; Petri Auvinen; Olli-Pekka Smolander; Metagenomic and metatranscriptomic analysis of the microbial community in Swiss-type Maasdam cheese during ripening. International journal of food microbiology 2018, 281, 10-22, 10.1016/j.ijfoodmicro.2018.05.017.
  3. Sandra Reitmeier; Thomas C. A. Hitch; Nicole Treichel; Nikolaos Fikas; Bela Hausmann; Amanda E. Ramer-Tait; Klaus Neuhaus; David Berry; Dirk Haller; Ilias Lagkouvardos; et al.Thomas Clavel Handling of spurious sequences affects the outcome of high-throughput 16S rRNA gene amplicon profiling. ISME Communications 2021, 1, 1-12, 10.1038/s43705-021-00033-z.
  4. Christian Quast; Elmar Pruesse; Pelin Yilmaz; Jan Gerken; Timmy Schweer; Pablo Yarza; Jörg Peplies; Frank Oliver Glöckner; The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Research 2012, 41, D590-D596, 10.1093/nar/gks1219.
  5. Henrik Nilsson; Karl-Henrik Larsson; Andy F S Taylor; Johan Bengtsson-Palme; Thomas Stjernegaard Jeppesen; Dmitry Schigel; Peter Kennedy; Kathryn Picard; Frank Oliver Glöckner; Leho Tedersoo; et al.Irja SaarUrmas KõljalgKessy Abarenkov The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Research 2018, 47, D259-D264, 10.1093/nar/gky1022.
  6. Orna Mizrahi-Man; Emily Davenport; Yoav Gilad; Taxonomic Classification of Bacterial 16S rRNA Genes Using Short Sequencing Reads: Evaluation of Effective Study Designs. Plos One 2013, 8, e53608, 10.1371/journal.pone.0053608.
  7. Kamimura, B.A.; De Filippis, F.; Sant’Ana, A.S.; Ercolini, D. Large-Scale Mapping of Microbial Diversity in Artisanal Brazilian Cheeses. Food Microbiol. 2019, 80, 40–49.
  8. Sant’Anna, F.M.; Wetzels, S.U.; Cicco, S.H.S.; Figueiredo, R.C.; Sales, G.A.; Figueiredo, N.C.; Nunes, C.A.; Schmitz-Esser, S.; Mann, E.; Wagner, M.; et al. Microbial Shifts in Minas Artisanal Cheeses from the Serra Do Salitre Region of Minas Gerais, Brazil throughout Ripening Time. Food Microbiol. 2019, 82, 349–362.
  9. Li, J.; Zheng, Y.; Xu, H.; Xi, X.; Hou, Q.; Feng, S.; Wuri, L.; Bian, Y.; Yu, Z.; Kwok, L.-Y.; et al. Bacterial Microbiota of Kazakhstan Cheese Revealed by Single Molecule Real Time (SMRT) Sequencing and Its Comparison with Belgian, Kalmykian and Italian Artisanal Cheeses. BMC Microbiol. 2017, 17, 13.
  10. Lucchini, R.; Cardazzo, B.; Carraro, L.; Negrinotti, M.; Balzan, S.; Novelli, E.; Fasolato, L.; Fasoli, F.; Farina, G. Contribution of Natural Milk Culture to Microbiota, Safety and Hygiene of Raw Milk Cheese Produced in Alpine Malga. Ital. J. Food Saf. 2018, 7, 6967.
  11. Zago, M.; Bardelli, T.; Rossetti, L.; Nazzicari, N.; Carminati, D.; Galli, A.; Giraffa, G. Evaluation of Bacterial Communities of Grana Padano Cheese by DNA Metabarcoding and DNA Fingerprinting Analysis. Food Microbiol. 2021, 93, 103613.
  12. Turri, F.; Cremonesi, P.; Battelli, G.; Severgnini, M.; Brasca, M.; Gandini, G.; Pizzi, F. High Biodiversity in a Limited Mountain Area Revealed in the Traditional Production of Historic Rebel Cheese by an Integrated Microbiota–Lipidomic Approach. Sci. Rep. 2021, 11, 10374.
  13. Dalmasso, A.; Soto del Rio, M.d.l.D.; Civera, T.; Pattono, D.; Cardazzo, B.; Bottero, M.T. Characterization of Microbiota in Plaisentif Cheese by High-Throughput Sequencing. LWT-Food Sci. Technol. 2016, 69, 490–496.
  14. Xue, J.; Yang, Y.; Wang, Z.; Guo, Y.; Shao, Y. Bacterial Diversity in Chinese Rushan Cheese From Different Geographical Origins. Front. Microbiol. 2018, 9, 1920.
  15. Piqueras, J.; Chassard, C.; Callon, C.; Rifa, E.; Theil, S.; Lebecque, A.; Delbès, C. Lactic Starter Dose Shapes S. Aureus and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses. Microorganisms 2021, 9, 1081.
  16. Johnson, J.; Curtin, C.; Waite-Cusic, J. The Cheese Production Facility Microbiome Exhibits Temporal and Spatial Variability. Front. Microbiol. 2021, 12, 644828.
  17. Kamilari, E.; Anagnostopoulos, D.A.; Papademas, P.; Kamilaris, A.; Tsaltas, D. Characterizing Halloumi Cheese’s Bacterial Communities through Metagenomic Analysis. LWT 2020, 126, 109298.
  18. Aldrete-Tapia, A.; Escobar-Ramírez, M.C.; Tamplin, M.L.; Hernández-Iturriaga, M. High-Throughput Sequencing of Microbial Communities in Poro Cheese, an Artisanal Mexican Cheese. Food Microbiol. 2014, 44, 136–141.
  19. Gobbetti, M.; Di Cagno, R.; Calasso, M.; Neviani, E.; Fox, P.F.; De Angelis, M. Drivers That Establish and Assembly the Lactic Acid Bacteria Biota in Cheeses. Trends Food Sci. Technol. 2018, 78, 244–254.
  20. Porcellato, D.; Skeie, S.B. Bacterial Dynamics and Functional Analysis of Microbial Metagenomes during Ripening of Dutch-Type Cheese. Int. Dairy J. 2016, 61, 182–188.
  21. Cardinali, F.; Ferrocino, I.; Milanović, V.; Belleggia, L.; Corvaglia, M.R.; Garofalo, C.; Foligni, R.; Mannozzi, C.; Mozzon, M.; Cocolin, L.; et al. Microbial Communities and Volatile Profile of Queijo de Azeitão PDO Cheese, a Traditional Mediterranean Thistle-Curdled Cheese from Portugal. Food Res. Int. 2021, 147, 110537.
  22. Ceugniez, A.; Taminiau, B.; Coucheney, F.; Jacques, P.; Delcenserie, V.; Daube, G.; Drider, D. Fungal Diversity of “Tomme d’Orchies” Cheese during the Ripening Process as Revealed by a Metagenomic Study. Int. J. Food Microbiol. 2017, 258, 89–93.
  23. Biolcati, F.; Bottero, M.T.; Dalmasso, A. Microbiological Analysis of the Robiola Di Roccaverano Cheese by Means of Traditional Culturing Methods. Ital. J. Food Saf. 2019, 8.
  24. Montel, M.-C.; Buchin, S.; Mallet, A.; Delbes-Paus, C.; Vuitton, D.A.; Desmasures, N.; Berthier, F. Traditional Cheeses: Rich and Diverse Microbiota with Associated Benefits. Int. J. Food Microbiol. 2014, 177, 136–154.
  25. Martín, J.F.; Coton, M. Blue cheese: Microbiota and fungal metabolites. In Fermented Foods in Health and Disease Prevention; Academic Press: San Diego, CA, USA, 2017; pp. 275–303.
  26. Coton, E.; Coton, M.; Hymery, N.; Mounier, J.; Jany, J.-L. Penicillium Roqueforti: An Overview of Its Genetics, Physiology, Metabolism and Biotechnological Applications. Fungal Biol. Rev. 2020, 34, 59–73.
  27. Gillot, G.; Jany, J.-L.; Poirier, E.; Maillard, M.-B.; Debaets, S.; Thierry, A.; Coton, E.; Coton, M. Functional Diversity within the Penicillium Roqueforti Species. Int. J. Food Microbiol. 2017, 241, 141–150.
  28. Garnier, L.; Valence, F.; Mounier, J. Diversity and Control of Spoilage Fungi in Dairy Products: An Update. Microorganisms 2017, 5, 42.
  29. Mounier, J.; Coton, M. Kluyveromyces spp. In Encyclopedia of Dairy Sciences; McSweeney, P.L.H., McNamara, J.P., Eds.; Elsevier: San Diego, CA, USA, 2022; Volume 4, pp. 569–574.
  30. Lundin, D.; Severin, I.; Logue, J.B.; Östman, Ö.; Andersson, A.F.; Lindström, E.S. Which Sequencing Depth Is Sufficient to Describe Patterns in Bacterial α- and β-Diversity? Sequencing Depth in Diversity Research. Environ. Microbiol. Rep. 2012, 4, 367–372.
  31. Rubin, B.E.R.; Sanders, J.G.; Hampton-Marcell, J.; Owens, S.M.; Gilbert, J.A.; Moreau, C.S. DNA Extraction Protocols Cause Differences in 16S RRNA Amplicon Sequencing Efficiency but Not in Community Profile Composition or Structure. MicrobiologyOpen 2014, 3, 910–921.
  32. Nearing, J.T.; Douglas, G.M.; Comeau, A.M.; Langille, M.G.I. Denoising the Denoisers: An Independent Evaluation of Microbiome Sequence Error-Correction Approaches. PeerJ 2018, 6, e5364.
  33. Kleerebezem, R.; Stouten, G.; Koehorst, J.; Langenhoff, A.; Schaap, P.; Smidt, H. Experimental Infrastructure Requirements for Quantitative Research on Microbial Communities. Curr. Opin. Biotechnol. 2021, 67, 158–165.
  34. Unno, R.; Suzuki, T.; Matsutani, M.; Ishikawa, M. Evaluation of the Relationships Between Microbiota and Metabolites in Soft-Type Ripened Cheese Using an Integrated Omics Approach. Front. Microbiol. 2021, 12, 681185.
  35. Hugerth, L.W.; Andersson, A.F. Analysing Microbial Community Composition through Amplicon Sequencing: From Sampling to Hypothesis Testing. Front. Microbiol. 2017, 8, 1561.
  36. Piwosz, K.; Shabarova, T.; Pernthaler, J.; Posch, T.; Šimek, K.; Porcal, P.; Salcher, M.M. Bacterial and Eukaryotic Small-Subunit Amplicon Data Do Not Provide a Quantitative Picture of Microbial Communities, but They Are Reliable in the Context of Ecological Interpretations. mSphere 2020, 5, e00052-20.
  37. Yeh, Y.-C.; Needham, D.M.; Sieradzki, E.T.; Fuhrman, J.A. Taxon Disappearance from Microbiome Analysis Reinforces the Value of Mock Communities as a Standard in Every Sequencing Run. mSystems 2018, 3, e00023-18.
  38. Bokulich, N.A.; Mills, D.A. Facility-Specific “House” Microbiome Drives Microbial Landscapes of Artisan Cheesemaking Plants. Appl. Environ. Microbiol. 2013, 79, 5214–5223.
  39. Calasso, M.; Ercolini, D.; Mancini, L.; Stellato, G.; Minervini, F.; Di Cagno, R.; De Angelis, M.; Gobbetti, M. Relationships among House, Rind and Core Microbiotas during Manufacture of Traditional Italian Cheeses at the Same Dairy Plant. Food Microbiol. 2016, 54, 115–126.
  40. Brumfield, K.D.; Huq, A.; Colwell, R.R.; Olds, J.L.; Leddy, M.B. Microbial Resolution of Whole Genome Shotgun and 16S Amplicon Metagenomic Sequencing Using Publicly Available NEON Data. PLoS ONE 2020, 15, e0228899.
  41. Johnson, J.S.; Spakowicz, D.J.; Hong, B.-Y.; Petersen, L.M.; Demkowicz, P.; Chen, L.; Leopold, S.R.; Hanson, B.M.; Agresta, H.O.; Gerstein, M.; et al. Evaluation of 16S RRNA Gene Sequencing for Species and Strain-Level Microbiome Analysis. Nat. Commun. 2019, 10, 5029.
  42. Afshari, R.; Pillidge, C.J.; Dias, D.A.; Osborn, A.M.; Gill, H. Cheesomics: The Future Pathway to Understanding Cheese Flavour and Quality. Crit. Rev. Food Sci. Nutr. 2020, 60, 33–47.
More
Video Production Service