Insight into Shiga toxin genes encoded by Escherichia coli O157 from whole genome sequencing

Genomic analysis of an outbreak of Shiga toxin-producing Escherichia coli O183:H18 in the United Kingdom, 2023

In June 2023, UKHSA surveillance systems detected an outbreak of severe gastrointestinal symptoms caused by a rare serotype of Shiga toxin-producing Escherichia coli, STEC O183:H18. There were 26 cases aged 6 months to 74 years (42 % cases were aged 0-9 years), distributed across the UK with onset dates range between 22 May 2023 and 4 July 2023. The epidemiological and food chain investigations were inconclusive, although meat products made from beef mince were implicated as a potential vehicle. The outbreak strain belonged to sequence type (ST) 657 and harboured a Shiga toxin (stx) subtype stx2a located on a prophage that was unique in the UKHSA stx-encoding bacteriophage database. Plasmid encoded, putative virulence genes subA, ehxA, saa, iha, lpfA and iss were detected, however, the established STEC virulence genes involved in attachment to the gut mucosa (eae and aggR) were absent. The acquisition of stx across the global population structure of ST657 appeared to correspond with the presence of subA, ehxA, saa, iha, lpfA and iss. During the outbreak investigation, we used long read sequencing to characterise the plasmid and prophage content of this atypical STEC, to look for evidence to explain its recent emergence. Although we were unable to determine source and transmission route of the outbreak strain, the genomic analysis revealed potential clues as to how novel strains for STEC evolve. With the implementation of PCR capable of detecting all STEC, and genome sequencing for typing and virulence profiling, we have the tools to enable us to monitor the changing landscape of STEC. Improvements in the standardised collection of epidemiological data and trace-back strategies within the food industry, will ensure we have a surveillance system capable of alerting us to emerging threats to public health.

A Comprehensive Review on Shiga Toxin Subtypes and Their Niche-Related Distribution Characteristics in Shiga-Toxin-Producing E. coli and Other Bacterial Hosts

Shiga toxin (Stx), the main virulence factor of Shiga-toxin-producing E. coli (STEC), was first discovered in Shigella dysenteriae strains. While several other bacterial species have since been reported to produce Stx, STEC poses the most significant risk to human health due to its widespread prevalence across various animal hosts that have close contact with human populations. Based on its biochemical and molecular characteristics, Shiga toxin can be grouped into two types, Stx1 and Stx2, among which a variety of variants and subtypes have been identified in various bacteria and host species. Interestingly, the different Stx subtypes appear to vary in their host distribution characteristics and in the severity of diseases that they are associated with. As such, this review provides a comprehensive overview on the bacterial species that have been recorded to possess stx genes to date, with a specific focus on the various Stx subtype variants discovered in STEC, their prevalence in certain host species, and their disease-related characteristics. This review provides a better understanding of the Stx subtypes and highlights the need for rapid and accurate approaches to toxin subtyping for the proper evaluation of the health risks associated with Shiga-toxin-related bacterial food contamination and human infections.

Pathogenomes and virulence profiles of representative big six non-O157 serogroup Shiga toxin-producing Escherichia coli

Shiga toxin (Stx)-producing Escherichia coli (STEC) of non-O157:H7 serotypes are responsible for global and widespread human food-borne disease. Among these serogroups, O26, O45, O103, O111, O121, and O145 account for the majority of clinical infections and are colloquially referred to as the "Big Six." The "Big Six" strain panel we sequenced and analyzed in this study are reference type cultures comprised of six strains representing each of the non-O157 STEC serogroups curated and distributed by the American Type Culture Collection (ATCC) as a resource to the research community under panel number ATCC MP-9. The application of long- and short-read hybrid sequencing yielded closed chromosomes and a total of 14 plasmids of diverse functions. Through high-resolution comparative phylogenomics, we cataloged the shared and strain-specific virulence and resistance gene content and established the close relationship of serogroup O26 and O103 strains featuring flagellar H-type 11. Virulence phenotyping revealed statistically significant differences in the Stx-production capabilities that we found to be correlated to the strain's individual stx-status. Among the carried Stx1a, Stx2a, and Stx2d phages, the Stx2a phage is by far the most responsive upon RecA-mediated phage mobilization, and in consequence, stx2a + isolates produced the highest-level of toxin in this panel. The availability of high-quality closed genomes for this "Big Six" reference set, including carried plasmids, along with the recorded genomic virulence profiles and Stx-production phenotypes will provide a valuable foundation to further explore the plasticity in evolutionary trajectories in these emerging non-O157 STEC lineages, which are major culprits of human food-borne disease.

Clinical and public health implications of increasing notifications of LEE-negative Shiga toxin-producing Escherichia coli in England, 2014-2022

Introduction. Shiga toxin-producing Escherichia coli (STEC) belong to a diverse group of gastrointestinal pathogens. The pathogenic potential of STEC is enhanced by the presence of the pathogenicity island called the Locus of Enterocyte Effacement (LEE), including the intimin encoding gene eae.Gap statement. STEC serotypes O128:H2 (Clonal Complex [CC]25), O91:H14 (CC33), and O146:H21 (CC442) are consistently in the top five STEC serotypes isolated from patients reporting gastrointestinal symptoms in England. However, they are eae/LEE-negative and perceived to be a low risk to public health, and we know little about their microbiology and epidemiology.Aim. We analysed clinical outcomes and genome sequencing data linked to patients infected with LEE-negative STEC belonging to CC25 (O128:H2, O21:H2), CC33 (O91:H14) and, and CC442 (O146:H21, O174:H21) in England to assess the risk to public health.Results. There was an almost ten-fold increase between 2014 and 2022 in the detection of all STEC belonging to CC25, CC33 and CC442 (2014 n=38, 2022 n=336), and a total of 1417 cases. There was a higher proportion of female cases (55-70 %) and more adults than children, with patients aged between 20-40 and >70 most at risk across the different serotypes. Symptoms were consistent across the three dominant serotypes O91:H14 (CC33), O146:H21 (CC442) and O128:H2 (CC25) (diarrhoea >75 %; bloody diarrhoea 25-32 %; abdominal pain 64-72 %; nausea 37-45 %; vomiting 10-24 %; and fever 27-30 %). Phylogenetic analyses revealed multiple events of acquisition and loss of different stx-encoding prophage. Additional putative virulence genes were identified including iha, agn43 and subA.Conclusions. Continued monitoring and surveillance of LEE-negative STEC infections is essential due to the increasing burden of infectious intestinal disease, and the risk that highly pathogenic strains may emerge following acquisition of the Shiga toxin subtypes associated with the most severe clinical outcomes.

Phylogenetic relationship and virulence composition of Escherichia coli O26:H11 cattle and human strain collections in Scotland; 2002-2020

O26 is the commonest non-O157 Shiga toxin (stx)-producing Escherichia coli serogroup reported in human infections worldwide. Ruminants, particularly cattle, are the primary reservoir source for human infection. In this study, we compared the whole genomes and virulence profiles of O26:H11 strains (n = 99) isolated from Scottish cattle with strains from human infections (n = 96) held by the Scottish Escherichia coli O157/STEC Reference Laboratory, isolated between 2002 and 2020. Bovine strains were from two national cross-sectional cattle surveys conducted between 2002-2004 and 2014-2015. A maximum likelihood phylogeny was constructed from a core-genome alignment with the O26:H11 strain 11368 reference genome. Genomes were screened against a panel of 2,710 virulence genes using the Virulence Finder Database. All stx-positive bovine O26:H11 strains belonged to the ST21 lineage and were grouped into three main clades. Bovine and human source strains were interspersed, and the stx subtype was relatively clade-specific. Highly pathogenic stx2a-only ST21 strains were identified in two herds sampled in the second cattle survey and in human clinical infections from 2010 onwards. The closest pairwise distance was 9 single-nucleotide polymorphisms (SNPs) between Scottish bovine and human strains and 69 SNPs between the two cattle surveys. Bovine O26:H11 was compared to public EnteroBase ST29 complex genomes and found to have the greatest commonality with O26:H11 strains from the rest of the UK, followed by France, Italy, and Belgium. Virulence profiles of stx-positive bovine and human strains were similar but more conserved for the stx2a subtype. O26:H11 stx-negative ST29 (n = 17) and ST396 strains (n = 5) were isolated from 19 cattle herds; all were eae-positive, and 10 of these herds yielded strains positive for ehxA, espK, and Z2098, gene markers suggestive of enterohaemorrhagic potential. There was a significant association (p < 0.001) between nucleotide sequence percent identity and stx status for the bacteriophage insertion site genes yecE for stx2 and yehV for stx1. Acquired antimicrobial resistance genes were identified in silico in 12.1% of bovine and 17.7% of human O26:H11 strains, with sul2, tet, aph(3″), and aph(6″) being most common. This study describes the diversity among Scottish bovine O26:H11 strains and investigates their relationship to human STEC infections.

Prevalence of Shiga Toxin-Producing Escherichia coli O157 in Wild Scottish Deer with High Human Pathogenic Potential

Shiga toxin-producing E. coli (STEC) infections associated with wildlife are increasing globally, highlighting many 'spillover' species as important reservoirs for these zoonotic pathogens. A human outbreak of STEC serogroup O157 in 2015 in Scotland, associated with the consumption of venison meat products, highlighted several knowledge gaps, including the prevalence of STEC O157 in Scottish wild deer and the potential risk to humans from wild deer isolates. In this study, we undertook a nationwide survey of wild deer in Scotland and determined that the prevalence of STEC O157 in wild deer is low 0.28% (95% confidence interval = 0.06-0.80). Despite the low prevalence of STEC O157 in Scottish wild deer, identified isolates were present in deer faeces at high levels (>104 colony forming units/g faeces) and had high human pathogenic potential based on whole genome sequencing and virulence gene profiling. A retrospective epidemiological investigation also identified one wild deer isolate from this study as a possible source of a Scottish human outbreak in 2017. These results emphasise the importance of food hygiene practices during the processing of wild deer carcasses for human consumption.

Analysis of Escherichia coli O157 strains in cattle and humans between Scotland and England & Wales: implications for human health

For the last two decades, the human infection frequency of Escherichia coli O157 (O157) in Scotland has been 2.5-fold higher than in England and Wales. Results from national cattle surveys conducted in Scotland and England and Wales in 2014/2015 were combined with data on reported human clinical cases from the same time frame to determine if strain differences in national populations of O157 in cattle could be associated with higher human infection rates in Scotland. Shiga toxin subtype (Stx) and phage type (PT) were examined within and between host (cattle vs human) and nation (Scotland vs England and Wales). For a subset of the strains, whole genome sequencing (WGS) provided further insights into geographical and host association. All three major O157 lineages (I, II, I/II) and most sub-lineages (Ia, Ib, Ic, IIa, IIb, IIc) were represented in cattle and humans in both nations. While the relative contribution of different reservoir hosts to human infection is unknown, WGS analysis indicated that the majority of O157 diversity in human cases was captured by isolates from cattle. Despite comparable cattle O157 prevalence between nations, strain types were localized. PT21/28 (sub-lineage Ic, Stx2a+) was significantly more prevalent in Scottish cattle [odds ratio (OR) 8.7 (2.3-33.7; P<0.001] and humans [OR 2.2 (1.5-3.2); P<0.001]. In England and Wales, cattle had a significantly higher association with sub-lineage IIa strains [PT54, Stx2c; OR 5.6 (1.27-33.3); P=0.011] while humans were significantly more closely associated with sub-lineage IIb [PT8, Stx1 and Stx2c; OR 29 (4.9-1161); P<0.001]. Therefore, cattle farms in Scotland were more likely to harbour Stx2a+O157 strains compared to farms in E and W (P<0.001). There was evidence of limited cattle strain migration between nations and clinical isolates from one nation were more similar to cattle isolates from the same nation, with sub-lineage Ic (mainly PT21/28) exhibiting clear national association and evidence of local transmission in Scotland. While we propose the higher rate of O157 clinical cases in Scotland, compared to England and Wales, is a consequence of the nationally higher level of Stx2a+O157 strains in Scottish cattle, we discuss the multiple additional factors that may also contribute to the different infection rates between these nations.

First molecular characterization of Escherichia coli O157:H7 isolates from clinical samples in Paraguay using whole-genome sequencing

Escherichia coli O157:H7 is a foodborne pathogen implicated in numerous outbreaks worldwide that has the ability to cause extra-intestinal complications in humans. The Enteropathogens Division of the Central Public Health Laboratory (CPHL) in Paraguay is working to improve the genomic characterization of Shiga toxin-producing E. coli (STEC) to enhance laboratory-based surveillance and investigation of foodborne disease outbreaks. Whole genome sequencing (WGS) is proposed worldwide to be used in the routine laboratory as a high-resolution tool that allows to have all the results in a single workflow. This study aimed to carry out for the first time, the genomic characterization by WGS of nine STEC O157:H7 strains isolated from human samples in Paraguay. We were able to identify virulence and resistance mechanisms, MLST subtype, and even establish the phylogenetic relationships between isolates. Furthermore, we detected the presence of strains belonging to hypervirulent clade 8 in most of the isolates studied.

Outbreak of STEC O157:H7 linked to a milk pasteurisation failure at a dairy farm in England, 2019

In November 2019, an outbreak of Shiga toxin-producing Escherichia coli O157:H7 was detected in South Yorkshire, England. Initial investigations established consumption of milk from a local dairy as a common exposure. A sample of pasteurised milk tested the next day failed the phosphatase test, indicating contamination of the pasteurised milk by unpasteurised (raw) milk. The dairy owner agreed to immediately cease production and initiate a recall. Inspection of the pasteuriser revealed a damaged seal on the flow divert valve. Ultimately, there were 21 confirmed cases linked to the outbreak, of which 11 (52%) were female, and 12/21 (57%) were either <15 or >65 years of age. Twelve (57%) patients were treated in hospital, and three cases developed haemolytic uraemic syndrome. Although the outbreak strain was not detected in the milk samples, it was detected in faecal samples from the cattle on the farm. Outbreaks of gastrointestinal disease caused by milk pasteurisation failures are rare in the UK. However, such outbreaks are a major public health concern as, unlike unpasteurised milk, pasteurised milk is marketed as ‘safe to drink’ and sold to a larger, and more dispersed, population. The rapid, co-ordinated multi-agency investigation initiated in response to this outbreak undoubtedly prevented further cases.

Characterisation of atypical Shiga toxin gene sequences and description of Stx2j, a new subtype

Shiga toxin (Stx) is the definitive virulence factor of Shiga toxin-producing Escherichia coli (STEC). Stx variants are currently organised into a taxonomic system of three Stx1 (a,c,d) and seven Stx2 (a,b,c,d,e,f,g) subtypes. In this study, seven STEC isolates from food and clinical samples possessing stx2 sequences that do not fit current Shiga toxin taxonomy were identified. Genome assemblies of the STEC strains was created from Oxford Nanopore and Illumina sequence data. The presence of atypical stx2 sequences were confirmed by Sanger sequencing, as were Stx2 expression and cytotoxicity. A strain of O157:H7 was found to possess stx1a and a truncated stx2a, which were originally misidentified as an atypical stx2. Two strains possessed unreported variants of Stx2a (O8:H28) and Stx2b (O146:H21). In four of the strains we found three Stx-subtypes that are not included in the current taxonomy. Stx2h (O170:H18) was identified in a Canadian sprout isolate; this subtype has only previously been reported in STEC from Tibetan Marmots. Stx2o (O85:H1) was identified in a clinical isolate. Finally, Stx2j (O158:H23 and O33:H14) was found in lettuce and clinical isolates. The results of this study expands the number of known Stx subtypes, the range of STEC serotypes, and isolation sources in which they may be found. The presence of the Stx2j and Stx2o in clinical isolates of STEC indicates that strains carrying these variants are potential human pathogens. Highlights Atypical Shiga toxin (stx) genes in Escherichia coli were sequenced. Two new variants of stx2a and stx2b are described. Two strains carried subtypes Stx2h and Stx2o, which have only one previous report. Two strains carried a previously undescribed subtype, Stx2j.

Inconsistent PCR detection of Shiga toxin-producing Escherichia coli: Insights from whole genome sequence analyses

Shiga toxin-producing Escherichia coli (STEC) have been linked to food-borne disease outbreaks. As PCR is routinely used to screen foods for STEC, it is important that factors leading to inconsistent detection of STEC by PCR are understood. This study used whole genome sequencing (WGS) to investigate causes of inconsistent PCR detection of stx1, stx2, and serogroup-specific genes. Fifty strains isolated from Alberta feedlot cattle from three different studies were selected with inconsistent or consistent detection of stx and serogroup by PCR. All isolates were initially classified as STEC by PCR. Sequencing was performed using Illumina MiSeq® with sample library by Nextera XT. Virtual PCRs were performed using Geneious and bacteriophage content was determined using PHASTER. Sequencing coverage ranged from 47 to 102x, averaging 74x, with sequences deposited in the NCBI database. Eleven strains were confirmed by WGS as STEC having complete stxA and stxB subunits. However, truncated stx fragments occurred in twenty-two other isolates, some having multiple stx fragments in the genome. Isolates with complete stx by WGS had consistent stx1 and stx2 detection by PCR, although one also having a stx2 fragment had inconsistent stx2 PCR. For all STEC and 18/39 non-STEC, serogroups determined by PCR agreed with those determined by WGS. An additional three WGS serotypes were inconclusive and two isolates were Citrobacter spp. Results demonstrate that stx fragments associated with stx-carrying bacteriophages in the E. coli genome may contribute to inconsistent detection of stx1 and stx2 by PCR. Fourteen isolates had integrated stx bacteriophage but lacked complete or fragmentary stx possibly due to partial bacteriophage excision after sub-cultivation or other unclear mechanisms. The majority of STEC isolates (7/11) did not have identifiable bacteriophage DNA in the contig(s) where stx was located, likely increasing the stability of stx in the bacterial genome and its detection by PCR.

Epidemiological investigations identified an outbreak of Shiga toxin-producing Escherichia coli serotype O26:H11 associated with pre-packed sandwiches

In October 2019, public health surveillance systems in Scotland identified an increase in the number of reported infections of Shiga toxin-producing Escherichia coli (STEC) O26:H11 involving bloody diarrhoea. Ultimately, across the United Kingdom (UK) 32 cases of STEC O26:H11 stx1a were identified, with the median age of 27 years and 64% were male; six cases were hospitalised. Among food exposures there was an association with consuming pre-packed sandwiches purchased at outlets belonging to a national food chain franchise (food outlet A) [odds ratio (OR) = 183.89, P < 0.001]. The common ingredient identified as a component of the majority of the sandwiches sold at food outlet A was a mixed salad of Apollo and Iceberg lettuce and spinach leaves. Microbiological testing of food and environmental samples were negative for STEC O26:H11, although STEC O36:H19 was isolated from a mixed salad sample taken from premises owned by food outlet A. Contamination of fresh produce is often due to a transient event and detection of the aetiological agent in food that has a short-shelf life is challenging. Robust, statistically significant epidemiological analysis should be sufficient evidence to direct timely and targeted on-farm investigations. A shift in focus from testing the microbiological quality of the produce to investigating the processes and practices through the supply chain and sampling the farm environment is recommended.

Outbreak of Shiga toxin-producing Escherichia coli O157 linked with consumption of a fast-food product containing imported cucumbers, United Kingdom, August 2020

BACKGROUND

In August 2020, an outbreak of Shiga toxin-producing Escherichia coli (STEC) O157:H7 occurred in the United Kingdom. Whole genome sequencing revealed that these cases formed a genetically distinct cluster.

METHODS

Hypotheses generated from case interviews were tested in analytical studies, and results informed environmental sampling and food chain analysis. A case-case study used non-outbreak 'comparison' STEC cases; a case-control study used a market research panel to recruit controls.

RESULTS

A total of 36 cases were identified; all cases reported symptom onset between August 3 and August 16, 2020. The majority of cases (83%) resided in the Midlands region of England and in Wales. A high proportion of cases reported eating out, with one fast-food restaurant chain mentioned by 64% (n = 23) of cases. Both the case-case study (adjusted odds ratio (aOR) 31.8, 95% confidence interval (CI) 1.6-624.9) and the case-control study (aOR 9.19, 95% CI 1.0-82.8) revealed statistically significant results, showing that the consumption of a specific fast-food product was independently associated with infection.

CONCLUSIONS

Consumption of a specific fast-food product was a likely cause of this outbreak. The only ingredient specific to the product was cucumbers. The supply of cucumbers was immediately halted, and no further cases have been identified.

Phylogenetic structure of Shiga toxin-producing Escherichia coli O157:H7 from sub-lineage to SNPs

Sequence similarity of pathogen genomes can infer the relatedness between isolates as the fewer genetic differences identified between pairs of isolates, the less time since divergence from a common ancestor. Clustering based on hierarchical single linkage clustering of pairwise SNP distances has been employed to detect and investigate outbreaks. Here, we evaluated the evidence-base for the interpretation of phylogenetic clusters of Shiga toxin-producing Escherichia coli (STEC) O157:H7. Whole genome sequences of 1193 isolates of STEC O157:H7 submitted to Public Health England between July 2015 and December 2016 were mapped to the Sakai reference strain. Hierarchical single linkage clustering was performed on the pairwise SNP difference between all isolates at descending distance thresholds. Cases with known epidemiological links fell within 5-SNP single linkage clusters. Five-SNP single linkage community clusters where an epidemiological link was not identified were more likely to be temporally and/or geographically related than sporadic cases. Ten-SNP single linkage clusters occurred infrequently and were challenging to investigate as cases were few, and temporally and/or geographically dispersed. A single linkage cluster threshold of 5-SNPs has utility for the detection of outbreaks linked to both persistent and point sources. Deeper phylogenetic analysis revealed that the distinction between domestic UK and imported isolates could be inferred at the sub-lineage level. Cases associated with domestically acquired infection that fall within clusters that are predominantly travel associated are likely to be caused by contaminated imported food.

MLST based serotype prediction for the accurate identification of non typhoidal Salmonella serovars

Traditional serotyping based on the phenotypic variation of O- and H-antigen remains as the gold-standard for the identification and classification of Salmonella isolates for last 70 years. Although this classification is a globally recognized nomenclature, huge diversity of Salmonella serotypes have made the serovar identification to be very complex. Seven gene multilocus sequence typing (MLST) on the other hand can provide serovar prediction as well as the evolutionary origin between the serovars. In this study non typhoidal Salmonella (NTS) strains (n = 45) isolated from clinical samples (blood, faeces and pus) were identified by traditional phenotypic serotyping and biochemical testing. All the tested Salmonella isolates were designated as serovar Typhimurium based on phenotyping. However, by MLST 60% (27/45) of the isolates were S. Typhimurium, 35.5% (16/45) were S. Agona (ST13), 2.2% (1/45) were S. Kentucky (ST198) and 2.2% (1/45) were S. Saintpaul (ST27). MLST analysis assigned S. Typhimurium isolates as ST36 (18/127), ST19 (7/27) and ST313 (2/27). Mismatches in serovar designation between MLST database and phenotypic serotyping can be due to the misinterpretation of phenotypic serotyping as the antigenic structures of S. Typhimurium, S. Agona differs by a surface antigen. MLST based phylogeny of study isolates showed clustering according to sequence types. Concordance between MLST based sequence type and phenotypic serotype is important to provide insights into genetic population structure of Salmonella.

Impact of whole genome sequencing on the investigation of food-borne outbreaks of Shiga toxin-producing Escherichia coli serogroup O157:H7, England, 2013 to 2017

We aim to provide insight and guidance on the utility of whole genome sequencing (WGS) data for investigating food-borne outbreaks of Shiga toxin-producing Escherichia coli (STEC) O157:H7 in England between 2013 and 2017. Analysis of WGS data delivered an unprecedented level of strain discrimination when compared with multilocus variable number tandem repeat analysis. The robustness of the WGS method ensured confidence in the microbiological identification of linked cases, even when epidemiological links were obscured. There was evidence that phylogeny derived from WGS data can be used to trace the geographical origin of an isolate. Further analysis of the phylogenetic data provided insight on the evolutionary context of emerging pathogenic strains. Publically available WGS data linked to the clinical, epidemiological and environmental context of the sequenced strain has improved trace back investigations during outbreaks. Expanding the use of WGS-based typing analysis globally will ensure the rapid implementation of interventions to protect public health, inform risk assessment and facilitate the management of national and international food-borne outbreaks of STEC O157:H7.

Outbreak of Shiga toxin-producing Escherichia coli O157:H7 linked to raw drinking milk resolved by rapid application of advanced pathogen characterisation methods, England, August to October 2017

An outbreak of Shiga toxin-producing Escherichia coli (STEC) O157:H7 occurred on the Isle of Wight between August and October 2017. Of the seven cases linked to the outbreak, five were identified through the statutory notification process and two were identified through national surveillance of whole genome sequencing data. Enhanced surveillance questionnaires established a common link to a farm, and link to the likely food vehicle, raw drinking milk (RDM). Microbiological investigations, including PCR, identified the presence of STEC O157:H7 in samples of RDM. Analysis of core genome single nucleotide polymorphism (SNP) data of STEC O157:H7 from human stool specimens, animal faecal samples and RDM demonstrated a one SNP difference between isolates, and therefore close genetic relatedness. Control measures that were put in place included suspension of sales and recall of RDM, as well as restrictions on public access to parts of the farm. Successful integration of traditional epidemiological surveillance and advanced laboratory methods for the detection and characterisation of STEC O157:H7 from human, animal and environmental samples enabled prompt identification of the outbreak vehicle and provided evidence to support the outbreak control team’s decision-making, leading to implementation of effective control measures in a timely manner.

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.

Evaluation of chromogenic selective agar (CHROMagar STEC) for the direct detection of Shiga toxin-producing Escherichia coli from faecal specimens

Shiga toxin-producing Escherichia coli (STEC) are zoonotic pathogens that cause symptoms of severe gastrointestinal disease, including haemolytic uraemic syndrome (HUS), in humans. Currently in England, STEC serotypes other than O157:H7 are not cultured at the local hospital laboratories. The aim of this study was to evaluate the utility of CHROMagar STEC for the direct detection of STEC from faecal specimens in a diagnostic setting, compared to the current reference laboratory method using PCR targeting the Shiga-toxin gene (stx) to test multiple colonies cultured on MacConkey agar. Of the 292 consecutive faecal specimens submitted to the Gastrointestinal Bacterial Reference Unit that tested positive for stx by PCR, STEC could not be cultured on MacConkey agar or CHROMagar STEC from 87/292 (29.8 %). Of the 205 that were cultured, 106 (51.7 %) were detected on both MacConkey agar and CHROMagar STEC and 99 (48.3 %) were detected on MacConkey agar only. All 106 (100 %) isolates that grew on CHROMagar STEC had the ter gene cassette, known to be associated with resistance to tellurite, compared to 13/99 (13.1 %) that were not detected on CHROMagar STEC. CHROMagar STEC supported the growth of 36/40 (90 %) isolates harbouring stx2a or stx2d, the subtypes most frequently associated with progression to HUS. Of the 92 isolates harbouring eae, an important STEC virulence marker, 77 (83.7 %) grew on CHROMagar STEC. CHROMagar STEC is a useful selective media for the rapid, near-patient detection of STEC that have the potential to cause HUS.

Contribution and Interaction of Shiga Toxin Genes to Escherichia coli O157:H7 Virulence

Escherichia coli O157:H7 is the predominant cause of diarrhea-associated hemolytic uremic syndrome (HUS) worldwide. Its cardinal virulence traits are Shiga toxins, which are encoded by stx genes, the most common of which are stx1a, stx2a, and stx2c. The toxins these genes encode differ in their in vitro and experimental phenotypes, but the human population-level impact of these differences is poorly understood. Using Shiga toxin-encoding bacteriophage insertion typing and real-time polymerase chain reaction, we genotyped isolates from 936 E. coli O157:H7 cases and verified HUS status via chart review. We compared the HUS risk between isolates with stx2a and those with stx2a and another gene and estimated additive interaction of the stx genes. Adjusted for age and symptoms, the HUS incidence of E. coli O157:H7 containing stx2a alone was 4.4% greater (95% confidence interval (CI) −0.3%, 9.1%) than when it occurred with stx1a. When stx1a and stx2a occur together, the risk of HUS was 27.1% lower (95% CI −87.8%, −2.3%) than would be expected if interaction were not present. At the population level, temporal or geographic shifts toward these genotypes should be monitored, and stx genotype may be an important consideration in clinically predicting HUS among E. coli O157:H7 cases.

Eleven High-Quality Reference Genome Sequences and 360 Draft Assemblies of Shiga Toxin-Producing Escherichia coli Isolates from Human, Food, Animal, and Environmental Sources in Canada

We report high-quality closed reference genomes for 1 bovine strain and 10 human Shiga toxin (Stx)-producing Escherichia coli (STEC) strains from serogroups O26, O45, O91, O103, O104, O111, O113, O121, O145, and O157. We also report draft assemblies, with standardized metadata, for 360 STEC strains isolated from watersheds, animals, farms, food, and human infections.

We report high-quality closed reference genomes for 1 bovine strain and 10 human Shiga toxin (Stx)-producing Escherichia coli (STEC) strains from serogroups O26, O45, O91, O103, O104, O111, O113, O121, O145, and O157. We also report draft assemblies, with standardized metadata, for 360 STEC strains isolated from watersheds, animals, farms, food, and human infections.

Shiga toxin sub-type 2a increases the efficiency of Escherichia coli O157 transmission between animals and restricts epithelial regeneration in bovine enteroids

Specific Escherichia coli isolates lysogenised with prophages that express Shiga toxin (Stx) can be a threat to human health, with cattle being an important natural reservoir. In many countries the most severe pathology is associated with enterohaemorrhagic E. coli (EHEC) serogroups that express Stx subtype 2a. In the United Kingdom, phage type (PT) 21/28 O157 strains have emerged as the predominant cause of life-threatening EHEC infections and this phage type commonly encodes both Stx2a and Stx2c toxin types. PT21/28 is also epidemiologically linked to super-shedding (>10³ cfu/g of faeces) which is significant for inter-animal transmission and human infection as demonstrated using modelling studies. We demonstrate that Stx2a is the main toxin produced by stx2a⁺/stx2c⁺ PT21/28 strains induced with mitomycin C and this is associated with more rapid induction of gene expression from the Stx2a-encoding prophage compared to that from the Stx2c-encoding prophage. Bacterial supernatants containing either Stx2a and/or Stx2c were demonstrated to restrict growth of bovine gastrointestinal organoids with no restriction when toxin production was not induced or prevented by mutation. Isogenic strains that differed in their capacity to produce Stx2a were selected for experimental oral colonisation of calves to assess the significance of Stx2a for both super-shedding and transmission between animals. Restoration of Stx2a expression in a PT21/28 background significantly increased animal-to-animal transmission and the number of sentinel animals that became super-shedders. We propose that while both Stx2a and Stx2c can restrict regeneration of the epithelium, it is the relatively rapid and higher levels of Stx2a induction, compared to Stx2c, that have contributed to the successful emergence of Stx2a+ E. coli isolates in cattle in the last 40 years. We propose a model in which Stx2a enhances E. coli O157 colonisation of in-contact animals by restricting regeneration and turnover of the colonised gastrointestinal epithelium.

Enterohaemorrhagic E. coli (EHEC) O157 strains are found in cattle where they are asymptomatic, while human exposure can lead to severe symptoms including bloody diarrhoea and kidney damage due to the activity of Shiga toxin (Stx). The most serious symptoms in humans are associated with isolates that encode Stx subtype 2a. The advantage of these toxins in the animal reservoir is still not clear, however there is experimental evidence implicating Stx with increased bacterial adherence, immune modulation and suppression of predatory protozoa. In this study, the hypothesis that Stx2a is important for super-shedding and calf-to-calf transmission was tested by comparing excretion and transmission dynamics of E. coli O157 strains with and without Stx2a. While Stx2a did not alter excretion levels when calfs were orally challenge, it enabled colonisation of more in contact ‘sentinel’ animals in our transmission model. We show that Stx2a is generally induced more rapidly than Stx2c, resulting in increased levels of Stx2a expression. Both Stx2a and Stx2c were able to restrict cellular proliferation of epithelial cells in cultured bovine enteroids. Taken together, we propose that rapid production of Stx2a and its role in establishing E. coli O157 colonisation in the bovine gastrointestinal tract facilitate effective transmission and have led to its expansion in the cattle E. coli O157 population.

Phylogenetic context of Shiga toxin-producing Escherichia coli serotype O26:H11 in England

The increasing use of PCR for the detection of gastrointestinal pathogens in hospital laboratories in England has improved the detection of Shiga toxin-producing Escherichia coli (STEC), and the diagnosis of haemolytic uraemic syndrome (HUS). We aimed to analyse the microbiological characteristics and phylogenetic relationships of STEC O26:H11, clonal complex (CC) 29, in England to inform surveillance, and to assess the threat to public health. There were 502 STEC belonging to CC29 isolated between 2014 and 2019, of which 416 were from individual cases. The majority of isolates belonged to one of three major sequence types (STs), ST16 (n=37), ST21 (n=350) and ST29 (n=24). ST16 and ST29 were mainly isolated from cases reporting recent travel abroad. Within ST21, there were three main clades associated with domestic acquisition. All three domestic clades had Shiga toxin subtype gene (stx) profiles associated with causing severe clinical outcomes including STEC-HUS, specifically either stx1a, stx2a or stx1a/stx2a. Isolates from the same patient, same household or same outbreak with an established source for the most part fell within 5-SNP single linkage clusters. There were 19 5-SNP community clusters, of which six were travel-associated and one was an outbreak of 16 cases caused by the consumption of contaminated salad leaves. Of the remaining 12 clusters, 9/12 were either temporally or geographically related or both. Exposure to foodborne STEC O26:H11 ST21 capable of causing severe clinical outcomes, including STEC-HUS, is an emerging risk to public health in England. The lack of comprehensive surveillance of this STEC serotype is a concern, and there is a need to expand the implementation of methods capable of detecting STEC in local hospital settings.

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.

Comparison of single-nucleotide variants identified by Illumina and Oxford Nanopore technologies in the context of a potential outbreak of Shiga toxin-producing Escherichia coli

BACKGROUND

We aimed to compare Illumina and Oxford Nanopore Technology sequencing data from the 2 isolates of Shiga toxin-producing Escherichia coli (STEC) O157:H7 to determine whether concordant single-nucleotide variants were identified and whether inference of relatedness was consistent with the 2 technologies.

RESULTS

For the Illumina workflow, the time from DNA extraction to availability of results was ∼40 hours, whereas with the ONT workflow serotyping and Shiga toxin subtyping variant identification were available within 7 hours. After optimization of the ONT variant filtering, on average 95% of the discrepant positions between the technologies were accounted for by methylated positions found in the described 5-methylcytosine motif sequences, CC(A/T)GG. Of the few discrepant variants (6 and 7 difference for the 2 isolates) identified by the 2 technologies, it is likely that both methodologies contain false calls.

CONCLUSIONS

Despite these discrepancies, Illumina and Oxford Nanopore Technology sequences from the same case were placed on the same phylogenetic location against a dense reference database of STEC O157:H7 genomes sequenced using the Illumina workflow. Robust single-nucleotide polymorphism typing using MinION-based variant calling is possible, and we provide evidence that the 2 technologies can be used interchangeably to type STEC O157:H7 in a public health setting.

SMRT sequencing reveals differential patterns of methylation in two O111:H- STEC isolates from a hemolytic uremic syndrome outbreak in Australia

In 1995 a severe haemolytic-uremic syndrome (HUS) outbreak in Adelaide occurred. A recent genomic analysis of Shiga toxigenic Escherichia coli (STEC) O111:H- strains 95JB1 and 95NR1 from this outbreak found that the more virulent isolate, 95NR1, harboured two additional copies of the Shiga toxin 2 (Stx2) genes encoded within prophage regions. The structure of the Stx2-converting prophages could not be fully resolved using short-read sequence data alone and it was not clear if there were other genomic differences between 95JB1 and 95NR1. In this study we have used Pacific Biosciences (PacBio) single molecule real-time (SMRT) sequencing to characterise the genome and methylome of 95JB1 and 95NR1. We completely resolved the structure of all prophages including two, tandemly inserted, Stx2-converting prophages in 95NR1 that were absent from 95JB1. Furthermore we defined all insertion sequences and found an additional IS1203 element in the chromosome of 95JB1. Our analysis of the methylome of 95NR1 and 95JB1 identified hemi-methylation of a novel motif (5′-CTGC^m6AG-3′) in more than 4000 sites in the 95NR1 genome. These sites were entirely unmethylated in the 95JB1 genome, and included at least 177 potential promoter regions that could contribute to regulatory differences between the strains. IS1203 mediated deactivation of a novel type IIG methyltransferase in 95JB1 is the likely cause of the observed differential patterns of methylation between 95NR1 and 95JB1. This study demonstrates the capability of PacBio SMRT sequencing to resolve complex prophage regions and reveal the genetic and epigenetic heterogeneity within a clonal population of bacteria.

Highly Pathogenic Clone of Shiga Toxin-Producing Escherichia coli O157:H7, England and Wales

We used whole-genome sequencing to investigate the evolutionary context of an emerging highly pathogenic strain of Shiga toxin–producing Escherichia coli (STEC) O157:H7 in England and Wales. A timed phylogeny of sublineage IIb revealed that the emerging clone evolved from a STEC O157:H7 stx-negative ancestor ≈10 years ago after acquisition of a bacteriophage encoding Shiga toxin (stx) 2a, which in turn had evolved from a stx2c progenitor ≈20 years ago. Infection with the stx2a clone was a significant risk factor for bloody diarrhea (OR 4.61, 95% CI 2.24–9.48; p<0.001), compared with infection with other strains within sublineage IIb. Clinical symptoms of cases infected with sublineage IIb stx2c and stx-negative clones were comparable, despite the loss of stx2c. Our analysis highlighted the highly dynamic nature of STEC O157:H7 Stx-encoding bacteriophages and revealed the evolutionary history of a highly pathogenic clone emerging within sublineage IIb, a sublineage not previously associated with severe clinical symptoms.

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.

Validation of Whole-Genome Sequencing for Identification and Characterization of Shiga Toxin-Producing Escherichia coli To Produce Standardized Data To Enable Data Sharing

Whole-genome sequencing (WGS) is rapidly becoming the method of choice for outbreak investigations and public health surveillance of microbial pathogens. The combination of improved cluster resolution and prediction of resistance and virulence phenotypes provided by a single tool is extremely advantageous. However, the data produced are complex, and standard bioinformatics pipelines are required to translate the output into easily interpreted epidemiologically relevant information for public health action. The main aim of this study was to validate the implementation of WGS at the Scottish Escherichia coli O157/STEC Reference Laboratory (SERL) using the Public Health England (PHE) bioinformatics pipeline to produce standardized data to enable interlaboratory comparison of results generated at two national reference laboratories. In addition, we evaluated the BioNumerics whole-genome multilocus sequence typing (wgMLST) and E. coli genotyping plug-in tools using the same data set. A panel of 150 well-characterized isolates of Shiga toxin-producing E. coli (STEC) that had been sequenced and analyzed at PHE using the PHE pipeline and database (SnapperDB) was assembled to provide identification and typing data, including serotype (O:H type), sequence type (ST), virulence genes (eae and Shiga toxin [stx] subtype), and a single-nucleotide polymorphism (SNP) address. To validate the implementation of sequencing at the SERL, DNA was reextracted from the isolates and sequenced and analyzed using the PHE pipeline, which had been installed at the SERL; the output was then compared with the PHE data. The results showed a very high correlation between the data, ranging from 93% to 100%, suggesting that the standardization of WGS between our reference laboratories is possible. We also found excellent correlation between the results obtained using the PHE pipeline and BioNumerics, except for the detection of stx_2a and stx_2c when these subtypes are both carried by strains.

An outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with contaminated salad leaves: epidemiological, genomic and food trace back investigations

In August 2015, Public Health England detected an outbreak of Shiga toxin-producing Escherichia coli (STEC) serotype O157:H7 caused by contaminated salad leaves in a mixed leaf prepacked salad product from a national retailer. The implicated leaves were cultivated at five different farms and the zoonotic source of the outbreak strain was not determined. In March 2016, additional isolates from new cases were identified that shared a recent common ancestor with the outbreak strain. A case-case study involving the cases identified in 2016 revealed that ovine exposures were associated with illness (n = 16; AOR 8·24; 95% CI 1·55-39·74). By mapping the recent movement of sheep and lambs across the United Kingdom, epidemiological links were established between the cases reporting ovine exposures. Given the close phylogenetic relationship between the outbreak strain and the isolates from cases with ovine exposures, it is plausible that ovine faeces may have contaminated the salad leaves via untreated irrigation water or run-off from fields nearby. Timely and targeted veterinary and environmental sampling should be considered during foodborne outbreaks of STEC, particularly where ready to eat vegetables and salads are implicated.

Shiga Toxin Subtypes of Non-O157 Escherichia coli Serogroups Isolated from Cattle Feces

Shiga toxin producing Escherichia coli (STEC) are important foodborne pathogens responsible for human illnesses. Cattle are a major reservoir that harbor the organism in the hindgut and shed in the feces. Shiga toxins (Stx) are the primary virulence factors associated with STEC illnesses. The two antigenically distinct Stx types, Stx1 and Stx2, encoded by stx1 and stx2 genes, share approximately 56% amino acid sequence identity. Genetic variants exist within Stx1 and Stx2 based on differences in amino acid composition and in cytotoxicity. The objective of our study was to identify the stx subtypes in strains of STEC serogroups, other than O157, isolated from cattle feces. Shiga toxin gene carrying E. coli strains (n = 192), spanning 27 serogroups originating from cattle (n = 170) and human (n = 22) sources, were utilized in the study. Shiga toxin genes were amplified by PCR, sequenced, and nucleotide sequences were translated into amino acid sequences using CLC main workbench software. Shiga toxin subtypes were identified based on the amino acid motifs that define each subtype. Shiga toxin genotypes were also identified at the nucleotide level by in silico restriction fragment length polymorphism (RFLP). Of the total 192 STEC strains, 93 (48.4%) were positive for stx1 only, 43 (22.4%) for stx2 only, and 56 (29.2%) for both stx1 and stx2. Among the 149 strains positive for stx1, 132 (88.6%) were stx1a and 17 (11.4%) were stx1c. Shiga toxin 1a was the most common subtype of stx1 among cattle (87.9%; 123/140) and human strains (100%; 9/9) of non-O157 serogroups. Of the total 99 strains positive for stx2, 79 were stx2a (79.8%), 11 (11.1%) were stx2c, 12 (12.1%) were stx2d. Of the 170 strains originating from cattle feces, 58 (34.1%) were stx2a subtype, 11 (6.5%) were stx2c subtype, and 11 were of subtype stx2d (6.5%). All but one of the human strains were positive for stx2a. Three strains of cattle origin were positive for both stx2a and stx2d. In conclusion, a number of non-O157 STEC serogroups harbored by cattle possess a wide variety of Shiga toxin subtypes, with stx1a and stx2a being the most predominant stx subtypes occurring individually or in combination. Cattle are a reservoir of a number of non-O157 STEC serogroups and information on the Shiga toxin subtypes is useful in assessing the potential risk as human pathogens.

Severe Outbreak of Sorbitol-Fermenting Escherichia coli O157 via Unpasteurized Milk and Farm Visits, Finland 2012

Shiga toxin-producing, sorbitol-fermenting Escherichia coli O157 (SF O157) has emerged as a cause of severe human illness. Despite frequent human findings, its transmission routes and reservoirs remain largely unknown. Foodborne transmission and reservoir in cattle have been suspected, but with limited supporting evidence. This study describes the outbreak of SF O157 that occurred in Finland in 2012. The outbreak originated from a recreational farm selling unpasteurized milk, as revealed by epidemiologic and microbiological investigations, and involved six hospitalized children and two asymptomatic adults with culture-confirmed infection. An identical strain of SF O157 was isolated from patients, cattle and the farm environment, and epidemiologic analysis suggested unpasteurized milk as the vehicle of transmission. This study reports the first milkborne outbreak of SF O157, provides supporting evidence of cattle as a reservoir and highlights the health risks related to the consumption of unpasteurized milk.

Comparative Evaluation of Genomic and Laboratory Approaches for Determination of Shiga Toxin Subtypes in Escherichia coli

The determination of Shiga toxin (ST) subtypes can be an important element in the risk characterization of foodborne ST-producing Escherichia coli (STEC) isolates for making risk management decisions. ST subtyping methods include PCR techniques based on electrophoretic or pyrosequencing analysis of amplicons and in silico techniques based on whole genome sequence analysis using algorithms that can be readily incorporated into bioinformatics analysis pipelines for characterization of isolates by their genetic composition. The choice of technique will depend on the performance characteristics of the method and an individual laboratory's access to specialized equipment or personnel. We developed two whole genome sequence-based ST subtyping tools: (i) an in silico PCR algorithm requiring genome assembly to replicate a reference PCR-based method developed by the Statens Serum Institut (SSI) and (ii) an assembly-independent routine in which raw sequencing results are mapped to a database of known ST subtype sequence variants (V-Typer). These tools were evaluated alongside the SSI reference PCR method and a recently described PCR-based pyrosequencing technique. The V-Typer method results corresponded closely with the reference method in the analysis of 67 STEC cultures obtained from a World Health Organization National Reference Laboratory. In contrast, the in silico PCR method failed to detect ST subtypes in several cases, a result which we attribute to assembly-induced errors typically encountered with repetitive gene sequences. The V-Typer can be readily integrated into bioinformatics protocols used in the identification and characterization of foodborne STEC isolates.

Subtyping Escherichia coli Virulence Genes Isolated from Feces of Beef Cattle and Clinical Cases in Alberta

Clinical outcomes of Shiga toxin (stx)-producing Escherichia coli infection are largely determined by virulence gene subtypes. This study used a polymerase chain reaction (PCR)-pyrosequencing assay to analyze single-nucleotide polymorphisms for subtyping three major virulence genes (stx₁, stx₂, eae) of pathogenic E. coli (O157, O26, O111, and O103) isolated from cattle over a 2-year interval (n = 465) and human clinical cases (n = 42) in western Canada. Most bovine isolates were PCR positive for at least one target virulence gene (367/465), whereas 100% of human isolates harbored eae in combination with at least one stx gene. Four Shiga toxin (1a, 2a, 2c, and 2e) and four eae (λ/γ1-eae, ɛ-eae, θ/γ2-eae, and β-eae) subtypes were identified in over 25 distinct virulence genotypes. Among cattle isolates, every serogroup, but O103, presented a dominant genotype (O157: stx_1a+stx_2a+λ/γ1-eae, O26: β-eae alone, and O111: stx_1a+θ/γ2-eae). Similar patterns were found in human isolates, although it was not possible to establish a clear genotypic association between the two sources. Many O157 and non-O157 cattle isolates lacked stx genes; the absence was greater in non-O157 (75/258) and O157:non-H7 (19/40) than in O157:H7 strains (1/164). In addition, there was a greater diversity of virulence genotypes of E. coli isolated from cattle than those of human diseases, which could be due to sample characteristics (e.g., source and clinical condition). However, the majority of cattle strains had virulence profiles identical to those of clinical cases. Consequently, determining the presence of certain stx (stx_1a and stx_2a) and eae (λ/γ1-eae) subtypes known to cause human disease would be a valuable tool for risk assessment and prediction of disease outcome along the farm-to-fork continuum.

Whole genome sequencing improved case ascertainment in an outbreak of Shiga toxin-producing Escherichia coli O157 associated with raw drinking milk

Five cases of STEC O157 phage type (PT) 21/28 reported consumption of raw cows' drinking milk (RDM) produced at a dairy farm in the South West of England. STEC O157 PT21/28 was isolated from faecal specimens from milking cows on the implicated farm. Whole genome sequencing (WGS) showed that human and cattle isolates were the same strain. Further analysis of WGS data confirmed that sequences of isolates from an additional four cases (who did not report consumption of RDM when first questioned) fell within the same five single nucleotide polymorphism cluster as the initial five cases epidemiologically linked to the consumption of RDM. These four additional cases identified by WGS were investigated further and were, ultimately, associated with the implicated farm. The RDM outbreak strain encoded stx2a, which is associated with increased pathogenicity and severity of symptoms. Further epidemiological analysis showed that 70% of isolates within a wider cluster containing the outbreak strain were from cases residing in, or linked to, the same geographical region of England. During this RDM outbreak, use of WGS improved case ascertainment and provided insights into the evolution of a highly pathogenic clade of STEC O157 PT21/28 stx2a associated with the South West of England.

Implementation of Whole Genome Sequencing (WGS) for Identification and Characterization of Shiga Toxin-Producing Escherichia coli (STEC) in the United States

Shiga toxin-producing Escherichia coli (STEC) is an important foodborne pathogen capable of causing severe disease in humans. Rapid and accurate identification and characterization techniques are essential during outbreak investigations. Current methods for characterization of STEC are expensive and time-consuming. With the advent of rapid and cheap whole genome sequencing (WGS) benchtop sequencers, the potential exists to replace traditional workflows with WGS. The aim of this study was to validate tools to do reference identification and characterization from WGS for STEC in a single workflow within an easy to use commercially available software platform. Publically available serotype, virulence, and antimicrobial resistance databases were downloaded from the Center for Genomic Epidemiology (CGE) (www.genomicepidemiology.org) and integrated into a genotyping plug-in with in silico PCR tools to confirm some of the virulence genes detected from WGS data. Additionally, down sampling experiments on the WGS sequence data were performed to determine a threshold for sequence coverage needed to accurately predict serotype and virulence genes using the established workflow. The serotype database was tested on a total of 228 genomes and correctly predicted from WGS for 96.1% of O serogroups and 96.5% of H serogroups identified by conventional testing techniques. A total of 59 genomes were evaluated to determine the threshold of coverage to detect the different WGS targets, 40 were evaluated for serotype and virulence gene detection and 19 for the stx gene subtypes. For serotype, 95% of the O and 100% of the H serogroups were detected at > 40x and ≥ 30x coverage, respectively. For virulence targets and stx gene subtypes, nearly all genes were detected at > 40x, though some targets were 100% detectable from genomes with coverage ≥20x. The resistance detection tool was 97% concordant with phenotypic testing results. With isolates sequenced to > 40x coverage, the different databases accurately predicted serotype, virulence, and resistance from WGS data, providing a fast and cheaper alternative to conventional typing techniques.

Phenotypic and Genotypic Characteristics of Shiga Toxin-Producing Escherichia coli Isolated from Surface Waters and Sediments in a Canadian Urban-Agricultural Landscape

A hydrophobic grid membrane filtration-Shiga toxin immunoblot method was used to examine the prevalence of Shiga toxin-producing Escherichia coli (STEC) in four watersheds located in the Lower Mainland of British Columbia, Canada, a region characterized by rapid urbanization and intensive agricultural activity. STEC were recovered from 21.6, 23.2, 19.5, and 9.2% of surface water samples collected monthly from five sites in each watershed over a period of 1 year. Overall prevalence was subject to seasonal variation however, ranging between 13.3% during fall months and 34.3% during winter months. STEC were also recovered from 23.8% of sediment samples collected in one randomly selected site. One hundred distinct STEC isolates distributed among 29 definitive and 4 ambiguous or indeterminate serotypes were recovered from water and sediments, including isolates from Canadian "priority" serogroups O157 (3), O26 (4), O103 (5), and O111 (7). Forty seven isolates were further characterized by analysis of whole genome sequences to detect Shiga toxin gene (stx 1 and stx 2), intimin gene (eaeA) allelic variants and acquired virulence factors. These analyses collectively showed that surface waters from the region support highly diverse STEC populations that include strains with virulence factors commonly associated with human pathotypes. The present work served to characterize the microbiological hazard implied by STEC to support future assessments of risks to public health arising from non-agricultural and agricultural uses of surface water resources in the region.

Whole Genome Sequencing for Public Health Surveillance of Shiga Toxin-Producing Escherichia coli Other than Serogroup O157

Shiga toxin-producing Escherichia coli (STEC) are considered to be a significant threat to public health due to the severity of gastrointestinal symptoms associated with human infection. In England STEC O157 is the most commonly detected STEC serogroup, however, the implementation of PCR at local hospital laboratories has resulted in an increase in the detection of non-O157 STEC. The aim of this study was to evaluate the use of whole genome sequencing (WGS) for routine public health surveillance of non-O157 STEC by comparing this approach to phenotypic serotyping and PCR for subtyping the stx-encoding genes. Of the 102 isolates where phenotypic and genotypic serotyping could be compared, 98 gave fully concordant results. The most common non-O157 STEC serogroups detected were O146 (22) and O26 (18). All but one of the 38 isolates that could not be phenotypically serotyped (designated O unidentifiable or O rough) were serotyped using the WGS data. Of the 73 isolates where a flagella type was available by traditional phenotypic typing, all results matched the H-type derived from the WGS data. Of the 140 sequenced non-O157 isolates, 52 (37.1%) harboured stx1 only, 42 (30.0%) had stx2 only, 46 (32.9%) carried stx1 and stx2. Of these, stx subtyping PCR results were available for 131 isolates and 121 of these had concordant results with the stx subtype derived from the WGS data. Of the 10 discordant results, non-specific primer binding during PCR amplification, due to the similarity of the stx2 subtype gene sequences was the most likely cause. The results of this study showed WGS provided a reliable and robust one-step process for characterization of STEC. Deriving the full serotype from WGS data in real time has enabled us to report a higher level of strain discrimination while stx subtyping provides data on the pathogenic potential of each isolate, enabling us to predict clinical outcome of each case and to monitor the emergence of hyper-virulent strains.

Subtilase contributes to the cytotoxicity of a Shiga toxin-producing Escherichia coli strain encoding three different toxins

Food-borne Shiga toxin-producing Escherichia coli (STEC) O113:H21 strain TS18/08, that has previously been isolated from mixed minced meat, harbors the Shiga toxin (Stx) encoding allele stx2a, the plasmid-located subtilase cytotoxin encoding allele subAB1 and the cytolethal distending toxin type V encoding gene cdt-V. In the current study, it could be shown that each of these toxin genes was transcribed with different transcription levels at different time points by RT real time PCR under laboratory batch conditions in LB-broth. The transcription maximum for cdt-V and subAB1 was observed after 3h while stx2a transcription was highest after 6h of incubation. During this time the mean relationship of the amount of stx2a:subAB1:cdt-V transcripts was 1:26:100. Furthermore, isogenic stx2a and cdt-V chromosomal deletion mutants were constructed to measure the contribution of SubAB1 to the overall cytotoxicity of this strain. In this context, a further copy of stx2 was detected in this strain and was also deleted. Comparing the cytotoxicity of supernatants of the resulting mutant strains TS18/08-3 (Δstx2-1Δstx2-2Δcdt-V) and TS18/08-4 (Δstx2-1Δstx2-2Δcdt-VΔsubAB1) on Vero cells demonstrated a contribution of SubAB1 to the overall cytotoxic effect while the 4-fold isogenic deletion mutant did not show any cytotoxic effect and that was comparable to the non-toxic laboratory E. coli strain C600. The cytotoxic effect could be restored by complementation with the recombinant low copy plasmid pWSK29 harboring subAB1 under the control of its own promoter. In addition, the cytotoxicity of wild type strain TS18/08 to Vero cells was in the same range as the EHEC O157:H7 strain EDL933. Therefore, food-borne STEC O113:H21 strain TS18/08 can be considered as a putative human pathogen.

Whole Genome Sequencing demonstrates that Geographic Variation of Escherichia coli O157 Genotypes Dominates Host Association

Genetic variation in an infectious disease pathogen can be driven by ecological niche dissimilarities arising from different host species and different geographical locations. Whole genome sequencing was used to compare E. coli O157 isolates from host reservoirs (cattle and sheep) from Scotland and to compare genetic variation of isolates (human, animal, environmental/food) obtained from Scotland, New Zealand, Netherlands, Canada and the USA. Nei’s genetic distance calculated from core genome single nucleotide polymorphisms (SNPs) demonstrated that the animal isolates were from the same population. Investigation of the Shiga toxin bacteriophage and their insertion sites (SBI typing) revealed that cattle and sheep isolates had statistically indistinguishable rarefaction profiles, diversity and genotypes. In contrast, isolates from different countries exhibited significant differences in Nei’s genetic distance and SBI typing. Hence, after successful international transmission, which has occurred on multiple occasions, local genetic variation occurs, resulting in a global patchwork of continental and trans-continental phylogeographic clades. These findings are important for three reasons: first, understanding transmission and evolution of infectious diseases associated with multiple host reservoirs and multi-geographic locations; second, highlighting the relevance of the sheep reservoir when considering farm based interventions; and third, improving our understanding of why human disease incidence varies across the world.

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.

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

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

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.