text:
Molecular Genomics
Molecular Evolution of a Pathogenicity Island from Enterohemorrhagic Escherichia coli O157:H7
Nicole T. Perna, George F. Mayhew, György Pósfai, Simon Elliott, Michael S. Donnenberg, James B. Kaper, Frederick R. Blattner
Nicole T. Perna
Laboratory of Genetics, University of Wisconsin-Madison, Madison Wisconsin 53706 1 ;
George F. Mayhew
György Pósfai
Institute of Biochemistry, Biological Research Center, H-6701 Szeged, Hungary 2 ; and
Simon Elliott
Center for Vaccine Development 3 and
Michael S. Donnenberg
Division of Infectious Disease, 4 University of Maryland School of Medicine, Baltimore, Maryland 21201
James B. Kaper
Frederick R. Blattner
We report the complete 43,359-bp sequence of the locus of enterocyte effacement (LEE) from EDL933, an enterohemorrhagicEscherichia coli O157:H7 serovar originally isolated from contaminated hamburger implicated in an outbreak of hemorrhagic colitis. The locus was isolated from the EDL933 chromosome with a homologous-recombination-driven targeting vector. Recent completion of the LEE sequence from enteropathogenic E. coli (EPEC) E2348/69 afforded the opportunity for a comparative analysis of the entire pathogenicity island. We have identified a total of 54 open reading frames in the EDL933 LEE. Of these, 13 fall within a putative P4 family prophage designated 933L. The prophage is not present in E2348/69 but is found in a closely related EPEC O55:H7 serovar and other O157:H7 isolates. The remaining 41 genes are shared by the two complete LEEs, and we describe the nature and extent of variation among the two strains for each gene. The rate of divergence is heterogeneous along the locus. Most genes show greater than 95% identity between the two strains, but other genes vary more than expected for clonal divergence among E. coli strains. Several of these highly divergent genes encode proteins that are known to be involved in interactions with the host cell. This pattern suggests recombinational divergence coupled with natural selection and has implications for our understanding of the interaction of both pathogens with their host, for the emergence of O157:H7, and for the evolutionary history of pathogens in general.
The locus of enterocyte effacement (LEE) is a 35-kb cluster of genes involved in the intimate adherence of pathogens to intestinal epithelial cells, the initiation of host signal transduction pathways, and the formation of attaching and effacing lesions (32, 33). Colony hybridization studies indicate that sequences homologous to the entire element are found in numerous enteropathogenic Escherichia coli (EPEC) and enterohemorrhagic E. coli (EHEC) strains and in other related bacteria (32). However, sequence data have been available for only a limited number of genes, and often from just one strain.
EPEC is an important cause of nonbloody infantile diarrhea in developing countries (8). EHEC O157:H7 causes both nonbloody diarrhea and a clinically distinct form of diarrheal disease known as hemorrhagic colitis that can lead to hemolytic uremic syndrome (17). These differences in clinical syndromes can be at least partly explained by the presence or absence of specific virulence factor genes. For example, the possession of bacteriophage-encoded Shiga toxins is a crucial distinction in the pathogenesis of disease due to these two pathogens, and acquisition of the phage-encoded toxins played a significant role in the evolution of EHEC from EPEC (48). Apart from the differences in virulence factors, EPEC and EHEC share the important intestinal histopathological phenotype known as attaching and effacing. This phenotype is distinguished by effacement of intestinal epithelial cell microvilli, intimate adherence of the bacteria to the epithelial cells, and marked changes in the host cell cytoskeleton. Both EPEC and EHEC have been shown to produce attaching and effacing lesions in tissue culture (26, 27) and animal models (15, 34, 45, 46). The genes responsible for this phenotype are contained in the LEE pathogenicity island, and sequence variation in one LEE gene has been implicated in different intestinal colonization sites (44, 50). How much of the similarities and differences in the pathogenesis of EPEC and EHEC disease can be attributed to sequence variation in the LEE is unknown.
The majority of information on genes contained in the LEE has been generated with EPEC strain E2348/69, and the complete LEE sequence for this strain has recently been reported (14). Theeae gene encodes a cell surface protein, intimin, involved with the intimate interaction of the pathogen and host epithelial cells (21). Three genes, espA (24),espB (12), and espD (29), encode secreted proteins involved with host signal transduction pathways. The esc (20) and sep(39) genes encode components of a type III secretory apparatus. The recently identified tir gene encodes a protein that is translocated from the bacterium to the host where it serves as a receptor for intimin (22).
Only three genes of the O157:H7 LEE have been previously sequenced,eae (3, 50), espB (13), and an open reading frame (ORF) of unknown function immediately upstream ofeae (orfU) (52). We now report the complete sequence of the LEE for EHEC O157:H7 strain EDL933. This information provides a rare opportunity for a comparative analysis of pathogenicity islands from two organisms that have important similarities as well as important differences in the pathogenesis of disease.
Bacterial strains.Eight strains of E. coli were used in this study. EDL933, an O157:H7 serotype obtained from the American Type Culture Collection (ATCC 43895), was originally isolated from contaminated hamburger implicated in an outbreak of hemorrhagic colitis (47). EPEC strain E2348/69 (O127:H6) and the five diarrheal E. coli (DEC) strains are from the University of Maryland Center for Vaccine Development collection. The DEC strains (48), kindly provided by Tom Whittam (The Pennsylvania State University) are C54-58 (O55:H6), F60-51 (O55:H7), 5625-50 (O55:H7), 3077-88 (O157:H7), and C374-83 (O157:H7). The K-12 strain MG1655, from the University of Wisconsin collection, was described by Guyer et al. (18).
LEE isolation and sequencing.Isolation of the LEE from EDL933 is discussed in detail elsewhere (37). In short, targeting vectors containing approximately 800 bp of known chromosomal sequence flanking the LEE were used to introduce a novel Flp recombinase target site on each side of the EDL933 LEE by homologous recombination. Expression of Flp from a helper plasmid promotes excision of the target as a plasmid. Random shotgun cloning into the Janus M13 vector (7) provided templates for automated sequencing (Applied Biosystems; model 377) with the ABI PRISM DyeTerminator Cycle Sequencing Ready Reaction Kit under standard conditions. Random clones were sequenced to provide an average of eightfold coverage of the entire element with a minimum quality of threefold coverage, including both strands, for all regions. The sequences were assembled with SeqManII software (DNASTAR) and edited manually. We note that the EDL933 sequence is presented in the conventional clockwise orientation of the K-12 genome. This representation is the reverse complement of that published for the EPEC E2384/69 LEE, which was oriented with respect to the transcription of the eae gene.
LEE gene identification and comparative analysis.ORFs were located with GeneQuest software (DNASTAR). Each ORF was annotated based on protein sequence searches (DeCypherII hardware/software system; TimeLogic, Inc.) against a combined SwissProt release 34 and trEMBLSP release 1 database. Each ORF in the sequence was assigned a unique identifier (L0001 to L0057), which appears in the GenBank submission. All sequence alignments of the EDL933 and E2348/69 LEEs were done with MegAlign or Align software (DNASTAR), and levels of divergence were assessed with the Molecular Evolutionary Genetic Analysis software package (28). The codon adaptation index was calculated by the method of Sharp and Li (41).
LEE phage PCR and sequencing.The selC-LEE junctions of EDL933, E2348/69, and five DEC strains were amplified by PCR. Genomic DNA was isolated in 2% Incert agarose (FMC BioProducts), as described by Kirkpatrick and Blattner (25). Agarose plugs were melted in equal volumes of 10 mM Tris-Cl (pH 8), and 1 μl was used as the template for a 50-μl PCR reaction mixture containing 2.5 U of TaKaRa LA Taq polymerase (PanVera), a 200 μM concentration of each dNTP, and a 2.5 μM concentration of each primer. The primers were designated leephage-f (located inselC) and leephage-r (in the nearest region of the LEE conserved between EDL933 and E2348/69) (Table1). These primers were chosen based on the EDL933 LEE sequence generated in this study. Strains with a LEE lacking a prophage generate a fragment of approximately 700 bp, while those with a full, intact sequence generate an approximately 8-kb fragment. An initial melting step of 94°C for 2 min was followed by 30 cycles of 94°C for 1 min and 68°C for 10 min, followed by a 72°C extension for 10 min. The PCR products were gel purified (Microcon/Micropure) according to the manufacturer's instructions and sequenced directly with the amplification primers.
Primers used in this study
Malate dehydrogenase (mdh) gene sequencing.Themdh genes from EDL933, E2384/69, and the DEC strains were amplified by PCR in a 50-μl PCR reaction mixture containing 1 μl of the genomic template described above, 2.5 U of TaKaRa Ex Taqpolymerase (PanVera), a 200 μM concentration of each dNTP, and a 2.5 μM concentration of each primer. The primers were designated mdh-a and mdh-c (Table 1). An initial melting step of 94°C for 1 min was followed by 25 cycles of 94°C for 15 s, 64°C for 15 s, and 72°C for 4 min. PCR products of the correct size (∼900 bp) were sequenced by primer walking. Primers were selected from the known MG1655 sequence for mdh (Table 1). Templates were sequenced under standard reaction conditions. Generated trace data files were then assembled and edited with SeqManII (DNASTAR). Consensus sequences were used to reconstruct an mdh phylogeny with the Kimura two-parameter distance matrix and the neighbor-joining algorithm of the Molecular Evolutionary Genetic Analysis software package (28).
Nucleotide sequence accession numbers.The sequences determined in this study have been assigned GenBank accession numbersAF071027 to AF071034.
Location and length of O157:H7 LEE.The region sequenced from EDL933 included chromosomal sequence flanking the LEE itself, allowing verification of the location by comparison with the K-12 sequence (4). A diagram of the EDL933 LEE and intact genes homologous to K-12 is shown in Fig. 1A. Near identity between the O157:H7 and K-12 chromosomes ends just beyond (16 bp) the 3′ end of the mature selC tRNA at bp 3833943 in the K-12 genome (GenBank accession no., U00096). In the K-12 chromosome,selC is followed by two hypothetical ORFs, yicKand yicL. In the EDL933 chromosome, the LEE is substituted for a segment (792 bp) of the K-12 chromosome including theselC-yicK intergenic region and 5′ end ofyicK. The substitution of the LEE is accompanied by at least one other local rearrangement relative to K-12, a 781-bp internal deletion in the yicK ORF (bp 3834812 to 3835592 in the K-12 genome). As inferred by a PCR assay, both junctions are conserved between EDL933 and EPEC strain E2348/69 (32).
Diagram of the EDL933 LEE. (A) ORFs are shown above and below the line to indicate the direction of transcription. Genes of the putative prophage are shown in black. Genes common to both the EDL933 and E2348/69 LEEs are shown in white. Genes homologous to those on the K-12 chromosome are hatched. The number of synonymous changes per synonymous site (dS) and (B) the number of nonsynonymous changes per nonsynonymous site (dN) (C) are shown for each ORF shared with the EPEC E2348/69 LEE (values are taken from Table 2). A scale (in base pairs) is shown along the bottom.
The EDL933 LEE is 43,359 bp (bp 903 to 44,262 of the sequence with GenBank accession no. AF071034), which is considerably greater than the 35,624-bp element sequenced from the EPEC strain (14). A 7.5-kb putative prophage near the selC end of the locus (Fig. 1A) accounts for most of the size difference. Although the 5′ ends of the two LEEs are nearly identical, they diverge abruptly after 118 bp. Excluding the 13 bp at the selC junction, this 105-bp region (bp 916 to 1020) is duplicated in the EDL933 LEE, 7,539 bp downstream (bp 8454 to 8558). An alignment of the two copies from EDL933 and the single copy from E2348/69 is shown in Fig.2. The sequence bounded by these repeats in EDL933 has no homolog in the E2348/69 LEE but contains several ORFs similar to ones found in retronphage phiR73 (42) and other members the CP4 family of cryptic prophage from K-12 previously described (4). The LEE prophage, with the direct repeats delimiting attR and attL, has been designated 933L.
Alignment of the two repeats flanking 933L, theselC end of the E2348/69 LEE, and the end of PAI-1 from uropathogenic E. coli (GenBank accession no., M13943). Residues that match the first copy of the repeat in EDL933 are represented by a dot. Coordinates shown to the left of the alignment correspond to the sequence position in each GenBank entry.
We have confirmed that the putative prophage is present in the genomic copy of the EDL933 LEE and absent from E2348/69 by PCR (Fig.3), with primers specific toselC and the nearest region of the locus shared by the two strains. The EDL933 PCR product is approximately 7.5 kb larger than the E2348/69 product. We also tested five additional strains, representing each major DEC clone (48) known to have its LEE adjacent toselC (49). Of these, two O157:H7 serovars exhibited PCR products identical in size to that of EDL933. Two others, both O55:H6 serovars, lack the LEE prophage and produce a band identical in size to that of E2348/69. One O55:H7 serovar yielded a product much larger than that of E2348/69 but roughly 1 kb smaller than that of EDL933 as a result of an internal deletion of a segment of the prophage. Single-read sequencing from the leephage-f and leephage-r primers was used to verify the identity of each product.
PCR to detect 933L in genomic DNA from MG1655, EDL933, 5624-50, and E2348/69 with primers leephage-f and leephage-r from Table1. One-half microliter of the reaction mixture was loaded on a 1% horizontal agarose gel in Tris-acetate-EDTA buffer and electrophoresed at 40 V for 1.5 h to visualize the product. The sizes of bands in Lambda DNA/HindIII markers (λHindIII) from Promega and a Low DNA MASS ladder (LDML) from GIBCO BRL in kilobase pairs are shown.
ORFs of the EDL933 LEE.We have identified a total of 54 ORFs in the extended EDL933 LEE (Fig. 1A). Of these, 13 fall within the putative prophage and 41 correspond to those previously described from the recently completed EPEC E2348/69 LEE sequence (14). The average codon adaptation index for the 54 ORFs, 0.219, is well below the average for the K-12 genome. As observed for EPEC E2348/69, the overall GC content (40.91%) is also below the K-12 average (50.80%). The prophage base composition is 51.72% G+C, while the remainder of the element is only 39.59% G+C.
The prophage ORF nearest to selC, L0003, encodes a 393-amino-acid (aa) protein with greater than 40% identity to integrases from a number of known bacteriophages including SF6 (SwissProt accession no., P37317), P4 (P08320), CP4-57 (P32053), and phiR73. L0004 encodes a 116-aa product 48% identical to a hypothetical Shigella dysenteriae IS911 protein (P39213) of similar size. A similar ORF is also found in CP4-6. L0005 encodes a short (60-aa) peptide with no significant matches in the combined SwissProt release 34 and trEMBLSP release 1 database or in the K-12 genome. The L0006 gene product resembles putative transposases from E. coli IS3 (P77673),Acinetobacter calcoaceticus (Q43916), and Erwinia amylovora (Q57113). The next two ORFs are similar to genes found only in the P4-like family of cryptic prophage from K-12 (4). The first of these two, L0007, encoding a 124-aa product, is shared by CP4-57, CP4-6, and CP4-44, all of which encode comparable-size peptides with >58% of the amino acids identical. L0008 has a match only in CP4-44, which encodes a considerably smaller product (163 versus 60 aa). The next four ORFs (L0009, L0010, L0011, and L0012) are complete unknowns. That is, there are no matches in the existing protein databases exceeding 30% identity alignable across at least 60% of both proteins. The remaining three ORFs within the boundaries of the putative prophage are most similar to hypothetical ORFs annotated in Agrobacterium and Rhizobiumplasmid sequences. The products of all three are similar in length to the database entries (133 aa versus 154 aa for Q52592, 115 aa versus 115 aa for P50359, and 512 aa versus 511 aa for P55504). The first gene in this group, L0013, is also similar to the 5′ end of a hypothetical insertion element IS2 gene roughly three times its size. The products of the latter two, L0014 and L0015, match proteins that are part of larger families of paralogous hypothetical proteins encoded by a Rhizobium plasmid.
The identities of the 41 genes common to both the EHEC EDL933 and EPEC 2348/69 LEEs have been presented in some detail elsewhere (14). Absolute conservation of gene order and number is observed between these elements. Table 2lists the corresponding ORFs from the EPEC strain and a brief description of the known or putative function. Further comment on individual ORFs will be made in the context of exploring the sequence variation among these strains. Additional small ORFs can be found in both the EHEC and EPEC LEEs but are not included in our final determinations.
ORFs shared by the EPEC and EHEC LEEs
Sequence comparison with EPEC.Each O157:H7 LEE gene and protein were aligned to their homologs from E2348/69. Protein alignments were used to refine nucleotide sequence alignments by eye. Summary data for each comparison are shown in Table 2. The average level of nucleotide identity is 93.9%. The number of gaps and total gap length in each alignment are shown. We calculated rates of synonymous and nonsynonymous substitution for each gene. These estimates are plotted in Fig. 1B and C to illustrate the distribution of variation across the length of the LEE element. It is clear from both Table 2 and Fig. 1 that there is considerable variation among genes in the level of sequence divergence. In general, variable ORFs are more variable in every way. They have more gaps and higher rates of both synonymous and nonsynonymous substitution.
One group of highly variable genes includes espB (25.99% difference), espD (19.64% difference), and espA(15.37% difference). All three encode secreted proteins implicated in the activation of host epithelial signal transduction, intimate adherence, and formation of attaching and effacing lesions (23,24, 29). espB, espD, and espAare adjacent in a tightly packed, presumably cotranscribed, cluster that also encodes a potential chaperone, one putative component of the secretion apparatus, and two proteins of unknown function (9). The first gene of this group, one of the unknowns, shows moderate sequence divergence (6.60%) in addition to a variable-number repeat structure described in more detail below. The putative secretion apparatus gene in this region, escF, is invariant between EDL933 and E2348/69. A gene of unknown function, L0023, is found on the other side of the variable cluster directly adjacent to espA. It varies at 5.68% of the nucleotide sites but has an absolutely conserved length in these two strains. L0023 is followed by a highly conserved secretion system gene,escD, contained on the opposite strand.
Another gene that shows marked differences is eae, which encodes intimin, perhaps the best characterized of the LEE proteins. This outer membrane protein is required (but not sufficient) for the intimate adherence to epithelial cells characteristic of attaching and effacing enteropathogens and for full virulence of EPEC in human volunteers (10, 11). Cell binding activity is known to reside within the C-terminal 192 aa of the intimin protein (16). Yu and Kaper (50) sequenced eaefrom both EDL933 and E2348/69 and suggested that the C-terminal variability between the corresponding intimins might indicate that the strains bind different eukaryotic receptors. Since then, eaesequences from numerous other strains have been reported, and comparative analyses of these have shown divergence patterns similar to that observed in the initial EDL933-E2348/69 comparison and have also shown that there are at least three groupings of eae(2, 16). In our study, the two eae genes were found to be 87.23% identical, or slightly less variable than the secreted proteins described above.
A second group of contiguous genes is highly variable. The most variable of this group is tir (33.52% difference), which was recently discovered to encode a product that is translocated from the bacterium to the host cell where it most likely serves as the intimin receptor (22). For the other three genes, L0028 (17.48% difference), L0029 (21.94% difference), and L0030 (25.30% difference), very little is known about their role in pathogenesis. The L0028 product resembles a hypothetical protein encoded in aShigella virulence-associated cluster (14). Without additional functional characterization, it is impossible to evaluate the adaptive significance of variation at these loci. Among the remaining LEE genes, only one, sepZ, exhibits elevated divergence (29.29%) between EDL933 and E2348/69. Interestingly, asepZ mutant of EPEC exhibits certain phenotypic similarities with EHEC, including reduced invasion efficiency and lack of tyrosine phosphorylation of Hp90, yet retains the ability to form attaching and effacing lesions (39).
Alignments of one gene, L0016, revealed a particularly interesting feature. The E2348/69 sequence requires a single 126-bp gap to accommodate the obviously homologous sequences at the 5′ and 3′ ends. The EDL933 sequence in this gap encodes a third copy of a proline-rich repeat noted by Donnenberg et al. (9) in the E2348/69 protein. Figure 4 is an alignment of the three complete copies and one partial copy of the repeat from EDL933 with the two complete copies and one partial copy from E2348/69. The first copies of the repeat in both strains, R1EDL933 and R1E2348/69, have several bases in common (positions 5, 6, 10, and 12) near the beginning of the alignment that are not conserved among the other complete copies, R2EDL933, R3EDL933, and R2E2348/69. Similarly, the two partial repeats near the ends of the genes, R4EDL933 and R3E2348/69, are most like each other in the central region of the alignment (positions 24, 55 to 67, 76, and 114). At other sites (positions 96, 99, 102, and 125 to 141) the copies within each strain appear to be more similar to each other than to the copies from the other strain. This suggests either gene conversion to homogenize the repeats within a sequence or separate expansion of repeats in each lineage by either unequal recombination or slipped-strand mispairing. Several other database proteins, mainly eukaryotic, exhibit similar repetitive structures, although there is no functional consensus.
Alignment of the repetitive motifs of EDL933 gene L0016 and the E2348/69 homolog. Coordinates shown to the left of the alignment correspond to the sequence position in each GenBank entry. Each sequence is labeled with the repeat number and subscripted strain designation (R1EDL933, R2EDL933, etc.). The repeat is in frame in both strains, and the putative peptide sequence is shown above the alignment with asterisks replacing residues that are not absolutely conserved. Dots indicate nucleotides identical to those in the R1EDL933 sequence. Dashes have been added to maximize the similarity among sequences.
Many of the intergenic regions are quite small, but several additional noncoding regions merit comparison. The large intergenic region between L0054 and L0055 in the EDL933 LEE and that between rorf2 andorf1 of the E2348/69 LEE include a member of the ERIC family of repeated elements and are 98.4% identical in the two strains. Previously, similarity between the selC end of the EPEC LEE and PAI-1 of uropathogenic E. coli has been noted (32). This corresponds to the repeated region flanking the prophage in EDL933. The EPEC copy is more like the second repeat in EHEC (97 versus 83% identical). All three LEE repeats are equally similar to PAI-1 (73% identical). An alignment of these sequences is shown in Fig. 2.
Comparisons with other published espA andespB genes.Our EDL933 espB sequence differs from a published one (13) by one synonymous change. Sequence data were available for espA and espBfrom a rabbit EPEC strain, RDEC-1 (1). The espAsequences of E2348/69 and RDEC-1 are more similar to each other (91.0% similarity) than either is to that of EDL933 (84.6 and 87.6% similarity, respectively). The RDEC-1 espB is more similar to both that of EDL933 (76.8% similarity) and E2348/69 (73.4% similarity) than they are to each other (68.6% similarity). In sharp contrast, the espB sequences of RDEC-1 and bovine EHEC O26 serovar 413/89-1 (13) are nearly identical (two synonymous differences in 945 bp) (1).
Phylogenetic relationship among strains.To interpet the observed variation between the two complete LEEs and the distribution of the phage, it is helpful to understand the relationship between the strains under consideration. The gene encoding malate dehydrogenase (mdh) has been used to describe the relationships amongE. coli strains (6, 38). We have sequenced PCR-amplified mdh from EDL933, E2348/69, and the DEC strains used in this study. None of the E. coli mdh sequences differ from another by more than 3% identity. These sequences were used to reconstruct an updated phylogeny (Fig.5). This tree is consistent with the relationships observed in a multilocus enzyme electrophoresis-based phylogeny of EPEC and EHEC strains (49) but additionally suggests that EDL933 and E2348/69 are both more closely related to strains from the nondiarrheagenic ECOR collection than to each other. As predicted by multilocus enzyme electrophoresis (48, 49), the O55:H7 serovar and all three O157:H7 strains, including EDL933, are very closely related. The distribution of the LEE prophage is indicated on the tree (Fig. 5).
Evolutionary relationships among strains used in this study and 20 E. coli reference strains based onmdh sequences. A phylogeny was reconstructed with published data from the ECOR collection (6, 38) (GenBank accession no., U04742 through U04758, U04770, and AF004201) and our EDL933, E2348/69, and DEC strain mdh sequences. The tree is based on a MegAlign (DNASTAR) multiple sequence alignment of positions 34 through 878 of the 916-bp mdh coding region, excluding positions 652 and 653. The phylogeny was reconstructed with a Kimura two-parameter distance matrix and the neighbor-joining algorithim. The percentages of 2,000 bootstrap replicates supporting each cluster are shown along the branches.
We have described the 43.359-kb LEE from EHEC EDL933 and compared it to the entire homologous locus from EPEC E2348/69. The EDL933 LEE has precisely the same chromosomal integration point as the LEE of EPEC E2348/69 and is highly similar for most genes. However, the O157:H7 LEE differs structurally from the 35.5-kb EPEC locus by the inclusion of P4 family cryptic prophage 933L, and the two LEEs exhibit high levels of divergence for some genes.
The base composition and codon usage of the LEE are atypical ofE. coli, suggesting that the element was horizontally transferred into the species. One very basic question is whether the LEE was acquired once or repeatedly. As predicted by McDaniel et al. (32), EDL933 and E2348/69 LEEs are found in homologous chromosomal positions and have very similar sequences at both junctions of the insertion point relative to the K-12 chromosome. An examination of the mdh phylogeny shows that EDL933 and E2348/69 are relatively distantly related E. coli, both of which are more closely related to nondiarrheagenic ECOR strains than to each other. In a colony hybridization survey, no nonpathogenic E. coliexamined exhibited a LEE element (32). If the LEE was acquired only once and is not found among commensal strains, it must have been subsequently lost independently by numerous lineages. Alternatively, the LEEs shared by EDL933, E2348/69, and a variety of other EPEC strains may be the result of lateral exchange. Under this model, LEE is introduced again and again into different clonal lineages by recombination either with a related species or with another strain of E. coli that already harbors the element. Wieler et al. (49) have recently shown that there are both EPEC and EHEC strains that contain the LEE but that have an intact selC, indicating that there is at least one other chromosomal location for the LEE among E. coli strains.
The only previously presented evidence suggesting autonomous mobility of the LEE is the slight similarity of the E2348/69 LEE ends and the IS600 family of insertion elements, including a small ORF encoding a peptide with similarity to part of the putative transposase (9). These regions are highly conserved between the EDL933 and E2348/69 sequences. If both LEEs were independently mobilized by a precursor of these transposase remnants, we would not expect them to differ from a functional transposon in such similar ways. Although we do not rule out the possibility that the LEE originally entered anE. coli selC locus by this mechanism, the conserved remnants suggest that the EDL933 and E2348/69 LEEs have a common ancestor with an already-defective transposase. One possible scenario is transposition of the LEE into one E. coli selC, loss of mobility, and subsequent lateral transfer distributing the LEE at the selC locus among E. coli lineages.
The cryptic P4-like prophage associated with the EDL933 LEE appears to integrate into the very end of the LEE itself rather than theselC locus. This is notable since selC acts as a recombinational hot spot, serving as the site of integration for the related retronphage phiR73, other phage, and another pathogenicity island, PAI-1, in uropathogenic E. coli. PCR and sequencing suggest that this phage integrated into the LEE prior to the divergence of the O55:H7 and O157:H7 strains shown in the mdh tree. Given the ability of prophage to be excised, leaving an intact integration site, it is not possible to tell whether there was ever a phage in the E2348/69 LEE. Regardless, it seems unlikely that the LEE was ever mobilized by this prophage since both putative attsites occur on one side of the LEE coding region.
Other than the putative prophage, gene content and organization are highly conserved between the EDL933 and E2348/69 LEEs. There is however, strikingly large and nonrandomly distributed sequence variation within the cluster (Fig. 2). Homologous E. coliloci derived from direct clonal descent generally differ at less than 5% of the nucleotide sites. mdh, with an average pairwise divergence of 2.1%, is a standard example (36). Many of the LEE genes exhibit comparable levels of sequence divergence, but eight genes vary at more than 15% of the alignable nucleotide sites.
Subdivision of nucleotide variation into synonymous and nonsynonymous changes sheds some additional light on the mechanisms of LEE evolution. Synonymous base substitutions are those changes that do not affect the amino acid sequence of the encoded protein. Nonsynonymous substitutions are reflected in the protein and are subject to natural selection. Several observations can be made about the pattern of divergence for the two complete LEEs. Synonymous substitutions do not occur at a constant rate across the element. High rates of nonsynonymous change are always accompanied by elevated synonymous divergence and in no case does the estimated nonsynonymous rate exceed the synonymous rate, indicating purifying selection on all proteins encoded in the locus.
In the absence of selection constraining the accumulation of synonymous differences, the rate should be constant across sequences that have been evolving independently as intact units. Instead, synonymous substitution rates for the 41 shared LEE proteins are quite heterogeneous. For many of the less variable genes, the estimated synonymous substitution rate is very close to zero indicating very recent divergence of at least some regions of the EDL933 and E2348/69 LEEs. All genes previously described as variable overall have notably higher synonymous substitution rates, as do L0023, L0026,sepQ, and L0038, all of which are adjacent to a hypervariable gene. It is difficult to envision a selection-based model of evolution that could lead to the observed disparity in the frequency of synonymous changes in the LEEs if the mutation rate is constant along the entire locus. The patchy distribution of synonymous variation is most easily explained by recombinational events that unite genes or clusters of genes with distinct mutational histories. The number of recombinations and donor sources involved may be discernable by comparisons of more LEE sequences from E. coli strains and related species.
Nonsynonymous changes and natural selection of the encoded proteins most likely determines the fate of new alleles. The most striking feature of the divergent genes is that they include all of those encoding proteins known to interact directly with the host:eae, espA, espB, espD, andtir. Intimin binds to the host cells. Intimin also appears to influence the site of intestinal colonization. Complementation of an O157:H7 eae mutant with plasmid-borne EPEC eaeresults in colonization of the distal half of the small intestine and the surface of the large intestine in gnotobiotic pigs (44). This is more typical of EPEC than EHEC, which normally colonizes only the lower bowel. Kenny and Finlay have observed the espA-,espB-, and espD-encoded proteins associated with HeLa cells after EPEC infection (23). The tirgene product is translocated into the host cell (22). In EPEC, the tir gene product is tyrosine phosphorylated, becoming the protein originally described as Hp90 (40), and is presented on the host cell surface where it serves as the intimin receptor. EHEC O157:H7 fails to induce tyrosine phosphorylation of the receptor protein in HEp-2 and T84 cells (19). Together, these observations suggest that the variability we observe between the EPEC and EHEC LEEs has phenotypic effects that could be subject to natural selection for adaptation to either host specificity or evasion of the host immune system. In contrast to the hypervariability seen in genes encoding proteins known to interact directly with the host, other LEE-encoded proteins that do not interact with the host are highly conserved between both EHEC and EPEC LEEs, notably theesc genes encoding the type III secretion system.
Hypervariability of genes directly involved with host interaction is emerging as a paradigm of molecular evolution in pathogens. Two studies (5, 30) have noted a relationship between cellular location of the gene product and rate of divergence of the inv-spainvasion gene complex in Salmonella enterica. Three adjacent genes in this cluster exhibit elevated nonsynonymous site variation relative to other genes at the locus, among which are components of a type III secretion apparatus similar to that of the LEE. Two of the three, spaO and spaN, encode secreted proteins, while the product of the third, spaM, remains unknown. Boyd et al. (5) ruled out selection for host specificity as a mechanism driving the observed variability because of the lack of polymorphism among strongly host-adapted serovars but suggested that selection for antigenic diversity may be responsible. High levels of polymorphism have been observed at loci involved with host-pathogen interactions in many different bacterial species, includingBorrelia burgdorferi (43), Streptococcus pyogenes (31), Neisseria spp. (35), and Mycoplasma hominis (51). Future genomic sequencing studies of these and other pathogens will likely yield additional insights into those proteins that are important in host-pathogen interactions.
This work was supported by NSF-Sloan Foundation Molecular Evolution Postdoctoral Fellowship BIR-9626042 (N.T.P.); NIH grants AI41329-01 (F.R.B.), AI32074 (M.S.D.), and AI21637 and AI41325 (J.B.K.); HHMI grant 75195-542102 (G.P.); and Ronald McDonald House Charities.
We thank Guy Plunkett, Val Burland, and Jeremy Glasner for advice and thoughtful consideration of both the data and manuscript. We are grateful for the expert technical assistance provided by Heather Kirkpatrick, Jason Gregor, Guy Peyrot, Pritin Soni, Mike Goeden, and Wayne Davis.
Received 29 December 1997.
Returned for modification 6 March 1998.
Accepted 27 May 1998.
Copyright © 1998 American Society for Microbiology
Abe A.,
Kenny B.,
Stein M.,
Finlay B. B.
(1997) Characterization of two virulence proteins secreted by rabbit enteropathogenic Escherichia coli, EspA and EspB, whose maximal expression is sensitive to host body temperature. Infect. Immun. 65:3547–3555.
Agin T. S.,
Wolf M. K.
(1997) Identification of a family of intimins common to Escherichia coli causing attaching-effacing lesions in rabbits, humans, and swine. Infect. Immun. 65:320–326.
Beebakhee G.,
Louie M.,
Azavedo J. D.,
Brunton J.
(1992) Cloning and nucleotide sequence of the eae gene homologue from enterohemorrhagic Escherichia coli serotype O157:H7. FEMS Microbiol. Lett. 91:63–68.
Blattner F. R.,
Plunkett III G.,
Bloch C. A.,
Perna N. T.,
Burland V.,
Riley M.,
Collado-Vides J.,
Glasner J. D.,
Rode C. K.,
Mayhew G.,
Gregor J.,
Davis N. W.,
Kirkpatrick H. A.,
Goeden M.,
Rose D.,
Mau B.,
Shao Y.
(1997) The complete sequence of Escherichia coli K-12. Science 277:1453–1462.
Boyd E. F.,
Li J.,
Ochman H.,
Selander R. K.
(1997) Comparative genetics of the inv-spa invasion gene complex of Salmonella enterica. J. Bacteriol. 179:1985–1991.
Nelson K.,
Wang F.-S.,
Whittam T. S.,
(1994) Molecular genetic basis of allelic polymorphism in malate dehydrogenase (mdh) in natural populations of Escherichia coli and Salmonella enterica. Proc. Natl. Acad. Sci. USA 91:1280–1284.
Daniels D. L.,
Blattner F. R.
(1993) Genome sequencing on both strands: the Janus strategy. Nucleic Acids Res. 21:3385–3390.
Donnenberg M. S.
(1995) Enteropathogenic Escherichia coli. in Infections of the gastroinstestinal tract. eds Blaser M. J., Smith P. D., Ravdin J. I., Greenberg H. B., Guerrant R. L. (Raven Press, New York, N.Y), pp 709–726.
Donnenberg M. S.,
Lai L.-C.,
Taylor K. A.
(1997) The locus of enterocyte effacement pathogenicity island of enteropathogenic E. coli encodes secretion functions and remnants of transposons at its extreme right end. Gene 184:107–114.
Tacket C. O.,
James S. P.,
Losonsky G.,
Nataro J. P.,
Wasserman S. S.,
Kaper J. B.,
Levine M. M.
(1993) Role of the eaeA gene in experimental enteropathogenic Escherichia coli infection. J. Clin. Invest. 92:1412–1417.
Tzipori S.,
McKee M. L.,
O'Brien A. D.,
Alroy J.,
Kaper J. B.
(1993) The role of the eae gene of enterohemorrhagic Escherichia coli in intimate attachment in vitro and in a porcine model. J. Clin. Invest. 92:1418–1424.
Yu J.,
(1993) A second chromosomal gene necessary for intimate attachment of enteropathogenic Escherichia coli to epithelial cells. J. Bacteriol. 175:4670–4680.
Ebel F.,
Deibel C.,
Kresse A. U.,
Guzman C. A.,
Chakraborty T.
(1996) Temperature- and medium-dependent secretion of proteins by Shiga toxin-producing Escherichia coli. Infect. Immun. 64:4472–4779.
Elliott S.,
Wainwright L. A.,
McDaniel T.,
MacNamara B.,
Donnenberg M.,
(1998) The complete sequence of the locus of enterocyte effacement (LEE) from enteropathogenic E. coli E2348/69. Mol. Microbiol. 28:1–4.
Francis D. H.,
Collins J. E.,
Duimstra J. R.
(1986) Infection of gnotobiotic pigs with an Escherichia coli O157:H7 strain associated with an outbreak of hemorrhagic colitis. Infect. Immun. 51:953–956.
Frankel G.,
Candy D. C. A.,
Fabiani E.,
Adu-Bobie J.,
Gil S.,
Novakova M.,
Phillips A. D.,
Dougan G.
(1995) Molecular characterization of a carboxy-terminal eukaryotic-cell-binding domain of intimin from enteropathogenic Escherichia coli. Infect. Immun. 63:4323–4328.
Griffin P. M.,
Tauxe R. V.
(1991) The epidemiology of infections caused by Escherichia coli O157:H7, and other enterohemorrhagic E. coli, and the associated hemolytic uremic syndrome. Epidemiol. Rev. 13:60–98.
Guyer M.,
Reed R. E.,
Steitz T.,
Low K. B.
(1981) Identification of a sex-factor-affinity site in E. coli as γδ. Cold Spring Harbor Symp. Quant. Biol. 45:135–140.
Ismaili A.,
Philpott D. J.,
Dytoc M. T.,
Sherman P. M.
(1995) Signal transduction responses following adhesion of verotoxin-producing Escherichia coli. Infect. Immun. 63:3316–3326.
Jarvis K. G.,
Girón J. A.,
Jerse A. E.,
McDaniel T. K.,
(1995) Enteropathogenic Escherichia coli contains a putative type III secretion system necessary for the export of proteins involved in attaching and effacing lesion formation. Proc. Natl. Acad. Sci. USA 92:7996–8000.
Tall B. D.,
(1990) A genetic locus of enteropathogenic Escherichia coli necessary for the production of attaching and effacing lesions on tissue culture cells. Proc. Natl. Acad. Sci. USA 87:7839–7843.
DeVinney R.,
Reinscheid D. J.,
Frey E. A.,
(1997) Enteropathogenic E. coli (EPEC) transfers its receptor for intimate adherence into mammalian cells. Cell 91:511–520.
(1995) Protein secretion by enteropathogenic Escherichia coli is essential for transducing signals to epithelial cells. Proc. Natl. Acad. Sci. USA 92:7991–7995.
Lai L. C.,
Finlay B. B.,
(1996) EspA, a protein secreted by enteropathogenic Escherichia coli, is required to induce signals in epthelial cells. Mol. Microbiol. 20:313–323.
(1997) Isolation of intact, high molecular weight DNA fragments for the E. coli genome project. Epicentre Forum 4:11–13.
Knutton S.,
Baldwin T.,
Williams P. H.,
McNeish A. S.
(1989) Actin accumulation at sites of bacterial adhesion to tissue culture cells: basis of a new diagnostic test for enteropathogenic and enterohemorrhagic Escherichia coli. Infect. Immun. 57:1290–1298.
Lloyd D. R.,
(1987) Adhesion of enteropathogenic Escherichia coli to human intestinal enterocytes and cultured human intestinal mucosa. Infect. Immun. 55:69–77.
Kumar S.,
Tamura K.,
Nei M.
(1993) MEGA: molecular evolutionary genetics analysis, 1.01 ed. (The Pennsylvania State University, University Park, Pa).
Stone K. D.,
(1997) A third secreted protein that is encoded by the enteropathogenic Escherichia coli pathogenicity island is required for transduction of signals and for attaching and effacing activities in host cells. Infect. Immun. 65:2211–2217.
Groisman E. A.,
Solomon F.,
(1995) Relationship between evolutionary rate and cellular location among the Inv/Spa invasion proteins of Salmonella enterica. Proc. Natl. Acad. Sci. USA 92:7252–7256.
Marciel A. M.,
Kapur V.,
Musser J. M.
(1997) Molecular population genetic analysis of a Streptococcus pyogenes bacteriophage-encoded hyaluronidase gene: recombination contributes to allelic variation. Microb. Pathog. 22:209–217.
(1995) A genetic locus of enterocyte effacement conserved among diverse enterobacterial pathogens. Proc. Natl. Acad. Sci. USA 92:1664–1668.
(1997) A cloned pathogenicity island from enteropathogenic Escherichia coli confers the attaching and effacing phenotype on E. coli K-12. Mol. Microbiol. 23:399–407.
Moon H. W.,
Whipp S. C.,
Argenzio M.,
Levine M.,
Gianella R. A.
(1983) Attaching and effacing activities of rabbit and human enteropathogenic Escherichia coli in pig and rabbit intestines. Infect. Immun. 41:1340–1351.
Nassif X.,
Lowy J.,
Stenberg P.,
O'Gaora P.,
Ganji A.,
So M.
(1993) Antigenic variation of pilin regulates adhesion of Neisseria meningitidis to human epilethial cells. Mol. Microbiol. 8:719–725.
(1994) Intergeneric transfer and recombination of the 6-phosphogluconate dehydrogenase gene (gnd) in enteric bacteria. Proc. Natl. Acad. Sci. USA 91:10227–10231.
Pósfai G.,
Koob M. D.,
(1997) Versatile insertion plasmids for targeted genome manipulations in bacteria: isolation, deletion, and rescue of the pathogenicity island LEE of the Escherichia coli O157:H7 genome. J. Bacteriol. 179:4426–4428.
Pupo G. M.,
Karaolis D. K. R.,
Lan R.,
Reeves P. R.
(1997) Evolutionary relationships among pathogenic and nonpathogenic Escherichia coli strains inferred from multilocus enzyme electrophoresis and mdh sequence studies. Infect. Immun. 65:2685–2692.
Rabinowitz R. P.,
Jarvis K.,
(1996) Attaching and effacing of host cells by enteropathogenic Escherichia coli: the absence of detectable tyrosine kinase mediated signal transduction. Microb. Pathog. 21:157–171.
Rosenshine I.,
(1992) Signal transduction between enteropathogenic Escherichia coli (EPEC) and epilethial cells: EPEC induces tyrosine phosphorylation of host cell proteins to initiate cytoskeletal rearrangement and bacterial uptake. EMBO J. 11:3551–3560.
Sharp P. M.,
Li W. H.
(1987) The codon adaptation index—a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 15:1281–1295.
Sun J.,
Inouye M.,
Inouye S.
(1991) Association of a retroelement with a P4-like cryptic prophage (retronphage φR73) integrated into the selenocystyl tRNA gene of Escherichia coli. J. Bacteriol. 173:4171–4181.
Theisen M.,
Borre M.,
Mathiesen M. J.,
Mikkelsen B.,
Lebech A.-M.,
Hansen K.
(1995) Evolution of the Borrelia burgdorferi outer surface protein OspC. J. Bacteriol. 177:3036–3044.
Gunzer F.,
DeMontigny L.,
Donohue-Rolfe A.
(1995) The role of the eaeA gene in diarrhea and neurological complications in a gnotobiotic piglet model of enterohemorrhagic Escherichia coli infection. Infect. Immun. 63:3621–3627.
Robins-Browne R. M.,
Gonis G.,
Hayes J.,
Whithers M.,
McCartney E.
(1985) Enteropathogenic Escherichia coli enteritis: evaluation of the gnotobiotic piglet as a model of human infection. Gut 26:570–578.
Wachsmuth I. K.,
Chapman C.,
Birden R.,
Brittingham J.,
Jackson C.,
Hogg J.
(1986) The pathogenesis of hemorrhagic colitis caused by Escherichia coli O175:H7 in gnotobiotic piglets. J. Infect. Dis. 154:712–716.
Wells J. G.,
Davis B. R.,
Riley L. W.,
Remis R. S.,
Sokolow R.,
Morris G. K.
(1983) Laboratory investigation of hemorrhagic colitis outbreaks associated with a rare Escherichia coli serotype. J. Clin. Microbiol. 18:512–520.
Wolfe M. L.,
Ørskov F.,
Ørskov I.,
Wilson R. A.
(1993) Clonal relationships among Escherichia coli strains that cause hemorrhagic colitis and infantile diarrhea. Infect. Immun. 61:1619–1629.
Wieler L. H.,
(1997) Insertion site of the locus of enterocyte effacement in enteropathogenic and enterohemorrhagic Escherichia coli differs in relation to the clonal phylogeny of the strains. FEMS Microbiol. Lett. 156:49–53.
(1992) Cloning and characterization of the eae gene of enterohaemorrhagic Escherichia coli O157:H7. Mol. Microbiol. 6:411–417.
Zhang Q.,
Wise K. S.
(1996) Molecular basis of size and antigenic variation of a Mycoplasma hominis adhesin encoded by divergent vaa genes. Infect. Immun. 64:2737–2744.
Zhao S.,
Mitchell S. E.,
Meng J.,
Doyle M. P.,
Kresovich S.
(1995) Cloning and nucleotide sequence of a gene upstream of the eaeA gene of enterohemorrhagic Escherichia coli O157:H7. FEMS Microbiol. Lett. 133:35–39.
Infection and Immunity Aug 1998, 66 (8) 3810-3817; DOI:
You are going to email the following Molecular Evolution of a Pathogenicity Island from Enterohemorrhagic Escherichia coli O157:H7