Microbiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Microbiology 152 (2006), 85-93; DOI  10.1099/mic.0.28482-0
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.
Agricola
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.
Microbiology 152 (2006), 85-93; DOI  10.1099/mic.0.28482-0
© 2006 Society for General Microbiology

Development of a multilocus sequence typing method for analysis of Lactobacillus plantarum strains

Blanca de las Rivas, Ángela Marcobal and Rosario Muñoz

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
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Lactobacillus plantarum is a species of considerable industrial and medical interest. To date, the lack of reliable molecular methods for definite identification at strain level has hindered studies of the population biology of this organism. Here, a multilocus sequence typing (MLST) system for this organism is described, which exploits the genetic variation present in six housekeeping loci to determine the genetic relationship among isolates. The MLST system was established using 16 L. plantarum strains that were also characterized by ribotyping and restriction fragment length polymorphism (RFLP) analysis of the PCR-amplified 16S–23S rDNA intergenic spacer region (ISR). Ribotyping grouped the strains into four groups; however, RFLP analysis of the ISRs showed no differences in the strains analysed. In contrast, MLST had a good discriminatory ability. The sequence analysis of the six genes showed 14 different allelic combinations, with 12 of them represented by only one strain. By using this MLST approach we were able to confirm the identity of two strains deposited in the Spanish Type Culture Collection as different strains. Phylogenetic analysis indicated a panmictic population structure of L. plantarum and split decomposition analysis indicated that recombination plays a role in creating genetic heterogeneity in L. plantarum. As MLST allows precise identification, and easy comparison and exchange of results obtained in different laboratories, the future application of this new molecular method could be useful for the identification of valuable L. plantarum strains.


Abbreviations: dN, non-synonymous substitutions; dS, synonymous substitutions; ISA, standardized index of association; ISR, intergenic spacer region; MLST, multilocus sequence typing; RFLP, restriction fragment length polymorphism; ST, sequence typing

The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are AJ966402–AJ966404 (pgm), AJ966405–AJ966409 (ddl), AJ966364–AJ966370 (gyrB), AJ966371–AJ966378 (purK1), AJ966379–AJ966388 (gdh), AJ966389–AJ966396 (mutS) and AJ966397–AJ966401 (tkt4).


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Lactic acid bacteria (LAB) are often responsible for food preservation and flavour development. The specific environmental conditions prevailing in a fermenting food substrate promote the growth of some of these bacteria. Lactobacillus plantarum is predominantly found in fermented food and feed products, and is implicated in the processing of food for human consumption, such as sauerkraut, dry fermented sausage, wine and green olive fermentations and cheese making (Ruiz-Barba et al., 1994Down; Oneca et al., 2003Down), as well as in animal nutrition, such as crop preservation (Merry et al., 1995Down), fish and crab waste fermentation (Abazinge et al., 1993Down) and poultry by-product fermentation (Urlings et al., 1993Down).

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., 1999Down), 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., 1995Down; Oneca et al., 2003Down; Elegado et al., 2004Down) or variations of restriction enzyme analysis, for example ribotyping (Yansanjav et al., 2003Down; Rodas et al., 2005Down) and pulsed-field gel electrophoresis (Sánchez et al., 2004Down). 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., 2001Down), Vibrio cholerae (Farfán et al., 2002Down) and Bacillus cereus (Helgason et al., 2004Down); recently, MLST was also applied to the non-pathogenic food production bacterium Oenococcus oeni (De las Rivas et al., 2004Down).

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 16S–23S rDNA intergenic spacer region (ISR) and MLST in this species; and (ii) use MLST to analyse L. plantarum population structure.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Bacterial strains.
A total of 16 strains of L. plantarum were used in this study. The sources of the strains are listed in Table 1Down. Eight strains were provided by the Spanish Type Culture Collection (CECT). L. plantarum strains were routinely grown in MRS medium at 30 °C without shaking. Chromosomal DNA was prepared as described previously (Vaquero et al., 2004Down).


View this table:
[in this window]
[in a new window]
 
Table 1. Properties of L. plantarum isolates

RT, Ribotype; ITS, 16S–23S ISR type; ST, sequence type; ND, no data available.

 
Ribotyping.
Chromosomal DNA was digested with EcoRI and the products were separated by electrophoresis in 0·7 % agarose gels in 1x Tris/acetate/EDTA buffer. Digested DNA was transferred onto positively charged nylon membranes by Southern blotting. Probe 16S rDNA was obtained from L. plantarum CECT 748T by PCR by using the eubacterial universal primer pair 63f and 1387r (Marchesi et al., 1998Down). The 16S rDNA probe was digoxigenin-labelled and detected by chemiluminescence by using a DIG-High Prime DNA Labelling and Detection Starter Kit (Roche), according to the manufacturer's instructions.

RFLP of the PCR-amplified 16S–23S 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., 1996Down). 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)Down. The amplified 16S–23S 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 16S–23S 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 2Down). 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., 2003Down).


View this table:
[in this window]
[in a new window]
 
Table 2. Primers used for MLST

 
PCR was performed to amplify gene fragments from chromosomal DNA of the L. plantarum strains by using oligonucleotides listed in Table 2Up. The conditions of PCR, purification and sequencing of DNA fragments were described previously (De las Rivas et al., 2004Down).

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 1Up). 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, 1999Down) 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., 1994Down), 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, 1998Down). 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., 1993Down) 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, 1986Down). To avoid dependence on the number of loci, we also calculated the standardized IA ({85equ1}) (Haubold & Hudson, 2000Down).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Ribotyping and ISR RFLP analysis
Ribopattern analysis with EcoRI revealed four bands for all L. plantarum strains. The results are shown in Fig. 1Down(a). Four different groups of strains were defined on this basis: ribotype 1, containing L. plantarum CECT 220; ribotype 2, including L. plantarum type strain (CECT 748T); ribotype 3, represented only by L. plantarum RM40; and ribotype 4, containing L. plantarum CECT 4645 and RM35. The assignment of ribotype groups to each strain is listed in Table 1Up. Eleven out of 16 strains belonged to ribotype group 2.



View larger version (51K):
[in this window]
[in a new window]
 
Fig. 1. Typing analysis of the 16 L. plantarum strains examined in this study. (a) Ribotyping patterns of an EcoRI digest of L. plantarum chromosomal DNAs. Lanes: 1, ribotype 1; 2, ribotype 2; 3, ribotype 3; 4, ribotype 4. The sizes of the labelled fragments are indicated on the left. (b) 16S–23S rDNA ISR patterns obtained after AluI (lane 1), CfoI (2), DdeI (3) and TaqI (4) digestion of the 550 bp DNA fragment PCR-amplified with primers 16S14f and 23S1R. The sizes of some fragments in the 50 bp molecular mass marker are indicated on the left.

 
Primers 16S14F and 23S1R, complementary to target sequences at about 140 nt from the 3' end of the 16S rRNA gene and about 120 nt from the 5' end of the 23S rRNA gene (Zavaleta et al., 1996Down), respectively, were used to amplify the ISR of the L. plantarum strains. A PCR product of approximately 550 bp was obtained for each strain, indicating that the ISR is highly similar in all strains analysed (data not shown). PCR products of the 16S–23S rDNA ISR were digested with AluI, CfoI, DdeI and TaqI. All the strains analysed showed identical RFLP patterns for each enzyme (Fig. 1bUp).

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 3Down). Between 3 (pgm) and 10 (gdh) alleles were found per locus (Table 3Down, Fig. 2Down). The sequence differences in each allele are shown in Fig. 2Down. 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., 2003Down), the gdh allele sequences showed an unexpectedly high G+C content.


View this table:
[in this window]
[in a new window]
 
Table 3. Sequence variation at six loci

 


View larger version (51K):
[in this window]
[in a new window]
 
Fig. 2. Polymorphic nucleotide sites in L. plantarum MLST genes. Only the variable sites are shown. The nucleotide at each site is shown for a putative consensus sequence; only those that differ from the nucleotide in the consensus sequence are shown for the alleles. Nucleotide sites are numbered in vertical format from the first nucleotide position of the corresponding gene. The number of strains possessing the allele is indicated in parentheses.

 
The proportion of nucleotide alterations that changed the amino acid sequence (non-synonymous substitutions, dN) and the proportion of silent changes (synonymous substitutions, dS) were calculated for each gene. With these data, the dN/dS ratios were calculated for all loci and all were <1 (Table 3Up).

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. 2Up, Table 1Up) were dominant in the population.

A UPGMA dendrogram showing the genetic relatedness among L. plantarum strains investigated in this study is shown in Fig. 3Down. 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 1Up, Fig. 3Down).



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 3. UPGMA dendrogram showing the genetic relatedness of the L. plantarum strains examined in this study. Phylogenetic and molecular evolutionary analyses were conducted with the MEGA program (version 2.1) by the UPGMA method. Bootstrap confidence intervals are shown.

 
Relationship of L. plantarum strains by sequence analysis of housekeeping loci
For all genes, multiple strains carried identical alleles. However, only two pairs of strains had identical sequences for all fragments, and all other strains could be distinguished from each other because they had unique combinations of alleles (Table 1Up). Allelic profiles (ST) were assigned (Table 1Up). All STs differed in various loci, except ST-5 and ST-6 that differed only in the gdh locus. Concatenated pgm, ddl, gyrB, purK1, gdh and mutS gene sequence fragments were analysed and graphically displayed with SPLITSTREE (Fig. 4Down). The algorithm used in this software is able to display conflicting results in the phylogenetic descent of sequences. A tree-like structure is created when the descent is clonal, but an interconnecting network or a parallelogram will appear whenever recombination has been involved in the evolution of the analysed gene. Fig. 5Down shows the split graphs for all alleles of the six fragments analysed. The fit parameter was 100, indicating that all phylogenetic information in the sequences could be visualized in the graphs. We observed parallelograms for two of the six genes examined. The purK1 and gdh loci present network-like or parallelogram structures, indicating the presence of homoplasies, probably evolved by intergenic recombination. The split graphs obtained with pgm, ddl, gyrB and mutS loci showed no evidence of network-like evolution.



View larger version (8K):
[in this window]
[in a new window]
 
Fig. 4. Split decomposition analysis based on the allelic profilesof the 16 L. plantarum strains examined in this study. The numbering refers to strain numbers.

 


View larger version (10K):
[in this window]
[in a new window]
 
Fig. 5. Split decomposition analysis of alleles obtained from 16 L. plantarum strains for six loci. The observation that in the purK1 and gdh graphs several alleles in the sample are connected to each other by multiple pathways, forming parallelogram structures, is suggestive of recombination. All branch lengths are drawn to scale. The numbering refers to allele numbers.

 
Sawyer's test revealed no detectable cases of intragenic recombination in the sample, except in the case of the gdh locus in which there was evidence against the null hypothesis of no recombination (P=0·02) (Table 4Down). This evidence disappeared when the maximum condensed fragment value (P=1·00) was considered.


View this table:
[in this window]
[in a new window]
 
Table 4. Sawyer's test analysis for evidence of intragenic recombination

The Sawyer's test analysis was carried out by using the START program.Results were obtained from 10 000 random resamplings. SSCF, Sum of the square of condensed fragments; MCF, maximum condensed fragment.

 
Linkage disequilibrium between alleles was estimated with {85equ2}, a statistic that is expected to be zero when the alleles are in linkage equilibrium (free recombination). The latter probably means that the distribution of alleles occurs independently from each other. The {85equ3} value was 0·139, indicating that the genes investigated are close to linkage equilibrium. This low {85equ4} value is also indicative of extensive recombination.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The identification of strains belonging to the same species is still difficult. So far, a suitable and precise L. plantarum typing method is not available and it is urgent that one becomes available since strain characterization is necessary prior to patenting and release of a valuable L. plantarum strain for commercial applications.

Ribotyping is a both an inter- and an intraspecific typing method (Manfreda et al., 2005Down; Zadoks et al., 2005Down) that has been successfully used to differentiate between strains of the same Lactobacillus species. Yansanjav et al. (2003)Down 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., 1998Down). 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 16S–23S 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 1Up and Fig. 1Up, 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 3Up). 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 3Up). Our results corroborate the previously described genetic heterogeneity of this species (Vescovo et al., 1993Down; Sánchez et al., 2004Down; Molenaar et al., 2005Down). 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 1Up). 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, 2003Down). 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 1Up).

Using this typing method we could not find an association between ST and the source from which the different strains were isolated (Table 1Up). 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)Down 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. 4Up) 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., 2000Down). 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. 2Up). 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, 1998Down). 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., 2005Down). 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 3Up). 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 {85equ5} 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
 
This work has been supported by grants AGL2005-00470 and RM03-002. We thank R. González and E. García for their critical reading of the manuscript. The technical assistance of M. V. Santamaría is greatly appreciated. We also thank A. Hexter for correcting the English version. B. de las Rivas is a recipient of a postdoctoral fellowship and A. Marcobal of a predoctoral fellowship both from the Comunidad de Madrid.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Abazinge, M. D. A., Fontenot, J. P., Allen, V. G. & Flick, G. J. (1993). Ensiling characteristics of crab waste and wheat straw treated with different additives. J Agric Food Chem 41, 657–661.[CrossRef]

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, 7210–7219.[Abstract/Free Full Text]

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, 14–23.[Abstract/Free Full Text]

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, 11–18.[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, 3049–3060.[Abstract]

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, 1304–1313.[Abstract/Free Full Text]

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, 465–469.[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, 95–98.

Haubold, B. & Hudson, R. R. (2000). LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinformatics 14, 68–73.

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, 191–201.[Abstract/Free Full Text]

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, 283–290.[CrossRef][Medline]

Huson, D. H. (1998). SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14, 68–73.[Abstract/Free Full Text]

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, 155–159.[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, 1990–1995.[Abstract/Free Full Text]

Kumar, S., Tamura, K. & Nei, M. (1994). MEGA: molecular evolutionary genetics analysis software for microcomputers. Comput Appl Biosci 10, 189–191.[Abstract/Free Full Text]

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, 287–293.[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, 795–799.[Abstract/Free Full Text]

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, 4384–4388.[Abstract/Free Full Text]

Merry, R. J., Dhanoa, M. S. & Theodorou, M. K. (1995). Use of freshly cultured lactic acid bacteria as silage inoculants. Grass For Sci 50, 112–123.[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, 6119–6127.[Abstract/Free Full Text]

Nei, M. & Gojobori, T. (1986). Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3, 418–426.[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, 271–277.[CrossRef][Medline]

Rodas, A. M., Ferrer, S. & Pardo, I. (2005). Polyphasic study of wine Lactobacillus strains: taxonomic implications. Int J Syst Evol Microbiol 55, 197–207.[Abstract/Free Full Text]

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, 2059–2064.[Abstract/Free Full Text]

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, 446–458.[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, 2420–2426.[Abstract]

Urwin, R. & Maiden, M. C. J. (2003). Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol 11, 479–487.[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, 199–204.[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, 261–284.

Yansanjav, A., Svec, P., Sedlacek, I., Hollerova, I. & Nemec, M. (2003). Ribotyping of lactobacilli isolated from spoiled beer. FEMS Microbiol Lett 229, 141–144.[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, 257–265.[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, 2105–2114.[Abstract]

Zhong, W., Millsap, K., Bialkowska-Hobrazanska, H. & Reid, G. (1998). Differentiation of Lactobacillus species by molecular typing. Appl Environ Microbiol 64, 2418–2423.[Abstract/Free Full Text]

Received 31 August 2005; revised 18 October 2005; accepted 19 October 2005.


This article has been cited by other articles:


Home page
Appl. Environ. Microbiol.Home page
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]


Home page
Appl. Environ. Microbiol.Home page
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]


Home page
MicrobiologyHome page
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]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.
Agricola
Right arrow Articles by de las Rivas, B.
Right arrow Articles by Muñoz, R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2006 Society for General Microbiology.