Plesiomonas shigelloides is an emerging pathogen that is widespread in the aquatic environment and is responsible for intestinal diseases and extraintestinal infections in humans and other animals. Virtually nothing is known about its genetic diversity, population structure, and evolution, which severely limits epidemiological control. We addressed these questions by developing a multilocus sequence typing (MLST) system based on five genes (fusA, leuS, pyrG, recG, and rpoB) and analyzing 77 epidemiologically unrelated strains from several countries and several ecological sources. The phylogenetic position of P. shigelloides within family Enterobacteriaceae was precisely defined by phylogenetic analysis of the same gene portions in other family members. Within P. shigelloides, high levels of nucleotide diversity (average percentage of nucleotide differences between strains, 1.49%) and genotypic diversity (64 distinct sequence types; Simpson's index, 99.7%) were found, with no salient internal phylogenetic structure. We estimated that homologous recombination in housekeeping genes affects P. shigelloides alleles and nucleotides 7 and 77 times more frequently than mutation, respectively. These ratios are similar to those observed in the naturally transformable species Streptococcus pneumoniae with a high rate of recombination. In contrast, recombination within Salmonella enterica, Escherichia coli, and Yersinia enterocolitica was much less frequent. P. shigelloides thus stands out among members of the Enterobacteriaceae. Its high rate of recombination results in a lack of association between genomic background and O and H antigenic factors, as observed for the 51 serotypes found in our sample. Given its robustness and discriminatory power, we recommend MLST as a reference method for population biology studies and epidemiological tracking of P. shigelloides strains.

Phylogenetic relations within the family Enterobacteriaceae were analyzed using partial dnaJ sequences of 165 strains belonging to 93 species from 27 enterobacterial genera. The dnaJ phylogeny was in relative agreement with that constructed by 16S rDNA sequences, but more monophyletic groups were obtained from the dnaJ tree than from the 16S rDNA tree. The degree of divergence of the dnaJ gene was approximately 6 times greater than that of 16S rDNA. Also, the dnaJ gene showed the most discriminatory power in comparison with tuf and atpD genes, facilitating clear differentiation of any 2 enterobacterial species by dnaJ sequence analysis. The application of dnaJ sequences to the identification was confirmed by assigning 72 clinical isolates to the correct enterobacterial species. Our data indicate that analysis of the dnaJ gene sequences can be used as a powerful marker for phylogenetic study and identification at the species level of the family Enterobacteriaceae.

Real-time polymerase chain reaction and sequencing were used to characterize a 506-bp-long DNA fragment internal to the gyrB gene (gyrBint). The sequences obtained from 32 Enterobacteriaceae-type strains and those available in the Genbank nucleotide sequence database (n = 24) were used as a database to identify 240 clinical enterobacteria isolates. Sequence analysis of the gyrBint fragment of 240 strains showed that gyrBint constitutes a discriminative target sequence to differentiate between Enterobacteriaceae species. Comparison of these identifications with those obtained by phenotypic methods (Vitek 1 system and/or Rapid ID 32E; bioM?rieux, Marcy l'Etoile, France) revealed discrepancies essentially with genera Citrobacter and Enterobacter. Most of the strains identified as Enterobacter cloacae by phenotypic methods were identified as Enterobacter hormaechei strains by gyrBint sequencing. The direct sequencing of gyrBint would be useful as a complementary tool in the identification of clinical Enterobacteriaceae isolates.