Next Article in Journal
Global Genomic Analysis of SARS-CoV-2 RNA Dependent RNA Polymerase Evolution and Antiviral Drug Resistance
Next Article in Special Issue
Polyphenols Content and In Vitro α-Glycosidase Activity of Different Italian Monofloral Honeys, and Their Effect on Selected Pathogenic and Probiotic Bacteria
Previous Article in Journal
Novel High-Throughput Multiplex qPCRs for the Detection of Canine Vector-Borne Pathogens in the Asia-Pacific
Previous Article in Special Issue
Characterization of Soil Bacteria with Potential to Degrade Benzoate and Antagonistic to Fungal and Bacterial Phytopathogens
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 9
Issue 5
10.3390/microorganisms9051093
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Yeasts and Lactic Acid Bacteria for Panettone Production: An Assessment of Candidate Strains

by
Luciana De Vero
1,
Giovanna Iosca
1,
Salvatore La China
1,
Fabio Licciardello
1,2,
Maria Gullo
1 and
Andrea Pulvirenti
1,2,*
1
Department of Life Sciences, University of Modena and Reggio Emilia, 42122 Reggio Emilia, Italy
2
Interdepartmental Research Centre for the Improvement of Agri-Food Biological Resources (BIOGEST-SITEIA), University of Modena and Reggio Emilia, 42124 Reggio Emilia, Italy
*
Author to whom correspondence should be addressed.
Microorganisms 2021, 9(5), 1093; https://doi.org/10.3390/microorganisms9051093
Submission received: 3 May 2021 / Revised: 18 May 2021 / Accepted: 18 May 2021 / Published: 19 May 2021
(This article belongs to the Special Issue Selected Papers from 1st International Electronic Conference on Microbiology (ECM 2020))
Browse Figures
Versions Notes

Abstract

:
The recovery of yeasts and lactic acid bacteria (LAB) involved in sourdough fermentation is the first step in the selection of starters with suitable technological aptitude and capable of producing desired aromas and/or aromatic precursors. In this work, two sourdoughs samples (MA and MB) and the derived doughs (samples A and B) were collected from a bakery during artisanal Panettone manufacture. Yeasts and bacteria were isolated at different fermentation steps on selective agar media. A total of 77 isolates were obtained and characterized. Representative strains of yeasts and LAB were identified by sequencing the D1/D2 domain of the 26S rRNA and the 16S rRNA genes, respectively. Moreover, the volatile organic compounds (VOCs) produced in the collected samples were detected and correlated to the species found in the same samples. The results highlighted the occurrence of Kazachstania humilis in both samples A and B, while Saccharomyces cerevisiae strains were detected only in samples B. Among LAB, Fructilactobacillus sanfranciscensis was the main species detected in both sourdoughs. Furthermore, strains belonging to the species Lactiplantibacillus plantarum, Furfurilactobacillus rossiae, Lactobacillus parabuchneri, Leuconostoc citreum, and Leuconostoc mesenteroides were assessed in the dough samples.

1. Introduction

Panettone is a traditional Italian baked cake produced from sourdough, also called mother dough (MD). Sourdough can be defined as a matrix of flour and water including yeasts and lactic acid bacteria (LAB) as main functional microorganisms [1,2].
Commonly, the type I sourdough is specifically used for the artisanal production of Panettone and other traditional Italian sweet baked products, such as Pandoro and Colomba [3,4,5]. This type of sourdough is characterized by daily refreshments, carried out every 4–16 h at 25–35 °C according to a back-slopping procedure, which consists of the addition of flour and water with an aliquot of previously fermented dough [6,7]. This practice maintains the microorganisms in an active state and promotes the selection of a stable and characteristic microflora, which is well adapted to the specific recipe and manufacturing procedure [8,9].
The main species occurring in type I sourdoughs are the LAB Fructilactobacillus sanfranciscensis (formerly Lactobacillus sanfranciscensis), Levilactobacillus brevis (formerly Lactobacillus brevis), and Lactiplantibacillus plantarum (formerly Lactobacillus plantarum) and the yeasts Saccharomyces cerevisiae, Kazachstania exigua, and Kazachstania humilis [7,10,11].
Technological parameters, such as temperature, flour composition, degree of dough hydration, and sodium chloride content significantly contribute to the selection of the microflora [12]. Moreover, the richness and diversity of yeasts in sourdoughs can differ considerably as they are influenced by many factors including tolerance to the organic acids produced by the LAB and the availability of carbon sources [12,13,14].
Generally, LAB activity causes the acidification of sourdough, while yeasts are mainly responsible for the dough leavening [4,6]. In addition, the metabolism of LAB and yeasts contributes to the aroma’s formation through the production of important flavor compounds such as diacetyl, other carbonyls, ethyl acetate, and isoalcohols [7,15,16].
Although sourdoughs have been widely studied, currently, there is an interest in exploring their microbial composition, due to the numerous scientific studies that highlight their positive influence on sensory, nutritional, and shelf-life characteristics of naturally leavened products [17,18,19]. Additionally, the continuous search for starter cultures with new attributes, aimed at enhancing the nutritional and nutraceutical traits of sourdoughs, is driving the interest in exploring their microbial diversity [20,21,22].
In the present work, we investigated the yeasts and LAB population of sourdoughs for Panettone production, as well as the volatile fraction of dough samples, with the aim to select candidate strains for developing single and/or multiple starter cultures.

2. Materials and Methods

2.1. Panettone Production and Sample Collection

Two sourdoughs, defined as MA and MB, were used to produce an artisanal Panettone in a bakery (Panificio Fantuzzi) of the Emilia Romagna region (Italy). They differed in the kind of storage after the periodical refreshments; specifically, MA was kept in a cotton bag, while MB in a bucket with water. Starting from these, the refreshment procedures were made using only wheat flour (type 0) and water. Three back-slopping steps were made for MA and only two for MB. In both cases, about 28% (w/w) of the final refreshed sourdough was used for the first dough preparation by adding the following ingredients: flour (4.0 kg), sugar (1.25 kg), yolks (1.0 kg), butter (1.5 kg), and water (2.1 L). After 18 h of incubation at 20–22°C, a second dough was prepared by mixing the first dough with flour (1.5 kg), sugar (1.0 kg), yolks (1,0 kg), butter (2.0 kg), barley malt (0.05 kg), salt (0.065 kg), and water (0.25 L). The final dough was left to rise for 5–6 h before baking. No aromatic ingredients were added.
During the entire process, six samples were taken concurrently in sterile containers for the laboratory analyses. Specifically, the samples were MA and MB; the first doughs were named IMPA and IMPB, collected after the leavening of 18 h, and the final doughs were called FINA and FINB, collected just before the cooking. A schematic representation of the process with the sampling points is reported in Figure 1.

2.2. Determination of Physicochemical Parameters

The following parameters pH, total titratable acidity (TTA), and water activity (aw) were evaluated on all the samples collected. In detail, pH was detected with a pH meter XS series pH 70 (Bormac srl, Carpi, Italy), and TTA was determined on the homogenized samples and expressed as the amount (mL) of 0.1 N NaOH necessary to achieve pH 8.3 [23]. The aw was measured with the AquaLab 4TE instrument (Meter Group, Pullman, WA, USA) following the manufacturer’s instructions. Each sample was analyzed in triplicate.

2.3. Yeasts and LAB Isolation

To obtain microbial isolates, 10 g of each sample was collected and mixed with 90 mL of physiological solution (9 g/L of NaCl) in a Stomacher bag and homogenized for 2 min. From the homogenate samples, serial dilutions from 10-1 to 10-7 were made and appropriate dilutions were plated in triplicate on different agar media. Specifically, Yeast Peptone Dextrose Agar (YPDA; 10 g/L yeasts extract, 10 g/L peptone, 20 g/L dextrose, and 20 g/L agar) supplemented with 0.1 g/L of chloramphenicol (Sigma-Aldrich, Milan, Italy) was used for yeasts isolation. Plates were incubated at 27 °C for 2 days. LAB isolation was made on de Man Rogosa Sharpe (MRS agar; Oxoid, Milan, Italy) and Sourdough Bacteria Agar (SDB; 6 g/L trypticase, 3 g/L yeast extract, 20 g/L maltose, 3 mL of 10% Tween 80, 1.5 % of fresh yeast extract, 15 g/L agar), both supplemented with 0.1 g/L cycloheximide (Sigma-Aldrich). Plates were incubated in jars with the Anaerogen system (Oxoid) at 30 °C for 3 days. After the colonies were counted on the plates of each medium, between 3 and 10 colonies were randomly selected at the highest dilutions and purified by streaking on the correspondent isolation medium. Basic phenotypic tests, such as Gram staining and catalases, were done for bacteria. Cell morphology of the isolates was observed using a Zeiss Axiolab microscope (Carl Zeiss Ltd., Cambridge, UK). After the molecular characterizations, yeasts and bacteria strains were preserved in the Unimore Microbial Culture Collection (UMCC) in accordance with standard procedures [24]. A seed lot of the strains was stored at −80 °C, in cryovials containing an aliquot of each culture, grown in liquid medium, and 25% (v/v) glycerol (Carlo Erba, Milan, Italy).

2.4. DNA Extraction and Molecular Characterization

Genomic DNA from yeast and bacterial cultures was extracted following the protocols of Hoffman and Winston [25] and Ausubel et al. [26], respectively. The quality and concentration of the extracted DNA were checked by NanoDrop™ 1000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA).
For a preliminary characterization at species level, PCR-RFLP of the internal transcribed spacer (ITS) region and 5.8S gene of rDNA was performed on yeast gDNA following the protocol described by Pulvirenti et al. [14]. ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′), ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) oligonucleotide primers (Eurofins MWG Operon, Ebersberg, Germany), and PCR amplification kit, which included TaKaRa Taq DNA Polymerase (Takara Bio, Inc., Otsu, Shiga, Japan), were used. PCR was made on a final volume of 50 µL by using the thermal cycler BioRAD T100TM (Bio-Rad Laboratories, Milan, Italy). Amplicons digestion was performed with HaeIII restriction endonucleases (Fermentas, Hanover, ND, USA), according to supplier’s recommendation. The fragments were separated on 2% (w/v) agarose gel electrophoresis in 0.5X Tris-Borate-EDTA (TBE) buffer at 90 volts for 2 h. The gels were stained with SYBR Safe reagent (Thermo Fisher Scientific) and visualized on the Transilluminator Safe View (Cleaver Scientific LTD, Rugby, UK).
A first rapid detection of LAB strains belonging to Fl. sanfranciscensis species was made based on PCR amplification of the 16S–23S rDNA intergenic spacer region (ISR). PCR was performed according to the protocol described by Valcheva et al. [27] by using the oligonucleotide primers 16S p2: (5′-CTTGTACACACCGCCCGTC-3′) and 23S p7: (5′-GGTACTTAGATGTTTCAGTTC-3′) (Eurofins MWG Operon). The fragments were separated on 1.5% (w/v) agarose gel electrophoresis in 0.5X Tris-Borate-EDTA (TBE) buffer at 90 volts for 1.5 h.
For strain characterization, repetitive element sequence-based (rep)-PCR genotyping [28] using (GTG)5 oligonucleotide primer (Eurofins MWG Operon) was performed on gDNA of all the isolates according to the protocol described by La China et al. [29]. Fingerprinting profiles were assessed after electrophoresis on 1.5% (w/v) agarose gel, carried out as previously described. Pattern band lengths were determined by comparison against a 100 bp plus DNA ladder (Thermo Fisher Scientific).
Digitalized images of the electrophoretic profiles were analyzed by using the GelCompare software v3.0 (Applied Maths, Kortrijk, Belgium). The similarity matrix of the bands’ patterns was computed by using Pearson’s correlation with optimization and curve smoothening values at 0.5%. Dendrograms were constructed by the unweighted-pair group method using arithmetic averages (UPGMA) clustering method.

2.5. Species Assignment and Phylogenetic Clustering

At least one strain of yeasts and bacteria isolated from each sample was subjected to sequencing analysis for the species assignment and the successive phylogenetic clustering. To this aim, targeted regions were amplified using specific primer sets. In particular, the D1/D2 region of 26 rRNA large subunit from the yeast gDNA was amplified using primers NL1 (5′-GCATATCAATAAGCGGAGGAAAAG-3′) and NL4 (5′-GGTCCGTGTTTCAAGACGG-3′), as reported by Kurtzman et al. [30].
Yeast Ab1 files obtained from the sequencing were processed using CodonCode aligner to trim the sequences based on phred score (phred 20). The sequence ends were trimmed according to primer length. High-quality sequences were aligned against YeastIP database [31]. Reference sequences were downloaded from YeastIP database, based on matches, e-value, and identity scores.
Regarding bacteria, amplification of 16S rRNA gene was performed using primers 27f (5′-CTGGGATCCATTTACTCGAGAGTTTGATCCTGGCTCAG-3′) and 1490r (5′-GGTTCCCCTAAGCTTACCTTGTTACGACTTC-3′), as reported by Sato et al. [32].
PCR amplicons were purified and cleaned using ExoSAP PCR cleanup reagent (Thermo Fisher Scientific) and sequenced by Eurofins MWG Biotech Company (Ebersberg, Germany). The strains sequenced were deposited in the UMCC Database (https://umcc.bio-aware.com/, accessed on 20 April 2021).
Bacteria Ab1 files were processed as described for the yeast dataset. High-quality sequences were aligned against 16S rRNA sequence database from NCBI using blast algorithm. The top hit with the highest percentage of identities (considering a minimum threshold of 97%) was downloaded and used as reference.
The sequences of selected strains and references were aligned all-vs.-all using Clustal W v2.1 [33]. For constructing the phylogenetic trees, the multiple sequences alignment was imported into MegaX v10.2 and trimmed to match the sequence length. The phylogeny was inferred computing a neighbor-joining phylogenetic tree based on 1000 bootstraps for both yeasts and LAB dataset. The Tamura-Nei DNA evolutionary model was used, applying a discrete Gamma distribution to model evolutionary rate differences among sites [34].
Newick trees of both bacteria and yeasts were visualized using Interactive Tree of Life (ITOL) v4 [35] and rooted at outgroup reference strains. The sequences of D1/D2 region and 16S rRNA gene were submitted to GenBank/EMBL/DDBJ under the accession numbers from MZ170795 to MZ170810 and from MZ170701 to MZ170718, respectively.

2.6. HS-SPME-GC/MS Profiling of Fermented Products

Volatile organic compounds (VOCs) developed during dough fermentation were determined by headspace solid-phase microextraction (HS-SPME) followed by gas-chromatography/mass spectrometry (GC-MS) analysis.
Seven grams of the six samples collected were weighted into 25 mL screw-cap glass vials provided with Mininert© valves.
Vials were conditioned at 60 °C for 30 min in a thermoblock (Falc Instruments, Treviglio, Italy), and then a divinylbenzene/carboxen/polydimethylsiloxane (DVB/CAR/PDMS) SPME fiber was exposed in the headspace for 30 min at the same temperature for the extraction of volatile compounds. Chromatographic separation of analytes was carried out by an Agilent (Paolo Alto, CA, USA) 6890N GC followed by mass spectrometer Agilent 5973 Network Mass Selective Detector. After extraction, fibers were desorbed for 3 min into the GC injector port set in splitless mode at 240 °C. The GC carrier gas used was helium at 1 mL/min and the detector temperature set at 240 °C. GC oven temperature program was: start at 50 °C for 3 min, 5 °C/min until 160 °C, hold at 160 °C for 2 min, 20 °C/min until 240 °C, hold at 240 °C for 2 min. Peaks identification was carried out by comparison with system libraries (Wiley, Nist). The analyses were performed in triplicate. A correlation analysis of VOCs and species detected in sourdough and dough samples was obtained using the Corrplot v. 0.84 [36] package implemented in R v 4.0.3 [37].

2.7. Statistical Analysis

pH, temperature, TTA, and aw data were expressed as mean value ± standard deviation (sd). Microbial counts were converted to log colony-forming units (CFU) per g of sample and expressed as mean ± sd. The data were subjected to one-way analysis of variance (ANOVA) followed by Tukey HSD post hoc test to establish significant differences between means (p < 0.05; n = 3). For the statistical analysis, the software GraphPad Prism version 8 was used (GraphPad Software, San Diego, CA, USA).
Total VOCs in sourdough and dough samples were expressed as a mean value of the total chromatographic area detected by GC-MS analysis, while the relative abundance of each VOC was reported as % of the total chromatographic area ± sd.

3. Results

3.1. Samples Analysis and Microbial Counts

The physicochemical parameters detected on the samples collected during the Panettone production are reported in Table 1. The pH values of sourdoughs MA and MB pH were 3.90 and 3.81, respectively. A significant increase in the pH was observed in all the other samples because of the dilution effect of the ingredients added during the production process. Accordingly, the TTA was higher in the two sourdoughs and decreased in doughs IMPA and IMPB, as well as in the final doughs. The same trend was observed also for the water activity (aw), which reached a final value of 0.9441 and 0.9411 in FINA and FINB, respectively. In these cases, above the dilution effect, the addition of ingredients, such as sugars, also contributed to lowering the aw values [8].
Regarding the colony counts of the yeasts assessed on YPDA, it ranged from 6.81 to 7.83 log cfu/g for samples A and from 7.56 to 7.89 log cfu/g for B. The viable bacterial population enumerated on MRS and SDB media was generally higher in MA and MB, with values around 6 log cfu/g; in all the other samples, the values were around 5 log cfu/g (Table 2). In total, 57 yeast isolates were collected and further investigated. Concerning the isolates from MRS and SDB media, 20 were selected as presumptive LAB after the basic microbiological screenings.

3.2. Yeasts’ Molecular Characterization and Phylogenetic Analysis

The analysis of the RFLP-PCR of ITS-5.8S rDNA and the yeasts’ amplicons digestion with HaeIII allowed us to detect two different profiles (Table 3). Specifically, all the isolates belonging to samples A showed the same profile characterized by two bands of 230 and 400 base pairs (bp). As reported in the literature, this profile was tentatively attributed to the species K. humilis (formerly Candida milleri) [13,14]. Regarding the isolates from samples B, 12 showed a profile with two bands, as before, whereas a profile with three bands of 150, 180, 230, 320 bp were detected in the remaining ones. These isolates were preliminarily attributed to S. cerevisiae species according to Esteve-Zarzoso et al. [38].
For strain clusterization, the (GTG)5 rep-PCR fingerprinting techniques were performed. Figure 2 reports the dendrogram obtained from UPGMA analysis using Pearson’s coefficient. The discrimination of biotypes was performed considering a similarity threshold of 89%, which allowed the grouping of yeast isolates into 25 clades. The resulting discrimination power, calculated using the Simpson’s index, was 0.92. The clades with the highest number of strains (from 5 to 13) included those preliminarily attributed to K. humilis. Considering a similarity threshold of 40%, all the tentative K. humilis strains (42 in total) were grouped closely in a major clade clearly distinct from a minor clade that included the remaining 15 strains preliminarily attributed to S. cerevisiae.
A total of 43 high-quality sequences were used to infer the phylogeny of yeast isolates. Reference sequences belong to Saccharomyces and Kazachstania genera, detected after alignment in the YeastIP database. After all-vs.-all sequence alignment, sequences were trimmed all at the same length, resulting in a total length of 599 bp. The neighbor-joining phylogenetic tree was reported in Figure 3. The bootstrap values higher than 50%, used to define the accuracy of branch prediction, were reported.
According to sequence clusterization, the phylogenetic tree resulted in three well-defined clades, distinguishing with high accuracy the detected species, as reported by bootstrap values. The total average of dissimilarities among all phylotypes considered was 7.04%. The two major clades included Saccharomyces and Kazachstania species, supported by a bootstrap value of 98%. A total of four isolates, coming from sample B (LFB1, LFB10, LIB7, and LMB9), were grouped in S. cerevisiae clade, having a percentage of identity of 100% with S. cerevisiae CBS 1171T. The remaining Saccharomyces species were grouped in three different clades, of which one includes S. uvarum CBS 395T and the hybrid species from S. uvarum, S. bayanus CBS 380T, and S. pastorianus CBS 1538T. Another cluster of Saccharomyces species includes S. kudriavzevii CBS 8840T and S. arboricolus CBS 10644T, which share 99.2% of sequence identitt. The remaining clade of Saccharomyces includes S. cariocanus CBS 8841T, S. paradoxus CBS 432T, S. jurei NCYC D5088T, and S. mikatae CBS 8839T. Species included in this clade were phylogenetically closest, except for S. mikatae, sharing a dissimilarity percentage of about 1% with the other clade members.
The Kazachstania clade is branched out into three subclades. A clade including isolates coming from sample A and sample B was represented by K. humilis. Three isolates, represented by LMB1, LIB5, and LFA1, share 100% identity with K. humilis CBS 5658T and K. humilis H38. The remaining strains in the K. humilis clade, including C. milleri CBS 6897 (currently a synonym of K. humilis), are phylogenetically quite distant from K. humilis CBS 5658T, which shares the 99.8% similarity, indicating a low divergence as also reported from different studies [39,40]. Given this consideration, K. humilis and C. milleri are considered conspecific [41]. The most closely related to K. humilis strains is K. pseudohumilis CBS 11404T, with divergence supported by a bootstrap of 60%. The remaining Kazachstania strains included in the dataset were in a different clade, supported by a bootstrap of 86%. This clade includes K. barnettii CBS 5648 and K. barnettii CBS 6946T, forming a single group and the stains K. rupicola CBS 12684T, K. serrabonitensis CBS 14236T, K. australis CBS 2141T, and K. exigua CBS 379T. The represented species are low in divergence since the sequence differs by 0.7% to 1.7%.

3.3. Bacteria Molecular Characterization and Phylogenetic Analysis

The analysis of the 16S-23S rDNA ISR allowed to detect Fl. sanfranciscensis species among the bacterial isolates (Table 4). In particular, PCR products from all the MA and MB isolates yielded similar profiles containing three bands of 600, 700, and 800 bp; these isolates were preliminarily attributed to Fl. sanfranciscensis [27]. All the remaining isolates showed a profile with only two bands of 600 and 800 bp. This result was consistent with previously reported data for Lactobacillus 16S–23S ISR organization [27,42,43].
The digitized patterns obtained from (GTG)5 rep-PCR, analyzed with the UPGMA method and using Pearson’s coefficient, allowed the construction of the dendrogram shown in Figure 4. The biotypes discrimination was performed considering a similarity threshold of 95%, which allowed the grouping of the strains into 19 different clades; only the strains BIA5 and BIB8 were included in the same clade. The calculated Simpson’s index was 0.99, indicating a high diversity among the isolated strains including those preliminarily attributed to the species Fl. sanfranciscensis. This is consistent with the evidence reported by several authors that highlighted the great variability among F. sanfranciscensis strains [7,27].
High-quality sequences obtained from the bacteria dataset were analyzed in order to determine the phylotypes of isolates. After all-vs.-all alignment, sequences were trimmed at the same length, resulting in a total length of 904 bp. Reference sequences were downloaded from NCBI 16S rRNA database, resulting in two genera including Leuconostoc and Lactobacillus. The sequence of Bacillus albidus MCCC 1A02146T was used as an outgroup.
The average phylogenetic distance among all sequences included in the dataset was 13.17%. The resulting phylogenetic tree was represented in Figure 5. A total of two major clades are depicted, represented by the clade of Leuconostoc genus and the clade including the genera Furfurilactobacillus, Lactiplantibacillus, Lentilactobacillus, and Fructilactobacillus. In particular, Leuconostoc clade includes five isolates and is branched in two subclades, Leuc. mesenteroides NBRC 100469T and Leuc. citreum ATCC 49370T.
Two strains, BIA7 and BFB6, were grouped with Leuc. mesenteroides NBRC 100469T, sharing 100% similarity.
Three strains, BFA2, BFB8, and BIA6, were clustered with Leuc. citreum ATCC 49370T, sharing 100% of similarity among them and 99.9% of similarity with the reference strain. The phylotype classification was supported by high bootstrap values (ranging between 77% and 100%).
The remaining isolated strains were included in the other major clade, represented by four species: Furfurilactobacillus rossiae (formerly Lactobacillus rossiae) DSM 15814T, Lacp. plantarum JCM 1149T, Lentilactobacillus parabuchneri (formerly Lactobacillus parabuchneri) LMG 11457T, and Fl. sanfranciscensis JCM 5669T. Branches in this clade, representing species differentiation, were supported by high bootstrap values, ranging between 63% and 100%.
Specifically, BIA10 and BIB7 strains were grouped with Furl. rossiae DSM 15814T. The two strains shared 100% of identity and were phylogenetic divergent with the reference strain of 1.51%. The strain BFA1 was identified as Lacp. Plantarum, having 100% sequence similarity to Lacp. plantarum JCM 1149T. The branching was supported by 100% of bootstrap value, meaning an inference with high confidence. Isolate BIA2 clustered with Lenl. parabuchneri LMG 11457T, supported by a 100% bootstrap value and sharing 100% identity. The remaining strains, represented by BMB5, BMA2, BMB7, BMA8, and BMA10, were included in the Fl. sanfranciscensis clade (bootstrap value 100%). All references and isolates showed low divergence, since the maximum phylogenetic distances observed was 0.2% in the case of BMB10 with BMA8. Given the high sequence similarities and the high confidence represented by the bootstrap values, it is reasonable to assume that all of the isolates belonging to Fl. sanfranciscensis clade could be assigned to this species. Our results are consistent with previous evidence that reported the remarkable polymorphism of the 16S rDNA within the Fl. sanfranciscensis species and placed this microorganism in the Lacp. plantarum group based on the 16S rDNA phylogenetic analysis [27].

3.4. Volatile Organic Compounds Profiling

The VOCs profiling of sourdoughs and following dough samples revealed significant quantitative and qualitative differences, which can be attributed to the different microflora of each of the two systems A and B, given that the ingredients used were the same.
Sourdoughs were characterized by much higher total VOCs, as can be inferred from Figure 6. Indeed, the volatiles produced in the sourdough undergo a dramatic dilution from the addition of ingredients. It is noteworthy that the significantly higher (1.5-fold) VOCs levels, expressed as total chromatographic area, observed in MB (4.76 × 107) compared with MA (3.16 × 107), are followed by significantly higher levels also recovered in IMPB (2.26 × 107) and FINB (1.88x107) samples compared with IMPA (1.62 × 107) and FINA (1.56 × 107), respectively, even if the difference is slightly reduced to 1.2-fold.
The relative abundance of chemical classes of VOCs in the volatile fraction of samples is shown in Figure 7. As much as >80% of the volatile fraction of sourdoughs was made up of alcohols and esters together, while aldehydes and acids each represented about 6-7% of the total chromatographic area.
Interestingly, the VOCs profile of MA was characterized by a balanced level of esters and alcohols, while in MB samples the alcohols prevailed; this finding agrees with the presence of S. cerevisiae in B samples. The relative composition of the volatile fraction changed after formulation and leavening, the most noteworthy changes being the increase in aldehydes from 6–7% to 17–23% of the total VOCs, a reduction in esters to 5–7%, and the development of ketones, represented by acetoin (3-hydroxy, 2-butanone), which was not detected in the sourdough samples.
The main aliphatic aldehydes found were acetaldehyde, characterizing the VOC profile of all samples and showing an evident increasing trend from MB to FINB samples (Table 5). To a lower extent, there were 2-butenal, only found in fermented dough samples, hexanal, the product of lipid oxidation by the lipoxygenase pathway [44], and heptenal, retrieved only in MA and MB samples; they accounted for about 4% and 0.3–0.7%, respectively. These aldehydes could not be detected in further steps, according to previous reports [16]; however, the hexanal loss was not followed by a corresponding 1-hexanol increase during maturation.
Among aromatic aldehydes, phenylacetaldehyde accounted for about 5.5–7.6% of total VOCs in dough samples, while benzaldehyde accounted for 0.5–2%. The B samples always showed higher relative abundances of these compounds, while none of them were detectable in sourdough samples, suggesting the contribution of dough ingredients in their formation.
Among esters, ethyl acetate represented about 40% and 21% of the total VOCs of sourdoughs MA and MB, respectively, followed by isoamyl acetate (about 2%), ethyl lactate (1.5–2%), and ethyl hexanoate, while their concentration was dramatically reduced in the leavened products. This finding agrees with previous data on the natural fermentation of bakery products [16], which report ethyl acetate relative abundances ranging from 16 to 27% in sourdough, falling to about 3% in leavened products.
Ethyl butanoate, which was not detected in sourdoughs, was determined in both leavened products, at about 0.8-0.9% of total VOCs.
VOCs belonging to the alcohols group were dominant, with ethanol being the most concentrated, followed by isoamyl alcohol and phenethyl alcohol. Isoamyl alcohol, ranging from 6.4 to 8.9% of total VOCs in the dough samples, was also present at comparable levels in the sourdough samples and always at higher relative abundance in B samples: this alcohol is considered the most impacting aroma component produced by yeasts, together with isobutanol [45]. Phenethyl alcohol represented only 0.5% of VOCs in MA, while it was more concentrated in MB (3.3%); the relative abundance of this alcohol increased with the production steps, representing 3.9 and 8.2% of total VOCs in FINA and FINB, respectively. As for ethanol, phenethyl alcohol was always more concentrated in B samples. Isobutanol was not detected in sourdough samples but represented 2.5% of total VOCs in FINB, while it was not retrieved in FINA.
Interestingly, MA and MB samples did not show ketones among VOCs, while IMPA and IMPB, as well as FINA and FINB, revealed the presence of acetoin (3-hydroxy-2-butanone) up to about 8% of total VOCs: this compound can be attributed to both fermentation by S. cerevisiae and homofermentative LAB [46] and to the contribution of butter in the dough formulation. According to Montanari et al. [16], acetoin was not retrieved in sourdough, but could only be detected in the leavened dough at higher levels in dough samples A.

3.5. Correlation of VOCs and Species Detected in Sourdough and Dough Samples

To understand the correlation among the detected VOCs and bacteria or yeasts isolated a Pearson’s correlation index was calculated. The correlation plot is reported in Figure 8. Considering bacteria isolates, the major contributor to the aroma composition was attributed to Fl. sanfranciscensis, which was shown to be positively correlated (R2 ≥ 0.6; p < 0.05) with a wide range of VOCs belonging to four classes, such as aldehydes, esters, alcohols, and acids, as well as with a minor compound as limonene. Two other species, specifically Leuc. citreum and Lacp. plantarum, contributed in a moderate way to the aroma production, resulting in each of them being positively correlated to four different compounds. In particular, Leuc. citreum was correlated to acetaldehyde, ethyl butanoate, ethyl hexanoate, and acetoin (R2 ≥ 0.6; p < 0.05). In the case of Lacp. plantarum, the major contribution was related to 3,3-dimethyl hexane, ethyl butanoate, 2-butenal, and 1-butyl-2-propyl cyclopentane. Interestingly, Furl. rossiae was shown to have a good correlation to ethyl alcohol production (R2 = 0.83; p < 0.05). The other two species detected, Lenl. parabuchneri and Leuc. mesenteroides, were significantly correlated with just one compound, precisely 2-hexanol (R2 = 0.87; p < 0.05) and isobutanol (R 2 = 0.71; p < 0.05), respectively. Regarding the yeasts detected, K. humilis and S. cerevisiae, their contribution in terms of VOCs was shown to be, generally, very low. In particular, S. cerevisiae was found to be positively correlated with isoamyl alcohol, ethyl hexanoate, phenethyl alcohol. No relevant correlation was highlighted for K. humilis.

4. Discussion

In the present work, the viable microbial population of different dough samples collected during the production process of Panettone was assessed.
The sourdough MA showed a strong dominance of K. humilis species, with only one sample identified as a group-species K. humilis/K. parahumilis. This homogeneous result was already observed in traditional sourdoughs [8,10,13,47]. In sourdough MB, four sequenced strains belonged to the S. cerevisiae species, while the other strains belonged to the K. humilis species. The interactions of these two yeast species were also usually found in several traditional sourdoughs [8,48,49]. In our case, the minor number of refreshments of MB sample and the reduced fermentation step probably could have favored the dominance of S. cerevisiae, thanks to its rapid metabolism and growth rate.
Concerning LAB, the results showed the prevalence of heterofermentative-obliged bacteria, with the exception of Lacp. plantarum, which is a facultative heterofermentative bacterium, normally colonizing sourdoughs characterized by continuous back-slopping and low incubation temperature [3]. Fl. sanfranciscensis was found in both sourdoughs with a major occurrence in the A samples. Strains of this species are commonly reported as the main components of sourdoughs due to their ability to adapt to the highly selective conditions of this environment [4,6,18,49]. Moreover, Fl. sanfranciscensis strains utilize preferentially maltose and are generally unable to ferment fructose; therefore, there is a non-competitive association of this species with other different LABs. For instance, a stable association is reported between Fl. sanfranciscensis and Lacp. Plantarum, as the latter species preferentially ferments glucose and fructose rather than maltose, whose metabolism is subject to carbon catabolite repression [3,50,51]. Regarding maltose utilization, another stable association is described between Fl. sanfranciscensis and maltose-negative yeasts such as K. humilis. On the other hand, the interactions of Fl. sanfranciscensis and S. cerevisiae are reported to be strain-specific and evidently competitive [52].
Regarding Furl. Rossiae, its presence in both sourdough samples is quite interesting. This species was described by Corsetti et al. [53] and subsequently studied for its fermentative properties by Di Cagno et al. [54]. Furl. Rossiae was found in both samples A and B, mostly in association with other LAB species, as shown in IMPA. The species Leuc. Citreum and Leuc. Mesenteroides, which were also found in both dough types, are typical of the last stages of the Panettone process [18]. Their presence, along with Lacp. Plantarum, Furl. Rossiae, and Fl. Sanfranciscensis, could be explained by the low temperature of the process [3]. The occurrence of Lacp. plantarum in sourdough is also widely described as a both sub-dominant and dominant species [18,55]. Less common but explainable is the presence of Lenl. parabuchneri, a species commonly found in beer and dairy products [55].
Concerning the VOCs detected and their correlation to the yeasts and LAB species found in the samples, the contribution of the microbial activity to the production of the different metabolites was evident.
Noteworthy are the VOCs belonging to the alcohols group. Ethanol, isoamyl alcohol, and phenethyl alcohol were much higher in MB compared to MA and, even though their amount in IMP and FIN samples decreased because of a dilution effect (due to addition of ingredients), they kept higher concentrations in the different stages of method B, with a final difference by around 37% between the two methods. This result may be ascribed to the presence of S. cerevisiae in the B samples. Hansen and Schieberle [56] reported that the content of some volatile compounds in baked leavened products such as bread is related to the concentration in the corresponding sourdough, where they are usually more concentrated, since dough recipes imply the use of a fraction of sourdough added to other ingredients. Ethanol represented about 31% and 44% in MA and MB, respectively, thus confirming itself as the principal alcohol in the volatile fraction. Similarly, Montanari et al. [16] reported a relative composition ranging from 37 to 42% in the volatile fraction of sourdough samples, which were of the same origin but obtained with different maturation conditions.
In the present study, the different origins and compositions of the two sourdoughs gave rise to higher differences and to a wider range of variability. In particular, the comparison between the two methods for sourdough preparation highlighted higher total esters content in MA samples; however, isoamyl acetate and ethyl lactate were significantly higher in MB samples (Table 5). Ethyl acetate was detected, though to a lower extent, in the leavened products, decreasing from IMP to FIN samples, while ethyl lactate and isoamyl acetate were not detected; this trend also agreed with the literature [16]. The decrease in ethyl acetate through the production steps was limited in method B, which revealed the occurrence of S. cerevisiae.
According to De Luca et al. [57], alcohols, acids, and aldehydes represent the most characteristic VOC families influencing the fruity, green, floral, sweet, alcoholic, and fatty acid odors in bread. The LAB species detected in our samples were differently correlated to the specific VOCs, showing a complementary behavior in some cases. This supports the use of a defined combination of selected strains to achieve the desired aroma production.

5. Conclusions

In the present work, the molecular characterization of the dough’s population revealed the presence of the dominant yeasts, K. humilis and S. cerevisiae, as well as the LAB species, namely Fl. sanfranciscensis, Lacp. plantarum, Furl. rossiae, Lenl. parabuchneri, Leuc. mesenteroides, and Leuc. citreum.
Results highlight the relevance of refreshment procedures as the main factor in maintaining a stable and active microflora in industrial conditions.
Moreover, the correlation among the selected strains and the VOCs produced in the different samples allowed us to estimate the species contribution to aroma formation.
Therefore, the research of sourdoughs’ selected strains, with known metabolic properties and high technological performances, is fundamental for their exploitation as starter cultures capable of controlling the manufacturing of sourdoughs. A pool of LAB and yeasts of this study deposited at the UMCC culture collection comprises candidate strains for the development of single and multiple selected starters for Panettone production.

Author Contributions

Conceptualization and original draft preparation: L.D.V.; investigation and visualization: L.D.V., G.I., S.L.C., F.L.; writing—review and editing: L.D.V., S.L.C., M.G., F.L., A.P.; supervision, project administration and funding acquisition: A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

We thank Panificio Fantuzzi s.r.l. (Reggio Emilia, Italy) for providing us the two sourdoughs samples in the different stages of the process. A special thanks to Katia D’Ambra and Andrea Feroce for their contribution in laboratory experiments. The JRU MIRRI-IT (http://www.mirri-it.it/) and IS_MIRRI21 (https://ismirri21.mirri.org/) are greatly acknowledged for scientific support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vernocchi, P.; Valmorri, S.; Gatto, V.; Torriani, S.; Gianotti, A.; Suzzi, G.; Guerzoni, M.E.; Gardini, F. A Survey on yeast microbiota associated with an Italian traditional sweet-leavened baked good fermentation. Food Res. Int. 2004, 37, 469–476. [Google Scholar] [CrossRef]
  2. Lattanzi, A.; Minervini, F.; Di Cagno, R.; Diviccaro, A.; Antonielli, L.; Cardinali, G.; Cappelle, S.; De Angelis, M.; Gobbetti, M. The lactic acid bacteria and yeast microbiota of eighteen sourdoughs used for the manufacture of traditional Italian sweet leavened baked goods. Int. J. Food Microbiol. 2013, 163, 71–79. [Google Scholar] [CrossRef]
  3. De Vuyst, L.; Vrancken, G.; Ravyts, F.; Rimaux, T.; Weckx, S. Biodiversity, ecological determinants, and metabolic exploitation of sourdough microbiota. Food Microbiol. 2009, 26, 666–675. [Google Scholar] [CrossRef] [PubMed]
  4. De Vuyst, L.; Neysens, P. The sourdough microflora: Biodiversity and metabolic interactions. Trends Food Sci. Technol. 2005, 16, 43–56. [Google Scholar] [CrossRef]
  5. Foschino, R.; Gallina, S.; Andrighetto, C.; Rossetti, L.; Galli, A. Comparison of Cultural Methods for the Identification and Molecular Investigation of Yeasts from Sourdoughs for Italian Sweet Baked Products. FEMS Yeast Res. 2004, 4, 609–618. [Google Scholar] [CrossRef] [Green Version]
  6. Garofalo, C.; Silvestri, G.; Aquilanti, L.; Clementi, F. PCR-DGGE Analysis of lactic acid bacteria and yeast dynamics during the production processes of three varieties of Panettone. J. Appl. Microbiol. 2008, 105, 243–254. [Google Scholar] [CrossRef] [PubMed]
  7. Rogalski, E.; Ehrmann, M.A.; Vogel, R.F. Role of Kazachstania humilis and Saccharomyces cerevisiae in the strain-specific assertiveness of Fructilactobacillus sanfranciscensis strains in rye sourdough. Eur. Food Res. Technol. 2020, 246, 1817–1827. [Google Scholar] [CrossRef]
  8. Venturi, M.; Guerrini, S.; Vincenzini, M. Stable and non-competitive association of Saccharomyces cerevisiae, Candida milleri and Lactobacillus sanfranciscensis during manufacture of two traditional sourdough baked goods. Food Microbiol. 2012, 31, 107–115. [Google Scholar] [CrossRef] [PubMed]
  9. Di Cagno, R.; Pontonio, E.; Buchin, S.; De Angelis, M.; Lattanzi, A.; Valerio, F.; Gobbetti, M.; Calasso, M. Diversity of the lactic acid bacterium and yeast microbiota in the switch from firm- to liquid-sourdough fermentation. Appl. Environ. Microbiol. 2014, 80, 3161–3172. [Google Scholar] [CrossRef] [Green Version]
  10. De Vuyst, L.; Van Kerrebroeck, S.; Harth, H.; Huys, G.; Daniel, H.M.; Weckx, S. Microbial ecology of sourdough fermentations: Diverse or uniform? Food Microbiol. 2014, 37, 11–29. [Google Scholar] [CrossRef]
  11. Zheng, J.; Wittouck, S.; Salvetti, E.; Franz, C.M.A.P.; Harris, H.M.B.; Mattarelli, P.; O’toole, P.W.; Pot, B.; Vandamme, P.; Walter, J.; et al. A Taxonomic Note on the Genus Lactobacillus: Description of 23 Novel Genera, Emended Description of the Genus Lactobacillus Beijerinck 1901, and Union of Lactobacillaceae and Leuconostocaceae. Int. J. Syst. Evol. Microbiol. 2020, 70, 2782–2858. [Google Scholar] [CrossRef]
  12. Gobbetti, M.; Corsetti, A.; Rossi, J.; La Rosa, F.; De Vincenzi, S. Identification and clustering of lactic acid bacteria and yeasts from wheat sourdoughs of central Italy. Ital. J. Food. Sci. 1994, 1, 85–94. [Google Scholar]
  13. Gullo, M.; Romano, A.D.; Pulvirenti, A.; Giudici, P. Candida humilis—Dominant species in sourdoughs for the production of durum wheat bran flour bread. Int. J. Food Microbiol. 2003, 80, 55–59. [Google Scholar] [CrossRef]
  14. Pulvirenti, A.; Solieri, L.; Gullo, M.; De Vero, L.; Giudici, P. Occurrence and dominance of yeast species in sourdough. Lett. Appl. Microbiol. 2004, 38, 113–117. [Google Scholar] [CrossRef] [Green Version]
  15. Chavan, R.S.; Chavan, S.R. Sourdough Technology-A Traditional Way for Wholesome Foods: A Review. Compr. Rev. Food Sci. Food Saf. 2011, 10, 169–182. [Google Scholar] [CrossRef]
  16. Montanari, C.; Bargossi, E.; Lanciotti, R.; Chinnici, F.; Gardini, F.; Tabanelli, G. Effects of two different sourdoughs on the characteristics of Pandoro, a typical Italian sweet leavened baked good. LWT 2014, 59, 289–299. [Google Scholar] [CrossRef]
  17. Gobbetti, M.; Rizzello, C.G.; Di Cagno, R.; De Angelis, M. How the sourdough may affect the functional features of leavened baked goods. Food Microbiol. 2014, 37, 30–40. [Google Scholar] [CrossRef] [PubMed]
  18. Reale, A.; Di Renzo, T.; Boscaino, F.; Nazzaro, F.; Fratianni, F.; Aponte, M. Lactic acid bacteria biota and aroma profile of italian traditional sourdoughs from the irpinian area in Italy. Front. Microbiol. 2019, 10, 1621. [Google Scholar] [CrossRef] [Green Version]
  19. De Vero, L.; Iosca, G.; Gullo, M.; Pulvirenti, A. Functional and Healthy Features of Conventional and Non-Conventional Sourdoughs. Appl. Sci. 2021, 11, 3694. [Google Scholar] [CrossRef]
  20. Pulvirenti, A.; Rainieri, S.; Boveri, S.; Giudici, P. Optimizing the selection process of yeast starter cultures by preselecting strains dominating spontaneous fermentations. Can. J. Microbiol. 2009, 55, 326–332. [Google Scholar] [CrossRef] [PubMed]
  21. Reale, A.; Zotta, T.; Ianniello, R.G.; Mamone, G.; Di Renzo, T. Selection criteria of lactic acid bacteria to be used as starter for sweet and salty leavened baked products. LTW 2020, 133, 110092. [Google Scholar] [CrossRef]
  22. Palla, M.; Blandino, M.; Grassi, A.; Giordano, D.; Sgerri, C.; Quartacci, M.F.; Reyneri, A.; Agnolucci, M.; Giovanetti, M. Characterization and selection of functional yeast strains during sourdough fermentation of different cereal wholegrain flours. Sci. Rep. 2020, 10, 12856. [Google Scholar] [CrossRef]
  23. Ercolini, D.; Pontonio, E.; De Filippis, F.; Minervini, F.; La Storia, A.; Gobbetti, M.; Di Cagno, R. Microbial ecology dynamics during rye and wheat sourdough preparation. Appl. Environ. Microbiol. 2013, 79, 7827–7836. [Google Scholar] [CrossRef] [Green Version]
  24. De Vero, L.; Boniotti, M.B.; Budroni, M.; Buzzini, P.; Cassanelli, S.; Comunian, R.; Gullo, M.; Logrieco, A.F.; Mannazzu, I.; Musumeci, R.; et al. Preservation, characterization and exploitation of microbial biodiversity: The perspective of the italian network of culture collections. Microorganisms 2019, 7, 685. [Google Scholar] [CrossRef] [Green Version]
  25. Hoffman, C.S.; Winston, F. A ten-minute DNA preparation from yeast efficiently releases autonomous plasmids for transformation of Escherichia coli. Gene 1987, 57, 267–272. [Google Scholar] [CrossRef]
  26. Ausubel, F.M.; Brent, R.; Kingston, R.E.; Moore, D.D.; Seidman, J.G.; Smith, J.A.; Struhl, K. Current Protocols in Molecular Biology: Preface, Ringbou ed.; John Wiley & Sons Inc: Hoboken, NJ, USA, 2003; ISBN 047150338X. [Google Scholar]
  27. Valcheva, R.; Kabadjova, P.; Rachman, C.; Ivanova, I.; Onno, B.; Prévost, H.; Dousset, X. A Rapid PCR procedure for the specific identification of Lactobacillus sanfranciscensis, based on the 16S–23S Intergenic Spacer Regions. J. Appl. Microbiol. 2007, 102, 290–302. [Google Scholar] [CrossRef]
  28. Versalovic, J.; Schneider, M.; De Bruijn, F.J.; Lupski, J.R. Genomic fingerprinting of bacteria using repetitive sequence-based polymerase chain reaction. Methods Mol. Cell. Biol. 1994, 5, 25–40. [Google Scholar]
  29. La China, S.; De Vero, L.; Anguluri, K.; Brugnoli, M.; Mamlouk, D.; Gullo, M. Kombucha tea as a reservoir of cellulose producing bacteria: Assessing diversity among Komagataeibacter isolates. Appl. Sci. 2021, 11, 1595. [Google Scholar] [CrossRef]
  30. Kurtzman, C.P. Identification of food and beverage spoilage yeasts from DNA sequence analyses. Int. J. Food Microbiol. 2015, 213, 71–78. [Google Scholar] [CrossRef]
  31. Weiss, S.; Samson, F.; Navarro, D.; Casaregola, S. YeastIP: A database for identification and phylogeny of Saccharomycotina yeasts. FEMS Yeast Res. 2013, 13, 117–125. [Google Scholar] [CrossRef] [Green Version]
  32. Sato, H.; Yanagida, F.; Shinohara, T.; Yokotsuka, K. Restriction Fragment Length Polymorphism Analysis of 16S RRNA genes in lactic acid aacteria isolated from Red Wine. J. Biosci. Bioeng. 2000, 90, 335–337. [Google Scholar] [CrossRef]
  33. Higgins, D.; Thompson, J.; Gibson TThompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar]
  34. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar]
  35. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v4: Recent updates and new developments. Nucleic Acids Res. 2019, 47, 256–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Taiyun, W.; Viliam, S. R Package “corrplot”: Visualization of a Correlation Matrix: Version 0.84. 2017. Available online: https://github.com/taiyun/corrplot (accessed on 2 April 2021).
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 2 April 2021).
  38. Esteve-Zarzoso, B.; Belloch, C.; Uruburu, F.; Querol, A. Identification of Yeasts by RFLP Analysis of the 5.8S rRNA Gene and the Two Ribosomal Internal Transcribed Spacers. Int. J. Syst. Bacteriol. 1999, 49, 329–337. [Google Scholar] [CrossRef] [Green Version]
  39. Kurtzman, C.P.; Robnett, C.J. Identification and Phylogeny of Ascomycetous Yeasts from Analysis of Nuclear Large Subu.Pdf. Antonie Van Leeuwenhoek 1998, 98, 331–371. [Google Scholar] [CrossRef] [PubMed]
  40. Lachance, M.A.; Wijayanayaka, T.M.; Bundus, J.D.; Wijayanayaka, D.N. Ribosomal DNA sequence polymorphism and the delineation of two ascosporic yeast species, Metschnikowia agaves and Starmerella bombicola. FEMS Yeast Res. 2011, 11, 324–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Jacques, N.; Sarilar, V.; Urien, C.; Lopes, M.R.; Morais, C.G.; Uetanabaro, A.P.T.; Tinsley, C.R.; Rosa, C.A.; Sicard, D.; Casaregola, S. Three novel ascomycetous yeast species of the Kazachstania clade, Kazachstania saulgeensis sp. nov., Kazachstania serrabonitensis sp. nov. and Kazachstania australis sp. nov. Reassignment of Candida humilis to Kazachstania humilis f.a. comb. nov. and Candida pseudohumilis to Kazachstania pseudohumilis f.a. comb. nov. Int. J. Syst. Evol. Microbiol. 2016, 66, 5192–5200. [Google Scholar]
  42. Coppola, S.; Blaiotta, G.; Ercolini, D.; Moschetti, G. Molecular evaluation of microbial diversity occurring in different types of Mozzarella cheese. J. Appl. Microbiol. 2001, 90, 414–420. [Google Scholar] [CrossRef] [Green Version]
  43. Rachman, C.N.; Kabadjova, P.; Prevost, H.; Dousset, X. Identification of Lactobacillus alimentarius and Lactobacillus farciminis with 16S-23S rDNA intergenic spacer region polymorphism and PCR amplification using species-specific oligonucleotide. J. Appl. Microbiol. 2003, 95, 1207–1216. [Google Scholar] [CrossRef] [Green Version]
  44. Martínez-Anaya, M.A. Enzymes and bread flavor. J. Agr. Food Chem. 1996, 44, 2469–2480. [Google Scholar] [CrossRef]
  45. Lund, B.; Hansen, A.; Lewis, M.J. The Influence of Yeast and Sourdough Fermentation on Volatile Compounds in Rye Bread Crumb. In Cereal Science and Technology, Proceedings from 23. Nordic Cereal Congress; Munck, L., Ed.; The Danish Cereal Society: Copenhagen, Denmark, 1987; pp. 185–195. [Google Scholar]
  46. Damiani, P.; Gobbetti, M.; Cossignani, L.; Corsetti, A.; Simonetti, M.S.; Rossi, J. The sourdough microflora. Characterization of hetero- and homo- fermentative lactic acid bacteria, yeasts and their interactions on the basis of the volatile compounds produced. LWT 1996, 29, 63–70. [Google Scholar] [CrossRef]
  47. De Vuyst, L.; Harth, H.; Van Kerrebroeck, S.; Leroy, F. Yeast Diversity of Sourdoughs and Associated Metabolic Properties and Functionalities. Int. J. Food Microbiol. 2016, 239, 26–34. [Google Scholar] [CrossRef]
  48. Carbonetto, B.; Nidelet, T.; Guezenec, S.; Perez, M.; Segond, D.; Sicard, D. Interactions between Kazachstania humilis yeast species and Lactic Acid Bacteria in Sourdough. Microorganisms 2020, 8, 240. [Google Scholar] [CrossRef] [Green Version]
  49. Fujimoto, A.; Ito, K.; Narushima, N.; Miyamoto, T. Identification of Lactic Acid Bacteria and Yeasts, and characterization of food components of Sourdoughs used in Japanese bakeries. J. Biosci. Bioeng. 2019, 127, 575–581. [Google Scholar] [CrossRef]
  50. Gobbetti, M. The sourdough microflora: Interactions of lactic acid bacteria and yeasts. Trends Food Sci. Technol. 1998, 9, 267–274. [Google Scholar] [CrossRef]
  51. Corsetti, A.; Lavermicocca, P.; Morea, M.; Baruzzi, F.; Tosti, N.; Gobbetti, M. Phenotypic and molecular identification and clustering of lactic acid bacteria and yeasts from wheat (species Triticum durum and Triticum aestivum) sourdoughs of Southern Italy. Int. J. Food Microbiol. 2001, 64, 95–104. [Google Scholar] [CrossRef]
  52. Rogalski, E.; Ehrmann, M.A.; Vogel, R.F. Strain-specific interaction of Fructilactobacillus sanfranciscensis with yeasts in the sourdough fermentation. Eur. Food Res. Technol. 2021, 247, 1437–1447. [Google Scholar] [CrossRef]
  53. Corsetti, A.; Settanni, L.; Van Sinderen, D.; Felis, G.E.; Dellaglio, F.; Gobbetti, M. Lactobacillus rossii sp. nov., isolated from Wheat Sourdough. Int. J. Syst. Evol. Microbiol. 2005, 55, 35–40. [Google Scholar] [CrossRef] [Green Version]
  54. Di Cagno, R.; De Angelis, M.; Gallo, G.; Settanni, L.; Berloco, M.G.; Siragusa, S.; Parente, E. Genotypic and phenotypic diversity of Lactobacillus rossiae strains isolated from Sourdough. J. Appl. Microbiol. 2007, 103, 821–835. [Google Scholar] [CrossRef]
  55. Arena, M.P.; Russo, P.; Spano, G.; Capozzi, V. From Microbial ecology to innovative applications in food quality improvements: The Case of Sourdough as a Model Matrix. J. Multidiscip. Sci. J. 2020, 3, 9–19. [Google Scholar] [CrossRef] [Green Version]
  56. Hansen, A.; Schieberle, P. Generation of aroma compounds during sourdough fermentation: Applied and fundamental aspects. Trends Food Sci. Technol. 2005, 16, 85–94. [Google Scholar] [CrossRef]
  57. De Luca, L.; Aiello, A.; Pizzolongo, F.; Blaiotta, G.; Aponte, M.; Romano, R. Volatile Organic Compounds in Breads Prepared with Different Sourdoughs. Appl. Sci. 2021, 11, 1330. [Google Scholar] [CrossRef]
Microorganisms 09 01093 g001 550
Figure 1. Artisanal Panettone production using two different sourdoughs (MA and MB). The samples highlighted were collected concurrently at the same stages for their investigation.
Figure 1. Artisanal Panettone production using two different sourdoughs (MA and MB). The samples highlighted were collected concurrently at the same stages for their investigation.
Microorganisms 09 01093 g001
Microorganisms 09 01093 g002 550
Figure 2. Dendrogram obtained from unweighted pair group method with arithmetic mean (UPGMA) analysis, using Pearson’s coefficient, of the (GTG)5 rep-PCR patterns. The similarity threshold for biotypes discrimination was 89% (blue line). The cophenetic coefficient is represented by numbers and dots colored as red, orange, yellow, or green, based on the branch quality. Two major clades are distinct with a threshold of 40% (red line).
Figure 2. Dendrogram obtained from unweighted pair group method with arithmetic mean (UPGMA) analysis, using Pearson’s coefficient, of the (GTG)5 rep-PCR patterns. The similarity threshold for biotypes discrimination was 89% (blue line). The cophenetic coefficient is represented by numbers and dots colored as red, orange, yellow, or green, based on the branch quality. Two major clades are distinct with a threshold of 40% (red line).
Microorganisms 09 01093 g002
Microorganisms 09 01093 g003 550
Figure 3. Neighbor-joining phylogenetic tree representing the phylogenetic distances of isolates from sourdough MA and MB and the strain types of identified yeasts species. The tree was made using Tamura-Nei model, and a discrete Gamma distribution was used to model evolutionary rate differences among sites. Node numbers indicate bootstrap values obtained using 1000 replicates.
Figure 3. Neighbor-joining phylogenetic tree representing the phylogenetic distances of isolates from sourdough MA and MB and the strain types of identified yeasts species. The tree was made using Tamura-Nei model, and a discrete Gamma distribution was used to model evolutionary rate differences among sites. Node numbers indicate bootstrap values obtained using 1000 replicates.
Microorganisms 09 01093 g003
Microorganisms 09 01093 g004 550
Figure 4. Dendrograms obtained from unweighted pair group method with arithmetic mean (UPGMA) analysis, using Pearson’s coefficient, of (GTG)5 rep-PCR. The similarity threshold for biotypes discrimination was 95% (blue line). The cophenetic coefficient is represented by numbers and dots colored as red, orange, yellow, or green, based on the branch quality.
Figure 4. Dendrograms obtained from unweighted pair group method with arithmetic mean (UPGMA) analysis, using Pearson’s coefficient, of (GTG)5 rep-PCR. The similarity threshold for biotypes discrimination was 95% (blue line). The cophenetic coefficient is represented by numbers and dots colored as red, orange, yellow, or green, based on the branch quality.
Microorganisms 09 01093 g004
Microorganisms 09 01093 g005 550
Figure 5. Neighbor-joining phylogenetic tree obtained by 16S multiple alignment. A gamma distribution was applied to model the rate variation among sites. The node numbers indicate the percentage of replicates in which taxa clustered at the same distance (bootstrap test). References species used in the dataset are reported, in the figure, using their former name and the associated NCBI id. Given the recent re-organization of the Lactobacillus genus, the different genera are highlighted with different colors and named using the new taxonomical nomenclature, according to Zheng et al. [11].
Figure 5. Neighbor-joining phylogenetic tree obtained by 16S multiple alignment. A gamma distribution was applied to model the rate variation among sites. The node numbers indicate the percentage of replicates in which taxa clustered at the same distance (bootstrap test). References species used in the dataset are reported, in the figure, using their former name and the associated NCBI id. Given the recent re-organization of the Lactobacillus genus, the different genera are highlighted with different colors and named using the new taxonomical nomenclature, according to Zheng et al. [11].
Microorganisms 09 01093 g005
Microorganisms 09 01093 g006 550
Figure 6. Mean total chromatographic area of the identified volatile organic compounds (VOCs) in sourdoughs and dough samples A and B.
Figure 6. Mean total chromatographic area of the identified volatile organic compounds (VOCs) in sourdoughs and dough samples A and B.
Microorganisms 09 01093 g006
Microorganisms 09 01093 g007 550
Figure 7. Relative abundance (%) of the chemical classes of volatile organic compounds (VOCs), determined in sourdough and dough samples at different steps of Panettone preparation.
Figure 7. Relative abundance (%) of the chemical classes of volatile organic compounds (VOCs), determined in sourdough and dough samples at different steps of Panettone preparation.
Microorganisms 09 01093 g007
Microorganisms 09 01093 g008 550
Figure 8. Correlation plot among discovered VOCs and isolated species, obtained using R package Corrplot and Pearson’s correlation index. The positive and negative correlation were considered based on R2 ≥ 0.6 and R2 ≤ −0.6, respectively.
Figure 8. Correlation plot among discovered VOCs and isolated species, obtained using R package Corrplot and Pearson’s correlation index. The positive and negative correlation were considered based on R2 ≥ 0.6 and R2 ≤ −0.6, respectively.
Microorganisms 09 01093 g008
Table
Table 1. Physicochemical parameters detected on the samples collected during the Panettone production.
Table 1. Physicochemical parameters detected on the samples collected during the Panettone production.
SamplespHWater Activity (aw)Total Titratable Acidity (mL NaOH 0.1N/10 g)
MA3.90 ± 0.03 d0.0.9855 ± 0.002 d10.85 ± 0.45 b
IMPA5.44 ± 0.05 b0.0.9636 ± 0.001 b3.35 ± 0.25 a
FINA5.35 ± 0.02 bc0.0.9441 ± 0.002 a3.55 ± 0.45 a
MB3.81 ± 0.02 d0.0.9894 ± 0.001 d10.45 ± 0.45 b
IMPB5.30 ± 0.07 ac0.0.9562 ± 0.002 c3.45 ± 0.30 a
FINB5.21 ± 0.03 a0.0.9411 ± 0.003 a4.05 ± 0.05 b
Values are means ± standard deviation of three replicates (n = 3). Values within a column with different superscript letters are significantly different (p < 0.05).
Table
Table 2. Colony counts on the growth media used for the isolation of yeasts and lactic acid bacteria from dough samples.
Table 2. Colony counts on the growth media used for the isolation of yeasts and lactic acid bacteria from dough samples.
SamplesMeans Counts * on the Isolation Media
YPDAMRSSDB
MA6.81 ± 0.25 b6.15 ± 0.14 ab6.17 ± 0.21 b
IMPA7.81 ± 0.05 a5.39 ± 0.12 a5.65 ± 0.21 ab
FINA7.83 ± 0.17 a5.46 ± 0.07 a5.45 ± 0.19 ab
MB7.56 ± 0.03 a6.87 ± 0.42 b6.16 ± 0.15 b
IMPB8.02 ± 0.01 a5.57 ± 0.06 a5.27 ± 0.09 a
FINB7.89 ± 0.17 a5.64 ± 0.04 a5.64 ± 0.02 ab
* expressed as log of colony forming unit (cfu)/ g of sample ± standard deviation. Values are means ± standard deviation of three replicates (n = 3). Values within a column with different superscript letters are significantly different (p < 0.05). YPDA medium (yeast extract, peptone, dextrose, agar); MRS medium (de Man Rogosa Sharpe); SDB medium (Sourdough Bacteria agar).
Table
Table 3. RFLP-PCR of ITS-5.8S rDNA.
Table 3. RFLP-PCR of ITS-5.8S rDNA.
SampleStrains *Size of Amplicon (bp)Amplicons Digestion with HaeIII (Size Fragments, bp)
MALMA2, LMA3 (UMCC 2973), LMA1, LMA4, LMA5 (UMCC 2974), LMA6, LMA7, LMA8, LMA9 (UMCC 2975), LMA10630400–230
IMPALIA1 (UMCC 2976), LIA9, LIA10, LIA2, LIA3, LIA4, LIA5 (UMCC 2977), LIA7, LIA6 (UMCC 2978), LIA8630400–230
FINALFA1 (UMCC 2979), LFA2, LFA3 (UMCC 2980), LFA4, LFA5, LFA6 (UMCC 2981), LFA7, LFA8, LFA9, LFA10630400–230
MBLMB1 (UMCC 2982), LMB5630400–230
LMB2, LMB3, LMB7, LMB9 (UMCC 2984), LMB10880320–230–180–150
IMPBLIB1, LIB4, LIB7 (UMCC 2986), LIB10880320–230–180–150
LIB2, LIB3, LIB5 (UMCC 2985), LIB6, LIB8, LIB9880400–230
FINBLFB1 (UMCC 2987), LFB2, LFB3, LFB4, LFB6, LFB10 (UMCC 2988)880320–230–180–150
LFB5, LFB7, LFB8, LFB9630400–230
* The strains sequenced are shown in bold and in brackets are indicated the correspondent UMCC code.
Table
Table 4. PCR amplification of the 16S-23S rDNA intergenic spacer region (ISR).
Table 4. PCR amplification of the 16S-23S rDNA intergenic spacer region (ISR).
SampleStrain *16S-23S rDNA ISR Size Fragments (bp)
MABMA2 (UMCC 2989), BMA8 (UMCC 2990), BMA10 (UMCC 2991)600–700–800
IMPABIA2 (UMCC 2992), BIA3, BIA5, BIA6 (UMCC 2993), BIA7 (UMCC 2994), BIA8, BIA9, BIA10 (UMCC 2995)600–800
FINABFA1 (UMCC 2996), BFA2 (UMCC 2997)600–800
MBBMB5 (UMCC 2998), BMB7 (UMCC 2999)600–700–800
IMPBBIB7 (UMCC 3002), BIB8600–800
FINBBFB3, BFB6 (UMCC 3004), BFB8 (UMCC 3006)600–800
* The strains sequenced are shown in bold and in brackets are indicated the correspondent UMCC code.
Table
Table 5. Relative volatile organic compound (VOC) composition of sourdough and dough samples.
Table 5. Relative volatile organic compound (VOC) composition of sourdough and dough samples.
Chemical CompoundsSamples *
MAMBImpAImpBFinAFinB
Acetaldehyde2.22 ± 0.111.62 ± 0.138.43 ± 0.193.28 ± 0.018.31 ± 0.106.27 ± 0.05
2-Butenal--2.44 ± 0.161.56 ± 0.238.61 ± 0.653.80 ± 0.26
Hexanal4.13 ± 0.173.86 ± 0.37----
Heptenal0.31 ± 0.040.66 ± 0.06----
Benzaldehyde--0.55 ± 0.061.96 ± 0.120.86 ± 0.211.37 ± 0.05
Phenylacetaldehyde--6.80 ±0.477.62 ± 0.185.52 ± 0.336.08 ± 0.03
Aldehydes6.346.1418.2214.4323.3017.52
Ethyl Acetate39.73 ± 1.5120.78 ± 0.437.17 ± 0.455.43 ± 0.234.08 ± 0.045.85 ± 0.15
Ethyl butanoate----0.86 ± 0.120.76 ± 0.13
Isoamyl acetate2.03 ± 0.142.17 ± 0.23---0.39 ± 0.12
Ethyl hexanoate-0.38 ± 0.070.31 ± 0.000.20 ± 0.030.32 ± 0.020.64 ± 0.07
Hexyl acetate-0.26±0.01----
Ethyl lactate1.44 ± 0.042.14 ± 0.01----
Esters43.1925.747.485.645.267.26
Ethyl alcohol31.19 ± 0.2244.05 ± 0.0050.77 ± 0.6653.19 ± 2.5041.09 ± 2.6540.51 ± 1.46
Isobutanol--1.40 ± 0.171.75 ± 0.11-2.51 ± 0.05
Isoamyl alcohol4.30 ± 0.256.80 ± 0.326.42 ± 0.568.92 ± 1.066.70 ± 1.327.77 ± 0.25
1-Hexanol1.96 ± 0.153.48 ± 0.421.60 ± 0.230.95 ± 0.121.12 ± 0.011.25 ± 0.06
2-Hexanol--0.64 ± 0.11---
1-Octen-3-ol0.30 ± 0.030.69 ± 0.03----
Phenethyl alcohol 0.51 ± 0.013.27 ± 0.452.70 ± 0.387.44 ± 1.093.90 ± 0.628.22 ± 1.96
Alcohols38.2758.3063.5272.2652.8260.25
Acetic acid5.34 ± 1.445.93 ± 0.091.23 ± 0.051.19 ± 0.091.11 ± 0.372.76 ± 0.35
Hexanoic acid0.43 ± 0.070.94 ± 0.04--0.45 ± 0.120.31 ± 0.01
Acids5.766.871.231.191.112.76
Limonene 2.49 ± 0.070.74 ± 0.011.08 ± 0.240.75 ± 0.110.96 ± 0.200.67 ± 0.06
Acetoin --7.42 ± 0.375.05 ± 0.339.39 ± 0.267.64 ±0.91
3,3-dimethyl hexane----2.48 ± 0.361.47 ± 0.19
1-butyl-2-propyl cyclopentane----3.24 ± 0.671.00 ± 0.40
2-Pentylfuran0.92 ± 0.000.72 ± 0.02----
Benzene,1,3-bis(1,1-dimethylethyl)3.02 ± 0.281.03 ± 0.140.38 ± 0.060.69 ± 0.121.44 ± 0.021.44 ± 0.04
Others6.432.498.896.4917.5112.21
* Data are expressed as % GC area ± standard deviation. The chemical classes of VOCs and their corresponding amounts are indicated in italics.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

De Vero, L.; Iosca, G.; La China, S.; Licciardello, F.; Gullo, M.; Pulvirenti, A. Yeasts and Lactic Acid Bacteria for Panettone Production: An Assessment of Candidate Strains. Microorganisms 2021, 9, 1093. https://doi.org/10.3390/microorganisms9051093

AMA Style

De Vero L, Iosca G, La China S, Licciardello F, Gullo M, Pulvirenti A. Yeasts and Lactic Acid Bacteria for Panettone Production: An Assessment of Candidate Strains. Microorganisms. 2021; 9(5):1093. https://doi.org/10.3390/microorganisms9051093

Chicago/Turabian Style

De Vero, Luciana, Giovanna Iosca, Salvatore La China, Fabio Licciardello, Maria Gullo, and Andrea Pulvirenti. 2021. "Yeasts and Lactic Acid Bacteria for Panettone Production: An Assessment of Candidate Strains" Microorganisms 9, no. 5: 1093. https://doi.org/10.3390/microorganisms9051093

APA Style

De Vero, L., Iosca, G., La China, S., Licciardello, F., Gullo, M., & Pulvirenti, A. (2021). Yeasts and Lactic Acid Bacteria for Panettone Production: An Assessment of Candidate Strains. Microorganisms, 9(5), 1093. https://doi.org/10.3390/microorganisms9051093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop