Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2005, p. 7872-7879, Vol. 71, No. 12
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.12.7872-7879.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Área de Microbiología, Departamento de Ciencias de la Salud, Facultad de Ciencias Experimentales, Universidad de Jaén, 23071 Jaén, Spain
Received 7 March 2005/ Accepted 15 August 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Conventional culture-dependent methods of strain characterization by physiological and biochemical tests and molecular techniques, such as ribotyping, amplified fragment polymorphism, and randomly amplified polymorphic DNA, can give significant insight on specific isolates and microbial populations occurring in caper berry manufacture. However, cultivation can over- or underestimate the microbial diversity, as media may not be sufficiently selective for monitoring population dynamics. Novel molecular approaches, using molecular methods such as PCR amplification of the 16S rRNA gene in combination with denaturing gradient gel electrophoresis (DGGE) or temperature gradient gel electrophoresis is commonly performed in microbial ecology (2-5, 15, 16). These methods have provided a new insight into microbial diversity and allow a more rapid, high-resolution description of microbial communities than that provided by the traditional approach of isolation of microorganisms. The use of 16S rRNA gene and PCR-DGGE or PCR-temperature gradient gel electrophoresis has often been combined with sequencing and subsequent identification of the species present in a sample.
In this study, the dynamics of the microbial community responsible for the fermentation of caper berries was investigated for the first time by using a polyphasic approach combining microbial enumerations with culture media, randomly amplified polymorphic DNA (RAPD) and temporal temperature gradient gel electrophoresis (TTGE) fingerprinting of total community DNA and sequencing of partial 16S rRNA genes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Enumeration of microorganisms.
Serial dilutions of caper samples were used for microbial enumerations with the following media: plate count agar (Scharlab, Barcelona, Spain) for estimation of total aerobic mesophilic bacteria; MRS agar (Scharlab) for lactic acid bacteria (LAB); McConkey agar (Scharlab) for enterobacteria; potato dextrose agar (Merck) containing (per liter) 14 mg of tartaric acid (Sigma), 50 mg of chloramphenicol (Sigma), and 50 mg of rose bengal (Sigma) for yeasts and molds. All enumerations were done by plate counting. Portions (0.1 ml) of appropriate dilutions were spread plated in triplicate. Counts were obtained after 48 h and 5 days of incubation at 30°C. Results were calculated as the means of three determinations.
For LAB isolation, the plates were incubated under anaerobic conditions at 32°C. A total of 133 colonies were randomly picked from MRS agar plates. The isolates were purified, checked for catalase activity and gram reaction, and microscopically examined prior to storage in liquid culture using 20% glycerol at 80°C.
Isolates with coccoid morphology were tested for the capacity to grow in brain heart infusion (BHI) broth (Scharlab) at temperatures of 10°C and 45°C. The capacity to grow at 37°C in BHI broth adjusted to pH 9.6 with NaOH or in BHI supplemented with NaCl to 6.5% final concentration as well as growth and esculin hydrolysis on bile esculin agar (Scharlab) were also tested.
Preparation of DNA from pure culture.
Total DNA was extracted from overnight cultures by the method of De los Reyes-Gavilan et al. (13). Two-milliliter samples were centrifuged at 7,000 x g for 10 min. The cell pellets were resuspended in 200 µl of a solution of lysozyme (20 mg/ml) in TES buffer (50 mM Tris-HCl, pH 8.0, 1 mM EDTA, 8.56% [wt/vol] sucrose) containing 10 µl of a mutanolysin solution (1 U/µl) and were incubated for 1 h at 37°C. Fifty microliters of 20% sodium dodecyl sulfate was then added to each preparation to lyse the cells. Total DNA was then purified by repeated extraction with phenol-chloroform-isoamyl alcohol (25:24:1), followed by final extraction with chloroform-isoamyl alcohol (24:1). The DNA was precipitated with isopropanol, washed with 70% ethanol, and vacuum dried. The resulting pellets were resuspended in 200-µl portions of a solution containing 50 µg/ml RNase and incubated for 15 min at 37°C. DNA concentration and quality were assessed by measuring optical density at 260 and 280 nm with a Beckman DU670 spectrophotometer.
DNA isolation from fermented capers.
Samples of the fermentation broth (1.5 ml each) were centrifuged for 5 min at 12,000 rpm in a microcentrifuge. The cell pellets were kept overnight at 80°C. Cells were resuspended in 400 µl TESL buffer (25 mM Tris-HCl, pH 8.0, 10 mM EDTA, 20% sucrose, and 20 mg/ml lysozyme), and 20 µl of 1-U/µl mutanolysin (Sigma) was added to each sample. Samples were vortexed for 1 min and incubated for 1 h at 37°C. Fifty microliters of a proteinase K solution (10 mg/ml) was added to each tube, and the tubes were incubated for 50 min at 50°C and then for 10 min at 65°C. To lyse the cell membranes and separate DNA from cell debris and protein, 100 µl of 20% sodium dodecyl sulfate, 700 µl of phenol, pH 8.0 (Sigma), and 100 µl of chloroform-isoamyl alcohol were then added to each sample. The samples were mixed by inversion followed by microcentrifugation for 10 min at 12,000 x g, 4°C. The aqueous phase was transferred to a new tube, and another extraction was performed with phenol and chloroform-isoamyl alcohol followed by one extraction with chloroform. The resulting aqueous phase was precipitated by the addition of isopropanol. Samples were held at 20°C for at least 1 h and centrifuged again. DNA pellets were washed twice with 70% ethanol, dried briefly under a vacuum, dissolved in 50 µl water, and stored at 80°C. The quality of the DNA recovered was routinely checked on agarose gels, and DNA was quantified spectrophotometrically as previously described.
RAPD-PCR amplification.
RAPD-PCR analysis was performed as described previously (6) using the primers M13 (5'-GAG GGT GGC GGT TCT-3') and B10 (5'-GGT GAC GCA G-3') in separate reactions. DNA was amplified with primer M13 for 35 cycles (94°C for 1 min; 40°C for 20 s, ramp to 72°C at 0.6°C/s for 2 min; 72°C for 2 min) and with primer B10 for 30 cycles (94°C for 1 min, 34°C for 1 min, ramp to 72°C at 0.6°C/s for 2 min; 72°C for 2 min). The gels were stained in ethidium bromide and photographed on a UV transilluminator. Photopositives were scanned into a computer and subsequently analyzed using the Bionumerics software, version 2.5 (Applied Maths, Kortrijk, Belgium). RAPD-PCR patterns were grouped by means of cluster analysis with the Pearson product moment correlation coefficient and the unweighted pair group method using arithmetic averages (UPGMA). The RAPD patterns obtained with both primers were analyzed together using Bionumerics software to obtain a single dendrogram.
Identification of the isolates.
Molecular identification was carried out by PCR. The isolates were identified by using species-specific primers for Enterococcus faecalis (E1/E2), Enterococcus faecium (F1/F2) (14), Lactobacillus plantarum (PLANT1/LOWLAC), and Lactobacillus fermentum (FERM1/LOWLAC) (10). Isolates that tested positive in PCR amplification for L. plantarum were then tested by multiplex PCR for identification of L. plantarum, Lactobacillus paraplantarum, and Lactobacillus pentosus with recA gene-based primers paraF (5'-GTC ACA GGC ATT ACG AAA AC-3'), pentF (5'-CAG TGG CGC GGT TGA TAT C-3'), planF (5'-CCG TTT ATG CGG AAC ACC TA-3'), and pREV (5'-TCG GGA TTA CCA AAC ATC AC-3') as described by Torriani et al. (32). The following strains were used as references: Enterococcus faecalis CECT 481T, E. faecium CECT 410T, L. plantarum CECT 220, L. paraplantarum CECT 5783, L. pentosus CECT 4023, and L. fermentum CECT 4007 (CECT strains are from the Colección Española de Cultivos Tipo [http://www.cect.org], Burjassot, Valencia, Spain). The PCR products were analyzed by electrophoresis on an agarose gel in Tris-borate-EDTA buffer at 100 V, and stained with ethidium bromide solution (0.5 µg/ml).
For 16S rRNA gene sequencing, genomic DNAs of strains were used to amplify the almost complete 16S rRNA gene. Amplification of the 16S rRNA gene sequence by PCR was performed using the following primers: 5'-AGAGTTTGATCMTGGCTC-3' (forward), corresponding to conserved Escherichia coli 16S rRNA positions 8 to 25 (8), and 5'-CNCGTCCTTCATCGCCT-3' (reverse), corresponding to conserved E. coli 23S rRNA positions 44 to 60 (9). PCR was carried out according to the method of Cibik et al. (11). The product of amplification was purified by using a Quantum Prep PCR Kleen Spin column (Bio-Rad, Madrid, Spain) and then sequenced by using the primers Sp3 (5'-TACGCATTTCACCKCTACA-3', position 684 reverse), Sp4 (5'-CTCGTTGCGGGACTTAAC-3', position 1089 reverse), and Sp5 (5'-GNTACCTTGTTACGACTT-3', position 1492 reverse) according to the method of Weisburg et al. (34) in a CEQ 2000 XL DNA analysis system (Beckman Coulter, California). The partial sequences of the 16S rRNA gene were determined by using CEQ 2000 dye terminator cycle sequencing with a quick start kit (Beckman Coulter, California) according to the manufacturer's instructions. The sequence data were analyzed with a CEQ DNA analysis system (version 4.0). The overlapping sequences obtained with SP3; SP4 and SP5 were merged using theLaser gene program, version 5.05 (DNASTAR, Inc., Madison, WI). A search for homology of the DNA sequence was done using the BLAST algorithm available at the National Center for Biotechnology Information.
PCR amplification of the microbial community 16S rRNA gene.
TTGE samples were prepared by performing two successive PCR amplifications (nested PCR) by using the primer pairs described elsewhere (24). First, a 700-bp fragment of the 16S rRNA gene including the V3 region was amplified. The reaction mixture (100 µl) contained reaction buffer (10 mM Tris-HCl, 1.5 mM MgCl2, 50 mM KCl [final concentrations]), each deoxynucleoside triphosphate at a concentration of 200 µM, 60 pmol of primer W01 (5'-AGA GTT TGA TC[A/C] TGG CTC-3'), 60 pmol of primer W012 (5'-TAC GCA TTT CAC C[G/T]C TAC A-3'),
50 ng of total DNA, and 2.5 U of Taq DNA polymerase (Amersham Biosciences). The amplification program was 96°C for 2 min; 30 cycles of 96°C for 1 min, 50°C for 30 s, and 72°C for 1 min; and finally, 72°C for 2 min. Second, the 700-bp fragment was used to amplify the V3 region with the following primers (24): HDA1-GC (5'-CGC CCG GGG CGC GCC CCG GGC GGG GCG GGG GCA CGG GGG GAC TCC TAC GGG AGG CAG CAG T-3', the GC clamp sequence is in bold) and HDA2 (5'-GTA TTA CCG CGG CTG CTG GCA-3'). The reaction mixture (100 µl) consisted of a reaction buffer (10 mM Tris-HCl, 1.5 mM MgCl2, 50 mM KCl [final concentrations]), each deoxynucleoside triphosphate at a concentration of 200 µM, 60 pmol of each primer, 1 µl of the amplified 700-bp fragment, and 2.5 U of Taq DNA polymerase. The amplification program was 94°C for 2 min; 30 cycles of 94°C for 1 min, 58°C for 30 s, and 72°C for 1 min; and finally, 72°C for 7 min. The sizes and quantities of PCR products were determined by 2% agarose gel electrophoresis.
TTGE analysis.
PCR products obtained from the V3 region amplification as described above were subjected to TTGE analyses. TTGE was performed by using the Dcode universal mutation detection system (Bio-Rad). Polyacrylamide (8%) gels were prepared and run with 1x TAE buffer (40 mM Tris base, 20 mM glacial acetic acid, 1 mM EDTA, pH 8.3) diluted from 50x TAE buffer. Gels were prepared with 8% (wt/vol) acrylamide stock solutions (acrylamide-bisacrylamide, 37.5:1) and a final urea concentration of 7 M. The final electrophoresis conditions were 70 V for 16 h with an initial temperature of 63°C and a final temperature of 70°C (the temperature was increased 0.4°C per h). Samples (20 µl) of PCR products were deposited in wells. A magnetic stirrer was used to mix the buffers and improve the temperature gradient homogeneity. After runs, gels were stained for 30 min with ethidium bromide (0.5 µg/ml of 1x TAE buffer), rinsed for 20 min in 1x TAE buffer, and photographed on a UV transillumination table.
DNA sequencing and analysis of sequence data.
TTGE fragments were excised with a sterile scalpel. The DNA of each fragment was eluted in 20 µl of sterile water overnight at 4°C. One microliter of the eluted DNA of each TTGE band was reamplified using primers HDA1-GC and HDA2. The success of this procedure was checked by electrophoresing 3-µl portions of the PCR products in TTGE gel as described above with caper-amplified DNA as control.
To determine the closet known relatives of the partial 16S rRNA gene sequences obtained, searches of public data libraries (RDP and GenBank) were performed by using the BLAST and RDP programs.
| RESULTS |
|---|
|
|
|---|
|
Identification and cluster analysis of bacilli.
By means of PCR amplification with species-specific probes having rRNA genes as targets, 6 isolates were identified as Lactobacillus fermentum (Table 1) and 43 isolates tested positive for Lactobacillus plantarum. The later were further identified by multiplex PCR amplification of recA gene-based primers as L. plantarum (37 isolates), L. paraplantarum (1 isolate), and L. pentosus (5 isolates) (Table 1). The rest of the isolates could not be amplified by PCR and were identified as Lactobacillus brevis by 16S rRNA gene sequencing.
|
|
Group G2 includes all six isolates of L. fermentum (P15g, P15p, P18g, Lb44p, Lb45, and Lb34). These are genetically distant and show a similarity coefficient of 65.51%.
Group G3 includes 23 isolates that cluster in three subgroups (G3-1 to G3-3). Subgroup G3-1 is composed of three isolates of L. brevis (L9, L13, and Lb14g) with an 80% similarity coefficient. Subgroup G3-2 includes 15 isolates genetically close that show a similarity coefficient of 83.13%: 11 isolates of L. plantarum (Lb17, Lb38p, Lb35, Lb31g, Lb32g, Lb14p, Lb25, Lb26, Lb27, Lb29, and Lb23), 3 isolates of L. pentosus (Lb19p, P19g, and P19p), and the single isolate L. paraplantarum Lb20g. Subgroup G3-3 is the most homogeneous and includes 5 isolates of L. plantarum (Lb1g, Lb3, Lb2g, Lb5p, and Lb4g) that are genetically very close and show a similarity coefficient of 91.89%.
Group G4 is more heterogeneous and includes six isolates of L. brevis (L16, Lb21, Lb13, Lb15, L15, and L17) with a similarity coefficient of 73%.
Identification and cluster analysis of cocci.
Preliminary physiological tests indicated that all isolates were able to grow at 45°C, but only two isolates (S1 and S11) were able to grow at 10°C as well as in media containing 6.5% NaCl or buffered at pH 9.6. They were also able to hydrolyze esculin. Both isolates were identified as Enterococcus faecium by PCR amplification with specific probes (Table 1). The remaining isolates were unable to grow at 10°C or at pH 9.6. Fourteen of them (S4, S5, S7, S8, S9, S10, S12, S13, S14, P1, P2, P3, P6, and P9) were able to hydrolyze esculin and to grow in media containing 6.5% NaCl. They were identified as Pediococcus pentosaceus by 16S rRNA gene sequencing (Table 1). Finally, isolate S15 was unable to hydrolyze esculin and showed a weak growth capacity in the presence of 6.5% NaCl. This isolate was identified as Pediococcus acidilactici by 16S rRNA gene sequencing.
Cluster analysis of RAPD-PCR patterns obtained for cocci revealed two main groups, G1 and G2 (Fig. 3). Group G1 includes all isolates of P. pentosaceus distributed in two subgroups plus the single isolate of P. acidilactici (S15). Isolates included in subgroup G1-1 (S4, P9, P2, P6, P1, P3, S7, S5, S10, S9, and S12) are genetically close and show a similarity coefficient of 89.55%, while isolates in subgroup G1-2 (S13, S14, and S8) are genetically very close and show a higher similarity coefficient of 93.77%. Group G2 includes E. faecium isolates (S1 and S11) with a similarity coefficient of 89.24%.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
A collection of bacterial isolates obtained during the fermentation period was used as starting material for this study. Biochemical characterization, although useful for identification purposes, has a poor discriminating capability, and molecular techniques must be used with typing purposes. Of the various typing methods available, those based on PCR product analysis, such as RAPD technique, are reported to be simple and rapid to perform and to provide good levels of discrimination; they are therefore considered the most suitable when a large number of strains must be analyzed. After comparison of RAPD profiles and elimination of repeated isolates coming from the same day of fermentation with identical RAPD profiles, a total of 58 bacilli and 17 cocci were selected for species identification and cluster analysis. Following molecular identification methods using species-specific PCR primers and 16S rRNA gene sequencing, the 75 isolates were identified as L. plantarum, L. paraplantarum, L. pentosus, L. brevis, L. fermentum, P. pentosaceus, P. acidilactici, and E. faecium. About 49% of isolates were L. plantarum. This species showed greatest genetic diversity by RAPD-PCR analysis. This greater genotypic differentiation could be the result of the higher number of isolates belonging to these species used for comparison, which would increase the likelihood of finding more distantly related taxonomic units.
Results from this study indicate that caper fruit fermentation shares many characteristics with other vegetable fermentations, where L. plantarum is usually the predominant species as well, both because of its higher acid tolerance and its ability to degrade sugars which are present in vegetables (12, 17, 19, 23, 30). However, while Leuconostoc species are also present in other vegetable fermentations, they were not detected in caper fermentations. A similar profile has been described in Almagro eggplant fermentation, where Lactobacillus species were also predominant (31). The absence of Leuconostoc could be attributed to the rapid decrease of pH, since this bacterium is not able to grow below pH 4.5 (23). It should be noticed as well that this process takes place during a very warm time in the year, and the exposure of the containers to the solar radiation allows the fermentation to take place at temperatures over 40°C. On the contrary, growth of Leuconostoc appears be favored in fermentations that take place at low temperatures, between 8 and 18°C (29), and therefore, the high temperatures in caper fermentations could very well inhibit their growth.
The application of molecular methods in the field of microbiology allows a better understanding of the ecology of food fermentations. Since the results are simple to interpret and mistakes in the identification of species due to the problems related to biochemical tests are avoided, dominant strains responsible for the main transformations during fermentations can be easily detected and identified. The use of culture-independent methods is a recent trend in the study of the microbial ecology of foods (18). PCR-DGGE fingerprinting has been used to study microbial communities in different fermented foods including pozol, cassava, fermented sausages, and dairy products (15). TTGE, which has been used to characterize other types of microbial populations (7, 24), is a technique that works on the same basic principle than DGGE but without the requirement for a chemical denaturing gradient, thus producing more reproducible data (35). In the present work, TTGE was used to monitor the population dynamics during natural fermentation of capers. DNA was amplified by nested PCR, and gels were visually inspected to identify the bands representing the populations involved in the fermentation. To circumvent the biases inherent in subjective interpretation, the presence of the bands was confirmed by direct sequencing. When the results obtained from both traditional plating and TTGE were analyzed, it became evident that the fermentation was characterized by strong LAB activity. The profiles obtained by TTGE allowed us to differentiate the main groups of lactic acid bacteria (L. plantarum, L. brevis, L. fermentum, Pediococcus sp., Enterococcus sp., E. faecium, and E. casseliflavus). These results also show that the predominant species in this fermentation process is L. plantarum, which was represented by a highly intense band throughout the fermentation.
Comparison of results from molecular and culture-dependent studies revealed similar species composition for caper fermentation, and all species of Lactobacillus were correctly identified in both cases. However, Pediococcus could only be identified at species level by culture-dependent methods, while TTGE provided a more accurate picture of the distribution of this bacterium through the fermentation. E. faecium was also identified in both cases, but once more, its distribution during fermentation could only be determined by TTGE. Another enterococcal species (E. casseliflavus) could only be detected by TTGE analysis. These results suggest that TTGE analysis can offer better results for detection of species that are present at lower concentrations. However, it should also be considered that amplification results may largely be dependent on the extraction efficacy for different bacterial species and food matrices. Results from other works also revealed significant differences in the microbial compositions of fermented foods depending on the use of culture-dependent or culture-independent methods (5, 16, 21, 22), suggesting that polyphasic studies should be used to better understand the microbial ecology of foods.
The results presented in this work revealed the microbial complexity of caper berry fermentation. This would be expected for a natural fermentation, which relies solely on the microbiota of the fruit surface, manufacturing plant, and other natural sources. Once the predominant species have been identified, bacterial isolates with selected properties could specifically be used for development of starter cultures to avoid problems of stuck fermentations and to achieve a more homogeneous flavor of the fermented product.
| ACKNOWLEDGMENTS |
|---|
We also acknowledge the Technical Services of this University for support with DNA sequencing.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|