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

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

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

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

Insertion Sequence (IS)-Excision Enhancer (IEE)-Mediated IS Excision from the lacZ Gene Restores the Lactose Utilization Defect of Shiga Toxin-Producing Escherichia coli O121:H19 Strains and Is Responsible for Their Delayed Lactose Utilization Phenotype

Lactose utilization is one of the general biochemical characteristics of Escherichia coli, and the lac operon is responsible for this phenotype, which can be detected on lactose-containing media, such as MacConkey agar, after 24 h of incubation. However, some Shiga toxin-producing E. coli (STEC) O121:H19 strains exhibit an unusual phenotype called delayed lactose utilization (DLU), in which lactose utilization can be detected after 48 h of cultivation but not after only 24 h of cultivation. Insertion of an insertion sequence (IS), IS600, into the lacZ gene appears to be responsible for the DLU phenotype, and exposure to lactose has been reported to be necessary to observe this phenotype, but the mechanism underlying these phenomena remains to be elucidated. Here, we performed detailed analyses of the lactose utilization abilities of a set of O121:H19 strains and their mutants and found that IS-excision enhancer (IEE)-mediated excision of IS600 reactivates the lacZ gene and that the selective proliferation of IS-cured subclones in lactose-supplemented culture medium is responsible for the expression of the DLU phenotype. In addition, we analyzed the patterns of IS insertion into the lacZ and iee genes in the global O121:H19 population and revealed that while there are O121:H19 strains or lineage/sublineages that contain the IS insertion into iee or intact lacZ and thus do not show the DLU phenotype, most currently circulating O121:H19 strains contain IS600-inserted lacZ and intact iee and thus exhibit this phenotype.

IMPORTANCE Insertion sequences (ISs) can modulate gene expression by gene inactivation or activation. While phenotypic changes due to IS insertion/transposition are frequently observed, gene reactivation by precise or simple IS excision rarely occurs. In this study, we show that IS600 is excised from the lacZ gene by IS-excision enhancer (IEE) during the cultivation of Shiga toxin-producing Escherichia coli (STEC) O121:H19 strains that show an unusual phenotype called delayed lactose utilization (DLU). This excision rescued their lactose utilization defect, and the subsequent selective proliferation of IS-cured subclones in lactose-containing medium resulted in the expression of the DLU phenotype. As we also show that most currently circulating O121:H19 strains exhibit this phenotype, this study not only provides information helpful for the isolation and identification of O121:H19 STEC but also offers novel insights into the roles of IS and IEE in the generation of phenotypic variation in bacterial populations.

Chronological set of E. coli O157:H7 bovine strains establishes a role for repeat sequences and mobile genetic elements in genome diversification

Background

Enterohemorrhagic Escherichia coli O157:H7 (EHEC) is a significant foodborne pathogen that resides asymptomatically within cattle and other ruminants. The EHEC genome harbors an extensive collection of mobile genetic elements (MGE), including multiple prophage, prophage-like elements, plasmids, and insertion sequence (IS) elements.

Results

A chronological collection of EHEC strains (FRIK804, FRIK1275, and FRIK1625) isolated from a Wisconsin dairy farm (farm X) comprised a closely related clade genetically differentiated by structural alterations to the chromosome. Comparison of the FRIK804 genome with a reference EHEC strain Sakai found a unique prophage like element (PLE, indel 1) and an inversion (1.15 Mb) situated symmetrically with respect to the terminus region. Detailed analysis determined the inversion was due to homologous recombination between repeat sequences in prophage. The three farm X strains were distinguished by the presence or absence of indel 3 (61 kbp) and indel 4 (48 kbp); FRIK804 contained both of these regions, FRIK1275 lacked indel 4, and indels 3 and 4 were both absent in FRIK1625. Indel 3 was the stx2 prophage and indel 4 involved a deletion between two adjacent prophage with shared repeat sequences. Both FRIK804 and FRIK1275 produced functional phage while FRIK1625 did not, which is consistent with indel 3. Due to their involvement in recombination events, direct and inverted repeat sequences were identified, and their locations mapped to the chromosome. FRIK804 had a greater number and overall length of repeat sequences than E. coli K12 strain MG1655. Repeat sequences were most commonly associated with MGE.

Conclusions

This research demonstrated that three EHEC strains from a Wisconsin dairy farm were closely related and distinguished by variability within prophage regions and other MGE. Chromosome alterations were associated with recombination events between repeat sequences. An inventory of direct and inverted repeat sequences found a greater abundance and total length of repeat sequences in the EHEC strains compared to E. coli strain MG1655. The locations of the repeat sequences were biased towards MGE. The findings from this study expand our understanding of the precise molecular events and elements that contributed to genetic diversification of wild-type EHEC in the bovine and farm environments.

Genetic and Functional Analyses of Virulence Potential of an Escherichia coli O157:H7 Strain Isolated From Super-Shedder Cattle

Shiga toxin (Stx)-producing Escherichia coli (STEC) O157:H7 is an enteric pathogen that causes life-threatening disease in humans, with cattle being major natural reservoirs. A group of STEC O157:H7 with a dramatic combination of high virulence potentials and super-shedder bovine origin have been isolated. Here, an STEC O157:H7 isolate, JEONG-1266, was analyzed by comparative genomics, stx genotyping, and phenotypic analyses. The phylogenetic typing and whole-genome comparison consistently showed that JEONG-1266 is genetically close to EC4115 (one of 2006 Spinach outbreak isolates) and SS17 (an isolate from super-shedder cattle) strains, all of which belong to lineage I/II and Clade 8. Both lineage I/II and Clade 8 are known to be mostly associated with clinical strains with high virulence and severe clinical symptoms. Further, JEONG-1266, like EC4115 and SS17, harbors stx2a/stx2c genes, and carries Stx-encoding prophages, specifically the φstx2a-γ subtype. Possession of the φstx2a-γ subtype of Stx-encoding prophages and production of Stx2a have been shown to be a key signature associated with hypervirulent STEC O157:H7 strains. In silico virulence typing elucidated JEONG-1266, EC4115, and SS17 shared a highly conserved profile of key virulence genes at the nucleotide sequence level. Consistently, phenotypic data showed that JEONG-1266 expressed a high level of Stx2 toxins and had the full capacity of adhesion in vitro. Taken together, our study suggests that JEONG-1266 may represent an emerging STEC O157:H7 group, which are hypervirulent strains that originate from super-shedders, that can be a threat to food safety and public health.

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.

Colonization of Beef Cattle by Shiga Toxin-Producing Escherichia coli during the First Year of Life: A Cohort Study

Each year Shiga toxin-producing Escherichia coli (STEC) are responsible for 2.8 million acute illnesses around the world and > 250,000 cases in the US. Lowering the prevalence of this pathogen in animal reservoirs has the potential to reduce STEC outbreaks in humans by controlling its entrance into the food chain. However, factors that modulate the colonization and persistence of STEC in beef cattle remain largely unidentified. This study evaluated if animal physiological factors such as age, breed, sex, and weight gain influenced the shedding of STEC in beef cattle. A cohort of beef calves (n = 260) from a multi-breed beef calf population was sampled every three months after birth to measure prevalence and concentration of STEC during the first year of life. Metagenomic analysis was also used to understand the association between the STEC colonization and the composition of gut microflora. This study identified that beef calves were more likely to shed STEC during the first 6 months and that STEC shedding decreased as the animal matured. Animal breed group, sex of the calf, and average weight gain were not significantly associated with STEC colonization. The metagenomic analysis revealed for the first time that STEC colonization was correlated with a lower diversity of gut microflora, which increases as the cattle matured. Given these findings, intervention strategies that segregate younger animals, more likely to be colonized by STEC from older animals that are ready to be harvested, could be investigated as a method to reduce zoonotic transmission of STEC from cattle to humans.

Escherichia coli O157:H7 strains harbor at least three distinct sequence types of Shiga toxin 2a-converting phages

Background

Shiga toxin-producing Escherichia coli O157:H7 is a foodborne pathogen that causes severe human diseases including hemolytic uremic syndrome (HUS). The virulence factor that mediates HUS, Shiga toxin (Stx), is encoded within the genome of a lambdoid prophage. Although draft sequences are publicly available for a large number of E. coli O157:H7 strains, the high sequence similarity of stx-converting bacteriophages with other lambdoid prophages poses challenges to accurately assess the organization and plasticity among stx-converting phages due to assembly difficulties.

Methods

To further explore genome plasticity of stx-converting prophages, we enriched phage DNA from 45 ciprofloxacin-induced cultures for subsequent 454 pyrosequencing to facilitate assembly of the complete phage genomes. In total, 22 stx2a-converting phage genomes were closed.

Results

Comparison of the genomes distinguished nine distinct phage sequence types (PSTs) delineated by variation in obtained sequences, such as single nucleotide polymorphisms (SNPs) and insertion sequence element prevalence and location. These nine PSTs formed three distinct clusters, designated as PST1, PST2 and PST3. The PST2 cluster, identified in two clade 8 strains, was related to stx2a-converting phages previously identified in non-O157 Shiga-toxin producing E. coli (STEC) strains associated with a high incidence of HUS. The PST1 cluster contained phages related to those from E. coli O157:H7 strain Sakai (lineage I, clade 1), and PST3 contained a single phage that was distinct from the rest but most related to the phage from E. coli O157:H7 strain EC4115 (lineage I/II, clade 8). Five strains carried identical stx2a-converting phages (PST1-1) integrated at the same chromosomal locus, but these strains produced different levels of Stx2.

Conclusion

The stx2a-converting phages of E. coli O157:H7 can be categorized into at least three phage types. Diversification within a phage type is mainly driven by IS629 and by a small number of SNPs. Polymorphisms between phage genomes may help explain differences in Stx2a production between strains, however our data indicates that genes encoded external to the phage affect toxin production as well.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1934-1) contains supplementary material, which is available to authorized users.

Simultaneous Presence of Insertion Sequence Excision Enhancer and Insertion Sequence IS629 Correlates with Increased Diversity and Virulence in Shiga Toxin-Producing Escherichia coli

Although new serotypes of enterohemorrhagic Escherichia coli (EHEC) emerge constantly, the mechanisms by which these new pathogens arise and the reasons emerging serotypes tend to carry more virulence genes than other E. coli are not understood. An insertion sequence (IS) excision enhancer (IEE) was discovered in EHEC O157:H7 that promoted the excision of IS3 family members and generating various genomic deletions. One IS3 family member, IS629, actively transposes and proliferates in EHEC O157:H7 and enterotoxigenic E. coli (ETEC) O139 and O149. The simultaneous presence of the IEE and IS629 (and other IS3 family members) may be part of a system promoting not only adaptation and genome diversification in E. coli O157:H7 but also contributing to the development of pathogenicity among predominant serotypes. Prevalence comparisons of these elements in 461 strains, representing 72 different serotypes and 5 preassigned seropathotypes (SPT) A to E, showed that the presence of these two elements simultaneously was serotype specific and associated with highly pathogenic serotypes (O157 and top non-O157 Shiga toxin-producing Escherichia coli [STEC]) implicated in outbreaks and sporadic cases of human illness (SPT A and B). Serotypes lacking one or both elements were less likely to have been isolated from clinical cases. Our comparisons of IEE sequences showed sequence variations that could be divided into at least three clusters. Interestingly, the IEE sequences from O157 and the top 10 non-O157 STEC serotypes fell into clusters I and II, while less commonly isolated serotypes O5 and O174 fell into cluster III. These results suggest that IS629 and IEE elements may be acting synergistically to promote genome plasticity and genetic diversity among STEC strains, enhancing their abilities to adapt to hostile environments and rapidly take up virulence factors.

Comparative analysis of super-shedder strains of Escherichia coli O157:H7 reveals distinctive genomic features and a strongly aggregative adherent phenotype on bovine rectoanal junction squamous epithelial cells

Shiga toxin-producing Escherichia coli O157:H7 (O157) are significant foodborne pathogens and pose a serious threat to public health worldwide. The major reservoirs of O157 are asymptomatic cattle which harbor the organism in the terminal recto-anal junction (RAJ). Some colonized animals, referred to as “super-shedders” (SS), are known to shed O157 in exceptionally large numbers (>104 CFU/g of feces). Recent studies suggest that SS cattle play a major role in the prevalence and transmission of O157, but little is known about the molecular mechanisms associated with super-shedding. Whole genome sequence analysis of an SS O157 strain (SS17) revealed a genome of 5,523,849 bp chromosome with 5,430 open reading frames and two plasmids, pO157 and pSS17, of 94,645 bp and 37,446 bp, respectively. Comparative analyses showed that SS17 is clustered with spinach-associated O157 outbreak strains, and belongs to the lineage I/II, clade 8, D group, and genotype 1, a subgroup of O157 with predicted hyper-virulence. A large number of non-synonymous SNPs and other polymorphisms were identified in SS17 as compared with other O157 strains (EC4115, EDL933, Sakai, TW14359), including in key adherence- and virulence-related loci. Phenotypic analyses revealed a distinctive and strongly adherent aggregative phenotype of SS17 on bovine RAJ stratified squamous epithelial (RSE) cells that was conserved amongst other SS isolates. Molecular genetic and functional analyses of defined mutants of SS17 suggested that the strongly adherent aggregative phenotype amongst SS isolates is LEE-independent, and likely results from a novel mechanism. Taken together, our study provides a rational framework for investigating the molecular mechanisms associated with SS, and strong evidence that SS O157 isolates have distinctive features and use a LEE-independent mechanism for hyper-adherence to bovine rectal epithelial cells.

Genetic diversity of Shiga toxin-producing Escherichia coli O157 : H7 recovered from human and food sources

The aim of this study was to identify an epidemiological association between Shiga toxin-producing Escherichia coli O157 : H7 strains associated with human infection and with food sources. Frequency distributions of different genetic markers of E. coli O157 : H7 strains recovered from human and food sources were compared using molecular assays to identify E. coli O157 : H7 genotypes associated with variation in pathogenic potential and host specificity. Genotypic characterization included: lineage-specific polymorphism assay (LSPA-6), clade typing, tir (A255T) polymorphism, Shiga toxin-encoding bacteriophage insertion site analysis and variant analysis of Shiga toxin 2 gene (stx2a and stx2c) and antiterminator Q genes (Q933 and Q21). The intermediate lineage (LI/II) dominated among both food and human strains. Compared to other clades, clades 7 and 8 were more frequent among food and human strains, respectively. The tir (255T) polymorphism occurred more frequently among human strains than food strains. Q21 and Q933 + Q21 were found at significantly higher frequencies among food and human strains, respectively. Moreover, stx2a and stx2a+c were detected at significantly higher frequencies among human strains compared to food strains. Bivariate analysis revealed significant concordance (P<0.05) between the LSPA-6 assay and the other typing methods. Multivariable regression analysis suggested that tir (255T) was the most distinctive genotype that can be used to detect bacterial clones with potential risk for human illness from food sources. This study supported previous reports of the existence of diversity in genetic markers among different isolation sources by including E. coli O157 : H7 strains from both food and human sources. This might enable tracking genotypes with potential risk for human illness from food sources.