Metabolic framework of spontaneous and synthetic sourdough metacommunities to reveal microbial players responsible for resilience and performance
Microbiome volume 10, Article number: 148 (2022)
In nature, microbial communities undergo changes in composition that threaten their resiliency. Here, we interrogated sourdough, a natural cereal-fermenting metacommunity, as a dynamic ecosystem in which players are subjected to continuous environmental and spatiotemporal stimuli.
The inspection of spontaneous sourdough metagenomes and transcriptomes revealed dominant, subdominant and satellite players that are engaged in different functional pathways. The highest microbial richness was associated with the highest number of gene copies per pathway. Based on meta-omics data collected from 8 spontaneous sourdoughs and their identified microbiota, we de novo reconstructed a synthetic microbial community SDG. We also reconstructed SMC-SD43 from scratch using the microbial composition of its spontaneous sourdough equivalent for comparison. The KEGG number of dominant players in the SDG was not affected by depletion of a single player, whereas the subdominant and satellite species fluctuated, revealing unique contributions. Compared to SMC-SD43, SDG exhibited broader transcriptome redundancy. The invariant volatilome profile of SDG after in situ long-term back slopping revealed its stability. In contrast, SMC-SD43 lost many taxon members. Dominant, subdominant and satellite players together ensured gene and transcript redundancy.
Our study demonstrates how, by starting from spontaneous sourdoughs and reconstructing these communities synthetically, it was possible to unravel the metabolic contributions of individual players. For resilience and good performance, the sourdough metacommunity must include dominant, subdominant and satellite players, which together ensure gene and transcript redundancy. Overall, our study changes the paradigm and introduces theoretical foundations for directing food fermentations.
Microbiomes are vital components of natural ecosystems whose functions are typically performed not by single species but by metacommunities . The structure, function, resilience, and stability of these metacommunities result from a dynamic interplay of natural selection, historical contingency, and chance events in a manner that remains poorly understood. The ability to predict and reconstruct large multispecies communities requires an understanding of how such microbiomes shape, behave and coexist in natural niches to ultimately design de novo functional microbiomes [2,3,4]. These concepts translate easily to food microbiomes, and specifically to the sourdough microbiome, or “fermentome” (fermenting metacommunity), as we prefer to define this highly variable ecosystem. Here, we experimentally interrogated sourdough as a model dynamic microbial ecosystem. As one of the oldest examples of a natural starter, sourdough is increasingly used for making baked goods, and almost 30 years of research has been carried out to understand and define its performance under abiotic and biotic stresses . The main secret to sourdough performance lies in its microbial diversity. Up to 59 bacterial genera were culturable in sourdoughs, with Lactobacillus, under its former taxonomic nomenclature, as the most abundant genus, comprising 82 species, mostly nomadic and heterofermentative . Yeasts are also irreplaceable players in the sourdough metacommunity, with Saccharomyces cerevisiae being the most frequently identified . The house microbiota, that of the flours, with resilient and resistant bacteria, also behave as cereal plant endophytes, and those of the other eventual ingredients in the dough serve as the main sources of microbial contamination and spontaneous sourdough fermentation . Frequent back slopping, as a self-renewing community-scale inoculum, and an acidic environment generate a complex microbiome of interacting players (dominants, subdominants and satellites), and the structure and metabolic network of this microbiome are constantly influenced by consistent biotic pressure and spatiotemporal changes. As a result, in most cases, sourdough is prone to instability such that the metabolic networks are constantly threatened. Metacommunity interplay is, therefore, critical in maintaining stability as a precondition for ensuring sourdough functionality. Consequently, sourdough composition and function cannot be estimated or predicted simply by summing its single parts or monitoring the most abundant species . Basic approaches are needed to anticipate the response of the metacommunity, to steer coexisting species to more desirable and resilient states and to design microbiomes de novo. The synthetic microbial communities (SMCs) approach constitutes the only route to confirm or invalidate the hypotheses of our study [1, 7]. SMCs allow a more detailed study of the resilience and performance of the metacommunity in simple and well-controlled habitats, where selective and nonselective forces shape the microbiome in mechanistic and quantitative ways. We assume that subdominant players cooperate synergistically with dominant species by expressing accessory genes that are critical for key metabolic pathways and that, although not encoding key functions, satellite members may be metabolically active and equipped with unique genes and enzymes. More generally, gene redundancy should prevent the loss of specific metabolic pathways in the case of microbial perturbance.
Using a multiomics workflow (Fig. 1), we first deciphered the taxon composition and the functional redundancy of spontaneous sourdoughs at both the metagenomics and metatranscriptomics levels; then, based on these data, we reconstructed the main metabolic pathways that are able to exert key functions within themselves. Subsequently, we used this backbone knowledge to design specifically depleted synthetic microbial communities (SMCs) (where one species at a time was subtracted from the consortium) that guided us in monitoring transcript profile changes during microbial species interactions. As milestones in the workflow, we reconstructed two additional de novo SMCs. The first mimicked a spontaneous sourdough, whereas the second resulted from the funnel selection of species based on omics results. The comparison of their two volatilomes at different back slopping times allowed us to validate the good performance of our selected communities, which could be measured in terms of resiliency, stability, and robustness.
The meta-omics inspection proposed here was useful in laying a theoretical foundation for steering the metacommunity by reconstructing a sourdough fermentome that is as resilient as possible to in situ perturbations.
Traditional and spontaneous sourdoughs
Traditional and spontaneous sourdoughs were obtained from the International Sourdough Library (Puratos, St. Vith, Belgium, https://www.puratos.com/commitments/next-generation/product-heritage/sourdough-library), the only one recognized worldwide. In particular, the 8 sourdoughs used were from Italy (SD1 and SD69), the USA (SD43, SD102 and SD104), Spain (SD88 and SD93) and France (SD44). SD69, SD93 and SD102 were produced with soft wheat flour, with final dough yields (DY) of 143, 160 and 200, respectively. SD43 (DY = 150) and SD104 (DY = 200) were produced with strong soft wheat flour; SD1 (DY = 209), with durum wheat flour; SD44 (DY = 200), with type 80 (specific residual of ash content) soft wheat flour with added honey; and SD88 (DY = 200), with type 80 soft wheat flour with added beer. All sourdoughs were considered mature (constant acidifying capacity) after two back slopping steps, in which the first step included 40% (w/w) inoculum, while the second step used a percentage unique to each sourdough, ranging from 8 to 50%. The back slopping temperature was also unique to each sourdough (from 8 °C up to 28 °C), and the fermentation time ranged from 2 to 24 h. These sourdoughs were the most representative of the International Sourdough Library because of their broad use for making leavened baked goods, and the different geographic provenances reflect countries with a long history and tradition of using sourdough. Supplementary Table S1, Additional File 1 summarizes their origin, ingredients, and technology parameters.
The eight mature spontaneous sourdoughs were subjected to analyses in triplicate using M17 (30 °C, 48 h), mostly for coccus-shaped lactic acid bacteria, and SDB, mMRS and MRS5 (30 °C, 48 h), mostly for dominant lactobacilli and Weissella. Slanetz & Bartley agar (spread plate method, 37 °C, 24–48 h) was used for enterococci, Baird-Parker agar (spread plate method, 30 °C, 48 h) for staphylococci and micrococci, VRBGA for total coliforms, and WA and SDA (30 °C, 48 h) for yeasts. For all media, colonies of bacteria and yeasts with different morphologies were picked and isolated from the penultimate dilution, with the exception of the subdominant culturable lactobacilli, which was isolated from 140 mm diameter plates and considering several decimal dilutions. Isolated colonies were cultivated in broth media and then restreaked onto agar medium (MRS for all lactic acid bacteria and SDA for yeasts) until they were purified. Bacterial and yeast isolates were identified by partial sequencing of the 16S rRNA and 26S rRNA genes, respectively [8, 9]. The identified strains comprising the sourdough biobank were further used upon selection to reconstruct the de novo synthetic microbial communities.
Shotgun metagenomics and metatranscriptomics sequencing
The DNeasy PowerFood Microbial Kit (Qiagen, Hilden, Germany) was used for DNA extraction . To improve the DNA yield, 100 μl lysozyme (10 mg/ml), 10 μl mutanolysin (20 U/μl  and 2 μl Zymolyase (5 U/μl)  were added. All library preparations, NGS and quality control steps were performed through RTL-Genomics (Lubbock, Texas). Sequencing was performed on the Illumina MiSeq platform (Illumina Inc., San Diego, CA) applying a depth of coverage of 70X and a fragment size of 300X2 bp paired end (PE) reads. The sequence quality was assessed by inspecting raw read data with FastQC software (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Starting from the fastq raw sequencing files, within the SqueezeMeta automated bioinformatics pipeline , we used Trimmomatic  for adapter removal, trimming and filtering the reads by quality. The details of the bioinformatics workflow are provided below in this section.
Total RNA was extracted using the RNeasy PowerMicrobiome Kit (Qiagen, Hilden, Germany). For each RNA sample, a directional library was prepared using the TruSeq Stranded Total RNA Sample Prep Kit with Ribo-Zero Plus rRNA Depletion Technology (Illumina, San Diego, CA, USA). In detail, RNA-Seq libraries were prepared from 100 ng of total RNA following the manufacturer’s instructions. The cDNA libraries obtained were validated using the High Sensitivity DNA assay (Agilent Technologies, Santa Clara, CA, USA) on a Bioanalyzer 2100 instrument (Agilent Technologies, Santa Clara, CA, USA) and quantified by fluorimetry with the Quant-iTTM PicoGreen® dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA) on a NanoDrop™ 3300 fluorospectrometer (Thermo Scientific, Waltham, MA, USA). Finally, RNA-Seq libraries were pooled in equimolar ratios and sequenced at a concentration of 1.3 pM on the Illumina NextSeq 500 platform (Illumina, San Diego, CA, USA), generating ca. 20–25 M paired 75 bp reads for each sample.
SqueezeMeta bioinformatics parameters and nested software
The sequencing data of the 8 sourdough metagenomes/metatranscriptomes were analyzed in silico by using the SqueezeMeta pipeline Version 1.0, July 2019 (https://github.com/jtamames/SqueezeMeta) plus other related ad hoc custom utilities developed to handle the assembly of metagenomes and to obtain genome bins. A step-by-step bioinformatics pipeline comprising software for assembly, annotation and bin statistics was used for sourdough gene and transcript annotation and quantification.
The inspection of annotation data was fundamental for the retrieval of unique/shared functions up to the species taxonomic level. As a dedicated option, we chose to run the metagenome and metatranscriptome assembly using the “coassembly” mode. This procedure allowed us to pool reads from all samples, thus performing a single assembly. Therefore, the derived per sample gene abundances come from mapping the singular read sample set against the coassembly. The assembly was performed using Megahit (v1.1.2) , which makes use of succinct de Bruijn graphs  to take full advantage of multiple k-mer sizes, optimizing both sensitivity and accuracy. Contig statistics were determined using prinseq , which performs rapid quality control and data preprocessing of genomic and metagenomic datasets. The locations of the ribosomal RNA genes in genomes were predicted using Barrnap (BAsic Rapid Ribosomal RNA Predictor) , which supports bacteria (5S, 23S, and 16S), archaea (5S, 5.8S, 23S, and 16S) and eukaryotes (5S, 5.8S, 28S, and 18S). The 16S rRNA sequences were taxonomically classified using the RDP classifier . tRNA/tmRNA sequences were predicted using Aragorn . ORFs were predicted using Prodigal . Similarity searches for GenBank , eggNOG , and KEGG  were performed using Diamond . HMM homology searches were performed by HMMER3  for the Pfam database . Read mapping against contigs was performed using Bowtie2 , and binning was performed using MaxBin2  and Metabat2 . The combination of binning results was performed using DAS Tool . Bin statistics were computed using CheckM . Pathway prediction for the KEGG  and MetaCyc databases was performed using MinPath .
Construction of metagenomics and metatranscriptomics project databases
SqueezeMeta outputs were used to build customized databases viewable in the R environment. Specifically, a set of integrated functions implemented within the SQMtools R package Version 0.3.3 (https://github.com/jtamames/SqueezeMeta/wiki/Using-R-to-analyze-your-SQM-results) allowed us to query specific KEGG output results. The whole metagenomics and metatranscriptomics projects were used to create the two relative databases and then inspected for gene and transcript content and abundances, respectively.
Synthetic microbial community (SMC) study design
A meta-omics approach was used to design SMCs from the 8 spontaneous sourdoughs. Quantitative data (cell densities) from culturomics were confirmed by the genomic potential and expression levels derived from the metagenomics and metatranscriptomics analyses. By inspecting and merging all these omics data, the sourdough members were classified as core dominant, subdominant (core or dispensable) or satellite species. Starting from this classification, we aimed to de novo reconstruct a potentially stable and resilient SMC (Sourdough Global, SDG). Here, dominant and subdominant species were those harboring at least 20 key genes belonging to carbohydrate/pyruvate and nitrogen metabolism and, at the same time, were shared by at least 50% (4 out of 8) of traditional sourdoughs. Thirteen species out of 49 dominant and subdominant lactic acid bacteria and yeasts met these criteria. Upon considering this filter, further considerations dealing with the overall metabolic capabilities, and worldwide frequency of isolation, SDG comprised 1 core (Sac. cerevisiae) and 1 dispensable dominant yeast (Pichia kudriavzevii), 2 core dominant (Lactiplantibacillus plantarum and Limosilactobacillus fermentum) and 2 core subdominant (Furfurilactobacillus rossiae and Pediococcus pentosaceus) lactic acid bacteria, and 1 satellite species (Staphylococcus epidermidis) (Table 1). The species selected to reconstruct the SDG originated from the sourdough biobank and more specifically from SD104, which was the traditional spontaneous sourdough with the highest species biodiversity and functional potential. Strains within the same species from SD104 were selected randomly. Another SMC (SMC-SD43), comprising all the species encompassed by sourdough SD43, was reconstructed with the aim of mimicking the potential of a spontaneous sourdough. The cell densities of all species used to reproduce the SMCs corresponded to those found in traditional spontaneous sourdoughs, as estimated by culturomics (Table 1). One experiment with these 2 SMCs involved their daily propagation at 30 °C for 30 days under in vivo sourdough conditions. Plate counts and isolation for RAPD-PCR profiling  were performed every 10 days. Analysis and comparison of biotypes were performed with BioNumerics software (v. 8.0, Applied Maths) using the reference profiles of each bacterial and yeast species belonging to the 2 SMCs. Metabolomics analyses allowed the characterization of volatile organic compounds (VOCs) after 1 and 30 days of propagation. To simulate the sourdough-like environment under sterile conditions, another experiment used WFH (wheat flour hydrolyzed) liquid medium . The 2 SMCs were inoculated with 7 log cfu/mL of core dominant bacteria, 5 log cfu/mL of core subdominant bacteria and satellites, and 6 log cfu/mL of yeasts. The fermentation lasted for 12 h at 30 °C. With the aim of assessing how the elimination of microbial members might affect the transcriptomic profile of SDG, 7 SMCs were reconstructed by depleting the consortium 1 member at a time (SDG1 to SDG7). Acidification and growth kinetics were recorded using a benchtop online pH meter with a food probe (Hanna Instruments, Woonsocket, RI, USA) and a Tecan Infinite M NANO + spectrophotometer (Tecan Ltd., Switzerland). The growth kinetics were modeled according to the Gompertz Eq .. Metatranscriptomes were sequenced at exponential (6 h) and stationary (12 h) phases of growth. The bioinformatics pipeline described in the previous section (SqueezeMeta bioinformatics parameters and nested software) and the bioinformatics approach described above allowed the annotation and quantification of transcripts. An ad hoc reconstructed database allowed us to query the transcriptome project data and to extract per-species contributions.
Metabolic networking of sourdough core, dispensable and satellite microbial members
Based on the importance of sourdough performance, we selected the following metabolic mega-pathways: carbohydrate metabolism, pyruvate/energy production-conversion and nitrogen metabolism. Core enzymatic pathways, including enzymes detected in all sourdough metagenomes, were identified. In addition, dispensable/accessory genes encoding for the enzymatic portfolio among sourdough metagenomes were identified by inspecting the ad hoc build up sourdough SqueezeMeta database. KEGG abundance values (expressed as transcripts per million (TPM) values after normalizing for both sequencing depth and gene length) were used for a two-sample G-test (w/Yates’ + Fisher test) with Benjamini–Hochberg multiple test correction in STAMP (statistical analysis of taxonomic and functional profiles) software . Metatranscriptomics was used together with metagenomics to find the active pathways within sourdoughs.
The VOC profile was analyzed by gas chromatography–mass spectrometry (GC–MS). A PAL COMBI-xt autosampler (CTC combiPAL, CTC Analysis AG, Zwingen, Switzerland) was used to standardize the headspace solid-phase microextraction (HS-SPME) procedure according to Liu et al. (2020) . A Clarus 680 (Perkin Elmer, Beaconsfield, UK) gas chromatograph equipped with a Rtx-WAX column (30 m × 0.25 mm i.d., 0.25 μm film thickness) (Restek Superchrom, Milano, Italy) was used to thermally desorb and separate the headspace VOC . Each chromatogram was analyzed for peak identification by comparing (i) the retention time (RT) of the detected compound with those of the provided pure standard for HPLC (Sigma–Aldrich, St. Louis, MO, USA) and (ii) experimental mass spectra with those of the National Institute of Standards and Technology database (NIST/EPA/NIH Mass Spectral Library with Search Program, data version NIST 05, software version 2.0d). A peak area threshold larger than 1 M and a match percentage higher than 85% were used as the criteria for identification. VOCs were quantified by using 2-methyl-4-pentanol (final concentration of 33 mg/l) as the internal standard.
Metagenomics and metatranscriptomics data belonging to eight spontaneous sourdoughs and sixteen metatranscriptomics datasets from seven depleted SMCs plus one complete reconstructed SMC (at the exponential and stationary phases of growth) were analyzed. Two de novo reconstructed SMCs (SDG and SMC-SD43) were propagated for 30 days, and the single contributions of dominant, subdominant and satellite species were used to populate the subset of the investigated metabolic pathway functions. In turn, these data have been compared with those of the reconstructed sourdough (SDG1-SDG7). Statistically significant gene and transcript comparisons were obtained using an ANOVA Tukey–Kramer post hoc comparison test (corrected for multiple tests with the Benjamini–Hochberg procedure) in STAMP (statistical analysis of taxonomic and functional profiles) v2.1.3 software . The input file matrices from the metagenomics and metatranscriptomics data and the relative metadata were customized to fit the STAMP software requirements. Culturomics analyses were carried out on 2 biological replicates for the same length of time (up to 30 days) as for VOC analyses.
Sourdough richness in culturable bacteria and yeasts
Our isolation, purification and characterization procedures resulted in the creation of a sourdough biobank comprising a total of 1,661 isolates (1,488 bacteria and 173 yeasts). Spontaneous sourdoughs harbored presumptive lactic acid bacteria at ca. 2.0 to ca. 9.9 log cfu/g. Presumptive staphylococci and micrococci ranged from ca. 2.0 to ca. 7.0 log cfu/g. Presumptive coliforms were identified only in SD1 (ca. 2.5 log cfu/g) and SD69 (ca. 4.5 log cfu/g). All sourdoughs harbored ca. 6.5 log cfu/g yeasts. Partial sequencing of the 16S rRNA and 26S rRNA genes of all 1,661 isolates allowed the taxonomic identification of 36 species. The sourdough richness in culturable bacteria and yeast species varied from 8 (SD93) to 15 (SD44 and SD104), with differences in cell densities and prevalence. The sourdough species prevalence was calculated as the number of isolates identified per species divided by their total number. This percentage varied from 1 to 67%, mostly depending on the abundance level of each species (Table 2). Following this culturomics approach, the sourdough biobank also allowed us to obtain an overview of the ratio among players. Based on the correlation between back slopping steps and LAB cell density until sourdough reaches its maturity, we defined a species to be dominant or subdominant based on whether its cell density was greater or lower than 6 log cfu/g . In addition, considering commensal interaction as a primary impacting factor, the natural selection of species until sourdough is mature is dependent on cell load and other technological parameters whose impacts can be even more important .
Moreover, we consider satellite species to be those nonlactic acid bacteria that are occasionally found in sourdough. Lac. plantarum was dominant or codominant in all sourdoughs. Leuconostoc citreum was also identified in all sourdoughs, with cell densities varying from ca. 3.2 to 9.0 log cfu/g. Except for SD1, enterococci were present at low cell density (ca. 2 log cfu/g) only in some sourdoughs. Except for SD104, staphylococcal species were mostly identified as satellite members. Sac. cerevisiae was the dominant yeast found in all spontaneous sourdoughs, with cell densities ranging from 6.1 to 6.4 log cfu/g. The only exception was SD88, which harbored Pic. kudriavzevii (ca. 6.5 log cfu/g).
Sourdough richness in bacteria and yeasts as estimated by metagenomics
Approximately 6.9 M contigs, for a total of 2.8 Gbp, led to the assembly combining all sourdoughs (Supplementary Table S2, Additional File 2). The longest contig was ca. 245 Kbp. The assembled N50 and N90 values were 813 and 516 base pairs, respectively. Looking at the rank assignation for the whole assembly, 195, 333 and 385 contigs were assigned to family, genus and species ranks, respectively. Among prokaryotes, sourdoughs harbored 202 genera and 275 species (Supplementary Table S3, Additional File 3). Eukaryotic taxa (yeasts) accounted for 30 genera and 31 species (Supplementary Table S4, Additional File 4). The high biodiversity we detected was ultimately reflected in the species relative abundances, ranging from 0.0001 to 80%. For metagenomics data, considering those species that were present in at least one sourdough, we defined dominant or subdominant species based on whether their percentage of relative abundance was greater or lower than 5%, as normalized to the total metagenome content. Similarly, in other metagenomics studies, dominance/subdominance was defined based on species richness and relative abundance level [41, 42]. Regardless of the relative abundance, we defined core or dispensable species according to whether they were shared by all sourdoughs or not (Fig. 2). Within the core microbiome set of species, Lac. plantarum, Fru. sanfranciscensis, Lim. fermentum, Leuc. citreum and Weissella confusa were dominant in 1 or more sourdoughs, which agreed with culturomics data (Table 2). Eighteen subdominant species of lactic acid bacteria were also found. The core sourdough microbiome also harbored other bacterial species. The core yeast microbiome (Saccharomycodes ludwigii and Sac. cerevisiae) was dominant in all sourdoughs except for SD1 and SD69, in which Saccharomyces eubayanus was dominant. Dispensable dominant yeasts included Pichia kudriavzevii, Clavispora lusitaniae and Saccharomyces bayanus, which were shared by 5 and 4 sourdoughs, respectively (Supplementary Figure S1, Additional File 5).
Inspecting gene annotation and deciphering the functional redundancy of sourdough genomes
To decipher the extent of functional redundancy in sourdough genomes, we retrieved the gene abundances in each sample by mapping a single set of reads (from each individual sample) against the total assembly obtained by pooling reads from all samples in a single bin (Supplementary Table S5, Additional File 6). The highest mapping percentage (84.47) was found for SD43. To detect missing genes and correct errors in gene prediction, we used a double-pass step procedure named “extrasensitive ORF prediction”, which consisted of gene prediction as a first step and BLASTX search as a second step.
The total predicted ORF number was nonhomogeneous (Supplementary Table S6, Additional File 7). Considering the Prodigal results, SD1, SD69 and SD93 showed comparable numbers (over 1 M ORFs, with SD1 having a peak of approximately 2 M ORFs). SD43, SD44, SD88, SD102 and SD104 exhibited fewer than 1 M total predicted ORFs. On the other hand, the number of predicted KEGG functions was more homogeneously distributed, with the lowest and highest numbers (ca. 61 K and 551 K) found in SD102 and SD104, respectively. The set of KEGG functions shared among all spontaneous sourdoughs (core metagenome) spanned 1163 genes (Supplementary Table S7, Additional File 8). The other accessory genes were heterogeneously distributed. A set of unique genes was found in each sourdough, ranging from 33 (SD43) to 684 (SD1) unique genes. Adhering to the concept that housekeeping gene functions were always maintained by every species, sourdoughs with a higher number of species showed a much higher number of copies per pathway for these functions (Supplementary Table S8, Additional File 9). The requirements for fermentation, in terms of biochemical pathways, for all sourdoughs are related to carbohydrates, pyruvate and energy production and conversion, nitrogen (proteolytic systems, transport and catabolism of amino acids) and stress response. The majority of genes that we found in our sourdough metagenomics data are related to these main pathways.
Players actively engaged in sourdough metabolism
With the aim of identifying metabolically active species, we performed a metatranscriptomics experiment. Transcript profiles and comparison with metagenomics data were performed only for those species with a relative abundance higher than 0.1% in at least 1 sourdough. As a result, all metagenomics data from Lactobacillus (former taxonomy) and Weissella species, were confirmed in terms of taxa viability and relative expression profiles by culturomics and metatranscriptomics experiments (Fig. 2). Notably, species of Lactococcus and Enterococcus, as identified by metagenomics, were culturable but not transcriptionally active. Moreover, all Leuconostoc, Oenococcus and Pediococcus species identified by metagenomics were transcriptionally active, and almost all were culturable. Conversely, most of the other nonlactic acid bacteria were transcriptionally inactive and nonculturable. Staphylococcus aureus, Staphylococcus epidermidis and Staphylococcus warneri were culturable but transcriptionally inactive. Except in SD88, Sac. cerevisiae was always transcriptionally active and cultivable. Some yeast species were transcriptionally inactive and/or nonculturable (Fig. 2).
Reconstruction of the main metabolic pathways underlying player functionality
We aimed to capture the potential contribution of those players involved in the metabolic pathways that have a key role in sourdough functionality and thus impact resilience and stability. The investigation of genes and transcripts related to the main pathways relevant to sourdough biotechnological properties allowed us to determine the contributions of species belonging to dominant, subdominant and satellite groups (Supplementary Figure S2, Additional File 10). Starch, nonstarch polysaccharide and sucrose pathway core genes (shared by all sourdoughs), together with those heterogeneously represented by a few species (accessory genes), are listed in Table 3. More specifically, higher gene and transcript numbers were harbored by core dominant species, such as Fru. sanfranciscensis, Lim. fermentum, Lac. plantarum, Leuc. citreum and Weissella confusa and Sac. cerevisiae, whereas subdominant and dispensable species contributed fewer genes and transcripts.
The pathways involved in hexose fermentation (Embden-Meyerhof-Parnas and phosphoketolase pathways), pentose metabolism (phosphoketolase pathways), and the alternative fate of pyruvate are shown in Table 4.. Additionally, in this case, the selected core metagenome encoded all transcripts needed to ferment hexoses and pentoses (Table 4.). With few exceptions, all the basic functions at both the genomic and transcriptomic levels are performed by the dominant species. Among subdominant species, most of the genes for the pentose phosphate pathway were found in Lev. brevis, Fur. rossiae, Ped. pentosaceus, and/or Pic. kudriavzevii genomes. Key genes, such as transaldolase (EC:18.104.22.168) and transketolase (EC:22.214.171.124), were detectable in all sourdoughs. The genes encoding enzymes to synthesize acetate were detected in the dominant species of each sourdough, except for Fru. sanfranciscensis. The acetate pathway was also present in subdominant species Fur. rossiae, Lev brevis, Ped. pentosaceus and Leu. mesenteroides. Some strains of Sac. cerevisiae harbored the gene encoding acetyl-CoA synthetase (EC:126.96.36.199). This activity was found in all sourdoughs except for SD102. Additionally, for this pathway, the core dominant group contributed the highest number of genes and relative transcripts (Table 4.).
All sourdough metagenomes contained the bacterial and yeast genes responsible for peptide catabolism (Supplementary Table S9, Additional File 11). As evidenced by the species abundance in each sourdough, gene redundancy markedly varied in terms of transcript number. The core metagenome comprises a pattern of peptidases with different substrate specificities, such as 5 endopeptidases, 4 aminopeptidases, a tripeptide aminopeptidase and 3 dipeptidases, whose genes were expressed in all sourdoughs. The only exception was found for lactocepin EC:188.8.131.52. The core metagenome of all sourdoughs also harbored several proline-specific peptidases. Various other accessory peptidases were identified. Overall, Sac. cerevisiae and Pic. kudriavzevii have few proline-associated peptidases. The genes and transcripts relevant to the catabolism of branched chain amino acids (BCAAs), aromatic amino acids (ArAs) and free amino acids (FAAs) are listed (Supplementary Table S10, Additional File 12). Core dominants showed the highest gene and transcript contributions. The branched chain amino transferase (BcaT EC:184.108.40.206) was encoded by Lac. plantarum, Lim. fermentum, Leuc. citreum and one yeast (Sac. cerevisiae). Aminotransferases for aromatic amino acids (ArAT EC:220.127.116.11, EC:18.104.22.168 and EC:22.214.171.124) were only related to Sac. cerevisiae and/or Pic. kudriavzevii. Genes encoding L-asparaginase (EC:126.96.36.199) were widely distributed and active within the genome of dominant and subdominant lactic acid bacteria. Genes encoding arginine deiminase [EC:188.8.131.52], which is relevant to arginine degradation, were found to be active within the genome of Lac. plantarum, Leuc. citreum, Wei. confusa and Ped. pentosaceus.
De novo design of synthetic microbial communities
We de novo reconstructed a synthetic microbial community (Sourdough Global, SDG) by merging the results from all meta-omics. The biobank was used as the resource of the dominant and subdominant species, which were eligible to be included if they harbored at least 20 key genes and transcripts for each of the above pathways and were present in 4 out of 8 spontaneous sourdoughs. Satellite members were selected based on their frequency of isolation. The SDG comprised 7 species (Table 1). The factorial approach we used allowed us to investigate whether the depletion of one species at a time affected the functions of the remaining members, implying sourdough performance and stability.
Players contributing to the metabolic resilience of de novo SDG
The dominant KEGG function number for Lac. plantarum and Lim. fermentum did not vary between the exponential and stationary phases of growth, and it was not influenced by the depletion of a single species within the SDG (Fig. 3). Although the absolute value of KEGG transcripts was lower than that of the dominant species, the same stability was also observed for the subdominant Ped. pentosaceus. In contrast, the KEGG numbers for other bacteria and yeasts decreased upon reaching the stationary phase of growth and showed differential expression depending on the SDG composition. Species interactions modified the KEGG profiles when any single species was removed from the pools. The number of KEGG functions of the subdominant Fur. rossiae was driven by dominant species: the highest number of transcripts was observed in the pool without Lim. fermentum and Lac. plantarum. The opposite trend was observed during the stationary phase of growth without Sta. epidermidis or Ped. pentosaceus. The satellite Sta. epidermidis had the lowest KEGG transcript number when Fur. rossiae was absent and, more generally, showed a decreased number of KEGG functions when entering the stationary phase of growth. Nevertheless, its contribution during the exponential phase of growth was not negligible. The core dominant Sac. cerevisiae and the dispensable dominant Pic. kudriavzevii had the highest number of KEGG transcripts during the exponential phase of growth and when Sta. epidermidis and Ped. pentosaceus were absent (Fig. 3). Sac. cerevisiae had a higher number of transcripts than Pic. kudriavzevii in all SMCs.
The SDG had the capability to mimic all the functions needed to express carbohydrate metabolism. The potential contribution of each species to the starch, nonstarch polysaccharide and sucrose pathways was reconstructed (Fig. 4a, b). Key enzymes for starch, cyclodextrin and maltodextrin hydrolysis were encoded by core and dispensable dominant species of SDG. Lac. plantarum was the only species involved in the pathways leading to maltose formation (e.g., alpha-amylase, EC:184.108.40.206 and cyclomaltodextrinase/neopullulanase, EC:220.127.116.11), while the other two core dominant species (Lim. fermentum and Sac. cerevisiae) showed the ability to modify the starch/amylose structure by releasing D-glucose-1P. Concomitant activity on the D-glucose-1-P pathway was also detected for the dispensable dominant Pic. kudriavzevii.
Sac. cerevisiae was the only species capable of directly hydrolyzing starch in D-glucose (glucoamylase, EC:18.104.22.168). The core subdominant Fur. rossiae did not possess starch-hydrolyzing enzymes but showed the capability to ferment dextrin/isomaltose and maltose through oligo-1,6-glucosidase (EC:22.214.171.124) and alpha-glucosidase (EC:126.96.36.199). β-D-Glucosidase (EC:188.8.131.52) gene transcripts were unique to Fur. rossiae. Arabinoxylan metabolism leading to D-xylose and arabinose was exclusively found in Fur. rossiae.
Maltose phosphorylase (EC:184.108.40.206), glucokinase (EC:220.127.116.11) and β-phosphoglucomutase (EC:18.104.22.168) were encoded by all core dominant and subdominant lactic acid bacteria. Concomitantly, Sac. cerevisiae contributed to the metabolism of dextrin via maltose phosphorylase and uniquely showed the capability to catabolize maltose by amylomaltase (E. C 22.214.171.124). Lac. plantarum and Ped. pentosaceus showed the capability to use alternative energy sources such as trehalose and cellobiose via trehalose phosphorylase (EC:126.96.36.199) (only Lac. plantarum), 6-phospho-beta-glucosidase (EC:188.8.131.52) and trehalose-6-phosphate hydrolase (EC:184.108.40.206). The gene encoding β-fructofuranosidase (EC:220.127.116.11), which participates in fructan metabolism, was detected in all the core dominant species and Sta. epidermidis. Fructokinase (EC:18.104.22.168), which metabolizes fructose, was shared by all lactic acid bacteria. Within the fructose pathway, supplementary activity was shown by yeasts through hexokinase (EC:22.214.171.124). As observed downstream of the UDP-glucose intermediate, the contribution of Sac. cerevisiae was often redundant with that of Pic. kudriavzevii. The synergistic activity of endo-xylanase (EC:126.96.36.199), β-d-xylosidase (EC:188.8.131.52), β-glucosidase (EC:184.108.40.206) and α-arabinofuranosidase (EC:220.127.116.11) was mainly expressed by Fur. rossiae (Fig. 4b). Regarding the fermentation of hexoses and pentoses and the alternative fate of pyruvate (Fig. 4b), the genes encoding members of the pyruvate dehydrogenase complex [pyruvate dehydrogenase (EC:18.104.22.168), dihydrolipoamide dehydrogenase (EC:22.214.171.124), pyruvate dihydrolipoamide acetyltransferase (EC:126.96.36.199)] were detected in core dominant and subdominant strains as well as in Sta. epidermidis. This latter also provided a supplementary contribution with D-lactate dehydrogenase (EC:188.8.131.52). L-lactate dehydrogenase (cytochrome) (EC:184.108.40.206), D-lactate dehydrogenase (cytochrome) (EC:220.127.116.11) and (R)-2-hydroxyglutarate-pyruvate transhydrogenase (EC:18.104.22.168) activities were detected only in yeasts. Yeasts and Sta. epidermidis encompassed the encoding gene for acetyl-CoA synthetase (EC:22.214.171.124), an enzyme using acetate to synthesize acetyl-CoA. Another alternative pathway leading to the synthesis of acetate was detected. This included pyruvate decarboxylase (EC:126.96.36.199) and aldehyde dehydrogenase NAD( +)/NADP( +) (EC:188.8.131.52/184.108.40.206). Sac. cerevisiae encodes a mannitol 2-dehydrogenase (EC:220.127.116.11) that reduces D-fructose to D-mannitol. An alternative mechanism for NAD+ regeneration was found in Lac. plantarum, involving mannitol-1-phosphate 5-dehydrogenase (EC:18.104.22.168), which reduces D-fructose 6-phosphate to D-mannitol 1-phosphate.
Transcriptomic evidence of unique species-specific contributions also emerged for nitrogen metabolism (Fig. 5). The contribution of Sta. epidermidis is mainly related to tripeptide aminopeptidase (EC:22.214.171.124) and Xaa-pro aminopeptidase (EC:126.96.36.199). Despite the few genes transcribed, this satellite species showed higher expression levels than the other SDG species (Fig. 5a and c), whereas its global contribution was lowest (Fig. 5b). The SDG analysis showed that most nitrogen metabolism functions were shared among bacteria, while this redundancy might be extended to yeasts for only some pathways. Analysis of the BCAA metabolism pathway revealed that various enzymes were actively transcribed only in Sac. cerevisiae and Pic. kudriavzevii, which, together with Sta. epidermidis, contributed to several indispensable pathways. This is the case for glutamate dehydrogenase (EC:188.8.131.52) and aldehyde dehydrogenase (NAD +) (EC:184.108.40.206) (Fig. 4c). Gdh was encoded only by Lac. plantarum, Sta. epidermidis and yeasts. Aspartate ammonia-lyase (EC:220.127.116.11), the enzyme that converts aspartate into fumarate, was uniquely transcribed by Lac. plantarum (Fig. 4d). The 2 routes for arginine degradation, the arginine-urease and ADI pathways, were detected. The gene encoding arginase is present in Lim. fermentum, Sac. cerevisiae and Pic. kudriavzevii. Clustered genes for the ADI pathway were found in Lim. fermentum, Fur. rossiae and Ped. pentosaceus.
The de novo synthetic community shows transcriptome redundancies
Together with de novo SDG, we reconstructed another synthetic microbial community (SMC-SD43) to compare transcriptomics data in exponential and stationary growth conditions. SMC-SD43 comprised all bacterial and yeast species that were detected by culturomics in the spontaneous sourdough SD43. Both SDG and SMC-SD43 were cultivated in WFH, and their levels of transcripts assigned to KEGG functions (including copies of each transcript) were assessed (Fig. 6 a and b). The data focused on the robustness of the SDG with respect to SMC-SD43. The number of total KEGG transcripts in SDG was 4 times higher than that in SMC-SD43. This overall trend was also confirmed for each of the main reconstructed pathways. The highest number of gene transcripts in SDG belonged to carbohydrate metabolism, which was represented by more than 2500 gene copies.
The de novo synthetic community shows steady state under in situ conditions
Next, we assessed SDG and SMC-SD43 during 30 days of daily back slopping under in situ sourdough conditions. Both stabilized at similar pH values (4.3 ± 0.03 and 4.2 ± 0.02) after 5 days of propagation. The kinetics of acidification remained almost constant throughout 30 days. The cell densities of presumptive lactic acid bacteria were also steady, reaching similar values (ca. 8.51 to 8.65 log cfu/g). Presumptive staphylococci and micrococci (ca. 3.0 log cfu/g) and enterobacteria (ca. 3.9 to 4.9 log cfu/g) did not vary significantly. Yeasts slightly differed, at ca. 6.82 to 6.86 ± 0.02 log cfu/g for SMC-SD43 and 7.3 to 7.70 ± 0.04 log cfu/g for SDG. Biotype profiling (Fig. 7b and c) showed the sudden disappearance of Fru. sanfranciscensis, Lac. lactis and Staphylococcus sp. strains from SMC-SD43. After 10 days of propagation, other biotype members, such as Lac. paracasei and Lac. rhamnosus, were no longer detectable. At 30 days, SMC-SD43 harbored only 4 of the initial strains. In contrast, all strains inoculated in SDG persisted throughout back sloping (Fig. 7a). The only exception was Pic. kudriavzevii, which was no longer detectable after 10 days of propagation. The GC–MS data showed a profile of fifty-four VOCs, which reflected the stability of the sourdough metacommunities. Indeed, the levels of VOCs varied substantially during the propagation of SMC-SD43, while they remained almost constant for SDG. The statistical comparison (Welch test with BH correction) between sampling revealed variations of 45 VOCs in SMC-SD43 and only 3 VOCs in SDG (Supplementary Table S11, Additional File 13). Based on the linear distance of the principal component analysis plot, the two sampling times of SDG almost overlapped in the same quadrant (third), while those from SMC-SD43 were plotted in the second and fourth quadrants (Fig. 8).
Deterministic drivers render the sourdough microbiome/fermentome highly variable, dynamic and often unpredictable in terms of evolution, stability and performance . As further evidence of this, the representative spontaneous sourdoughs used in this study had a heterogeneous richness in culturable bacterial and yeast species that could not be appreciated simply by evaluating the cell densities of dominant members and the ratio between lactic acid bacteria and yeasts .
As we discussed in an earlier systematic review , lactic acid bacteria usually act as dominant players in sourdoughs, but we have also been able to detect several subdominant species, which together accounted for the majority of genes/transcripts. Deterministic drivers during sourdough back slopping promote the dominance of competitive species with ecological fitness over a pool of less abundant members, the role of which must not be neglected . Overall, we confirmed the constant dominance of Sac. cerevisiae among the yeasts . Acting as satellite members, other bacterial genera (e.g., staphylococci and few enterobacteria) completed the sourdough metacommunity. Repeated back slopping provides the opportunity for even unlikely contaminants to become significant members of the sourdough microbiome .
Deciphering the functional redundancy of sourdough genomes and the main players actively engaged in the main metabolic pathways
Sourdoughs shared a very high number of genes in common, although numerous unique genes were also assigned to unique sourdoughs. Biochemical pathways relevant to the main properties of sourdoughs, such as carbohydrate, pyruvate, energy and nitrogen metabolism, accounted for the highest number of genes. Based on metatranscriptomics, we might assert that higher microbial richness is accompanied by higher sourdough gene redundancy. The contribution of taxonomically close members was redundant in terms of functions . This evidence confirmed the findings from other microbial ecosystems, such as those of the human gut , animals , foods  and plants . We succeeded in assigning almost all functions to reconstructed pathways, which were fundamental for the resilience and performance of the sourdough. At the same time, we underlined the individual contributions of many players within the metacommunity that are actively part of dominant, subdominant and satellite groups. As already reported when describing the interactions between Lac. plantarum and Sac. cerevisiae , we demonstrated, in a more complex metacommunity, how single enzyme activities might rely on unique or multiple contributions from detectable species. All Weissella spp. and Lactobacillus spp. (former taxonomic nomenclature) detectable at the metagenomics level were also transcriptionally active in sourdoughs, whereas Lactococcus and Enterococcus were not. Although frequently harbored by spontaneous sourdoughs, these last two mentioned genera might have been displaced by lactobacilli during back slopping . Some other culturable bacteria (e.g., Enterobacteriaceae) also did not exhibit transcripts, probably because of acid inhibition and/or their inability to transcribe at low cell densities. Staphylococcus spp. are commonly isolated from cereal flours, but they do not show sufficient competitiveness for maintenance during sourdough back slopping . Nevertheless, when Sta. epidermidis was inoculated at higher cell densities and grown in the WFH model medium, it expressed gene transcripts, even within a heterogeneous metacommunity. To complete the critical analysis of this finding, we must acknowledge that by capturing single snapshots at specific timepoints, metatranscriptomics may suffer from specific biases relevant to transcript detection and quantification. Overall, the dominant and subdominant taxa were distinguished by the difference in numbers of genes and transcripts. Fru. sanfranciscensis, Lac. plantarum, Lim. fermentum, Ped. pentosaceus, Leu. citreum and Wei. confusa actively contributed to almost all reactions of all reconstructed pathways. This set of species was shared by several spontaneous sourdoughs, in particular those with the highest genetic potential and numbers of transcripts. Key enzymes such as fructose-1,6-bisphosphate aldolase (EMP pathway) and phosphoketolase (phosphogluconate pathway) are broadly distributed within the sourdough metatranscriptomes [54, 55]. Phosphoketolase was invariably present, which also indicated a role for anabolic reactions. NADH-dependent mannitol dehydrogenase activity contributing to additional cofactor regeneration was found only in Sac. cerevisiae. Since this enzyme is also pivotal for lactic acid bacteria competitiveness, the lack of the encoding gene might be ascribed to its high sequence similarity with alcohol- and L-threonine dehydrogenases . In addition, all sourdoughs harbored the gene encoding mannitol-1-phosphate 5-dehydrogenase, which was exclusively ascribed to Lac. plantarum . Specific reactions leading to the synthesis or release of amylose, 1–3 β-glucans and trehalose were exclusively contributed by Sac. cerevisiae and P. kudriavzevii. As confirmed by constructing the synthetic microbial community, Sac. cerevisiae was also unique in contributing to starch breakdown through amyloglucosidase activity (EC: 18.104.22.168). Another parallel pathway leading to the synthesis of D-glucose from β-D-glucoside revealed the unique contribution of Fur. rossiae.
De novo design of synthetic microbial communities and their players contributing to metabolic resilience
Based on the taxonomic classification, functional annotation and pathway inspection, we established the groundwork for selecting those species most suitable to reconstruction of a synthetic resilient and high-performing sourdough metacommunity. According to our hypotheses, subdominant genes synergistically cooperate with dominant players to express key genes, and satellite members are equipped with unique genes that are also important for the stability of the metacommunity. In summary, gene complementarity and redundancy are the main prerequisites to guarantee sourdough efficiency. The engineering of a multispecies consortium relies on a simplified and controlled environment , and therefore, WFH  was chosen as the model system. To unravel the contribution of each player to the de novo SDGs, the metatranscriptomics assessment was repeated. All the pathways reconstructed from spontaneous sourdoughs were confirmed, and the roles of individual players were deepened. Growth phases and depletion of single members from SDG did not influence the number of transcripts of Lac. plantarum, Lim. fermentum and Ped. pentosaceus, which contributed redundant functions, particularly for carbohydrate pathways. These species constitutively expressed enzymes for maltose utilization (e.g., maltose phosphorylase, β-phosphoglucomutase and glucokinase), which indicated an adaptation to the sourdough ecosystem, as maltose is the most abundant fermentable carbohydrate in cereal flours, but also a physiology and ecology that matched their phylogenetic placement . The comparison of exponential and stationary phases of growth revealed an adaptation in terms of species contribution. Specifically, when entering the stationary phase of growth, the number of transcripts from Sac. cerevisiae and Pic. kudriavzevii decreased; this behavior might be related to a quiescence state because of the lack of carbon sources . Moreover, slightly influenced by the phase of growth, the number of transcripts of Fur. rossiae increased when the dominant players were depleted from SDG. Fur. rossiae was the only species contributing to some pathways (e.g., arabinoxylans) and enzyme activities (e.g., β-D-glucosidase), which are fundamental for the rheologic, sensory and nutritional attributes of sourdough [60, 61]. Even though staphylococci are not competitive enough to grow under sourdough conditions , we highlighted some fundamental contributions from these bacteria, mainly to nitrogen metabolism. Staphylococci have already been used as alternatives to lactobacilli  to leverage their nitrogen metabolism (e.g., converting arginine into ornithine)  and, in general, for the uptake and release of free amino acids.
The de novo synthetic community shows steady state under in situ conditions
Compared to another reconstructed synthetic microbial community (SMC-SD43), which mimicked the composition of a spontaneous sourdough and had a higher number of bacterial and yeast species (9 vs. 7), the SDG more reliably maintained all the functions needed to resiliently express carbohydrate and nitrogen metabolism. The robustness of the SDG was attributable to its higher gene expression in terms of total KEGG pathways and copy numbers. The last proof to confirm our hypotheses was based on the in situ resilience and performance, which we monitored under perturbations that occurred during the daily sourdough back slopping. As expected, the two synthetic microbial communities did not show variation in terms of acidification or numbers of presumptive lactic acid bacteria and yeasts. It has been widely demonstrated that contaminant species may replace the initial ones and that these qualitative changes within the metacommunity are not observable using simple determinations such as pH and plating . Here, much more finely tuned biotype profiling revealed how the metacommunity of SDG remained stable during 30 days of propagation, while that of SMC-SD43 lost most of its species members. These changes might affect sourdough performance in terms of flavor, where VOC synthesis may be a suitable indicator. In fact, the VOC profiles of SDG overlapped during 30 days of propagation, while those of SMC-SD43 drastically changed.
In our complex workflow, we envisioned the sourdough as a specialized social structure from which members colonize the new flour environment in a complementary manner that is orchestrated strategically to cover all accompanying metabolic pathways. Our study demonstrates how, by starting from spontaneous sourdoughs and reconstructing the synthetic community, it was possible to unravel the metabolic contributions of individual players. For resilience and good performance, the sourdough metacommunity needed to include dominant, subdominant and satellite players, which together ensured gene and transcript redundancy. Overall, our study changes the paradigm and introduced theoretical foundations for directing food fermentations.
The BioProject, including metagenomics and metatranscriptomics samples, is under submission.
Availability of data and materials
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Friedman J, Higgins LM, Gore J. Community structure follows simple assembly rules in microbial microcosms. Nat Ecol Evol. 2017;1:109.
novel avenues for bio-based processes. Diender, M., Parera Olm, I. & Sousa, D. Z. Synthetic co-cultures. Curr Opin Biotechnol. 2021;67:72–9.
Antoniewicz MR. A guide to deciphering microbial interactions and metabolic fluxes in microbiome communities. Curr Opin Biotechnol. 2020;64:230–7.
Comasio A, Verce M, Van Kerrebroeck S, De Vuyst L. Diverse Microbial Composition of Sourdoughs From Different Origins. Front Microbiol. 2020;11:1212.
Arora K, et al. Thirty years of knowledge on sourdough fermentation: A systematic review. Trends Food Sci Technol. 2021;108:71–83.
Minervini F, Dinardo FR, Celano G, De Angelis M, Gobbetti M. Lactic Acid Bacterium Population Dynamics in Artisan Sourdoughs Over One Year of Daily Propagations Is Mainly Driven by Flour Microbiota and Nutrients. Front Microbiol. 2018;9:1984.
Grosskopf T, Soyer OS. Synthetic microbial communities. Curr Opin Microbiol. 2014;18:72–7.
De Angelis M, et al. Selection of potential probiotic lactobacilli from pig feces to be used as additives in pelleted feeding. Res Microbiol. 2006;157:792–801.
Kurtzman CP, Robnett CJ. Identification and phylogeny of ascomycetous yeasts from analysis of nuclear large subunit (26S) ribosomal DNA partial sequences. Antonie Van Leeuwenhoek. 1998;73:331–71.
Lugli GA, et al. Compositional assessment of bacterial communities in probiotic supplements by means of metagenomic techniques. Int J Food Microbiol. 2019;294:1–9.
Wagner Mackenzie B, Waite DW, Taylor MW. Evaluating variation in human gut microbiota profiles due to DNA extraction method and inter-subject differences. Front Microbiol. 2015;6:130.
Fuxman Bass, J. I., Reece-Hoyes, J. S. & Walhout, A. J. M. Zymolyase-Treatment and Polymerase Chain Reaction Amplification from Genomic and Plasmid Templates from Yeast. Cold Spring Harb Protoc 2016, (2016).
Tamames J, Puente-Sánchez F. SqueezeMeta, A Highly Portable, Fully Automatic Metagenomic Analysis Pipeline. Front Microbiol. 2018;9:3349.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.
Muggli MD, et al. Succinct colored de Bruijn graphs. Bioinformatics. 2017;33:3181–7.
Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
Laslett D, Canback B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res. 2004;32:11–6.
Hyatt D, et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.
Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2016;44:D67-72.
Huerta-Cepas, J. et al. eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res 44, D286–293 (2016).
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.
A new generation of homology search tools based on probabilistic inference - PubMed. https://pubmed.ncbi.nlm.nih.gov/20180275/.
Finn RD, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44:D279-285.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Wu, Y.-W., Simmons, B. A. & Singer, S. W. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 32, 605–607 (2016).
Kang DD, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7: e7359.
Sieber CMK, et al. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat Microbiol. 2018;3:836–43.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.
Ye Y, Doak TG. A parsimony approach to biological pathway reconstruction/inference for genomes and metagenomes. PLoS Comput Biol. 2009;5: e1000465.
Zapparoli G, Torriani S, Pesente P, Dellaglio F. Design and evaluation of malolactic enzyme gene targeted primers for rapid identification and detection of Oenococcus oeni in wine. Lett Appl Microbiol. 1998;27:243–6.
Di Cagno R, et al. Cell-cell communication in sourdough lactic acid bacteria: a proteomic study in Lactobacillus sanfranciscensis CB1. Proteomics. 2007;7:2430–46.
Zwietering, M. H., Jongenburger, I., Rombouts, F. M. & van ’t Riet, K. Modeling of the bacterial growth curve. Appl Environ Microbiol 56, 1875–1881 (1990).
Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30:3123–4.
Liu G, et al. Headspace solid-phase microextraction of semi-volatile ultraviolet filters based on a superhydrophobic metal-organic framework stable in high-temperature steam. Talanta. 2020;219: 121175.
Perri G, et al. Bioprocessing of Barley and Lentil Grains to Obtain In Situ Synthesis of Exopolysaccharides and Composite Wheat Bread with Improved Texture and Health Properties. Foods. 2021;10:1489.
Ercolini D, et al. Microbial Ecology Dynamics during Rye and Wheat Sourdough Preparation. Appl Environ Microbiol. 2013;79:7827–36.
Michel E, et al. Characterization of relative abundance of lactic acid bacteria species in French organic sourdough by cultural, qPCR and MiSeq high-throughput sequencing methods. Int J Food Microbiol. 2016;239:35–43.
Menezes, L. a. A. et al. Sourdough bacterial dynamics revealed by metagenomic analysis in Brazil. Food Microbiol 85, 103302 (2020).
Oshiro M, Tanaka M, Momoda R, Zendo T, Nakayama J. Mechanistic Insight into Yeast Bloom in a Lactic Acid Bacteria Relaying-Community in the Start of Sourdough Microbiota Evolution. Microbiol Spectr. 2021;9: e0066221.
Gänzle M, Ripari V. Composition and function of sourdough microbiota: From ecological theory to bread quality. Int J Food Microbiol. 2016;239:19–25.
Lau SW, Chong AQ, Chin NL, Talib RA, Basha RK. Sourdough Microbiome Comparison and Benefits Microorganisms. 2021;9:1355.
Gänzle MG, Zheng J. Lifestyles of sourdough lactobacilli - Do they matter for microbial ecology and bread quality? Int J Food Microbiol. 2019;302:15–23.
Louca S, et al. Function and functional redundancy in microbial systems. Nat Ecol Evol. 2018;2:936–43.
Moya A, Ferrer M. Functional Redundancy-Induced Stability of Gut Microbiota Subjected to Disturbance. Trends Microbiol. 2016;24:402–13.
Estrada-Peña A, Cabezas-Cruz A, Obregón D. Behind Taxonomic Variability: The Functional Redundancy in the Tick Microbiome. Microorganisms. 2020;8:1829.
Kothe CI, et al. Unraveling the world of halophilic and halotolerant bacteria in cheese by combining cultural, genomic and metagenomic approaches. Int J Food Microbiol. 2021;358: 109312.
Harbort CJ, et al. Root-Secreted Coumarins and the Microbiota Interact to Improve Iron Nutrition in Arabidopsis. Cell Host Microbe. 2020;28:825-837.e6.
Zhang G, et al. Proteomic Analysis Explores Interactions between Lactiplantibacillus plantarum and Saccharomyces cerevisiae during Sourdough Fermentation. Microorganisms. 2021;9:2353.
De Vuyst L, Van Kerrebroeck S, Leroy F. Microbial Ecology and Process Technology of Sourdough Fermentation. Adv Appl Microbiol. 2017;100:49–160.
Weckx S, et al. Community dynamics of bacteria in sourdough fermentations as revealed by their metatranscriptome. Appl Environ Microbiol. 2010;76:5402–8.
De Vuyst, L., Comasio, A. & Kerrebroeck, S. V. Sourdough production: fermentation strategies, microbial ecology, and use of non-flour ingredients. Crit Rev Food Sci Nutr 1–33 (2021) doi:https://doi.org/10.1080/10408398.2021.1976100.
Weckx S, Van Kerrebroeck S, De Vuyst L. Omics approaches to understand sourdough fermentation processes. Int J Food Microbiol. 2019;302:90–102.
Zhang Y, Kastman EK, Guasto JS, Wolfe BE. Fungal networks shape dynamics of bacterial dispersal and community assembly in cheese rind microbiomes. Nat Commun. 2018;9:336.
Zheng J, Ruan L, Sun M, Gänzle M. A Genomic View of Lactobacilli and Pediococci Demonstrates that Phylogeny Matches Ecology and Physiology. Appl Environ Microbiol. 2015;81:7233–43.
Martinez MJ, et al. Genomic analysis of stationary-phase and exit in Saccharomyces cerevisiae: gene expression and identification of novel essential genes. Mol Biol Cell. 2004;15:5295–305.
Pontonio E, et al. Cloning, expression and characterization of a β-D-xylosidase from Lactobacillus rossiae DSM 15814(T). Microb Cell Fact. 2016;15:72.
Son S-H, et al. Probiotic lactic acid bacteria isolated from traditional Korean fermented foods based on β-glucosidase activity. Food Sci Biotechnol. 2018;27:123–9.
De Vuyst L, Van Kerrebroeck S, Leroy F. Microbial Ecology and Process Technology of Sourdough Fermentation. Adv Appl Microbiol. 2017;100:49–160.
Sánchez Mainar, M., Weckx, S. & Leroy, F. Coagulase-negative Staphylococci favor conversion of arginine into ornithine despite a widespread genetic potential for nitric oxide synthase activity. Appl Environ Microbiol 80, 7741–7751 (2014).
Gobbetti M, Minervini F, Pontonio E, Di Cagno R, De Angelis M. Drivers for the establishment and composition of the sourdough lactic acid bacteria biota. Int J Food Microbiol. 2016;239:3–18.
This work was supported by the Open Access Publishing Fund of the Free University of Bozen-Bolzano. We thank the International Sourdough Library (Puratos, St. Vith, Belgium) for providing us with sourdough samples.
F.M.C.: Methodology, Investigation, Formal analysis, Data analysis, Writing – original draft. H.A.: Methodology, Investigation, Formal analysis, Data analysis, Writing – original draft. O.N.: Methodology, Investigation, Formal analysis, Writing – original draft. G.C.: Formal analysis, Writing – original draft. M.V.: Formal analysis. W.J.F.L.J.: Data analysis. C.M.: Formal analysis. F.V.: review & editing. R.D.C.: Supervision, Writing – review & editing. G.P.: review & editing. M.D.A.: Supervision, Writing – review & editing. M.G.: Funding acquisition, Conceptualization, Methodology, Supervision, Project administration, Writing – original draft, Writing – review & editing. All authors read and approved the final manuscript.
This work was supported by the Open Access Publishing Fund of the Free University of Bozen-Bolzano.
Ethics approval and consent to participate
Consent for publication
Fabienne Vertè is employed by Puratos NV. The other authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Origin, ingredients and technology parameters of sourdoughs.
Assembly statistics relative to the eight spontaneous sourdough assembled metagenomes.
Total number of genera and species of bacteria detected in the eight spontaneous sourdoughs by shotgun metagenomics.
Total number of genera and species of yeasts detected in the eight spontaneous sourdoughs by shotgun metagenome analysis.
Core and dispensable species found in the eight spontaneous sourdoughs by culturing, and shotgun meta-genomics and -transcriptomics. Pseudo-heatmap displaying the comparison between the core (species shared by the eight sourdoughs) and dispensable (species with relative abundance >0.1% variously detected within one or more sourdoughs) transcriptionally active microbiome species under sourdough conditions, and cultivable microbiota. Colour bar on the right describes species prevalence. Among omics C stands for culturomics, SM and MT stand for shotgun metagenomics and metatranscriptomics, respectively.
Total number of mapped reads in the eight spontaneous sourdoughs by metatranscriptomics analysis.
Metagenomic ORFs predicted by using Aragorn, Prodigal and Barrnap and relative ORF annotations based on the KEGG, COG and Pfam databases for each sourdough (SD).
Shared and unique KEGG genes from metagenomes in the eight spontaneous sourdoughs.
List of KEGG ID, gene name, function, transcript level (transcripts per million = TPM), copy number, total bases and coverage of each gene found in sourdough metagenomes.
Gene and transcript class content (core and dispensable dominants and subdominants) for aminotransferase a) deaminase/lyase b), carbohydrate c) and pyruvate d) microbial pathways in the 8 spontaneous sourdoughs.
Core and accessory genes encoding the enzymatic portfolio related to protein/peptide degradation in the eight spontaneous sourdoughs as estimated by metagenomic analyses.
Distribution of genes encoding the catabolism of branched chain amino acids (BCAAs), aromatic amino acids (ArAAs) and free amino acids. The expressed genes, evaluated through metatranscriptomic analyses, are marked by asterisks.
Statistically significant VOCs in SDGlobal and SMC-SD43. VOC standardized abundance whole matrices were used to compute a Welch corrected test (Benjamini–Hochberg). Mean rel. freq. = mean relative frequency; std.dev. = standard deviation; CI = confidence interval.
About this article
Cite this article
Calabrese, F.M., Ameur, H., Nikoloudaki, O. et al. Metabolic framework of spontaneous and synthetic sourdough metacommunities to reveal microbial players responsible for resilience and performance. Microbiome 10, 148 (2022). https://doi.org/10.1186/s40168-022-01301-3