Abstract

Intimin is a highly polymorphic protein encoded by the eae gene and plays a crucial role in the attaching-effacing phenotype of diarrheagenic Escherichia coli and related pathogens. We have developed a method to quickly and accurately uncover allelic variation at the eae locus through the use of fluorescent RFLP (fRFLP). Application of fRFLP to 151 eae-positive strains (including the newly described Escherichia albertii) revealed 26 different fRFLP types that correspond to 20 of the 28 previously described eae alleles. Two sequence variants of the γ, ι, κ, and ζ alleles and three variants of ɛ were also observed. In addition to being reliable and accurate, the method can be easily adapted to accommodate new eae allelic sequences, as they become known.

Introduction

Strains of attaching and effacing Escherichia coli (AEEC) are capable of intimately attaching to intestinal epithelial cells (Finlayet al,1992; Kaper, 1998). Once the bacteria attach to the host cells, they manipulate the host cytoskeleton to form pedestal structures beneath the bacterial cells and, in the process, efface the microvilli of the intestinal mucosa (Hartlandet al,2000; Clearyet al,2004). This process creates a characteristic intestinal histopathology, termed attaching effacing (AE) lesions. The ability of pathogenic E. coli to form AE lesions is encoded in a pathogenicity island referred to as the locus of enterocyte effacement or LEE (McDanielet al,1995; McDaniel & Kaper, 1997). The LEE island is ∼35kb in length and can be inserted into one of several chromosomal locations in AEEC clonal groups (Wieleret al,1997; Elliottet al,1998; Pernaet al,1998; Sperandioet al,1998). The LEE island comprises ∼40 genes including those that encode the structural components of a type III secretion system (TTSS) (Elliottet al,1998; Pernaet al,1998). The TTSS acts as a molecular delivery system to translocate bacterial effector proteins, many of which are also LEE-encoded, into the host cell (Denget al,2005).

Intimin plays a crucial role in the AE phenotype and is encoded by eae, one of the most highly polymorphic genes of the LEE island (Castilloet al,2005). Intimin is composed of six domains: a periplasmic domain, a transmembrane domain, three extracellular immunoglobulin-like domains, and an extracellular lectin-like domain (Luoet al,2000). Much of the transmembrane domain is homologous to the invasins of pathogenic Yersinia and this part of the molecule has been termed the central conserved domain (McGrawet al,1999). The extracellular domains of intimin interact with cellular receptors (Sinclairet al,2006) including Tir, the translocated intimin receptor that is encoded in LEE and moves to the eukaryotic cell via the LEE-encoded TTSS (Kennyet al,1997). An analysis of 27 intimin alleles from GenBank indicates that the four C-terminal extracellular domains contain more than 75% of the nucleotide variation present within those sequences (unpublished data). The AE strains represent a variety of enteric pathotypes of both humans and animals, including both typical and atypical enteropathogenic E. coli (EPEC), enterhemorrhagic E. coli (EHEC), Escherichia albertii, and Citrobacter rodentium (Donnenberg, 2002; Waleset al,2005). Among the AEEC, there is a striking association of different eae alleles with specific clones or clonal groups of pathogenic E. coli (Wieleret al,1997; Adu-Bobieet al,1998; Tarr & Whittam, 2002). This association of eae alleles with pathogenic lineages has been used to help classify strains into pathotypes, so a rapid and reliable intimin typing scheme is a valuable tool.

Here we describe a method to identify allelic variants of the eae locus through the use of fluorescent RFLP (fRFLP; Lazzaroet al,2002). In this method, the entire highly variable 3′ half of eae is amplified in a standard PCR reaction using primers that are located in the conserved central domain of eae and in the conserved downstream gene escD. The exact size of the amplicon depends on the specific eae allele, but is typically about 2kb in length. The PCR amplicon is then digested with one or more restriction enzymes that leave a 5′ overhang that acts as a template for the incorporation of a fluorescent dye-terminator nucleotide from a standard cycle sequencing kit. The multiple labeled restriction fragments are then separated on a capillary-based sequencer and their sizes estimated to within a few base pairs in length. This fRFLP system provides a rapid method for uncovering allelic variation in eae and for classifying eae subtypes based on complex restriction digests.

Materials and methods

Bacterial strains and DNA isolation

The strains in this study included 144 AEEC representing 38 O-serogroups, many of which comprise EPEC serotypes, and seven strains from the Escherichia albertii and Shigella boydii 13 clonal lineage (Huyset al,2003; Hymaet al,2005). All 151 AE strains were grown overnight at 37°C in 10mL of LB broth with moderate shaking. Genomic DNA was isolated using the Puregene DNA isolation kit (Gentra Systems Incet al, Minneapolis, MN). DNA concentrations were determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies Incet al, Rockland, DE). Working template concentrations of 50ngμL−1 were used for PCR.

PCR and eae-escD primer design

PCR primers were designed in the central conserved domain of eae (eae-F1: 5′-ACT CCG ATT CCT CTG GTG AC-3′) and the conserved downstream gene escD (escD-R1: 5′-GTA TCA ACA TCT CCC GCC CA-3′) based on available sequences for the α, β, γ, ɛ, ζ, and θ intimin alleles. The eae-F1 and escD-R1 primers are located at positions 25986–26005 and 27918–27937, respectively, of the complete LEE sequence from strain E2348/69 (GenBank accession number AF022236). Each 25μL reaction contained 2.5μL 10 × buffer II (Applied Biosystems, Foster City, CA), 2.5μL 2mM dNTP, 2.0μL 25mM MgCl2, 0.5μL 10μM eae-F1 primer, 0.5μL 10μM escD-R1 primer, 1.5U AmpliTaq Gold (Applied Biosystems), 1μL 50ngμL−1 genomic DNA template, and 15.7μL ddH2O. Amplification utilized an initial denaturing step at 94°C for 10min, followed by 35 cycles of 92°C for 1min, 55°C for 1min, and 72°C for 2min. A final step of 72°C for 5min was used for final completion of any partially extended product. PCR products (5μL) were visualized on ethidium bromide-stained 0.8% agarose gels by illumination with UV light.

fRFLP development

In silico digestions were performed to find suitable restriction endonucleases for fRFLP. To be considered for fRFLP, the restriction enzymes had to (1) leave a 5′ overhang, (2) mainly produce fragments in the range of the labeled size standard (60–640bp), and (3) produce fragments that could be labeled with fluorescent tags that were different than that of the size standard. Out of 532 restriction enzymes tested in silico, MseI was chosen for fRFLP because it meets all of the above criteria and was found to have the greatest power to discriminate known eae alleles (Figs 1 and 2).

1

Location of PCR primers and predicted MseI restriction sites for four major intimin alleles associated with human disease. The forward primer (eae-F1) is located in the central conserved domain of eae and the reverse primer (escD-R1) is located in the conserved downstream gene escD. The resulting amplicon is ∼1800–2100bp depending on the allele.

1

Location of PCR primers and predicted MseI restriction sites for four major intimin alleles associated with human disease. The forward primer (eae-F1) is located in the central conserved domain of eae and the reverse primer (escD-R1) is located in the conserved downstream gene escD. The resulting amplicon is ∼1800–2100bp depending on the allele.

2

Fragment sizes of DNA resulting from MseI digestion of the eae-escD amplicon based on in silico analysis of known sequence of the eae locus. The number of restriction fragments ranges from 16 to 22 with an average of 19.1 fragments per amplicon. Two additional patterns were resolved for lane 5 (γ1.1 and γ1.2) and lane 7 (ɛ1.2 and ɛ1.3) by AseI/DdeI/SalI restriction digestions. The gray line denotes the 60nt length that is the lower limit of fragment size determined by fRFLP.

2

Fragment sizes of DNA resulting from MseI digestion of the eae-escD amplicon based on in silico analysis of known sequence of the eae locus. The number of restriction fragments ranges from 16 to 22 with an average of 19.1 fragments per amplicon. Two additional patterns were resolved for lane 5 (γ1.1 and γ1.2) and lane 7 (ɛ1.2 and ɛ1.3) by AseI/DdeI/SalI restriction digestions. The gray line denotes the 60nt length that is the lower limit of fragment size determined by fRFLP.

To confirm the eae allele assignments based on the MseI findings, a second restriction digestion was designed. Multiple enzymes were needed to give a similar amount of pattern variation as seen with MseI. For simplicity, restriction enzymes were selected that use the same reaction conditions (reaction buffer, incubation temperature, and labeled ddNTP). Three sets of triple-enzyme digests were found. A combination of AseI, DdeI, and SalI was chosen over the others based on lower enzyme cost and the ddNTP used is the same as that for the MseI digest. The two digests (MseI and AseI/DdeI/SalI) were used to subtype the intimin alleles from a diverse set of 151 AE strains.

PCR clean-up and restriction enzyme digestion

Samples that were positive for the eae-escD PCR were treated with ExoSAP-IT (USB Corporation, Cleveland, OH) to remove unincorporated dNTPs and PCR primers (5μL of PCR product and 2μL of ExoSAP-IT). All restriction endonucleases, reaction buffers, and bovine serum albumin (BSA) solutions were obtained from New England Biolabs Inc. (Ipswich, MA). MseI digests were set up so that each 30μL reaction contained 7.0μL ExoSAP-treated PCR product, 0.3μL 100 × BSA, 0.5μL 10UμL−1MseI, 3.0μL 10 × NEB2 buffer, and 19.2μL ddH2O. AseI/DdeI/SalI digests were set up so that each 30μL reaction contained 7.0μL ExoSAP-treated PCR product, 0.3μL 100 × BSA, 0.5μL 10UμL−1AseI, 0.5L 10UμL−1DdeI, 0.25μL 20UμL−1SalI, 3.0μL 10 × NEB3 buffer, and 18.45μL ddH2O. Reactions were incubated overnight at 37°C.

fRFLP

fRFLP was performed using the CEQ DTCS standard kit (Beckman Coulter Inc, Fullerton, CA). Each reaction contained 2.0μL of unpurified restriction enzyme digest, 1.5μL 10 × reaction buffer, 0.1μL ddUTP, 0.1μL Taq DNA polymerase, and 11.3μL ddH2O. Samples were incubated at 60°C for 1h, purified with Sephadex G-50 Fine columns (Amersham Pharmacia Biotech Incet al, Piscataway, NJ), dried under vacuum centrifugation (Savant Instruments Incet al, Holbrook, NY), and suspended in 10μL of deionized formamide. Of this, 2μL were mixed with 0.6μL CEQ DNA Size Standard 600 (Beckman Coulter Incet al, Fullerton, CA), and 39.4μL deionized formamide, and run on a CEQ2000XL (Beckman Coulter Inc.) using a capillary temperature of 50°C, a denature step at 90°C for 2min, injection at 2.0kV for 30s, and separation at 4.8kV for 65min (MseI digest) or 90min (AseI/DdeI/SalI digest). Fragment sizes were determined with the CEQ2000XL software, version 4.3.9.

DNA sequencing

The complete eae gene was sequenced in at least one representative strain for each fRFLP pattern identified. The 5′ half of eae was amplified using primers cesT-F9 (5′-TCA GGG AAT AAC ATT AGA AA-3′) and eae-R3 (5′-TCT TGT GCG CTT TGG CTT-3′) using the same PCR conditions described above. PCR products were purified using the QIAquick PCR purification kit (QIAGEN Incet al, Valencia, CA) and quantified with a NanoDrop ND-1000 spectrophotometer. Cycle sequencing reactions contained 6.0μL CEQ DTCS Quick Start premix (Beckman Coulter Inc.), 1.5μL 20μM primer, c. 180ng of cesT/eae product or 250ng of eae/escD product, and ddH2O to 15μL. Amplification utilized an initial denaturing step at 94°C for 1min, followed by 35 cycles of 96°C for 30s, 52°C for 30s, and 60°C for 2min. Upon completion of cycle sequencing, samples were purified with Sephadex G-50 Fine columns, dried under vacuum centrifugation, suspended in 40μL of deionized formamide, and sequenced using a Beckman CEQ2000XL DNA sequencer. Samples were analyzed using the CEQ2000XL software and then exported for further analysis with the SeqMan and MegAlign modules of the Lasergene software (DNASTAR Incet al, Madison, WI). Internal sequencing primers were designed as new sequence data were generated.

Results

The ability of the eae fRFLP method to resolve intimin alleles was first tested on 53 strains based on the expected patterns from the in silico digestion of available eae sequences from GenBank. This initial study included 34 AE strains, most of which had previously been known to have the α, β, or γ intimin alleles, most often found in strains associated with human infection (Table 1). Nineteen additional AE strains were selected and tested by the fRFLP method because their intimin allele had been determined by other researchers or by previous work in our laboratory based on DNA sequencing or conventional RFLP (Table 2). The intimin alleles of all 53 strains were confirmed by both the MseI and AseI/DdeI/SalI digests.

1

Allelic variation in eae based on fRFLP

Accession number Strain name Serotype eae fRFLP pattern Source and pathotype (reference) 
TW00588 DEC 1a (572-56) O55:H6 α 1.1 Human EPEC (Rodrigueset al,1996
TW07884 E851/71 O142:H6 α 1.1 Human EPEC (Craviotoet al,1979
TW07923 RN587/1 O157:[h45] α 1.1 Human EPEC (Blanket al,2003
TW04262 TB269C O145:[h34] α 2.1 Human atypical EPEC (Boketeet al,1997
TW01120 B170 O111:[h2] β 1.1 Human EPEC (Moyenuddinet al,1989
TW05355 13180-25 O111:H11 β 1.1 EHEC from food (Fenget al,1998
TW00148 3448-87 O114:H2 β 1.1 Human EPEC 
TW00389 29315 (AC-C12) O119:H2 β 1.1 Human EPEC (Goncalveset al,1997
TW07099 LT119-80 O119:H2 β 1.1 Human EPEC (Goncalveset al,1997
TW01266 C342-62 O126:H2 β 1.1 Human EPEC (Orskovet al,1990
TW07896 E56/54 O128:H2 β 1.1 Human EPEC (Robins-Browneet al,1993
TW01664 DEC 10i (87-1713) O145:H16 β 1.1 Human EHEC 
TW05149 BCL73 O145:[h-] β 1.1 Bovine STEC 
TW07860 314-S O145:[h16] β 1.1 Bovine STEC 
TW08894 02-3422 O145:[h2] β 1.1 Rabbit EPEC (Garciaet al,2002
TW09153 IH 16 O145:[h-] β 1.1 Human STEC 
TW07924 Z188-93 O110:H6 β 2.1 Avian EPEC (Blanket al,2000
TW01225 1396/69 O119:H6 β 2.1 Human EPEC (Craviotoet al,1979
TW03293 ECOR-37 O-:[h7] γ 1.1 Marmoset atypical EPEC (McGrawet al,1999
TW00947 DEC 5d (C586-65) O55:H7 γ 1.1 Human atypical EPEC (McGrawet al,1999
TW03064 B6820-C1 O145:[h28] γ 1.1 Bovine STEC 
TW07596 GS G5578620 O145:[h28] γ 1.1 Human STEC (Feyet al,2000
TW07865 IHIT0304 O145:H28 γ 1.1 Bovine STEC 
TW08087 MT#66 O145:[h28] γ 1.1 Human STEC (Jelacicet al,2003
TW09356 4865/96 O145:[h28] γ 1.1 Human STEC (Unkmeir & Schmidt, 2000
TW08881 3556-77 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW08887 3557-77 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW08889 3053-94 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW07618 98ST607 O110:H28 ζ 1.1 Human STEC 
TW00964 75-83 O145:[h25] ζ 1.1 Human STEC (Boppet al,1987
TW05307 LTO55-43 O55:H7 θ 1.1 Human atypical EPEC (Rodrigueset al,1996
TW07960 DA-34 O103:[h25] θ 1.1 Human STEC 
TW00970 DEC 8b (3030A-86) O111:H8 θ 1.1 Human EHEC (Tarr & Whittam, 2002
TW07888 010-311082 O76:H51 μ 1.1 Human EPEC (Blanket al,2000
Accession number Strain name Serotype eae fRFLP pattern Source and pathotype (reference) 
TW00588 DEC 1a (572-56) O55:H6 α 1.1 Human EPEC (Rodrigueset al,1996
TW07884 E851/71 O142:H6 α 1.1 Human EPEC (Craviotoet al,1979
TW07923 RN587/1 O157:[h45] α 1.1 Human EPEC (Blanket al,2003
TW04262 TB269C O145:[h34] α 2.1 Human atypical EPEC (Boketeet al,1997
TW01120 B170 O111:[h2] β 1.1 Human EPEC (Moyenuddinet al,1989
TW05355 13180-25 O111:H11 β 1.1 EHEC from food (Fenget al,1998
TW00148 3448-87 O114:H2 β 1.1 Human EPEC 
TW00389 29315 (AC-C12) O119:H2 β 1.1 Human EPEC (Goncalveset al,1997
TW07099 LT119-80 O119:H2 β 1.1 Human EPEC (Goncalveset al,1997
TW01266 C342-62 O126:H2 β 1.1 Human EPEC (Orskovet al,1990
TW07896 E56/54 O128:H2 β 1.1 Human EPEC (Robins-Browneet al,1993
TW01664 DEC 10i (87-1713) O145:H16 β 1.1 Human EHEC 
TW05149 BCL73 O145:[h-] β 1.1 Bovine STEC 
TW07860 314-S O145:[h16] β 1.1 Bovine STEC 
TW08894 02-3422 O145:[h2] β 1.1 Rabbit EPEC (Garciaet al,2002
TW09153 IH 16 O145:[h-] β 1.1 Human STEC 
TW07924 Z188-93 O110:H6 β 2.1 Avian EPEC (Blanket al,2000
TW01225 1396/69 O119:H6 β 2.1 Human EPEC (Craviotoet al,1979
TW03293 ECOR-37 O-:[h7] γ 1.1 Marmoset atypical EPEC (McGrawet al,1999
TW00947 DEC 5d (C586-65) O55:H7 γ 1.1 Human atypical EPEC (McGrawet al,1999
TW03064 B6820-C1 O145:[h28] γ 1.1 Bovine STEC 
TW07596 GS G5578620 O145:[h28] γ 1.1 Human STEC (Feyet al,2000
TW07865 IHIT0304 O145:H28 γ 1.1 Bovine STEC 
TW08087 MT#66 O145:[h28] γ 1.1 Human STEC (Jelacicet al,2003
TW09356 4865/96 O145:[h28] γ 1.1 Human STEC (Unkmeir & Schmidt, 2000
TW08881 3556-77 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW08887 3557-77 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW08889 3053-94 Boydii 13 γ 1.2 Human atypical B13 (Hymaet al,2005
TW07618 98ST607 O110:H28 ζ 1.1 Human STEC 
TW00964 75-83 O145:[h25] ζ 1.1 Human STEC (Boppet al,1987
TW05307 LTO55-43 O55:H7 θ 1.1 Human atypical EPEC (Rodrigueset al,1996
TW07960 DA-34 O103:[h25] θ 1.1 Human STEC 
TW00970 DEC 8b (3030A-86) O111:H8 θ 1.1 Human EHEC (Tarr & Whittam, 2002
TW07888 010-311082 O76:H51 μ 1.1 Human EPEC (Blanket al,2000
*

Lower case H-types in brackets were inferred from fliC allele determination.

fRFLP, fluorescent restriction fragment length polymorphism.

2

 Reference strains for defined intimin alleles and patterns determined by fRFLP and confirmed by DNA sequencing

Pattern no. Intimin allele eae fRFLP pattern Reference strain Species and serotype GenBank number (reference) 
α α 1.1 TW06375 (E2348/69) E. coli O127:H6 AF022236 (Elliottet al,1998
α2 α 2.1 TW01270 (C712-65) E. coli O125:H6 DQ523600 (this study) 
β β 1.1 TW07862 (413/89-1) E. coli O26:[h11] AJ277443 
β2 β 2.1 TW07894 (0659-79) E. coli O119:H6 DQ523605 (this study) 
γ γ 1.1 TW08264 (Sakai) E. coli O157:H7 NC002695 (Hayashiet al,2001
γ γ 1.2 TW08888 (3092-94) S. boydii type 13 DQ523608 (this study) 
ɛ ɛ 1.1 TW08101 (MT#80) E. coli O103:H2 DQ523606 (this study) 
ɛ ɛ 1.2 TW08023 (MT#2) E. coli O121:H19 AY186750 (Tarret al,2002
ɛ ɛ 1.3 TW10363 (83F4) E. coli O-:[h8] DQ523612 (this study) 
10 ɛ2 ɛ 2.1 TW10371 (98B3) E. coli O116:[h9] DQ523614 (this study) 
11 ζ ζ 1.1 TW07863 (537/89) E. coli O84:[h2] AJ298279 (Joreset al,2003
12 ζ ζ 1.2 TW04892 (921) E. coli O111:H9 AF449417 (Tarr & Whittam, 2002
13 η η 1.1 TW07892 (012-050982) E. coli O142:[h21] DQ523604 (this study) 
14 θ θ 1.1 TW01387 (CL-37) E. coli O111:H8 AF449418 (Tarr & Whittam, 2002
15 ι ι 1.1 TW01933 (1252-59) E. coli O55:[h34] DQ523601 (this study) 
16 ι ι 1.2 TW04174 (TB227C) E. coli O86:[h8] DQ523602 (this study) 
17 ι2 ι 2.1 TW08839 (C-425) S. boydii type 13 AY696842 (Hymaet al,2005
18 κ κ 1.1 TW06584 (C295-53) E. coli O86:H34 DQ523603 (this study) 
19 κ κ 1.2 TW10337 (64B4) E. coli O49:[h10] DQ523611 (this study) 
20 λ λ 1.1 TW10327 (57A1) E. coli O33:[h34] DQ523609 (this study) 
21 μ μ 1.1 TW08260 (MA551/1) E. coli O55:[h51] DQ523607 (this study) 
22 ν ν 1.1 TW10376 (106A5) E. albertii DQ523615 (this study) 
23 ξ ξ 1.1 TW10334 (60A3) E. coli O5:[h2] DQ523610 (this study) 
24 o 1.1 TW07627 (Albert 19982) E. albertii AY696838 (Hymaet al,2005
25 ρ ρ 1.1 TW10366 (93I4) E. coli O21:[h5] DQ523613 (this study) 
26 τ τ 1.1 TW08933 (K-1) S. boydii type 7 AY696839 (Hymaet al,2005
Pattern no. Intimin allele eae fRFLP pattern Reference strain Species and serotype GenBank number (reference) 
α α 1.1 TW06375 (E2348/69) E. coli O127:H6 AF022236 (Elliottet al,1998
α2 α 2.1 TW01270 (C712-65) E. coli O125:H6 DQ523600 (this study) 
β β 1.1 TW07862 (413/89-1) E. coli O26:[h11] AJ277443 
β2 β 2.1 TW07894 (0659-79) E. coli O119:H6 DQ523605 (this study) 
γ γ 1.1 TW08264 (Sakai) E. coli O157:H7 NC002695 (Hayashiet al,2001
γ γ 1.2 TW08888 (3092-94) S. boydii type 13 DQ523608 (this study) 
ɛ ɛ 1.1 TW08101 (MT#80) E. coli O103:H2 DQ523606 (this study) 
ɛ ɛ 1.2 TW08023 (MT#2) E. coli O121:H19 AY186750 (Tarret al,2002
ɛ ɛ 1.3 TW10363 (83F4) E. coli O-:[h8] DQ523612 (this study) 
10 ɛ2 ɛ 2.1 TW10371 (98B3) E. coli O116:[h9] DQ523614 (this study) 
11 ζ ζ 1.1 TW07863 (537/89) E. coli O84:[h2] AJ298279 (Joreset al,2003
12 ζ ζ 1.2 TW04892 (921) E. coli O111:H9 AF449417 (Tarr & Whittam, 2002
13 η η 1.1 TW07892 (012-050982) E. coli O142:[h21] DQ523604 (this study) 
14 θ θ 1.1 TW01387 (CL-37) E. coli O111:H8 AF449418 (Tarr & Whittam, 2002
15 ι ι 1.1 TW01933 (1252-59) E. coli O55:[h34] DQ523601 (this study) 
16 ι ι 1.2 TW04174 (TB227C) E. coli O86:[h8] DQ523602 (this study) 
17 ι2 ι 2.1 TW08839 (C-425) S. boydii type 13 AY696842 (Hymaet al,2005
18 κ κ 1.1 TW06584 (C295-53) E. coli O86:H34 DQ523603 (this study) 
19 κ κ 1.2 TW10337 (64B4) E. coli O49:[h10] DQ523611 (this study) 
20 λ λ 1.1 TW10327 (57A1) E. coli O33:[h34] DQ523609 (this study) 
21 μ μ 1.1 TW08260 (MA551/1) E. coli O55:[h51] DQ523607 (this study) 
22 ν ν 1.1 TW10376 (106A5) E. albertii DQ523615 (this study) 
23 ξ ξ 1.1 TW10334 (60A3) E. coli O5:[h2] DQ523610 (this study) 
24 o 1.1 TW07627 (Albert 19982) E. albertii AY696838 (Hymaet al,2005
25 ρ ρ 1.1 TW10366 (93I4) E. coli O21:[h5] DQ523613 (this study) 
26 τ τ 1.1 TW08933 (K-1) S. boydii type 7 AY696839 (Hymaet al,2005
*

Lower case H-types in brackets were inferred from fliC allele determination.

fRFLP pattern discovered in AEEC strains from Guinea-Bissau study (Valentiner-Branthet al,2003).

Variants of an intimin allele that are indistinguishable by fRFLP with the MseI digest.

§

Variants of an intimin allele that are indistinguishable by fRFLP with the AseI/DdeI/SalI triple digest.

Benkel P. and Chakraborty T. (unpublished).

fRFLP, fluorescent restriction fragment length polymorphism; E. coli, Escherichia coli; S. boydii, Shigella boydii; E. albertii, Escherichia albertii; AEEC, attaching and effacing E. coli.

To further evaluate the new method, we then tested 98 eae-positive strains for which no allelic subtyping data existed. These strains were originally recovered from two separate populations: a pediatric population in Seattle, Washington (Dennoet al,2005) and a cohort study of childhood diarrheal disease in Guinea-Bissau, West Africa (Valentiner-Branthet al,2003). Among these 98 eae-positive strains, most strains (85%) exhibited known digestion patterns with both MseI and the triple digest, and therefore could be easily subtyped to an allele. The 15 strains for which the eae allele could not be determined had one of seven previously unobserved digestion patterns. A representative of each pattern was sequenced and the alleles were either previously described or variants of previously described alleles (ɛ1.3, ɛ2.1, κ1.2, λ1.1, ν1.1, ρ1.1, and ξ1.1).

Among the 151 strains examined, there were a total of 24 different fRFLP patterns observed for the MseI digests (Fig. 2). The AseI/DdeI/SalI digest also resolved 24 distinct fRFLP patterns (data not shown). For two alleles, ɛ1.2/ɛ1.3 and γ1.2/γ1.3 variants are resolved by the AseI/DdeI/SalI digests and are indistinguishable with MseI digestion. In contrast, two variants that are resolved by MseI (ι1.1/ι1.2 and ζ1.1/ζ1.2) are indistinguishable with the triple digest. Combined there are 26 distinct fRFLP patterns (see Table S1). In all, the new fRFLP method is able to identify 20 of the 28 previously described alleles, and differentiated two variants of the γ, ι, κ and ζ alleles, and three variants of ɛ.

To gauge the accuracy of the fragment size estimation, we compared the expected fragment sizes based on in silico digestions to those observed from the capillary sequencer. Overall the fragment scoring was accurate with most fragments over 100nt in length less than 2% different in size from their expected values (Fig. 3). Examination of the plot reveals that for fragments less than ∼100nt and greater than ∼550nt, the estimated fragments size deviates from the expected (Fig. 3). For the 83 fragments observed in the MseI digest, the average deviation from expected size is 1.76%. The triple digest was more accurate with an average deviation of 1.01% across the 70 distinct fragments.

3

Accuracy of fragments sizes estimated by fRFLP. The percentage observed difference from the expected fragment size is plotted against the expected size. The lines mark the percentage deviation for 1, 2, and 5nt, respectively.

3

Accuracy of fragments sizes estimated by fRFLP. The percentage observed difference from the expected fragment size is plotted against the expected size. The lines mark the percentage deviation for 1, 2, and 5nt, respectively.

Discussion

Current eae typing schemes either focus on allele-specific PCR amplification (Adu-Bobieet al,1998; Reidet al,1999; Oswaldet al,2000; Zhanget al,2002; Blancoet al,2003; Joreset al,2003; Ordenet al,2003) or conventional RFLP analysis (Schmidtet al,1993; Oswaldet al,2000; Jenkinset al,2003; Ramachandranet al,2003). Sequence analysis of intimin alleles has revealed that many eae alleles are mosaics with segments having different evolutionary histories (McGrawet al,1999; Tarr & Whittam, 2002). Therefore, allele-specific PCR amplification can lead to erroneous typing results. For example, the μ allele would be erroneously typed as γ by the Reid method (Reidet al,1999) because the γ allele primer is located in a region that is shared between these two alleles (data not shown). In addition, as new alleles are discovered, new primers are necessary to amplify these alleles and the subtyping scheme becomes more complex.

Conventional RFLP analysis also has its limitations. The 5′ half of eae is relatively conserved among the different alleles, so there may not be sufficient variation in the amplicons to accurately and reliably differentiate the alleles. For example, EPEC strain 1396/69 (O119:H6) was typed by Jenkins (2003) as possessing the γ allele of eae instead of the β2 allele. In silico analysis revealed that the β2, γ, and λ alleles of eae are virtually indistinguishable under the Jenkins system. Another limitation of conventional RFLP is that a system based on the highly variable 3′ half of the gene may be difficult to score since small differences in the banding patterns of different alleles may not be easily discernible under standard electrophoretic conditions.

The fRFLP method described here addresses many of the limitations of the existing eae typing methods. A drawback of the devised method, however, is that it requires the location of escD relative to eae to be conserved. If this location changes such that escD is either upstream or far downstream of eae, PCR amplification will not occur and the strain will be nontypeable by this method. However, most of the known alleles for eae have been amplified and observed, so this situation does not appear to be common. We have tested the new subtyping method against a diverse panel of 87 eae-positive isolates originally recovered from children in West Africa and also uncovered seven additional variants. The remaining alleles (β3, ɛ3, ɛ4, η2, π, σ) still need to be tested. Another limitation is that some of the eae alleles in GenBank do not contain their associated downstream escD sequence, so the expected amplicon size and fRFLP profile cannot be fully determined in all cases. When new fRFLP patterns are observed, eae needs to be completely sequenced to verify the allele, but then the new pattern can be added to the fRFLP database for future reference.

The fRFLP method of eae allelic subtyping is easy to perform, cost-effective, and sensitive. The principal advantages are that it makes use of existing technology that is widely available, and can reveal previously unknown allelic variants. In the future, we expect more sensitive, high throughput methods for detecting nucleotide polymorphisms, such as those based on liquid microsphere suspension (Gilmouret al,2006), will provide quantitative technologies for rapid allelic identification of intimin and other crucial virulence factors.

Acknowledgements

This project has been funded in part with Federal funds from the NIAID, NIH, DHHS, under NIH Research Contract # N01-AI-30058.

References

Adu-Bobie
J
Frankel
G
Bain
C
Goncalves
AG
Trabulsi
LR
Douce
G
Knutton
S
Dougan
G
(
1998
)
Detection of intimins α, β, γ, and δ, four intimin derivatives expressed by attaching and effacing microbial pathogens
.
J Clin Microbiol
 
36
:
662
668
.
Blanco
M
Blanco
JE
Mora
A
, et al.  . (
2003
)
Serotypes, virulence genes, and intimin types of Shiga toxin (verotoxin)-producing Escherichia coli isolates from healthy sheep in Spain
.
J Clin Microbiol
 
41
:
1351
1356
.
Blank
TE
Zhong
H
Bell
AL
Whittam
TS
Donnenberg
MS
(
2000
)
Molecular variation among type IV pilin (bfpA) genes from diverse enteropathogenic Escherichia coli strains
.
Infect Immun
 
68
:
7028
7038
.
Blank
TE
Lacher
DW
Scaletsky
IC
Zhong
H
Whittam
TS
Donnenberg
MS
(
2003
)
Enteropathogenic Escherichia coli O157 strains from Brazil
.
Emerg Infect Dis
 
9
:
113
115
.
Bokete
TN
Whittam
TS
Wilson
RA
Clausen
CR
O'Callahan
CM
Moseley
SL
Fritsche
TR
Tarr
PI
(
1997
)
Genetic and phenotypic analysis of Escherichia coli with enteropathogenic characteristics isolated from Seattle children
.
J Infect Dis
 
175
:
1382
1389
.
Bopp
CA
Greene
KD
Downes
FP
Sowers
EG
Wells
JG
Wachsmuth
IK
(
1987
)
Unusual verotoxin-producing Escherichia coli associated with hemorrhagic colitis
.
J Clin Microbiol
 
25
:
1486
1489
.
Castillo
A
Eguiarte
LE
Souza
V
(
2005
)
A genomic population genetics analysis of the pathogenic enterocyte effacement island in Escherichia coli: the search for the unit of selection
.
Proc Natl Acad Sci USA
 
102
:
1542
1547
.
Cleary
J
Lai
LC
Shaw
RK
Straatman-Iwanowska
A
Donnenberg
MS
Frankel
G
Knutton
S
(
2004
)
Enteropathogenic Escherichia coli (EPEC) adhesion to intestinal epithelial cells: role of bundle-forming pili (BFP), EspA filaments and intimin
.
Microbiology
 
150
:
527
538
.
Cravioto
A
Gross
RJ
Scotland
SM
Rowe
B
(
1979
)
An adhesive factor found in strains of Escherichia coli belonging to the traditional infantile enteropathogenic serotypes
.
Curr Microbiol
 
3
:
95
99
.
Deng
W
Li
Y
Hardwidge
PR
Frey
EA
Pfuetzner
RA
Lee
S
Gruenheid
S
Strynakda
NC
Puente
JL
Finlay
BB
(
2005
)
Regulation of type III secretion hierarchy of translocators and effectors in attaching and effacing bacterial pathogens
.
Infect Immun
 
73
:
2135
2146
.
Donnenberg
MS
(
2002
)
Escherichia coli: Virulence Mechanism of a Versatile Pathogens
(
Donnenberg
MS
, ed), pp.
xxi
xxv
.
Academic Press
, New York.
Denno
DM
Stapp
JR
Boster
DR
Qin
X
Clausen
CR
Del Beccaro
KH
Swerdlow
DL
Braden
CR
Tarr
PI
(
2005
)
Etiology of diarrhea in pediatric outpatient settings
.
Pediatr Infect Dis J
 
24
:
142
148
.
Elliott
SJ
Wainwright
LA
McDaniel
TK
Jarvis
KG
Deng
YK
Lai
LC
McNamara
BP
Donnenberg
MS
Kaper
JB
(
1998
)
The complete sequence of the locus of enterocyte effacement (LEE) from enteropathogenic Escherichia coli E2348/6
.
Mol Microbiol
 
28
:
1
4
.
Feng
P
Lampel
KA
Karch
H
Whittam
TS
(
1998
)
Genotypic and phenotypic changes in the emergence of Escherichia coli O157:H7
.
J Infect Dis
 
177
:
1750
1753
.
Fey
PD
Wickert
RS
Rupp
ME
Safranek
TJ
Hinrichs
SH
(
2000
)
Prevalence of non-O157:H7 Shiga toxin-producing Escherichia coli in diarrheal stool samples from Nebraska
.
Emerg Infect Dis
 
6
:
530
533
.
Finlay
BB
Rosenshine
I
Donnenberg
MS
Kaper
JB
(
1992
)
Cytoskeletal composition of attaching and effacing lesions associated with enteropathogenic Escherichia coli adherence to HeLa cells
.
Infect Immun
 
60
:
2541
2543
.
Garcia
A
Marini
RP
Feng
Y
Vitsky
A
Knox
KA
Taylor
NS
Schauer
DB
Fox
JG
(
2002
)
A naturally occurring rabbit model of enterohemorrhagic Escherichia coli-induced disease
.
J Infect Dis
 
186
:
1682
1686
.
Gilmour
MW
Tracz
DM
Andrysiak
AK
Clark
CG
Tyson
S
Severini
A
Ng
LK
(
2006
)
Use of the espZ gene encoded in the locus of enterocyte effacement for molecular typing of Shiga toxin-producing Escherichia coli
.
J Clin Microbiol
 
44
:
449
458
.
Goncalves
AG
Campos
LC
Gomes
TA
Rodrigues
J
Sperandio
V
Whittam
TS
Trabulsi
LR
(
1997
)
Virulence properties and clonal structures of strains of Escherichia coli O119 serotypes
.
Infect Immun
 
65
:
2034
2040
.
Hartland
EL
Daniell
SJ
Delahay
RM
Neves
BC
Wallis
T
Shaw
RK
Hale
C
Knutton
S
Frankel
G
(
2000
)
The type III protein translocation system of enteropathogenic Escherichia coli involves EspA–EspB protein interactions
.
Mol Microbiol
 
35
:
1483
1492
.
Hayashi
T
Makino
K
Ohnishi
M
, et al.  . (
2001
)
Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12
.
DNA Res
 
8
:
11
22
.
Huys
G
Cnockaert
M
Janda
JM
Swings
J
(
2003
)
Escherichia albertii sp. nov., a diarrhoeagenic species isolated from stool specimens of Bangladeshi children
.
Int J Syst Evol Microbiol
 
53
:
807
810
.
Hyma
KE
Lacher
DW
Nelson
AM
Bumbaugh
AC
Janda
JM
Strockbine
NA
Young
VB
Whittam
TS
(
2005
)
Evolutionary genetics of a new pathogenic Escherichia species: Escherichia albertii and related Shigella boydii strains
.
J Bacteriol
 
187
:
619
628
.
Jelacic
JK
Damrow
T
Chen
GS
Jelacic
S
Bielaszewska
M
Ciol
M
Carvalho
HM
Melton-Celsa
AR
O'Brien
AD
Tarr
PI
(
2003
)
Shiga toxin-producing Escherichia coli in Montana: bacterial genotypes and clinical profiles
.
J Infect Dis
 
188
:
719
729
.
Jenkins
C
Lawson
AJ
Cheasty
T
Willshaw
GA
Wright
P
Dougan
G
Frankel
G
Smith
HR
(
2003
)
Subtyping intimin genes from enteropathogenic Escherichia coli associated with outbreaks and sporadic cases in the United Kingdom and Eire
.
Mol Cell Probes
 
17
:
149
156
.
Jores
J
Zehmke
K
Eichberg
J
Rumer
L
Wieler
LH
(
2003
)
Description of a novel intimin variant (type zeta) in the bovine O84:NM verotoxin-producing Escherichia coli strain 537/89 and the diagnostic value of intimin typing
.
Exp Biol Med (Maywood)
 
228
:
370
376
.
Kaper
JB
(
1998
)
EPEC delivers the goods
.
Trends Microbiol
 
6
:
169
172
; discussion 172–163.
Kenny
B
DeVinney
R
Stein
M
Reinscheid
DJ
Frey
EA
Finlay
BB
(
1997
)
Enteropathogenic E. coli (EPEC) transfers its receptor for intimate adherence into mammalian cells
.
Cell
 
91
:
511
520
.
Lazzaro
BP
Sceurman
BK
Carney
SL
Clark
AG
(
2002
)
fRFLP and fAFLP: medium-throughput genotyping by fluorescently post-labeling restriction digestion
.
Biotechniques
 
33
:
539
540
, 542, 545–546.
Luo
Y
Frey
EA
Pfuetzner
RA
Creagh
AL
Knoechel
DG
Haynes
CA
Finlay
BB
Strynadka
NC
(
2000
)
Crystal structure of enteropathogenic Escherichia coli intimin-receptor complex
.
Nature
 
405
:
1073
1077
.
McDaniel
TK
Kaper
JB
(
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
.
McDaniel
TK
Jarvis
KG
Donnenberg
MS
Kaper
JB
(
1995
)
A genetic locus of enterocyte effacement conserved among diverse enterobacterial pathogens
.
Proc Natl Acad Sci USA
 
92
:
1664
1668
.
McGraw
EA
Li
J
Selander
RK
Whittam
TS
(
1999
)
Molecular evolution and mosaic structure of α, β, and γ intimins of pathogenic Escherichia coli
.
Mol Biol Evol
 
16
:
12
22
.
Moyenuddin
M
Wachsmuth
IK
Moseley
SL
Bopp
CA
Blake
PA
(
1989
)
Serotype, antimicrobial resistance, and adherence properties of Escherichia coli strains associated with outbreaks of diarrheal illness in children in the United States
.
J Clin Microbiol
 
27
:
2234
2239
.
Orden
JA
Yuste
M
Cid
D
Piacesi
T
Martinez
S
Ruiz-Santa-Quiteria
JA
De la Fuente
R
(
2003
)
Typing of the eae and espB genes of attaching and effacing Escherichia coli isolates from ruminants
.
Vet Microbiol
 
96
:
203
215
.
Orskov
F
Whittam
TS
Cravioto
A
Orskov
I
(
1990
)
Clonal relationships among classic enteropathogenic Escherichia coli (EPEC) belong to different O groups
.
J Infect Dis
 
162
:
76
81
.
Oswald
E
Schmidt
H
Morabito
S
Karch
H
Marches
O
Caprioli
A
(
2000
)
Typing of intimin genes in human and animal enterohemorrhagic and enteropathogenic Escherichia coli: characterization of a new intimin variant
.
Infect Immun
 
68
:
64
71
.
Perna
NT
Mayhew
GF
Posfai
G
Elliott
S
Donnenberg
MS
Kaper
JB
Blattner
FR
(
1998
)
Molecular evolution of a pathogenicity island from enterohemorrhagic Escherichia coli O157:H7
.
Infect Immun
 
66
:
3810
3817
.
Ramachandran
V
Brett
K
Hornitzky
MA
Dowton
M
Bettelheim
KA
Walker
MJ
Djordjevic
SP
(
2003
)
Distribution of intimin subtypes among Escherichia coli isolates from ruminant and human sources
.
J Clin Microbiol
 
41
:
5022
5032
.
Reid
SD
Betting
DJ
Whittam
TS
(
1999
)
Molecular detection and identification of intimin alleles in pathogenic Escherichia coli by multiplex PCR
.
J Clin Microbiol
 
37
:
2719
2722
.
Robins-Browne
RM
Yam
WC
O'Gorman
LE
Bettelheim
KA
(
1993
)
Examination of archetypal strains of enteropathogenic Escherichia coli for properties associated with bacterial virulence
.
J Med Microbiol
 
38
:
222
226
.
Rodrigues
J
Scaletsky
IC
Campos
LC
Gomes
TA
Whittam
TS
Trabulsi
LR
(
1996
)
Clonal structure and virulence factors in strains of Escherichia coli of the classic serogroup O55
.
Infect Immun
 
64
:
2680
2686
.
Schmidt
H
Russmann
H
Karch
H
(
1993
)
Virulence determinants in nontoxinogenic Escherichia coli O157 strains that cause infantile diarrhea
.
Infect Immun
 
61
:
4894
4898
.
Sinclair
JF
Dean-Nystrom
EA
O'Brien
AD
(
2006
)
The established intimin receptor Tir and the putative eucaryotic intimin receptors nucleolin and β1 integrin localize at or near the site of enterohemorrhagic Escherichia coli O157:H7 adherence to enterocytes in vivo
.
Infect Immun
 
74
:
1255
1265
.
Sperandio
V
Kaper
JB
Bortolini
MR
Neves
BC
Keller
R
Trabulsi
LR
(
1998
)
Characterization of the locus of enterocyte effacement (LEE) in different enteropathogenic Escherichia coli (EPEC) and Shiga-toxin producing Escherichia coli (STEC) serotypes
.
FEMS Microbiol Lett
 
164
:
133
139
.
Tarr
CL
Whittam
TS
(
2002
)
Molecular evolution of the intimin gene in O111 clones of pathogenic Escherichia coli
.
J Bacteriol
 
184
:
479
487
.
Tarr
CL
Large
TM
Moeller
CL
Lacher
DW
Tarr
PI
Acheson
DW
Whittam
TS
(
2002
)
Molecular characterization of a serotype O121:H19 clone, a distinct Shiga toxin-producing clone of pathogenic Escherichia coli
.
Infect Immun
 
70
:
6853
6859
.
Unkmeir
A
Schmidt
H
(
2000
)
Structural analysis of phage-borne stx genes and their flanking sequences in shiga toxin-producing Escherichia coli and Shigella dysenteriae type 1 strains
.
Infect Immun
 
68
:
4856
4864
.
Valentiner-Branth
P
Steinsland
H
Fischer
TK
Perch
M
Scheutz
F
Dias
F
Aaby
P
Molbak
K
Sommerfelt
H
(
2003
)
Cohort study of Guinean children: incidence, pathogenicity, conferred protection, and attributable risk for enteropathogens during the first 2 years of life
.
J Clin Microbiol
 
41
:
4238
4245
.
Wales
AD
Woodward
MJ
Pearson
GR
(
2005
)
Attaching-effacing bacteria in animals
.
J Comp Pathol
 
132
:
1
26
.
Wieler
LH
McDaniel
TK
Whittam
TS
Kaper
JB
(
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
.
Zhang
WL
Kohler
B
Oswald
E
Beutin
L
Karch
H
Morabito
S
Caprioli
A
Suerbaum
S
Schmidt
H
(
2002
)
Genetic diversity of intimin genes of attaching and effacing Escherichia coli strains
.
J Clin Microbiol
 
40
:
4486
4492
.