Genotypic characterization and prevalence of virulence factors among Canadian Escherichia coli O157:H7 strains

Commensal Escherichia coli Strains of Bovine Origin Competitively Mitigated Escherichia coli O157:H7 in a Gnotobiotic Murine Intestinal Colonization Model with or without Physiological Stress

Cattle are a primary reservoir of enterohemorrhagic Escherichia coli (EHEC) O157:H7. Currently, there are no effective methods of eliminating this important zoonotic pathogen from cattle, and colonization resistance in relation to EHEC O157:H7 in cattle is poorly understood. We developed a gnotobiotic EHEC O157:H7 murine model to examine aspects of the cattle pathogen-microbiota interaction, and to investigate competitive suppression of EHEC O157:H7 by 18 phylogenetically distinct commensal E. coli strains of bovine origin. As stress has been suggested to influence enteric colonization by EHEC O157:H7 in cattle, corticosterone administration (±) to incite a physiological stress response was included as an experimental variable. Colonization of the intestinal tract (IT) of mice by the bovine EHEC O157:H7 strain, FRIK-2001, mimicked characteristics of bovine IT colonization. In this regard, FRIK-2001 successfully colonized the IT and temporally incited minimal impacts on the host relative to other EHEC O157:H7 strains, including on the renal metabolome. The presence of the commensal E. coli strains decreased EHEC O157:H7 densities in the cecum, proximal colon, and distal colon. Moreover, histopathologic changes and inflammation markers were reduced in the distal colon of mice inoculated with commensal E. coli strains (both propagated separately and communally). Although stress induction affected the behavior of mice, it did not influence EHEC O157:H7 densities or disease. These findings support the use of a gnotobiotic murine model of enteric bovine EHEC O157:H7 colonization to better understand pathogen-host-microbiota interactions toward the development of effective on-farm mitigations for EHEC O157:H7 in cattle, including the identification of bacteria capable of competitively colonizing the IT.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Genome Sequences of 14 Escherichia coli O157:H7 Strains Isolated before and during the Time Frame of the 2018 Multistate Outbreak Associated with Romaine Lettuce

Several outbreaks of Escherichia coli O157:H7 associated with contaminated leafy green vegetables have been documented. Here, we report the draft genome sequences of 14 strains isolated from human patients in the state of Wisconsin during a multistate outbreak in early 2018 that was linked to consumption of romaine lettuce.

Epidemiology of Shiga Toxin-Producing Escherichia coli O157 in the Province of Alberta, Canada, 2009-2016

Shiga toxin-producing Escherichia coli (STEC) infections are the product of the interaction between bacteria, phages, animals, humans, and the environment. In the late 1980s, Alberta had one of the highest incidences of STEC infections in North America. Herein, we revisit and contextualize the epidemiology of STEC O157 human infections in Alberta for the period 2009–2016. STEC O157 infections were concentrated in large urban centers, but also in rural areas with high cattle density. Hospitalization was often required when the Shiga toxin genotype stx2a stx2c was involved, however, only those aged 60 years or older and infection during spring months (April to June) independently predicted that need. Since the late 1980s, the rate of STEC O157-associated hemolytic uremic syndrome (HUS) in Alberta has remained unchanged at 5.1%, despite a marked drop in the overall incidence of the infection. While Shiga toxin genotypes stx1a stx2c and stx2a stx2c seemed associated with HUS, only those aged under 10 years and infection during spring months were independently predictive of that complication. The complexity of the current epidemiology of STEC O157 in Alberta highlights the need for a One Health approach for further progress to be made in mitigating STEC morbidity.

Comparative genomics reveals structural and functional features specific to the genome of a foodborne Escherichia coli O157:H7

Background

Escherichia coli O157:H7 (O157) has been linked to numerous foodborne disease outbreaks. The ability to rapidly sequence and analyze genomes is important for understanding epidemiology, virulence, survival, and evolution of outbreak strains. In the current study, we performed comparative genomics to determine structural and functional features of the genome of a foodborne O157 isolate NADC 6564 and infer its evolutionary relationship to other O157 strains.

Results

The chromosome of NADC 6564 contained 5466 kb compared to reference strains Sakai (5498 kb) and EDL933 (5547 kb) and shared 41 of its 43 Linear Conserved Blocks (LCB) with the reference strains. However, 18 of 41 LCB had inverse orientation in NADC 6564 compared to the reference strains. NADC 6564 shared 18 of 19 bacteriophages with reference strains except that the chromosomal positioning of some of the phages differed among these strains. The additional phage (P19) of NADC 6564 was located on a 39-kb insertion element (IE) encoding several hypothetical proteins, an integrase, transposases, transcriptional regulators, an adhesin, and a phosphoethanolamine transferase (PEA). The complete homologs of the 39-kb IE were found in E. coli PCN061 of porcine origin. The IE-encoded PEA showed low homology (32–33%) to four other PEA in NADC 6564 and PEA linked to mobilizable colistin resistance in E. coli but was highly homologous (95%) to a PEA of uropathogenic, avian pathogenic, and enteroaggregative E. coli. NADC 6564 showed slightly higher minimum inhibitory concentration of colistin compared to the reference strains. The 39-kb IE also contained dndBCDE and dptFGH operons encoding DNA S-modification and a restriction pathway, linked to oxidative stress tolerance and self-defense against foreign DNA, respectively. Evolutionary tree analysis grouped NADC 6564 with lineage I O157 strains.

Conclusions

These results indicated that differential phage counts and different chromosomal positioning of many bacteriophages and genomic islands might have resulted in recombination events causing altered chromosomal organization in NADC 6564. Evolutionary analysis grouped NADC 6564 with lineage I strains and suggested its earlier divergence from these strains. The ability to perform S-DNA modification might affect tolerance of NADC 6564 to various stressors.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5568-6) contains supplementary material, which is available to authorized users.

Bacteriophages Contribute to the Spread of Antibiotic Resistance Genes among Foodborne Pathogens of the Enterobacteriaceae Family - A Review

Foodborne illnesses continue to have an economic impact on global health care systems. There is a growing concern regarding the increasing frequency of antibiotic resistance in foodborne bacterial pathogens and how such resistance may affect treatment outcomes. In an effort to better understand how to reduce the spread of resistance, many research studies have been conducted regarding the methods by which antibiotic resistance genes are mobilized and spread between bacteria. Transduction by bacteriophages (phages) is one of many horizontal gene transfer mechanisms, and recent findings have shown phage-mediated transduction to be a significant contributor to dissemination of antibiotic resistance genes. Here, we review the viability of transduction as a contributing factor to the dissemination of antibiotic resistance genes in foodborne pathogens of the Enterobacteriaceae family, including non-typhoidal Salmonella and Shiga toxin-producing Escherichia coli, as well as environmental factors that increase transduction of antibiotic resistance genes.

Diversity of Survival Patterns among Escherichia coli O157:H7 Genotypes Subjected to Food-Related Stress Conditions

The purpose of this study was to evaluate the resistance patterns to food-related stresses of Shiga toxin producing Escherichia coli O157:H7 strains belonging to specific genotypes. A total of 33 E. coli O157:H7 strains were exposed to seven different stress conditions acting as potential selective pressures affecting the transmission of E. coli O157:H7 to humans through the food chain. These stress conditions included cold, oxidative, osmotic, acid, heat, freeze-thaw, and starvation stresses. The genotypes used for comparison included lineage-specific polymorphism, Shiga-toxin-encoding bacteriophage insertion sites, clade type, tir (A255T) polymorphism, Shiga toxin 2 subtype, and antiterminator Q gene allele. Bacterial resistance to different stressors was calculated by determining D-values (times required for inactivation of 90% of the bacterial population), which were then subjected to univariate and multivariate analyses. In addition, a relative stress resistance value, integrating resistance values to all tested stressors, was calculated for each bacterial strain and allowed for a ranking-type classification of E. coli O157:H7 strains according to their environmental robustness. Lineage I/II strains were found to be significantly more resistant to acid, cold, and starvation stress than lineage II strains. Similarly, tir (255T) and clade 8 encoding strains were significantly more resistant to acid, heat, cold, and starvation stress than tir (255A) and non-clade 8 strains. Principal component analysis, which allows grouping of strains with similar stress survival characteristics, separated strains of lineage I and I/II from strains of lineage II, which in general showed reduced survival abilities. Results obtained suggest that lineage I/II, tir (255T), and clade 8 strains, which have been previously reported to be more frequently associated with human disease cases, have greater multiple stress resistance than strains of other genotypes. The results from this study provide a better insight into how selective pressures encountered through the food chain may play a role in the epidemiology of STEC O157:H7 through controlling the transmission of highly adapted strains to humans.

Genetic features of human and bovine Escherichia coli O157:H7 strains isolated in Argentina

Shiga toxin-producing Escherichia coli (STEC) are important food-borne pathogens associated with human diseases. In Argentina, O157:H7 is the dominant serotype in hemolytic uremic syndrome (HUS) cases. Previously, we have described the almost exclusive circulation of human E. coli O157 strains belonging to the hypervirulent clade 8 in Neuquén Province. The aim of the present study was to investigate, by a broad molecular characterization, if this particular distribution of E. coli O157 clades in Neuquén is similar to the situation in other regions of the country and if it may be originated in a similar profile in cattle, its main reservoir. Two-hundred and eighty O157 strains (54 bovine and 226 human) isolated between 2006 and 2008 in different regions of Argentina were studied. All strains harbored rfbO157, fliCH7, eae, and ehxA genes. The predominant genotype was stx2a/stx2c in human (76.1%) and bovine (55.5%) strains. All human isolates tested by Lineage-Specific Polymorphism Assay (LSPA-6), were lineage I/II; among bovine strains, 94.1% belonged to lineage I/II and 5.9% to lineage I. No LSPA-6 lineage II isolates were detected. Single nucleotide polymorphism (SNP) analysis has revealed the existence of nine clade phylogenetic groups. In our clinical strains collection, 87.6% belonged to the hypervirulent clade 8, and 12.4% were classified as clade 4/5. In bovine isolates, 59.3% strains were clade 8, 33.3% clade 4/5 and 7.4% clade 3. More than 80% of human strains showed the presence of 6 of the 7 virulence determinants described in the TW14359 O157 strain associated with the raw spinach outbreak in the U.S. in 2006. More than 80% of bovine strains showed the presence of 3 of these factors. The q933 allele, which has been related to high toxin production, was present in 98.2% of clinical strains and 75.9% of the bovine isolates. The molecular characterization of human STEC O157 strains allows us to conclude that the particular situation previously described for Neuquén Province, may actually be a characteristic of the whole country. These genetic features are quite similar to those observed in the bovine reservoir and may be derived from it. This data confirms that, unlike the rest of the world, in Argentina most of the STEC O157 strains present in cattle may cause human infections of varying severity and the marked virulence described for these strains may be related to the high incidence of HUS in our country.

Escherichia coli O157:H7 strains harbor at least three distinct sequence types of Shiga toxin 2a-converting phages

Background

Shiga toxin-producing Escherichia coli O157:H7 is a foodborne pathogen that causes severe human diseases including hemolytic uremic syndrome (HUS). The virulence factor that mediates HUS, Shiga toxin (Stx), is encoded within the genome of a lambdoid prophage. Although draft sequences are publicly available for a large number of E. coli O157:H7 strains, the high sequence similarity of stx-converting bacteriophages with other lambdoid prophages poses challenges to accurately assess the organization and plasticity among stx-converting phages due to assembly difficulties.

Methods

To further explore genome plasticity of stx-converting prophages, we enriched phage DNA from 45 ciprofloxacin-induced cultures for subsequent 454 pyrosequencing to facilitate assembly of the complete phage genomes. In total, 22 stx2a-converting phage genomes were closed.

Results

Comparison of the genomes distinguished nine distinct phage sequence types (PSTs) delineated by variation in obtained sequences, such as single nucleotide polymorphisms (SNPs) and insertion sequence element prevalence and location. These nine PSTs formed three distinct clusters, designated as PST1, PST2 and PST3. The PST2 cluster, identified in two clade 8 strains, was related to stx2a-converting phages previously identified in non-O157 Shiga-toxin producing E. coli (STEC) strains associated with a high incidence of HUS. The PST1 cluster contained phages related to those from E. coli O157:H7 strain Sakai (lineage I, clade 1), and PST3 contained a single phage that was distinct from the rest but most related to the phage from E. coli O157:H7 strain EC4115 (lineage I/II, clade 8). Five strains carried identical stx2a-converting phages (PST1-1) integrated at the same chromosomal locus, but these strains produced different levels of Stx2.

Conclusion

The stx2a-converting phages of E. coli O157:H7 can be categorized into at least three phage types. Diversification within a phage type is mainly driven by IS629 and by a small number of SNPs. Polymorphisms between phage genomes may help explain differences in Stx2a production between strains, however our data indicates that genes encoded external to the phage affect toxin production as well.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1934-1) contains supplementary material, which is available to authorized users.

Utility of Whole-Genome Sequencing of Escherichia coli O157 for Outbreak Detection and Epidemiological Surveillance

Detailed laboratory characterization of Escherichia coli O157 is essential to inform epidemiological investigations. This study assessed the utility of whole-genome sequencing (WGS) for outbreak detection and epidemiological surveillance of E. coli O157, and the data were used to identify discernible associations between genotypes and clinical outcomes. One hundred five E. coli O157 strains isolated over a 5-year period from human fecal samples in Lothian, Scotland, were sequenced with the Ion Torrent Personal Genome Machine. A total of 8,721 variable sites in the core genome were identified among the 105 isolates; 47% of the single nucleotide polymorphisms (SNPs) were attributable to six “atypical” E. coli O157 strains and included recombinant regions. Phylogenetic analyses showed that WGS correlated well with the epidemiological data. Epidemiological links existed between cases whose isolates differed by three or fewer SNPs. WGS also correlated well with multilocus variable-number tandem repeat analysis (MLVA) typing data, with only three discordant results observed, all among isolates from cases not known to be epidemiologically related. WGS produced a better-supported, higher-resolution phylogeny than MLVA, confirming that the method is more suitable for epidemiological surveillance of E. coli O157. A combination of insilico analyses (VirulenceFinder, ResFinder, and local BLAST searches) were used to determine stx subtypes, multilocus sequence types (15 loci), and the presence of virulence and acquired antimicrobial resistance genes. There was a high level of correlation between the WGS data and our routine typing methods, although some discordant results were observed, mostly related to the limitation of short sequence read assembly. The data were used to identify sublineages and clades of E. coli O157, and when they were correlated with the clinical outcome data, they showed that one clade, Ic3, was significantly associated with severe disease. Together, the results show that WGS data can provide higher resolution of the relationships between E. coli O157 isolates than that provided by MLVA. The method has the potential to streamline the laboratory workflow and provide detailed information for the clinical management of patients and public health interventions.

Model-based clustering of Escherichia coli O157:H7 genotypes and their potential association with clinical outcome in human infections

This study addresses the potential association of Escherichia coli O157:H7 genetic clusters with severe clinical manifestations in humans. The genotypes used in this model-based clustering had been delineated on the basis of lineage-specific polymorphism assay, Shiga toxin-encoding bacteriophage insertion site assay, clade typing, tir (A255T) polymorphism, variant analysis of Shiga toxin 2 gene, and antiterminator Q genes. Based on this model, the distribution of genotypes among tested strains suggested the presence of 6 main genetic clusters of E. coli O157:H7 strains. Clusters 1 and 3 were observed to be more frequent among E. coli O157:H7 strains isolated from bloody diarrhea and hemolytic uremic syndrome, respectively. Consequently, our findings supported the growing evidence of the existence of distinct genotypes of E. coli O157:H7 that differ in their virulence levels to human.

Geographically distinct Escherichia coli O157 isolates differ by lineage, Shiga toxin genotype, and total shiga toxin production

While the differential association of Escherichia coli O157 genotypes with animal and human hosts has recently been well documented, little is known about their distribution between countries and how this might affect regional disease rates. Here, we used a 48-plex single nucleotide polymorphism (SNP) assay to segregate 148 E. coli O157 isolates from Australia, Argentina, and the United States into 11 SNP lineages. We also investigated the relationship between SNP lineages, Shiga toxin (Stx) gene profiles, and total Stx production. E. coli O157 isolates clearly segregated into SNP lineages that were differentially associated with each country. Of the 11 SNP lineages, seven were detected among isolates from a single country, two were detected among isolates from all three countries, and another two were detected only among U.S. and Argentinean isolates. A number of Australian (30%) and Argentinean (14%) isolates were associated with novel, previously undescribed SNP lineages that were unique to each country. Isolates within SNP lineages that were strongly associated with the carriage of stx2a produced comparatively more Stx on average than did those lacking the stx2a subtype. Furthermore, the proportion of isolates in stx2a-associated SNP lineages was significantly higher in Argentina and the United States than Australia (P < 0.05). This study provides evidence for the geographic divergence of E. coli O157 and for a prominent role of stx2a in total Stx production. These results also highlight the need for more comprehensive studies of the global distribution of E. coli O157 lineages and the impacts of regionally predominant E. coli O157 lineages on the prevalence and severity of disease.

Multilocus genotypic characterization of Escherichia coli O157:H7 recovered from food sources

Escherichia coli O157:H7 strains (n = 33) recovered from different food sources in Egypt were characterized using molecular assays to identify strain genotypes associated with various levels of pathogenic potential. Genotypic characterization included: lineage-specific polymorphism assay (LSPA-6), Shiga-toxin-encoding bacteriophage insertion site assay (SBI), clade 8 typing, Tir (A255 T) polymorphism, and variant analysis of Shiga toxin 2 gene (Stx 2a and Stx 2c), and anti-terminator Q genes (Q 933 and Q 21). Genotypes LI/II (76%), SBI 1 (60·6%), clade 8 (69·7%), Tir (255 T) (72·7%) and Stx 2c (45·5%) were found to be significantly more frequent compared to other genetic markers in the strains analysed. Multivariable analysis revealed a significant association between LPSA-6 and clade types as well as Tir (A255 T). To the best of our knowledge, this is the first study to report the characterization of these genetic markers in E. coli O157:H7 strains in the Middle East and Africa.

Genetic diversity of Shiga toxin-producing Escherichia coli O157 : H7 recovered from human and food sources

The aim of this study was to identify an epidemiological association between Shiga toxin-producing Escherichia coli O157 : H7 strains associated with human infection and with food sources. Frequency distributions of different genetic markers of E. coli O157 : H7 strains recovered from human and food sources were compared using molecular assays to identify E. coli O157 : H7 genotypes associated with variation in pathogenic potential and host specificity. Genotypic characterization included: lineage-specific polymorphism assay (LSPA-6), clade typing, tir (A255T) polymorphism, Shiga toxin-encoding bacteriophage insertion site analysis and variant analysis of Shiga toxin 2 gene (stx2a and stx2c) and antiterminator Q genes (Q933 and Q21). The intermediate lineage (LI/II) dominated among both food and human strains. Compared to other clades, clades 7 and 8 were more frequent among food and human strains, respectively. The tir (255T) polymorphism occurred more frequently among human strains than food strains. Q21 and Q933 + Q21 were found at significantly higher frequencies among food and human strains, respectively. Moreover, stx2a and stx2a+c were detected at significantly higher frequencies among human strains compared to food strains. Bivariate analysis revealed significant concordance (P<0.05) between the LSPA-6 assay and the other typing methods. Multivariable regression analysis suggested that tir (255T) was the most distinctive genotype that can be used to detect bacterial clones with potential risk for human illness from food sources. This study supported previous reports of the existence of diversity in genetic markers among different isolation sources by including E. coli O157 : H7 strains from both food and human sources. This might enable tracking genotypes with potential risk for human illness from food sources.

Phylogenetic characterization of Escherichia coli O157 : H7 based on IS629 distribution and Shiga toxin genotype

Shiga toxin (stx)-producing Escherichia coli O157 : H7 is a prominent food-borne pathogen. Symptoms in human infections range from asymptomatic to haemorrhagic colitis and haemolytic uraemic syndrome, and there is a need for methods that yield information that can be used to better predict clinical and epidemiological outcomes. IS629 is an insertion sequence notable for its prevalence and variable distribution in the chromosome of E. coli O157 : H7, which has been exploited for subtyping and strain characterization. In particular, IS629 distribution is closely aligned with the major phylogenetic lineages that are known to be distinctive in their genome structure and virulence potential. In the present study, a comprehensive subtyping method in which IS629-typing was combined with stx genotyping was developed using a conventional PCR approach. This method consisted of a set of 32 markers based on the unique distribution of IS629 in the three major phylogenetic lineages of E. coli O157 : H7 and six additional markers to determine the stx genotype, a key virulence signature associated with each lineage. The analysis of IS629 loci variation with the 32 markers allowed us to determine the IS629 distribution profile (IDP), phylogenetic lineage and genetic relatedness of the 31 E. coli O157 : H7 strains examined. An association between IDP typing and stx genotype was observed. The use of both IDP and the stx genotype for strain characterization provided confirmative and complementary data in support of lineage placement of closely related strains. In addition, IS629/stx profiles were in agreement with strain segregation based on LSPA-6 (lineage-specific polymorphism assay) and PFGE subtyping, demonstrating its potential as a subtyping and strain tracking method.

Characteristics of clinical Shiga toxin-producing Escherichia coli isolated from British Columbia

Shiga toxin-producing Escherichia coli (STEC) are significant public health threats. Although STEC O157 are recognized foodborne pathogens, non-O157 STEC are also important causes of human disease. We characterized 10 O157:H7 and 15 non-O157 clinical STEC derived from British Columbia (BC). Eae, hlyA, and stx were more frequently observed in STEC O157, and 80 and 100% of isolates possessed stx₁ and stx₂, respectively. In contrast, stx₁ and stx₂ occurred in 80 and 40% of non-O157 STEC, respectively. Comparative genomic fingerprinting (CGF) revealed three distinct clusters (C). STEC O157 was identified as lineage I (LI; LSPA-6 111111) and clustered as a single group (C1). The cdi gene previously observed only in LII was seen in two LI O157 isolates. CGF C2 strains consisted of diverse non-O157 STEC while C3 included only O103:H25, O118, and O165 serogroup isolates. With the exception of O121 and O165 isolates which were similar in virulence gene complement to STEC O157, C1 O157 STEC produced more Stx2 than non-O157 STEC. Antimicrobial resistance (AMR) screening revealed resistance or reduced sensitivity in all strains, with higher levels occurring in non-O157 STEC. One STEC O157 isolate possessed a mobile bla_CMY-2 gene transferrable across genre via conjugation.

Multilocus genotype analysis of Escherichia coli O157 isolates from Australia and the United States provides evidence of geographic divergence

Escherichia coli O157 is a food-borne pathogen whose major reservoir has been identified as cattle. Recent genetic information has indicated that populations of E. coli O157 from cattle and humans can differ genetically and that this variation may have an impact on their ability to cause severe human disease. In addition, there is emerging evidence that E. coli O157 strains from different geographical regions may also be genetically divergent. To investigate the extent of this variation, we used Shiga toxin bacteriophage insertion sites (SBI), lineage-specific polymorphisms (LSPA-6), multilocus variable-number tandem-repeat analysis (MLVA), and a tir 255T>A polymorphism to examine 606 isolates representing both Australian and U.S. cattle and human populations. Both uni- and multivariate analyses of these data show a strong association between the country of origin and multilocus genotypes (P < 0.0001). In addition, our results identify factors that may play a role in virulence that also differed in isolates from each country, including the carriage of stx1 in the argW locus uniquely observed in Australian isolates and the much higher frequency of stx2-positive (also referred to as stx2a) strains in the U.S. isolates (4% of Australian isolates versus 72% of U.S. isolates). LSPA-6 lineages differed between the two continents, with the majority of Australian isolates belonging to lineage I/II (LI/II) (LI, 2%; LI/II, 85%; LII, 13%) and the majority of U.S. isolates belonging to LI (LI, 60%; LI/II, 16%; LII, 25%). The results of this study provide strong evidence of phylogeographic structuring of E. coli O157 populations, suggesting divergent evolution of enterohemorrhagic E. coli O157 in Australia and the United States.

Characterization of Escherichia coli O157:H7 strains isolated from supershedding cattle

Previous reports have indicated that a small proportion of cattle shedding high levels of Escherichia coli O157:H7 is the main source for transmission of this organism between animals. Cattle achieving a fecal shedding status of 10(4) CFU of E. coli O157:H7/gram or greater are now referred to as supershedders. The aim of this study was to investigate the contribution of E. coli O157:H7 strain type to supershedding and to determine if supershedding was restricted to a specific set of E. coli O157:H7 strains. Fecal swabs (n = 5,086) were collected from cattle at feedlots or during harvest. Supershedders constituted 2.0% of the bovine population tested. Supershedder isolates were characterized by pulsed-field gel electrophoresis (PFGE), phage typing, lineage-specific polymorphism assay (LSPA), Stx-associated bacteriophage insertion (SBI) site determination, and variant analysis of Shiga toxin, tir, and antiterminator Q genes. Isolates representing 52 unique PFGE patterns, 19 phage types, and 12 SBI clusters were obtained from supershedding cattle, indicating that there is no clustering to E. coli O157:H7 genotypes responsible for supershedding. While being isolated directly from cattle, this strain set tended to have higher frequencies of traits associated with human clinical isolates than previously collected bovine isolates with respect to lineage and tir allele, but not for SBI cluster and Q type. We conclude that no exclusive genotype was identified that was common to all supershedder isolates.

Molecular and antimicrobial susceptibility analyses distinguish clinical from bovine Escherichia coli O157 strains

A population-based study combining (i) antimicrobial, (ii) genetic, and (iii) virulence analyses with molecular evolutionary analyses revealed segregative characteristics distinguishing human clinical and bovine Escherichia coli O157 strains from western Canada. Human (n = 50) and bovine (n = 50) strains of E. coli O157 were collected from Saskatchewan and Manitoba in 2006 and were analyzed by using the six-marker lineage-specific polymorphism assay (LSPA6), antimicrobial susceptibility analysis, the colicin assay, plasmid and virulence profiling including the eae, ehxA, espA, iha, stx1, stx2, stx2c, stx2d, stx2d-activatable, stx2e, and stx2f virulence-associated genes, and structure analyses. Multivariate logistic regression and Fisher's exact test strongly suggested that antimicrobial susceptibility was the most distinctive characteristic (P = 0.00487) associated with human strains. Among all genetic, virulence, and antimicrobial determinants, resistance to tetracycline (P < 0.000) and to sulfisoxazole (P < 0.009) were the most strongly associated segregative characteristics of bovine E. coli O157 strains. Among 11 virulence-associated genes, stx2c showed the strongest association with E. coli O157 strains of bovine origin. LSPA6 genotyping showed the dominance of the lineage I genotype among clinical (90%) and bovine (70%) strains, indicating the importance of lineage I in O157 epidemiology and ecology. Population structure analysis revealed that the more-diverse bovine strains came from a unique group of strains characterized by a high degree of antimicrobial resistance and high frequencies of lineage II genotypes and stx2c variants. These findings imply that antimicrobial resistance generated among bovine strains of E. coli O157 has a large impact on the population of this human pathogen.

Evolution of the Stx2-encoding prophage in persistent bovine Escherichia coli O157:H7 strains

Escherichia coli O157:H7 is a human pathogen that resides asymptomatically in its bovine host. The level of Shiga toxin (Stx) produced is variable in bovine-derived strains in contrast to human isolates that mostly produce high levels of Stx. To understand the genetic basis for varied Stx production, chronological collections of bovine isolates from Wisconsin dairy farms, R and X, were analyzed for multilocus prophage polymorphisms, stx(2) subtypes, and the levels of stx(2) transcript and toxin. The E. coli O157:H7 that persisted on both farms were phylogenetically distinct and yet produced little to no Stx2 due to gene deletions in Stx2c-encoding prophage (farm R) or insertional inactivation of stx(2a) by IS1203v (farm X). Loss of key regulatory and lysis genes in Stx2c-encoding prophage abolished stx(2c) transcription and induction of the prophage and stx(2a)::IS1203v in Stx2a-encoding prophage generated a truncated stx(2a) mRNA without affecting phage production. Stx2-producing strains were transiently present (farm R) and became Stx2 negative on farm X (i.e., stx(2a)::IS1203v). To our knowledge, this is the first study that details the evolution of E. coli O157:H7 and its Stx2-encoding prophage in a chronological collection of natural isolates. The data suggest the bovine and farm environments can be niches where Stx2-negative E. coli O157:H7 emerge and persist, which explains the Stx variability in bovine isolates and may be part of an evolutionary step toward becoming bovine specialists.

Carriage of stx2a differentiates clinical and bovine-biased strains of Escherichia coli O157

Background

Shiga toxin (Stx) are cardinal virulence factors of enterohemorrhagic E. coli O157:H7 (EHEC O157). The gene content and genomic insertion sites of Stx-associated bacteriophages differentiate clinical genotypes of EHEC O157 (CG, typical of clinical isolates) from bovine-biased genotypes (BBG, rarely identified among clinical isolates). This project was designed to identify bacteriophage-mediated differences that may affect the virulence of CG and BBG.

Methods

Stx-associated bacteriophage differences were identified by whole genome optical scans and characterized among >400 EHEC O157 clinical and cattle isolates by PCR.

Results

Optical restriction maps of BBG strains consistently differed from those of CG strains only in the chromosomal insertion sites of Stx2-associated bacteriophages. Multiplex PCRs (stx1, stx2a, and stx2c as well as Stx-associated bacteriophage - chromosomal insertion site junctions) revealed four CG and three BBG that accounted for >90% of isolates. All BBG contained stx2c and Stx2c-associated bacteriophage – sbcB junctions. All CG contained stx2a and Stx2a-associated bacteriophage junctions in wrbA or argW.

Conclusions

Presence or absence of stx2a (or another product encoded by the Stx2a-associated bacteriophage) is a parsimonious explanation for differential virulence of BBG and CG, as reflected in the distributions of these genotypes in humans and in the cattle reservoir.

A comparison of Shiga-toxin 2 bacteriophage from classical enterohemorrhagic Escherichia coli serotypes and the German E. coli O104:H4 outbreak strain

Escherichia coli O104:H4 was associated with a severe foodborne disease outbreak originating in Germany in May 2011. More than 4000 illnesses and 50 deaths were reported. The outbreak strain was a typical enteroaggregative E. coli (EAEC) that acquired an antibiotic resistance plasmid and a Shiga-toxin 2 (Stx2)-encoding bacteriophage. Based on whole-genome phylogenies, the O104:H4 strain was most closely related to other EAEC strains; however, Stx2-bacteriophage are mobile, and do not necessarily share an evolutionary history with their bacterial host. In this study, we analyzed Stx2-bacteriophage from the E. coli O104:H4 outbreak isolates and compared them to all available Stx2-bacteriophage sequences. We also compared Stx2 production by an E. coli O104:H4 outbreak-associated isolate (ON-2011) to that of E. coli O157:H7 strains EDL933 and Sakai. Among the E. coli Stx2-phage sequences studied, that from O111:H- strain JB1-95 was most closely related phylogenetically to the Stx2-phage from the O104:H4 outbreak isolates. The phylogeny of most other Stx2-phage was largely concordant with their bacterial host genomes. Finally, O104:H4 strain ON-2011 produced less Stx2 than E. coli O157:H7 strains EDL933 and Sakai in culture; however, when mitomycin C was added, ON-2011 produced significantly more toxin than the E. coli O157:H7 strains. The Stx2-phage from the E. coli O104:H4 outbreak strain and the Stx2-phage from O111:H- strain JB1-95 likely share a common ancestor. Incongruence between the phylogenies of the Stx2-phage and their host genomes suggest the recent Stx2-phage acquisition by E. coli O104:H4. The increase in Stx2-production by ON-2011 following mitomycin C treatment may or may not be related to the high rates of hemolytic uremic syndrome associated with the German outbreak strain. Further studies are required to determine whether the elevated Stx2-production levels are due to bacteriophage or E. coli O104:H4 host related factors.

Phylogenetically related Argentinean and Australian Escherichia coli O157 isolates are distinguished by virulence clades and alternative Shiga toxin 1 and 2 prophages

Shiga toxigenic Escherichia coli O157 is the leading cause of hemolytic uremic syndrome (HUS) worldwide. The frequencies of stx genotypes and the incidences of O157-related illness and HUS vary significantly between Argentina and Australia. Locus-specific polymorphism analysis revealed that lineage I/II (LI/II) E. coli O157 isolates were most prevalent in Argentina (90%) and Australia (88%). Argentinean LI/II isolates were shown to belong to clades 4 (28%) and 8 (72%), while Australian LI/II isolates were identified as clades 6 (15%), 7 (83%), and 8 (2%). Clade 8 was significantly associated with Shiga toxin bacteriophage insertion (SBI) type stx(2) (locus of insertion, argW) in Argentinean isolates (P < 0.0001). In Argentinean LI/II strains, stx(2) is carried by a prophage inserted at argW, whereas in Australian LI/II strains the argW locus is occupied by the novel stx(1) prophage. In both Argentinean and Australian LI/II strains, stx(2c) is almost exclusively carried by a prophage inserted at sbcB. However, alternative q(933)- or q(21)-related alleles were identified in the Australian stx(2c) prophage. Argentinean LI/II isolates were also distinguished from Australian isolates by the presence of the putative virulence determinant ECSP_3286 and the predominance of motile O157:H7 strains. Characteristics common to both Argentinean and Australian LI/II O157 strains included the presence of putative virulence determinants (ECSP_3620, ECSP_0242, ECSP_2687, ECSP_2870, and ECSP_2872) and the predominance of the tir255T allele. These data support further understanding of O157 phylogeny and may foster greater insight into the differential virulence of O157 lineages.

Enterohemorrhagic Escherichia coli O157:H7 Shiga toxins inhibit gamma interferon-mediated cellular activation

Enterohemorrhagic Escherichia coli (EHEC) serotype O157:H7 is a food-borne pathogen that causes significant morbidity and mortality in developing and industrialized nations. EHEC infection of host epithelial cells is capable of inhibiting the gamma interferon (IFN-γ) proinflammatory pathway through the inhibition of Stat-1 phosphorylation, which is important for host defense against microbial pathogens. The aim of this study was to determine the bacterial factors involved in the inhibition of Stat-1 tyrosine phosphorylation. Human HEp-2 and Caco-2 epithelial cells were challenged directly with either EHEC or bacterial culture supernatants and stimulated with IFN-γ, and then the protein extracts were analyzed by immunoblotting. The data showed that IFN-γ-mediated Stat-1 tyrosine phosphorylation was inhibited by EHEC secreted proteins. Using two-dimensional difference gel electrophoresis, EHEC Shiga toxins were identified as candidate inhibitory factors. EHEC Shiga toxin mutants were then generated and complemented in trans, and mutant culture supernatant was supplemented with purified Stx to confirm their ability to subvert IFN-γ-mediated cell activation. We conclude that while other factors are likely involved in the suppression of IFN-γ-mediated Stat-1 tyrosine phosphorylation, E. coli-derived Shiga toxins represent a novel mechanism by which EHEC evades the host immune system.

Variation in stress resistance patterns among stx genotypes and genetic lineages of shiga toxin-producing Escherichia coli O157

To evaluate the relationship between bacterial genotypes and stress resistance patterns, we exposed 57 strains of Shiga toxin-producing Escherichia coli (STEC) O157 to acid, freeze-thaw, heat, osmotic, oxidative, and starvation stresses. Inactivation rates were calculated in each assay and subjected to univariate and multivariate analyses, including principal component analysis (PCA) and cluster analysis. The stx genotype was determined for each strain as was the lineage-specific polymorphism assay (LSPA6) genotype. In univariate analyses, strains of the stx(1) stx(2) genotype showed greater resistance to heat than strains of the stx(1) stx(2c) genotype; moreover, strains of the stx(1) stx(2) genotype showed greater resistance to starvation than strains of the stx(2) or stx(2c) genotypes. LSPA6 lineage I (LI) strains showed greater resistance to heat and starvation than LSPA6 lineage II (LII) strains. PCA revealed a general trend that a strain with greater resistance to one type of stress tended to have greater resistance to other types of stresses. In cluster analysis, STEC O157 strains were grouped into stress-resistant, stress-sensitive, and intermediate clusters. In stx genotypes, all strains of the stx(1) stx(2) genotype were grouped with the stress-resistant cluster, whereas 72.7% (8/11) of strains of the stx(1) stx(2c) genotype grouped with the stress-sensitive cluster. In LI strains, 77.8% (14/18) of the strains were grouped with the stress-resistant cluster, whereas 64.7% (11/17) of LII strains were grouped with the stress-sensitive cluster. These results indicate that the genotypes of STEC O157 that are frequently associated with human illness, i.e., LI or the stx(1) stx(2) genotype, have greater multiple stress resistance than do strains of other genotypes.

Genetic features differentiating bovine, food, and human isolates of shiga toxin-producing Escherichia coli O157 in The Netherlands

The frequency of Escherichia coli O157 genotypes among bovine, food, and human clinical isolates from The Netherlands was studied. Genotyping included the lineage-specific polymorphism assay (LSPA6), the Shiga-toxin-encoding bacteriophage insertion site assay (SBI), and PCR detection and/or subtyping of virulence factors and markers [stx1, stx(2a)/stx(2c), q21/Q933, tir(A255T), and rhsA(C3468G)]. LSPA6 lineage II dominated among bovine isolates (63%), followed by lineage I/II (35.6%) and lineage I (1.4%). In contrast, the majority of the human isolates were typed as lineage I/II (77.6%), followed by lineage I (14.1%) and lineage II (8.2%). Multivariate analysis revealed that the tir(A255T) SNP and the stx(2a)/stx(2c) gene variants were the genetic features most differentiating human from bovine isolates. Bovine and food isolates were dominated by stx(2c) (86.4% and 65.5%, respectively). Among human isolates, the frequency of stx(2c) was 36.5%, while the frequencies of stx(2a) and stx(2a) plus stx(2c) were 41.2% and 22.4%, respectively. Bovine isolates showed equal distribution of tir(255A) (54.8%) and tir(255T) (45.2%), while human isolates were dominated by the tir(255T) genotype (92.9%). LSPA6 lineage I isolates were all genotype stx(2c) and tir(255T), while LSPA6 lineage II was dominated by tir(255A) (86.4%) and stx(2c) (90.9%). LSPA6 lineage I/II isolates were all genotype tir(255T) but showed more variation in stx(2) types. The results support the hypothesis that in The Netherlands, the genotypes primarily associated with human disease form a minor subpopulation in the bovine reservoir. Comparison with published data revealed that the distribution of LSPA6 lineages among bovine and human clinical isolates differs considerably between The Netherlands and North America.

Genome signatures of Escherichia coli O157:H7 isolates from the bovine host reservoir

Cattle comprise a main reservoir of Shiga toxin-producing Escherichia coli O157:H7 (STEC). The significant differences in host prevalence, transmissibility, and virulence phenotypes among strains from bovine and human sources are of major interest to the public health community and livestock industry. Genomic analysis revealed divergence into three lineages: lineage I and lineage I/II strains are commonly associated with human disease, while lineage II strains are overrepresented in the asymptomatic bovine host reservoir. Growing evidence suggests that genotypic differences between these lineages, such as polymorphisms in Shiga toxin subtypes and synergistically acting virulence factors, are correlated with phenotypic differences in virulence, host ecology, and epidemiology. To assess the genomic plasticity on a genome-wide scale, we have sequenced the whole genome of strain EC869, a bovine-associated E. coli O157:H7 isolate. Comparative phylogenomic analysis of this key isolate enabled us to place accurately bovine lineage II strains within the genetically homogenous E. coli O157:H7 clade. Identification of polymorphic loci that are anchored both in the chromosomal backbone and horizontally acquired regions allowed us to associate bovine genotypes with altered virulence phenotypes and host prevalence. This study catalogued numerous novel lineage II-specific genome signatures, some of which appear to be associated intimately with the altered pathogenic potential and niche adaptation within the bovine rumen. The presented extended list of polymorphic markers is valuable in the development of a robust typing system critical for forensic, diagnostic, and epidemiological studies of this emerging human pathogen.

Characterization of the relationship between integrase, excisionase and antirepressor activities associated with a superinfecting Shiga toxin encoding bacteriophage

Shigatoxigenic Escherichia coli emerged as new food borne pathogens in the early 1980s, primarily driven by the dispersal of Shiga toxin-encoding lambdoid bacteriophages. At least some of these Stx phages display superinfection phenotypes, which differ significantly from lambda phage itself, driving through in situ recombination further phage evolution, increasing host range and potentially increasing the host's pathogenic profile. Here, increasing levels of Stx phage Φ24(B) integrase expression in multiple lysogen cultures are demonstrated along with apparently negligible repression of integrase expression by the cognate CI repressor. The Φ24(B) int transcription start site and promoter region were identified and found to differ from in silico predictions. The unidirectional activity of this integrase was determined in an in situ, inducible tri-partite reaction. This indicated that Φ24(B) must encode a novel directionality factor that is controlling excision events during prophage induction. This excisionase was subsequently identified and characterized through complementation experiments. In addition, the previous proposal that a putative antirepressor was responsible for the lack of immunity to superinfection through inactivation of CI has been revisited and a new hypothesis involving the role of this protein in promoting efficient induction of the Φ24(B) prophage is proposed.

Multivariate analyses revealed distinctive features differentiating human and cattle isolates of Shiga toxin-producing Escherichia coli O157 in Japan

Genotypes of Shiga toxin-producing Escherichia coli (STEC) O157 isolated from humans and cattle were analyzed by uni- and multivariable logistic regression, and population structure methods, to gain insight into transmission and the nature of human infection. Eleven genotyping assays, including PCR typing of five virulence factors (stx(1), stx(2), stx(2c), eae, and ehxA) and a lineage-specific polymorphism assay using six markers (LSPA6), were considered in the analyses. The prevalence of the stx(1), stx(2), and stx(2c) virulence factors was significantly different between human and cattle isolates. However, multivariable regression revealed that the presence of only the stx(2) gene was significantly associated with human isolates after controlling for confounding effects. LSPA6 typing demonstrated an apparent difference in the distribution of LSPA6 lineages between human and cattle isolates and a strong association between stx genotypes and LSPA6 genotypes. Population genetics tools identified three genetically distinct clusters of STEC O157. Each cluster was characterized by stx genotypes and LSPA6 genotypes. The human isolates typically comprised LSPA6 lineage I with stx(1) stx(2) strains and LSPA6 lineage I/II with stx(2c) or stx(2) stx(2c) strains [corrected]. In contrast, the cattle isolates comprised LSPA6 lineage II strains withstx(2c) or stx(1) stx(2c) strains [corrected] in addition to the clusters identified for the human isolates. Our analyses provide new evidence that the stx(2) gene is the most distinctive feature in human isolates compared to cattle isolates in Japan, and only a subset of the genetically diverse population isolated from cattle is involved in human illnesses. Our results may contribute to international comparisons and risk assessments of STEC O157.

Gene cluster conferring streptomycin, sulfonamide, and tetracycline resistance in Escherichia coli O157:H7 phage types 23, 45, and 67

Multidrug resistance to streptomycin, sulfonamide, and tetracycline (AMR-SSuT) was identified in 156 of 171 isolates of Escherichia coli O157:H7 of phage types 23, 45, and 67. In 154 AMR-SSuT isolates, resistance was encoded by strA, strB, sul2, and tet(B), which in 59 of 63 tested isolates were found clustered together on the chromosome within the cdiA locus.

Lineage and host source are both correlated with levels of Shiga toxin 2 production by Escherichia coli O157:H7 strains

Escherichia coli O157:H7 strains fall into three major genetic lineages that differ in their distribution among humans and cattle. Several recent studies have reported differences in the expression of virulence factors between E. coli O157:H7 strains from these two host species. In this study, we wished to determine if important virulence-associated "mobile genetic elements" such as Shiga toxin 2 (Stx2)-encoding prophage are lineage restricted or are host source related and acquired independently of the pathogen genotype. DNA sequencing of the stx(2) flanking region from a lineage II (LII) strain, EC970520, revealed that the transcriptional activator gene Q in LI strain EDL933 (upstream of stx(2)) is replaced by a pphA (serine/threonine phosphatase) homologue and an altered Q gene in this and all other LII strains tested. In addition, nearly all LI strains carried stx(2), whereas all LII strains carried variant stx(2c) and 4 of 14 LI/II strains had copies of both stx(2) and variant stx(2c). Real-time PCR (RT-PCR) and enzyme-linked immunosorbent assay (ELISA) demonstrated that LI and LI/II strains produce significantly more stx(2) mRNA and Stx2 than LII strains. However, among LI strains significantly more Stx2 is also produced by strains from humans than from cattle. Therefore, lineage-associated differences among E. coli O157:H7 strains such as prophage content, toxin type, and toxin expression may contribute to host isolation bias. However, the level of Stx2 production alone may also play an important role in the within-lineage association of E. coli O157:H7 strains with human clinical disease.

rhs genes are potential markers for multilocus sequence typing of Escherichia coli O157:H7 strains

DNA sequence-based molecular subtyping methods such as multilocus sequence typing (MLST) are commonly used to generate phylogenetic inferences for monomorphic pathogens. The development of an effective MLST scheme for subtyping Escherichia coli O157:H7 has been hindered in the past due to the lack of sequence variation found within analyzed housekeeping and virulence genes. A recent study suggested that rhs genes are under strong positive selection pressure, and therefore in this study we analyzed these genes within a diverse collection of E. coli O157:H7 strains for sequence variability. Eighteen O157:H7 strains from lineages I and II and 15 O157:H7 strains from eight clades were included. Examination of these rhs genes revealed 44 polymorphic loci (PL) and 10 sequence types (STs) among the 18 lineage strains and 280 PL and 12 STs among the 15 clade strains. Phylogenetic analysis using rhs genes generally grouped strains according to their known lineage and clade classifications. These findings also suggested that O157:H7 strains from clades 6 and 8 fall into lineage I/II and that strains of clades 1, 2, 3, and 4 fall into lineage I. Additionally, unique markers were found in rhsA and rhsJ that might be used to define clade 8 and clade 6. Therefore, rhs genes may be useful markers for phylogenetic analysis of E. coli O157:H7.

In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence

Background

Many approaches have been used to study the evolution, population structure and genetic diversity of Escherichia coli O157:H7; however, observations made with different genotyping systems are not easily relatable to each other. Three genetic lineages of E. coli O157:H7 designated I, II and I/II have been identified using octamer-based genome scanning and microarray comparative genomic hybridization (mCGH). Each lineage contains significant phenotypic differences, with lineage I strains being the most commonly associated with human infections. Similarly, a clade of hyper-virulent O157:H7 strains implicated in the 2006 spinach and lettuce outbreaks has been defined using single-nucleotide polymorphism (SNP) typing. In this study an in silico comparison of six different genotyping approaches was performed on 19 E. coli genome sequences from 17 O157:H7 strains and single O145:NM and K12 MG1655 strains to provide an overall picture of diversity of the E. coli O157:H7 population, and to compare genotyping methods for O157:H7 strains.

Results

In silico determination of lineage, Shiga-toxin bacteriophage integration site, comparative genomic fingerprint, mCGH profile, novel region distribution profile, SNP type and multi-locus variable number tandem repeat analysis type was performed and a supernetwork based on the combination of these methods was produced. This supernetwork showed three distinct clusters of strains that were O157:H7 lineage-specific, with the SNP-based hyper-virulent clade 8 synonymous with O157:H7 lineage I/II. Lineage I/II/clade 8 strains clustered closest on the supernetwork to E. coli K12 and E. coli O55:H7, O145:NM and sorbitol-fermenting O157 strains.

Conclusion

The results of this study highlight the similarities in relationships derived from multi-locus genome sampling methods and suggest a "common genotyping language" may be devised for population genetics and epidemiological studies. Future genotyping methods should provide data that can be stored centrally and accessed locally in an easily transferable, informative and extensible format based on comparative genomic analyses.

Escherichia coli O157:H7 strain origin, lineage, and Shiga toxin 2 expression affect colonization of cattle

Enterohemorrhagic Escherichia coli O157:H7 has evolved into an important human pathogen with cattle as the main reservoir. The recent discovery of E. coli O157:H7-induced pathologies in challenged cattle has suggested that previously discounted bacterial virulence factors may contribute to the colonization of cattle. The objective of the present study was to examine the impact of lineage type, cytotoxin activity, and cytotoxin expression on the amount of E. coli O157:H7 colonization of cattle tissue and cells in vitro. Using selected bovine- and human-origin strains, we determined that lineage type predicted the amount of E. coli O157:H7 strain colonization: lineage I > intermediate lineages > lineage II. All E. coli O157:H7 strain colonization was dose dependent, with threshold colonization at 10(3) to 10(5) CFU and maximum colonization at 10(7) CFU. We also determined that an as-yet-unknown factor of strain origin was the most dominant predictor of the amount of strain colonization in vitro. The amount of E. coli O157:H7 colonization was also influenced by strain cytotoxin activity and the inclusion of cytotoxins from lineage I or intermediate lineage strains increased colonization of a lineage II strain. There was a higher level of expression of the Shiga toxin 1 gene (stx(1)) in human-origin strains than in bovine-origin strains. In addition, lineage I strains expressed higher levels of the Shiga toxin 2 gene (stx(2)). The present study supports a role for strain origin, lineage type, cytotoxin activity, and stx(2) expression in modulating the amount of E. coli O157:H7 colonization of cattle.

Genomic regions conserved in lineage II Escherichia coli O157:H7 strains

Populations of the food- and waterborne pathogen Escherichia coli O157:H7 are comprised of two major lineages. Recent studies have shown that specific genotypes within these lineages differ substantially in the frequencies with which they are associated with human clinical disease. While the nucleotide sequences of the genomes of lineage I strains E. coli O157 Sakai and EDL9333 have been determined, much less is known about the genomes of lineage II strains. In this study, suppression subtractive hybridization (SSH) was used to identify genomic features that define lineage II populations. Three SSH experiments were performed, yielding 1,085 genomic fragments consisting of 811 contigs. Bacteriophage sequences were identified in 11.3% of the contigs, 9% showed insertions and 2.3% deletions with respect to E. coli O157:H7 Sakai, and 23.2% did not have significant identity to annotated sequences in GenBank. In order to test for the presence of these novel loci in lineage I and II strains, 27 PCR primer sets were designed based on sequences from these contigs. All but two of these PCR targets were found in the majority (51.9% to 100%) of 27 lineage II strains but in no more than one (<6%) of the 17 lineage I strains. Several of these lineage II-related fragments contain insertions/deletions that may play an important role in virulence. These lineage II-related loci were also shown to be useful markers for genotyping of E. coli O157:H7 strains isolated from human and animal sources.

Super-shedding and the link between human infection and livestock carriage of Escherichia coli O157

Cattle that excrete more Escherichia coli O157 than others are known as super-shedders. Super-shedding has important consequences for the epidemiology of E. coli O157 in cattle--its main reservoir--and for the risk of human infection, particularly owing to environmental exposure. Ultimately, control measures targeted at super-shedders may prove to be highly effective. We currently have only a limited understanding of both the nature and the determinants of super-shedding. However, super-shedding has been observed to be associated with colonization at the terminal rectum and might also occur more often with certain pathogen phage types. More generally, epidemiological evidence suggests that super-shedding might be important in other bacterial and viral infections.

Bacteriophage 2851 is a prototype phage for dissemination of the Shiga toxin variant gene 2c in Escherichia coli O157:H7

The production of Shiga toxin (Stx) (verocytotoxin) is a major virulence factor of Escherichia coli O157:H7 strains (Shiga toxin-producing E. coli [STEC] O157). Two types of Shiga toxins, designated Stx1 and Stx2, are produced in STEC O157. Variants of the Stx2 type (Stx2, Stx2c) are associated with high virulences of these strains for humans. A bacteriophage designated 2851 from a human STEC O157 encoding the Stx2c variant was described previously. Nucleotide sequence analysis of the phage 2851 genome revealed 75 predicted coding sequences and indicated a mosaic structure typical for lambdoid phages. Analyses of free phages and K-12 phage 2851 lysogens revealed that upon excision from the bacterial chromosome, the loss of a phage-encoded IS629 element leads to fusion of phage antA and antB genes, with the generation of a recombined antAB gene encoding a strong antirepressor. In wild-type E. coli O157 as well as in K-12 strains, phage 2851 was found to be integrated in the sbcB locus. Additionally, phage 2851 carries an open reading frame which encodes an OspB-like type III effector similar to that found in Shigella spp. Investigation of 39 stx(2c) E. coli O157 strains revealed that all except 1 were positive for most phage 2851-specific genes and possessed a prophage with the same border sequences integrated into the sbcB locus. Phage 2851-specific sequences were absent from most stx(2c)-negative E. coli O157 strains, and we suggest that phage 2851-like phages contributed significantly to the dissemination of the Stx2c variant toxin within this group of E. coli.

Rapid determination of Escherichia coli O157:H7 lineage types and molecular subtypes by using comparative genomic fingerprinting

In this study, variably absent or present (VAP) regions discovered through comparative genomics experiments were targeted for the development of a rapid, PCR-based method to subtype and fingerprint Escherichia coli O157:H7. Forty-four VAP loci were analyzed for discriminatory power among 79 E. coli O157:H7 strains of 13 phage types (PT). Twenty-three loci were found to maximize resolution among strains, generating 54 separate fingerprints, each of which contained strains of unique PT. Strains from the three previously identified major E. coli O157:H7 lineages, LSPA6-LI, LSPA6-LI/II, and LSPA6-LII, formed distinct branches on a dendrogram obtained by hierarchical clustering of comparative genomic fingerprinting (CGF) data. By contrast, pulsed-field gel electrophoresis (PFGE) typing generated 52 XbaI digestion profiles that were not unique to PT and did not cluster according to O157:H7 lineage. Our analysis identified a subpopulation comprised of 25 strains from a closed herd of cattle, all of which were of PT87 and formed a cluster distinct from all other E. coli O157:H7 strains examined. CGF found five related but unique fingerprints among the highly clonal herd strains, with two dominant subtypes characterized by a shift from the presence of locus fprn33 to its absence. CGF had equal resolution to PFGE typing but with greater specificity, generating fingerprints that were unique among phenotypically related E. coli O157:H7 lineages and PT. As a comparative genomics typing method that is amenable for use in high-throughput platforms, CGF may be a valuable tool in outbreak investigations and strain characterization.