Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2007, p. 6262-6269, Vol. 73, No. 19
0099-2240/07/$08.00+0 doi:10.1128/AEM.00894-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Laboratory of Microbiology,1 BCCM/LMG Bacteria Collection, Department of Biochemistry, Physiology, and Microbiology, Ghent University, K. L. Ledeganckstraat 35, B-9000 Ghent, Belgium,2 Research Group of Industrial Microbiology and Food Biotechnology, Department of Applied Biological Sciences and Engineering, Vrije Universiteit Brussel, Pleinlaan 2, B-1050 Brussels, Belgium3
Received 20 April 2007/ Accepted 27 July 2007
|
|
|---|
|
|
|---|
The increasing interest in starter cultures for sourdough fermentations requires better insights into the genetic and phenotypic diversity of strains for exploitation in technological processes. Although still important for the description of new species, conventional methods such as carbohydrate fermentation pattern analysis and cell wall analyses are not reliable for the accurate identification of LAB at the species level. The use of 16S rRNA gene sequencing is generally regarded as a more reliable solution for the classification and identification of LAB, although the differentiation of closely related Lactobacillus species is not always straightforward, due to a high degree of conservation within this genus (12). Alternatively, sequence-based identification using one or more protein-coding genes has been explored as an alternative to determine the genomic relatedness between LAB strains (4). In this regard, the phenylalanyl-tRNA synthase gene (pheS) has proven to be a valuable tool for the identification of Lactobacillus species and the delineation of novel taxa (24, 24a, 28, 33).
The objective of the present study was to assess the taxonomic diversity of LAB isolated from traditional sourdoughs collected from local bakeries throughout different provinces of Belgium. The assessment of the LAB biodiversity was based entirely on genotypic methods, such as repetitive DNA element PCR (rep-PCR) fingerprinting and sequencing of the pheS gene. The rep-PCR method was used to reduce the initial collection of LAB isolates to a set of genotypically unique isolates and allowed the tentative identification of the majority of sourdough LAB isolates. These identifications were further confirmed by pheS sequence analysis, and the taxonomic position of one pheS cluster that could not be allocated to any of the LAB species included in the pheS reference database was characterized by 16S rRNA gene sequence analysis and DNA-DNA hybridization experiments. Using statistical analysis, the influence of geographical origin, flour type, bakery environment, and technological characteristics on the LAB biodiversity in Belgian sourdoughs was assessed.
|
|
|---|
Sourdough sampling and isolation of LAB and yeasts.
A total of 11 different bakeries from seven provinces throughout Belgium were selected for sampling (Table 1). Samples from 21 natural sourdoughs were collected aseptically and kept at 4°C. Microbiological analyses were performed within 24 h after sampling. For microbiological analyses, sourdough samples (10 g) were homogenized in 90 ml peptone physiological solution (0.1% [wt/vol] bacteriological peptone [L37; Oxoid, Basingstoke, Hampshire, United Kingdom] and 0.85% [wt/vol] NaCl) by using a Stomacher 80 lab blender (Seward Medical, London, United Kingdom), after which serial dilutions in peptone physiological solution were prepared. LAB and yeasts were isolated on MRS-5 agar containing 0.1 g of cycloheximide (C-7698; Sigma-Aldrich Co., St. Louis, MO) and yeast-glucose (YG) agar (21a) supplemented with 0.1 g of chloramphenicol (C-0378; Sigma-Aldrich Co.) per liter, respectively. Routinely, 50 µl of the 10–4 to 10–7 dilutions was plated in triplicate on a series of MRS-5 agar plates that were incubated at four different incubation conditions: aerobic at 30°C, aerobic at 37°C, anaerobic at 30°C, and anaerobic at 37°C. Additionally, on YG plates, 50 µl of the 10–1 to 10–6 dilutions was plated on YG agar and aerobically incubated at 30°C (21a). Following incubation, an average of 40 colonies per sourdough sample was randomly selected from MRS-5 agar plates containing between 20 and 200 colonies. The colonies were further purified and stored in Microbank vials at –80°C. From YG agar plates, 10 colonies were picked and stored at 4°C on YG slants overlaid with sterile tap water. After the isolation, the remainder of the sourdough samples were stored at –20°C. A uniform code was assigned to each sourdough sample, encoding information on the depositor (D01 to D11), the type of flour used to produce the sourdough (WW, wheat flour; RR, rye flour; SS, spelt flour; WR, mixture of wheat and rye flour; and ME, mixture of wheat, rye, and spelt flour), and the time of sampling (T). The main technological characteristics of the 21 sourdough samples are shown in Table 1.
|
View this table: [in a new window] |
TABLE 1. Main characteristics of 21 Belgian traditional sourdoughs included in this studya
|
rep-PCR genomic fingerprinting of LAB isolates.
For the establishment of a taxonomical reference framework, a total of 63 type and reference strains belonging to 36 LAB species associated with sourdough were subjected to (GTG)5-PCR fingerprinting. The genomic DNA of sourdough LAB was extracted from single colonies grown on MRS-5 agar plates by alkaline lysis. For this purpose, a small portion of one colony was resuspended in 20 µl lysis buffer (2.5 ml 10% sodium dodecyl sulfate, 5.0 ml 1 N NaOH, and 92.5 ml Milli-Q water). The mixture was heated at 95°C for 15 min and cooled immediately on ice, and 180 µl sterile Milli-Q water was added after a short centrifugation (1 s) at 13,000 x g. Subsequently, the mixture was centrifuged at 13,000 x g for 5 min and stored at –20°C. If the alkaline lysis method did not produce high quality (GTG)5-PCR profiles, a phenol-chloroform method was used (13) and DNA work solutions of 50 ng/µl were prepared. For rep-PCR, the oligonucleotide primer (GTG)5 (5'-GTGGTGGTGGTGGTG-3') was used under the PCR conditions previously described (35). PCR amplifications were performed in a Perkin-Elmer 9600 PCR apparatus as previously described (19). PCR products were separated in a 1.5% agarose gel (15 by 20 cm) for 16 h at a constant voltage of 55 V in 1x TAE (40 mM Tris-acetate, 1 mM EDTA, pH 8.0) at 4°C. (GTG)5-PCR profiles were visualized under UV light after being stained with ethidium bromide, followed by digital image capturing using a charge-coupled-device camera. The resulting fingerprints were analyzed by using the BioNumerics version 4.61 software package (Applied Maths, Sint-Martens-Latem, Belgium). The similarities among the digitized profiles were calculated by using the Pearson correlation (expressed for convenience as a percentage-of-similarity value) and a dendrogram was derived from the profiles by using the unweighted pair-group method using arithmetic averages.
Phenylalanyl-tRNA synthase (pheS) gene sequencing of selected LAB isolates.
For pheS sequence analysis, genomic DNA was obtained by alkaline lysis (see above). If the alkaline lysis method did not produce good-quality pheS sequence data, a phenol-chloroform method was used (13). The amplification conditions and sequencing reactions, using PheS-21-F and PheS-23-R, were those described by Naser et al. (24). The taxonomic reference framework of pheS sequences consisted of sequences of 255 Lactobacillus strains (representing 94 species), 51 Enterococcus strains (representing 9 species), 3 Lactococcus strains (representing 2 species), 41 Leuconostoc strains (representing 12 species), 26 Pediococcus strains (representing 9 species), and 20 Weissella strains (representing 10 species). Each of these species was represented by its type strain and, generally, by two to five additional reference strains. Newly obtained pheS sequences were imported into the BioNumerics version 4.61 software (Applied Maths), aligned, and compared using the neighbor-joining method with sequences available from the pheS database.
16S rRNA gene sequencing.
Genomic DNA was extracted using the alkaline lysis method (see above) or a phenol-chloroform method (13). 16S rRNA gene amplification, purification, and sequencing were performed as described by Vancanneyt et al. (33). Phylogenetic analysis was performed using the software package BioNumerics version 4.61 (Applied Maths). To determine the closest relatives of 16S rRNA gene sequences, a search was performed in GenBank using the BLAST algorithm.
DNA-DNA hybridization experiments.
For the extraction of high-molecular-mass DNA, strains were cultivated in MRS-5 broth. DNA-DNA hybridizations were performed as described previously (28), at a temperature of 37°C in the presence of 50% formamide.
Data analysis.
All data processing and statistical analysis were performed using the software package BioNumerics version 4.61 (Applied Maths). The numerical data input was handled as an open character type, meaning that a nonfixed number of characters (i.e., the bacterial species) were allowed. For each sourdough sample, a microbial community profile was composed that reflects the qualitative (number of species) and quantitative (number of isolates per species) diversity of sourdough LAB in the corresponding sample. Similarities were expressed using the Pearson product-moment correlation coefficient, and a dendrogram was constructed using the unweighted pair-group method using arithmetic averages. In addition, principal component analysis (PCA) was carried out to visualize the possible impact of several technological characteristics on the sourdough LAB biodiversity.
|
|
|---|
Polyphasic identification of sourdough LAB isolates.
(GTG)5-PCR fingerprinting was used to screen 714 presumptive LAB isolates. Numerical analysis of the digitized (GTG)5-PCR fingerprints with the reference database resulted in the delineation of 24 clusters and six single isolates at 50% Pearson similarity. Visual inspection of (GTG)5-PCR fingerprints indicated that the variations in the band patterns within one rep-PCR cluster were mainly due to differences in band intensities. However, minor qualitative differences revealed heterogeneity among isolates from sourdoughs produced at different bakeries, suggesting that such isolates represent different strains. Isolates were tentatively assigned to a given species when they belonged to a (GTG)5-PCR cluster that contained the type strain and/or one or more additional reference strains. This allowed the tentative identification of 514 LAB isolates (72% of the total) belonging to 10 different (GTG)5-PCR clusters. The reproducibility of each run was evaluated by the inclusion of control strain Lactobacillus alimentarius LMG 21683. No marked qualitative variations between runs were observed for this strain (data not shown).
Subsequently, one or more representative strains from each (GTG)5-PCR cluster were selected and subjected to pheS sequence analysis (a total of 96 isolates). The GenBank/EMBL/DDBJ accession numbers for the pheS gene sequences of the selected isolates are AM745639 to AM745731, AM285025, AM285026, and AM259120. For the identification of Lactobacillus and Enterococcus species, the pheS gene generally offers an intraspecies interval with a variation of up to 3% for both genera and interspecies gaps exceeding 10% and 16% divergence, respectively (90% confidence interval) (24; Naser et al., unpublished results). Data on the pheS intra- and interspecies variations for Pediococcus, Leuconostoc, and Weissella species are currently not available. Alignment and cluster analysis of the pheS sequences of the selected LAB sourdough isolates and reference strains resulted in the delineation of 25 pheS clusters at a similarity level of 97% (see Fig. S1 in the supplemental material). In this way, 91 LAB sourdough isolates (95%) belonging to 24 different pheS clusters could be identified at the species level. One pheS cluster that showed the highest pheS sequence similarity (94.2%) with the type strain of Lactobacillus rossiae (LMG 22972T) could not be unequivocally allocated to a LAB species included in the reference database. Based on 16S rRNA gene sequencing, isolate D03WW01T01-35 of this "unknown" pheS cluster showed 99.5% 16S rRNA gene sequence similarity with L. rossiae LMG 22972T. A DNA-DNA binding value of 83% was found between the type strain of L. rossiae (LMG 22972T) and strain D03WW01T01-35, which is well above the proposed threshold of 70% for species delineation (30). It thus indicates that strain D03WW01T01-35 and all other members of the unidentified pheS cluster belong to L. rossiae.
Taxonomic distribution and prevalence of sourdough LAB species.
Seven hundred and fourteen LAB were randomly isolated from 21 Belgian traditional sourdoughs. By polyphasic identification, 82% of the isolates were assigned to seven species: Lactobacillus paralimentarius (21%), Lactobacillus sanfranciscensis (17%), Lactobacillus plantarum (11%), Lactobacillus pontis (10%), Lactobacillus fermentum (9%), Lactobacillus hammesii (8%), and Weissella cibaria (6%). The heterofermentative species L. paralimentarius, L. sanfranciscensis, L. plantarum, and L. pontis occurred in 10, 9, 9, and 9 sourdoughs, respectively, and are therefore considered the dominant cultured microbiota in Belgian sourdoughs. However, all four species were never found together in any of the sourdoughs analyzed, and only in samples D01WW01T01 and D06WW01T01 was an association of three of the species observed. Furthermore, the occurrence of the sourdough species L. sanfranciscensis and L. paralimentarius seemed to be negatively correlated (Spearman correlation value,
= –0.656; P = 0.000625). L. fermentum, L. hammesii, and W. cibaria occurred in three, five, and two sourdoughs, respectively.
The remaining 129 isolates (18%) were assigned to LAB species Lactobacillus sakei (n = 26), Lactobacillus brevis (n = 15), Lactobacillus helveticus (n = 15), L. rossiae (n = 11), Lactobacillus curvatus (n = 11), Pediococcus pentosaceus (n = 9), Lactobacillus parabuchneri (n = 8), Lactobacillus spicheri (n = 8), Lactobacillus nantensis (n = 7), Lactobacillus namurensis (n = 6), Lactobacillus crustorum (n = 3), Pediococcus acidilactici (n = 3), Leuconostoc mesenteroides (n = 2), Weissella confusa (n = 2), Enterococcus mundtii (n = 1), Lactobacillus buchneri (n = 1), and Lactobacillus paracasei (n = 1). The microbiota of the sampled sourdoughs consisted of two to seven different LAB species, with the exception of D02WR01T02, from which only L. sanfranciscensis was isolated. The recently described Lactobacillus species L. namurensis (28) and L. crustorum (29) were isolated from one (D09ME01T01) and two (D01WW01T01 and D05WW01T01) sourdoughs, respectively. Notably, the occurrence of L. crustorum and L. helveticus seemed to be positively correlated in the sampled sourdoughs (Spearman correlation value,
= 0.995; P < 0.01).
To study the influence of descriptive and technological characteristics on the LAB diversity and the relationships between dominant LAB species in the sourdough samples, cluster and PCA analyses were performed. Seven bacterial species, corresponding to the seven most-dominant species that together represent 82% of the total LAB population, were included in the analyses (Fig. 1). Except for sample D07WR02T01 from depositor 7, all sourdough samples produced at one particular bakery were found together in separate groups and tended to show a similar qualitative and quantitative microbial diversity. As an example, all sourdough samples originating from depositor 11 and produced with different types of flour (rye, spelt, and two varieties of wheat flour) clearly grouped together and are all dominated by the species L. paralimentarius, L. plantarum, and L. hammesii (Fig. 1). On the other hand, a number of sourdoughs from different bakeries (e.g., those from depositor 5 and 10) displayed highly similar LAB compositions.
![]() View larger version (30K): [in a new window] |
FIG. 1. Biodiversity data were imported into BioNumerics version 4.61 as an open character type, resulting in microbial community profiles reflecting the diversity of dominant LAB species in 21 Belgian traditional sourdoughs analyzed. Descriptive (depositor and type of flour) characteristics are shown; technological characteristics (pH and fermentation time and temperature) are in Table 1. The dendrogram was constructed by using the unweighted-pair group method using arithmetic averages, with correlation levels expressed as percentage-of-similarity values of the Pearson correlation coefficient. Bootstrap values (based on 1,000 replications) at or above 50% are indicated at the branch points.
|
![]() View larger version (84K): [in a new window] |
FIG. 2. Dendrogram showing the genetic similarities between selected strains of the four most frequently isolated species from Belgian artisan sourdoughs (L. paralimentarius, L. sanfranciscenis, L. plantarum, and L. pontis) based on (GTG)5-PCR fingerprints. One representative fingerprint pattern per species isolated from sourdough samples originating from depositors D2, D6, D7, D10, and D11 and the pattern of the corresponding type strain (in bold face) were included for cluster analysis. The dendrogram was constructed using the unweighted-pair group method using arithmetic averages, with correlation levels expressed as percentage-of-similarity values of the Pearson correlation coefficient.
|
![]() View larger version (57K): [in a new window] |
FIG. 3. Three-dimensional score plot of the PCA of the data obtained from 21 Belgian artisan sourdough samples. 1.1, D01WW01T01; 2.1, D02WR01T01; 2.2, D02WR01T02; 2.3, D02WW01T01, 3.1, D03WW01T01; 4.1, D04WW01T01; 5.1, D05WW01T01; 6.1, D01SS01T01; 6.2, D01WW01T01; 7.1, D07WR01T01; 7.2, D07WR01T02; 7.3, D07WR02T01; 8.1, D08WW01T01; 9.1, D09ME01T01; 10.1, D10WW01T01; 10.2, D10SS01T01; 10.3, D10WR01T01; 11.1; D11RR01T01, 11.2, D11WW01T01; 11.3, D11WW02T01; 11.4, D11SS01T01.
|
|
|
|---|
Consistent with the results of other European sourdough biodiversity studies, heterofermentative species appear to dominate the microbiota of Belgian type I sourdoughs (6, 9, 14, 36), although complex associations of homo- and heterofermentative strains were found in 12 of the 21 sourdoughs. All Belgian sourdough samples were dominated by one or more of the heterofermentative species L. paralimentarius, L. sanfranciscensis, L. plantarum, and L. pontis. The relative dominance of L. sanfranciscensis in type I sourdoughs has been ascribed to environmental conditions during sourdough fermentation (6, 8, 27, 37). In contrast to these studies, in which the frequent association of L. sanfranciscensis and L. plantarum is reported, we detected both species in only 2 of the 21 Belgian sourdoughs. Likewise, Corsetti et al. (6) suggested an association of L. sanfranciscensis and L. alimentarius in Italian wheat sourdoughs, but no members of the latter species were found in the present study. However, the closely related L. paralimentarius (3) is common to traditional Greek wheat sourdoughs (9) and was isolated from 10 Belgian sourdoughs produced at six different bakeries. In addition, our results suggest that the copresence of L. sanfranciscensis and L. paralimentarius in Belgian sourdough samples is negatively correlated. Interestingly, the negative effect of L. paralimentarius on L. sanfranciscensis growth has been previously reported (25). Another characteristic of Belgian sourdoughs is the presence of the obligate heterofermentative species L. pontis, which is frequently isolated from European sourdoughs (37, 38) and was found in 9 of the 21 Belgian sourdoughs sampled. However, L. pontis has been isolated from pig intestines as well, indicating that L. pontis is not associated only with sourdough (18). The rare presence of the recently described LAB species L. namurensis (28) and L. crustorum (29) in one and two Belgium sourdoughs, respectively, may be the result of random colony picking or nonoptimal culture conditions (8). Hence, culture-independent studies using denaturing gradient gel electrophoresis and/or real-time PCR are expected to yield more information about the relative abundance and ecology of these and other species in Belgian sourdoughs.
Although it is known that the type of flour, process technology, and other factors strongly influence the composition of the sourdough microbiota (9, 21), cluster analysis and PCA indicate that the typical microbial composition of Belgian sourdoughs is influenced by the bakery environment, rather than the type of flour used to produce the sourdough. Individual exceptions to this conclusion may be the result of arbitrary colony selection. Although it was beyond the scope of this study to investigate the clonal relatedness of the sourdough isolates in depth, (GTG)5-PCR fingerprint analysis gave evidence of high infraspecific similarity among isolates from sourdoughs produced at the same bakery. These results may indicate that the "in-house microbiota" of the bakery setting largely determines the microbial diversity in sourdoughs.
In conclusion, the combined use of (GTG)5-PCR fingerprinting and pheS sequence analysis proved to be a valuable method for the screening and identification of a large set of LAB sourdough isolates. This polyphasic approach was applied for the analysis of 21 Belgian sourdough samples originating from 11 different artisan bakeries located throughout Belgium. Large differences in species composition were found, except in sourdoughs from the same bakery, which demonstrated a similar qualitative and quantitative LAB population. These findings suggest that the bakery environment significantly influences the LAB community structure of Belgian sourdoughs, irrespective of the type of flour used to prepare the sourdough.
We thank the owners and staff of the bakeries for providing the sourdough samples used in this study. We also thank Bruno Pot for helpful discussions on statistical analysis.
Published ahead of print on 3 August 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»