The nucleotide sequence of Shiga toxin (Stx) 2e-encoding phage phiP27 is not related to other Stx phage genomes, but the modular genetic structure is conserved

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.

Phage tailspike modularity and horizontal gene transfer reveals specificity towards E. coli O-antigen serogroups

BACKGROUND

The interaction between bacteriophages and their hosts is intricate and highly specific. Receptor-binding proteins (RBPs) of phages such as tail fibers and tailspikes initiate the infection process. These RBPs bind to diverse outer membrane structures, including the O-antigen, which is a serogroup-specific sugar-based component of the outer lipopolysaccharide layer of Gram-negative bacteria. Among the most virulent Escherichia coli strains is the Shiga toxin-producing E. coli (STEC) pathotype dominated by a subset of O-antigen serogroups.

METHODS

Extensive phylogenetic and structural analyses were used to identify and validate specificity correlations between phage RBP subtypes and STEC O-antigen serogroups, relying on the principle of horizontal gene transfer as main driver for RBP evolution.

RESULTS

We identified O-antigen specific RBP subtypes for seven out of nine most prevalent STEC serogroups (O26, O45, O103, O104, O111, O145 and O157) and seven additional E. coli serogroups (O2, O8, O16, O18, 4s/O22, O77 and O78). Eight phage genera (Gamaleya-, Justusliebig-, Kaguna-, Kayfuna-, Kutter-, Lederberg-, Nouzilly- and Uetakeviruses) emerged for their high proportion of serogroup-specific RBPs. Additionally, we reveal sequence motifs in the RBP region, potentially serving as recombination hotspots between lytic phages.

CONCLUSION

The results contribute to a better understanding of mosaicism of phage RBPs, but also demonstrate a method to identify and validate new RBP subtypes for current and future emerging serogroups.

Characteristics and Comparative Genomic Analysis of a Novel Virus, VarioGold, the First Bacteriophage of Variovorax

Variovorax represents a widespread and ecologically significant genus of soil bacteria. Despite the ecological importance of these bacteria, our knowledge about the viruses infecting Variovorax spp. is quite poor. This study describes the isolation and characterization of the mitomycin-induced phage, named VarioGold. To the best of our knowledge, VarioGold represents the first characterized virus for this genus. Comparative genomic analyses suggested that VarioGold is distinct from currently known bacteriophages at both the nucleotide and protein levels; thus, it could be considered a new virus genus. In addition, another 37 prophages were distinguished in silico within the complete genomic sequences of Variovorax spp. that are available in public databases. The similarity networking analysis highlighted their general high diversity, which, despite clustering with previously described phages, shows their unique genetic load. Therefore, the novelty of Variovorax phages warrants the great enrichment of databases, which could, in turn, improve bioinformatic strategies for finding (pro)phages.

Comparative Genomics of Escherichia coli Serotype O55:H7 Using Complete Closed Genomes

Escherichia coli O55:H7 is a human foodborne pathogen and is recognized as the progenitor strain of E. coli O157:H7. While this strain is important from a food safety and genomic evolution standpoint, much of the genomic diversity of E. coli O55:H7 has been demonstrated using draft genomes. Here, we combine the four publicly available E. coli O55:H7 closed genomes with six newly sequenced closed genomes to provide context to this strain’s genomic diversity. We found significant diversity within the 10 E. coli O55:H7 strains that belonged to three different sequence types. The prophage content was about 10% of the genome, with three prophages common to all strains and seven unique to one strain. Overall, there were 492 insertion sequences identified within the six new sequence strains, with each strain on average containing 75 insertions (range 55 to 114). A total of 31 plasmids were identified between all isolates (range 1 to 6), with one plasmid (pO55) having an identical phylogenetic tree as the chromosome. The release and comparison of these closed genomes provides new insight into E. coli O55:H7 diversity and its ability to cause disease in humans.

Characterization of Shiga Toxin 2a Encoding Bacteriophages Isolated From High-Virulent O145:H25 Shiga Toxin-Producing Escherichia coli

Shiga toxin-producing Escherichia coli (STEC) may cause severe disease mainly due to the ability to produce Shiga toxins (Stx) encoded on bacteriophages. In Norway, more than 30% of the reported cases with STEC O145:H25 develop hemolytic uremic syndrome (HUS), and most cases, with known travel history, acquired the infection domestically. To describe phage characteristics associated with high virulence, we extracted the Stx2a phage sequences from eight clinical Norwegian O145:H25 STEC to conduct in-depth molecular characterization using long and short read sequencing. The Stx2a phages were annotated, characterized, and compared with previously published Stx2a phages isolated from STEC of different serotypes. The Norwegian O145:H25 Stx2a phages showed high sequence identity (>99%) with 100% coverage. The Stx2a phages were located at the integration site yciD, were approximately 45 kbp long, and harbored several virulence-associated genes, in addition to stx2a, such as nanS and nleC. We observed high sequence identity (>98%) and coverage (≥94%) between Norwegian O145:H25 Stx2a phages and publicly available Stx2a phages from O145:H25 and O145:H28 STEC, isolated from HUS cases in the USA and a hemorrhagic diarrhea case from Japan, respectively. However, low similarity was seen when comparing the Norwegian O145:H25 Stx2a phage to Stx2a phages from STEC of other serotypes. In all the Norwegian O145:H25 STEC, we identified a second phage or remnants of a phage (a shadow phage, 61 kbp) inserted at the same integration site as the Stx2a phage. The shadow phage shared similarity with the Stx2a phage, but lacked stx2a and harbored effector genes not present in the Stx2a phage. We identified a conserved Stx2a phage among the Norwegian O145:H25 STEC that shared integration site with a shadow phage in all isolates. Both phage and shadow phage harbored several virulence-associated genes that may contribute to the increased pathogenicity of O145:H25 STEC.

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.

F4- and F18-Positive Enterotoxigenic Escherichia coli Isolates from Diarrhea of Postweaning Pigs: Genomic Characterization

This study aimed to characterize in silico enterotoxigenic Escherichia coli F4- and F18-positive isolates (n = 90) causing swine postweaning diarrhea, including pathogenic potential, phylogenetic relationship, antimicrobial and biocide resistance, prophage content, and metal tolerance rates. F4 strains belonged mostly to the O149 and O6 serogroups and ST100 and ST48 sequence types (STs). F18 strains were mainly assigned to the O8 and O147 serogroups and ST10, ST23, and ST42. The highest rates of antimicrobial resistance were found against streptomycin, sulfamethoxazole, tetracycline, trimethoprim, and ampicillin. No resistance was found toward ciprofloxacin, cefotaxime, ceftiofur, and colistin. Genes conferring tolerance to copper (showing the highest diversity), cadmium, silver, and zinc were predicted in all genomes. Enterotoxin genes (ltcA, 100% F4, 62% F18; astA, 100% F4, 38.1% F18; sta, 18.8% F4, 38.1% F18; stb, 100% F4, 76.2% F18) and fimbria-encoding genes typed as F4ac and F18ac were detected in all strains, in addition to up to 16 other virulence genes in individual strains. Phage analysis predicted between 7 and 20 different prophage regions in each strain. A highly diverse variety of plasmids was found; IncFII, IncFIB, and IncFIC were prevalent among F4 isolates, while IncI1 and IncX1 were dominant among F18 strains. Interestingly, F4 isolates from the early 1990s belonged to the same clonal group detected for most of the F4 strains from 2018 to 2019 (ONT:H10-A-ST100-CH27-0). The small number of single-nucleotide polymorphism differences between the oldest and recent F4 ST100 isolates suggests a relatively stable genome. Overall, the isolates analyzed in this study showed remarkably different genetic traits depending on the fimbria type.IMPORTANCE Diarrhea in the postweaning period due to enterotoxigenic E. coli (ETEC) is an economically relevant disease in pig production worldwide. In Denmark, prevention is mainly achieved by zinc oxide administration (to be discontinued by 2022). In addition, a breeding program has been implemented that aims to reduce the prevalence of this illness. Treatment with antimicrobials contributes to the problem of antimicrobial resistance (AMR) development. As a novelty, this study aims to deeply understand the genetic population structure and variation among diarrhea-associated isolates by whole-genome sequencing characterization. ST100-F4ac is the dominant clonal group circulating in Danish herds and showed high similarity to ETEC ST100 isolates from China, the United States, and Spain. High rates of AMR and high diversity of virulence genes were detected. The characterization of diarrhea-related ETEC is important for understanding the disease epidemiology and pathogenesis and for implementation of new strategies aiming to reduce the impact of the disease in pig production.

Investigation on the Evolution of Shiga Toxin-Converting Phages Based on Whole Genome Sequencing

Bacteriophages are pivotal elements in the dissemination of virulence genes. The main virulence determinants of Shiga Toxin producing E. coli, Shiga Toxins (Stx), are encoded by genes localized in the genome of lambdoid bacteriophages. Stx comprise two antigenically different types, Stx1 and Stx2, further divided into subtypes. Among these, certain Stx2 subtypes appear to be more commonly occurring in the most severe forms of the STEC disease, haemorrhagic colitis and haemolytic uremic syndrome (HUS). This study aimed at obtaining insights on the evolution of Stx2 bacteriophages, due to their relevance in public health, and we report here on the analysis of the genomic structure of Stx2 converting phages in relation with the known reservoir of the E. coli strains harboring them. Stx2-converting phages conveying the genes encoding different stx2 subtypes have been isolated from STEC strains and their whole genomes have been sequenced, analyzed and compared to those of other Stx2 phages available in the public domain. The phages' regions containing the stx2 genes have been analyzed in depth allowing to make inference on the possible mechanisms of selection and maintenance of certain Stx2 phages in the reservoir. The "stx regions" of different stx2 gene subtypes grouped into three different evolutionary lines in the comparative analysis, reflecting the frequency with which these subtypes are found in different animal niches, suggesting that the colonization of specific reservoir by STEC strains could be influenced by the Stx phage that they carry. Noteworthy, we could identify the presence of nanS-p gene exclusively in the "stx regions" of the phages identified in STEC strains commonly found in cattle. As a matter of fact, this gene encodes an esterase capable of metabolizing sialic acids produced by submaxillary glands of bovines and present in great quantities in their gastrointestinal tract.

Ea22 Proteins from Lambda and Shiga Toxin-Producing Bacteriophages Balance Structural Diversity with Functional Similarity

Enterohemorrhagic Escherichia coli (EHEC) outbreaks are commonly associated with contaminated food sources. Unlike normal intestinal bacteria, EHEC are lysogens of lambdoid bacteriophages that also carry a gene for Shiga toxin. Oxidative attack by the immune system or other stressors on the bacterial host can activate the lytic pathway of the latent phage genome to produce phage progeny and the release of Shiga toxin into the surrounding tissues. Within the genomes of bacteriophage λ and Shiga toxin-expressing (Stx⁺) phages such as φ24_B and φP27, there is a conserved set of open reading frames that is located between the exo and xis genes that influences the lysogenic-lytic decision. In this report, we have focused on the largest exo-xis region open reading frame termed ea22 that has been shown previously to have prolysogenic properties. Using a variety of biophysical and bioinformatic methods, we demonstrate that λ and φP27 Ea22 proteins are tetrameric in solution and can be considered in terms of an amino-terminal region, a central coiled-coil region, and a carboxy-terminal region. The carboxy-terminal regions of λ and φ24_B Ea22, expressed on their own, form dimers with exceptional thermostability. Limited proteolysis of φP27 Ea22 also identified a C-terminal region along the predicted boundaries. While the three Ea22 proteins all appear to have the hallmarks of a domain in their respective C-terminal regions, each sequence is remarkably dissimilar. To reconcile this difference among Ea22 proteins from λ and Stx⁺ phages alike, we speculate that each Ea22 may achieve the same function by targeting different components of the same regulatory process in the host.

Emerging Multidrug-Resistant Hybrid Pathotype Shiga Toxin-Producing Escherichia coli O80 and Related Strains of Clonal Complex 165, Europe

Enterohemorrhagic Escherichia coli serogroup O80, involved in hemolytic uremic syndrome associated with extraintestinal infections, has emerged in France. We obtained circularized sequences of the O80 strain RDEx444, responsible for hemolytic uremic syndrome with bacteremia, and noncircularized sequences of 35 O80 E. coli isolated from humans and animals in Europe with or without Shiga toxin genes. RDEx444 harbored a mosaic plasmid, pR444_A, combining extraintestinal virulence determinants and a multidrug resistance-encoding island. All strains belonged to clonal complex 165, which is distantly related to other major enterohemorrhagic E. coli lineages. All stx-positive strains contained eae-ξ, ehxA, and genes characteristic of pR444_A. Among stx-negative strains, 1 produced extended-spectrum β-lactamase, 1 harbored the colistin-resistance gene mcr1, and 2 possessed genes characteristic of enteropathogenic and pyelonephritis E. coli. Because O80-clonal complex 165 strains can integrate intestinal and extraintestinal virulence factors in combination with diverse drug-resistance genes, they constitute dangerous and versatile multidrug-resistant pathogens.

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.

Multidrug resistance and multivirulence plasmids in enterotoxigenic and hybrid Shiga toxin-producing/enterotoxigenic Escherichia coli isolated from diarrheic pigs in Switzerland

Enterovirulent Escherichia coli infections cause significant losses in the pig industry. However, information about the structures of the virulence and multidrug resistance (MDR) plasmids harboured by these strains is sparse. In this study, we used whole-genome sequencing with PacBio and Illumina platforms to analyse the molecular features of the multidrug-resistant enterotoxigenic E. coli (ETEC) strain 14OD0056 and the multidrug-resistant hybrid Shiga toxin-producing/enterotoxigenic E. coli (STEC/ETEC) strain 15OD0495 isolated from diarrheic pigs in Switzerland. Strain 14OD0056 possessed three virulence plasmids similar to others previously found in ETEC strains, while 15OD0495 harboured a 119-kb multivirulence IncFII/IncX1 hybrid STEC/ETEC plasmid (p15ODTXV) that co-carried virulence genes of both ETEC and STEC pathotypes, confirming the key role of plasmids in the emergence of hybrid pathotypes. All resistance genes of 14OD0056 that conferred resistance to ampicillin (bla_TEM-1b), gentamicin (aac(3)-IIa), kanamycin (aph(3')-Ia), sulfonamide (sul1 and sul2), streptomycin (aph(3″)-Ib, aph(6)-Id), tetracycline (tet(B)) and trimethoprim (dfrA1) were identified on a single 207-kb conjugative MDR plasmid of incompatibility group (Inc) IncHI1/IncFIA (p14ODMR). Strain 15OD0495 carried two antimicrobial resistance plasmids (p15ODAR and p15ODMR). The 99-kb IncI1 plasmid p15ODAR harboured only aminoglycoside resistance genes (aac(3)-IIa, aph(3″)-Ib, aph(6)-Id, aph(4)-Ia), whilst the 49-kb IncN MDR plasmid p15ODMR carried genes conferring resistance to ampicillin (bla_TEM-1b), sulfonamide (sul2), streptomycin (aph(6)-Id), tetracycline (tet(A)) and trimethoprim (dfrA14). Filter mating assays showed that p14ODMR, p15ODMR and p15ODAR were conjugative at room temperature and 37°C. The co-localization of multiple resistance genes on MDR conjugative plasmids such as p14ODMR and p15ODMR poses the risk of simultaneous selection of several resistance traits during empirical treatment. Thus, preventive strategies and targeted therapy following antibiotic susceptibility testing should be encouraged to avoid further dissemination of such plasmids.

Lysogenic conversion of atypical enteropathogenic Escherichia coli (aEPEC) from human, murine, and bovine origin with bacteriophage Φ3538 Δstx2::cat proves their enterohemorrhagic E. coli (EHEC) progeny

Bacteriophages play an important role in the evolution of bacterial pathogens. A phage-mediated transfer of stx-genes to atypical enteropathogenic E. coli (aEPEC) which are prevalent in different hosts, would convert them to enterohemorrhagic E. coli (EHEC). We decided to confirm this hypothesis experimentally to provide conclusive evidence that aEPEC isolated from different mammalian hosts are indeed progenitors of typical EHEC which gain the ability to produce Shiga-Toxin by lysogeny with stx-converting bacteriophages, utilizing the model phage Φ3538 Δstx₂::cat. We applied a modified in vitro plaque-assay, using a high titer of a bacteriophage carrying a deletion in the stx₂ gene (Φ3538 Δstx₂::cat) to increase the detection of lysogenic conversion events. Three wild-type aEPEC strains were chosen as acceptor strains: the murine aEPEC-strain IMT14505 (sequence type (ST)28, serotype Ont:H6), isolated from a striped field mouse (Apodemus agrarius) in the surrounding of a cattle shed, and the human aEPEC-strain 910#00 (ST28, Ont:H6). The close genomic relationship of both strains implies a high zoonotic potential. A third strain, the bovine aEPEC IMT19981, was of serotype O26:H11 and ST21 (STC29). All three aEPEC were successfully lysogenized with phage Φ3538 Δstx₂::cat. Integration of the bacteriophage DNA into the aEPEC host genomes was confirmed by amplification of chloramphenicol transferase (cat) marker gene and by Southern-Blot hybridization. Analysis of the whole genome sequence of each of the three lysogens showed that the bacteriophage was integrated into the known tRNA integration site argW, which is highly variable among E. coli. In conclusion, the successful lysogenic conversion of aEPEC with a stx-phage in vitro underlines the important role of aEPEC as progenitors of EHEC. Given the high prevalence and the wide host range of aEPEC acceptors, their high risk of zoonotic transmission should be recognized in infection control measures.

Why Bacteriophage Encode Exotoxins and other Virulence Factors

This study considers gene location within bacteria as a function of genetic element mobility. Our emphasis is on prophage encoding of bacterial virulence factors (VFs). At least four mechanisms potentially contribute to phage encoding of bacterial VFs: (i) Enhanced gene mobility could result in greater VF gene representation within bacterial populations. We question, though, why certain genes but not others might benefit from this mobility. (ii) Epistatic interactions—between VF genes and phage genes that enhance VF utility to bacteria—could maintain phage genes via selection acting on individual, VF-expressing bacteria. However, is this mechanism sufficient to maintain the rest of phage genomes or, without gene co-regulation, even genetic linkage between phage and VF genes? (iii) Phage could amplify VFs during disease progression by carrying them to otherwise commensal bacteria colocated within the same environment. However, lytic phage kill bacteria, thus requiring assumptions of inclusive fitness within bacterial populations to explain retention of phage-mediated VF amplification for the sake of bacterial utility. Finally, (iv) phage-encoded VFs could enhance phage Darwinian fitness, particularly by acting as ecosystem-modifying agents. That is, VF-supplied nutrients could enhance phage growth by increasing the density or by improving the physiology of phage-susceptible bacteria. Alternatively, VF-mediated break down of diffusion-inhibiting spatial structure found within the multicellular bodies of host organisms could augment phage dissemination to new bacteria or to environments. Such phage-fitness enhancing mechanisms could apply particularly given VF expression within microbiologically heterogeneous environments, ie, ones where phage have some reasonable potential to acquire phage-susceptible bacteria.

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.

Global occurrence and heterogeneity of the Roseobacter-clade species Ruegeria mobilis

Tropodithietic acid (TDA)-producing Ruegeria mobilis strains of the Roseobacter clade have primarily been isolated from marine aquaculture and have probiotic potential due to inhibition of fish pathogens. We hypothesized that TDA producers with additional novel features are present in the oceanic environment. We isolated 42 TDA-producing R. mobilis strains during a global marine research cruise. While highly similar on the 16S ribosomal RNA gene level (99-100% identity), the strains separated into four sub-clusters in a multilocus sequence analysis. They were further differentiated to the strain level by average nucleotide identity using pairwise genome comparison. The four sub-clusters could not be associated with a specific environmental niche, however, correlated with the pattern of sub-typing using co-isolated phages, the number of prophages in the genomes and the distribution in ocean provinces. Major genomic differences within the sub-clusters include prophages and toxin-antitoxin systems. In general, the genome of R. mobilis revealed adaptation to a particle-associated life style and querying TARA ocean data confirmed that R. mobilis is more abundant in the particle-associated fraction than in the free-living fraction occurring in 40% and 6% of the samples, respectively. Our data and the TARA data, although lacking sufficient data from the polar regions, demonstrate that R. mobilis is a globally distributed marine bacterial species found primarily in the upper open oceans. It has preserved key phenotypic behaviors such as the production of TDA, but contains diverse sub-clusters, which could provide new capabilities for utilization in aquaculture.

Postinfectious Hemolytic Uremic Syndrome

Post-infectious hemolytic uremic syndrome (HUS) is caused by specific pathogens in patients with no identifiable HUS-associated genetic mutation or autoantibody. The majority of episodes is due to infections by Shiga toxin (Stx) producing Escherichia coli (STEC). This chapter reviews the epidemiology and pathogenesis of STEC-HUS, including bacterial-derived factors and host responses. STEC disease is characterized by hematological (microangiopathic hemolytic anemia), renal (acute kidney injury) and extrarenal organ involvement. Clinicians should always strive for an etiological diagnosis through the microbiological or molecular identification of Stx-producing bacteria and Stx or, if negative, serological assays. Treatment of STEC-HUS is supportive; more investigations are needed to evaluate the efficacy of putative preventive and therapeutic measures, such as non-phage-inducing antibiotics, volume expansion and anti-complement agents. The outcome of STEC-HUS is generally favorable, but chronic kidney disease, permanent extrarenal, mainly cerebral complication and death (in less than 5 %) occur and long-term follow-up is recommended. The remainder of this chapter highlights rarer forms of (post-infectious) HUS due to S. dysenteriae, S. pneumoniae, influenza A and HIV and discusses potential interactions between these pathogens and the complement system.

Isolation, genome sequencing and functional analysis of two T7-like coliphages of avian pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC) causes colibacillosis, which results in significant economic losses to the poultry industry worldwide. Due to the drug residues and increased antibiotic resistance caused by antibiotic use, bacteriophages and other alternative therapeutic agents are expected to control APEC infection in poultry. Two APEC phages, named P483 and P694, were isolated from the feces from the farmers market in China. We then studied their biological properties, and carried out high-throughput genome sequencing and homology analyses of these phages. Assembly results of high-throughput sequencing showed that the structures of both P483 and P694 genomes consist of linear and double-stranded DNA. Results of the electron microscopy and homology analysis revealed that both P483 and P694 belong to T7-like virus which is a member of the Podoviridae family of the Caudovirales order. Comparative genomic analysis showed that most of the predicted proteins of these two phages showed strongest sequence similarity to the Enterobacteria phages BA14 and 285P, Erwinia phage FE44, and Kluyvera phage Kvp1; however, some proteins such as gp0.6a, gp1.7 and gp17 showed lower similarity (<85%) with the homologs of other phages in the T7 subgroup. We also found some unique characteristics of P483 and P694, such as the two types of the genes of P694 and no lytic activity of P694 against its host bacteria in liquid medium. Our results serve to further our understanding of phage evolution of T7-like coliphages and provide the potential application of the phages as therapeutic agents for the treatment of diseases.

Mobile effector proteins on phage genomes

Bacteriophage genomes found in a range of bacterial pathogens encode a diverse array of virulence factors ranging from superantigens or pore forming lysins to numerous exotoxins. Recent studies have uncovered an entirely new class of bacterial virulence factors, called effector proteins or effector toxins, which are encoded within phage genomes that reside among several pathovars of Escherichia coli and Salmonella enterica. These effector proteins have multiple domains resulting in proteins that can be multifunctional. The effector proteins encoded within phage genomes are translocated directly from the bacterial cytosol into their eukaryotic target cells by specialized bacterial type three secretion systems (T3SSs). In this review, we will give an overview of the different types of effector proteins encoded within phage genomes and examine their roles in bacterial pathogenesis.

Taxonomy Meets Public Health: The Case of Shiga Toxin-Producing Escherichia coli

To help assess the clinical and public health risks associated with different Shiga toxin-producing Escherichia coli (STEC) strains, an empirical classification scheme was used to classify STEC into five “seropathotypes” (seropathotype A [high risk] to seropathotypes D and E [minimal risk]). This definition is of considerable value in cases of human infection but is also problematic because not all STEC infections are fully characterized and coupled to reliable clinical information. Outbreaks with emerging hybrid strains continuously challenge our understanding of virulence potential and may result in incorrect classification of specific pathotypes; an example is the hybrid strain that caused the 2011 outbreak in Germany, STEC/EAggEC O104:H4, which may deserve an alternative seropathotype designation. The integration of mobile virulence factors in the stepwise and parallel evolution of pathogenic lineages of STEC collides with the requirements of a good taxonomy, which separates elements of each group into subgroups that are mutually exclusive, unambiguous, and, together, include all possibilities. The concept of (sero)-pathotypes is therefore challenged, and the need to identify factors of STEC that absolutely predict the potential to cause human disease is obvious. Because the definition of hemolytic-uremic syndrome (HUS) is distinct, a basic and primary definition of HUS-associated E. coli (HUSEC) for first-line public health action is proposed: stx2 in a background of an eae- or aggR -positive E. coli followed by a second-line subtyping of stx genes that refines the definition of HUSEC to include only stx2a and stx2d . All other STEC strains are considered “low-risk” STEC.

Draft Genome Sequence of an Escherichia coli O157:H43 Strain Isolated from Cattle

Here we report the draft genome sequence of an Escherichia coli O157:H43 strain, designated T22, with an atypical virulence gene profile and isolated from healthy cattle. T22 produces cytolethal distending toxin V (CDT-V) and belongs to phylogenetic group B1 and sequence type 155 (ST155).

Phylogenetic and molecular analysis of food-borne shiga toxin-producing Escherichia coli

Seventy-five food-associated Shiga toxin-producing Escherichia coli (STEC) strains were analyzed by molecular and phylogenetic methods to describe their pathogenic potential. The presence of the locus of proteolysis activity (LPA), the chromosomal pathogenicity island (PAI) PAI ICL3, and the autotransporter-encoding gene sabA was examined by PCR. Furthermore, the occupation of the chromosomal integration sites of the locus of enterocyte effacement (LEE), selC, pheU, and pheV, as well as the Stx phage integration sites yehV, yecE, wrbA, z2577, and ssrA, was analyzed. Moreover, the antibiotic resistance phenotypes of all STEC strains were determined. Multilocus sequence typing (MLST) was performed, and sequence types (STs) and sequence type complexes (STCs) were compared with those of 42 hemolytic-uremic syndrome (HUS)-associated enterohemorrhagic E. coli (HUSEC) strains. Besides 59 STs and 4 STCs, three larger clusters were defined in this strain collection. Clusters A and C consist mostly of highly pathogenic eae-positive HUSEC strains and some related food-borne STEC strains. A member of a new O26 HUS-associated clone and the 2011 outbreak strain E. coli O104:H4 were found in cluster A. Cluster B comprises only eae-negative food-borne STEC strains as well as mainly eae-negative HUSEC strains. Although food-borne strains of cluster B were not clearly associated with disease, serotypes of important pathogens, such as O91:H21 and O113:H21, were in this cluster and closely related to the food-borne strains. Clonal analysis demonstrated eight closely related genetic groups of food-borne STEC and HUSEC strains that shared the same ST and were similar in their virulence gene composition. These groups should be considered with respect to their potential for human infection.

Carriage of stx2a differentiates clinical and bovine-biased strains of Escherichia coli O157

Background

Shiga toxin (Stx) are cardinal virulence factors of enterohemorrhagic E. coli O157:H7 (EHEC O157). The gene content and genomic insertion sites of Stx-associated bacteriophages differentiate clinical genotypes of EHEC O157 (CG, typical of clinical isolates) from bovine-biased genotypes (BBG, rarely identified among clinical isolates). This project was designed to identify bacteriophage-mediated differences that may affect the virulence of CG and BBG.

Methods

Stx-associated bacteriophage differences were identified by whole genome optical scans and characterized among >400 EHEC O157 clinical and cattle isolates by PCR.

Results

Optical restriction maps of BBG strains consistently differed from those of CG strains only in the chromosomal insertion sites of Stx2-associated bacteriophages. Multiplex PCRs (stx1, stx2a, and stx2c as well as Stx-associated bacteriophage - chromosomal insertion site junctions) revealed four CG and three BBG that accounted for >90% of isolates. All BBG contained stx2c and Stx2c-associated bacteriophage – sbcB junctions. All CG contained stx2a and Stx2a-associated bacteriophage junctions in wrbA or argW.

Conclusions

Presence or absence of stx2a (or another product encoded by the Stx2a-associated bacteriophage) is a parsimonious explanation for differential virulence of BBG and CG, as reflected in the distributions of these genotypes in humans and in the cattle reservoir.

Comparative genomics and stx phage characterization of LEE-negative Shiga toxin-producing Escherichia coli

Infection by Escherichia coli and Shigella species are among the leading causes of death due to diarrheal disease in the world. Shiga toxin-producing E. coli (STEC) that do not encode the locus of enterocyte effacement (LEE-negative STEC) often possess Shiga toxin gene variants and have been isolated from humans and a variety of animal sources. In this study, we compare the genomes of nine LEE-negative STEC harboring various stx alleles with four complete reference LEE-positive STEC isolates. Compared to a representative collection of prototype E. coli and Shigella isolates representing each of the pathotypes, the whole genome phylogeny demonstrated that these isolates are diverse. Whole genome comparative analysis of the 13 genomes revealed that in addition to the absence of the LEE pathogenicity island, phage-encoded genes including non-LEE encoded effectors, were absent from all nine LEE-negative STEC genomes. Several plasmid-encoded virulence factors reportedly identified in LEE-negative STEC isolates were identified in only a subset of the nine LEE-negative isolates further confirming the diversity of this group. In combination with whole genome analysis, we characterized the lambdoid phages harboring the various stx alleles and determined their genomic insertion sites. Although the integrase gene sequence corresponded with genomic location, it was not correlated with stx variant, further highlighting the mosaic nature of these phages. The transcription of these phages in different genomic backgrounds was examined. Expression of the Shiga toxin genes, stx(1) and/or stx(2), as well as the Q genes, were examined with quantitative reverse transcriptase polymerase chain reaction assays. A wide range of basal and induced toxin induction was observed. Overall, this is a first significant foray into the genome space of this unexplored group of emerging and divergent pathogens.

Development and application of a method for the purification of free shigatoxigenic bacteriophage from environmental samples

Shiga toxin producing Escherichia coli (STEC) strains are foodborne pathogens whose ability to produce Shiga toxin (Stx) is due to the integration of Stx-encoding lambdoid bacteriophage (Stx phage). Circulating, infective Stx phages are very difficult to isolate, purify and propagate such that there is no information on their genetic composition and properties. Here we describe a novel approach that exploits the phage's ability to infect their host and form a lysogen, thus enabling purification of Stx phages by a series of sequential lysogen isolation and induction steps. A total of 15 Stx phages were rigorously purified from water samples in this way, classified by TEM and genotyped using a PCR-based multi-loci characterisation system. Each phage possessed only one variant of each target gene type, thus confirming its purity, with 9 of the 15 phages possessing a short tail-spike gene and identified by TEM as Podoviridae. The remaining 6 phages possessed long tails, four of which appeared to be contractile in nature (Myoviridae) and two of which were morphologically very similar to bacteriophage lambda (Siphoviridae).

Comparative genomics of Shiga toxin encoding bacteriophages

Background

Stx bacteriophages are responsible for driving the dissemination of Stx toxin genes (stx) across their bacterial host range. Lysogens carrying Stx phages can cause severe, life-threatening disease and Stx toxin is an integral virulence factor. The Stx-bacteriophage vB_EcoP-24_B, commonly referred to as Ф24_B, is capable of multiply infecting a single bacterial host cell at a high frequency, with secondary infection increasing the rate at which subsequent bacteriophage infections can occur. This is biologically unusual, therefore determining the genomic content and context of Ф24_B compared to other lambdoid Stx phages is important to understanding the factors controlling this phenomenon and determining whether they occur in other Stx phages.

Results

The genome of the Stx2 encoding phage, Ф24_B was sequenced and annotated. The genomic organisation and general features are similar to other sequenced Stx bacteriophages induced from Enterohaemorrhagic Escherichia coli (EHEC), however Ф24_B possesses significant regions of heterogeneity, with implications for phage biology and behaviour. The Ф24_B genome was compared to other sequenced Stx phages and the archetypal lambdoid phage, lambda, using the Circos genome comparison tool and a PCR-based multi-loci comparison system.

Conclusions

The data support the hypothesis that Stx phages are mosaic, and recombination events between the host, phages and their remnants within the same infected bacterial cell will continue to drive the evolution of Stx phage variants and the subsequent dissemination of shigatoxigenic potential.

Escherichia coli serotype O55:H7 diversity supports parallel acquisition of bacteriophage at Shiga toxin phage insertion sites during evolution of the O157:H7 lineage

Enteropathogenic Escherichia coli (EPEC) continues to be a leading cause of mortality and morbidity in children around the world. Two EPEC genomes have been fully sequenced: those of EPEC O127:H6 strain E2348/69 (United Kingdom, 1969) and EPEC O55:H7 strain CB9615 (Germany, 2003). The O55:H7 serotype is a recent precursor to the virulent enterohemorrhagic E. coli O157:H7. To explore the diversity of O55:H7 and better understand the clonal evolution of O157:H7, we fully sequenced EPEC O55:H7 strain RM12579 (California, 1974), which was collected 1 year before the first U.S. isolate of O157:H7 was identified in California. Phage-related sequences accounted for nearly all differences between the two O55:H7 strains. Additionally, O55:H7 and O157:H7 strains were tested for the presence and insertion sites of Shiga toxin gene (stx)-containing bacteriophages. Analysis of non-phage-associated genes supported core elements of previous O157:H7 stepwise evolutionary models, whereas phage composition and insertion analyses suggested a key refinement. Specifically, the placement and presence of lambda-like bacteriophages (including those containing stx) should not be considered stable evolutionary markers or be required in placing O55:H7 and O157:H7 strains within the stepwise evolutionary models. Additionally, we suggest that a 10.9-kb region (block 172) previously believed unique to O55:H7 strains can be used to identify early O157:H7 strains. Finally, we defined two subsets of O55:H7 strains that share an as-yet-unobserved or extinct common ancestor with O157:H7 strains. Exploration of O55:H7 diversity improved our understanding of the evolution of E. coli O157:H7 and suggested a key revision to accommodate existing and future configurations of stx-containing bacteriophages into current models.

Reticulate classification of mosaic microbial genomes using NeAT website

The tree of life is the classical representation of the evolutionary relationships between existent species. A tree is appropriate to display the divergence of species through mutation, i.e., by vertical descent. However, lateral gene transfer (LGT) is excluded from such representations. When LGT contribution to genome evolution cannot be neglected (e.g., for prokaryotes and mobile genetic elements), the tree becomes misleading. Networks appear as an intuitive way to represent both vertical and horizontal relationships, while overlapping groups within such graphs are more suitable for their classification. Here, we describe a method to represent both vertical and horizontal relationships. We start with a set of genomes whose coded proteins have been grouped into families based on sequence similarity. Next, all pairs of genomes are compared, counting the number of proteins classified into the same family. From this comparison, we derive a weighted graph where genomes with a significant number of similar proteins are linked. Finally, we apply a two-step clustering of this graph to produce a classification where nodes can be assigned to multiple clusters. The procedure can be performed using the Network Analysis Tools (NeAT) website.

Differential virulence of clinical and bovine-biased enterohemorrhagic Escherichia coli O157:H7 genotypes in piglet and Dutch belted rabbit models

Enterohemorrhagic Escherichia coli O157:H7 (EHEC O157) is an important cause of food and waterborne illness in the developed countries. Cattle are a reservoir host of EHEC O157 and a major source of human exposure through contaminated meat products. Shiga toxins (Stxs) are an important pathogenicity trait of EHEC O157. The insertion sites of the Stx-encoding bacteriophages differentiate EHEC O157 isolates into genogroups commonly isolated from cattle but rarely from sick humans (bovine-biased genotypes [BBG]) and those commonly isolated from both cattle and human patients (clinical genotypes [CG]). Since BBG and CG share the cardinal virulence factors of EHEC O157 and are carried by cattle at similar prevalences, the infrequent occurrence of BBG among human disease isolates suggests that they may be less virulent than CG. We compared the virulence potentials of human and bovine isolates of CG and BBG in newborn conventional pig and weaned Dutch Belted rabbit models. CG-challenged piglets experienced severe disease accompanied by early and high mortality compared to BBG-challenged piglets. Similarly, CG-challenged rabbits were likely to develop lesions in kidney and intestine compared with the BBG-challenged rabbits. The CG strains used in this study carried stx2 and produced significantly higher amounts of Stx, whereas the BBG strains carried the stx2c gene variant only. These results suggest that BBG are less virulent than CG and that this difference in virulence potential is associated with the Stx2 subtype(s) carried and/or the amount of Stx produced.

Genetic background and mobility of variants of the gene nleA in attaching and effacing Escherichia coli

In this study, we characterized the genetic background of various nleA variants in 106 Shiga toxin-producing Escherichia coli (STEC) and enteropathogenic Escherichia coli (EPEC) strains. The flanking regions of eight nleA variants were analyzed by DNA sequencing and compared with the corresponding regions of five previously described NleA-encoding prophages. The analyzed nleA variants were all located downstream of the DNA region responsible for phage morphogenesis. In particular, the type III effector genes avrA, ospB, nleH, and nleG and IS elements were detected in the neighborhood of nleA. The structure of the eight analyzed regions flanking nleA primarily resembled the corresponding region of the NleA₄₇₉₅-encoding prophage BP-4795. Using PCR, the gene order flanking 13 nleA variants in strains of different serogroups was compared to the respective regions in reference strains. The analyses showed that strains which harbor prophages with conserved flanking regions of a particular nleA variant predominantly occurred, and IS elements were additionally detected in these regions. We were able to mobilize nleA by transduction in 20% of strains determined, which comprised in particular EPEC strains harboring an nleA variant, the gene encoding the protein known as "EspI-like." Plaque hybridization was used to identify phages that harbor the genes stx and nleA. However, only two strains harbored variant nleA₄₇₉₅ in the genome of an Stx1 prophage.

The use of genomic signature distance between bacteriophages and their hosts displays evolutionary relationships and phage growth cycle determination

BACKGROUND

Bacteriophage classification is mainly based on morphological traits and genome characteristics combined with host information and in some cases on phage growth lifestyle. A lack of molecular tools can impede more precise studies on phylogenetic relationships or even a taxonomic classification. The use of methods to analyze genome sequences without the requirement for homology has allowed advances in classification.

RESULTS

Here, we proposed to use genome sequence signature to characterize bacteriophages and to compare them to their host genome signature in order to obtain host-phage relationships and information on their lifestyle. We analyze the host-phage relationships in the four most representative groups of Caudoviridae, the dsDNA group of phages. We demonstrate that the use of phage genomic signature and its comparison with that of the host allows a grouping of phages and is also able to predict the host-phage relationships (lytic vs. temperate).

CONCLUSIONS

We can thus condense, in relatively simple figures, this phage information dispersed over many publications.

Non-O157 verotoxigenic Escherichia coli and beef: a Canadian perspective

Verotoxigenic Escherichia coli (VTEC) are important foodborne pathogens in Canada. VTEC of the O157:H7 serogroup have been the focus of regulatory action and surveillance in both Canada and the USA, due to their role in a number of high profile outbreaks. However, there is increasing evidence that other VTEC serogroups cause a substantial proportion of human illness. This issue is of particular importance to the cattle industry due to the role of beef as a vehicle for VTEC transmission. In this review, the evidence for non-O157 VTEC as cause of human illness in Canada and the potential for Canadian beef and cattle to serve as a source of VTEC are presented. In addition, the available strategies for the control of VTEC in cattle and beef are discussed.

Comparative analyses of prophage-like elements present in bifidobacterial genomes

So far, only scarce and rather diffuse information is available on bacteriophages infecting members of the genus Bifidobacterium. In the current study, we investigated the genetic organization, phylogenetic relationships, and, in some cases, transcription profiles and inducibility of 19 prophage-like elements present on the individual chromosomes of nine bifidobacterial strains, which represent six different species.

Genomic diversity of pathogenic Escherichia coli of the EHEC 2 clonal complex

BACKGROUND

Evolutionary analyses of enterohemorrhagic Escherichia coli (EHEC) have identified two distantly related clonal groups: EHEC 1, including serotype O157:H7 and its inferred ancestor O55:H7; and EHEC 2, comprised of several serogroups (O26, O111, O118, etc.). These two clonal groups differ in their virulence and global distribution. Although several fully annotated genomic sequences exist for strains of serotype O157:H7, much less is known about the genomic composition of EHEC 2. In this study, we analyzed a set of 24 clinical EHEC 2 strains representing serotypes O26:H11, O111:H8/H11, O118:H16, O153:H11 and O15:H11 from humans and animals by comparative genomic hybridization (CGH) on an oligoarray based on the O157:H7 Sakai genome.

RESULTS

Backbone genes, defined as genes shared by Sakai and K-12, were highly conserved in EHEC 2. The proportion of Sakai phage genes in EHEC 2 was substantially greater than that of Sakai-specific bacterial (non-phage) genes. This proportion was inverted in O55:H7, reiterating that a subset of Sakai bacterial genes is specific to EHEC 1. Split decomposition analysis of gene content revealed that O111:H8 was more genetically uniform and distinct from other EHEC 2 strains, with respect to the Sakai O157:H7 gene distribution. Serotype O26:H11 was the most heterogeneous EHEC 2 subpopulation, comprised of strains with the highest as well as the lowest levels of Sakai gene content conservation. Of the 979 parsimoniously informative genes, 15% were found to be compatible and their distribution in EHEC 2 clustered O111:H8 and O118:H16 strains by serotype. CGH data suggested divergence of the LEE island from the LEE1 to the LEE4 operon, and also between animal and human isolates irrespective of serotype. No correlation was found between gene contents and geographic locations of EHEC 2 strains.

CONCLUSION

The gene content variation of phage-related genes in EHEC 2 strains supports the hypothesis that extensive modular shuffling of mobile DNA elements has occurred among EHEC strains. These results suggest that EHEC 2 is a multiform pathogenic clonal complex, characterized by substantial intra-serotype genetic variation. The heterogeneous distribution of mobile elements has impacted the diversification of O26:H11 more than other EHEC 2 serotypes.

Diverse phage-encoded toxins in a protective insect endosymbiont

The lysogenic bacteriophage APSE infects "Candidatus Hamiltonella defensa," a facultative endosymbiont of aphids and other sap-feeding insects. This endosymbiont has established a beneficial association with aphids, increasing survivorship following attack by parasitoid wasps. Although APSE and "Ca. Hamiltonella defensa" are effectively maternally transmitted between aphid generations, they can also be horizontally transferred among insect hosts, which results in genetically distinct "Ca. Hamiltonella defensa" strains infecting the same aphid species and sporadic distributions of both APSE and "Ca. Hamiltonella defensa" among hosts. Aphids infected only with "Ca. Hamiltonella defensa" have significantly less protection than those infected with both "Ca. Hamiltonella defensa" and APSE. This protection has been proposed to be connected to eukaryote-targeted toxins previously discovered in the genomes of two characterized APSE strains. In this study, we have sequenced partial genomes from seven additional APSE strains to address the evolution and extent of toxin variation in this phage. The APSE lysis region has been a hot spot for nonhomologous recombination of novel virulence cassettes. We identified four new toxins from three protein families, Shiga-like toxin, cytolethal distending toxin, and YD-repeat toxins. These recombination events have also resulted in reassortment of the downstream lysozyme and holin genes. Analysis of the conserved APSE genes flanking the variable toxin cassettes reveals a close phylogenetic association with phage sequences from two other facultative endosymbionts of insects. Thus, phage may act as a conduit for ongoing gene exchange among heritable endosymbionts.

Evaluation of major types of Shiga toxin 2E-producing Escherichia coli bacteria present in food, pigs, and the environment as potential pathogens for humans

Shiga toxin 2e (Stx2e)-producing strains from food (n = 36), slaughtered pigs (n = 25), the environment (n = 21), diseased pigs (n = 19), and humans (n = 9) were investigated for production of Stx2e by enzyme-linked immunosorbent assay, for virulence markers by PCR, and for their serotypes to evaluate their role as potential human pathogens. Stx2e production was low in 64% of all 110 strains. Stx2e production was inducible by mitomycin C but differed considerably between strains. Analysis by nucleotide sequencing and transcription of stx(2e) genes in high- and low-Stx2e-producing strains showed that toxin production correlated with transcription rates of stx(2e) genes. DNA sequences specific for the int, Q, dam, and S genes of the stx(2e) bacteriophage P27 were found in 109 strains, indicating cryptic P27-like prophages, although 102 of these were not complete for all genes tested. Genes encoding intimin (eae), enterohemorrhagic Escherichia coli hemolysin (ehx), or other stx(1) or stx(2) variants were not found, whereas genes for heat-stable enterotoxins STI, STII, or EAST1 were present in 54.5% of the strains. Seven major serotypes that were associated with diseased pigs (O138:H14, O139:H1, and O141:H4) or with slaughter pigs, food, and the environment (O8:H4, O8:H9, O100:H30, and O101:H9) accounted for 60% of all Stx2e strains. The human Stx2e isolates did not belong to these major serotypes of Stx2e strains, and high production of Stx2e in human strains was not related to diarrheal disease. The results from this study and other studies do not point to Stx2e as a pathogenicity factor for diarrhea and hemolytic uremic syndrome in humans.

Genome classification by gene distribution: an overlapping subspace clustering approach

BACKGROUND

Genomes of lower organisms have been observed with a large amount of horizontal gene transfers, which cause difficulties in their evolutionary study. Bacteriophage genomes are a typical example. One recent approach that addresses this problem is the unsupervised clustering of genomes based on gene order and genome position, which helps to reveal species relationships that may not be apparent from traditional phylogenetic methods.

RESULTS

We propose the use of an overlapping subspace clustering algorithm for such genome classification problems. The advantage of subspace clustering over traditional clustering is that it can associate clusters with gene arrangement patterns, preserving genomic information in the clusters produced. Additionally, overlapping capability is desirable for the discovery of multiple conserved patterns within a single genome, such as those acquired from different species via horizontal gene transfers. The proposed method involves a novel strategy to vectorize genomes based on their gene distribution. A number of existing subspace clustering and biclustering algorithms were evaluated to identify the best framework upon which to develop our algorithm; we extended a generic subspace clustering algorithm called HARP to incorporate overlapping capability. The proposed algorithm was assessed and applied on bacteriophage genomes. The phage grouping results are consistent overall with the Phage Proteomic Tree and showed common genomic characteristics among the TP901-like, Sfi21-like and sk1-like phage groups. Among 441 phage genomes, we identified four significantly conserved distribution patterns structured by the terminase, portal, integrase, holin and lysin genes. We also observed a subgroup of Sfi21-like phages comprising a distinctive divergent genome organization and identified nine new phage members to the Sfi21-like genus: Staphylococcus 71, phiPVL108, Listeria A118, 2389, Lactobacillus phi AT3, A2, Clostridium phi3626, Geobacillus GBSV1, and Listeria monocytogenes PSA.

CONCLUSION

The method described in this paper can assist evolutionary study through objectively classifying genomes based on their resemblance in gene order, gene content and gene positions. The method is suitable for application to genomes with high genetic exchange and various conserved gene arrangement, as demonstrated through our application on phages.

Stx-phages: drivers and mediators of the evolution of STEC and STEC-like pathogens

Bacteriophages, also known as phages, are viruses that infect bacteria. Until recently they have been ignored by most of the scientific community, but their impact upon our world is enormous. They are the most abundant lifeform on the globe and drive the diversity and abundance of bacteria around us, including, in many instances, the pathogenic profiles of many of mankind's most feared bacterial pathogens. This article focuses on how a group of bacteriophages, Stx-phages, which carry the genes encoding Shiga toxin, have driven and are driving the emergence of Shiga toxin-producing pathogens such as the infamous Escherichia coli O157:H7. Since the emergence of this foodborne pathogen as a cause of significant human disease in 1982, more than 500 different serogroups of E. coli have been reported to produce Shiga toxin, as well as a few other organisms. These events and many more are all controlled by the biology of Stx-phages.

Multilocus characterization scheme for shiga toxin-encoding bacteriophages

Shiga toxin-producing Escherichia coli (STEC) strains are food-borne pathogens whose ability to produce Shiga toxin (Stx) is due to integration of Stx-encoding lambdoid bacteriophages. These Stx phages are both genetically and morphologically heterogeneous, and here we report the design and validation of a PCR-based multilocus typing scheme. PCR primer sets were designed for database variants of a range of key lambdoid bacteriophage genes and applied to control phages and 70 stx(+) phage preparations induced from a collection of STEC isolates. The genetic diversity residing within these populations could be described, and observations were made on the heterogeneity of individual gene targets, including the unexpected predominance of short-tailed phages with a highly conserved tail spike protein gene. Purified Stx phages can be profiled using this scheme, and the lambdoid phage-borne genes in induced STEC preparations can be identified as well as those residing in the noninducible prophage complement. The ultimate goal is to enable robust and realistically applicable epidemiological studies of Stx phages and their traits. The impact of Stx phage on STEC epidemiology is currently unknown.

Short-tailed stx phages exploit the conserved YaeT protein to disseminate Shiga toxin genes among enterobacteria

Infection of Escherichia coli by Shiga toxin-encoding bacteriophages (Stx phages) was the pivotal event in the evolution of the deadly Shiga toxin-encoding E. coli (STEC), of which serotype O157:H7 is the most notorious. The number of different bacterial species and strains reported to produce Shiga toxin is now more than 500, since the first reported STEC infection outbreak in 1982. Clearly, Stx phages are spreading rapidly, but the underlying mechanism for this dissemination has not been explained. Here we show that an essential and highly conserved gene product, YaeT, which has an essential role in the insertion of proteins in the gram-negative bacterial outer membrane, is the surface molecule recognized by the majority (ca. 70%) of Stx phages via conserved tail spike proteins associated with a short-tailed morphology. The yaeT gene was initially identified through complementation, and its role was confirmed in phage binding assays with and without anti-YaeT antiserum. Heterologous cloning of E. coli yaeT to enable Stx phage adsorption to Erwinia carotovora and the phage adsorption patterns of bacterial species possessing natural yaeT variants further supported this conclusion. The use of an essential and highly conserved protein by the majority of Stx phages is a strategy that has enabled and promoted the rapid spread of shigatoxigenic potential throughout multiple E. coli serogroups and related bacterial species. Infection of commensal bacteria in the mammalian gut has been shown to amplify Shiga toxin production in vivo, and the data from this study provide a platform for the development of a therapeutic strategy to limit this YaeT-mediated infection of the commensal flora.

Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes

BACKGROUND

Enterohemorrhagic Escherichia coli (EHEC) O157 causes severe food-borne illness in humans. The chromosome of O157 consists of 4.1 Mb backbone sequences shared by benign E. coli K-12, and 1.4 Mb O157-specific sequences encoding many virulence determinants, such as Shiga toxin genes (stx genes) and the locus of enterocyte effacement (LEE). Non-O157 EHECs belonging to distinct clonal lineages from O157 also cause similar illness in humans. According to the 'parallel' evolution model, they have independently acquired the major virulence determinants, the stx genes and LEE. However, the genomic differences between O157 and non-O157 EHECs have not yet been systematically analyzed.

RESULTS

Using microarray and whole genome PCR scanning analyses, we performed a whole genome comparison of 20 EHEC strains of O26, O111, and O103 serotypes with O157. In non-O157 EHEC strains, although genome sizes were similar with or rather larger than O157 and the backbone regions were well conserved, O157-specific regions were very poorly conserved. Around only 20% of the O157-specific genes were fully conserved in each non-O157 serotype. However, the non-O157 EHECs contained a significant number of virulence genes that are found on prophages and plasmids in O157, and also multiple prophages similar to, but significantly divergent from, those in O157.

CONCLUSION

Although O157 and non-O157 EHECs have independently acquired a huge amount of serotype- or strain-specific genes by lateral gene transfer, they share an unexpectedly large number of virulence genes. Independent infections of similar but distinct bacteriophages carrying these virulence determinants are deeply involved in the evolution of O157 and non-O157 EHECs.

Distribution, functional expression, and genetic organization of Cif, a phage-encoded type III-secreted effector from enteropathogenic and enterohemorrhagic Escherichia coli

Enteropathogenic Escherichia coli (EPEC) and enterohemorrhagic E. coli (EHEC) inject effector proteins into host cells via a type III secretion system encoded by the locus of enterocyte effacement (LEE). One of these effectors is Cif, encoded outside the LEE by a lambdoid prophage. In this study, we demonstrated that the Cif-encoding prophage of EPEC strain E22 is inducible and produces infectious phage particles. We investigated the distribution and functional expression of Cif in 5,049 E. coli strains of human, animal, and environmental origins. A total of 115 E. coli isolates from diverse origins and geographic locations carried cif. The presence of cif was tightly associated with the LEE, since all the cif-positive isolates were positive for the LEE. These results suggested that the Cif-encoding prophages have been widely disseminated within the natural population of E. coli but positively selected within the population of LEE-positive strains. Nonetheless, 66% of cif-positive E. coli strains did not induce a typical Cif-related phenotype in eukaryotic cells due to frameshift mutations or insertion of an IS element in the cif gene. The passenger region of the prophages carrying cif was highly variable and showed various combinations of IS elements and genes coding for other effectors such as nleB, nleC, nleH, nleG, espJ, and nleA/espI (some of which were also truncated). This diversity and the presence of nonfunctional effectors should be taken into account to assess EPEC and EHEC pathogenicity and tropism.

Curated list of prokaryote viruses with fully sequenced genomes

Genome sequencing is of enormous importance for classification of prokaryote viruses and for understanding the evolution of these viruses. This survey covers 284 sequenced viruses for which a full description has been published and for which the morphology is known. This corresponds to 219 (4%) of tailed and 75 (36%) of tailless viruses of prokaryotes. The number of sequenced tailless viruses almost doubles if viruses of unknown morphology are counted. The sequences are from representatives of 15 virus families and three groups without family status, including eight taxa of archaeal viruses. Tailed phages, especially those with large genomes and hosts other than enterobacteria or lactococci, mycobacteria and pseudomonads, are vastly under investigated.

Insertion site occupancy by stx2 bacteriophages depends on the locus availability of the host strain chromosome

Shiga toxin-producing Escherichia coli (STEC) is an emergent pathogen characterized by the expression of Shiga toxins, which are encoded in the genomes of lambdoid phages. These phages are infectious for other members of the Enterobacteriaceae and establish lysogeny when they integrate into the host chromosome. Five insertion sites, used mainly by these prophages, have been described to date. In the present study, the insertion of stx(2) prophages in these sites was analyzed in 168 STEC strains isolated from cattle. Additionally, insertion sites were determined for stx(2) phages which (i) converted diverse laboratory host strains, (ii) coexisted with another stx(2) prophage, and (iii) infected a recombinant host strain lacking the most commonly used insertion site. Results show that depending on the host strain, phages preferentially use one insertion site. For the most part, yehV is occupied in STEC strains while wrbA is preferentially selected by the same stx phages in E. coli laboratory strains. If this primary insertion site is unavailable, then a secondary insertion site is selected. It can be concluded that insertion site occupancy by stx phages depends on the host strain and on the availability of the preferred locus in the host strain.

Shiga toxin gene loss and transfer in vitro and in vivo during enterohemorrhagic Escherichia coli O26 infection in humans

Escherichia coli serogroup O26 consists of enterohemorrhagic E. coli (EHEC) and atypical enteropathogenic E. coli (aEPEC). The former produces Shiga toxins (Stx), major determinants of EHEC pathogenicity, encoded by bacteriophages; the latter is Stx negative. We have isolated EHEC O26 from patient stools early in illness and aEPEC O26 from stools later in illness, and vice versa. Intrapatient EHEC and aEPEC isolates had quite similar pulsed-field gel electrophoresis (PFGE) patterns, suggesting that they might have arisen by conversion between the EHEC and aEPEC pathotypes during infection. To test this hypothesis, we asked whether EHEC O26 can lose stx genes and whether aEPEC O26 can be lysogenized with Stx-encoding phages from EHEC O26 in vitro. The stx2 loss associated with the loss of Stx2-encoding phages occurred in 10% to 14% of colonies tested. Conversely, Stx2- and, to a lesser extent, Stx1-encoding bacteriophages from EHEC O26 lysogenized aEPEC O26 isolates, converting them to EHEC strains. In the lysogens and EHEC O26 donors, Stx2-converting bacteriophages integrated in yecE or wrbA. The loss and gain of Stx-converting bacteriophages diversifies PFGE patterns; this parallels findings of similar but not identical PFGE patterns in the intrapatient EHEC and aEPEC O26 isolates. EHEC O26 and aEPEC O26 thus exist as a dynamic system whose members undergo ephemeral interconversions via loss and gain of Stx-encoding phages to yield different pathotypes. The suggested occurrence of this process in the human intestine has diagnostic, clinical, epidemiological, and evolutionary implications.

Why bacteriophage encode exotoxins and other virulence factors

This study considers gene location within bacteria as a function of genetic element mobility. Our emphasis is on prophage encoding of bacterial virulence factors (VFs). At least four mechanisms potentially contribute to phage encoding of bacterial VFs: (i) Enhanced gene mobility could result in greater VF gene representation within bacterial populations. We question, though, why certain genes but not others might benefit from this mobility. (ii) Epistatic interactions-between VF genes and phage genes that enhance VF utility to bacteria-could maintain phage genes via selection acting on individual, VF-expressing bacteria. However, is this mechanism sufficient to maintain the rest of phage genomes or, without gene co-regulation, even genetic linkage between phage and VF genes? (iii) Phage could amplify VFs during disease progression by carrying them to otherwise commensal bacteria colocated within the same environment. However, lytic phage kill bacteria, thus requiring assumptions of inclusive fitness within bacterial populations to explain retention of phage-mediated VF amplification for the sake of bacterial utility. Finally, (iv) phage-encoded VFs could enhance phage Darwinian fitness, particularly by acting as ecosystem-modifying agents. That is, VF-supplied nutrients could enhance phage growth by increasing the density or by improving the physiology of phage-susceptible bacteria. Alternatively, VF-mediated break down of diffusion-inhibiting spatial structure found within the multicellular bodies of host organisms could augment phage dissemination to new bacteria or to environments. Such phage-fitness enhancing mechanisms could apply particularly given VF expression within microbiologically heterogeneous environments, ie, ones where phage have some reasonable potential to acquire phage-susceptible bacteria.