Variable number of tandem repeats and pulsed-field gel electrophoresis cluster analysis of enterohemorrhagic Escherichia coli serovar O157 strains

Application of whole genome sequence data in analyzing the molecular epidemiology of Shiga toxin-producing Escherichia coli O157:H7/H

Seventeen clusters of Shiga toxin-producing Escherichia coli O157:H7/- (O157) strains, determined by cluster analysis of pulsed-field gel electrophoresis patterns, were analyzed using whole genome sequence (WGS) data to investigate this pathogen's molecular epidemiology. The 17 clusters included 136 strains containing strains from nine outbreaks, with each outbreak caused by a single source contaminated with the organism, as shown by epidemiological contact surveys. WGS data of these strains were used to identify single nucleotide polymorphisms (SNPs) by two methods: short read data were directly mapped to a reference genome (mapping derived SNPs) and common SNPs between the mapping derived SNPs and SNPs in assembled data of short read data (common SNPs). Among both SNPs, those that were detected in genes with a gap were excluded to remove ambiguous SNPs from further analysis. The effectiveness of both SNPs was investigated among all the concatenated SNPs that were detected (whole SNP set); SNPs were divided into three categories based on the genes in which they were located (i.e., backbone SNP set, O-island SNP set, and mobile element SNP set); and SNPs in non-coding regions (intergenic region SNP set). When SNPs from strains isolated from the nine single source derived outbreaks were analyzed using an unweighted pair group method with arithmetic mean tree (UPGMA) and a minimum spanning tree (MST), the maximum pair-wise distances of the backbone SNP set of the mapping derived SNPs were significantly smaller than those of the whole and intergenic region SNP set on both UPGMAs and MSTs. This significant difference was also observed when the backbone SNP set of the common SNPs were examined (Steel-Dwass test, P≤0.01). When the maximum pair-wise distances were compared between the mapping derived and common SNPs, significant differences were observed in those of the whole, mobile element, and intergenic region SNP set (Wilcoxon signed rank test, P≤0.01). When all the strains included in one complex on an MST or one cluster on a UPGMA were designated as the same genotype, the values of the Hunter-Gaston Discriminatory Power Index for the backbone SNP set of the mapping derived and common SNPs were higher than those of other SNP sets. In contrast, the mobile element SNP set could not robustly subdivide lineage I strains of tested O157 strains using both the mapping derived and common SNPs. These results suggested that the backbone SNP set were the most effective for analysis of WGS data for O157 in enabling an appropriation of its molecular epidemiology.

Emergence of enterohemorrhagic Escherichia coli serovar O157 strains in clade 8 with highly similar pulsed-field gel electrophoresis patterns

Enterohemorrhagic Escherichia coli serovar O157 (O157) strains with highly similar pulsed-field gel electrophoresis (PFGE) patterns were isolated in Japan during 2007 and 2008. Several genetic features related to O157 evolution were investigated to indicate whether homoplasy might have contributed to the highly similar PFGE patterns in these strains. The O157 strains were classified in lineage I/II, as defined by a lineage-specific polymorphism assay-6 with an atypical allele in Z5935 (code: 231111). Analysis of the insertion sites of stx(2) phage in these strains showed that the sites were "occupied" in yehV and "intact" in wrbA, indicating that the strains were derived from "Cluster 1" of "Subgroup C." When a specific single-nucleotide polymorphism in ECs2357 in clade 8 strains was investigated, all of the strains in the present study were confirmed to be clade 8 strains. These results indicated that the O157 strains in this study had common genetic features, suggesting that the highly similar PFGE patterns of these strains were not due to homoplasy. Because no common source of these strains could be identified in 2007 to 2008 in Japan, these strains may have emerged from a unique O157 clade 8 clone and then spread by dissemination in Japan.