Skip to main content
Download PDF

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

BMC Genomics volume 10, Article number: 287 (2009) Cite this article

Abstract

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.

Background

Escherichia coli O157:H7 is the most commonly implicated serotype of Shiga-toxin producing E. coli (STEC) or enterohemorrhagic E. coli (EHEC) associated with hemorrhagic colitis and the hemolytic uremic syndrome (HUS) [1, 2]. It is recognized world-wide as an important cause of both sporadic cases and outbreaks of food- and waterborne disease. Genomic diversity within populations of this pathogen is extensive and genome comparisons have revealed many DNA insertion/deletion and recombination events which are thought to be driven primarily through bacteriophage-mediated and other mechanisms of lateral gene transfer [3–5].

Three genetic lineages of E. coli O157:H7 have been described using octamer-based genome scanning (OBGS) and microarray-based comparative genomic hybridization (mCGH). Lineage I strains are commonly isolated from both cattle and clinically ill humans, lineage II strains are isolated primarily from cattle and intermediate lineage I/II strains are less well characterized with respect to phenotype and host distribution [6–8]. Additional studies on the evolution, population structure and genetic diversity of Escherichia coli O157:H7 have been carried out using a number of different genotyping approaches, each of which is based on targeting polymorphisms in a particular locus or set of loci. For example, the presence of Shiga-toxin (Stx) containing bacteriophage integration sites has been used to describe sixteen E. coli O157:H7 genotypes (with the majority associated with strains of bovine origin) [9]. Multi-locus variable number tandem repeat analysis (MLVA), which targets the number of repeats at nine genetic loci in the method proposed for PulseNet, has recently been used for epidemiological typing of O157:H7 strains [10]. mCGH, which examines the presence or absence of every gene in the genome of a reference strain or strains, has been used to elucidate the step-wise emergence of O157:H7 strains from an O55:H7-like ancestor [11] and by our laboratory to first characterize lineage I/II strains [7]. Single-nucleotide polymorphism (SNP) typing examines the single nucleotide changes throughout the E. coli O157:H7 genome; SNPs in 96 loci have been used to delineate 39 genotypes in over 500 strains, which have been partitioned into nine evolutionary clades [12]. One of these clades (clade 8) was found to contain putative hyper-virulent strains, including those implicated in the 2006 spinach- and lettuce-associated outbreaks in the United States. Comparative genomic fingerprinting (CGF), which examines the presence or absence of the most common variable loci within the genome, has recently been used with 23 loci to analyze 79 O157:H7 strains, and to group them into lineages and epidemiologically and phenotypically related clusters [13].

The various methods have proved useful on their own but in most cases the strain groupings from one method are not easily relatable to another, and none of these methods has yet achieved the "gold standard" status to be used as a common genotyping method. Because each of these genotyping approaches can provide important information on isolates with different pathogenic and virulence characteristics, the need to replicate each of these methods using a common group of reference strains has arisen. However, present restrictions on international strain exchange make the acquisition of certain reference strains nearly impossible for researchers outside the country of origin of the strains. While the stepwise emergence of E. coli O157:H7 from an O55:H7-like serotype has been advanced [11, 14], the relationship among E. coli O157:H7 lineages and clades incorporating data from standard and more advanced genotyping methods has not been evaluated. Thus, discoveries such as the recently described hyper-virulent clade of E. coli O157:H7 strains [12] are framed only in the context of the individual study.

Recently, three complete E. coli O157:H7 strain sequences and 11 whole genome shotgun sequences have become available in Genbank. With the additional sequencing of one O157:H7 and one O145:NM strain by our group, the availability of multiple whole-genome sequences has allowed us to apply several molecular typing schemes in silico to a group of E. coli O157:H7 and related strains. We have applied six molecular typing methods shown to differentiate among E. coli O157:H7 strains to a panel of nineteen genome-sequenced strains: Stx-bacteriophage insertion site typing [9], CGF [13], mCGH [7, 11], SNP genotyping [12], genomic in silico subtractive hybridization (GISSH) (Laing et al. in preparation) and MLVA[10].

In this study we provide a view of the relationships among E. coli O157:H7 strains based on a supernetwork representation of these six separate molecular typing methods. We also show how hyper-virulent clade 8 E. coli O157:H7 strains fit into this scheme and suggest the use of genotyping approaches that are relatable as part of a "common genotyping language" until whole genome sequencing becomes routine in bacterial strain genotyping.

Results

Lineage typing

The is- PCR LSPA6 typing of the E. coli O157:H7 strains in Table 1 designated strains EC4501 and TW14588 as lineage I, eleven of the 14 GenBank strains and the in-house strain EC71074 as lineage I/II, and strain EC869 as lineage II.

Table 1 The 19 E. coli strains analyzed in silico in this study.
Full size table

Stx-phage insertion site typing

Besser et al. [9] found greater diversity within O157:H7 strains isolated from cattle than those from human clinical illness, but the relationship between their typing scheme and lineage groupings was not explored. In Table 1 it is evident from in silico analysis that all genotypes based on the Stx-phage integration site typing can be related to mCGH lineage.

All lineage I strains were of genotype 3 in that they possessed both stx1 and stx2 and both bacteriophage integration sites yehV and wrbA were occupied. Additionally, they all possessed the "K" form of the N135K polymorphism in the FimH mannose binding pocket, lacked stx2c, had the "T" form of the T238A polymorphism in tir and were positive for the Q-stx2 junction of phage 933W.

All lineage I/II strains were of genotype 1 in that they contained stx2 but not stx1 and contained an occupied yehV and an intact wrbA bacteriophage integration site. Additionally, they all possessed the "N" form of the N135K polymorphism in the FimH mannose binding pocket and the "T" form of the T238A polymorphism in tir. However, the Q-stx2 junction of phage 933W and stx2c were variably present among the lineage I/II strains.

Lineage II strain EC869 was of genotype 6. It contained both stx1 and stx2 had the variant-R form of the yehV bacteriophage right junction and an intact wrbA integration site. EC869 possessed stx2c, the "A" form of the T238A polymorphism in tir, the "N" form of the N135K polymorphism in the FimH mannose binding pocket and was negative for the Q-stx2 junction of phage 933W. The raw data for this and all other in silico typing methods are provided [see Additional file 1].

SNP genotyping

The clade of each strain based on SNP typing was determined as described by Riordan et al. [15] where SNPs in the four genes ECs2357, ECs2521, ECs3881 and ECs4130 allow the delineation of the four most common clades of E. coli O157:H7. Our analysis typed EDL933 as clade 3 (SNP profile in the gene order given above as CGCT), Sakai as clade 1 (SNP profile CTTT), TW14588 and EC4501 as clade 2 (SNP profile CGTT), EC869 as one of clades 4–7 (SNP profile CGCC) and the remaining O157:H7 strains as members of the hyper-virulent clade 8 (SNP profile AGCC) (Table 1).

GISSH-based novel region distribution

The microarray used in the study by Zhang et al. [7] was based on lineage I strains EDL933 and Sakai and the K12 strain MG1655, so it is not surprising that the lineage I strains were found to contain more of the lineage I markers than strains from the other lineages. To account for novel genomic regions present in other O157:H7 strains, but not found in EDL933 or Sakai, an in silico subtractive hybridization was performed on every O157:H7 sequence in Table 1 against EDL933 and Sakai (Laing et al., in preparation). The study found 417 separate regions comprising 1456 segments of approximately 500 bp in length (~0.8 Mbp of novel DNA sequence). The distribution of these segments in silico is shown in Figure 1 and highlights the fact that lineage I strains possess genomic regions not represented in the original microarray probe set, as well as the fact that there are other lineage I/II and lineage II specific genomic regions.

Figure 1
figure 1

Novel regions distribution among the E. coli strains in Table 1. The distribution of 1456 regions ~500 bp in size among 17 E. coli O157:H7 strains and the O145:NM strain EC33264 and K12 strain MG1655. Regions are not necessarily contiguous and are defined as novel based on less than 80% sequence identity to the genome of either EDL933 or Sakai. Black indicates the presence of a region and white indicates the absence of a region.

Full size image

Combined analysis

The maximum parsimony trees from each method shared a number of commonalities [see Additional file 2]. All three lineages grouped distinctly with the following methods: CGF, SNP genotyping, mCGH and GISSH-based novel region typing. Stx-phage insertion site typing distinguished lineage I strains as a separate cluster but lineage II strain EC869 did not form a distinct cluster from the lineage I/II strains. Similarly, MLVA put lineage II strain EC869 on its own branch, whereas the lineage I/II and lineage I strains did not group separately with this method. All of the genotyping methods except mCGH identified lineage I/II strains EC508 and EC71074 as a separate cluster, which was placed into a group close to the lineage I strains in the CGF and Stx-phage insertion site based maximum parsimony trees and was placed between lineage I strain Sakai and the other three lineage I strains in the maximum parsimony tree based on MLVA data.

In order to obtain a more complete understanding of the relationships among E. coli O157:H7 strains, the maximum parsimony trees derived using in silico data from the six separate molecular typing methods included in this study were combined to form a supernetwork. We accomplished this by combing the six maximum parsimony trees generated using PAUP* instead of concatenating the data into a single matrix and generating a single most parsimonious tree. This ensured a method with many data points such as mCGH did not outweigh a method like MLVA that contains few data points. The supernetwork based on the six trees given equal weighting (Figure 2) shows competing signals, rather than an inferred absolute tree [16]. As can be seen, the 19 strains were distributed into four main groups. One contained strains K12 MG1655 and O145:NM EC33264 and another the single lineage II strain EC869. The remaining two groups of the supernetwork were lineage-specific, comprising one cluster of lineage I strains and another of lineage I/II strains. The distribution of the four lineage I strains was consistent with the clade breakdown of Manning et al. [12], as the clade 2 strains TW14588 and EC4501 were more closely related to each other than clade 3 strain EDL933 or clade 1 strain Sakai.

Figure 2
figure 2

Supernetwork constructed from in silico data of six typing methods. The supernetwork created from the combination of each maximum parsimony tree from the following typing methods: Stx-phage insertion site typing, MLVA, CGF, SNP genotyping, mCGH, and GISSH-based novel region distribution typing. Maximum parsimony trees were combined using the un-weighted mean distance and Z-closure with 1000 iterations; the resulting supernetwork was displayed using the equal angle method.

Full size image

Because only a relatively small number of strains have been sequenced, we repeated the analysis including experimental microarray [7, 11] and CGF data [13] for which lineage information was available to provide a larger context for our in silico analysis. This added 83 strains to the 19 strain in silico dataset and included data for three E. coli O55:H7 strains and two sorbitol-fermenting E. coli O157:H-strains (Table 2), allowing us to frame the resulting supernetwork within the context of the proposed stepwise emergence of E. coli O157:H7 [11]. The supernetwork in Figure 3 again showed four main groups. All lineage I and II strains formed discrete branches in the supernetwork, while the lineage I/II strains formed a cluster that was closest to E. coli K12 and strains of O55:H7, O145:NM and sorbitol-fermenting O157:H-.

Table 2 The 102 E. coli O157:H7 strains used in the construction of the supernetwork in Figure 3.
Full size table
Figure 3
figure 3

Supernetwork constructed from in silico data of six typing methods and experimental data from two methods. The supernetwork created from the combination of each maximum parsimony tree from the following typing methods: Stx-phage insertion site typing, MLVA, CGF, SNP genotyping, mCGH, and GISSH-based novel region distribution typing. Both mCGH and CGF datasets included the in silico data from the strains in Table 1 in addition to experimental data from the remaining strains in Table 2. Maximum parsimony trees were combined using the un-weighted mean distance and Z-closure with 1000 iterations; the resulting supernetwork was displayed using the equal angle method.

Full size image

Discussion

The acquisition and loss of genetic elements in E. coli O157:H7 is thought to affect the virulence of this pathogen in humans [17–19]. While attributes like the Stxs [20–22] and the locus of enterocyte attachment and effacement (LEE) [23] are well known, other less well characterized elements such as non-LEE effectors (NLE)s [24] may contribute to the spectrum of virulence that has been captured for STEC within the seropathotype classification [25].

This genomic diversity can largely be attributed to bacteriophages [4, 26] and other mobile elements [27] that cause DNA segment insertion and deletion events, rather than single-nucleotide changes [28]. The best characterization of diversity within E. coli O157:H7 must therefore take both single-nucleotide changes and large region turnover events into consideration. Whole genome sequencing is the best method for such analysis, but until the number of completely sequenced strains increases and whole genome sequencing becomes both routine and cost-effective, an estimation of whole genomic change based on sampling of polymorphisms or variability at a number of specific loci will have to suffice. A number of molecular genotyping methods have recently emerged that target different regions and types of variation, so the combination of data from these methods can offer a picture that is greater than the sum of their individual views. It can be argued that lateral gene transfer obscures the relationship among bacterial strains; however, once mobile elements are acquired and if they are stably maintained, they can be especially valuable in assessing strain relationships. We therefore advocate the use of genotyping methods that rely on multiple spatially distinct loci to provide a robust view of the O157:H7 population structure. Each individual multi-locus method in our study pointed to a similar tree structure and the combination of methods in the supernetwork (Figure 3) showed three distinct lineages, two of which were originally proposed by Kim et al. [6] and the third by Zhang et al. [7], adding confidence to the conclusion that these lineages exist.

The results of SNP genotyping over 500 O157:H7 isolates found hyper-virulent clade 8 strains, the causative agents in the spinach- and lettuce-related outbreaks of 2006 in the United States, to be most closely related to sorbitol fermenting E. coli O157:NM strains [12]. Based on our in silico analyses, these clade 8 strains, which appear to be increasing in prevalence and have a greater association with HUS than other strains [12], are members of lineage I/II (Table 1). Interestingly, E. coli O157:H7 strain TW14588 was designated as clade 8 in the initial publication on variation in virulence among the clades [12], but all in silico analyses conducted during this study on the whole-genome shotgun sequence available from GenBank (NZ_ABKY) suggests it belongs to clade 2 (lineage I).

It must be recognized that the publicly available genome sequences are not error-free and that this could have affected the architecture of trees that are highly dependent on single nucleotide changes in sequences of target genes, such as SNP-genotyping. However, the in silico tree architecture was very similar to that described by Manning et al. [12] based on experimental data; therefore we suspect that this was not a significant source of error in this study.

Differences between lineage I and lineage II strains have been described [29–31], but much less is known about lineage I/II strains with respect to host/disease association or expression of virulence attributes. We have recently demonstrated in a study of E. coli O157:H7 strains in Canada that certain phage types (PT)s are specific to O157:H7 lineages [29]. In the latter study, PT2 strains were shown to belong exclusively to lineage I/II, however, others strains in lineage I/II belonged to PTs that were not lineage-restricted e.g. PT23, PT8 and PT1, suggesting that this lineage is widespread and also diverse. Such diversity was also apparent in the examination of the novel regions of O157:H7 DNA, presented in Figure 1.

The Stxs, which are bacteriophage-encoded and the primary virulence factors of E. coli O157:H7 [20] show differential distribution among the lineages. The stx2 gene was found in nearly all O157:H7 lineage I and lineage I/II strains [12, 29], while the stx1 gene was absent in lineage I/II strains but present in nearly all other strains studied. Additionally, Ziebell et al. [29] found the stx2c gene in 96.7% of lineage II strains, 50.0% of lineage I/II strains and 1.8% of lineage I strains, while Manning et al. [12] found stx2c to be present in 57.6% of clade 8 strains. This SNP genotyping study also found a significant relationship between the presence of stx2 in conjunction with stx2c among strains of clade 8, in that no other clade associated with human illness displayed this combination. This is not surprising as lineage I strains only rarely contain stx2c despite their high association with human disease and lineage II strains nearly always contain stx2c despite a rare association with human disease [6]. It is unlikely that the combination of stx2 and stx2c alone is the reason for the hyper-virulence of lineage I/II strains, as the presence of stx2c is nearly ubiquitous among bovine-associated lineage II strains. The SNP genotyping study cited above only considered isolates associated with human disease; therefore it is likely that few lineage II strains were included in the study. A study by Friedrich et al. [32] examining stx2 subtypes and their association with clinical symptoms found stx2c to be the only subtype besides stx2 present in strains isolated from cases of HUS, but found no correlation between the presence of stx2c and the development of HUS. It has recently been shown that the level of Stx2 production is greater in lineage I strains than lineage II strains [33] so it may be that lineage I/II strains implicated in cases of HUS simply produce more toxin than other O157:H7 strains. However, it is possible that other factors possessed by the hyper-virulent lineage I/II group strains are responsible for their greater virulence in humans, and remain to be discovered.

The findings in this study highlight the need for a common genotyping approach, as it is evident that the same groups of genetically related strains have been given multiple designations based on the use of different comparative genotyping methods. This need exists for epidemiological studies of outbreak strains, where strain discrimination is the primary focus, as well as for population genetic studies where genomic information is of central importance. The value of being able to compare a pattern produced in a particular laboratory to one found in a national central database has made PFGE and PulseNet very useful in tracking outbreaks that are widely disseminated [34], despite the fact that PFGE is labour intensive and difficult to standardize [35]. As this common approach in identifying outbreak strains has been useful for a pattern-based method such as PFGE, approaches based on a multi-locus sampling of the genome could take this centralized database concept and extend it to contain presence/absence data for specific loci, SNPs, and other measures of heterogeneity. In this way, whether the goal is identifying an outbreak source or using information for a population-based study, the data would be available in a central repository. This type of system would be important in monitoring and identifying the emergence of new clones of O157:H7, such as the hyper-virulent lineage I/II/clade 8 strains and recognizing other changes in the population when genotyping E. coli O157:H7 strains associated with disease outbreaks.

The availability of whole genome sequence information has led to the development of new genotyping methods that are easier to perform than traditional methods, more discriminatory and more informative with respect to genotype and phenotype. It is interesting that the approaches targeting genetic polymorphisms within conserved genes and those targeting genetic changes based on gene insertion/deletion events converge to give a similar picture of E. coli O157:H7 strain relationships. Such concordance of methods has been previously demonstrated with mCGH and MLST for Campylobacter jejuni [36] and Streptococcus pneumoniae [37]. However, given that methods differ greatly in terms of the time required for the analysis, labour and equipment required, need for expertise, freedom from subjectivity in interpretation of the data and portability of the genotyping results from one laboratory to another, there is considerable advantage in selecting a method that is simple, extensible and easily portable. While some typing methods are better than others in these aspects, most of the multi-locus typing methods examined produced a similar tree architecture. This suggests a "common typing language" is possible at least in the context that genotypes derived using different methods can be integrated and communicated in a broader framework. While it was shown that multiple methods converge to provide a similar picture of the O157:H7 population structure, we are not advocating the routine use of multiple genotyping methods but rather the use of methods based on comparative genomics.

With the genomic sequencing revolution well under way, the ability to harvest novel sequence information in a timely fashion from new genomic sequences will become increasingly important and the ability to include and compare in silico results to those from traditional laboratory experiments will become necessary.

The results of this study suggest that genotyping approaches based on common comparative genomic data are likely to form the basis for the next-generation of analytical tools used for both population-based comparative genotyping and epidemiological studies.

Conclusion

A combined analysis of several genotyping methods and supernetwork construction confirm that the E. coli O157:H7 population is distributed among three major lineages. Recent work using SNP genotyping has identified an emergent hyper-virulent clade of public health importance among nine SNP-based clades of O157:H7. However, due to the lack of a common genotyping language, this designation has not been relatable to other genotype classification schemes. Using in silico analyses of multiple genome-sequenced strains we were able to integrate the overall diversity derived from individual methods into a more complete picture of the O157:H7 population structure and show that based on currently available O157:H7 genome sequence data, lineage I/II and clade 8 strains are synonymous. Further, some factor or factors other than the simple presence of stx2c are responsible for hyper-virulence among strains of this group. We have also previously shown that all PT2 strains examined belonged to lineage I/II. Future molecular typing and epidemiological studies should be aimed at providing genotyping data based on comparative genomic analyses. These data should be stored centrally and accessed locally in an easily transferable, relatable and extensible format that will form the basis of a common genotyping language.

Methods

Genome-sequenced strains used in the analyses

Our analyses utilized the complete and fully annotated genomes of K12 MG1655 and E. coli O157:H7 strains EDL933 and Sakai. Additionally, 14 E. coli O157:H7 whole-genome shotgun sequences from GenBank were included in the study, along with two whole-genome shotgun sequences that were sequenced in-house: EC71074, an E. coli O157:H7 strain and EC33264, an E. coli O145:NM strain (Table 1).

In silico lineage typing

The LSPA6 primer sequences [8] were utilized for in silico polymerase chain reaction (is-PCR) experiments to determine the lineage types of the strains in Table 1. The is-PCR method involved BLASTN searches [38] of the primer sequences combined with data processing in Microsoft Excel to determine the expected band size; any primer sequence with less than 80% sequence identity in a BLASTN search was considered absent.

Shiga-toxin encoding bacteriophage insertion site typing

As described by Besser et al. [9] the primer sequences targeting regions of bacteriophage insertion were subjected to is-PCR, using a sequence identity of 80% as a positive threshold. Data were given discrete character scores, with 0 representing absence and 1 representing an amplified band. In cases where bands of more than one size were found for the same primer set, the digits 2 through 5 were used to denote bands of additional size.

Comparative genomic fingerprinting (CGF)

Twenty-three of the most variable loci among E. coli O157:H7 strains were targeted by the comparative genomic fingerprinting method of Laing et al. [13]. The 80 strain binary dataset from the Laing et al. study was used, along with is-PCR results from the strains in Table 1, using a sequence identity of 80% as a positive threshold. The is- PCR data were converted to binary characters, with 0 for absence and 1 for presence.

Genomic in silico subtractive hybridization (GISSH)

The method for determining the existence of 417 genomic regions larger than 500 bp that are absent from the genomes of EDL933 and Sakai but present in other E. coli O157:H7 whole genome shotgun sequences will be described in detail elsewhere (Laing et al. in preparation). Each of the novel regions was split into contiguous segments of 500 bp, creating 1456 segments representing the original set of novel regions.

BLASTN was used to determine the distribution and these 1456 novel regions among the strains in Table 1. A segment containing less than 80% total sequence identity was considered absent and denoted by a 0, while a region containing greater than or equal to 80% total sequence identity was considered present and denoted by a 1.

Single nucleotide polymorphism (SNP) analysis

Manning et al. [12] recently described the existence of 96 SNPs among E. coli O157:H7 strains. The sequence of each SNP-containing locus in E. coli O157:H7 strain Sakai was used to conduct BLASTN comparisons of all 96 SNPs in each genomic sequence. Data were scored as single-letter nucleotides for either version of the SNP as put forth by Manning et al. with the following exception: in the Manning et al. study, E. coli O157:H7 strain Sakai locus ECs0606 is listed as having a SNP of A117G, whereas our data analysis found the SNP to be A117C.

Microarray comparative genomic hybridization (mCGH)

Zhang et al. [7] conducted a comprehensive microarray analysis of 31 E. coli O157:H7 strains using probe hybridization data for 6057 open reading frames. Analysis of those data found that the 50 mer probes used in the microarray experiments shared sequence identity of at least 80% to one or more of the three control strains (K12 MG1655, EDL933 and Sakai) in 6021 of the 6057 probes. The positive threshold used for in silico microarray hybridization was 80% and therefore the 36 ORFs with probes having less than 80% homology in all three control strains were not used. The in silico microarray hybridization used BLASTN to test the presence of each ORF in every strain in Table 1. Probes with identity greater than or equal to 80% were scored as present and everything else as divergent. Results were converted to a binary format for subsequent analyses with 1 representing a present ORF and 0 a divergent one.

Zhang et al. used the used the same MWG probe set as Wick et al. [11], which was used to determine the stepwise emergence of O157:H7 strains from an O55:H7-like ancestor. In their study, Wick et al. included two sorbitol-fermenting strains (CB2755 and 493/89) and three O55:H7 strains (Dec 5d, TB182A and 5905) that were absent from the study by Zhang et al. To expand the scope of our study, binary data for these five strains were included, with any loci not included in the Wick et al. study scored as missing and denoted by "?".

Multi-locus variable number tandem repeat analysis (MLVA)

The MLVA PCR scheme described by Hyattia-Trees et al. [10] was used in silico to identify the number of repeats at each locus, which is based on the number of repeats at nine variable loci within the O157:H7 genome. In accordance with the reference, partial repeats were rounded to the nearest whole number to generate the final set of data.

Construction of dendrograms

The character data for each of the six methods were converted to Nexus format [39] and imported to PAUP* v4.0 [40], where the data were analyzed using maximum parsimony. The tree space was thoroughly sampled using 1000 random sequence additions, with 10 trees kept per search using the following PAUP* command: "hsearch enforce = no start = stepwise addseq = random nreps = 1000 nchuck = 5 chuckscore = 1". The maximum parsimony tree for each set of data was saved and visualized in SplitsTree v4.0, where K12 was used as the outgroup.

Construction of supernetworks

The six maximum parsimony trees representing the six datasets comprised of the 19 strains in Table 1 were imported to Splitstree v4.0 and a supernetwork was created using Z-closure with 1000 iterations and the un-weighted mean distance of each tree to generate an equal angle representation of the data. A second supernetwork was created in an analogous fashion; however the CGF and microarray trees included data from the previously published studies of Laing et al [13], Wick et al. [11] and Zhang et al. [7].

References

  1. Siegler RL: The hemolytic uremic syndrome. Pediatr Clin North Am. 1995, 42: 1505-29.

    CAS  PubMed  Google Scholar 

  2. Yoon JW, Hovde CJ: All blood, No stool: enterohemorrhagic Escherichia coli O157:H7 infection. J Vet Sci. 2008, 9: 219-31.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Ohnishi M, Kurokawa K, Hayashi T: Diversification of Escherichia coli genomes: are bacteriophages the major contributors?. Trends Microbiol. 2001, 9: 481-5. 10.1016/S0966-842X(01)02173-4.

    Article  CAS  PubMed  Google Scholar 

  4. Allison HE: Stx-phages: drivers and mediators of the evolution of STEC and STEC-like pathogens. Future Microbiol. 2007, 2: 165-74. 10.2217/17460913.2.2.165.

    Article  CAS  PubMed  Google Scholar 

  5. Miao EA, Miller SI: Bacteriophages in the evolution of pathogen-host interactions. Proc Natl Acad Sci USA. 1999, 96: 9452-4. 10.1073/pnas.96.17.9452.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Kim J, Nietfeldt J, Benson AK: Octamer-based genome scanning distinguishes a unique subpopulation of Escherichia coli O157:H7 strains in cattle. Proc Natl Acad Sci USA. 1999, 96: 13288-93. 10.1073/pnas.96.23.13288.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Zhang Y, Laing C, Steele M, Ziebell K, Johnson R, Benson A, Taboada E, Gannon V: Genome evolution in major Escherichia coli O157:H7 lineages. BMC Genomics. 2007, 8: 121-10.1186/1471-2164-8-121.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Yang Z, Kovar J, Kim J, Nietfeldt J, Smith DR, Moxley RA, Olson ME, Fey PD, Benson AK: Identification of common subpopulations of non-sorbitol-fermenting, beta-glucuronidase-negative Escherichia coli O157:H7 from bovine production environments and human clinical samples. Appl Environ Microbiol. 2004, 70: 6846-54. 10.1128/AEM.70.11.6846-6854.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Besser TE, Shaikh N, Holt NJ, Tarr PI, Konkel ME, Malik-Kale P, Walsh CW, Whittam TS, Bono JL: Greater diversity of Shiga toxin-encoding bacteriophage insertion sites among Escherichia coli O157:H7 isolates from cattle than in those from humans. Appl Environ Microbiol. 2007, 73: 671-679. 10.1128/AEM.01035-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Hyytiä-Trees E, Smole SC, Fields PA, Swaminathan B, Ribot EM: Second generation subtyping: a proposed PulseNet protocol for multiple-locus variable-number tandem repeat analysis of Shiga toxin-producing Escherichia coli O157 (STEC O157). Foodborne Pathog Dis. 2006, 3: 118-31. 10.1089/fpd.2006.3.118.

    Article  PubMed  Google Scholar 

  11. Wick LM, Qi W, Lacher DW, Whittam TS: Evolution of genomic content in the stepwise emergence of Escherichia coli O157:H7. J Bacteriol. 2005, 187: 1783-1791. 10.1128/JB.187.5.1783-1791.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Manning SD, Motiwala AS, Springman AC, Qi W, Lacher DW, Ouellette LM, Mladonicky JM, Somsel P, Rudrik JT, Dietrich SE, Zhang W, Swaminathan B, Alland D, Whittam TS: Variation in virulence among clades of Escherichia coli O157:H7 associated with disease outbreaks. Proc Natl Acad Sci USA. 2008, 105: 4868-73. 10.1073/pnas.0710834105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Laing C, Pegg C, Yawney D, Ziebell K, Steele M, Johnson R, Thomas JE, Taboada EN, Zhang Y, Gannon VPJ: Rapid determination of Escherichia coli O157:H7 lineage types and molecular subtypes by using comparative genomic fingerprinting. Appl Environ Microbiol. 2008, 74: 6606-15. 10.1128/AEM.00985-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Feng PCH, Monday SR, Lacher DW, Allison L, Siitonen A, Keys C, Eklund M, Nagano H, Karch H, Keen J, Whittam TS: Genetic diversity among clonal lineages within Escherichia coli O157:H7 stepwise evolutionary model. Emerging Infect Dis. 2007, 13: 1701-6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Riordan JT, Viswanath SB, Manning SD, Whittam TS: Genetic differentiation of Escherichia coli O157:H7 clades associated with human disease by real-time PCR. J Clin Microbiol. 2008, 46: 2070-2073. 10.1128/JCM.00203-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-67. 10.1093/molbev/msj030.

    Article  CAS  PubMed  Google Scholar 

  17. Abe H, Miyahara A, Oshima T, Tashiro K, Ogura Y, Kuhara S, Ogasawara N, Hayashi T, Tobe T: Global regulation by horizontally transferred regulators establishes the pathogenicity of Escherichia coli. DNA Res. 2008, 15: 25-38. 10.1093/dnares/dsm033.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Hacker J, Kaper JB: Pathogenicity islands and the evolution of microbes. Annu Rev Microbiol. 2000, 54: 641-79. 10.1146/annurev.micro.54.1.641.

    Article  CAS  PubMed  Google Scholar 

  19. Ogura Y, Ooka T, Asadulghani , Terajima J, Nougayrède J, Kurokawa K, Tashiro K, Tobe T, Nakayama K, Kuhara S, Oswald E, Watanabe H, Hayashi T: Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes. Genome Biol. 2007, 8: R138-10.1186/gb-2007-8-7-r138.

    Article  PubMed Central  PubMed  Google Scholar 

  20. O'Brien AD, Newland JW, Miller SF, Holmes RK, Smith HW, Formal SB: Shiga-like toxin-converting phages from Escherichia coli strains that cause hemorrhagic colitis or infantile diarrhea. Science. 1984, 226: 694-6. 10.1126/science.6387911.

    Article  PubMed  Google Scholar 

  21. Plunkett G, Rose DJ, Durfee TJ, Blattner FR: Sequence of Shiga toxin 2 phage 933W from Escherichia coli O157:H7: Shiga toxin as a phage late-gene product. J Bacteriol. 1999, 181: 1767-78.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Neely MN, Friedman DI: Functional and genetic analysis of regulatory regions of coliphage H-19B: location of shiga-like toxin and lysis genes suggest a role for phage functions in toxin release. Mol Microbiol. 1998, 28: 1255-67. 10.1046/j.1365-2958.1998.00890.x.

    Article  CAS  PubMed  Google Scholar 

  23. McDaniel TK, Jarvis KG, Donnenberg MS, Kaper JB: A genetic locus of enterocyte effacement conserved among diverse enterobacterial pathogens. Proc Natl Acad Sci USA. 1995, 92: 1664-8. 10.1073/pnas.92.5.1664.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Tobe T, Beatson SA, Taniguchi H, Abe H, Bailey CM, Fivian A, Younis R, Matthews S, Marches O, Frankel G, Hayashi T, Pallen MJ: An extensive repertoire of type III secretion effectors in Escherichia coli O157 and the role of lambdoid phages in their dissemination. Proc Natl Acad Sci USA. 2006, 103: 14941-6. 10.1073/pnas.0604891103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Karmali MA, Mascarenhas M, Shen S, Ziebell K, Johnson S, Reid-Smith R, Isaac-Renton J, Clark C, Rahn K, Kaper JB: Association of genomic O island 122 of Escherichia coli EDL 933 with verocytotoxin-producing Escherichia coli seropathotypes that are linked to epidemic and/or serious disease. J Clin Microbiol. 2003, 41: 4930-40. 10.1128/JCM.41.11.4930-4940.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Shaikh N, Tarr PI: Escherichia coli O157:H7 Shiga toxin-encoding bacteriophages: integrations, excisions, truncations, and evolutionary implications. J Bacteriol. 2003, 185: 3596-605. 10.1128/JB.185.12.3596-3605.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Kelly BG, Vespermann A, Bolton DJ: The role of horizontal gene transfer in the evolution of selected foodborne bacterial pathogens. Food Chem Toxicol. 2009, 47: 951-968. 10.1016/j.fct.2008.02.006.

    Article  CAS  PubMed  Google Scholar 

  28. Kudva IT, Evans PS, Perna NT, Barrett TJ, Ausubel FM, Blattner FR, Calderwood SB: Strains of Escherichia coli O157:H7 differ primarily by insertions or deletions, not single-nucleotide polymorphisms. J Bacteriol. 2002, 184: 1873-9. 10.1128/JB.184.7.1873-1879.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Ziebell K, Steele M, Zhang Y, Benson A, Taboada EN, Laing C, McEwen S, Ciebin B, Johnson R, Gannon V: Genotypic characterization and prevalence of virulence factors among Canadian Escherichia coli O157:H7 strains. Appl Environ Microbiol. 2008, 74: 4314-23. 10.1128/AEM.02821-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Steele M, Ziebell K, Zhang Y, Benson A, Konczy P, Johnson R, Gannon V: Identification of Escherichia coli O157:H7 genomic regions conserved in strains with a genotype associated with human infection. Appl Environ Microbiol. 2007, 73: 22-31. 10.1128/AEM.00982-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Dowd S, Ishizaki H: Microarray based comparison of two Escherichia coli O157:H7 lineages. BMC Microbiology. 2006, 6: 30-10.1186/1471-2180-6-30.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Friedrich AW, Bielaszewska M, Zhang W, Pulz M, Kuczius T, Ammon A, Karch H: Escherichia coli harboring Shiga toxin 2 gene variants: frequency and association with clinical symptoms. J Infect Dis. 2002, 185: 74-84. 10.1086/338115.

    Article  CAS  PubMed  Google Scholar 

  33. Dowd SE, Williams JB: Comparison of Shiga-like toxin II expression between two genetically diverse lineages of Escherichia coli O157:H7. J Food Prot. 2008, 71: 1673-8.

    CAS  PubMed  Google Scholar 

  34. Cooley M, Carychao D, Crawford-Miksza L, Jay MT, Myers C, Rose C, Keys C, Farrar J, Mandrell RE: Incidence and tracking of Escherichia coli O157:H7 in a major produce production region in California. PLoS ONE. 2007, 2: e1159-10.1371/journal.pone.0001159.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Murchan S, Kaufmann ME, Deplano A, de Ryck R, Struelens M, Zinn CE, Fussing V, Salmenlinna S, Vuopio-Varkila J, El Solh N, Cuny C, Witte W, Tassios PT, Legakis N, van Leeuwen W, van Belkum A, Vindel A, Laconcha I, Garaizar J, Haeggman S, Olsson-Liljequist B, Ransjo U, Coombes G, Cookson B: Harmonization of pulsed-field gel electrophoresis protocols for epidemiological typing of strains of methicillin-resistant Staphylococcus aureus : a single approach developed by consensus in 10 European laboratories and its application for tracing the spread of related strains. J Clin Microbiol. 2003, 41: 1574-85. 10.1128/JCM.41.4.1574-1585.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Taboada EN, Mackinnon JM, Luebbert CC, Gannon VPJ, Nash JHE, Rahn K: Comparative genomic assessment of Multi-Locus Sequence Typing: rapid accumulation of genomic heterogeneity among clonal isolates of Campylobacter jejuni. BMC Evol Biol. 2008, 8: 229-10.1186/1471-2148-8-229.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Dagerhamn J, Blomberg C, Browall S, Sjöström K, Morfeldt E, Henriques-Normark B: Determination of accessory gene patterns predicts the same relatedness among strains of Streptococcus pneumoniae as sequencing of housekeeping genes does and represents a novel approach in molecular epidemiology. J Clin Microbiol. 2008, 46: 863-8. 10.1128/JCM.01438-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Maddison DR, Swofford DL, Maddison WP: NEXUS: an extensible file format for systematic information. Syst Biol. 1997, 46: 590-621.

    Article  CAS  PubMed  Google Scholar 

  40. Wilgenbusch JC, Swofford D: Inferring evolutionary trees with PAUP*. Curr Protoc Bioinformatics. 2003, Chapter 6 (Unit 6.4):

Download references

Acknowledgements

We would like to thank the Canadian Food Inspection Agency (CFIA) for allowing this research to be conducted at the Animal Diseases Research Institute (ADRI), Lethbridge, Alberta, Canada. This research was funded by the Public Health Agency of Canada.

Author information

Authors and Affiliations

  1. Laboratory for Foodborne Zoonoses, Public Health Agency of Canada, Lethbridge, AB, Canada

    Chad R Laing, Cody Buchanan, Eduardo N Taboada, Yongxiang Zhang & Victor PJ Gannon

  2. Laboratory for Foodborne Zoonoses, Public Health Agency of Canada, Guelph, ON, Canada

    Mohamed A Karmali

  3. Faculty of Biological Sciences, University of Lethbridge, Lethbridge, AB, Canada

    James E Thomas

Authors
  1. Chad R Laing
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Cody Buchanan
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Eduardo N Taboada
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yongxiang Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mohamed A Karmali
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. James E Thomas
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Victor PJ Gannon
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Victor PJ Gannon.

Additional information

Authors' contributions

CRL planned and carried out experimental work, data analysis and writing of the manuscript. CB carried out experimental work. ENT helped with data analysis and with revision of the manuscript. YZ planned experimental work and helped with revision of the manuscript. MAK assisted with project planning and revision of the manuscript. JET helped with project planning and revision of the manuscript. VPJG planned experimental work, contributed to data analysis and with preparation of the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

12864_2008_2171_MOESM1_ESM.xls

Additional file 1: In silico data files. A Microsoft Excel workbook containing the dataset for each of the six in silico typing methods. (XLS 1 MB)

12864_2008_2171_MOESM2_ESM.pdf

Additional file 2: Maximum parsimony trees from the in silico analyses. The maximum parsimony trees created from the in silico datasets, which were used in the creation of Figures 2 and 3. (PDF 265 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Laing, C.R., Buchanan, C., Taboada, E.N. et al. In silico genomic analyses reveal three distinct lineages of Escherichia coli O157:H7, one of which is associated with hyper-virulence. BMC Genomics 10, 287 (2009). https://doi.org/10.1186/1471-2164-10-287

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-10-287

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us