Genomic diversity of pathogenic Escherichia coli of the EHEC 2 clonal complex

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

The global population structure and evolutionary history of the acquisition of major virulence factor-encoding genetic elements in Shiga toxin-producing Escherichia coli O121:H19

Shiga toxin (Stx)-producing Escherichia coli (STEC) are foodborne pathogens causing serious diseases, such as haemorrhagic colitis and haemolytic uraemic syndrome. Although O157:H7 STEC strains have been the most prevalent, incidences of STEC infections by several other serotypes have recently increased. O121:H19 STEC is one of these major non-O157 STECs, but systematic whole genome sequence (WGS) analyses have not yet been conducted on this STEC. Here, we performed a global WGS analysis of 638 O121:H19 strains, including 143 sequenced in this study, and a detailed comparison of 11 complete genomes, including four obtained in this study. By serotype-wide WGS analysis, we found that O121:H19 strains were divided into four lineages, including major and second major lineages (named L1 and L3, respectively), and that the locus of enterocyte effacement (LEE) encoding a type III secretion system (T3SS) was acquired by the common ancestor of O121:H19. Analyses of 11 complete genomes belonging to L1 or L3 revealed remarkable interlineage differences in the prophage pool and prophage-encoded T3SS effector repertoire, independent acquisition of virulence plasmids by the two lineages, and high conservation in the prophage repertoire, including that for Stx2a phages in lineage L1. Further sequence determination of complete Stx2a phage genomes of 49 strains confirmed that Stx2a phages in lineage L1 are highly conserved short-tailed phages, while those in lineage L3 are long-tailed lambda-like phages with notable genomic diversity, suggesting that an Stx2a phage was acquired by the common ancestor of L1 and has been stably maintained. Consistent with these genomic features of Stx2a phages, most lineage L1 strains produced much higher levels of Stx2a than lineage L3 strains. Altogether, this study provides a global phylogenetic overview of O121:H19 STEC and shows the interlineage genomic differences and the highly conserved genomic features of the major lineage within this serotype of STEC.

Population structure of Escherichia coli O26 : H11 with recent and repeated stx2 acquisition in multiple lineages

A key virulence factor of enterohaemorrhagic Escherichia coli (EHEC) is the bacteriophage-encoded Shiga toxin (Stx). Stxs are classified into two types, Stx1 and Stx2, and Stx2-producing strains are thought to cause more severe infections than strains producing only Stx1. Although O26 : H11 is the second most prevalent EHEC following O157 : H7, the majority of O26 : H11 strains produce Stx1 alone. However, Stx2-producing O26 strains have increasingly been detected worldwide. Through a large-scale genome analysis, we present a global phylogenetic overview and evolutionary timescale for E. coli O26 : H11. The origin of O26 has been estimated to be 415 years ago. Sequence type 21C1 (ST21C1), one of the two sublineages of ST21, the most predominant O26 : H11 lineage worldwide, emerged 213 years ago from one of the three ST29 sublineages (ST29C2). The other ST21 lineage (ST21C2) emerged 95 years ago from ST21C1. Increases in population size occurred in the late 20th century for all of the O26 lineages, but most remarkably for ST21C2. Analysis of the distribution of stx2-positive strains revealed the recent and repeated acquisition of the stx2 gene in multiple lineages of O26, both in ST21 and ST29. Other major EHEC virulence genes, such as type III secretion system effector genes and plasmid-encoded virulence genes, were well conserved in ST21 compared to ST29. In addition, more antimicrobial-resistance genes have accumulated in the ST21C1 lineage. Although current attention is focused on several highly virulent ST29 clones that have acquired the stx2 gene, there is also a considerable risk that the ST21 lineage could yield highly virulent clones.

Intestinal Microbial Community Dynamics of White-Tailed Deer (Odocoileus virginianus) in an Agroecosystem

The intestinal microbiota has important functions that contribute to host health. The compositional dynamics of microbial communities are affected by many factors, including diet and presence of pathogens. In contrast to humans and domestic mammals, the composition and seasonal dynamics of intestinal microbiota of wildlife species remain comparatively understudied. White-tailed deer (Odocoileus virginianus) is an ecologically and economically important wildlife species that inhabits agricultural ecosystems and is known to be a reservoir of enteric pathogens. Nevertheless, there is a lack of knowledge of white-tailed deer intestinal microbiota diversity and taxonomic composition. This study's first objective was to characterize and compare the intestinal microbiota of 66 fecal samples from white-tailed deer collected during two sampling periods (March and June) using 16S rDNA pyrosequencing. Associations between community diversity and composition and factors including season, sex, host genetic relatedness, and spatial location were quantified. Results revealed that white-tailed deer intestinal microbiota was predominantly comprised of phyla Firmicutes and Proteobacteria, whose relative frequencies varied significantly between sampling periods. The second objective was to examine the associations between the presence of Escherichia coli and Salmonella, and microbiota composition and diversity. Results indicated that relative abundance of some microbial taxa varied when a pathogen was present. This study provides insights into microbial compositional dynamics of a wildlife species inhabiting coupled natural and agricultural landscapes. Data focus attention on the high prevalence of Proteobacteria particularly during the summer and highlight the need for future research regarding the role of white-tailed deer as a natural pathogen reservoir in agroecosystems.

Genetic Relatedness and Novel Sequence Types of Non-O157 Shiga Toxin-Producing Escherichia coli Strains Isolated in Argentina

Shiga toxin-producing Escherichia coli (STEC) is a foodborne pathogen responsible for severe disease in humans such as hemolytic uremic syndrome (HUS) and cattle, the principal reservoir. Identification of the clones/lineages is important as several characteristics, among them propensity to cause disease varies with STEC phylogenetic origin. At present, we do not know what STEC clones, especially of non-O157:H7, are circulating in Argentina. To fill this knowledge gap we assessed the genetic diversity of STEC strains isolated in Argentina from various sources, mostly cattle and food, using multilocus sequence typing (MLST). Our objectives were to determine the phylogenetic relationships among strains and to compare them with strains from different geographic origins, especially with those from clinical human cases, in order to evaluate their potential health risk. A total of 59 STEC isolates from 41 serotypes were characterized by MLST. Analysis using EcMLST database identified 38 sequence types (ST), 17 (45%) of which were new STs detected in 18 serotypes. Fifteen out of 38 STs identified were grouped into 11 clonal groups (CGs) and, 23 not grouped in any of the defined CGs. Different STs were found in the same serotype. Results highlighted a high degree of phylogenetic heterogeneity among Argentinean strains and they showed that several cattle and food isolates belonged to the same STs that are commonly associated with clinical human cases in several geographical areas. STEC is a significant public health concern. Argentina has the highest incidence of HUS in the world and this study provides the first data about which STEC clones are circulating. Data showed that most of them might pose a serious zoonotic risk and this information is important for developing public health initiatives. However, the actual potential risk will be defined by the virulence profiles, which may differ among isolates belonging to the same ST.

Whole Genome Sequencing for Genomics-Guided Investigations of Escherichia coli O157:H7 Outbreaks

Multi isolate whole genome sequencing (WGS) and typing for outbreak investigations has become a reality in the post-genomics era. We applied this technology to strains from Escherichia coli O157:H7 outbreaks. These include isolates from seven North America outbreaks, as well as multiple isolates from the same patient and from different infected individuals in the same household. Customized high-resolution bioinformatics sequence typing strategies were developed to assess the core genome and mobilome plasticity. Sequence typing was performed using an in-house single nucleotide polymorphism (SNP) discovery and validation pipeline. Discriminatory power becomes of particular importance for the investigation of isolates from outbreaks in which macrogenomic techniques such as pulse-field gel electrophoresis or multiple locus variable number tandem repeat analysis do not differentiate closely related organisms. We also characterized differences in the phage inventory, allowing us to identify plasticity among outbreak strains that is not detectable at the core genome level. Our comprehensive analysis of the mobilome identified multiple plasmids that have not previously been associated with this lineage. Applied phylogenomics approaches provide strong molecular evidence for exceptionally little heterogeneity of strains within outbreaks and demonstrate the value of intra-cluster comparisons, rather than basing the analysis on archetypal reference strains. Next generation sequencing and whole genome typing strategies provide the technological foundation for genomic epidemiology outbreak investigation utilizing its significantly higher sample throughput, cost efficiency, and phylogenetic relatedness accuracy. These phylogenomics approaches have major public health relevance in translating information from the sequence-based survey to support timely and informed countermeasures. Polymorphisms identified in this work offer robust phylogenetic signals that index both short- and long-term evolution and can complement currently employed typing schemes for outbreak ex- and inclusion, diagnostics, surveillance, and forensic studies.

Characterization of Shiga Toxin Subtypes and Virulence Genes in Porcine Shiga Toxin-Producing Escherichia coli

Similar to ruminants, swine have been shown to be a reservoir for Shiga toxin-producing Escherichia coli (STEC), and pork products have been linked with outbreaks associated with STEC O157 and O111:H-. STEC strains, isolated in a previous study from fecal samples of late-finisher pigs, belonged to a total of 56 serotypes, including O15:H27, O91:H14, and other serogroups previously associated with human illness. The isolates were tested by polymerase chain reaction (PCR) and a high-throughput real-time PCR system to determine the Shiga toxin (Stx) subtype and virulence-associated and putative virulence-associated genes they carried. Select STEC strains were further analyzed using a Minimal Signature E. coli Array Strip. As expected, stx_2e (81%) was the most common Stx variant, followed by stx_1a (14%), stx_2d (3%), and stx_1c (1%). The STEC serogroups that carried stx_2d were O15:H27, O159:H16 and O159:H-. Similar to stx_2a and stx_2c, the stx_2d variant is associated with development of hemorrhagic colitis and hemolytic uremic syndrome, and reports on the presence of this variant in STEC strains isolated from swine are lacking. Moreover, the genes encoding heat stable toxin (estIa) and enteroaggregative E. coli heat stable enterotoxin-1 (astA) were commonly found in 50 and 44% of isolates, respectively. The hemolysin genes, hlyA and ehxA, were both detected in 7% of the swine STEC strains. Although the eae gene was not found, other genes involved in host cell adhesion, including lpfA_O113 and paa were detected in more than 50% of swine STEC strains, and a number of strains also carried iha, lpfA_O26, lpfA_O157, fedA, orfA, and orfB. The present work provides new insights on the distribution of virulence factors among swine STEC strains and shows that swine may carry Stx1a-, Stx2e-, or Stx2d-producing E. coli with virulence gene profiles associated with human infections.

Detection of Shiga Toxin-Producing Escherichia coli from Nonhuman Sources and Strain Typing

Shiga toxin-producing Escherichia coli (STEC) strains are commonly found in the intestine of ruminant species of wild and domestic animals. Excretion of STEC with animal feces results in a broad contamination of food and the environment. Humans get infected with STEC through ingestion of contaminated food, by contact with the environment, and from STEC-excreting animals and humans. STEC strains can behave as human pathogens, and some of them, called enterohemorrhagic E. coli (EHEC), may cause hemorrhagic colitis (HC) and hemolytic-uremic syndrome (HUS). Because of the diversity of STEC types, detection strategies for STEC and EHEC are based on the identification of Shiga toxins or the underlying genes. Cultural enrichment of STEC from test samples is needed for identification, and different protocols were developed for this purpose. Multiplex real-time PCR protocols (ISO/CEN TS13136 and USDA/FSIS MLG5B.01) have been developed to specifically identify EHEC by targeting the LEE (locus of enterocyte effacement)-encoded eae gene and genes for EHEC-associated O groups. The employment of more genetic markers (nle and CRISPR) is a future challenge for better identification of EHEC from any kinds of samples. The isolation of STEC or EHEC from a sample is required for confirmation, and different cultivation protocols and media for this purpose have been developed. Most STEC strains present in food, animals, and the environment are eae negative, but some of these strains can cause HC and HUS in humans as well. Phenotypic assays and molecular tools for typing EHEC and STEC strains are used to detect and characterize human pathogenic strains among members of the STEC group.

Genetic Diversity and Pathogenic Potential of Attaching and Effacing Escherichia coli O26:H11 Strains Recovered from Bovine Feces in the United States

Escherichia coli O26 has been identified as the most common non-O157 Shiga toxin-producing E. coli (STEC) serogroup to cause human illnesses in the United States and has been implicated in outbreaks around the world. E. coli has high genomic plasticity, which facilitates the loss or acquisition of virulence genes. Attaching and effacing E. coli (AEEC) O26 strains have frequently been isolated from bovine feces, and there is a need to better characterize the relatedness of these strains to defined molecular pathotypes and to describe the extent of their genetic diversity. High-throughput real-time PCR was used to screen 178 E. coli O26 isolates from a single U.S. cattle feedlot, collected from May to July 2011, for the presence or absence of 25 O26 serogroup-specific and virulence-associated markers. The selected markers were capable of distinguishing these strains into molecularly defined groups (yielding 18 unique marker combinations). Analysis of the clustered regularly interspaced short palindromic repeat 1 (CRISPR1) and CRISPR2a loci further discriminated isolates into 24 CRISPR types. The combination of molecular markers and CRISPR typing provided 20.8% diversity. The recent CRISPR PCR target SP_O26-E, which was previously identified only in stx2-positive O26:H11 human clinical strains, was identified in 96.4% (161/167 [95% confidence interval, 99.2 to 93.6%]) of the stx-negative AEEC O26:H11 bovine fecal strains. This supports that these stx-negative strains may have previously contained a prophage carrying stx or could acquire this prophage, thus possibly giving them the potential to become pathogenic to humans. These results show that investigation of specific genetic markers may further elucidate our understanding of the genetic diversity of AEEC O26 strains in bovine feces.

Diverse virulence gene content of Shiga toxin-producing Escherichia coli from finishing swine

Shiga toxin-producing Escherichia coli (STEC) infections are a critical public health concern because they can cause severe clinical outcomes, such as hemolytic uremic syndrome, in humans. Determining the presence or absence of virulence genes is essential in assessing the potential pathogenicity of STEC strains. Currently, there is limited information about the virulence genes carried by swine STEC strains; therefore, this study was conducted to examine the presence and absence of 69 virulence genes in STEC strains recovered previously from finishing swine in a longitudinal study. A subset of STEC strains was analyzed by pulsed-field gel electrophoresis (PFGE) to examine their genetic relatedness. Swine STEC strains (n = 150) were analyzed by the use of a high-throughput real-time PCR array system, which included 69 virulence gene targets. Three major pathotypes consisted of 16 different combinations of virulence gene profiles, and serotypes were determined in the swine STEC strains. The majority of the swine STEC strains (n = 120) belonged to serotype O59:H21 and carried the same virulence gene profile, which consisted of 9 virulence genes: stx2e, iha, ecs1763, lpfAO113, estIa (STa), ehaA, paa, terE, and ureD. The eae, nleF, and nleH1-2 genes were detected in one swine STEC strain (O49:H21). Other genes encoding adhesins, including iha, were identified (n = 149). The PFGE results demonstrated that swine STEC strains from pigs raised in the same finishing barn were closely related. Our results revealed diverse virulence gene contents among the members of the swine STEC population and enhance understanding of the dynamics of transmission of STEC strains among pigs housed in the same barn.

Taxonomy Meets Public Health: The Case of Shiga Toxin-Producing Escherichia coli

To help assess the clinical and public health risks associated with different Shiga toxin-producing Escherichia coli (STEC) strains, an empirical classification scheme was used to classify STEC into five “seropathotypes” (seropathotype A [high risk] to seropathotypes D and E [minimal risk]). This definition is of considerable value in cases of human infection but is also problematic because not all STEC infections are fully characterized and coupled to reliable clinical information. Outbreaks with emerging hybrid strains continuously challenge our understanding of virulence potential and may result in incorrect classification of specific pathotypes; an example is the hybrid strain that caused the 2011 outbreak in Germany, STEC/EAggEC O104:H4, which may deserve an alternative seropathotype designation. The integration of mobile virulence factors in the stepwise and parallel evolution of pathogenic lineages of STEC collides with the requirements of a good taxonomy, which separates elements of each group into subgroups that are mutually exclusive, unambiguous, and, together, include all possibilities. The concept of (sero)-pathotypes is therefore challenged, and the need to identify factors of STEC that absolutely predict the potential to cause human disease is obvious. Because the definition of hemolytic-uremic syndrome (HUS) is distinct, a basic and primary definition of HUS-associated E. coli (HUSEC) for first-line public health action is proposed: stx2 in a background of an eae- or aggR -positive E. coli followed by a second-line subtyping of stx genes that refines the definition of HUSEC to include only stx2a and stx2d . All other STEC strains are considered “low-risk” STEC.

Comparative genomics of enterohemorrhagic Escherichia coli O145:H28 demonstrates a common evolutionary lineage with Escherichia coli O157:H7

Background

Although serotype O157:H7 is the predominant enterohemorrhagic Escherichia coli (EHEC), outbreaks of non-O157 EHEC that cause severe foodborne illness, including hemolytic uremic syndrome have increased worldwide. In fact, non-O157 serotypes are now estimated to cause over half of all the Shiga toxin-producing Escherichia coli (STEC) cases, and outbreaks of non-O157 EHEC infections are frequently associated with serotypes O26, O45, O103, O111, O121, and O145. Currently, there are no complete genomes for O145 in public databases.

Results

We determined the complete genome sequences of two O145 strains (EcO145), one linked to a US lettuce-associated outbreak (RM13514) and one to a Belgium ice-cream-associated outbreak (RM13516). Both strains contain one chromosome and two large plasmids, with genome sizes of 5,737,294 bp for RM13514 and 5,559,008 bp for RM13516. Comparative analysis of the two EcO145 genomes revealed a large core (5,173 genes) and a considerable amount of strain-specific genes. Additionally, the two EcO145 genomes display distinct chromosomal architecture, virulence gene profile, phylogenetic origin of Stx2a prophage, and methylation profile (methylome). Comparative analysis of EcO145 genomes to other completely sequenced STEC and other E. coli and Shigella genomes revealed that, unlike any other known non-O157 EHEC strain, EcO145 ascended from a common lineage with EcO157/EcO55. This evolutionary relationship was further supported by the pangenome analysis of the 10 EHEC str ains. Of the 4,192 EHEC core genes, EcO145 shares more genes with EcO157 than with the any other non-O157 EHEC strains.

Conclusions

Our data provide evidence that EcO145 and EcO157 evolved from a common lineage, but ultimately each serotype evolves via a lineage-independent nature to EHEC by acquisition of the core set of EHEC virulence factors, including the genes encoding Shiga toxin and the large virulence plasmid. The large variation between the two EcO145 genomes suggests a distinctive evolutionary path between the two outbreak strains. The distinct methylome between the two EcO145 strains is likely due to the presence of a BsuBI/PstI methyltransferase gene cassette in the Stx2a prophage of the strain RM13514, suggesting a role of horizontal gene transfer-mediated epigenetic alteration in the evolution of individual EHEC strains.

Genomic diversity of Escherichia isolates from diverse habitats

Our understanding of the Escherichia genus is heavily biased toward pathogenic or commensal isolates from human or animal hosts. Recent studies have recovered Escherichia isolates that persist, and even grow, outside these hosts. Although the environmental isolates are typically phylogenetically distinct, they are highly related to and phenotypically indistinguishable from their human counterparts, including for the coliform test. To gain insights into the genomic diversity of Escherichia isolates from diverse habitats, including freshwater, soil, animal, and human sources, we carried out comparative DNA-DNA hybridizations using a multi-genome E. coli DNA microarray. The microarray was validated based on hybridizations with selected strains whose genome sequences were available and used to assess the frequency of microarray false positive and negative signals. Our results showed that human fecal isolates share two sets of genes (n>90) that are rarely found among environmental isolates, including genes presumably important for evading host immune mechanisms (e.g., a multi-drug transporter for acids and antimicrobials) and adhering to epithelial cells (e.g., hemolysin E and fimbrial-like adhesin protein). These results imply that environmental isolates are characterized by decreased ability to colonize host cells relative to human isolates. Our study also provides gene markers that can distinguish human isolates from those of warm-blooded animal and environmental origins, and thus can be used to more reliably assess fecal contamination in natural ecosystems.

Tarr C, Gerner-Smidt P, Datta AR. Genomic characterization of Listeria monocytogenes strains involved in a multistate listeriosis outbreak associated with cantaloupe in US

A multistate listeriosis outbreak associated with cantaloupe consumption was reported in the United States in September, 2011. The outbreak investigation recorded a total of 146 invasive illnesses, 30 deaths and one miscarriage. Subtyping of the outbreak associated clinical, food and environmental isolates revealed two serotypes (1/2a and 1/2b) and four pulsed-field gel electrophoresis two-enzyme pattern combinations I, II, III, and IV, including one rarely seen before this outbreak. A DNA-microarray, Listeria GeneChip®, developed by FDA from 24 Listeria monocytogenes genome sequences, was used to further characterize a representative sample of the outbreak isolates. The microarray data (in the form of present or absent calls of specific DNA sequences) separated the isolates into two distinct groups as per their serotypes. The gene content of the outbreak-associated isolates was distinct from that of the previously-reported outbreak strains belonging to the same serotypes. Although the 1/2b outbreak associated isolates are closely related to each other, the 1/2a isolates could be further divided into two distinct genomic groups, one represented by pattern combination I strains and the other represented by highly similar pattern combinations III and IV strains. Gene content analysis of these groups revealed unique genomic sequences associated with these two 1/2a genovars. This work underscores the utility of multiple approaches, such as serotyping, PFGE and DNA microarray analysis to characterize the composition of complex polyclonal listeriosis outbreaks.

High density microarray analysis reveals new insights into genetic footprints of Listeria monocytogenes strains involved in listeriosis outbreaks

Listeria monocytogenes, a foodborne bacterial pathogen, causes invasive and febrile gastroenteritis forms of listeriosis in humans. Both invasive and febrile gastroenteritis listeriosis is caused mostly by serotypes 1/2a, 1/2b and 4b strains. The outbreak strains of serotype 1/2a and 4b could be further classified into several epidemic clones but the genetic bases for the diverse pathophysiology have been unsuccessful. DNA microarray provides an important tool to scan the entire genome for genetic signatures that may distinguish the L. monocytogenes strains belonging to different outbreaks. We have designed a pan-genomic microarray chip (Listeria GeneChip) containing sequences from 24 L. monocytogenes strains. The chip was designed to identify the presence/absence of genomic sequences, analyze transcription profiles and identify SNPs. Analysis of the genomic profiles of 38 outbreak strains representing 1/2a, 1/2b and 4b serotypes, revealed that the strains formed distinct genetic clusters adhering to their serotypes and epidemic clone types. Although serologically 1/2a and 1/b strains share common antigenic markers microarray analysis revealed that 1/2a strains are further apart from the closely related 1/2b and 4b strains. Within any given serotype and epidemic clone type the febrile gastroenteritis and invasive strains can be further distinguished based on several genetic markers including large numbers of phage genome, and intergenic sequences. Our results showed that the microarray-based data can be an important tool in characterization of L. monocytogenes strains involved in both invasive and gastroenteritis outbreaks. The results for the first time showed that the serotypes and epidemic clones are based on extensive pan-genomic variability and the 1/2b and 4bstrains are more closely related to each other than the 1/2a strains. The data also supported the hypothesis that the strains causing these two diverse outbreaks are genotypically different and this finding might be important in understanding the pathophysiology of this organism.

Virulence gene profiling of enterohemorrhagic (EHEC) and enteropathogenic (EPEC) Escherichia coli strains: a basis for molecular risk assessment of typical and atypical EPEC strains

Background

Enterohaemorrhagic E. coli (EHEC) can cause severe disease such as bloody diarrhoea and haemolytic uraemic syndrome in humans. Besides production of Shiga toxins, the presence of LEE (eae-gene) and non-LEE (nle) encoded effector genes harboured on O-islands OI-122, OI-71 and OI-57 is associated with EHEC virulence and their frequency in outbreaks. Genes encoded by the EHEC-plasmid are putative virulence markers of EHEC. EHEC-plasmids, LEE and non-LEE effector genes have also been detected in some strains of enteropathogenic E. coli (EPEC). The objective of this study was to analyze the relationship between EHEC and EPEC for virulence genes encoded by genomic O-islands and by the EHEC-plasmids.

Results

Nle genes ent/espL2, nleB and nleE (OI-122), nleA, nleF and nleH1-2 (OI-71), nleG5-2 and nleG6-2 (OI-57), espK (CP-933N) and the EHEC-plasmid encoded genes ehxA, espP, etpD and katP were searched in 73 typical and in 235 atypical enteropathogenic E. coli (EPEC) strains. Typical and atypical EPEC each fall into two clusters. Cluster 1 typical (n = 46) and atypical (n = 129) EPEC strains were characterized by the presence of OI-122 encoded genes and grouped together with 64 investigated EHEC strains. Cluster 2 typical (n = 27) and atypical (n = 106) strains grouped together with 52 LEE-negative, Shiga toxin-producing E. coli (STEC) and with 21 apathogenic E. coli strains. Typical EPEC Cluster 1 strains belonged to serotypes frequently involved in severe illness and outbreaks in children (O111:H2, O114:H2, O55:H6, O127:H6 and O142:H6). Atypical EPEC Cluster 1 strains were characterized by serotypes related to EHEC (O26:H11, O55:H7, O145:H28, O103:H2 and O103:H25).

Conclusion

The OI-122 encoded nleB gene was found to be most closely associated with Cluster 1 strains and may serve as a diagnostic tool for the identification of virulent EHEC and EPEC seropathotypes. OI-71 encoded genes nleA, nleF and nleH1-2 are less associated with Cluster 1 strains. EHEC-plasmid, OI-57 and CP-933 associated genes showed only weak similarities with virulent Cluster 1 EHEC and EPEC strains.