16S-rRNA-Based Metagenomic Profiling of the Bacterial Communities in Traditional Bulgarian Sourdoughs

Abstract

Sourdoughs (SDs) are spontaneously formed microbial ecosystems composed of various species of lactic acid bacteria (LAB) and acid-tolerant yeasts in food matrices of cereal flours mixed with water. To date, more than 90 LAB species have been isolated, significantly impacting the organoleptic characteristics, shelf life, and health properties of bakery products. To learn more about the unique bacterial communities involved in creating regional Bulgarian sourdoughs, we examined the metacommunities of five sourdoughs produced by spontaneous fermentation and maintained by backslopping in bakeries from three geographic locations. The 16S rRNA gene amplicon sequencing showed that the former genus Lactobacillus was predominant in the studied sourdoughs (51.0–78.9%). Weissella (0.9–42.8%), Herbaspirillum (1.6–3.8%), Serratia (0.1–11.7%), Pediococcus (0.2–7.5%), Bacteroides (0.1–1.3%), and Sphingomonas (0.1–0.5%) were also found in all 5 samples. Genera Leuconostoc, Enterococcus, Bacillus, and Asaia were sample-specific. It is interesting to note that the genus Weissella was more abundant in wholegrain samples. The greatest diversity at the species level was found in the former genus Lactobacillus, presented in the sourdough samples with 13 species. The UPGMA cluster analysis clearly demonstrated similarity in species’ relative abundance between samples from the same location. In addition, we can conclude that the presence of two main clusters—one including samples from mountainous places (the cities of Smolyan and Bansko) and the other including samples from the city of Ruse (the banks of the Danube River)—may indicate the impact of climate and geographic location (e.g., terrain, elevation, land use, and nearby water bodies and their streams) on the abundance of microbiome taxa. As the bacterial population is crucial for bread standardization, we expect the local bakery sector to be interested in the relationship between process variables and their effect on bacterial dynamics described in this research study.