Evaluation of bacterial communities of Grana Padano cheese by DNA metabarcoding and DNA fingerprinting analysis

Cultivable microbial diversity, peptide profiles, and bio-functional properties in Parmigiano Reggiano cheese

Introduction

Lactic acid bacteria (LAB) communities shape the sensorial and functional properties of artisanal hard-cooked and long-ripened cheeses made with raw bovine milk like Parmigiano Reggiano (PR) cheese. While patterns of microbial evolution have been well studied in PR cheese, there is a lack of information about how this microbial diversity affects the metabolic and functional properties of PR cheese.

Methods

To fill this information gap, we characterized the cultivable fraction of natural whey starter (NWS) and PR cheeses at different ripening times, both at the species and strain level, and investigated the possible correlation between microbial composition and the evolution of peptide profiles over cheese ripening.

Results and discussion

The results showed that NWS was a complex community of several biotypes belonging to a few species, namely, Streptococcus thermophilus, Lactobacillus helveticus, and Lactobacillus delbrueckii subsp. lactis. A new species-specific PCR assay was successful in discriminating the cheese-associated species Lacticaseibacillus casei, Lacticaseibacillus paracasei, Lacticaseibacillus rhamnosus, and Lacticaseibacillus zeae. Based on the resolved patterns of species and biotype distribution, Lcb. paracasei and Lcb. zeae were most frequently isolated after 24 and 30 months of ripening, while the number of biotypes was inversely related to the ripening time. Peptidomics analysis revealed more than 520 peptides in cheese samples. To the best of our knowledge, this is the most comprehensive survey of peptides in PR cheese. Most of them were from β-caseins, which represent the best substrate for LAB cell-envelope proteases. The abundance of peptides from β-casein 38-88 region continuously increased during ripening. Remarkably, this region contains precursors for the anti-hypertensive lactotripeptides VPP and IPP, as well as for β-casomorphins. We found that the ripening time strongly affects bioactive peptide profiles and that the occurrence of Lcb. zeae species is positively linked to the incidence of eight anti-hypertensive peptides. This result highlighted how the presence of specific LAB species is likely a pivotal factor in determining PR functional properties.

Assessment of residual plant DNA in bulk milk for Grana Padano PDO production by a metabarcoding approach

The aim of this study was to evaluate the ability of DNA metabarcoding, by rbcl as barcode marker, to identify and classify the small traces of plant DNA isolated from raw milk used to produce Grana Padano (GP) cheese. GP is one of the most popular Italian PDO (Protected Designation of Origin) produced in Italy in accordance with the GP PDO specification rules that define which forage can be used for feeding cows. A total of 42 GP bulk tank milk samples were collected from 14 dairies located in the Grana Padano production area. For the taxonomic classification, a local database with the rbcL sequences available in NCBI on September 2020/March 2021 for the Italian flora was generated. A total of 8,399,591 reads were produced with an average of 204,868 per sample (range 37,002-408,724) resulting in 16, 31 and 28 dominant OTUs at family, genus and species level, respectively. The taxonomic analysis of plant species in milk samples identified 7 families, 14 genera and 14 species, the statistical analysis conducted using alpha and beta diversity approaches, did not highlight differences among the investigated samples. However, the milk samples are featured by a high plant variability and the lack of differences at multiple taxonomic levels could be due to the standardisation of the feed rationing, as requested by the GP rules. The results suggest that DNA metabarcoding is a valuable resource to explore plant DNA traces in a complex matrix such as milk.

Low-Level Clostridial Spores' Milk to Limit the Onset of Late Blowing Defect in Lysozyme-Free, Grana-Type Cheese

The growth of clostridial spores during ripening leads to late blowing (LB), which is the main cause of spoilage in Grana Padano Protected Designation of Origin (PDO) cheese and other hard, long-ripened cheeses such as Provolone, Comté, and similar cheeses. This study aimed to verify the cause-effect relationship between the level of clostridial butyric spores (BCS) in milk and the onset of the LB defect. To this end, experimental Grana-type cheeses were produced without lysozyme, using bulk milk with different average BCS content. The vat milk from the so-called "virtuous" farms (L1) contained average levels of BCS of 1.93 ± 0.61 log most probable number (MPN) L^-1, while the vat milk from farms with the highest load of spores (L2), were in the order of 2.99 ± 0.69 log MPN L^-1. Cheeses after seven months of ripening evidenced a strong connection between BCS level in vat milk and the occurrence of LB defect. In L2 cheeses, which showed an average BCS content of 3.53 ± 1.44 log MPN g^-1 (range 1.36-5.04 log MPN g^-1), significantly higher than that found in L1 cheeses (p < 0.01), the defect of LB was always present, with Clostridium tyrobutyricum as the only clostridial species identified by species-specific PCR from MPN-positive samples. The L1 cheeses produced in the cold season (C-L1) were free of defects whereas those produced in the warm season (W-L1) showed textural defects, such as slits and cracks, rather than irregular eyes. A further analysis of the data, considering the subset of the cheesemaking trials (W-L1 and W-L2), carried out in the warm season, confirmed the presence of a climate effect that, often in addition to the BCS load in the respective bulk milks (L1 vs. L2), may contribute to explain the significant differences in the chemical composition and some technological parameters between the two series of cheeses. Metagenomic analysis showed that it is not the overall structure of the microbial community that differentiates L1 from L2 cheeses but rather the relative distribution of the species between them. The results of our trials on experimental cheeses suggest that a low-level BCS in vat milk (<200 L^-1) could prevent, or limit, the onset of LB in Grana-type and similar cheeses produced without lysozyme.

The rotation of primary starter culture mixtures results in batch-to-batch variations during Gouda cheese production

Industrial production of Gouda cheeses mostly relies on a rotated use of different mixed-strain lactic acid bacteria starter cultures to avoid phage infections. However, it is unknown how the application of these different starter culture mixtures affect the organoleptic properties of the final cheeses. Therefore, the present study assessed the impact of three different starter culture mixtures on the batch-to-batch variations among Gouda cheeses from 23 different batch productions in the same dairy company. Both the cores and rinds of all these cheeses were investigated after 36, 45, 75, and 100 weeks of ripening by metagenetics based on high-throughput full-length 16S rRNA gene sequencing accompanied with an amplicon sequence variant (ASV) approach as well as metabolite target analysis of non-volatile and volatile organic compounds. Up to 75 weeks of ripening, the acidifying Lactococcus cremoris and Lactococcus lactis were the most abundant bacterial species in the cheese cores. The relative abundance of Leuconostoc pseudomesenteroides was significantly different for each starter culture mixture. This impacted the concentrations of some key metabolites, such as acetoin produced from citrate, and the relative abundance of non-starter lactic acid bacteria (NSLAB). Cheeses with the least Leuc. pseudomesenteroides contained more NSLAB, such as Lacticaseibacillus paracasei that was taken over by Tetragenococcus halophilus and Loigolactobacillus rennini upon ripening time. Taken together, the results indicated a minor role of leuconostocs in aroma formation but a major impact on the growth of NSLAB. The relative abundance of T. halophilus (high) and Loil. rennini (low) increased with ripening time from rind to core. Two main ASV clusters of T. halophilus could be distinguished, which were differently correlated with some metabolites, both beneficial (regarding aroma formation) and undesirable ones (biogenic amines). A well-chosen T. halophilus strain could be a candidate adjunct culture for Gouda cheese production.

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.

Artisanal Household Milk Pasteurization Is Not a Determining Factor in Structuring the Microbial Communities of Labneh Ambaris: A Pilot Study

Labneh Ambaris is a traditional Lebanese dairy product traditionally made using raw goat's milk in earthenware jars, but recently the use of artisanally pasteurized milk was introduced for safety reasons. In this study, 12 samples of labneh Ambaris were studied, six made using raw goat's milk and six others using artisanally pasteurized goat's milk. These samples were collected during fermentation and their microbial compositions were analyzed. The 16S V3-V4 and the ITS2 regions of the rDNA were sequenced by DNA metabarcoding analyses for the identification and comparison of bacterial and fungal communities, respectively. The samples had high microbial diversity but differences in samples microbiota were unrelated to whether or not milk was pasteurized. The samples were consequently clustered on the basis of their dominant bacterial or fungal species, regardless of the milk used. Concerning bacterial communities, samples were clustered into 3 groups, one with a higher abundance of Lactobacillus helveticus, another with Lactobacillus kefiranofaciens as the dominant bacterial species, and the third with Lentilactobacillus sp. as the most abundant species. Species belonging to the Enterobacteriaceae family were detected in higher abundance in all raw milk samples than in artisanally pasteurized milk samples. As for fungal communities, the samples were clustered into two groups, one dominated by Geotrichum candidum and the other by Pichia kudriavzevii.

Authenticity and Typicity of Traditional Cheeses: A Review on Geographical Origin Authentication Methods

Food fraud, corresponding to any intentional action to deceive purchasers and gain an undue economical advantage, is estimated to result in a 10 to 65 billion US dollars/year economical cost worldwide. Dairy products, such as cheese, in particular cheeses with protected land- and tradition-related labels, have been listed as among the most impacted as consumers are ready to pay a premium price for traditional and typical products. In this context, efficient food authentication methods are needed to counteract current and emerging frauds. This review reports the available authentication methods, either chemical, physical, or DNA-based methods, currently used for origin authentication, highlighting their principle, reported application to cheese geographical origin authentication, performance, and respective advantages and limits. Isotope and elemental fingerprinting showed consistent accuracy in origin authentication. Other chemical and physical methods, such as near-infrared spectroscopy and nuclear magnetic resonance, require more studies and larger sampling to assess their discriminative power. Emerging DNA-based methods, such as metabarcoding, showed good potential for origin authentication. However, metagenomics, providing a more in-depth view of the cheese microbiota (up to the strain level), but also the combination of methods relying on different targets, can be of interest for this field.

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.

Effect of using mycotoxin-detoxifying agents in dairy cattle feed on natural whey starter biodiversity

Natural whey cultures (NWC) are undefined multiple-strain bacterial starter communities that can be affected by even small changes along the entire dairy chain. We applied a multidisciplinary approach to investigate how the addition of 2 mycotoxin-detoxifying agents [sodium smectite and lignocellulose-based material (B1); leonardite and betaine (B2)] to cow diets modified the microbiota of the NWC in manufacture of a Grana-like cheese. Microbiological and flow cytometry analyses showed that the content and viability of lactic acid bacteria (LAB) and the total whey microbiota were not affected by the detoxifying agents, and Streptococcus thermophilus, Lactobacillus helveticus, and Limosilactobacillus fermentum were the dominant taxa. Random amplified polymorphic DNA-PCR fingerprinting and metagenomic analysis highlighted differences in the bacterial community of the NWC and in the relative abundance of Bacteroidetes that increased when B1 and B2 were included in the diet. Two of 6 St. thermophilus biotypes were detected only in control samples; conversely, none of the Lb. helveticus biotypes found in control samples were isolated from B1 and B2. In vitro tests showed that the 2 binders did not significantly affect the development of St. thermophilus, but they stimulated the growth of Lb. helveticus strains recovered only from B1 and B2 NWC. The addition of binders in cow feed can affect the LAB biotypes present in NWC.

Dairy streptococcal cell wall and exopolysaccharide genome diversity

The large-scale and high-intensity application of Streptococcus thermophilus species in milk fermentation processes is associated with a persistent threat of (bacterio)phage infection. Phage infection of starter cultures may cause inconsistent, slow or even failed fermentations with consequent diminished product quality and/or output. The phage life cycle commences with the recognition of, and binding to, a specific host-encoded and surface-exposed receptor, which in the case of S. thermophilus can be the rhamnose-glucose polysaccharide (RGP; specified by the rgp gene cluster) or exopolysaccharide (EPS; specified by the eps gene cluster). The genomic diversity of 23 S . thermophilus strains isolated from unpasteurized dairy products was evaluated, including a detailed analysis of the rgp and eps loci. In the present study, five novel eps genotypes were identified while variations of currently recognized rgp gene cluster types were also observed. Furthermore, the diversity of rgp genotypes amongst retrieved isolates positively correlated with phage diversity based on phageome analysis of eight representative dairy products. Our findings therefore substantially expand our knowledge on S. thermophilus’ strain and phage diversity in (artisanal) dairy products and highlight the merit of phageome analysis of artisanal and traditional fermented foods as a sensitive marker of dominant microbiota involved in the fermentation.

High-throughput qPCR and 16S rRNA gene amplicon sequencing as complementary methods for the investigation of the cheese microbiota

Background

Next-generation sequencing (NGS) methods and especially 16S rRNA gene amplicon sequencing have become indispensable tools in microbial ecology. While they have opened up new possibilities for studying microbial communities, they also have one drawback, namely providing only relative abundances and thus compositional data. Quantitative PCR (qPCR) has been used for years for the quantification of bacteria. However, this method requires the development of specific primers and has a low throughput. The constraint of low throughput has recently been overcome by the development of high-throughput qPCR (HT-qPCR), which allows for the simultaneous detection of the most prevalent bacteria in moderately complex systems, such as cheese and other fermented dairy foods. In the present study, the performance of the two approaches, NGS and HT-qPCR, was compared by analyzing the same DNA samples from 21 Raclette du Valais protected designation of origin (PDO) cheeses. Based on the results obtained, the differences, accuracy, and usefulness of the two approaches were studied in detail.

Results

The results obtained using NGS (non-targeted) and HT-qPCR (targeted) show considerable agreement in determining the microbial composition of the cheese DNA samples studied, albeit the fundamentally different nature of these two approaches. A few inconsistencies in species detection were observed, particularly for less abundant ones. The detailed comparison of the results for 15 bacterial species/groups measured by both methods revealed a considerable bias for certain bacterial species in the measurements of the amplicon sequencing approach. We identified as probable origin to this PCR bias due to primer mismatches, variations in the number of copies for the 16S rRNA gene, and bias introduced in the bioinformatics analysis.

Conclusion

As the normalized microbial composition results of NGS and HT-qPCR agreed for most of the 21 cheese samples analyzed, both methods can be considered as complementary and reliable for studying the microbial composition of cheese. Their combined application proved to be very helpful in identifying potential biases and overcoming methodological limitations in the quantitative analysis of the cheese microbiota.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02451-y.

Assessment of the microbiological origin of blowing defects in Grana Padano Protected Designation of Origin cheese

Recognized worldwide for its history, flavor, and high nutritional quality, Grana Padano (GP) is one of the most traditional Italian raw-milk, hard-cooked, long-ripened cheese. Throughout GP manufacturing, some well-known and undesired bacterial species, such as clostridia, can proliferate and lead to spoilage defects that mischaracterize the final product; however, little is known about the development of late-blowing defects in hard cheese samples without clostridia. Therefore, in this study we aimed to use metataxonomic analysis to identify bacterial taxa associated with the development of late-blowing defect in GP samples. Furthermore, the presence of several heterofermentative lactobacilli species in defective zones were verified by primer-specific PCR assay. Considering α- and β-diversity analyses, no statistically significant differences were detected between cheese samples with and without blowing defect. Following taxonomic assignment, Lactobacillus and Streptococcus were the dominant genera, whereas clostridia-related taxa were not detected in any of the 20 analyzed samples. Using EdgeR, the genera Propionibacterium and Acinetobacter were found to be prevalently more abundant in samples categorized as having "big regular holes." In samples with "small regular holes," multiplex PCR amplification revealed differences in terms of Lactobacillus population composition, mainly obligate homofermentative lactobacilli, between defective and non-defective zones of the same cheese wheel. This study demonstrated that GP samples with blowing defects not caused by clostridial development share similar biodiversity indices with GP collected from control zones, but an imbalance of obligate homofermentative lactobacilli was noticed between samples, which requires further analysis to better comprehend the exact mechanism involved in this process.

Advantages and disadvantages of non-starter lactic acid bacteria from traditional fermented foods: Potential use as starters or probiotics

Traditional fermented foods are a significant source of starter and/or non-starter lactic acid bacteria (nsLAB). Moreover, these microorganisms are also known for their role as probiotics. The potential of nsLAB is huge; however, there are still challenges to be overcome with respect to characterization and application. In the present review, the most important steps that autochthonous lactic acid bacteria isolated from fermented foods need to overcome, to qualify as novel starter cultures, or as probiotics, in food technology and biotechnology, are considered. These different characterization steps include precise identification, detection of health-promoting properties, and safety evaluation. Each of these features is strain specific and needs to be accurately determined. This review highlights the advantages and disadvantages of nsLAB, isolated from traditional fermented foods, discussing safety aspects and sensory impact.

The Microbiota of Grana Padano Cheese. A Review

Grana Padano (GP) is the most appreciated and marketed cheese with Protected Designation of Origin in the world. The use of raw milk, the addition of undefined cultures (defined as 'sieroinnesto naturale'), the peculiar manufacturing proces, and the long ripening make the cheese microbiota play a decisive role in defining the quality and the organoleptic properties of the product. The knowledge on the microbial diversity associated with GP has been the subject, in recent years, of several studies aimed at understanding its composition and characteristics in order, on the one hand, to improve its technological performances and, on the other hand, to indirectly enhance the nutritional quality of the product. This review aims to briefly illustrate the main available knowledge on the composition and properties of the GP microbiota, inferred from dozens of studies carried out by both classical microbiology techniques and metagenomic analysis. The paper will essentially, but not exclusively, be focused on the lactic acid bacteria (LAB) derived from starter (SLAB) and the non-starter bacteria, both lactic (NSLAB) and non-lactic, of milk origin.

Biodiversity and succession of lactic microbiota involved in Brazilian buffalo mozzarella cheese production

The biodiversity and succession of lactic acid bacteria (LAB) involved in the production and storage of Brazilian buffalo mozzarella cheese were evaluated. The isolates were characterized by Gram staining and catalase test, by the ability to grow at different conditions: temperatures, pH, concentrations of NaCl, and production of CO₂ from glucose. The biodiversity and succession of 152 LAB isolated during cheese production were evaluated by 16S rRNA gene sequencing, Random Amplified Polymorphic DNA (RAPD-PCR), and Restriction Fragment Length Polymorphism (RFLP-PCR) techniques. Most of the strains grow well at 30 °C and are tolerant to 6.5% of NaCl, and in general, the best pH for growing was 9.6. Leuconostoc mesenteroides, Lacticaseibacillus casei, Limosilactobacillus fermentum, and Enterococcus sp. were prevalent and present in almost all steps of production. The LAB strains are typically found in the traditional Italian cheese, except the Leuconostoc citreum species. Sixty clusters were obtained by RAPD-PCR with 85% of similarity (114 isolates) while most of the LAB was clustered with 100% of similarity by the RFLP-PCR technique. The applied techniques enabled a valuable elucidation of the LAB biodiversity and succession, contributing to a better understanding of the specific microbial cultures with a technological aptitude of this cheese.

Bacterial Community of Grana Padano PDO Cheese and Generical Hard Cheeses: DNA Metabarcoding and DNA Metafingerprinting Analysis to Assess Similarities and Differences

The microbiota of Protected Designation of Origin (PDO) cheeses plays an essential role in defining their quality and typicity and could be applied to protect these products from counterfeiting. To study the possible role of cheese microbiota in distinguishing Grana Padano (GP) cheese from generical hard cheeses (HC), the microbial structure of 119 GP cheese samples was studied by DNA metabarcoding and DNA metafingerprinting and compared with 49 samples of generical hard cheeses taken from retail. DNA metabarcoding highlighted the presence, as dominant taxa, of Lacticaseibacillus rhamnosus, Lactobacillus helveticus, Streptococcus thermophilus, Limosilactobacillus fermentum, Lactobacillus delbrueckii, Lactobacillus spp., and Lactococcus spp. in both GP cheese and HC. Differential multivariate statistical analysis of metataxonomic and metafingerprinting data highlighted significant differences in the Shannon index, bacterial composition, and species abundance within both dominant and subdominant taxa between the two cheese groups. A supervised Neural Network (NN) classification tool, trained by metagenotypic data, was implemented, allowing to correctly classify GP cheese and HC samples. Further implementation and validation to increase the robustness and improve the predictive capacity of the NN classifier will be needed. Nonetheless, the proposed tool opens interesting perspectives in helping protection and valorization of GP and other PDO cheeses.

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.

Microbial Ecology of Greek Wheat Sourdoughs, Identified by a Culture-Dependent and a Culture-Independent Approach

The aim of the present study was to assess the microecosystem of 13 homemade spontaneously fermented wheat sourdoughs from different regions of Greece, through the combined use of culture-dependent (classical approach; clustering by Random Amplified Polymorphic DNA-Polymerase Chain Reaction (RAPD-PCR) and identification by PCR species-specific for Lactiplantibacillus plantarum, and sequencing of the 16S-rRNA and 26S-rRNA gene, for Lactic Acid Bacteria (LAB) and yeasts, respectively) and independent approaches [DNA- and RNA-based PCR-Denaturing Gradient Gel Electrophoresis (DGGE)]. The pH and Total Titratable Acidity (TTA) values ranged from 3.64–5.05 and from 0.50–1.59% lactic acid, respectively. Yeast and lactic acid bacteria populations ranged within 4.60–6.32 and 6.28–9.20 log CFU/g, respectively. The yeast: LAB ratio varied from 1:23–1:10,000. A total of 207 bacterial and 195 yeast isolates were obtained and a culture-dependent assessment of their taxonomic affiliation revealed dominance of Lb. plantarum in three sourdoughs, Levilactobacillus brevis in four sourdoughs and co-dominance of these species in two sourdoughs. In addition, Companilactobacillusparalimentarius dominated in two sourdoughs and Fructilactobacillussanfranciscensis and Latilactobacillus sakei in one sourdough each. Lactococcus lactis, Lb. curvatus, Leuconostoc citreum, Ln. mesenteroides and Lb. zymae were also recovered from some samples. Regarding the yeast microbiota, it was dominated by Saccharomyces cerevisiae in 11 sourdoughs and Pichia membranifaciens and P. fermentans in one sourdough each. Wickerhamomyces anomalus and Kazachstania humilis were also recovered from one sample. RNA-based PCR-DGGE provided with nearly identical results with DNA-based one; in only one sample the latter provided an additional band. In general, the limitations of this approach, namely co-migration of amplicons from different species to the same electrophoretic position and multiband profile of specific isolates, greatly reduced resolution capacity, which resulted in only partial verification of the microbial ecology detected by culture-dependent approach in the majority of sourdough samples. Our knowledge regarding the microecosystem of spontaneously fermented Greek wheat-based sourdoughs was expanded, through the study of sourdoughs originating from regions of Greece that were not previously assessed.