MULTICOMP: a program for preparing sequence data for phylogenetic analysis

Population structure and evolution of non-O1/non-O139 Vibrio cholerae by multilocus sequence typing

Pathogenic non-O1/non-O139 Vibrio cholerae strains can cause sporadic outbreaks of cholera worldwide. In this study, multilocus sequence typing (MLST) of seven housekeeping genes was applied to 55 non-O1/non-O139 isolates from clinical and environmental sources. Data from five published O1 isolates and 17 genomes were also included, giving a total of 77 isolates available for analysis. There were 66 sequence types (STs), with the majority being unique, and only three clonal complexes. The V. cholerae strains can be divided into four subpopulations with evidence of recombination among the subpopulations. Subpopulations I and III contained predominantly clinical strains. PCR screening for virulence factors including Vibrio pathogenicity island (VPI), cholera toxin prophage (CTXΦ), type III secretion system (T3SS), and enterotoxin genes (rtxA and sto/stn) showed that combinations of these factors were present in the clinical isolates with 85.7% having rtxA, 51.4% T3SS, 31.4% VPI, 31.4% sto/stn (NAG-ST) and 11.4% CTXΦ. These factors were also present in environmental isolates but at a lower frequency. Five strains previously mis-identified as V. cholerae serogroups O114 to O117 were also analysed and formed a separate population with V. mimicus. The MLST scheme developed in this study provides a framework to identify sporadic cholera isolates by genetic identity.

Investigation of the enteric pathogenic potential of oral Campylobacter concisus strains isolated from patients with inflammatory bowel disease

Background

Campylobacter concisus, a bacterium colonizing the human oral cavity, has been shown to be associated with inflammatory bowel disease (IBD). This study investigated if patients with IBD are colonized with specific oral C. concisus strains that have potential to cause enteric diseases.

Methodology

Seventy oral and enteric C. concisus isolates obtained from eight patients with IBD and six controls were examined for housekeeping genes by multilocus sequence typing (MLST), Caco2 cell invasion by gentamicin-protection-assay, protein analysis by mass spectrometry and SDS-PAGE, and morphology by scanning electron microscopy. The whole genome sequenced C. concisus strain 13826 which was isolated from an individual with bloody diarrhea was included in MLST analysis.

Principal Findings

MLST analysis showed that 87.5% of individuals whose C. concisus belonged to Cluster I had inflammatory enteric diseases (six IBD and one with bloody diarrhea), which was significantly higher than that in the remaining individuals (28.6%) (P<0.05). Enteric invasive C. concisus (EICC) oral strain was detected in 50% of patients with IBD and none of the controls. All EICC strains were in Cluster 1. The C. concisus strain colonizing intestinal tissues of patient No. 1 was closely related to the oral C. concisus strain from patient No. 6 and had gene recombination with the patient’s own oral C. concisus. The oral and intestinal C. concisus strains of patient No. 3 were the same strain. Some individuals were colonized with multiple oral C. concisus strains that have undergone natural recombination.

Conclusions

This study provides the first evidence that patients with IBD are colonized with specific oral C. concisus strains, with some being EICC strains. C. concisus colonizing intestinal tissues of patients with IBD at least in some instances results from an endogenous colonization of the patient’s oral C. concisus and that C. concisus strains undergo natural recombination.

Population structure and evolution of pathogenicity of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is an enteric human pathogen but is widespread in the environment. Pathogenicity is determined by a number of virulence factors, including the virulence plasmid pYV, the high-pathogenicity island (HPI), and the Y. pseudotuberculosis-derived mitogen (YPM), a superantigen. The presence of the 3 virulence factors varies among Y. pseudotuberculosis isolates. We developed a multilocus sequence typing (MLST) scheme to address the population structure of Y. pseudotuberculosis and the evolution of its pathogenicity. The seven housekeeping genes selected for MLST were mdh, recA, sucA, fumC, aroC, pgi, and gyrB. An MLST analysis of 83 isolates of Y. pseudotuberculosis, representing 19 different serotypes and six different genetic groups, identified 61 sequence types (STs) and 12 clonal complexes. Out of 26 allelic changes that occurred in the 12 clonal complexes, 13 were mutational events while 13 were recombinational events, indicating that recombination and mutation contributed equally to the diversification of the clonal complexes. The isolates were separated into 2 distinctive clusters, A and B. Cluster A is the major cluster, with 53 STs (including Y. pestis strains), and is distributed worldwide, while cluster B is restricted to the Far East. The YPM gene is widely distributed on the phylogenetic tree, with ypmA in cluster A and ypmB in cluster B. pYV is present in cluster A only but is sporadically absent in some cluster A isolates. In contrast, an HPI is present only in a limited number of lineages and must be gained by lateral transfer. Three STs carry all 3 virulence factors and can be regarded as high-pathogenicity clones. Isolates from the same ST may not carry all 3 virulence factors, indicating frequent gain or loss of these factors. The differences in pathogenicity among Y. pseudotuberculosis strains are likely due to the variable presence and instability of the virulence factors.

The internal transcribed spacer region, a new tool for use in species differentiation and delineation of systematic relationships within the Campylobacter genus

The Campylobacter genus consists of a number of important human and animal pathogens. Although the 16S rRNA gene has been used extensively for detection and identification of Campylobacter species, there is currently limited information on the 23S rRNA gene and the internal transcribed spacer (ITS) region that lies between the 16S and 23S rRNA genes. We examined the potential of the 23S rRNA gene and the ITS region to be used in species differentiation and delineation of systematic relationships for 30 taxa within the Campylobacter genus. The ITS region produced the highest mean pairwise percentage difference (35.94%) compared to the 16S (5.34%) and 23S (7.29%) rRNA genes. The discriminatory power for each region was further validated using Simpson's index of diversity (D value). The D values were 0.968, 0.995, and 0.766 for the ITS region and the 23S and 16S rRNA genes, respectively. A closer examination of the ITS region revealed that Campylobacter concisus, Campylobacter showae, and Campylobacter fetus subsp. fetus harbored tRNA configurations not previously reported for other members of the Campylobacter genus. We also observed the presence of strain-dependent intervening sequences in the 23S rRNA genes. Neighbor-joining trees using the ITS region revealed that Campylobacter jejuni and Campylobacter coli strains clustered in subgroups, which was not observed in trees derived from the 16S or 23S rRNA gene. Of the three regions examined, the ITS region is by far the most cost-effective region for the differentiation and delineation of systematic relationships within the Campylobacter genus.

Population structure of Helicobacter pylori among ethnic groups in Malaysia: recent acquisition of the bacterium by the Malay population

BACKGROUND

Helicobacter pylori is a major gastric bacterial pathogen. This pathogen has been shown to follow the routes of human migration by their geographical origin and currently the global H. pylori population has been divided into six ancestral populations, three from Africa, two from Asia and one from Europe. Malaysia is made up of three major ethnic populations, Malay, Chinese and Indian, providing a good population for studying recent H. pylori migration and admixture.

RESULTS

Seventy eight H. pylori isolates, including 27 Chinese, 35 Indian and 16 Malay isolates from Malaysia were analysed by multilocus sequence typing (MLST) of seven housekeeping genes and compared with the global MLST data. STRUCTURE analysis assigned the isolates to previously identified H. pylori ancestral populations, hpEastAsia, hpAsia2 and hpEurope, and revealed a new subpopulation, hspIndia, within hpAsia2. Statistical analysis allowed us to identify population segregation sites that divide the H. pylori populations and the subpopulations. The majority of Malay isolates were found to be grouped together with Indian isolates.

CONCLUSION

The majority of the Malay and Indian H. pylori isolates share the same origin while the Malaysian Chinese H. pylori is distinctive. The Malay population, known to have a low infection rate of H. pylori, was likely to be initially H. pylori free and gained the pathogen only recently from cross infection from other populations.

Importation of the major pilin TcpA gene and frequent recombination drive the divergence of the Vibrio pathogenicity island in Vibrio cholerae

The Vibrio pathogenicity island (VPI) encodes the toxin-coregulated pilus and other virulence factors for Vibrio cholerae to colonize the human intestine to cause cholera. We assessed the level of genetic variation of VPI in nine nonpandemic isolates, and compared them with the sixth and seventh pandemic strains by sequencing c. 5 kb each from the start, middle and end regions of the VPI. Variation is similar among the three regions at around 2%, except for the tcpA gene, which has a much higher level of variation (23%). Numerous recombination segments were identified with sizes up to 2177 bp. Nearly all VPI genes sequenced have a ratio of synonymous to nonsynonymous substitutions considerably lower than that for housekeeping genes, suggesting that VPI genes are under positive selection pressure for change. The tagA gene was deleted or damaged in six isolates, which is likely to affect the efficiency of colonization of the human intestine. Two genes, orf2 and acfD, previously found to be translated differently in the sixth and seventh pandemic strains, were determined to be mutant in the seventh and sixth pandemic strains, respectively. These findings enhance our understanding of variation in the VPI, and of the pathogenic potential of VPI-positive environmental isolates.

Adaptation of multilocus sequencing for studying variation within a major clone: evolutionary relationships of Salmonella enterica serovar Typhimurium

Serovar Typhimurium of Salmonella enterica is a model organism for studies of pathogenesis that exhibits phage-type variation and variation in host range and virulence, but in a recent study showed no sequence variation in four genes, indicating the clonal nature of this serovar. We determined the relationships of 46 Typhimurium isolates of nine phage types using mutational changes detected either by matching AFLP (amplified fragment length polymorphism) fragments to computer-modeled LT2 AFLP fragments or by sequencing intergenic regions. Fifty-one polymorphic sites were detected, which gave a single phylogenetic tree. Comparison with genome sequences of five other serovars, Typhi, Paratyphi A, Gallinarum, Enteritidis, and Pullorum, enabled determination of the root of the tree. Only two parallel events were observed, giving high confidence in the tree branching order. The mutation-based tree provided a high level of consistency and a clear lineage for the Typhimurium isolates studied. This enabled us to show that for seven of the nine phage types used, the isolates studied have a single origin, but that two phage types clearly have more than one independent origin. We found that sequencing intergenic regions provides a good strategy for detection of mutational polymorphisms and study of phylogenetic relationships of closely related isolates and would be applicable to many other species.

Molecular basis of the indole-negative reaction in Shigella strains: extensive damages to the tna operon by insertion sequences

The molecular basis of the loss of tryptophan utilization (indole-negative phenotype) of Shigella strains, in effect clones of Escherichia coli, was investigated. Analysis of the tna operon of 23 Shigella strains representing each of the indole-negative serotypes revealed that insertion sequence-mediated insertion and/or deletions damaged the tna operon, leading to inability to convert tryptophan to indole. These events differ for cluster 1, cluster 3, and the outlier Shigella strains, confirming our previous observation of independent origins of these lineages from within E. coli. Parallel loss of the trait and prevalence of indole-negative strains suggest that the trait is deleterious in Shigella strains and advantages those without it.

Molecular evolutionary relationships of enteroinvasive Escherichia coli and Shigella spp

Enteroinvasive Escherichia coli (EIEC), a distinctive pathogenic form of E. coli causing dysentery, is similar in many properties to bacteria placed in the four species of Shigella. Shigella has been separated as a genus but in fact comprises several clones of E. coli. The evolutionary relationships of 32 EIEC strains of 12 serotypes have been determined by sequencing of four housekeeping genes and two plasmid genes which were used previously to determine the relationships of Shigella strains. The EIEC strains were grouped in four clusters with one outlier strain, indicating independent derivation of EIEC several times. Three of the four clusters contain more than one O antigen type. One EIEC strain (an O112ac:H- strain) was found in Shigella cluster 3 but is not identical to the Shigella cluster 3 D2 and B15 strains with the same O antigen. Two forms of the virulence plasmid pINV have been identified in Shigella strains by using the sequences of ipgD and mxiA genes, and all but two of our EIEC strains have pINV A. The EIEC strains were grouped in two subclusters with a very low level of variation, generally not intermingled with Shigella pINV A strains. The EIEC clusters based on housekeeping genes were reflected in the plasmid gene sequences, with some exceptions. Two strains were found in the pINV B form by using the ipgD sequence, with one strain having an mxiA sequence similar to the divergent sequence of D1. Clearly, EIEC and Shigella spp. form a pathovar of E. coli.

Comparison of two major forms of the Shigella virulence plasmid pINV: positive selection is a major force driving the divergence

All Shigella and enteroinvasive Escherichia coli (EIEC) strains carry a 230-kb virulence plasmid (pINV) which is essential for their invasiveness. There are two sequence forms, pINV A and pINV B, of the plasmid (R. Lan, B. Lumb, D. Ryan, and P. R. Reeves, Infect. Immun. 69:6303-6309, 2001), and the recently sequenced pINV plasmid from Shigella flexneri serotype 5 is a pINV B form. In this study we sequenced the majority of the coding region of the pINV A form from S. flexneri serotype 6 other than insertion sequence or related sequences and compared it with the pINV B form. More than half of the genes sequenced appear to be under positive selection based on their low ratio of synonymous to nonsynonymous substitutions. This high proportion of selected differences indicates that the two pINV forms have functional differences, and comparative studies of pathogenicity in different Shigella-EIEC strains could be informative. There are also genes absent in the S. flexneri serotype 6 plasmid, including the sepA gene encoding serine protease, the major secreted protein of S. flexneri serotype 2a, and the stbAB genes, which encode one of the two partition systems found in S. flexneri serotype 5. The incompatibility of the two pINV forms appears to be due to either small differences in the mvpAT postsegregational killing system or the presence of an unknown system in pINVA.

Extensive variation in the O-antigen gene cluster within one Salmonella enterica serogroup reveals an unexpected complex history

The 46 serogroups of Salmonella enterica have different O-antigens, and each is thought to have a specific form of the O-antigen cluster. Comparison of the 145 serovars of serogroup B revealed much more intraserogroup genetic diversity than expected. The O27 factor, due to an alpha 1-6 linkage between O units in place of the more common alpha 1-2 linkage and previously thought to be due to a converting bacteriophage, is now shown to be due to a wzy(alpha(1-6)) gene located within the major gene cluster. Surprisingly a remnant of this gene in all O27(-) serovars shows that the ancestor was O27(+). There are six distinct gene cluster forms, five apparently derived by a series of deletions and one by an insertion from an ancestral O27(+) form present in 57 serovars. The history of the gene cluster and movement between subspecies I and II can be traced. Two of the derivative forms still have a functional wzy(alpha(1-6)) gene, while in three it has been inactivated by deletion or insertion. Two of the forms lacking a functional wzy(alpha(1-6)) gene have the wzy(alpha(1-2)) gene first described for strain LT2 as rfc, whereas for the third the wzy gene has not been located.

Sequence analysis of four Shigella boydii O-antigen loci: implication for Escherichia coli and Shigella relationships

Shigella strains are in reality clones of Escherichia coli and are believed to have emerged relatively recently (G. M. Pupo, R. Lan, and P. R. Reeves, Proc. Natl. Acad. Sci. USA 97:10567-10572, 2000). There are 33 O-antigen forms in these Shigella clones, of which 12 are identical to O antigens of other E. coli strains. We sequenced O-antigen gene clusters from Shigella boydii serotypes 4, 5, 6, and 9 and also studied the O53- and O79-antigen gene clusters of E. coli, encoding O antigens identical to those of S. boydii serotype 4 and S. boydii serotype 5, respectively. In both cases the S. boydii and E. coli O-antigen gene clusters have the same genes and organization. The clusters of both S. boydii 6 and S. boydii 9 O antigens have atypical features, with a functional insertion sequence and a wzx gene located in the orientation opposite to that of all other genes in S. boydii serotype 9 and an rmlC gene located away from other rml genes in S. boydii serotype 6. Sequences of O-antigen gene clusters from another three Shigella clones have been published, and two of them also have abnormal structures, with either the entire cluster or one gene being located on a plasmid in Shigella sonnei or Shigella dysenteriae, respectively. It appears that a high proportion of clusters coding for O antigens specific to Shigella clones have atypical features, perhaps indicating recent formation of these gene clusters.

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

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.

Comparison of Vibrio cholerae pathogenicity islands in sixth and seventh pandemic strains

Epidemic Vibrio cholerae strains possess a large cluster of essential virulence genes on the chromosome called the Vibrio pathogenicity island (VPI). The VPI contains the tcp gene cluster encoding the type IV pilus toxin-coregulated pilus colonization factor which can act as the cholera toxin bacteriophage (CTXPhi) receptor. The VPI also contains genes that regulate virulence factor expression. We have fully sequenced and compared the VPI of the seventh-pandemic (El Tor biotype) strain N16961 and the sixth-pandemic (classical biotype) strain 395 and found that the N16961 VPI is 41,272 bp and encodes 29 predicted proteins, whereas the 395 VPI is 41,290 bp. In addition to various nucleotide and amino acid polymorphisms, there were several proteins whose predicted size differed greatly between the strains as a result of frameshift mutations. We hypothesize that these VPI sequence differences provide preliminary evidence to help explain the differences in virulence factor expression between epidemic strains (i.e., the biotypes) of V. cholerae.

Molecular characterization of Streptococcus pneumoniae type 4, 6B, 8, and 18C capsular polysaccharide gene clusters

Capsular polysaccharide (CPS) is a major virulence factor in Streptococcus pneumoniae. CPS gene clusters of S. pneumoniae types 4, 6B, 8, and 18C were sequenced and compared with those of CPS types 1, 2, 14, 19F, 19A, 23F, and 33F. All have the same four genes at the 5' end, encoding proteins thought to be involved in regulation and export. Sequences of these genes can be divided into two classes, and evidence of recombination between them was observed. Next is the gene encoding the transferase for the first step in the synthesis of CPS. The predicted amino acid sequences of these first sugar transferases have multiple transmembrane segments, a feature lacking in other transferases. Sugar pathway genes are located at the 3' end of the gene cluster. Comparison of the four dTDP-L-rhamnose pathway genes (rml genes) of CPS types 1, 2, 6B, 18C, 19F, 19A, and 23F shows that they have the same gene order and are highly conserved. There is a gradient in the nature of the variation of rml genes, the average pairwise difference for those close to the central region being higher than that for those close to the end of the gene cluster and, again, recombination sites can be observed in these genes. This is similar to the situation we observed for rml genes of O-antigen gene clusters of Salmonella enterica. Our data indicate that the conserved first four genes at the 5' ends and the relatively conserved rml genes at the 3' ends of the CPS gene clusters were sites for recombination events involved in forming new forms of CPS. We have also identified wzx and wzy genes for all sequenced CPS gene clusters by use of motifs.

Comparison of O-antigen gene clusters of Escherichia coli (Shigella) sonnei and Plesiomonas shigelloides O17: sonnei gained its current plasmid-borne O-antigen genes from P. shigelloides in a recent event

Escherichia coli Sonnei has an O antigen identical to that of Plesiomonas shigelloides O17, and its O-antigen gene cluster is located on a plasmid. By sequencing the chromosomal O-antigen gene cluster of P. shigelloides O17 and comparing it with that of Sonnei, we showed that Sonnei gained its O-antigen genes recently.

Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics

The evolutionary relationships of 46 Shigella strains representing each of the serotypes belonging to the four traditional Shigella species (subgroups), Dysenteriae, Flexneri, Boydii, and Sonnei, were determined by sequencing of eight housekeeping genes in four regions of the chromosome. Analysis revealed a very similar evolutionary pattern for each region. Three clusters of strains were identified, each including strains from different subgroups. Cluster 1 contains the majority of Boydii and Dysenteriae strains (B1-4, B6, B8, B10, B14, and B18; and D3-7, D9, and D11-13) plus Flexneri 6 and 6A. Cluster 2 contains seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and Dysenteriae 2. Cluster 3 contains one Boydii strain (B12) and the Flexneri serotypes 1-5 strains. Sonnei and three Dysenteriae strains (D1, D8, and D10) are outside of the three main clusters but, nonetheless, are clearly within Escherichia coli. Boydii 13 was found to be distantly related to E. coli. Shigella strains, like the other pathogenic forms of E. coli, do not have a single evolutionary origin, indicating convergent evolution of Shigella phenotypic properties. We estimate the three main Shigella clusters to have evolved within the last 35,000 to 270,000 years, suggesting that shigellosis was one of the early infectious diseases of humans.

Sequence diversity of the Escherichia coli H7 fliC genes: implication for a DNA-based typing scheme for E. coli O157:H7

Flagellar (H) antigens are mostly encoded by genes at the fliC locus in E. coli. We have sequenced 11 H7 fliC genes from Escherichia coli strains that belong to seven O serotypes. These sequences, together with those of nine other H7 fliC genes (from strains of three different O serotypes) sequenced recently (S. D. Reid, R. K. Selander, and T. S. Whittam, J. Bacteriol. 181:153-160, 1999), include 10 different sequences. The differences between these 10 sequences range from 0.06 to 3.12%. By comparison with other E. coli flagellin genes, we have identified primer length sequences specific for H7 genes in general and others specific for H7 genes of O157 and O55 strains: the specificity was confirmed by PCR testing the type strains for all 53 E. coli H types. We have previously identified genes specific for the E. coli O157 antigen, and use of the combination of O157- and H7-specific primers allows the sensitive and rapid detection of O157:H7 E. coli strains, which cause the majority of hemorrhagic colitis cases.

The wzz (cld) protein in Escherichia coli: amino acid sequence variation determines O-antigen chain length specificity

The O antigen is a polymer with a repeated unit. The chain length in most Escherichia coli strains has a modal value of 10 to 18 O units, but other strains have higher or lower modal values. wzz (cld/rol) mutants have a random chain length distribution, showing that the modal distribution is determined by the Wzz protein. Cloned wzz genes from E. coli strains with short (7 to 16), intermediate (10 to 18), and long (16 to 25) modal chain lengths were transferred to a model system, and their effects on O111 antigen were studied. The O111 chain length closely resembled that of the parent strains. We present data based on the construction of chimeric wzz genes and site-directed mutagenesis of the wzz gene to show that the modal value of O-antigen chain length of E. coli O1, O2, O7, and O157 strains can be changed by specific amino acid substitutions in wzz. It is concluded that the O-antigen chain length heterogeneity in E. coli strains is the result of amino acid sequence variation of the Wzz protein.

Evolutionary relationships among pathogenic and nonpathogenic Escherichia coli strains inferred from multilocus enzyme electrophoresis and mdh sequence studies

Within the species Escherichia coli, there are commensal strains and a variety of pathogenic strains, including enteropathogenic E. coli (EPEC), enterohemorrhagic E. coli (EHEC), enterotoxigenic E. coli (ETEC), enteroinvasive E. coli (EIEC), and urinary tract infection (UTI) strains. The pathogenic strains are identified by serotype and by possession of specific virulence determinants (toxins and adhesions, etc.) encoded by either monocistronic genes, plasmids, or pathogenicity islands. Although there are studies on the relationships between selected pathogenic strains, the relatedness among the majority of the pathogenic forms to each other, to commensal E. coli, and to the genus Shigella (which has often been suggested to be part of E. coli) has not been determined. We used multilocus enzyme electrophoresis (MLEE) at 10 enzyme loci and the sequence of the mdh housekeeping gene to study the genetic relationships of pathogenic E. coli strains (including Shigella clones), namely, 5 EPEC strains (serotypes O111 and O55), 3 EHEC strains (serotype O157), 6 ETEC strains (serotypes O78, O159, and O148), 5 EIEC strains (serotypes O124, O28, and O112), and 13 Shigella strains representing clones Flexneri, Dysenteriae, Boydii, and Sonnei, to commensal E. coli strains. Both the MLEE and mdh sequence trees reveal that EPEC, EHEC, ETEC, EIEC, and UTI strains are distributed among the ECOR set groups, with no overall clustering of EPEC, ETEC, EIEC, or UTI strains. The genus Shigella is shown to comprise a group of closely related pathogenic E. coli strains. Six pathogenic strains, i.e., M502 (EIEC; O112ac:NM), M503 (EPEC; O111:H12), M526 (ETEC; O159:H4), M522 (EPEC; O111ac:H12), M524 (ETEC; O78:H11), and M506 (ETEC; O78:H11), were found to have mdh sequences identical to those of five ECOR group A strains (ECOR5, ECOR10, ECOR14, ECOR6, and K-12). All 11 strains are closely related by MLEE. The results indicate that pathogenic strains of E. coli do not have a single evolutionary origin within E. coli but have arisen many times. The results also suggest the possibility that any E. coli strain acquiring the appropriate virulence factors may give rise to a pathogenic form.

The sixth and seventh cholera pandemics are due to independent clones separately derived from environmental, nontoxigenic, non-O1 Vibrio cholerae

The DNA sequences of the asd genes from 45 isolates of Vibrio cholerae (19 clinical O1 isolates, 2 environmental nontoxigenic O1 isolates, and 24 isolates with different non-O1 antigens) were determined. No differences were found within either sixth- or seventh-pandemic isolates; however, variation was found between the two forms and among the non-O1 isolates. O139 isolates had sequences identical to those of seventh-pandemic isolates. Phylogenetic trees with Vibrio mimicus as the outgroup suggest that the sixth-pandemic, seventh-pandemic, and U.S. Gulf isolates are three clones that have evolved independently from different lineages of environmental, nontoxigenic, non-O1 V. cholerae isolates. There is evidence for horizontal transfer of O antigen, since isolates with nearly identical asd sequences had different O antigens, and isolates with the O1 antigen did not cluster together but were found in different lineages. We also found evidence for recombination events within the asd gene of V. cholerae. V. cholerae may have a higher level of genetic exchange and a lower level of clonality than species such as Salmonella enterica and Escherichia coli.