Molecular evolution of large virulence plasmid in Shigella clones and enteroinvasive Escherichia coli

Abstract

Three genes, ipgD, mxiC, and mxiA, all in the invasion region of the Shigella virulence plasmid, were sequenced from strains representing a range of Shigella serotypes and from two enteroinvasive Escherichia coli (EIEC) isolates. The plasmids can be classified into two relatively homogeneous sequence forms which are quite distinct. pINV A plasmids are found in Shigella flexneri strains F6 and F6A, S. boydii strains B1, B4, B9, B10, B14, and B15, S. dysenteriae strains D3, D4, D6, D8, D9, D10, and D13, and the two EIEC strains (M519 and M520). pINV B plasmids are present in S. flexneri strains F1A, F2A, F3A, F3C, F4A, and FY, two S. boydii strains (B11 and B12), and S. sonnei. The D1 pINV plasmid is a recombinant with ipgD gene more closely related to those of pINV A but with mxiA and mxiC genes more closely related to those of pINV B. The phylogenetic relationships of the plasmid and those of the chromosomal genes of Shigella strains are largely consistent. The cluster 1 and cluster 3 strains tested (G.M. Pupo, R. Lan, and P. R. Reeves, Proc. Natl. Acad. Sci. USA 97:10567-10572, 2000) have pINV A and pINV B plasmids, respectively. However, of the three cluster 2 strains (B9, B11, and B15), B9 and B15 have pINV A while B11 has a pINV B plasmid. Those Shigella (D8 and D10 and S. sonnei) and EIEC strains which do not group with the main body of Shigella strains based on chromosomal genes were found to have plasmids belonging to one or the other of the two types and must have acquired these by lateral transfer.

FIG. 1

Alignment of the polymorphic sites for the three genes sequenced. Note that for some strains for which no sequence was obtained, the corresponding region appears blank. The numbers at the top of the figure, reading vertically, are base positions. Numbers (1, 2, and 3) at bottom are codon positions of that base (0 indicates a noncoding base), and an asterisk under a base indicates a phylogenetically informative site.

FIG. 2

Phylogenetic trees generated by the neighbor-joining method for the three genes sequenced. The number after the underscore gives the cluster number identified by chromosomal gene trees (26). Bootstrap values are percentages of 1,000 replications and are indicated at the nodes if the value is greater than 50%. Only strains with sequences for all three genes are included in the combined tree.