|
|
||||||||
Departamento de Microbiología, Instituto de Fermentaciones Industriales, CSIC, Juan de la Cierva 3, 28006 Madrid, Spain
Correspondence
Rosario Muñoz
rmunoz{at}ifi.csic.es
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are AJ966402AJ966404 (pgm), AJ966405AJ966409 (ddl), AJ966364AJ966370 (gyrB), AJ966371AJ966378 (purK1), AJ966379AJ966388 (gdh), AJ966389AJ966396 (mutS) and AJ966397AJ966401 (tkt4).
| INTRODUCTION |
|---|
|
|
|---|
Identification of L. plantarum strains is essential in both basic and applied research. Intraspecific differentiation is an important preliminary step for the selection of starter cultures, because technological, probiotic, antimicrobial and sensorial attributes are strain-specific and it may help to distinguish strains with particular technological properties. Currently, a great number of mostly molecular techniques are available for the identification of LAB, for industrial processes and food products. For each specific type of research or analysis, a well considered choice has to be made regarding the methodology to be applied in relation to taxonomic resolution, workload and cost. It is important to bear in mind that not every technique can be used for every purpose. In the course of safety assessments, it is crucial to use techniques working at the strain level to obtain a detailed fingerprint of individual isolates. The increasing interest in some L. plantarum properties, for example probiotic activities of specific strains (Herias et al., 1999
), contributes to the need for a reliable molecular method for the definite identification of L. plantarum at the strain level. The need for positive identification of different isolates is also acknowledged by research workers in the field, since many strains from diverse origins are often exchanged between laboratories and no reliable phenotypic method for certifying their identities is available. The identification of L. plantarum at the strain level is also important for industrial use. The biotechnological industry needs tools for monitoring, for example, the use of patented strains or to distinguish probiotic strains from natural isolates in the host gastrointestinal tract.
Genotypic methods used for L. plantarum strain typing are typically PCR-based methods, for example randomly amplified polymorphic DNA (Johansson et al., 1995
; Oneca et al., 2003
; Elegado et al., 2004
) or variations of restriction enzyme analysis, for example ribotyping (Yansanjav et al., 2003
; Rodas et al., 2005
) and pulsed-field gel electrophoresis (Sánchez et al., 2004
). Multilocus sequence typing (MLST) has recently been shown to be a powerful technique for bacterial typing. MLST makes use of automated DNA sequencing to characterize the alleles present at different housekeeping gene loci. As it is based on nucleotide sequence, it is highly discriminatory and provides unambiguous results that are directly comparable between laboratories. The MLST method was first described in 1998 and since then it has been applied to important bacterial pathogens, including several food-borne human pathogens such as Campylobacter jejuni (Dingle et al., 2001
), Vibrio cholerae (Farfán et al., 2002
) and Bacillus cereus (Helgason et al., 2004
); recently, MLST was also applied to the non-pathogenic food production bacterium Oenococcus oeni (De las Rivas et al., 2004
).
The aims of the present study were to: (i) develop an MLST method for L. plantarum; (ii) compare the discriminatory power of ribotyping, restriction fragment length polymorphism (RFLP) of the 16S23S rDNA intergenic spacer region (ISR) and MLST in this species; and (ii) use MLST to analyse L. plantarum population structure.
| METHODS |
|---|
|
|
|---|
|
RFLP of the PCR-amplified 16S23S rDNA ISR.
RFLP analysis of the ISRs was performed by using primers 16S14F and 23S1R based on conserved areas of aligned rRNA bacterial sequences (Zavaleta et al., 1996
). These primers amplified a 550 bp fragment in all the L. plantarum strains tested. PCR was performed in a volume of 50 µl as described by Zavaleta et al. (1996)
. The amplified 16S23S ISRs from L. plantarum strains were digested with restriction enzymes AluI, CfoI, DdeI and TaqI, since, by searching in the L. plantarum strain WCFS1 sequence (NC_004567), recognition sites for these enzymes were found in the 16S23S ISR. The digested products were separated by electrophoresis in 4·5 % MS-8 agarose gels.
Amplification and nucleotide sequencing.
The housekeeping genes encoding the following proteins were chosen for analysis: phosphoglucomutase (pgm), D-alanine-D-alanine ligase (ddl), B subunit of DNA gyrase (gyrB), ATPase subunit of phosphoribosylaminoimidazole carboxylase (purK1), glutamate dehydrogenase (gdh), DNA mismatch repair protein (mutS) and transketolase (tkt4). The DNA sequences of these candidate loci are available from GenBank (Table 2
). These housekeeping genes, which are not subject to any unusual selective forces and which diversify slowly by the random accumulation of neutral or nearly neutral variation, should provide reliable information about the relationships among isolates. To avoid the chance association of alleles, these genes were also selected on the criterion that they are widely separated on the chromosome, except purK1 and gdh which are only 28·5 kb apart (Kleerebezem et al., 2003
).
|
Data analysis.
For each locus, the sequences obtained for all isolates were compared and allele numbers were assigned to each unique sequence. Each isolate was defined by an allele profile or sequence type (ST) derived from the combination of numbers corresponding to the alleles at the loci analysed (Table 1
). Sequences which were different even at a single nucleotide site were considered distinct alleles.
Sequence alignments and comparison were done with the program BioEdit (http://jwbrown.mbio.ncsu.edu/BioEdit/bioedit.html) (Hall, 1999
) and converted into MEGA and NEXUS files with START (Sequence Type Analysis and Recombinatorial Tests; http://outbreak.ceid.ox.ac.uk/software.html). Phylogenetic trees were compiled with MEGA version 2.1 software (http://www.megasoftware.net). Dendrograms were constructed by the unweighted pair group method with arithmetic averages (UPGMA) with the Kimura two-parameter model. The percentage of bootstrap confidence levels for internal branches, as defined by the MEGA program (Kumar et al., 1994
), was calculated from 1000 random resamplings.
The method of split decomposition was used to assess the degree of tree-like structure present in the alleles found for each locus in the complete set of 16 isolates (Huson, 1998
). The sequence alignments were converted to NEXUS files and the split decomposition was performed with SPLITSTREE 2.0 (http://bibiserv.techfak.uni-bielefeld.de/splits/). Using the START program we performed two recombination tests, the index of association (IA) (Maynard Smith et al., 1993
) and the Sawyer's run test, and a test for selection, the ratio of mean non-synonymous substitutions per non-synonymous site/mean synonymous substitutions per synonymous site, the dN/dS ratio (Nei & Gojobori, 1986
). To avoid dependence on the number of loci, we also calculated the standardized IA (
) (Haubold & Hudson, 2000
).
| RESULTS |
|---|
|
|
|---|
|
Variation at the MLST loci
The sequences of the seven chosen loci were determined for the 16 strains, with the exception of the tkt4 locus which could not be amplified from strains CECT 223, 224 and 4645. In the strains where the tkt4 locus could be amplified, we sequenced a 567 bp fragment having six polymorphic sites that generate five different alleles (data not shown). Since the tkt4 locus could not be amplified from all the L. plantarum strains, it was discarded from the MLST scheme.
The MLST scheme defined alleles between 414 (gdh) and 704 bp (gyrB) in length (Table 3
). Between 3 (pgm) and 10 (gdh) alleles were found per locus (Table 3
, Fig. 2
). The sequence differences in each allele are shown in Fig. 2
. The proportion of variable sites present in the MLST alleles ranged from 1·03 (ddl) to 7·72 % (gdh). The mean G+C content of the different gene fragments varied from 42·7 (ddl) to 50·5 mol% (gdh). Since the G+C content of the L. plantarum chromosome is 44·5 mol% (Kleerebezem et al., 2003
), the gdh allele sequences showed an unexpectedly high G+C content.
|
|
The allele frequencies showed that for most of the loci, one allele (in gyrB, mutS and purK1) or two alleles (in ddl and pgm) (Fig. 2
, Table 1
) were dominant in the population.
A UPGMA dendrogram showing the genetic relatedness among L. plantarum strains investigated in this study is shown in Fig. 3
. There were no significant clusters that correlated with the geographic origin of the strains. All strains differed in various loci, except strains RM71 and RM72, and CECT 223 and CECT 224, that shared the same ST (Table 1
, Fig. 3
).
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Ribotyping is a both an inter- and an intraspecific typing method (Manfreda et al., 2005
; Zadoks et al., 2005
) that has been successfully used to differentiate between strains of the same Lactobacillus species. Yansanjav et al. (2003)
used ribotyping in a study of five L. plantarum strains. These five strains formed two ribotype patterns and the authors found that strains isolated from different breweries shared the same ribotype. In other studies, based on these ribotypes, apparently identical pairs of L. plantarum strains were isolated in very different settings: L. plantarum ATCC 14917T (CECT 748T) was purchased by the authors from ATCC and L. plantarum UH 2153 was isolated from the vagina of a woman in Ontario (Zhong et al., 1998
). In our study we found four different ribotypes; however, most of the strains (11/15) belonged to the same ribotype, ribotype 2. These ribotype 2 strains are not source-, time- or geographically related; therefore, more specific identification methods, such as MLST would be preferable.
In addition, the variation in length and sequence of the 16S23S rDNA ISRs has been used to type strains. ISRs usually show variations which make it possible to discriminate among strains within some species. As shown in Table 1
and Fig. 1
, all the L. plantarum strains analysed in this study shared the same RFLP ISR pattern. Therefore, the RFLP ISR results support the idea of the need for a more discriminating method able to differentiate L. plantarum strains.
As a first step for developing an MLST typing method, we analysed the sequence diversity of seven housekeeping genes to ascertain if they were sufficient to provide enough typing discrimination. The internal fragments of six loci that were selected (pgm, ddl, gyrB, purK1, gdh and mutS) could be amplified from all the strains examined. The amplified internal fragments were sequenced and from these sequences we were able to use fragments between 414 and 704 bp for analysis (Table 3
). However, internal tkt4 fragments could not be amplified from strains CECT 223, 224 and 4645. From these data, we cannot exclude the possibility of a tkt4 gene deletion or the presence of a non-homologous gene copy in these strains.
The number of alleles of the six housekeeping loci ranged from 3 to 10 (Table 3
). Our results corroborate the previously described genetic heterogeneity of this species (Vescovo et al., 1993
; Sánchez et al., 2004
; Molenaar et al., 2005
). The six loci were polymorphic and most of the types were represented by a single strain. We only found two exceptions. One of them was found in strains CECT 223 and 224 that were deposited in the CECT by D. M. Alvarez Marques in 1987 and isolated in Pamplona, Spain (Table 1
). In the three typing methods used in this study, both isolates were found to be identical; therefore, it is tempting to assume that both isolates could belong to the same original strain. A similar situation was observed in strains RM71 and RM72, since both strains were isolated from the same winery in the same year and they shared identical types in all the typing methods used in this study. The number of loci used can be increased to improve resolution, but there comes a point when it is not worth the reward because little additional information will be obtained, for epidemiological purposes, compared to the cost and effort involved (Urwin & Maiden, 2003
). The proposed MLST scheme showed highly discriminating power since it was able to differentiate between highly similar strains. This was the case for strains CECT 748T and 749 that were identical in five out of the six genes analysed, differing only in the gdh locus (Table 1
).
Using this typing method we could not find an association between ST and the source from which the different strains were isolated (Table 1
). The absence of such an association may be due to the versatility of individual L. plantarum strains, which can survive or even grow in different environments. Similar results were found by Molenaar et al. (2005)
when they explored L. plantarum genome diversity by using microarrays.
Examination of the sequences of housekeeping genes from biosynthetic pathways can provide evidence for the significance of recombination, since the variation within these genes is likely to be selectively neutral. Recombination can be detected by, for example, the appearance of a network of relationships among sequences rather than a bifurcating tree-like phylogeny. The split decomposition analysis of L. plantarum strains (Fig. 4
) reveals two uncentred edges, suggesting that the evolution of these strains stems from a couple of strains from which single branches radiate.
The most simple method to detect recombination in aligned sequences is to look for mosaic structures by eye. Significant mosaic structure is indicative of recombinatorial exchange, usually among isolates of the same species (Feil et al., 2000
). In our study, four possible examples of recombination events were found: allele 3 of gyrB, allele 6 of gdh, allele 8 of purK1 and allele 3 of pgm. The mean divergence for allele 3 of gyrB (3·97 %) is much higher than the mean diversity within the other gyrB alleles (0·57 %) (Fig. 2
). This divergence is similar to the divergence observed in a possible example of an intergenic recombinatorial event in the gdh gene between Streptococcus pneumoniae and Streptococcus mitis (Enright & Spratt, 1998
). The alleles in the other loci showed a divergence of 4·11 versus 1·45 % (allele 6 of gdh), 1·90 versus 0·33 % (allele 8 of purK1) and 0·29 versus 0·09 % (allele 3 of pgm). The wide range of variable sites found in the L. plantarum loci analysed (1 % in ddl to 7·7 % in gdh) reflects variation due to point mutations through to possible horizontal transfer gene events. Recently, and mainly based on unusual base composition, horizontal transfer events from closely related organisms or even the same species have been proposed for L. plantarum strains (Molenaar et al., 2005
). The same authors have postulated that lifestyle adaptation would encompass the frequent acquisition and loss of genes involved in adaptation to niche functions.
The utility of MLST for the analysis of the genetic structure of bacterial populations is mainly based on the characteristic of housekeeping genes to have selectively neutral variability. Analysis of synonymous and non-synonymous changes in the allele sequences of a locus can be used to determine if they are subject to positive selection, so a dN/dS ratio of greater than 1 implies selection for amino acid changes. In our genetic analysis, the six housekeeping loci had dN/dS ratios significantly lower than 1 (Table 3
). Another important characteristic in relation to this is the location of the loci on the chromosome. With the exception of the purK1 and gdh loci, the loci in this study were distant enough to make joint horizontal transfer of two loci unlikely.
Three types of populations structure are known in bacteria: clonal, panmictic and epidemic. Panmictic populations may be so variable that identical strains are only found among isolates from direct contacts. In our study, only strains isolated from the same source, such as strains CECT 223 and 224, and strains RM71 and 72 share the same ST. The analysis of the L. plantarum population structure presented here suggests a panmictic, non-clonal population structure with substantial recombination. The split decomposition analysis provides strong evidence that intraspecies recombination occurs frequently in L. plantarum and plays a major role in generating genetic heterogeneity among strains. A low
value (0·138) is indicative of a weak clonal population and also confirms the importance of recombination in L. plantarum and supports the estimation that the genes investigated in L. plantarum are close to linkage equilibrium. The extension of the present analysis to a larger number of isolates could contribute to improved knowledge about the structure of L. plantarum populations.
In conclusion, the MLST scheme presented here will be a useful new tool for the precise and unambiguous characterization of L. plantarum isolates and appears to have sufficient discriminatory power for population investigations. It should allow comparison of the isolates of this species by a precise typing procedure.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
De las Rivas, B., Marcobal, A. & Muñoz, R. (2004). Allelic diversity and population structure in Oenococcus oeni as determined from sequence analysis of housekeeping genes. Appl Environ Microbiol 70, 72107219.
Dingle, K. E., Colles, F. M., Wareing, D. R. A. & 7 other authors (2001). Multilocus sequence typing system for Campylobacter jejuni. J Clin Microbiol 39, 1423.
Elegado, F. B., Guerra, M. A. R. V., Macayan, R. A., Mendoza, H. A. & Lizaran, M. B. (2004). Spectrum of bacteriocin activity of Lactobacillus plantarum BS and fingerprinting by RAPD-PCR. Int J Food Microbiol 95, 1118.[CrossRef][Medline]
Enright, M. C. & Spratt, B. G. (1998). A multilocus sequence typing scheme for Streptococcus pneumoniae: identification of clones associated with serious invasive disease. Microbiology 144, 30493060.
Farfán, M., Miñana-Galbis, D., Fusté, M. C. & Lorén, J. G. (2002). Allelic diversity and population structure in Vibrio cholerae O139 Bengal based on nucleotide sequence analysis. J Bacteriol 184, 13041313.
Feil, E. J., Enright, M. C. & Spratt, B. G. (2000). Estimating the relative contributions of mutation and recombination to clonal diversification: a comparison between Neisseria meningitidis and Streptococcus pneumoniae. Res Microbiol 151, 465469.[Medline]
Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98NT. Nucleic Acids Symp Ser 41, 9598.
Haubold, B. & Hudson, R. R. (2000). LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics 14, 6873.[CrossRef]
Helgason, E., Tourasse, N. J., Meisal, R., Caugant, D. A. & Kolsto, A. B. (2004). Multilocus sequence typing scheme for bacteria of the Bacillus cereus group. Appl Environ Microbiol 70, 191201.
Herias, M. V., Hessle, C., Telemo, E., Midtvedt, T., Hanson, L. A. & Wold, E. (1999). Immunomodulatory effects of Lactobacillus plantarum colonizing the intestine of gnotobiotic rats. Clin Exp Immunol 116, 283290.[CrossRef][Medline]
Huson, D. H. (1998). SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14, 6873.
Johansson, M.-L., Quednau, M., Molin, G. & Ahrné, S. (1995). Randomly amplified polymorphic DNA (RAPD) for rapid typing of Lactobacillus plantarum strains. Lett Appl Microbiol 21, 155159.[Medline]
Kleerebezem, M., Boekhorst, J., van Kranenburg, R. & 17 other authors (2003). Complete genome sequence of Lactobacillus plantarum WCFS1. Proc Natl Acad Sci U S A 100, 19901995.
Kumar, S., Tamura, K. & Nei, M. (1994). MEGA: molecular evolutionary genetics analysis software for microcomputers. Comput Appl Biosci 10, 189191.
Manfreda, G., de Cesare, A., Stella, S., Cozzi, M. & Cantoni, C. (2005). Occurrence and ribotypes of Listeria monocytogenes in Gorgonzola cheeses. Int J Food Microbiol 102, 287293.[CrossRef][Medline]
Marchesi, J. R., Sato, T., Weightman, A. J., Martin, T. A., Fry, J. C., Hiom, S. J. & Wade, W. G. (1998). Design and evaluation of useful bacterium-specific PCR primers that amplify genes coding for bacterial 16S rRNA. Appl Environ Microbiol 64, 795799.
Maynard Smith, J., Smith, N. H., O'Rourke, M. & Spratt, B. G. (1993). How clonal are bacteria? Proc Natl Acad Sci U S A 90, 43844388.
Merry, R. J., Dhanoa, M. S. & Theodorou, M. K. (1995). Use of freshly cultured lactic acid bacteria as silage inoculants. Grass For Sci 50, 112123.[CrossRef]
Molenaar, D., Bringel, F., Schuren, F. H., de Vos, W. M., Siezen, R. J. & Kleerebezem, M. (2005). Exploring Lactobacillus plantarum genome diversity by using microarrays. J Bacteriol 187, 61196127.
Nei, M. & Gojobori, T. (1986). Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3, 418426.[Abstract]
Oneca, M., Irigoyen, A., Ortigosa, M. & Torre, P. (2003). PCR and RAPD identification of L. plantarum strains isolated from ovine milk and cheese. Geographical distribution of strains. FEMS Microbiol Lett 227, 271277.[CrossRef][Medline]
Rodas, A. M., Ferrer, S. & Pardo, I. (2005). Polyphasic study of wine Lactobacillus strains: taxonomic implications. Int J Syst Evol Microbiol 55, 197207.
Ruiz-Barba, J. L., Cathcart, D. P., Warner, P. J. & Jiménez-Díaz, R. (1994). Use of Lactobacillus plantarum LPCO10, a bacteriocin producer, as a starter culture in spanish-style green olive fermentations. Appl Environ Microbiol 60, 20592064.
Sánchez, I., Seseña, S. & Palop, L. L. (2004). Polyphasic study of the genetic diversity of lactobacilli associated with Almagro eggplants spontaneous fermentation, based on combined numerical analysis of randomly amplified polymorphic DNA and pulse-field gel electrophoresis patterns. J Appl Microbiol 97, 446458.[CrossRef][Medline]
Urlings, H. A. P., Bijker, P. G. H. & van Logtestijn, J. G. (1993). Fermentation of raw poultry by-products for animal nutrition. J Anim Sci 71, 24202426.[Abstract]
Urwin, R. & Maiden, M. C. J. (2003). Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol 11, 479487.[CrossRef][Medline]
Vaquero, I., Marcobal, A. & Muñoz, R. (2004). Tannase activity by lactic acid bacteria isolated from grape must and wine. Int J Food Microbiol 96, 199204.[CrossRef][Medline]
Vescovo, M., Torriani, S., Dellaglio, F. & Bottazi, V. (1993). Basic characteristics, ecology and application of Lactobacillus plantarum: a review. Ann Microbiol Enzymol 43, 261284.
Yansanjav, A., Svec, P., Sedlacek, I., Hollerova, I. & Nemec, M. (2003). Ribotyping of lactobacilli isolated from spoiled beer. FEMS Microbiol Lett 229, 141144.[CrossRef][Medline]
Zadoks, R. N., Tikofsky, L. L. & Boor, K. J. (2005). Ribotyping of Streptococcus uberis from a dairy's environment, bovine feces and milk. Vet Microbiol 109, 257265.[CrossRef][Medline]
Zavaleta, A. I., Martínez-Murcia, A. J. & Rodríguez-Varela, F. (1996). 16S-23S rDNA intergenic sequences indicate that Leuconostoc oenos is phylogenetically homogeneous. Microbiology 142, 21052114.
Zhong, W., Millsap, K., Bialkowska-Hobrazanska, H. & Reid, G. (1998). Differentiation of Lactobacillus species by molecular typing. Appl Environ Microbiol 64, 24182423.
Received 31 August 2005;
revised 18 October 2005;
accepted 19 October 2005.
This article has been cited by other articles:
![]() |
E. Bilhere, P. M. Lucas, O. Claisse, and A. Lonvaud-Funel Multilocus Sequence Typing of Oenococcus oeni: Detection of Two Subpopulations Shaped by Intergenic Recombination Appl. Envir. Microbiol., March 1, 2009; 75(5): 1291 - 1300. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. L. W. Rademaker, H. Herbet, M. J. C. Starrenburg, S. M. Naser, D. Gevers, W. J. Kelly, J. Hugenholtz, J. Swings, and J. E. T. van Hylckama Vlieg Diversity Analysis of Dairy and Nondairy Lactococcus lactis Isolates, Using a Novel Multilocus Sequence Analysis Scheme and (GTG)5-PCR Fingerprinting Appl. Envir. Microbiol., November 15, 2007; 73(22): 7128 - 7137. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Diancourt, V. Passet, C. Chervaux, P. Garault, T. Smokvina, and S. Brisse Multilocus Sequence Typing of Lactobacillus casei Reveals a Clonal Population Structure with Low Levels of Homologous Recombination Appl. Envir. Microbiol., October 15, 2007; 73(20): 6601 - 6611. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Cai, B. T. Rodriguez, W. Zhang, J. R. Broadbent, and J. L. Steele Genotypic and phenotypic characterization of Lactobacillus casei strains isolated from different ecological niches suggests frequent recombination and niche specificity Microbiology, August 1, 2007; 153(8): 2655 - 2665. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |