Comparison of the Microbiome of Artisanal Homemade and Industrial Feta Cheese through Amplicon Sequencing and Shotgun Metagenomics

Feta is the most renowned protected designation of origin (PDO) white brined cheese produced in Greece. The fine organoleptic characteristics and the quality of Feta rely on, among other factors, its overall microbial ecosystem. In this study, we employed 16S rDNA and internal transcribed spacer (ITS) amplicon sequencing, as well as shotgun metagenomics, to investigate the microbiome of artisanal homemade and industrial Feta cheese samples from different regions of Greece, which has very rarely been investigated. 16S rDNA data suggested the prevalence of the Lactococcus genus in the homemade samples, while Streptococcus and Lactobacillus genera prevailed in the industrial control samples. Species identification deriving from shotgun metagenomics corroborated these findings, as Lactococcus lactis dominated two homemade samples while Streptococcus thermophilus and Lactobacillus delbrueckii subsp. bulgaricus were found to be dominating one industrial sample. ITS data revealed a complex diversity of the yeast population among the samples analyzed. Debaryomyces, Kluyveromyces, Cutaneotrichosporon, Pichia, Candida, and Rhodotorula were the major genera identified, which were distributed in a rather arbitrary manner among the different samples. Furthermore, a number of potential metagenome-assembled genomes (MAGs) could be detected among assembled shotgun bins. The overall analysis of the shotgun metagenomics supported the presence of different foodborne pathogens in homemade samples (e.g., Staphylococcus aureus, Listeria monocytogenes, Enterobacter cloacae, and Streptococcus suis), but with low to very low abundances. Concluding, the combination of both amplicon sequencing and shotgun metagenomics allowed us to obtain an in-depth profile of the artisanal homemade Feta cheese microbiome.