Clade analysis of enterohemorrhagic Escherichia coli serotype O157:H7/H- strains and hierarchy of their phylogenetic relationships

Pathogenomes of Shiga Toxin Positive and Negative Escherichia coli O157:H7 Strains TT12A and TT12B: Comprehensive Phylogenomic Analysis Using Closed Genomes

Shiga toxin-producing Escherichia coli are zoonotic pathogens that cause food-borne human disease. Among these, the O157:H7 serotype has evolved from an enteropathogenic O55:H7 ancestor through the displacement of the somatic gene cluster and recurrent toxigenic conversion by Shiga toxin-converting bacteriophages. However, atypical strains that lack the Shiga toxin, the characteristic virulence hallmark, are circulating in this lineage. For this study, we analyzed the pathogenome and virulence inventories of the stx+ strain, TT12A, isolated from a patient with hemorrhagic colitis, and its respective co-isolated stx- strain, TT12B. Sequencing the genomes to closure proved critical to the cataloguing of subtle strain differentiating sequence and structural polymorphisms at a high-level of phylogenetic accuracy and resolution. Phylogenomic profiling revealed SNP and MLST profiles similar to the near clonal outbreak isolates. Their prophage inventories, however, were notably different. The attenuated atypical non-shigatoxigenic status of TT12B is explained by the absence of both the ΦStx1a- and ΦStx2a-prophages carried by TT12A, and we also recorded further alterations in the non-Stx prophage complement. Phenotypic characterization indicated that culture growth was directly impacted by the strains' distinct lytic phage complement. Altogether, our phylogenomic and phenotypic analyses show that these intimately related isogenic strains are on divergent Stx(+/stx-) evolutionary paths.

Escherichia coli O157:H7 tir 255 T > A allele strains differ in chromosomal and plasmid composition

Shiga toxin-producing Escherichia coli (STEC) O157:H7 strains with the T allele in the translocated intimin receptor polymorphism (tir) 255 A > T gene associate with human disease more than strains with an A allele; however, the allele is not thought to be the direct cause of this difference. We sequenced a diverse set of STEC O157:H7 strains (26% A allele, 74% T allele) to identify linked differences that might underlie disease association. The average chromosome and pO157 plasmid size and gene content were significantly greater within the tir 255 A allele strains. Eighteen coding sequences were unique to tir 255 A allele chromosomes, and three were unique to tir 255 T allele chromosomes. There also were non-pO157 plasmids that were unique to each tir 255 allele variant. The overall average number of prophages did not differ between tir 255 allele strains; however, there were different types between the strains. Genomic and mobile element variation linked to the tir 255 polymorphism may account for the increased frequency of the T allele isolates in human disease.

Clarification of relationship between single-nucleotide polymorphism panels of Shiga toxin-producing Escherichia coli O157:H7/H- strains

Eighty strains of enterohemorrhagic Escherichia coli O157:H7/H- were analyzed by three single-nucleotide polymorphism (SNP) panels using whole-genome sequencing data. The partial concordance of SNP types among the different SNP panels was observed on minimum spanning trees reconstructed with SNP data. As for lineage I/II strains, some of the clade 7 strains belonged to one unique SNP type as determined by three panels, suggesting that clade 7 should be divided into at least two genotypes, namely, the unique type and the rest. In addition, clade 8 contained two unique genotypes, which was consistent with the previous prediction. Similarly, for lineage II, clade 12 should be divided into three genotype strains. In contrast, many strains of several clades belonging to lineage I were clustered into the same node on each minimum spanning tree upon testing with the three SNP panels. Previous studies reported that lineage I diverged more recently than lineages I/II and II. Such low diversity in lineage I in this study may have arisen because this lineage has not accumulated SNPs because of its relatively recent divergence. Based on the concordance observed in this study, some of the previously published O157 genotype distribution data were successfully interpreted to clarify the clade distribution, which was well supported by previous literature.

Added Value of Genomic Surveillance of Virulence Factors in Shiga Toxin-Producing Escherichia coli in New South Wales, Australia

The disease caused by Shiga toxin-producing Escherichia coli (STEC) remains a significant public health challenge globally, but the incidence of human STEC infections in Australia remains relatively low. This study examined the virulence characteristics and diversity of STEC isolates in the state of New South Wales between December 2017 and May 2020. Utilisation of both whole and core genome multi-locus sequence typing (MLST) allowed for the inference of genomic diversity and detection of isolates that were likely to be epidemiologically linked. The most common STEC serotype and stx subtype detected in this study were O157:H7 and stx_1a, respectively. A genomic scan of other virulence factors present in STEC suggested interplay between iron uptake system and virulence factors that mediate either iron release or countermeasures against host defence that could result in a reduction of stx_1a expression. This reduced expression of the dominant stx genotype could contribute to the reduced incidence of STEC-related illness in Australia. Genomic surveillance of STEC becomes an important part of public health response and ongoing interrogation of virulence factors in STEC offers additional insights for the public health risk assessment.

Whole-genome sequence analysis of Shiga toxin-producing Escherichia coli O157 strains isolated from wild deer and boar in Japan

The prevalence of Shiga toxin-producing Escherichia coli O157 (STEC O157) strains in wild deer and boar in Japan was investigated. STEC O157 strains were isolated from 1.9% (9/474) of the wild deer and 0.7% (3/426) of the wild boar examined. Pulsed-field gel electrophoresis (PFGE) analysis classified the wild deer and boar strains into four and three PFGE patterns, respectively. The PFGE pattern of one wild boar strain was similar to that of a cattle strain that had been isolated from a farm in the same area the wild boar was caught, suggesting that a STEC O157 strain may have been transmitted between wild boar and cattle. Clade analysis indicated that, although most of the strains were classified in clade 12, two strains were classified in clade 7. Whole-genome sequence (WGS) analysis indicated that all the strains carried mdfA, a drug resistance gene for macrolide antibiotics, and also pathogenicity-related genes similar to those in the Sakai strain. In conclusion, our study emphasized the importance of food hygiene in processing meat from Japanese wild animals for human consumption.

Phylogeny, Prevalence, and Shiga Toxin (Stx) Production of Clinical Escherichia coli O157 Clade 2 Strains Isolated in Shimane Prefecture, Japan

This study investigated the genetic and pathogenic variation of the subgroups of clade 2 strains of Shiga toxin (Stx)-producing Escherichia coli (STEC) O157. A total of 111 strains of STEC O157 isolated in Shimane prefecture, Japan, were classified in clade 2 (n = 39), clade 3 (n = 16), clade 4/5 (n = 3), clade 7 (n = 14), clade 8 (n = 17), and clade 12 (n = 22) by single-nucleotide polymorphism analysis and lineage-specific polymorphism assay-6. These results showed a distinct difference from our previous study in which clade 3 strains were the most prevalent strains in three other prefectures in Japan, indicating that the clade distribution of O157 strains was different in different geographic areas in Japan. Phylogenetic analysis using insertion sequence (IS) 629 distribution data showed that clade 2 strains formed two clusters, designated 2a and 2b. Stx2 production by cluster 2b strains was significantly higher than by cluster 2a strains (P < 0.01). In addition, population genetic analysis of the clade 2 strains showed significant linkage disequilibrium in the IS629 distribution of the strains in clusters 2a and 2b (P < 0.05). The Φ_PT values calculated using the IS629 distribution data indicated that strains in clusters 2a and 2b were genetically different (P < 0.001). Cluster 2b strains are a highly pathogenic phylogenetic group and their geographic spread may be a serious public health concern.

Pathogenomes of Atypical Non-shigatoxigenic Escherichia coli NSF/SF O157:H7/NM: Comprehensive Phylogenomic Analysis Using Closed Genomes

The toxigenic conversion of Escherichia coli strains by Shiga toxin-converting (Stx) bacteriophages were prominent and recurring events in the stepwise evolution of enterohemorrhagic E. coli (EHEC) O157:H7 from an enteropathogenic (EPEC) O55:H7 ancestor. Atypical, attenuated isolates have been described for both non-sorbitol fermenting (NSF) O157:H7 and SF O157:NM serotypes, which are distinguished by the absence of Stx, the characteristic virulence hallmark of Stx-producing E. coli (STEC). Such atypical isolates either never acquired Stx-phages or may have secondarily lost stx during the course of infection, isolation, or routine subculture; the latter are commonly referred to as LST (Lost Shiga Toxin)-isolates. In this study we analyzed the genomes of 15 NSF O157:H7 and SF O157:NM strains from North America, Europe, and Asia that are characterized by the absence of stx, the virulence hallmark of STEC. The individual genomic basis of the Stx (-) phenotype has remained largely undetermined as the majority of STEC genomes in public genome repositories were generated using short read technology and are in draft stage, posing a major obstacle for the high-resolution whole genome sequence typing (WGST). The application of LRT (long-read technology) sequencing provided us with closed genomes, which proved critical to put the atypical non-shigatoxigenic NSF O157:H7 and SF O157:NM strains into the phylogenomic context of the stepwise evolutionary model. Availability of closed chromosomes for representative Stx (-) NSF O157:H7 and SF O157:NM strains allowed to describe the genomic basis and individual evolutionary trajectories underlying the absence of Stx at high accuracy and resolution. The ability of LRT to recover and accurately assemble plasmids revealed a strong correlation between the strains' featured plasmid genotype and chromosomally inferred clade, which suggests the coevolution of the chromosome and accessory plasmids. The identified ancestral traits in the pSFO157 plasmid of NSF O157:H7 strain LSU-61 provided additional evidence for its intermediate status. Taken together, these observations highlight the utility of LRTs for advancing our understanding of EHEC O157:H7/NM pathogenome evolution. Insights into the genomic and phenotypic plasticity of STEC on a lineage- and genome-wide scale are foundational to improve and inform risk assessment, biosurveillance, and prevention strategies.

Molecular subtyping and clonal relatedness of human and cattle verotoxin-producing Escherichia coli O157:H7 isolates

Verotoxin-producing Escherichia coli O157:H7 is the dominant serotype isolated from patients with hemolytic-uremic syndrome (HUS) and, Argentina has the highest rate of HUS in the world. However, not all O157:H7 isolates have the same ability to infect and cause disease in humans. It has been postulated that O157:H7 strains integrate subpopulations related to the origin and virulence. In order to study the population structure and genetic diversity of VTEC O157:H7 from Argentina, a combination of molecular subtyping methods such as multiple loci VNTR analysis (MLVA), single nucleotide polymorphisms (SNP) and phylogroups assignment were used. According to MLVA, high genetic diversity was found among strains isolated from cattle, humans and food. On the other hand, 92% of the isolates presented the allele tir 255 T > A T and 95% were assigned to phylogroup E. We did not find a significant association between the isolates origin and the allele T presence (P > 0,05) postulated as significantly overrepresented in human isolates. Our results show that human and cattle VTEC O157:H7 isolates from Argentina are a homogeneous group and, although it presents high genetic diversity in relation to their MLVA and virulence profiles, it is not possible to distinguish divergent populations. The presence in all the strains of a high number of T3SS effectors genes and the no association of genetic subtypes with strain source, is an alert about the potential risk in public health that VTEC O157:H7 cattle strains possess and, at less, a partial explication about the high incidence of HUS in Argentina.

Phage Predation Shapes the Population Structure of Shiga-Toxigenic Escherichia coli O157:H7 in the UK: An Evolutionary Perspective

Bacterial–host interactions are non-linear and actually threefold, involving significant selection through predatory lytic bacteriophages in the host environment. In studies of human and animal gut microbiome bacteria, it is important to consider phage in all host–pathogen interactions. We use an important zoonotic pathogen, Shiga toxigenic Escherichia coli (STEC) O157:H7, to investigate this. Our study provides evidence that phage resistance profiles are well maintained at the sub-lineage level with variation in profiles within sub-lineages uncommon. This indicates that phage resistance heterogeneity happened early on in the STEC O157:H7 natural history and that occasional “wobbles” do not often outcompete the stable lineage unless combined with a competitive advantage. We discuss an example of this in the acquisition of stx2a that, while an important virulence factor, also conveys increased phage cross-resistance. We also discuss the role of phage resistance in co-occurrence of the three stable lineages worldwide and whether differing phage resistance is maintaining diversity.

Emergence and divergence of major lineages of Shiga-toxin-producing Escherichia coli in Australia

Shiga-toxin-producing Escherichia coli (STEC) infection is an important global cause of foodborne disease. To date however, genomics-based studies of STEC have been predominately focused upon STEC collected in the Northern Hemisphere. Here, we demonstrate the population structure of 485 STEC isolates in Australia, and show that several clonal groups (CGs) common to Australia were infrequently detected in a representative selection of contemporary STEC genomes from around the globe. Further, phylogenetic analysis demonstrated that lineage II of the global O157:H7 STEC was most prevalent in Australia, and was characterized by a frameshift mutation in flgF, resulting in the H-non-motile phenotype. Strong concordance between in silico and phenotypic serotyping was observed, along with concordance between in silico and conventional detection of stx genes. These data represent the most comprehensive STEC analysis from the Southern Hemisphere, and provide a framework for future national genomics-based surveillance of STEC in Australia.

Longitudinal Study of Shiga Toxin-Producing Escherichia coli and Campylobacter jejuni on Finnish Dairy Farms and in Raw Milk

Shiga toxin-producing Escherichia coli (STEC) and Campylobacter jejuni are notable health hazards associated with the consumption of raw milk. These bacteria may colonize the intestines of asymptomatic cattle and enter bulk tank milk via fecal contamination during milking. We studied the frequency of STEC O157:H7 and C. jejuni contamination in tank milk (n = 785) and the in-line milk filters of milking machines (n = 631) versus the frequency of isolation from cattle feces (n = 257) on three Finnish dairy farms for 1 year. Despite simultaneous isolation of STEC O157:H7 (17%) or C. jejuni (53%) from cattle, these bacteria were rarely isolated from milk filters (2% or <1%, respectively) and milk (0%). As revealed by phylogenomics, one STEC O157:H7 strain at a time was detected on each farm and persisted for ≤12 months despite rigorous hygienic measures. C. jejuni strains of a generalist sequence type (ST-883 and ST-1080) persisted in the herds for ≥11 months, and several other C. jejuni types were detected sporadically. The stx gene carried by STEC was detected more frequently from milk filters (37%) than from milk (7%), suggesting that milk filters are more suitable sampling targets for monitoring than milk. A questionnaire of on-farm practices suggested lower stx contamination of milk when major cleansing in the barn, culling, or pasturing of dairy cows was applied, while a higher average outdoor temperature was associated with higher stx contamination. Because pathogen contamination occurred despite good hygiene and because pathogen detection from milk and milk filters proved challenging, we recommend heat treatment for raw milk before consumption.IMPORTANCE The increased popularity of raw milk consumption has created demand for relaxing legislation, despite the risk of contamination by pathogenic bacteria, notably STEC and C. jejuni However, the epidemiology of these milk-borne pathogens on the herd level is still poorly understood, and data are lacking on the frequency of milk contamination on farms with cattle shedding these bacteria in their feces. This study suggests (i) that STEC contamination in milk can be reduced, but not prevented, by on-farm hygienic measures while fecal shedding is observable, (ii) that milk filters are more suitable sampling targets for monitoring than milk although pathogen detection from both sample matrices may be challenging, and (iii) that STEC and C. jejuni genotypes may persist in cattle herds for several months. The results can be utilized in developing and targeting pathogen monitoring and risk management on the farm level and contributed to the revision of Finnish legislation in 2017.

Molecular Typing of Escherichia coli O157 Isolates from Romanian Human Cases

Verocytotoxin-producing Escherichia coli (VTEC) of serogroup O157 are among the most important causes of severe cases of foodborne disease and outbreaks worldwide. As little is known about the characteristic of these strains in Romania, we aimed to provide reference information on the virulence gene content, phylogenetic background, and genetic diversity of 7 autochthonous O157 strains collected during 2016 and 2017 from epidemiologically non-related cases. These strains were typed by a combination of phenotypic and molecular methods routinely used by the national reference laboratory. Additionally, 4 of them were subjected to whole-genome sequencing (WGS), and public web-based tools were used to extract information on virulence gene profiles, multilocus sequence types (MLST), and single nucleotide polymorphism (SNP)-based phylogenetic relatedness. Molecular typing provided evidence of the circulation of a polyclonal population while distinguishing a cluster of non-sorbitol-fermenting, glucuronidase-negative, phylogenetic group E, MLST 1804 strains, representing lineage II and clade 7, which harbored vtx2c, eae-gamma, and ehxA genes. A good correlation between the routine typing methods and WGS data was observed. However, SNP-based genotyping provided a higher resolution in depicting the relationships between the O157:H7 strains than that provided by Pulse-field gel electrophoresis. This study should be a catalyst for improved laboratory-based surveillance of autochthonous VTEC.

Genotypic Features of Clinical and Bovine Escherichia coli O157 Strains Isolated in Countries with Different Associated-Disease Incidences

There is great geographical variation in the frequency of Escherichia coli O157 infections that correlates with important differences in the bovine reservoir of each country. Our group carried out a broad molecular characterization of human and bovine E. coli O157 strains circulating in Argentina using different methodologies. Our data allows us to conclude that in Argentina, a high homogeneity is observed in both cattle and human strains, with almost exclusive circulation of strains belonging to the hypervirulent clade 8 described by Manning. The aim of this review was to compare the genetic background of E. coli O157 strains isolated in countries that have conducted similar studies, to try to correlate specific O157 genotypes with the incidence and severity of E. coli O157 associated diseases. The characteristics of the strains that cause disease in humans reflect the predominant genotypes in cattle in each of the countries analyzed. The main features clearly linked to high incidence or severity of E. coli O157 infections are lineage-specific polymorphism assay-6 lineage I/II, clade 8 strains and probably, clade 6 strains, the stx_2a/stx_2c genotype, the presence of q933 and q21 simultaneously, and putative virulence factor EC_3286. In countries with an absence of these features in O157 strains, the overall incidence of O157 disease is low. Argentina, where these characteristics are detected in most strains, shows the highest incidence of hemolytic uremic syndrome (HUS) worldwide.

Enterohemorrhagic Escherichia coli O157 subclade 8b strains in Chiba Prefecture, Japan, produced larger amounts of Shiga toxin 2 than strains in subclade 8a and other clades

Enterohemorrhagic Escherichia coli O157 (O157) strains can be classified into clades (one of several phylogenetic groups) by single nucleotide polymorphisms (SNPs): these are clade 1, clade 2, clade 3, descendant and ancestral clades 4/5, clade 6, clade 7, clade 8, clade 9, and clade 12. Some recent studies showed that some O157 strains in clade 8 produced a larger amount of Shiga toxin (Stx) 2 than other strains. In this study, 1121 epidemiologically unlinked strains of O157 isolated in Chiba Prefecture, Japan were classified into clades during 1996–2014. Clade 8 strains were further classified into subclade 8a (67 strains) and subclade 8b (48 strains) using SNP analysis. In the absence of mitomycin C (MMC), subclade 8a strains in this study produced significantly greater amounts of Stx2 than subclade 8b strains. However, in the presence of MMC, the levels of Stx2 production in subclade 8b strains were significantly greater than subclade 8a strains. On the other hand, a recent study reported that the Stx2 production level in O157 strains was determined mainly by the subtypes of Stx2a phage (ϕStx2_α, β, γ, δ, ε, and ζ). Using O157 strains in this study, the Stx2a phages were classified into these subtypes. In this study, all strains of subclades 8a and 8b carried ϕStx2a_γ and ϕStx2a_δ, respectively. Some strains in clade 6 also carried ϕStx2a_δ. In the presence of MMC, subclade 8b strains produced significantly greater amounts of Stx2 than clade 6 strains carrying ϕStx2_δ. In this study, we propose that Stx2 production in subclade 8b strains in the presence of MMC might be enhanced due to genetic factors other than ϕStx2_δ.

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.

Applying phylogenomics to understand the emergence of Shiga-toxin-producing Escherichia coli O157:H7 strains causing severe human disease in the UK

Shiga-toxin-producing Escherichia coli (STEC) O157:H7 is a recently emerged zoonotic pathogen with considerable morbidity. Since the emergence of this serotype in the 1980s, research has focussed on unravelling the evolutionary events from the E. coli O55:H7 ancestor to the contemporaneous globally dispersed strains observed today. In this study, the genomes of over 1000 isolates from both human clinical cases and cattle, spanning the history of STEC O157:H7 in the UK, were sequenced. Phylogenetic analysis revealed the ancestry, key acquisition events and global context of the strains. Dated phylogenies estimated the time to evolution of the most recent common ancestor of the current circulating global clone to be 175 years ago. This event was followed by rapid diversification. We show the acquisition of specific virulence determinates has occurred relatively recently and coincides with its recent detection in the human population. We used clinical outcome data from 493 cases of STEC O157:H7 to assess the relative risk of severe disease including haemolytic uraemic syndrome from each of the defined clades in the population and show the dramatic effect Shiga toxin repertoire has on virulence. We describe two strain replacement events that have occurred in the cattle population in the UK over the last 30 years, one resulting in a highly virulent strain that has accounted for the majority of clinical cases in the UK over the last decade. There is a need to understand the selection pressures maintaining Shiga-toxin-encoding bacteriophages in the ruminant reservoir and the study affirms the requirement for close surveillance of this pathogen in both ruminant and human populations.

Putative classification of clades of enterohemorrhagic Escherichia coli O157 using an IS-printing system

Enterohemorrhagic Escherichia coli O157 (O157) strains can be classified in clades by single nucleotide polymorphisms (SNPs), but this analysis requires significant laboratory effort. As the distribution of insertion sequence (IS) 629 insertions has been reported to be biased among different clades, O157 isolates can be putatively classified in clades by comparison with an IS629 distribution database. A database of the IS629 distribution in O157 strains isolated in Chiba Prefecture and their classification in clades was determined by SNP analysis and IS-printing, an easy and quick analytical tool for IS629 in the O157 genome. The IS629 distribution in O157 strains isolated in Fukuoka and Yamagata Prefectures was determined by IS-printing. These strains were putatively classified in clades by Relative Likelihood calculations that compared the IS-printing data and the IS629 distribution database. Concordance Ratios were calculated, which compared the number of strains putatively classified in a clade by Relative Likelihood to the number of strains classified in that clade by SNP analysis. For the Fukuoka and Yamagata strains, the Concordance Ratios for clades 3, 6 and 8 were 97-100%, for clade 7 about 88%, and for clades 2 and 12 over 90%. In conclusion, O157 clade 2, 3, 6, 7, 8 and 12 strains could be putatively classified by IS-printing.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study demonstrated that enterohemorrhagic E. coli O157 (O157) strains could be putatively classified in clades using an IS-printing system. IS-printing was previously developed as a relatively quick and easy tool for analysis of insertion sequence 629 in the O157 genome. Since most local government public health institutes in Japan carry out IS-printing for early detection of O157 outbreaks, these data should be useful for putative classification of O157 strains in each area.

Comparative analysis of Mycobacterium tuberculosis Beijing strains isolated in three remote areas of Japan

A quantitative and qualitative comparison was carried out of Mycobacterium tuberculosis Beijing strains isolated in three remote areas of Japan. A total of 452 strains from Chiba Prefecture, 75 from Yamagata Prefecture, and 315 from Kobe City were analyzed for 24 loci by variable number of tandem repeats typing (24(Beijing)-VNTR). All strains were classified in six Beijing subgroups (B(SUB)), B1 to B5 and T, based on a minimum spanning tree reconstructed using data of a standard set of 15 VNTR loci. No significant difference was found in the distribution of strains in the B(SUB) in the three areas, with one exception due to a B5 outbreak in Yamagata, indicating no significant quantitative difference in the B(SUB) in the three areas (P<0.01, Chi-square test). In addition, when strains in each B(SUB) isolated in the three areas were mixed and standardized index of association (I(A)(s)) and variance (Φ(PT)) values were calculated, no significant qualitative difference in the B(SUB) in the three areas was found. These results suggested that the B(SUB) diverged prior to the introduction of M. tuberculosis Beijing strains into Japan. Differences in the distribution of strains in each B(SUB) between Japan and continental Asian countries suggested there had been genetic drift in the continental Asian countries in which B4 had been dominant.

An evolutionary analysis of nitric oxide reductase gene norV in enterohemorrhagic Escherichia coli O157

A novel virulence gene, norV, that encodes nitric oxide (NO) reductase, was examined to investigate the emergence of enterohemorrhagic Escherichia coli (EHEC) O157 subgroup C clusters 2 and 3 from subgroup C cluster 1. Deletion of norV occurred at a point between cluster 1 and cluster 2 just after or at the same time that an stx2 bacteriophage, which retains Shiga toxin 2 gene, was inserted into wrbA, which encodes a novel multimeric flavodoxin-like protein, in EHEC O157. Sensitivity of NO to anaerobic growth was correlated with the deletion of norV in all EHEC O157 individuals tested. The C467A mutation of fimH, which encodes minor component of type 1 fimbriae, occurred within cluster 1, not as a transition from cluster 1 to cluster 2, indicating that there is a cluster 1 minority branch that leads to cluster 2. These data refine the evolutionary history of an emerging EHEC O157.

Public Health Investigation of Two Outbreaks of Shiga Toxin-Producing Escherichia coli O157 Associated with Consumption of Watercress

An increase in the number of cases of Shiga toxin-producing Escherichia coli (STEC) O157 phage type 2 (PT2) in England in September 2013 was epidemiologically linked to watercress consumption. Whole-genome sequencing (WGS) identified a phylogenetically related cluster of 22 cases (outbreak 1). The isolates comprising this cluster were not closely related to any other United Kingdom strain in the Public Health England WGS database, suggesting a possible imported source. A second outbreak of STEC O157 PT2 (outbreak 2) was identified epidemiologically following the detection of outbreak 1. Isolates associated with outbreak 2 were phylogenetically distinct from those in outbreak 1. Epidemiologically unrelated isolates on the same branch as the outbreak 2 cluster included those from human cases in England with domestically acquired infection and United Kingdom domestic cattle. Environmental sampling using PCR resulted in the isolation of STEC O157 PT2 from irrigation water at one implicated watercress farm, and WGS showed this isolate belonged to the same phylogenetic cluster as outbreak 2 isolates. Cattle were in close proximity to the watercress bed and were potentially the source of the second outbreak. Transfer of STEC from the field to the watercress bed may have occurred through wildlife entering the watercress farm or via runoff water. During this complex outbreak investigation, epidemiological studies, comprehensive testing of environmental samples, and the use of novel molecular methods proved invaluable in demonstrating that two simultaneous outbreaks of STEC O157 PT2 were both linked to the consumption of watercress but were associated with different sources of contamination.

Analysis of the population genetics of clades of enterohaemorrhagic Escherichia coli O157:H7/H- isolated in three areas in Japan

AIMS

The genetic differences of enterohaemorrhagic Escherichia coli O157 (O157) strains isolated from humans in three widely-separated areas in Japan were analysed to provide information on possible geographic aspects of O157 pathogenicity.

METHODS AND RESULTS

Epidemiologically unlinked O157 strains were isolated in Chiba (300 strains), Fukuoka (260 strains) and Yamagata (81 strains) prefectures. These strains were classified in clades by single nucleotide polymorphism in seven loci and lineage-specific polymorphism assay-6, and differences between the strains in each clade were compared by population genetic analyses using the IS-printing system. Analysis of the clades from the three areas showed linkage disequilibrium of the strains in each clade. Comparison of the genetic differences of strains from the three areas in each clade, from calculated ΦPT values, indicated that the strains in each clade were the same population in all three areas, except possibly the clade 12 strains.

CONCLUSIONS

Population genetics analyses confirmed that the distribution of O157 strains in the clades isolated in three areas in Japan were similar and stable.

SIGNIFICANCE AND IMPACT OF THE STUDY

The pathogenicity of O157 strains infecting humans was comparable due to the similar, stable geographic distribution of O157 clades.

Phylogenetic and population genetic analysis of Salmonella enterica subsp. enterica serovar Infantis strains isolated in Japan using whole genome sequence data

Salmonella enterica subsp. enterica serovar Infantis has been reported to be carried asymptomatically in chickens and isolated from some human patients with diarrhea. The aim of this study was to investigate the phylogeny of S. Infantis strains isolated in Japan from chicken meat, chicken egg shells, environmental samples from a grading and packing center for chicken eggs (GP center), diarrhea patients, and asymptomatic carriers based on whole genome sequence data. The S. Infantis strains were in five clusters in a phylogenetic tree reconstructed by the maximum likelihood method. The five clusters were confirmed by neighbor-net and Bayesian cluster analyses. Two of the five clusters formed a group containing all of the strains isolated from chicken meat and some of the strains isolated from diarrhea patients and asymptomatic carriers. The median-joining network reconstructed in this study showed that strains in one of these two clusters diverged from one node with similar relatively short branches, suggesting clonal dissemination of these strains. The other three clusters formed a group containing all of the strains isolated from chicken egg shells and the GP center, and the remaining strains from diarrhea patients and asymptomatic carriers. Interestingly, strains isolated from patients did not cluster in only one group, indicating that none of the S. Infantis strains in this study had significantly higher human pathogenicity. The population genetic analyses in this study showed the separation of the five clusters into two groups was concordant with the sources where the strains in the clusters were isolated. These results suggested that evolutionary groups with higher hierarchy than the clusters identified in this study may be present, although such groups could not be determined by phylogenetic, neighbor-net, and Bayesian analyses in this study. Determination of higher level S. Infantis evolutionary groups should be investigated using other types of genetic markers.

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.

Linkage disequilibrium of the IS629 insertion among different clades of enterohemorrhagic Escherichia coli O157:H7/H-strains

The distribution of insertion sequence (IS) 629 was investigated among enterohemorrhagic Escherichia coli O157:H7/H-(O157) strains in different clades. Minimum spanning tree analysis showed that most strains in each clade clustered in a separate branch, indicating biased distribution of the IS629 insertion in different clades. The standardized index of association of the IS629 distribution data showed linkage disequilibrium in the strains in every clade, indicating that IS629 distribution data could be used for population genetic analysis. For this population genetic analysis, the Φ(PT) value, an analogue of F(ST), was calculated and indicated that clade 7 strains could be split into two clades based on their lineages. The degree of pairwise linkage disequilibrium was quite different among these two groups. The clade 7 split was in agreement with the model of O157 paraphyletic evolution and a new designation of the lineage II clades was proposed. The prevalence of strains with an IS629 insertion in certain loci was significantly different in different clades. Some of these significant differences were in loci in strains in branches of clades that were close in the O157 phylogenetic model, suggesting that IS629 insertion/deletion was not directly related to the divergence of O157 clades.