Abstract
Robiola di Roccaverano is an artisanal Protected Designation of Origin (PDO) soft cheese made with raw goat’s milk and by the addition of Natural Milk Culture (NMC) to drive the fermentation process. Cheeses collected from five different dairy plants were analyzed for their bacterial and fungal microbiota diversity. Lactococcus lactis and Leuconostoc mesenteroides were the main bacterial population, while Galactomyces candidum and Kluyveromyces marxianus constituted the core mycobiota but many other minor taxa were observed, suggesting a high level of complexity in fungal composition by these cheeses compared to bacteria population.
Introduction
Robiola di Roccaverano is a Protected Designation of Origin (PDO) soft cheese produced with raw goat’s milk, in the Roccaverano area of the Piedmont region (North–West Italy) by a consortium of 15 dairies following strict procedural guidelines. During the cheese manufacturing process, raw milk from two subsequent milking sessions (evening and morning) is inoculated with Natural Milk Cultures (NMCs), obtained by back-slopping process in which the previous days’ fermented material is inoculated into the fresh milk to start the production process. Then, the material is also inoculated with animal rennet and left to coagulate for 8–36 h at room temperature. Molding procedure takes place in perforated plastic containers for 24–48 h, where the cheeses are periodically turned over to allow the whey to release and manually dry salted in both sides (). Robiola di Roccaverano can be consumed as fresh cheese after 5 days from the beginning of production or ripened for 15 days. During the production of Robiola di Roccaverano, the cheesemakers naturally select the best performing NMC that are well adapted to the dairy environment, able to ensure reproducibility and standardization of the final product without losing its typicality. Therefore, the sensory characteristics of Robiola di Roccaverano are highly influenced by the action of the microorganism deriving from raw milk and NMC (, , ), as well as by contaminating microorganisms often deriving from sources such as operators, dairy equipment, and environment (). Despite understudied, the fungal composition of the dairy products, commonly referred to as mycobiota, could strongly influence the composition of the final cheese and must be considered for a better characterization of PDO cheese.
In this work, the microbiome, both bacterial and fungal composition, of Robiola di Roccaverano has been studied by means of High Throughput Sequencing (HTS) approach. Cheese samples from different cheesemakers have been collected since the bacterial microbiota and the mycobiota may be associated with the geographical area, season, dairy environment, or raw materials used for the production. Nevertheless, due to the strict manufacture guidelines, it might be possible to define a common core microbial population, typical of Robiola di Roccaverano PDO cheese.
Materials and Methods
Samples Collection
Cheeses after 5 and 15 days of ripening were collected from five dairy plants (A, B, C, D, and E) located in different Roccaverano area of the Piedmont region. Six different batches were collected from each dairy plant (3 for 5-day and 15-day ripened cheese, respectively). For each batch, sampling took place on the same days and at the same time of year (January–February). All the producers used only goat’s milk for the production with the addition of NMC made by each cheesemaker. The 30 samples were collected aseptically and transported in laboratory at 4–6°C and kept at –20°C before the analysis.
DNA Extraction
Ten grams of cheese were homogenized in a sterile stomacher bag with 90 milliliter of sterile Ringer solution (Oxoid, Milan, Italy) and mixed in a Stomacher 400 Circulator (Seward, United Kingdom) at 300 rpm for 3 min. Due to the soft texture of Robiola di Roccaverano cheese, no distinction between core and rind was made during sampling. One milliliter of the first decimal dilution was transferred into 1.5 ml micro-tube, and DNA was extracted by using the DNeasy Blood & Tissue kit (Qiagen, Hilden, Germany) following manufacturer’s instruction. DNA was extracted and quantified by using the QUBIT dsHS kit (Thermo Fisher, Milan, Italy), and purity was checked by absorbance measure using the spectrophotometer (Nanodrop 2000, Thermo Fisher Scientific, Waltham, MA, United States).
High-Throughput Sequencing Analysis
The bacterial community was analyzed by amplifying the V3–V4 hypervariable region of the 16S rRNA gene using primers and procedure described by . The fungal population was studied by amplification of the D1 domain of the 26S rRNA gene using primers and conditions described by . 26S rRNA target was chosen because it showed greater taxonomic resolution and robustness compared to the ITS region. PCR amplicons were purified following the Illumina metagenomic pipeline (Illumina Inc., San Diego, CA, United States). Sequencing was performed with a MiSeq platform (Illumina), generating 250 bp paired-end reads.
Bioinformatic Analysis
Paired-end reads were imported in QIIME2 software for chimeric join and quality filtering (; ; ; ). Amplicon sequence variants (ASVs) generated through DADA2 were used for taxonomic assignment against the Greengenes 16S rRNA gene database for bacteria and using the in-house database of for fungi (; ). ASV taxonomy assignment was then confirmed by manual blast. The ASV table obtained with QIIME2 was rarefied at the lowest number of sequences, and when displayed, the higher taxonomy resolution was reached. When the taxonomy assignment was not able to reach species level, the genus or family name was displayed. PICRUSt () was used to predict the abundance of the KEGG gene families, based on 16S rRNA sequence data (; ). KEGG orthologs were then collapsed at a hierarchy level of 3, and the table was imported into R.
Sequencing data were deposited in the Sequence Read Archive of the National Center for Biotechnology Information (NCBI)1 under Bioproject accession numbers PRJNA704571 and PRJNA626588.
Statistical Analysis
Alpha diversity indices were calculated with vegan package () in R environment. ASV tables obtained through QIIME2 were imported in R2, using made4 function to obtain the Principal Component Analysis (PCA) plots. Venn diagrams were obtained by using Venn Diagram Maker. Inferred metagenome data were analyzed by the made4 package, and hierarchical Ward-linkage clustering, based on the Spearman correlation coefficients of the proportion of the activities belonging to lipid carbohydrates and amino acid metabolism pathways, was used to produce a heatplot.
Results
Bacterial Microbiota Composition of Robiola di Roccaverano Cheese
After sequencing, two samples were excluded from 16S rRNA dataset (one sample from dairy B and one sample from dairy D) due to lower number of sequences. The PCA (Figure 1A) showed that cheese samples were grouped according to the dairy plants, and no separation was observed between 5- and 15-day ripened cheese collected from the same dairy plant (Figure 1B). ASVs with a relative abundance of at least 0.5% in two samples are shown in Figure 2 and Supplementary Table 1. ASVs shared across different cheese type and dairy plant were analyzed by a Venn diagram analysis in order to detect unique ASVs or shared ones (Figure 3). In samples after 5 days of ripening, L. lactis, L. mesenteroides, Serratia sp., and the family Enterobacteriaceae were common to all dairy plants. Staphylococcus sciuri was unique of dairy C, while Ruminococcaceae were observed only in dairy plant A. After 15 days of ripening, Enterococcus sp. and Acinetobacter johnsonii were shared in the dataset, while Staphylococcus equorum was detected only in plant A, while Lactobacillus helveticus and Lactobacillus zeae were found mainly in plant D. Finally, S. sciuri and Ruminococcaceae were found unique in plant C (Figures 2, 3 and Supplementary Table 1).
FIGURE 1
FIGURE 2
FIGURE 3
Lactococcus lactis was the observed species with the highest relative abundance across all the samples, although it displayed some statistical differences between plants. Dairy E exhibit the highest relative frequency of L. lactis (from 96% to > 99%), followed by dairy plant C (from 91 to 99%), A and B showed similar relative abundance (from 76 to 99%), and D had the lowest relative abundance (from 53 to 84%). L. mesenteroides was found in all cheese samples with a relative abundance which reached 11%. In addition, Enterococcus spp. was observed as subdominant genus across all the samples.
Other taxa detected belonged to psychrotrophic genera such as Acinetobacter, Chryseobacterium, Corynebacterium, Pseudomonas, and Serratia. Additionally, few ASVs belonging to Staphylococcus sp. and the undesirable species S. aureus were observed in some samples at 1% of the relative abundance. The inferred metagenome was depicted in the heatplot, and two main cluster were identified (Supplementary Figure 1). Cheese samples collected from dairy plant A and B belong to cluster 1, dairy plant E and C belong to the cluster 2, while dairy D was quite spread between the two group. No specific distribution of cheeses at different ripening time was then observed. The 5-day cheeses of the dairy plant E and C showed a higher abundance of KEGG genes related to the lipid metabolism such as fatty acid and sphingolipid metabolism, biosynthesis of unsaturated fatty acids, and KEGG genes involved in citrate cycle (Supplementary Figure 1).
Mycobiota Composition of Robiola di Roccaverano Cheese
Five samples from the 26S rDNA sequencing were excluded from the analysis (two samples from dairy A; one sample from dairy B; one sample from dairy C; one sample from dairy D) due to poor DNA quality. The PCA analysis of the 26S data showed a cluster of cheese samples according to the plant (Figure 4A) independently of the time of ripening (Figure 4B). Venn diagram showed (Figure 3) 17 common ASVs in 5- and 15-day cheese. Galactomyces candidum and its telomorphic counterpart Geotrichum candidum (with a relative abundance accounting from 26 to 85%) and Kluyveromyces marxianus (with a relative abundance accounting from 4 to 49%) were the most abundant (Figure 5 and Supplementary Table 2). Several minor fraction of ASVs were in common in the dataset belonging to Yarrowia lipolytica, Debaryomyces hansenii, Trichosporon coremiiforme, Penicillium roqueforti, Kurtzmaniella anglica, Pichia fermentans, Geotrichum bryndzae, Kurtzmaniella santamariae, Wickerhamiella pararugosa, Cladosporium cladosporioides, Trichosporon, Galactomyces, and Yarrowia (for all details, see Supplementary Table 2). Some ASVs were detected in a specific dairy plant (Figures 3, 5). In 5-day ripened cheese, Candida parapsilosis, Rhodotorula mucilaginosa, Candida sake, and Fusarium spp. were detected in dairy A. Dairy plant C were characterized by the presence of Meyerozyma guilliermondii and Skvortzovia furfurella. Finally, some ASVs were found only in dairy plant E such as Actinomucor kuwaitiensis and Debaryomyces vindobonensis (Figures 3, 5). In 15-day ripened cheese, some differences have been observed. Dairy plant A was characterized by R. mucilaginosa and S. furfurella, while Torulospora delbrueckii (Figures 3, 5) was found only in dairy E.
FIGURE 4
FIGURE 5
Discussion
The present study aimed to characterize bacterial and fungal population of Robiola di Roccaverano cheese collected from five dairy plants. L. lactis was the dominant ASVs in all cheeses collected from the different plants probably due to the addition of natural ferments during the manufacturing process (). L. lactis is an important component of the defined and undefined starters worldwide () selected mainly for its ability to rapidly acidify milk during the cheese-making process. In addition to its starter activity, L. lactis has also been recognized for its ability to produce important aromatic compounds during cheese production (). In this study, samples from dairy plants E and C showed the highest frequency of L. lactis and the corresponding samples were mostly associated to inferred KEGG genes involved in lipid and citrate metabolism (Supplementary Figure 1). The principal compounds obtained through the citrate route are acetate, diacetyl, acetoin, and 2,3-butanediol, which has already been associated with L. lactis in miniaturized cheese (). In addition, it was reported that L. lactis displays the capability to metabolized carbohydrates through citrate route (). Strains belonging to the same species can have a different metabolic behavior; thus, we can hypothesize that the dairy environment selects specific strains that possess a strong lipolytic activity essential for cheese flavor development. The biological feature of the bacterial microbiota in this cheese can be due to the metabolic activity of specific L. lactis from dairy plant E and C. Among lactic acid bacteria, other ASVs were L. mesenteroides, Lb. helveticus, and Lcb. zeae. Lb. helveticus is commonly used as starter in Swiss-type cheese like Emmental but also to produce long-time ripened Italian cheeses like Provolone where its intense proteolytic activity is important during ripening process (). In contrast, the presence of Lcb. zeae has not been previously found in cheese made with NMC like Montasio and the so-called Traditional Specialty Guaranteed (TSG) Mozzarella cheese (; ).
Excluding plant E, which was probably more able to standardize the cheese production preventing the growth of undesirable species, the other dairies harbored several ASVs belonging to contaminant taxa. The presence of psychrotrophic bacteria like Pseudomonas spp., Acinetobacter spp., and Enterobacteriaceae, among which the genus Serratia, is common in raw milk microbiota (). Gram-negative bacteria are commonly present in milk samples as they are linked to animal and environmental factors and induce spoilage of dairy foods (; ; ). Although the results concerning the safety aspects could be considered promising, the viable bacterial flora was not evaluated in this study. The presence of toxin-producing Staphylococcus aureus, which is a common contaminant of raw milk and raw milk cheese, was not observed in Robiola di Roccaverano cheeses. Therefore, it could be postulated that the technology used in this study for the microbial characterization does not allow the detection of subdominant microbial population, among which Staphylococcus aureus.
Fungi possess an active role in the technological process related to the cheese production; however, they are also involved in spoilage phenomena. Indeed, some yeast species are recognized as primary drivers of defects during cheese making process, causing economic losses and leading to food safety issue (). Many studies have revealed that a great variety of fungal species are strictly associated with specific types of cheese, but few species tend to dominate the cheese systems (; ; ). In the present study, the core mycobiota consisted mainly of G. candidum followed by K. marxianus, which were found in high abundance in all the samples. G. candidum has been found in various environmental sources like milk and soft cheeses and recently also in fermented milk like Gioddu and Gwell (; ; ). The yeast K. marxianus could be a component of whey starter culture and NMC because of its ability to ferment lactose through β-galactosidase activity. Recently, this species has been isolated from the typical hard-cheese Fiore Sardo () where it influenced the ripening process through lactic acid utilization and by its proteolytic and lipolytic activity.
Other relevant yeasts observed in minor percentage were Y. lipolytica, D. hansenii, and T. coremiiforme. Debaryomyces spp. is frequently detected in milk and dairy products due to its resistance to low pH and high salt content ().
Yarrowia lipolytica represents one of the most common yeast species, but it cannot be added deliberately. In fact, it does not belong to the species commonly used to produce commercial starter culture used in the dairy industry. It is well known for its strong proteolytic and lipolytic activity which influence aroma and texture formation; thus, it is therefore particularly involved in ripening (). This yeast has already been detected by HTS in short-ripened cheese like Tomme d’Orchies ().
Some detected ASVs were characteristic of specific dairy plant such as P. roqueforti in dairy D. This mold is often present in artisanal blue or white mold cheeses (). It is known to possess the ability to produce volatile aromatic compound due to its strong lipolytic and proteolytic activity, which could have a positive impact on flavor formation during blue-cheese production (). Moreover, Saturnispora silvae, particularly detected in dairy E, was already isolated and identified from Robiola di Roccaverano cheese by culture-dependent methods and could be ascribed to an environmental source of contamination (, ). Finally, Kurtzmaniella anglica and Trichosporon coremiiforme, which were found mainly in dairy plants A, B, and E, respectively, were already detected in a previous work on Robiola di Roccaverano cheese ().
Some of the genera detected in low abundance in this study, namely, Penicillium, Cladosporium, and Fusarium, are undesirable fungi responsible for spoilage. Among them, Penicillium spp. constituted the largest group of spoilage molds in food and some species could be responsible for the production of mycotoxin (). Cladosporium spp. growing poorly is considered quite common in air system (). Mold contamination through ambient air is normal, more than for the yeast. Therefore, it is of notable importance, especially for small producers of artisanal cheeses such as Robiola di Roccaverano diaries, to follow accurate hygienic practice during the manufacturing and packaging to prevent the contamination.
Conclusion
In this study, metataxonomic analysis allowed the characterization of the microbial and fungal population of cheese samples collected from different Robiola di Roccaverano diaries. The differences detected suggested that, despite the localization of the dairy plants in different areas and the raw materials used, i.e., raw goat’s milk and NMC, it is possible to define a core microbial population constituted by L. lactis, G. candidum, and K. marxianus. However, some dairy-specific microbial and fungal traits have been identified and may help in the definition of specific hallmark.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
Sequencing data were deposited in the Sequence Read Archive of the National Center for Biotechnology Information (NCBI; https://www.ncbi.nlm.nih.gov/) under Bioproject accession numbers PRJNA704571 and PRJNA626588.
Author contributions
FB: sample collection and preparation, and manuscript drafting. FB and IF: sequencing, data analysis, and interpretation. MB and AD: project administration, conceptualization, study design, and funding acquisition. All authors contributed to manuscript revision and read and approved the submitted version.
Funding
This work was supported by the University of Turin (Italy; 60% grant).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.776862/full#supplementary-material
Supplementary Figure 1Heatplot of the abundance of KEGG gene belonging to carbohydrates, amino acid, and lipid metabolism pathways (left bar) in 5- and 15-day ripened cheese (upper bar from the top) and in the different dairy plant (lower bar from the top). The intensity of the color represents the degree of correlation between samples and inferred metabolic pathways.
References
1
AnelliP.HaidukowskiM.EpifaniF.CimmarustiM. T.MorettiA.LogriecoA.et al (2019). Fungal mycobiota and mycotoxin risk for traditional artisan Italian cave cheese.Food Microbiol.7862–72. 10.1016/j.fm.2018.09.014
2
BiolcatiF.AndrighettoC.BotteroM. T.DalmassoA. (2020). Microbial characterization of an artisanal production of Robiola di Roccaverano cheese.J. Dairy Sci.1034056–4067. 10.3168/jds.2019-17451
3
BiolcatiF.BotteroM. T.DalmassoA. (2019). Microbiological analysis of the Robiola di Roccaverano cheese by means of traditional culturing methods.Ital. J. Food Saf.8209–212. 10.4081/ijfs.2019.8574
4
BiolcatiF.FerrocinoI.BotteroM. T.DalmassoA. (2021). Mycobiota composition of Robiola di Roccaverano cheese along the production chain.Foods10:1859. 10.3390/foods10081859
5
BokulichN. A.AmiranashviliL.ChitchyanK.GhazanchyanN.DarbinyanK.GagelidzeN.et al (2015). Microbial biogeography of the transnational fermented milk matsoni.Food Microbiol.5012–19. 10.1016/j.fm.2015.01.018
6
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2.Nat. Biotechnol.37852–857. 10.1038/s41587-019-0209-9
7
BoutrouR.GuéguenM. (2005). Interests in Geotrichum candidum for cheese technology.Int. J. Food Microbiol.1021–20. 10.1016/j.ijfoodmicro.2004.12.028
8
CallahanB. J.McMurdieP. J.RosenM. J.HanA. W.JohnsonA. J. A.HolmesS. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data.Nat. Methods13581–583. 10.1038/nmeth.3869
9
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). Correspondence QIIME allows analysis of high- throughput community sequencing data Intensity normalization improves color calling in SOLiD sequencing.Nat. Publ. Gr.7335–336. 10.1038/nmeth0510-335
10
CaronT.PiverM.Le PéronA. C.LiebenP.LavigneR.BrunelS.et al (2021). Strong effect of Penicillium roqueforti populations on volatile and metabolic compounds responsible for aromas, flavor and texture in blue cheeses.Int. J. Food Microbiol.354:109174. 10.1016/j.ijfoodmicro.2021.109174
11
CeugniezA.TaminiauB.CoucheneyF.JacquesP.DelcenserieV.DaubeG.et al (2017). Fungal diversity of “Tomme d’Orchies” cheese during the ripening process as revealed by a metagenomic study.Int. J. Food Microbiol.25889–93. 10.1016/j.ijfoodmicro.2017.07.015
12
DixonP. (2003). Computer program review VEGAN, a package of R functions for community ecology.J. Veg. Sci.14927–930. 10.1111/j.1654-1103.2003.tb02228.x
13
EdgarR. C.HaasB. J.ClementeJ. C.QuinceC.KnightR. (2011). UCHIME improves sensitivity and speed of chimera detection.Bioinformatics272194–2200. 10.1093/bioinformatics/btr381
14
ErcoliniD.RussoF.FerrocinoI.VillaniF. (2009). Molecular identification of mesophilic and psychrotrophic bacteria from raw cow’s milk.Food Microbiol.26228–231. 10.1016/j.fm.2008.09.005
15
ErkusO.De JagerV. C. L.SpusM.Van Alen-BoerrigterI. J.Van RijswijckI. M. H.HazelwoodL.et al (2013). Multifactorial diversity sustains microbial community stability.ISME J.72126–2136. 10.1038/ismej.2013.108
16
FaddaM. E.MossaV.DeplanoM.PisanoM. B.CosentinoS. (2017). In vitro screening of Kluyveromyces strains isolated from Fiore Sardo cheese for potential use as probiotics.LWT Food Sci. Technol.75100–106. 10.1016/j.lwt.2016.08.020
17
FerrocinoI.GreppiA.La StoriaA.RantsiouK.ErcoliniD.CocolinL. (2016). Impact of nisin-activated packaging on microbiota of beef burgers during storage.Appl. Environ. Microbiol.82549–559. 10.1128/AEM.03093-15
18
FerrocinoI.RantsiouK.CocolinL.MicrobiologyA.SectorF. T. (2020). Microbiota of Milk and Dairy Foods?: Structure and Function by -omics Approaches.Amsterdam: Elsevier, 10.1016/B978-0-08-100596-5.22973-9
19
GarnierL.ValenceF.MounierJ. (2017). Diversity and control of spoilage fungi in dairy products: an update.Microorganisms5:42. 10.3390/microorganisms5030042
20
Gazzetta Ufficiale della Repubblica Italiana (2013). Della Repubblica Italiana, Vol. 154. 52–54. Available online at: https://www.gazzettaufficiale.it/eli/gu/2013/07/10/160/sg/pdf(accessed August 30, 2021).
21
GriffithsM. W.TellezA. M. (2013). Lactobacillus helveticus: the proteolytic system.Front. Microbiol.4:30. 10.3389/fmicb.2013.00030
22
GuidoneA.RicciardiA.RomanielloA.BonomoM. G.MoroneG.ZottaT.et al (2016). Microbial changes of natural milk cultures for mozzarella cheese during repeated propagation cycles.LWT Food Sci. Technol.65572–579. 10.1016/j.lwt.2015.08.031
23
Hantsis-ZacharovE.ShakeT.HalpernM. (2008). Chryseobacterium oranimense sp . nov ., a psychrotolerant, proteolytic and lipolytic bacterium isolated from raw cow’ s milk.Int. J. Syst. Evol. Microbiol.42635–2639. 10.1099/ijs.0.65819-0
24
KlindworthA.PruesseE.SchweerT.PepliesJ.QuastC.HornM.et al (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies.Nucleic Acids Res.411–11. 10.1093/nar/gks808
25
KureC. F.SkaarI. (2019). The fungal problem in cheese industry.Curr. Opin. Food Sci.2914–19. 10.1016/j.cofs.2019.07.003
26
LangilleM. G. I.ZaneveldJ.CaporasoJ. G.McDonaldD.KnightsD.ReyesJ. A.et al (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences.Nat. Biotechnol.31814–821. 10.1038/nbt.2676
27
MagočT.SalzbergS. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies.Bioinformatics272957–2963. 10.1093/bioinformatics/btr507
28
MaoloniA.BlaiottaG.FerrocinoI.MangiaN. P.OsimaniA.MilanovićV.et al (2020). Microbiological characterization of Gioddu, an Italian fermented milk.Int. J. Food Microbiol.323:108610. 10.1016/j.ijfoodmicro.2020.108610
29
MarínP.PalmeroD.JuradoM. (2015). Occurrence of moulds associated with ovine raw milk and cheeses of the Spanish region of Castilla La Mancha.Int. J. Dairy Technol.68565–572. 10.1111/1471-0307.12208
30
MarinoM.MaifreniM.RondininiG. (2003). Microbiological characterization of artisanal montasio cheese: analysis of its indigenous lactic acid bacteria.FEMS Microbiol. Lett.229133–140. 10.1016/S0378-1097(03)00816-4
31
McSweeneyP. L. H.SousaM. J. (2000). Biochemical pathways for the production of flavour compounds in cheeses during ripening: a review.Lait80293–324. 10.1051/lait:2000127
32
MontelM. C.BuchinS.MalletA.Delbes-PausC.VuittonD. A.DesmasuresN.et al (2014). Traditional cheeses: rich and diverse microbiota with associated benefits.Int. J. Food Microbiol.177136–154. 10.1016/j.ijfoodmicro.2014.02.019
33
Mota-GutierrezJ.FerrocinoI.RantsiouK.CocolinL. (2019). Metataxonomic comparison between internal transcribed spacer and 26S ribosomal large subunit (LSU) rDNA gene.Int. J. Food Microbiol.290132–140. 10.1016/j.ijfoodmicro.2018.10.010
34
RuggirelloM.GiordanoM.BertolinoM.FerrocinoI.CocolinL.DolciP. (2018). Study of Lactococcus lactis during advanced ripening stages of model cheeses characterized by GC-MS.Food Microbiol.74132–142. 10.1016/j.fm.2018.03.012
35
SuzziG.LanorteM. T.GalganoF.AndrighettoC.LombardiA.LanciottiR.et al (2001). Proteolytic, lipolytic and molecular characterisation of Yarrowia lipolytica isolated from cheese.Int. J. Food Microbiol.6969–77. 10.1016/S0168-1605(01)00574-8
36
von GastrowL.MadecM. N.ChuatV.LubacS.MorinièreC.LéS.et al (2020). Microbial diversity associated with Gwell, a traditional French mesophilic fermented milk inoculated with a natural starter.Microorganisms81–21. 10.3390/microorganisms8070982
37
WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy.Appl. Environ. Microbiol.735261–5267. 10.1128/AEM.00062-07
38
WoutersJ. T. M.AyadE. H. E.HugenholtzJ.SmitG. (2002). Microbes from raw milk for fermented dairy products.Int. Dairy J.1291–109. 10.1016/S0958-6946(01)00151-0
Summary
Keywords
cheese, mycobiota, metagenetic high-throughput sequencing, Robiola di Roccaverano, bacterial microbiota
Citation
Biolcati F, Ferrocino I, Bottero MT and Dalmasso A (2022) The Bacterial and Fungal Microbiota of “Robiola di Roccaverano” Protected Designation of Origin Raw Milk Cheese. Front. Microbiol. 12:776862. doi: 10.3389/fmicb.2021.776862
Received
14 September 2021
Accepted
21 December 2021
Published
31 January 2022
Volume
12 - 2021
Edited by
Rosanna Tofalo, University of Teramo, Italy
Reviewed by
Anne Thierry, INRA UMR 1253 Science & Technologie du Lait & de l’Oeuf, France; Giacomo Zara, University of Sassari, Italy
Updates
Copyright
© 2022 Biolcati, Ferrocino, Bottero and Dalmasso.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Federica Biolcati, federica.biolcati@unito.it
This article was submitted to Food Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.