Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read 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.

Re-analysis of an outbreak of Shiga toxin-producing Escherichia coli O157:H7 associated with raw drinking milk using Nanopore sequencing

The aim of this study was to compare Illumina and Oxford Nanopore Technology (ONT) sequencing data to quantify genetic variation to assess within-outbreak strain relatedness and characterise microevolutionary events in the accessory genomes of a cluster of 23 genetically and epidemiologically linked isolates related to an outbreak of Shiga toxin-producing Escherichia coli O157:H7 caused by the consumption of raw drinking milk. There were seven discrepant variants called between the two technologies, five were false-negative or false-positive variants in the Illumina data and two were false-negative calls in ONT data. After masking horizontally acquired sequences such as prophages, analysis of both short and long-read sequences revealed the 20 isolates linked to the outbreak in 2017 had a maximum SNP distance of one SNP between each other, and a maximum of five SNPs when including three additional strains identified in 2019. Analysis of the ONT data revealed a 47 kbp deletion event in a terminal compound prophage within one sample relative to the remaining samples, and a 0.65 Mbp large chromosomal rearrangement (inversion), within one sample relative to the remaining samples. Furthermore, we detected two bacteriophages encoding the highly pathogenic Shiga toxin (Stx) subtype, Stx2a. One was typical of Stx2a-phage in this sub-lineage (Ic), the other was atypical and inserted into a site usually occupied by Stx2c-encoding phage. Finally, we observed an increase in the size of the pO157 IncFIB plasmid (1.6 kbp) in isolates from 2019 compared to those from 2017, due to the duplication of insertion elements within the plasmids from the more recently isolated strains. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the genome plasticity and virulence between strains of this zoonotic, foodborne pathogen.

Characterization of the flanking region of the Shiga toxin operon in Stx2a bacteriophages reveals a diversity of the NanS-p sialate O-acetylesterase gene

Shiga toxin-producing E. coli (STEC) are diarrheagenic strains that can cause bloody diarrhea and hemolytic-uremic syndrome. Their main virulence factor, the Shiga toxin (Stx), is encoded by phages integrated into the bacterial chromosome. Stx phages are widely diverse and carry many genes with limited or unknown function. As the toxin subtype Stx2a is associated with highly pathogenic strains, this study was mainly focused on the characterization of the stx flanking region of Stx2a phages. Of particular interest was a sialate O-acetylesterase (NanS-p), which has been described previously to be encoded downstream stx in some phage genomes and may confer a growth advantage for STEC. Complete DNA sequences of Stx2a phages and prophages were retrieved from the GenBank database, and the genomic regions from anti-terminator Q to holin S genes were bioinformatically analyzed. Predicted NanSp sequences from phages encoding other Stx subtypes were also studied. Additionally, expression of nanS-p was quantified by qPCR in strains selected from our laboratory collection. The analysis of Stx2a phage genomes showed that all carried the Q, stx2a, nanS-p and S genes, but with allele diversity and other sequence differences. In particular, sequence differences were detected in each of the three domains of NanS-p esterases encoded by Stx2a phages and other Stx phages; however, nanS-p was not identified in the Stx2e, Stx2f and Stx2g phages analyzed. The expression of nanS-p increased in most stx2a-positive strains under phage inducing conditions, as was previously shown for stx2a. As the present work showed diversity at the Q-S region among Stx phages, and particularly in the encoded NanS-p enzyme, future studies will be necessary to evaluate if NanS-p variants differ in their activity and to assess the impact of the absence of nanS-p in certain Stx phages.

Dynamic changes in Shiga toxin (Stx) 1 transducing phage throughout the evolution of O26:H11 Stx-producing Escherichia coli

Shiga toxin (Stx) is the key virulence factor of Stx-producing Escherichia coli (STEC). All known Stxs (Stx1 and Stx2) are encoded by bacteriophages (Stx phages). Although the genetic diversity of Stx phages has frequently been described, systematic analyses of Stx phages in a single STEC lineage are limited. In this study, focusing on the O26:H11 STEC sequence type 21 (ST21) lineage, where the stx1a gene is highly conserved, we analysed the Stx1a phages in 39 strains representative of the entire ST21 lineage and found a high level of variation in Stx1a phage genomes caused by various mechanisms, including replacement by a different Stx1a phage at the same or different locus. The evolutionary timescale of events changing Stx1a phages in ST21 was also determined. Furthermore, by using an Stx1 quantification system developed in this study, we found notable variations in the efficiency of Stx1 production upon prophage induction, which sharply contrasted with the conserved iron regulated Stx1 production. These variations were associated with the Stx1a phage alteration in some cases but not in other cases; thus, Stx1 production in this STEC lineage was determined by differences not only in Stx1 phages but also in host-encoded factors.

Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli

Among the nine clades of Shiga toxin (Stx)-producing Escherichia coli O157:H7, clade 8 is thought to be highly pathogenic, as it causes severe disease more often than other clades. Two subclades have been proposed, but there are conflicting reports on intersubclade differences in Stx2 levels, although Stx2 production is a risk factor for severe disease development. The global population structure of clade 8 has also yet to be fully elucidated. Here, we present genome analyses of a global clade 8 strain set (n=510), including 147 Japanese strains sequenced in this study. The complete genome sequences of 18 of the 147 strains were determined to perform detailed clade-wide genome analyses together with 17 publicly available closed genomes. Intraclade variations in Stx2 production level and disease severity were also re-evaluated within the phylogenetic context. Based on phylogenomic analysis, clade 8 was divided into four lineages corresponding to the previously proposed SNP genotypes (SGs): SG8_30, SG8_31A, SG8_31B and SG8_32. SG8_30 and the common ancestor of the other SGs were first separated, with SG8_31A and SG8_31B emerging from the latter and SG8_32 emerging from SG8_31B. Comparison of 35 closed genomes revealed the overall structure of chromosomes and pO157 virulence plasmids and the prophage contents to be well conserved. However, Stx2a phages exhibit notable genomic diversity, even though all are integrated into the argW locus, indicating that subtype changes in Stx2a phage occurred from the γ subtype to its variant (γ_v1) in SG8_31A and from γ to δ in SG8_31B and SG8_32 via replacement of parts or almost entire phage genomes, respectively. We further show that SG8_30 strains (all carrying γ Stx2a phages) produce significantly higher levels of Stx2 and cause severe disease more frequently than SG8_32 strains (all carrying δ Stx2a phages). Clear conclusions on SG8_31A and SG8_31B cannot be made due to the small number of strains available, but as SG8_31A (carrying γ_v1 Stx2a phages) contains strains that produce much more Stx2 than SG8_30 strains, attention should also be paid to this SG.

Genomic Diversity, Virulence Gene, and Prophage Arrays of Bovine and Human Shiga Toxigenic and Enteropathogenic Escherichia coli Strains Isolated in Hungary

Escherichia coli belonging to the enterohemorrhagic (EHEC), Shiga toxin-producing (STEC) and atypical enteropathogenic (aEPEC) pathotypes are significant foodborne zoonotic pathogens posing serious health risks, with healthy cattle as their main reservoir. A representative sampling of Hungarian cattle farms during 2017–2018 yielded a prevalence of 6.5 and 5.8% for STEC and aEPEC out of 309 samples. The draft genomes of twelve STEC (of them 9 EHEC) and four aEPEC of bovine origin were determined. For comparative purposes, we also included 3 EHEC and 2 aEPEC strains of human origin, as well four commensal isolates and one extraintestinal pathogenic E. coli (ExPEC) obtained from animals in a final set of 26 strains for a WGS-based analysis. Apart from key virulence genes, these isolates harbored several additional virulence genes with arrays characteristic for the site of isolation. The most frequent insertion site of Shiga toxin (stx) encoding prophages was yehV for the Stx1 prophage and wrbA and sbcB for Stx2. For O157:H7 strains, the locus of enterocyte effacement pathogenicity island was present at the selC site, with integration at pheV for other serotypes, and pheU in the case of O26:H11 strains. Several LEE-negative STEC and aEPEC as well as commensal isolates carried additional prophages, with an average of ten prophage regions per isolate. Comparative phylogenomic analysis showed no clear separation between bovine and human lineages among the isolates characterized in the current study. Similarities in virulence gene arrays and close phylogenetic relations of bovine and human isolates underline the zoonotic potential of bovine aEPEC and STEC and emphasize the need for frequent monitoring of these pathogens in livestock.

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

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

Genome structural variation in Escherichia coli O157:H7

The human zoonotic pathogen Escherichia coli O157:H7 is defined by its extensive prophage repertoire including those that encode Shiga toxin, the factor responsible for inducing life-threatening pathology in humans. As well as introducing genes that can contribute to the virulence of a strain, prophage can enable the generation of large-chromosomal rearrangements (LCRs) by homologous recombination. This work examines the types and frequencies of LCRs across the major lineages of the O157:H7 serotype. We demonstrate that LCRs are a major source of genomic variation across all lineages of E. coli O157:H7 and by using both optical mapping and Oxford Nanopore long-read sequencing prove that LCRs are generated in laboratory cultures started from a single colony and that these variants can be recovered from colonized cattle. LCRs are biased towards the terminus region of the genome and are bounded by specific prophages that share large regions of sequence homology associated with the recombinational activity. RNA transcriptional profiling and phenotyping of specific structural variants indicated that important virulence phenotypes such as Shiga-toxin production, type-3 secretion and motility can be affected by LCRs. In summary, E. coli O157:H7 has acquired multiple prophage regions over time that act to continually produce structural variants of the genome. These findings raise important questions about the significance of this prophage-mediated genome contingency to enhance adaptability between environments.

Genomic analysis of shiga toxin-containing Escherichia coli O157:H7 isolated from Argentinean cattle

Cattle are the main reservoir of Enterohemorrhagic Escherichia coli (EHEC), with O157:H7 the distinctive serotype. EHEC is the main causative agent of a severe systemic disease, Hemolytic Uremic Syndrome (HUS). Argentina has the highest pediatric HUS incidence worldwide with 12–14 cases per 100,000 children. Herein, we assessed the genomes of EHEC O157:H7 isolates recovered from cattle in the humid Pampas of Argentina. According to phylogenetic studies, EHEC O157 can be divided into clades. Clade 8 strains that were classified as hypervirulent. Most of the strains of this clade have a Shiga toxin stx2a-stx2c genotype. To better understand the molecular bases related to virulence, pathogenicity and evolution of EHEC O157:H7, we performed a comparative genomic analysis of these isolates through whole genome sequencing. The isolates classified as clade 8 (four strains) and clade 6 (four strains) contained 13 to 16 lambdoid prophages per genome, and the observed variability of prophages was analysed. An inter strain comparison show that while some prophages are highly related and can be grouped into families, other are unique. Prophages encoding for stx2a were highly diverse, while those encoding for stx2c were conserved. A cluster of genes exclusively found in clade 8 contained 13 genes that mostly encoded for DNA binding proteins. In the studied strains, polymorphisms in Q antiterminator, the Q-stx2A intergenic region and the O and P γ alleles of prophage replication proteins are associated with different levels of Stx2a production. As expected, all strains had the pO157 plasmid that was highly conserved, although one strain displayed a transposon interruption in the protease EspP gene. This genomic analysis may contribute to the understanding of the genetic basis of the hypervirulence of EHEC O157:H7 strains circulating in Argentine cattle. This work aligns with other studies of O157 strain variation in other populations that shows key differences in Stx2a-encoding prophages.

Variability in the Occupancy of Escherichia coli O157 Integration Sites by Shiga Toxin-Encoding Prophages

Escherichia coli O157:H7 strains often produce Shiga toxins encoded by genes on lambdoid bacteriophages that insert into multiple loci as prophages. O157 strains were classified into distinct clades that vary in virulence. Herein, we used PCR assays to examine Shiga toxin (Stx) prophage occupancy in yehV, argW, wrbA, and sbcB among 346 O157 strains representing nine clades. Overall, yehV was occupied in most strains (n = 334, 96.5%), followed by wrbA (n = 213, 61.6%), argW (n = 103, 29.8%), and sbcB (n = 93, 26.9%). Twelve occupancy profiles were identified that varied in frequency and differed across clades. Strains belonging to clade 8 were more likely to have occupied sbcB and argW sites compared to other clades (p < 0.0001), while clade 2 strains were more likely to have occupied wrbA sites (p < 0.0001). Clade 8 strains also had more than the expected number of occupied sites based on the presence of stx variants (p < 0.0001). Deletion of a 20 kb non-Stx prophage occupying yehV in a clade 8 strain resulted in an ~18-fold decrease in stx2 expression. These data highlight the complexity of Stx prophage integration and demonstrate that clade 8 strains, which were previously linked to hemolytic uremic syndrome, have unique Stx prophage occupancy profiles that can impact stx2 expression.

Bacteriophages of Shiga Toxin-Producing Escherichia coli and Their Contribution to Pathogenicity

Shiga toxins (Stx) of Shiga toxin-producing Escherichia coli (STEC) are generally encoded in the genome of lambdoid bacteriophages, which spend the most time of their life cycle integrated as prophages in specific sites of the bacterial chromosome. Upon spontaneous induction or induction by chemical or physical stimuli, the stx genes are co-transcribed together with the late phase genes of the prophages. After being assembled in the cytoplasm, and after host cell lysis, mature bacteriophage particles are released into the environment, together with Stx. As members of the group of lambdoid phages, Stx phages share many genetic features with the archetypical temperate phage Lambda, but are heterogeneous in their DNA sequences due to frequent recombination events. In addition to Stx phages, the genome of pathogenic STEC bacteria may contain numerous prophages, which are either cryptic or functional. These prophages may carry foreign genes, some of them related to virulence, besides those necessary for the phage life cycle. Since the production of one or more Stx is considered the major pathogenicity factor of STEC, we aim to highlight the new insights on the contribution of Stx phages and other STEC phages to pathogenicity.

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.

Analysis of a small outbreak of Shiga toxin-producing Escherichia coli O157:H7 using long-read sequencing

Compared to short-read sequencing data, long-read sequencing facilitates single contiguous de novo assemblies and characterization of the prophage region of the genome. Here, we describe our methodological approach to using Oxford Nanopore Technology (ONT) sequencing data to quantify genetic relatedness and to look for microevolutionary events in the core and accessory genomes to assess the within-outbreak variation of four genetically and epidemiologically linked isolates. Analysis of both Illumina and ONT sequencing data detected one SNP between the four sequences of the outbreak isolates. The variant calling procedure highlighted the importance of masking homologous sequences in the reference genome regardless of the sequencing technology used. Variant calling also highlighted the systemic errors in ONT base-calling and ambiguous mapping of Illumina reads that results in variations in the genetic distance when comparing one technology to the other. The prophage component of the outbreak strain was analysed, and nine of the 16 prophages showed some similarity to the prophage in the Sakai reference genome, including the stx2a-encoding phage. Prophage comparison between the outbreak isolates identified minor genome rearrangements in one of the isolates, including an inversion and a deletion event. The ability to characterize the accessory genome in this way is the first step to understanding the significance of these microevolutionary events and their impact on the evolutionary history, virulence and potentially the likely source and transmission of this zoonotic, foodborne pathogen.

A Shiga Toxin-Encoding Prophage Recombination Event Confounds the Phylogenetic Relationship Between Two Isolates of Escherichia coli O157:H7 From the Same Patient

We compared genomes from multiple isolations of Shiga toxin-producing Escherichia coli (STEC) O157:H7 from the same patient, in cases notified to Public Health England (PHE) between 2015 and 2019. There were 261 cases where multiple isolates were sequenced from the same patient comprising 589 isolates. Serial isolates from the same patient fell within five single nucleotide polymorphisms (SNPs) of each other for 260/261 (99.6%) of the cases, indicating that there was little evidence of within host variation. The investigation into the 13 SNP discrepancy between one isolate pair revealed the cause to be a recombination event within a stx2a-encoding prophage resulting in the insertion/deletion of a fragment of the genome. This 50 kbp prophage fragment was homologous to a prophage in the reference genome, and the short reads from the isolate that had the 50 kbp fragment, mapped unambiguously to this region. The discrepant variants in the isolate without the 50 kbp fragment were attributed to ambiguous mapping of the short reads from other prophage regions to the 50 kbp fragment in the reference genome. Identification of such recombination events in this dataset appeared to be rare, most likely because the majority of prophage regions in the Sakai reference genome are masked during the analysis. Identification of SNPs under neutral selection, and masking recombination events, is a requirement for phylogenetic analysis used for public health surveillance, and for the detection of point source outbreaks. However, assaying the accessory genome by combining the use of short and long read technologies for public health surveillance may provide insight into how recombination events impact on the evolutionary course of STEC O157:H7.

Analysis Shiga Toxin-Encoding Bacteriophage in a Rare Strain of Shiga Toxin-Producing Escherichia coli O157:H7 stx2a/stx2c

In December 2015, six cases of Shiga toxin (Stx)-producing Escherichia coli (STEC) O157:H7 stx2a/stx2c phage type (PT) 24 were identified by the national gastrointestinal disease surveillance system at Public Health England (PHE). Frozen grated coconut imported from India was implicated as the vehicle of infection. Short and long read sequencing data were interrogated for genomic markers to provide evidence that the outbreak strain was from an imported source. The outbreak strain belonged to a sub-lineage (IIa) rare in domestically acquired infection in the United Kingdom, and indicative of an imported strain. Phylogenetic analysis identified the most closely related isolates to the outbreak strain were from cases reporting recent travel not to India, but to Uganda. Phylo-geographical signals based on travel data may be confounded by the failure of local and/or global monitoring systems to capture the full diversity of strains in a given country. This may be due to low prevalence strains circulating in-country under the surveillance radar, or a recent importation event involving the migration of animals and/or people. Comparison of stx2a-encoding prophage harbored by the outbreak strain with publicly available stx2a-encoding prophage sequences revealed that it was most closely related to stx2a-encoding prophage acquired by STEC O157:H7 that caused the first outbreak of STEC-hemolytic uremic syndrome (HUS) in England in 1982–83. Animal and people migration events may facilitate the transfer of stx2a-encoding prophage from indigenous STEC O157:H7 to recently imported strains, or vice versa. Monitoring the global transmission of STEC O157:H7 and tracking the exchange of stx2a-encoding phage between imported and indigenous strains may provide an early warning of emerging threats to public health.

Phages in Anaerobic Systems

Since the discovery of phages in 1915, these viruses have been studied mostly in aerobic systems, or without considering the availability of oxygen as a variable that may affect the interaction between the virus and its host. However, with such great abundance of anaerobic environments on the planet, the effect that a lack of oxygen can have on the phage-bacteria relationship is an important consideration. There are few studies on obligate anaerobes that investigate the role of anoxia in causing infection. In the case of facultative anaerobes, it is a well-known fact that their shifting from an aerobic environment to an anaerobic one involves metabolic changes in the bacteria. As the phage infection process depends on the metabolic state of the host bacteria, these changes are also expected to affect the phage infection cycle. This review summarizes the available information on phages active on facultative and obligate anaerobes and discusses how anaerobiosis can be an important parameter in phage infection, especially among facultative anaerobes.

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.

Comparison of Shiga toxin-encoding bacteriophages in highly pathogenic strains of Shiga toxin-producing Escherichia coli O157:H7 in the UK

Over the last 35 years in the UK, the burden of Shiga toxin-producing Escherichia coli (STEC) O157:H7 infection has, during different periods of time, been associated with five different sub-lineages (1983-1995, Ia, I/IIa and I/IIb; 1996-2014, Ic; and 2015-2018, IIb). The acquisition of a stx2a-encoding bacteriophage by these five sub-lineages appears to have coincided with their respective emergences. The Oxford Nanopore Technologies (ONT) system was used to sequence, characterize and compare the stx-encoding prophages harboured by each sub-lineage to investigate the integration of this key virulence factor. The stx2a-encoding prophages from each of the lineages causing clinical disease in the UK were all different, including the two UK sub-lineages (Ia and I/IIa) circulating concurrently and causing severe disease in the early 1980s. Comparisons between the stx2a-encoding prophage in sub-lineages I/IIb and IIb revealed similarity to the prophage commonly found to encode stx2c, and the same site of bacteriophage integration (sbcB) as stx2c-encoding prophage. These data suggest independent acquisition of previously unobserved stx2a-encoding phage is more likely to have contributed to the emergence of STEC O157:H7 sub-lineages in the UK than intra-UK lineage to lineage phage transmission. In contrast, the stx2c-encoding prophage showed a high level of similarity across lineages and time, consistent with the model of stx2c being present in the common ancestor to extant STEC O157:H7 and maintained by vertical inheritance in the majority of the population. Studying the nature of the stx-encoding bacteriophage contributes to our understanding of the emergence of highly pathogenic strains of STEC O157:H7.

Differential dynamics and impacts of prophages and plasmids on the pangenome and virulence factor repertoires of Shiga toxin-producing Escherichia coli O145:H28

Phages and plasmids play important roles in bacterial evolution and diversification. Although many draft genomes have been generated, phage and plasmid genomes are usually fragmented, limiting our understanding of their dynamics. Here, we performed a systematic analysis of 239 draft genomes and 7 complete genomes of Shiga toxin (Stx)-producing Escherichia coli O145:H28, the major virulence factors of which are encoded by prophages (PPs) or plasmids. The results indicated that PPs are more stably maintained than plasmids. A set of ancestrally acquired PPs was well conserved, while various PPs, including Stx phages, were acquired by multiple sublineages. In contrast, gains and losses of a wide range of plasmids have frequently occurred across the O145:H28 lineage, and only the virulence plasmid was well conserved. The different dynamics of PPs and plasmids have differentially impacted the pangenome of O145:H28, with high proportions of PP- and plasmid-associated genes in the variably present and rare gene fractions, respectively. The dynamics of PPs and plasmids have also strongly impacted virulence gene repertoires, such as the highly variable distribution of stx genes and the high conservation of a set of type III secretion effectors, which probably represents the core effectors of O145:H28 and the genes on the virulence plasmid in the entire O145:H28 population. These results provide detailed insights into the dynamics of PPs and plasmids, and show the application of genomic analyses using a large set of draft genomes and appropriately selected complete genomes.

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.

Genomic Characterization of β-Glucuronidase-Positive Escherichia coli O157:H7 Producing Stx2a

Among Shiga toxin (Stx)–producing Escherichia coli (STEC) O157:H7 strains, those producing Stx2a cause more severe diseases. Atypical STEC O157:H7 strains showing a β-glucuronidase–positive phenotype (GP STEC O157:H7) have rarely been isolated from humans, mostly from persons with asymptomatic or mild infections; Stx2a-producing strains have not been reported. We isolated, from a patient with bloody diarrhea, a GP STEC O157:H7 strain (PV15-279) that produces Stx2a in addition to Stx1a and Stx2c. Genomic comparison with other STEC O157 strains revealed that PV15-279 recently emerged from the stx1a/stx2c-positive GP STEC O157:H7 clone circulating in Japan. Major virulence genes are shared between typical (β-glucuronidase–negative) and GP STEC O157:H7 strains, and the Stx2-producing ability of PV15-279 is comparable to that of typical STEC O157:H7 strains; therefore, PV15-279 presents a virulence potential similar to that of typical STEC O157:H7. This study reveals the importance of GP O157:H7 as a source of highly pathogenic STEC clones.

Phylogenetic relationship of prophages is affected by CRISPR selection in Group A Streptococcus

BACKGROUND

Group A Streptococcus (GAS) is a major human pathogen, which is associated with a wide spectrum of invasive diseases, such as pharyngitis, scarlet fever, rheumatic fever, and streptococcal toxic shock syndrome (STSS). It is hypothesized that differences in GAS pathogenicity are related to the acquisition of diverse bacteriophages (phages). Nevertheless, the GAS genome also harbors clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (cas) genes, which play an important role in eliminating foreign DNA, including those of phages. However, the structure of prophages in GAS strains is mosaic, and the phylogenetic relationship between prophages and CRISPR is not clear. In this study, we analyzed CRISPR and prophage structure using 118 complete genome sequences of GAS strains to elucidate the relationship between two genomic elements. Additionally, phylogenetic and M-type analyses were performed.

RESULTS

Of the 118 GAS strains, 80 harbored type I-C and/or II-A CRISPR/cas loci. A total of 553 spacer sequences were identified from CRISPR/cas loci and sorted into 229 patterns. We identified and classified 373 prophages into 14 groups. Some prophage groups shared a common integration site, and were related to M-type. We further investigated the correlation between spacer sequences and prophages. Of the 229 spacer sequence patterns, 203 were similar to that of other GAS prophages. No spacer showed similarity with that of a specific prophage group with mutL integration site. Moreover, the average number of prophages in strains with type II-A CRISPR was significantly less than that in type I-C CRISPR and non-CRISPR strains. However, there was no statistical difference between the average number of prophages in type I-C strains and that in non-CRISPR strains.

CONCLUSIONS

Our results indicated that type II-A CRISPR may play an important role in eliminating phages and that the prophage integration site may be an important criterion for the acceptance of foreign DNA by GAS. M type, spacer sequence, and prophage group data were correlated with the phylogenetic relationships of GAS. Therefore, we hypothesize that genetic characteristics and/or phylogenetic relationships of GAS may be estimated by analyzing its spacer sequences.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

Transposon Insertion Sequencing Elucidates Novel Gene Involvement in Susceptibility and Resistance to Phages T4 and T7 in Escherichia coli O157

Experiments using bacteriophage (phage) to infect bacterial strains have helped define some basic genetic concepts in microbiology, but our understanding of the complexity of bacterium-phage interactions is still limited. As the global threat of antibiotic resistance continues to increase, phage therapy has reemerged as an attractive alternative or supplement to treating antibiotic-resistant bacterial infections. Further, the long-used method of phage typing to classify bacterial strains is being replaced by molecular genetic techniques. Thus, there is a growing need for a complete understanding of the precise molecular mechanisms underpinning phage-bacterium interactions to optimize phage therapy for the clinic as well as for retrospectively interpreting phage typing data on the molecular level. In this study, a genomics-based fitness assay (TraDIS) was used to identify all host genes involved in phage susceptibility and resistance for a T4 phage infecting Shiga-toxigenic Escherichia coli O157. The TraDIS results identified both established and previously unidentified genes involved in phage infection, and a subset were confirmed by site-directed mutagenesis and phenotypic testing of 14 T4 and 2 T7 phages. For the first time, the entire sap operon was implicated in phage susceptibility and, conversely, the stringent starvation protein A gene (sspA) was shown to provide phage resistance. Identifying genes involved in phage infection and replication should facilitate the selection of bespoke phage combinations to target specific bacterial pathogens.IMPORTANCE Antibiotic resistance has diminished treatment options for many common bacterial infections. Phage therapy is an alternative option that was once popularly used across Europe to kill bacteria within humans. Phage therapy acts by using highly specific viruses (called phages) that infect and lyse certain bacterial species to treat the infection. Whole-genome sequencing has allowed modernization of the investigations into phage-bacterium interactions. Here, using E. coli O157 and T4 bacteriophage as a model, we have exploited a genome-wide fitness assay to investigate all genes involved in defining phage resistance or susceptibility. This knowledge of the genetic determinants of phage resistance and susceptibility can be used to design bespoke phage combinations targeted to specific bacterial infections for successful infection eradication.

Evolutionary Context of Non-Sorbitol-Fermenting Shiga Toxin-Producing Escherichia coli O55:H7

In July 2014, an outbreak of Shiga toxin–producing Escherichia coli (STEC) O55:H7 in England involved 31 patients, 13 (42%) of whom had hemolytic uremic syndrome. Isolates were sequenced, and the sequences were compared with publicly available sequences of E. coli O55:H7 and O157:H7. A core-genome phylogeny of the evolutionary history of the STEC O55:H7 outbreak strain revealed that the most parsimonious model was a progenitor enteropathogenic O55:H7 sorbitol-fermenting strain, lysogenized by a Shiga toxin (Stx) 2a–encoding phage, followed by loss of the ability to ferment sorbitol because of a non-sense mutation in srlA. The parallel, convergent evolutionary histories of STEC O157:H7 and STEC O55:H7 may indicate a common driver in the evolutionary process. Because emergence of STEC O157:H7 as a clinically significant pathogen was associated with acquisition of the Stx2a-encoding phage, the emergence of STEC O55:H7 harboring the stx2a gene is of public health concern.

Enterohaemorrhagic and other Shiga toxin-producing Escherichia coli (STEC): Where are we now regarding diagnostics and control strategies?

Escherichia coli comprises a highly diverse group of Gram-negative bacteria and is a common member of the intestinal microflora of humans and animals. Generally, such colonization is asymptomatic; however, some E. coli strains have evolved to become pathogenic and thus cause clinical disease in susceptible hosts. One pathotype, the Shiga toxigenic E. coli (STEC) comprising strains expressing a Shiga-like toxin is an important foodborne pathogen. A subset of STEC are the enterohaemorrhagic E. coli (EHEC), which can cause serious human disease, including haemolytic uraemic syndrome (HUS). The diagnosis of EHEC infections and the surveillance of STEC in the food chain and the environment require accurate, cost-effective and timely tests. In this review, we describe and evaluate tests now in routine use, as well as upcoming test technologies for pathogen detection, including loop-mediated isothermal amplification (LAMP) and whole-genome sequencing (WGS). We have considered the need for improved diagnostic tools in current strategies for the control and prevention of these pathogens in humans, the food chain and the environment. We conclude that although significant progress has been made, STEC still remains an important zoonotic issue worldwide. Substantial reductions in the public health burden due to this infection will require a multipronged approach, including ongoing surveillance with high-resolution diagnostic techniques currently being developed and integrated into the routine investigations of public health laboratories. However, additional research requirements may be needed before such high-resolution diagnostic tools can be used to enable the development of appropriate interventions, such as vaccines and decontamination strategies.