Whole genome sequence comparison of vtx2-converting phages from Enteroaggregative Haemorrhagic Escherichia coli strains

Transduction of stx2a mediated by phage (Φ11-3088) from Escherichia coli O104:H4 in vitro and in situ during sprouting of mung beans

Escherichia coli O104:H4 strain 11-3088 encoding Stx2a is epidemiologically related to the foodborne outbreak associated with sprouts in Germany, 2011. Sprouting provides suitable conditions for bacterial growth and may lead to transduction of non-pathogenic strains of E. coli with Stx phages. Although transduction of E. coli by Stx phages in food has been documented, data on the phages from E. coli O104:H4 is limited. This study determined the host range of the bacteriophage Φ11-3088 from E. coli O104:H4 using E. coli O104:H4 ∆stx2::gfp::amp^r and demonstrated phage transduction during sprouting. The Φ11-3088∆stx transduced 5/45 strains, including generic E. coli, pap-positive E. coli O103:H2, ETEC, and S. sonnei. The expression level of Φ11-3088∆stx differed among lysogens upon induction. Of the 3 highly induced lysogens, the lytic cycle was induced in E. coli O104:H4∆stx2::gfp::amp^r and O103:H2 but not in S. sonnei. E. coli DH5α was the only strain susceptible to lytic infection by Φ11-3088∆stx. To explore the effect of drying and rehydration during seed storage and sprouting on phage induction and transduction, mung beans inoculated with the phage donor E. coli O104:H4∆stx2::gfp::amp^r (8 log CFU/g) were dried, rehydrated, and incubated with the phage recipient E. coli DH5α (7 log CFU/g) for 96 h. Sprouted seeds harbored about 3 log CFU/g of putative lysogens that acquired ampicillin resistance. At the end of sprouting, 71 % of putative lysogens encoded gfp, confirming phage transduction. Overall, stx transfer by phages may increase the cell counts of STEC during sprouting by converting generic E. coli to STEC.

Replication Region Analysis Reveals Non-lambdoid Shiga Toxin Converting Bacteriophages

Shiga toxin is the major virulence factor of enterohemorrhagic Escherichia coli (EHEC), and the gene encoding it is carried within the genome of Shiga toxin-converting phages (Stx phages). Numerous Stx phages have been sequenced to gain a better understanding of their contribution to the virulence potential of EHEC. The Stx phages are classified into the lambdoid phage family based on similarities in lifestyle, gene arrangement, and nucleotide sequence to the lambda phages. This study explores the replication regions of non-lambdoid Stx phages that completely lack the O and P genes encoding the proteins involved in initiating replication in the lambdoid phage genome. Instead, they carry sequences encoding replication proteins that have not been described earlier, here referred to as eru genes (after EHEC phage replication unit genes). This study identified three different types of Eru-phages, where the Eru1-type is carried by the highly pathogenic EHEC strains that caused the Norwegian O103:H25 outbreak in 2006 and the O104:H4 strain that caused the large outbreak in Europe in 2011. We show that Eru1-phages exhibit a less stable lysogenic state than the classical lambdoid Stx phages. As production of phage particles is accompanied by production of Stx toxin, the Eru1-phage could be associated with a high-virulence phenotype of the host EHEC strain. This finding emphasizes the importance of classifying Stx phages according to their replication regions in addition to their Stx-type and could be used to develop a novel strategy to identify highly virulent EHEC strains for improved risk assessment and management.

Isolation and Characterization of Shiga Toxin Bacteriophages

Shiga toxin (Stx) phages can be induced from Stx-producing Escherichia coli strains (STEC) or can be isolated as free virions from different samples. Here we describe methods used for the detection, enumeration, and isolation of Stx bacteriophages. Stx phages are temperate phages located in the genome of STEC. Their induction from the host strain cultures is achieved by different inducing agents, mitomycin C being one of the most commonly used. Detection of infectious Stx phages requires the production of visible plaques in a confluent lawn of the host strain using a double agar layer method. However, as the plaques produced by Stx phages are often barely visible and there is a possibility that non-Stx phages can also be induced from the strain, a hybridization step should be added to recognize and properly enumerate the lysis plaques generated after induction. Molecular methods can also be used to identify and enumerate Stx phages. Real-time quantitative PCR (qPCR) is the most accurate method for absolute quantification, although it cannot determine the infectivity of Stx phages. qPCR can also be useful for the detection of free Stx phage virions in different sample types.Stx phages induced from lysogenic bacterial strains can be purified by cesium chloride density gradients; this protocol also helps to specifically discriminate Stx phages from other prophages present in the genome of the host strain by selecting the phages expressing the Stx gene. High titer suspensions of Stx phages obtained after induction of large volumes of bacterial cultures and lysate concentration permits phage characterization by electron microscopy studies and genomic analysis.

Genomic Characterization of Two Shiga Toxin-Converting Bacteriophages Induced From Environmental Shiga Toxin-Producing Escherichia coli

Shiga toxin (Stx), encoded by stx genes located in prophage sequences, is the major agent responsible for the pathogenicity of Shiga toxin-producing Escherichia coli (STEC) and is closely associated with the development of hemolytic uremic syndrome (HUS). Although numerous Stx prophage sequences have been reported as part of STEC bacterial genomes, the information about the genomic characterization of Stx-converting bacteriophages induced from STEC strains is relatively scarce. The objectives of this study were to genomically characterize two Stx-converting phages induced from environmental STEC strains and to evaluate their correlations with published Stx-converting phages and STEC strains of different origins. The Stx1-converting phage Lys8385Vzw and the Stx2-converting phage Lys19259Vzw were induced from E. coli O103:H11 (RM8385) and E. coli O157:H7 (RM19259), respectively. Whole-genome sequencing of these phages was conducted on a MiSeq sequencer for genomic characterization. Phylogenetic analysis and comparative genomics were performed to determine the correlations between these two Stx-converting phages, 13 reference Stx-converting phages, and 10 reference STEC genomes carrying closely related Stx prophages. Both Stx-converting phages Lys8385Vzw and Lys19259Vzw had double-stranded DNA, with genome sizes of 50,953 and 61,072 bp, respectively. Approximately 40% of the annotated coding DNA sequences with the predicted functions were likely associated with the fitness for both phages and their bacterial hosts. The whole-genome–based phylogenetic analysis of these two Stx-converting phages and 13 reference Stx-converting phages revealed that the 15 Stx-converting phages were divided into three distinct clusters, and those from E. coli O157:H7, in particular, were distributed in each cluster, demonstrating the high genomic diversity of these Stx-converting phages. The genomes of Stx-converting phage Lys8385Vzw and Lys19259Vzw shared a high-nucleotide similarity with the prophage sequences of the selected STEC isolates from the clinical and environmental origin. The findings demonstrate the genomic diversity of Stx-converting phages induced from different STEC strains and provide valuable insights into the dissemination of stx genes among E. coli population via the lysogenization of Stx-converting phages.

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

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

Characterization and comparative genomic analysis of virulent and temperate Bacillus megaterium bacteriophages

Next-Generation Sequencing (NGS) technologies provide unique possibilities for the comprehensive assessment of the environmental diversity of bacteriophages. Several Bacillus bacteriophages have been isolated, but very few Bacillus megaterium bacteriophages have been characterized. In this study, we describe the biological characteristics, whole genome sequences, and annotations for two new isolates of the B. megaterium bacteriophages (BM5 and BM10), which were isolated from Egyptian soil samples. Growth analyses indicated that the phages BM5 and BM10 have a shorter latent period (25 and 30 min, respectively) and a smaller burst size (103 and 117 PFU, respectively), in comparison to what is typical for Bacillus phages. The genome sizes of the phages BM5 and BM10 were 165,031 bp and 165,213 bp, respectively, with modular organization. Bioinformatic analyses of these genomes enabled the assignment of putative functions to 97 and 65 putative ORFs, respectively. Comparative analysis of the BM5 and BM10 genome structures, in conjunction with other B. megaterium bacteriophages, revealed relatively high levels of sequence and organizational identity. Both genomic comparisons and phylogenetic analyses support the conclusion that the sequenced phages (BM5 and BM10) belong to different sub-clusters (L5 and L7, respectively), within the L-cluster, and display different lifestyles (lysogenic and lytic, respectively). Moreover, sequenced phages encode proteins associated with Bacillus pathogenesis. In addition, BM5 does not contain any tRNA sequences, whereas BM10 genome codes for 17 tRNAs.

ArgO145, a Stx2a prophage of a bovine O145:H- STEC strain, is closely related to phages of virulent human strains

Shiga toxins (Stx) are the main virulence factor of a pathogroup of Escherichia coli strains that cause severe human diseases. These toxins are encoded in prophages (Stx prophages), and generally their expression depends on prophage induction. Several studies have reported high diversity among both Stx prophages and Stx. In particular, the toxin subtype Stx2a is associated with high virulence and HUS. Here, we report the genome of ArgO145, an inducible Stx2a prophage identified in a bovine O145:H- strain which produced high levels of Shiga toxin and Stx phage particles. The ArgO145 genome shared lambda phage organization, with recombination, regulation, replication, lysis, and head and tail structural gene regions, although some lambda genes encoding regulatory proteins could not be identified. Remarkably, some Stx2a phages of strains isolated from patients in other countries showed high similarity to ArgO145.

Functional and comparative genome analysis of novel virulent actinophages belonging to Streptomyces flavovirens

Background

Next Generation Sequencing (NGS) technologies provide exciting possibilities for whole genome sequencing of a plethora of organisms including bacterial strains and phages, with many possible applications in research and diagnostics. No Streptomyces flavovirens phages have been sequenced to date; there is therefore a lack in available information about S. flavovirens phage genomics. We report biological and physiochemical features and use NGS to provide the complete annotated genomes for two new strains (Sf1 and Sf3) of the virulent phage Streptomyces flavovirens, isolated from Egyptian soil samples.

Results

The S. flavovirens phages (Sf1 and Sf3) examined in this study show higher adsorption rates (82 and 85%, respectively) than other actinophages, indicating a strong specificity to their host, and latent periods (15 and 30 min.), followed by rise periods of 45 and 30 min. As expected for actinophages, their burst sizes were 1.95 and 2.49 virions per mL. Both phages were stable and, as reported in previous experiments, showed a significant increase in their activity after sodium chloride (NaCl) and magnesium chloride (MgCl₂.6H₂O) treatments, whereas after zinc chloride (ZnCl2) application both phages showed a significant decrease in infection.

The sequenced phage genomes are parts of a singleton cluster with sizes of 43,150 bp and 60,934 bp, respectively. Bioinformatics analyses and functional characterizations enabled the assignment of possible functions to 19 and 28 putative identified ORFs, which included phage structural proteins, lysis components and metabolic proteins.

Thirty phams were identified in both phages, 10 (33.3%) of them with known function, which can be used in cluster prediction. Comparative genomic analysis revealed significant homology between the two phages, showing the highest hits among Sf1, Sf3 and the closest Streptomyces phage (VWB phages) in a specific 13Kb region. However, the phylogenetic analysis using the Major Capsid Protein (MCP) sequences highlighted that the isolated phages belong to the BG Streptomyces phage group but are clearly separated, representing a novel sub-cluster.

Conclusion

The results of this study provide the first physiological and genomic information for S. flavovirens phages and will be useful for pharmaceutical industries based on S. flavovirens and future phage evolution studies.

Multiple-locus variable-number tandem repeat analysis for strain discrimination of non-O157 Shiga toxin-producing Escherichia coli

Non-O157 Shiga toxin-producing Escherichia coli (STEC) are foodborne pathogens of growing concern worldwide that have been associated with several recent multistate and multinational outbreaks of foodborne illness. Rapid and sensitive molecular-based bacterial strain discrimination methods are critical for timely outbreak identification and contaminated food source traceback. One such method, multiple-locus variable-number tandem repeat analysis (MLVA), is being used with increasing frequency in foodborne illness outbreak investigations to augment the current gold standard bacterial subtyping technique, pulsed-field gel electrophoresis (PFGE). The objective of this study was to develop a MLVA assay for intra- and inter-serogroup discrimination of six major non-O157 STEC serogroups-O26, O111, O103, O121, O45, and O145-and perform a preliminary internal validation of the method on a limited number of clinical isolates. The resultant MLVA scheme consists of ten variable number tandem repeat (VNTR) loci amplified in three multiplex PCR reactions. Sixty-five unique MLVA types were obtained among 84 clinical non-O157 STEC strains comprised of geographically diverse sporadic and outbreak related isolates. Compared to PFGE, the developed MLVA scheme allowed similar discrimination among serogroups O26, O111, O103, and O121 but not among O145 and O45. To more fully compare the discriminatory power of this preliminary MLVA method to PFGE and to determine its epidemiological congruence, a thorough internal and external validation needs to be performed on a carefully selected large panel of strains, including multiple isolates from single outbreaks.

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.

Genome sequencing and comparative genomics provides insights on the evolutionary dynamics and pathogenic potential of different H-serotypes of Shiga toxin-producing Escherichia coli O104

BACKGROUND

Various H-serotypes of the Shiga toxin-producing Escherichia coli (STEC) O104, including H4, H7, H21, and H¯, have been associated with sporadic cases of illness and have caused food-borne outbreaks globally. In the U.S., STEC O104:H21 caused an outbreak associated with milk in 1994. However, there is little known on the evolutionary origins of STEC O104 strains, and how genotypic diversity contributes to pathogenic potential of various O104 H-antigen serotypes isolated from different ecological niches and/or geographical regions.

RESULTS

Two STEC O104:H21 (milk outbreak strain) and O104:H7 (cattle isolate) strains were shot-gun sequenced, and the genomes were closed. The intimin (eae) gene, involved in the attaching-effacing phenotype of diarrheagenic E. coli, was not found in either strain. Examining various O104 genome sequences, we found that two "complete" left and right end portions of the locus of enterocyte effacement (LEE) pathogenicity island were present in 13 O104 strains; however, the central portion of LEE was missing, where the eae gene is located. In O104:H4 strains, the missing central portion of the LEE locus was replaced by a pathogenicity island carrying the aidA (adhesin involved in diffuse adherence) gene and antibiotic resistance genes commonly carried on plasmids. Enteroaggregative E. coli-specific virulence genes and European outbreak O104:H4-specific stx2-encoding Escherichia P13374 or Escherichia TL-2011c bacteriophages were missing in some of the O104:H4 genome sequences available from public databases. Most of the genomic variations in the strains examined were due to the presence of different mobile genetic elements, including prophages and genomic island regions. The presence of plasmids carrying virulence-associated genes may play a role in the pathogenic potential of O104 strains.

CONCLUSIONS

The two strains sequenced in this study (O104:H21 and O104:H7) are genetically more similar to each other than to the O104:H4 strains that caused an outbreak in Germany in 2011 and strains found in Central Africa. A hypothesis on strain evolution and pathogenic potential of various H-serotypes of E. coli O104 strains is proposed.