Microbiological and Metagenomic Analysis to Assess the Effect of Container Material on the Microbiota of Feta Cheese during Ripening

Investigation of the Microbiome of Industrial PDO Sfela Cheese and Its Artisanal Variants Using 16S rDNA Amplicon Sequencing and Shotgun Metagenomics

Sfela is a white brined Greek cheese of protected designation of origin (PDO) produced in the Peloponnese region from ovine, caprine milk, or a mixture of the two. Despite the PDO status of Sfela, very few studies have addressed its properties, including its microbiology. For this reason, we decided to investigate the microbiome of two PDO industrial Sfela cheese samples along with two non-PDO variants, namely Sfela touloumotiri and Xerosfeli. Matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS), 16S rDNA amplicon sequencing and shotgun metagenomics analysis were used to identify the microbiome of these traditional cheeses. Cultured-based analysis showed that the most frequent species that could be isolated from Sfela cheese were Enterococcus faecium, Lactiplantibacillus plantarum, Levilactobacillus brevis, Pediococcus pentosaceus and Streptococcus thermophilus. Shotgun analysis suggested that in industrial Sfela 1, Str. thermophilus dominated, while industrial Sfela 2 contained high levels of Lactococcus lactis. The two artisanal samples, Sfela touloumotiri and Xerosfeli, were dominated by Tetragenococcus halophilus and Str. thermophilus, respectively. Debaryomyces hansenii was the only yeast species with abundance > 1% present exclusively in the Sfela touloumotiri sample. Identifying additional yeast species in the shotgun data was challenging, possibly due to their low abundance. Sfela cheese appears to contain a rather complex microbial ecosystem and thus needs to be further studied and understood. This might be crucial for improving and standardizing both its production and safety measures.

Pilot-Scale Production of Traditional Galotyri PDO Cheese from Boiled Ewes’ Milk Fermented with the Aid of Greek Indigenous Lactococcus lactis subsp. cremoris Starter and Lactiplantibacillus plantarum Adjunct Strains

The performance of a mixed thermophilic and mesophilic starter culture consisting of Streptococcus thermophilus ST1 and the Greek indigenous nisin-A-producing Lactococcus lactis subsp. cremoris M78 was evaluated in the absence (A: ST1+M78) or presence (B: ST1+M78+H25) of Lactiplantibacillus plantarum H25—another indigenous ripening strain—under real cheesemaking conditions. Three pilot-scale trials of fresh (6-day-old) Galotyri PDO cheese were made from boiled milk by an artisanal method using simple equipment, followed by cold ripening of the A1–A3 and B1–B3 cheeses at 4 °C for 30 days. All of the cheeses were analyzed microbiologically and for pH, gross composition, proteolysis, sugar and organic acid contents, and sensorial attributes before and after ripening. The artisanal (PDO) Galotyri manufacturing method did not ensure optimal growth of the ST1+M78 starter as regards the constant ability of the thermophilic strain ST1 to act as the primary milk acidifier under ambient (20–30 °C) fermentation conditions. Consequently, major trial-dependent microbial and biochemical differences between the Acheeses, and generally extended to the Bcheeses, were found. However, high-quality Galotyri was produced when either starter strain predominated in the fresh cheeses; only trial A1 had microbiological and sensory defects due to an outgrowth of post-thermal Gram-negative bacterial contaminants in the acidified curd. The H25 adjunct strain, which grew above 7 to 9 log CFU/g depending on the trial, had minor effects on the cheese’s pH, gross composition, and proteolysis, but it improved the texture, flavor, and the bacteriological quality of the Bcheeses during processing, and it exerted antifungal effects in the ripened cheeses.

An Insight into Goat Cheese: The Tales of Artisanal and Industrial Gidotyri Microbiota

The purpose of this study was to determine for the first time the microbiota in artisanal-type and industrial-type Gidotyri cheeses and investigate the influence of the cheese-making practices on their composition using culture-independent techniques. The microbiota present in artisanal with commercial starters (Artisanal_CS, n = 15), artisanal with in-house starters (Artisanal_IHS, n = 10) and industrial (Ind., n = 9) Gidotyri cheese samples were analyzed using a targeted metagenomic approach (16S rRNA gene). The Ind. Gidotyri cheese microbiota were less complex, dominated by the Streptococcaceae family (91%) that was more abundant compared to the artisanal Gidotyri cheeses (p < 0.05). Artisanal cheeses were more diverse compositionally with specific bacterial species being prevalent to each subtype. Particularly, Loigolactobacillus coryniformis (OTU 175), Secundilactobacillus malefermentans (OTU 48), and Streptococcus parauberis (OTU 50) were more prevalent in Artisanal_IHS cheeses compared to Artisanal_CS (p ≤ 0.001) and Ind. (p < 0.01) Gidotyri cheeses. Carnobacterium maltaromaticum (OTU 23) and Enterobacter hormaechei subsp. hoffmannii (OTU 268) were more prevalent in Artisanal_CS cheeses compared to Artisanal_IHS cheeses (p < 0.05) and Ind. cheeses (p < 0.05). Hafnia alvei (OTU 13) and Acinetobacter colistiniresistens (OTU 111) tended to be more prevalent in Artisanal_CS compared to the other two cheese groups (p < 0.10). In conclusion, higher microbial diversity was observed in the artisanal-type Gidotyri cheeses, with possible bacterial markers specific to each subtype identified with potential application to traceability of the manufacturing processes’ authenticity and cheese quality.

The Impact of Different Inoculation Schemes on the Microbiota, Physicochemical and Sensory Characteristics of Greek Kopanisti Cheese throughout Production and Ripening

Kopanisti is a Greek PDO cheese, which is traditionally produced by the addition of an amount of over-mature Kopanisti, called Mana Kopanisti, to initiate cheese ripening. The aim of this study was the production of four types of Kopanisti cheese (A-D) using pasteurized cow milk, and a combination of the following starters/adjuncts in order to test their ability to be used in Kopanisti cheese production: A: Lactococcus lactis subsp. lactis and Lacticaseibacillus paracasei, B: L. lactis and Lc. paracasei/Mana Kopanisti, C: L. lactis and Lc. paracasei/Ligilactobacillus acidipiscis and Loigolactobacillus rennini, D: Lig. acidipiscis and Loig. rennini. Throughout production and ripening, classical microbiological, metataxonomics and physicochemical analyses were employed, while the final products (Day 35) were subjected to sensory analysis as well. Most interestingly, beta-diversity analysis of the metataxonomics data revealed the clusters constructed among the Kopanisti types based on the different inoculation schemes. On day 35, Kopanisti A-C types clustered together due to their similar 16S microbiota, while Kopanisti D was highly differentiated. On the contrary, ITS data clustered Kopanisti B and C together, while Kopanisti A and D were grouped seperately. Finally, based on the sensory evaluation, Kopanisti C appeared to have the most suitable bacteria cocktail for the Kopanisti cheese production. Therefore, not only were the conventional starters used, but also the Lig. acidipiscis and Loig. rennini strains could be used in a standardized Kopanisti cheese production that could lead to final products of high quality and safety.

Microbiological Characterization of Greek Galotyri Cheese PDO Products Relative to Whether They Are Marketed Fresh or Ripened

Galotyri is the most popular traditional Greek PDO soft acid-curd cheese. This study compared the microbial numbers and types and characterized the lactic acid bacteria (LAB) biota of two artisan-type Galotyri PDO cheese varieties, one marketed fresh (Brand-K) and the other ripened (Brand-Z). Two retail batches of each cheese variety were analyzed, and a total of 102 LAB isolates were biochemically identified. LAB (7.2–9.3 log CFU/g) prevailed in all cheeses, followed by yeasts (5.8–6.8 log CFU/g). Typical starter strains of Streptococcus thermophilus and Lactobacillus delbrueckii were the most abundant species in all batches. However, the fresh Brand-K cheeses had 1–3 log units higher thermophilic starter LAB counts than the ripened Brand-Z cheeses, which contained a more diverse viable LAB biota comprising Lacticaseibacillus paracasei, Leuconostocmesenteroides, Lentilactobacillus (L. diolivorans, L. kefiri, L. hilgardii), Pediococcusinopinatus/parvulus, few spontaneous nonstarter thermophilic streptococci and lactobacilli, and Enterococcus faecium and E. faecalis at higher subdominant levels.Conversely, the fresh Brand-K cheeses were enriched in members of the Lactiplantibacillus plantarum group; other LAB species were sporadically isolated, including Lactococcus lactis. All retail cheeses were safe (pH 3.9–4.0). No Salmonella spp. or Listeria monocytogenes were detected in 25-g samples by culture enrichment; however, Listeria innocua and coagulase-positive staphylococci (850 CFU/g) survived in one ripened batch. Gram-negative bacteria were <100 CFU/g in all cheeses. In conclusion, ripening reduced the starter LAB viability but increased the nonstarter LAB species diversity in the present Galotyri PDO market cheeses.

Metataxonomic Mapping of the Microbial Diversity of Irish and Eastern Mediterranean Cheeses

The distinct sensorial characteristics of local cheeses influence consumer preferences, and make an essential contribution to the local economy. Microbial diversity in cheese is among the fundamental contributors to sensorial and qualitative characteristics. However, knowledge regarding the existence of microbial patterns associated with regional production practices in ripened cheeses remains limited. The present research was conducted to test the hypothesis that the background metagenome of cheeses could be used as a marker of their origin. We compared Irish versus Eastern Mediterranean cheeses-namely Greek and Cypriot-using High Throughput Sequencing (HTS). The study identified a significantly distinct separation among cheeses originating from the three different countries, in terms of the total microbial community composition. The use of machine learning and biomarkers discovery algorithms defined key microbes that differentiate each geographic region. Finally, the development of interaction networks revealed that the key species developed mostly negative interactions with the other members of the communities, highlighting their dominance in the community. The findings of the present research demonstrate that metagenome could indeed be used as a biological marker of the origin of mature cheeses, and could provide further insight into the dynamics of microbial community composition in ripened cheeses.

Metataxonomic sequencing to assess microbial safety of Turkish white cheeses

High-throughput sequencing has provided a way to monitor the large diversity of microorganisms in fermented foods that have complex microbiota. Up to date, many kinds of cheese have been characterized with the metataxonomic approach, but the safety of unpacked Turkish white cheeses, which are widely consumed in Turkey, has not been assessed. In this study, fifteen unpacked white cheeses sold in public bazaars in Ankara province have been collected and subjected to microbial enumeration as well as physicochemical analysis. Five white cheeses, which have relatively the highest foodborne pathogens, out of fifteen white cheeses, have been analyzed by next-generation sequencing and metataxonomic analysis. According to the results, abundant families were Lactobacillaceae, Oceanospirillaceae, Enterococcaceae, Pseudomonadaceae, and Vibrionaceae. Staphylococcus aureus, E. coli, and Salmonella, which are indicators of bad hygiene and sanitation conditions, were found in cheeses. In conclusion, culture-independent methods such as metataxonomic can be important to evaluate the safety of foods.

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.

Distinctive Traits of Four Apulian Traditional Agri-Food Product (TAP) Cheeses Manufactured at the Same Dairy Plant

This study aimed to highlight the distinctive features of four Traditional Agri-food Products (TAP), namely, Caprino, Pecorino, Vaccino, and Cacioricotta cheeses produced at the same dairy plant to reveal any possible relationships between their microbiological and biochemical characteristics. Two distinct natural whey starter (NWS) cultures were used during Caprino and Vaccino cheesemaking, whereas no starter was used for the other cheeses. Cacioricotta retained the highest concentrations of salt and residual carbohydrates. Lactic acid bacteria dominated the microbiota of the cheeses. Furthermore, staphylococci represented an additional dominant microbial population in Cacioricotta. Although culture-dependent analysis showed that the use of NWS cultures only slightly affected the microbial community of cheeses, 16S metagenetic analysis showed that Lactobacillus helveticus dominated both the NWS cultures and the corresponding Caprino and Vaccino cheeses. This analysis indicated that Staphylococcus equorum and Streptococcus thermophilus dominated Cacioricotta and Pecorino cheeses, respectively. The highest peptidase activities were found in either Caprino or Vaccino. Enzymes involved in the catabolism of free amino acids and esterase showed the highest activity in Pecorino cheese. Each cheese showed a distinct profile of volatile organic compounds, with Pecorino being the richest cheese in carboxylic acids, ketones, and esters, related to lipolysis. The results of this study contribute to valorizing and safeguarding these TAP cheeses, sustaining local farming.

Characterization of Bacterial Microbiota of P.D.O. Feta Cheese by 16S Metagenomic Analysis

Background: The identification of bacterial species in fermented PDO (protected designation of origin) cheese is important since they contribute significantly to the final organoleptic properties, the ripening process, the shelf life, the safety and the overall quality of cheese. Methods: Ten commercial PDO feta cheeses from two geographic regions of Greece, Epirus and Thessaly, were analyzed by 16S metagenomic analysis. Results: The biodiversity of all the tested feta cheese samples consisted of five phyla, 17 families, 38 genera and 59 bacterial species. The dominant phylum identified was Firmicutes (49% of the species), followed by Proteobacteria (39% of the species), Bacteroidetes (7% of the species), Actinobacteria (4% of the species) and Tenericutes (1% of the species). Streptococcaceae and Lactobacillaceae were the most abundant families, in which starter cultures of lactic acid bacteria (LAB) belonged, but also 21 nonstarter lactic acid bacteria (NSLAB) were identified. Both geographical areas showed a distinctive microbiota fingerprint, which was ultimately overlapped by the application of starter cultures. In the rare biosphere of the feta cheese, Zobellella taiwanensis and Vibrio diazotrophicus, two Gram-negative bacteria which were not previously reported in dairy samples, were identified. Conclusions: The application of high-throughput DNA sequencing may provide a detailed microbial profile of commercial feta cheese produced with pasteurized milk.

Serra da Estrela PDO Cheese Microbiome as Revealed by Next Generation Sequencing

Serra da Estrela PDO cheese is the oldest traditional cheese manufactured in Portugal. In this work, its microbiome as well as the main raw materials used in cheese production, raw ewes’ milk and thistle flowers (Cynara cardunculus L.), were characterized using next generation sequencing. Samples were accordingly retrieved from a local producer over two consecutive production campaigns and at different time periods within each campaign. The bacterial and fungi communities associated with each matrix were accessed through sequencing of V3–V4 and Internal Transcribed Spacer 2 regions of rRNA gene amplicons, respectively. A high microbial diversity was found associated to each matrix, differing significantly (p < 0.05) from each other. Over 500 taxa were identified in each analyzed matrix, ranging from dominant (relative abundance > 1%), sub-dominant (0.01–1%) and rare taxa (<0.01%). Specifically, in cheese, 30 taxa were present in all analyzed samples (core taxa), including species of Leuconostoc spp. and Lactococcus spp. for bacteria and Candida spp., Debaryomyces spp. and Yarrowia spp. for fungi, that were cumulatively the most prevalent genera in Serra da Estrela PDO cheese (average relative abundance ≥10%). Ultimately, this characterization study may contribute to a better understanding of the microbial dynamics of this traditional PDO product, namely the influence of raw materials on cheese microbiome, and could assist producers interested in preserving the identity, quality and safety of Serra da Estrela PDO cheese.

Data Science in the Food Industry

Food safety is one of the main challenges of the agri-food industry that is expected to be addressed in the current environment of tremendous technological progress, where consumers' lifestyles and preferences are in a constant state of flux. Food chain transparency and trust are drivers for food integrity control and for improvements in efficiency and economic growth. Similarly, the circular economy has great potential to reduce wastage and improve the efficiency of operations in multi-stakeholder ecosystems. Throughout the food chain cycle, all food commodities are exposed to multiple hazards, resulting in a high likelihood of contamination. Such biological or chemical hazards may be naturally present at any stage of food production, whether accidentally introduced or fraudulently imposed, risking consumers' health and their faith in the food industry. Nowadays, a massive amount of data is generated, not only from the next generation of food safety monitoring systems and along the entire food chain (primary production included) but also from the Internet of things, media, and other devices. These data should be used for the benefit of society, and the scientific field of data science should be a vital player in helping to make this possible.

Microbial profiles of Greek PDO cheeses assessed with amplicon metabarcoding

Greece is a country possessing many cheese products granted with a PDO (Protected Designation of Origin) certificate, with high exporting activities. In this study, we analyzed six popular cheese PDO products purchased from different industries to assess their microbial communities using amplicon metabarcoding analysis. To this end, using Next Generation Sequencing technology, we sequenced the 16S rRNA gene and the ITS spacer for prokaryotes and fungi, respectively. Alpha diversity indices revealed higher bacterial species richness for some cheeses (Kopanisti, Batzos) and poor for others (Feta, Galotiri). Kopanisti, together with Kalathaki and Anevato, also presented increased species diversity concerning fungal populations. Results showed that lactic acid bacteria (LAB) prevailed the bacterial populations in all samples (Lactococcus, Lactobacillus, Streptococcus, Leuconostoc), whereas for fungi, members of the Saccharomycetaceae, Dipodascaceae and Debaryomycetaceae families prevailed the fungal populations. Several other genera were identified that make up each product's microbiome leading to the creation of the unique organoleptic attributes of Greek PDO cheeses. However, the identified species could not be directly linked to certain cheese types, assuming that starter and adjunct cultures, combined with the raw material used during production greatly impact the microbial communities in cheeses. Our data, produced for the first time for six Greek PDO cheeses, can be exploited in the process of creating a core microbial signature within each cheese type, supporting the Greek brand name and valorizing cheese products.

Microbiological and Metagenomic Characterization of a Retail Delicatessen Galotyri-Like Fresh Acid-Curd Cheese Product

This study evaluated the microbial quality, safety, and ecology of a retail delicatessen Galotyri-like fresh acid-curd cheese traditionally produced by mixing fresh natural Greek yogurt with ‘Myzithrenio’, a naturally fermented and ripened whey cheese variety. Five retail cheese batches (mean pH 4.1) were analyzed for total and selective microbial counts, and 150 presumptive isolates of lactic acid bacteria (LAB) were characterized biochemically. Additionally, the most and the least diversified batches were subjected to a culture-independent 16S rRNA gene sequencing analysis. LAB prevailed in all cheeses followed by yeasts. Enterobacteria, pseudomonads, and staphylococci were present as <100 viable cells/g of cheese. The yogurt starters Streptococcus thermophilus and Lactobacillus delbrueckii were the most abundant LAB isolates, followed by nonstarter strains of Lactiplantibacillus, Lacticaseibacillus, Enterococcus faecium, E. faecalis, and Leuconostoc mesenteroides, whose isolation frequency was batch-dependent. Lactococcus lactis isolates were sporadic, except for one cheese batch. However, Lactococcus lactis, Enterobacteriaceae, Vibrionaceae, Salinivibrio, and Shewanellaceae were detected at fairly high relative abundances culture-independently, despite the fact that their viable counts in the cheeses were low or undetectable. Metagenomics confirmed the prevalence of S. thermophilus and Lb. delbrueckii. Overall, this delicatessen Galotyri-like cheese product was shown to be a rich pool of indigenous nonstarter LAB strains, which deserve further biotechnological investigation.

Microbiota of Chicken Breast and Thigh Fillets Stored under Different Refrigeration Temperatures Assessed by Next-Generation Sequencing

Chicken is one of the most widely consumed meats worldwide. The exploration of the bacterial diversity of chicken meat may provide new insights into the chicken-associated microbiome that will lead to moderation of food spoilage or safety. This study was undertaken to explore the bacterial communities of chicken breast and thigh fillets stored at refrigeration (0 °C and 5 °C) and slightly abuse (10 °C) temperatures for 5 days through conventional cultural methods along with next-generation sequencing (NGS) analysis. Total viable counts (TVC), Brochothrix thermosphacta, Pseudomonas spp., and lactic acid bacteria (LAB) were enumerated, while the bacterial communities were mapped through 16S rRNA gene amplicon sequencing. Chicken breast and thigh fillets possessed a complex bacterial structure that incorporated a total of >200 Operational Taxonomic Units (OTUs) at the genus level. The core microbiota of fresh samples consisted of Acinetobacter, Brochothrix, Flavobacterium, Pseudomonas, Psychrobacter, and Vibrionaceae (family). These genera persisted until the end of storage in >80% of samples, except Psychrobacter and Flavobacterium, while Photobacterium was also identified. Hierarchical clustering showed a distinction of samples based on storage time and chicken part. Conventional plate counting with growth media commonly used in spoilage studies did not always correspond to the microbial community profiles derived from NGS analysis, especially in Pseudomonas, Acinetobacter, Photobacterium, and Vibrionaceae. Results of the present study highlight Photobacterium and Vibrionaceae, in general, as potent chicken meat spoilers and suggest the necessity to combine classical microbiological methods along with NGS technologies to characterize chicken meat spoilage microbiota.

Exploring the Bacterial Communities of the Kaiafas Thermal Spring Anigrides Nymphes in Greece Prior to Rehabilitation Actions

Anigrides Nymphes of Lake Kaiafas is a thermal spring that is well known for its therapeutical properties, as the hot water (32-34 °C) is rich in sulfur compounds and minerals. Nowadays, efforts are made from the Hellenic Republic to modernize the existing facilities and infrastructure networks of the area. To study the complex ecosystem of the thermal spring, we collected water from four sampling points (Lake, and Caves 1, 2, and 3). Filtration method was used for microbial enumeration. In parallel, total bacterial DNA was extracted and subjected to next-generation sequencing (NGS). A total of 166 different bacterial families were detected. Differences in families, genera, and species abundances were detected between the different sampling points. Specifically, Comamonadaceae was the most common family detected in Lake and Cave 3. Similarly, in Caves 1 and 2, Rhodobacteraceae was detected at a higher percentage compared to the rest of the families. Moreover, the detection of sequences assigned to waterborne or opportunistic pathogens, i.e., Enterobacteriaceae, Legionellaceae, Coxiellaceae, and Clostridiaceae, as well as Enterococcus and Vibrio, is of great importance. Although the presence of pathogens was not examined by quantitative PCR, the detection of their sequences strengthens the need of the planned rehabilitation actions of this natural environment in order to allow human swimming.