A Geographically Widespread Outbreak Investigation and Development of a Rapid Screening Method Using Whole Genome Sequences of Enterohemorrhagic Escherichia coli O121

Increased resistance against tellurite is conferred by a mutation in the promoter region of uncommon tellurite resistance gene tehB in the ter-negative Shiga toxin-producing Escherichia coli O157:H7

Resistance to potassium tellurite (PT) is an important indicator in isolating Shiga toxin-producing Escherichia coli (STEC) O157:H7 and other major STEC serogroups. Common resistance determinant genes are encoded in the ter gene cluster. We found an O157:H7 isolate that does not harbor ter but is resistant to PT. One nonsynonymous mutation was found in another PT resistance gene, tehA, through whole-genome sequence analyses. To elucidate the contribution of this mutation to PT resistance, complementation of tehA and the related gene tehB in isogenic strains and quantitative RT‒PCR were performed. The results indicated that the point mutation not only changed an amino acid of tehA, but also was positioned on a putative internal promoter of tehB and increased PT resistance by elevating tehB mRNA expression. Meanwhile, the amino acid change in tehA had negligible impact on the PT resistance. Comprehensive screening revealed that 2.3% of O157:H7 isolates in Japan did not harbor the ter gene cluster, but the same mutation in tehA was not found. These results suggested that PT resistance in E. coli can be enhanced through one mutational event even in ter-negative strains.

IMPORTANCE

Selective agents are important for isolating Shiga toxin-producing Escherichia coli (STEC) because the undesirable growth of microflora should be inhibited. Potassium tellurite (PT) is a common selective agent for major STEC serotypes. In this study, we found a novel variant of PT resistance genes, tehAB, in STEC O157:H7. Molecular experiments clearly showed that one point mutation in a predicted internal promoter region of tehB upregulated the expression of the gene and consequently led to increased resistance to PT. Because tehAB genes are ubiquitous across E. coli, these results provide universal insight into PT resistance in this species.

Combined usage of serodiagnosis and O antigen typing to isolate Shiga toxin-producing Escherichia coli O76:H7 from a hemolytic uremic syndrome case and genomic insights from the isolate

IMPORTANCE

Hemolytic uremic syndrome (HUS) is a life-threatening disease caused by Shiga toxin-producing Escherichia coli (STEC) infection. The treatment approaches for STEC-mediated typical HUS and atypical HUS differ, underscoring the importance of rapid and accurate diagnosis. However, specific detection methods for STECs other than major serogroups, such as O157, O26, and O111, are limited. This study focuses on the utility of PCR-based O-serotyping, serum agglutination tests utilizing antibodies against the identified Og type, and isolation techniques employing antibody-conjugated immunomagnetic beads for STEC isolation. By employing these methods, we successfully isolated a STEC strain of a minor serotype, O76:H7, from a HUS patient.

Comparative genomics of Nocardia seriolae reveals recent importation and subsequent widespread dissemination in mariculture farms in the South Central Coast region, Vietnam

Between 2010 and 2015, nocardiosis outbreaks caused by Nocardia seriolae affected many permit farms throughout Vietnam, causing mass fish mortalities. To understand the biology, origin and epidemiology of these outbreaks, 20 N. seriolae strains collected from farms in four provinces in the South Central Coast region of Vietnam, along with two Taiwanese strains, were analysed using genetics and genomics. PFGE identified a single cluster amongst all Vietnamese strains that was distinct from the Taiwanese strains. Like the PFGE findings, phylogenomic and SNP genotyping analyses revealed that all Vietnamese N. seriolae strains belonged to a single, unique clade. Strains fell into two subclades that differed by 103 SNPs, with almost no diversity within clades (0–5 SNPs). There was no association between geographical origin and subclade placement, suggesting frequent N. seriolae transmission between Vietnamese mariculture facilities during the outbreaks. The Vietnamese strains shared a common ancestor with strains from Japan and China, with the closest strain, UTF1 from Japan, differing by just 220 SNPs from the Vietnamese ancestral node. Draft Vietnamese genomes range from 7.55 to 7.96 Mbp in size, have an average G+C content of 68.2 % and encode 7 602–7958 predicted genes. Several putative virulence factors were identified, including genes associated with host cell adhesion, invasion, intracellular survival, antibiotic and toxic compound resistance, and haemolysin biosynthesis. Our findings provide important new insights into the epidemiology and pathogenicity of N. seriolae and will aid future vaccine development and disease management strategies, with the ultimate goal of nocardiosis-free aquaculture.

Genome Sequence Analysis and Characterization of Shiga Toxin 2 Production by Escherichia coli O157:H7 Strains Associated With a Laboratory Infection

A laboratory-acquired E. coli O157:H7 infection with associated severe sequelae including hemolytic uremic syndrome occurred in an individual working in the laboratory with a mixture of nalidixic acid-resistant (Nal^R) O157:H7 mutant strains in a soil-biochar blend. The patient was hospitalized and treated with an intravenous combination of metronidazole and levofloxacin. The present study investigated the source of this severe laboratory acquired infection and further examined the influence of the antibiotics used during treatment on the expression and production of Shiga toxin. Genomes of two Stx_2a-and eae-positive O157:H7 strains isolated from the patient's stool were sequenced along with two pairs of the wt strains and their derived Nal^R mutants used in the laboratory experiments. High-resolution SNP typing determined the strains' individual genetic relatedness and unambiguously identified the two laboratory-derived Nal^R mutant strains as the source of the researcher's life-threatening disease, rather than a conceivable ingestion of unrelated O157:H7 isolates circulating at the same time. It was further confirmed that in sublethal doses, the antibiotics increased toxin expression and production. Our results support a simultaneous co-infection with clinical strains in the laboratory, which were the causative agents of previous O157:H7 outbreaks, and further that the administration of antibiotics may have impacted the outcome of the infection.

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

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

Whole Genome Analysis Detects the Emergence of a Single Salmonella enterica Serovar Chester Clone in Japan's Kanto Region

In Japan's Kanto region, the number of Salmonella enterica serovar Chester infections increased temporarily between 2014 and 2016. Concurrently with this temporal increase in the Kanto region, S. Chester isolates belonging to one clonal group were causing repetitive outbreaks in Europe. A recent study reported that the European outbreaks were associated with travelers who had been exposed to contaminated food in Morocco, possibly seafood. Because Japan imports a large amount of seafood from Morocco, we aimed to establish whether the temporal increase in S. Chester infections in the Kanto region was associated with imported Moroccan seafood. Short sequence reads from the whole-genome sequencing of 47 S. Chester isolates from people in the Kanto region (2014-2016), and the additional genome sequences from 58 isolates from the European outbreaks, were analyzed. The reads were compared with the complete genome sequence from a S. Chester reference strain, and 347 single nucleotide polymorphisms (SNPs) were identified. These SNPs were used in this study. Cluster and Bayesian cluster analyses showed that the Japanese and European isolates fell into two different clusters. Therefore, Φ PT and I A S values were calculated to evaluate genetic differences between these clusters. The results revealed that the Japanese and European isolates were genetically distinct populations. Our root-to-tip analysis showed that the Japanese isolates originating from one clone had accumulated mutations, suggesting that an emergence of this organism occurred. A minimum spanning tree analysis demonstrated no correlation between genetic and geographical distances in the Japanese isolates, suggesting that the emergence of the serovar in the Kanto region did not involve person-to-person contact; rather, it occurred through food consumption. The d N /d S ratio indicated that the Japanese strain has evolved under positive selection pressure. Generally, a population of bacterial clones in a reservoir faces negative selection pressure. Therefore, the Japanese strain must have existed outside of any reservoir during its emergence. In conclusion, S. Chester isolates originating from one clone probably emerged in the Kanto region via the consumption of contaminated foods other than imported Moroccan seafood. The emerging strain may have not established a reservoir for survival in the food supply chain resulting in its disappearance after 2017.

Whole-Genome Sequencing of Shiga Toxin-Producing Escherichia coli OX18 from a Fatal Hemolytic Uremic Syndrome Case

We report a fatal case of hemolytic uremic syndrome with urinary tract infection in Japan caused by Shiga toxin–producing Escherichia coli. We genotypically identified the isolate as OX18:H2. Whole-genome sequencing revealed 3 potentially pathogenic lineages (OX18:H2, H19, and H34) that have been continuously isolated in Japan.

Global distribution of epidemic-related Shiga toxin 2 encoding phages among enteroaggregative Escherichia coli

Since the Shiga toxin-producing enteroaggregative Escherichia coli (Stx-EAEC) O104:H4 strain caused a massive outbreak across Europe in 2011, the importance of Stx-EAEC has attracted attention from a public health perspective. Two Stx-EAEC O86 isolates were obtained from patients with severe symptoms in Japan in 1999 and 2015. To characterize the phylogeny and pathogenic potential of these Stx-EAEC O86 isolates, whole-genome sequence analyses were performed by short-and long-read sequencing. Among genetically diverse E. coli O86, the Stx-EAEC O86 isolates were clustered with the EAEC O86:H27 ST3570 subgroup. Strikingly, there were only two loci with single nucleotide polymorphisms (SNPs) between the Stx2a phage of a Japanese O86:H27 isolate and that of the European epidemic-related Stx-EAEC O104:H4 isolate. These results provide evidence of global distribution of epidemic-related Stx2a phages among various lineages of E. coli with few mutations.

Complete Genome Sequence of an Enterohemorrhagic Escherichia coli O111:H8 Strain Recovered from a Large Outbreak in Japan Associated with Consumption of Raw Beef

We present the complete genome sequence of an enterohemorrhagic Escherichia coli O111:H8 strain. This strain was isolated from a hemolytic-uremic syndrome patient and was responsible for a large outbreak associated with the consumption of raw beef in 2011.

We present the complete genome sequence of an enterohemorrhagic Escherichia coli O111:H8 strain. This strain was isolated from a hemolytic-uremic syndrome patient and was responsible for a large outbreak associated with the consumption of raw beef in 2011.

Effective Surveillance Using Multilocus Variable-Number Tandem-Repeat Analysis and Whole-Genome Sequencing for Enterohemorrhagic Escherichia coli O157

Due to the potential of enterohemorrhagic Escherichia coli (EHEC) serogroup O157 to cause large food borne outbreaks, national and international surveillance is necessary. For developing an effective method of molecular surveillance, a conventional method, multilocus variable-number tandem-repeat analysis (MLVA), and whole-genome sequencing (WGS) analysis were compared. WGS of 369 isolates of EHEC O157 belonging to 7 major MLVA types and their relatives were subjected to comprehensive in silico typing, core genome single nucleotide polymorphism (cgSNP), and core genome multilocus sequence typing (cgMLST) analyses. The typing resolution was the highest in cgSNP analysis. However, determination of the sequence of the mismatch repair protein gene mutS is necessary because spontaneous deletion of the gene could lead to a hypermutator phenotype. MLVA had sufficient typing resolution for a short-term outbreak investigation and had advantages in rapidity and high throughput. cgMLST showed less typing resolution than cgSNP, but it is less time-consuming and does not require as much computer power. Therefore, cgMLST is suitable for comparisons using large data sets (e.g., international comparison using public databases). In conclusion, screening using MLVA followed by cgMLST and cgSNP analyses would provide the highest typing resolution and improve the accuracy and cost-effectiveness of EHEC O157 surveillance.IMPORTANCE Intensive surveillance for enterohemorrhagic Escherichia coli (EHEC) serogroup O157 is important to detect outbreaks and to prevent the spread of the bacterium. Recent advances in sequencing technology made molecular surveillance using whole-genome sequence (WGS) realistic. To develop rapid, high-throughput, and cost-effective typing methods for real-time surveillance, typing resolution of WGS and a conventional typing method, multilocus variable-number tandem-repeat analysis (MLVA), was evaluated. Nation-level systematic comparison of MLVA, core genome single nucleotide polymorphism (cgSNP), and core genome multilocus sequence typing (cgMLST) indicated that a combination of WGS and MLVA is a realistic approach to improve EHEC O157 surveillance.

Enterohaemorrhagic Escherichia coli O121:H19 acquired an extended-spectrum β-lactamase gene during the development of an outbreak in two nurseries

Enterohaemorrhagic Escherichia coli (EHEC) is an important human pathogen worldwide. Although serotype O157 is currently the most dominant and important EHEC strain, serotypes O26, O111, O91, O103 and O121 are also recognized as serious pathogens that affect public health. EHEC outbreaks often occur in nurseries and elderly care facilities. In 2012, a nursery outbreak of EHEC O121 occurred during which the bacterium acquired a plasmid-borne extended-spectrum β-lactamase (ESBL) gene. ESBL-producing E. coli O86 was concurrently isolated from one of the EHEC patients. Therefore, we investigated the isolates by whole-genome sequence (WGS) analysis to elucidate the transmission dynamics of the EHEC strains and the ESBL plasmid. According to WGS-based phylogeny, all 17 EHEC O121 isolates were clonal, while E. coli O86 was genetically distant from the EHEC O121 isolates. The complete sequence of an ESBL plasmid encoding the CTX-M-55 β-lactamase was determined using S1-PFGE bands, and subsequent mapping of the WGS reads confirmed that the plasmid sequences from EHEC O121 and E. coli O86 were identical. Furthermore, conjugation experiments showed that the plasmid was capable of conjugative transfer. These results support the hypothesis that EHEC O121 acquired an ESBL-producing plasmid from E. coli O86 during the outbreak. This report demonstrates the importance of implementing preventive measures during EHEC outbreaks to control both secondary infection and the spread of antimicrobial resistance factors.

Polymorphism of Type I-F CRISPR/Cas system in Escherichia coli of phylogenetic group B2 and its application in genotyping

E. coli of phylogenetic group B2 is responsible for many extraintestinal infections, posing a great threat to health. The relatively polymorphic nature of CRISPR in phylogenetically related E. coli strains makes them potential markers for bacterial typing and evolutionary studies. In the current work, we investigated the occurrence and diversity of CRISPR/Cas system and explored its potential for genotyping. Type I-F CRISPR/Cas systems were found in 413 of 1190 strains of E. coli and exhibited the clustering within certain CCs and STs. And CRISPR spacer contents correlated well with MLST types. The divergence analysis of CRISPR showed stronger discriminatory power than MLST, and CRISPR polymorphism was instrumental for differentiating highly closely related strains. The timeline of spacer acquisition and deletion provided important information for inferring the evolution model between distinct serotypes. Identical spacer sequences were shared by strains with the same H-antigen type but not strains with the same O-antigen type. The homology between spacers and antibiotic-resistant plasmids demonstrated the role of Type I-F system in limiting the acquisition of antimicrobial resistance. Collectively, our data presents the dynamic nature of Type I-F CRISPR in E. coli of phylogenetic group B2 and provides new insights into the application of CRISPR-based typing in the species.

Characterization of OXA-48-like carbapenemase producers in Canada, 2011-14

Objectives

Since the first identification of the OXA-48 carbapenemase in 2001, Enterobacteriaceae harbouring OXA-48-like enzymes have been reported globally. Here, we applied WGS to characterize the molecular epidemiology of these bacterial isolates.

Methods

Enterobacteriaceae non-susceptible to carbapenems isolated from patients between 2011 and 2014 were voluntarily submitted to the Canadian National Microbiology Laboratory where they were screened for carbapenemase genes. WGS was conducted on OXA-48-like producers using the Illumina MiSeq platform. WGS data were used for single nucleotide variant (SNV) analysis, MLST analysis, detection of resistance genes and partial plasmid characterization. Susceptibilities were determined using Vitek2 and Etest. Patient data provided from sites were reviewed.

Results

Sixty-seven non-duplicated cases were identified among Escherichia coli (n = 21) and Klebsiella pneumoniae (n = 46). Recent international travel was observed in 40.4% of cases. OXA-181 (52.2%) and OXA-48 (31.3%) were the most common variants, one E. coli OXA-48 producer was found to harbour the acquired colistin resistance gene mcr-1. The dominant STs were ST38 and ST410 in E. coli and ST14 in K. pneumoniae. Three common plasmid types were observed among isolates: IncL/M associated with OXA-48 producers, and ColKP3 and IncX3 associated with OXA-181/232 producers.

Conclusions

Enterobacteriaceae with OXA-48-like carbapenemases are emerging in Canada. This study highlights the complexity of OXA-48-types identified in Canada owing to travel and the successful clones and plasmids harbouring the OXA-48-like enzyme.

Escherichia coli Sequence Type 410 Is Causing New International High-Risk Clones

Extraintestinal pathogenic Escherichia coli (ExPEC) is the main cause of urinary tract infections and septicemia. Significant attention has been given to the ExPEC sequence type ST131, which has been categorized as a “high-risk” clone. High-risk clones are globally distributed clones associated with various antimicrobial resistance determinants, ease of transmission, persistence in hosts, and effective transmission between hosts. The high-risk clones have enhanced pathogenicity and cause severe and/or recurrent infections. We show that clones of the E. coli ST410 lineage persist and/or cause recurrent infections in humans, including bloodstream infections. We found evidence of ST410 being a highly resistant globally distributed lineage, capable of patient-to-patient transmission causing hospital outbreaks. Our analysis suggests that the ST410 lineage should be classified with the potential to cause new high-risk clones. Thus, with the clonal expansion over the past decades and increased antimicrobial resistance to last-resort treatment options, ST410 needs to be monitored prospectively.

Escherichia coli sequence type 410 (ST410) has been reported worldwide as an extraintestinal pathogen associated with resistance to fluoroquinolones, third-generation cephalosporins, and carbapenems. In the present study, we investigated national epidemiology of ST410 E. coli isolates from Danish patients. Furthermore, E. coli ST410 was investigated in a global context to provide further insight into the acquisition of the carbapenemase genes bla_OXA-181 and bla_NDM-5 of this successful lineage. From 127 whole-genome-sequenced isolates, we reconstructed an evolutionary framework of E. coli ST410 which portrays the antimicrobial-resistant clades B2/H24R, B3/H24Rx, and B4/H24RxC. The B2/H24R and B3/H24Rx clades emerged around 1987, concurrently with the C1/H30R and C2/H30Rx clades in E. coli ST131. B3/H24Rx appears to have evolved by the acquisition of the extended-spectrum β-lactamase (ESBL)-encoding gene bla_CTX-M-15 and an IncFII plasmid, encoding IncFIA and IncFIB. Around 2003, the carbapenem-resistant clade B4/H24RxC emerged when ST410 acquired an IncX3 plasmid carrying a bla_OXA-181 carbapenemase gene. Around 2014, the clade B4/H24RxC acquired a second carbapenemase gene, bla_NDM-5, on a conserved IncFII plasmid. From an epidemiological investigation of 49 E. coli ST410 isolates from Danish patients, we identified five possible regional outbreaks, of which one outbreak involved nine patients with bla_OXA-181- and bla_NDM-5-carrying B4/H24RxC isolates. The accumulated multidrug resistance in E. coli ST410 over the past two decades, together with its proven potential of transmission between patients, poses a high risk in clinical settings, and thus, E. coli ST410 should be considered a lineage with emerging “high-risk” clones, which should be monitored closely in the future.

IMPORTANCE Extraintestinal pathogenic Escherichia coli (ExPEC) is the main cause of urinary tract infections and septicemia. Significant attention has been given to the ExPEC sequence type ST131, which has been categorized as a “high-risk” clone. High-risk clones are globally distributed clones associated with various antimicrobial resistance determinants, ease of transmission, persistence in hosts, and effective transmission between hosts. The high-risk clones have enhanced pathogenicity and cause severe and/or recurrent infections. We show that clones of the E. coli ST410 lineage persist and/or cause recurrent infections in humans, including bloodstream infections. We found evidence of ST410 being a highly resistant globally distributed lineage, capable of patient-to-patient transmission causing hospital outbreaks. Our analysis suggests that the ST410 lineage should be classified with the potential to cause new high-risk clones. Thus, with the clonal expansion over the past decades and increased antimicrobial resistance to last-resort treatment options, ST410 needs to be monitored prospectively.

Identification of Novel Biomarkers for Priority Serotypes of Shiga Toxin-Producing Escherichia coli and the Development of Multiplex PCR for Their Detection

It would be desirable to have an unambiguous scheme for the typing of Shiga toxin-producing Escherichia coli (STEC) isolates to subpopulations. Such a scheme should take the high genomic plasticity of E. coli into account and utilize the stratification of STEC into subgroups, based on serotype or phylogeny. Therefore, our goal was to identify specific marker combinations for improved classification of STEC subtypes. We developed and evaluated two bioinformatic pipelines for genomic marker identification from larger sets of bacterial genome sequences. Pipeline A performed all-against-all BLASTp analyses of gene products predicted in STEC genome test sets against a set of control genomes. Pipeline B identified STEC marker genes by comparing the STEC core proteome and the "pan proteome" of a non-STEC control group. Both pipelines defined an overlapping, but not identical set of discriminative markers for different STEC subgroups. Differential marker prediction resulted from differences in genome assembly, ORF finding and inclusion cut-offs in both workflows. Based on the output of the pipelines, we defined new specific markers for STEC serogroups and phylogenetic groups frequently associated with outbreaks and cases of foodborne illnesses. These included STEC serogroups O157, O26, O45, O103, O111, O121, and O145, Shiga toxin-positive enteroaggregative E. coli O104:H4, and HUS-associated sequence type (ST)306. We evaluated these STEC marker genes for their presence in whole genome sequence data sets. Based on the identified discriminative markers, we developed a multiplex PCR (mPCR) approach for detection and typing of the targeted STEC. The specificity of the mPCR primer pairs was verified using well-defined clinical STEC isolates as well as isolates from the ECOR, DEC, and HUSEC collections. The application of the STEC mPCR for food analysis was tested with inoculated milk. In summary, we evaluated two different strategies to screen large genome sequence data sets for discriminative markers and implemented novel marker genes found in this genome-wide approach into a DNA-based typing tool for STEC that can be used for the characterization of STEC from clinical and food samples.

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.