Diversity of phage integrases in Enterobacteriaceae: development of markers for environmental analysis of temperate phages

Multiomic spatial analysis reveals a distinct mucosa-associated virome

The human gut virome has been increasingly explored in recent years. However, nearly all virome-sequencing efforts rely solely on fecal samples and few studies leverage multiomic approaches to investigate phage-host relationships. Here, we combine metagenomics, metaviromics, and metatranscriptomics to study virome-bacteriome interactions at the colonic mucosal-luminal interface in a cohort of three individuals with inflammatory bowel disease; non-IBD controls were not included in this study. We show that the mucosal viral population is distinct from the stool virome and houses abundant crAss-like phages that are undetectable by fecal sampling. Through viral protein prediction and metatranscriptomic analysis, we explore viral gene transcription, prophage activation, and the relationship between the presence of integrase and temperate phages in IBD subjects. We also show the impact of deep sequencing on virus recovery and offer guidelines for selecting optimal sequencing depths in future metaviromic studies. Systems biology approaches such as those presented in this report will enhance our understanding of the human virome and its interactions with our microbiome and our health.

Impact of Shiga-toxin encoding gene transduction from O80:H2 Shiga toxigenic Escherichia coli (STEC) on non-STEC strains

Shiga toxin-producing Escherichia coli (STEC) are major foodborne pathogens that cause human diseases ranging from diarrhea to life-threatening complications including hemolytic-uremic syndrome. Virulence of STEC strains and their ability to cause severe diseases are associated with the activity of prophage-encoded Shiga toxins (Stxs). The first objective of this work was to isolate and characterize the Stx2d phage from STEC O80:H2 and to study the transfer of this phage in non-STEC strains. The second objective was to assess the survival of Galleria mellonella larvae inoculated with these transduced strains. Firstly, one bacteriophage isolated from a STEC O80:H2 strain was used to infect six non-STEC strains, resulting in the conversion of three strains. Then, stability assays were performed, showing that this phage was stable in the new STEC strains after three successive subculturing steps, as confirmed by a combination of short and long read genome sequencing approaches. This phage, vB_EcoS_ULI-O80_Stx2d, is resistant to moderate temperature and pH. It belongs to a currently unclassified genus and family within the Caudoviricetes class, shares 98% identity with Stx2_112808 phage and encodes several proteins involved in the lysogenic cycle. The yecE gene was identified at the insertion site. Finally, G. mellonella experiments showed that the transduced strains caused significantly higher mortality rates than the corresponding non-STEC strains. In conclusion, this study showed that stx2d gene from O80:H2 E. coli can be transferred to non-STEC strains and contributes to their virulence.

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.

The Site-Specific Recombination System of the Escherichia coli Bacteriophage Φ24B

Stx bacteriophages are members of the lambdoid group of phages and are responsible for Shiga toxin (Stx) production and the dissemination of Shiga toxin genes (stx) across shigatoxigenic Escherichia coli (STEC). These toxigenic bacteriophage hosts can cause life-threatening illnesses, and Stx is the virulence determinant responsible for the severe nature of infection with enterohemorrhagic E. coli, a subset of pathogenic STEC. Stx phages are temperate, and in the present study, the identification of what is actually required for Stx phage Φ24_B and bacterial DNA recombination was tested using both in vitro and in situ recombination assays. It is well established that phage λ, which underpins most of what we understand about lambdoid phage biology, requires its own encoded phage attachment site (attP) of 250 bp, a host-encoded attachment site (attB) of 21 bp, and a host-encoded DNA binding protein known as integration host factor (IHF). The assays applied in this study enabled the manipulation of the phage attachment site (attP) and the bacterial attachment site (attB) sequences and the inclusion or exclusion of a host-encoded accessory element known as integration host factor. We were able to demonstrate that the minimal attP sequence required by Φ24_B phage is between 350 and 427 bp. Unlike phage λ, the minimal necessary flanking sequences for the attB site do not appear to be equal in size, with a total length between 62 and 93 bp. Furthermore, we identified that the Φ24_B integrase does not require IHF to drive the integration and the recombination process. Understanding how this unusual Stx phage integrase works may enable exploitation of its promiscuous nature in the context of genetic engineering.

Prediction, Diversity, and Genomic Analysis of Temperate Phages Induced From Shiga Toxin-Producing Escherichia coli Strains

Shiga toxin-producing Escherichia coli (STEC) is a notorious foodborne pathogen containing stx genes located in the sequence region of Shiga toxin (Stx) prophages. Stx prophages, as one of the mobile elements, are involved in the transfer of virulence genes to other strains. However, little is known about the diversity of prophages among STEC strains. The objectives of this study were to predict various prophages from different STEC genomes and to evaluate the effect of different stress factors on Stx prophage induction. Forty bacterial whole-genome sequences of STEC strains obtained from National Center for Biotechnology Information (NCBI) were used for the prophage prediction using PHASTER webserver. Eight of the STEC strains from different serotypes were subsequently selected to quantify the induction of Stx prophages by various treatments, including antibiotics, temperature, irradiation, and antimicrobial agents. After induction, Stx1-converting phage Lys8385Vzw and Stx2-converting phage Lys12581Vzw were isolated and further confirmed for the presence of stx genes using conventional PCR. Phage morphology was observed by transmission electron microscopy. The prediction results showed an average of 8–22 prophages, with one or more encoding stx, were predicted from each STEC genome obtained in this study. Additionally, the phylogenetic analysis revealed high genetic diversity of Stx prophages among the 40 STEC genomes. However, the sequences of Stx prophages in the genomes of STEC O45, O111, and O121 strains, in general, shared higher genetic homology than those in other serotypes. Interestingly, most STEC strains with two or more stx genes carried at least one each of Stx1 and Stx2 prophages. The induction results indicated EDTA and UV were the most effective inducers of Stx1 and Stx2 prophages of the 8 selected STECs, respectively. Additionally, both Stx-converting phages could infect non-pathogenic E. coli (WG5, DH5α, and MG1655) and form new lysogens. The findings of this study confirm that Stx prophages can be induced by environmental stress, such as exposure to solar radiation, and lysogenize other commensal E. coli strains.

Development and evaluation of taxon-specific primers for the selected Caudovirales taxa

The phage taxonomy is primarily based on the morphology derived from Transmission Electron Microscopic (TEM) studies. TEM based characterization is authentic and accepted by scientific community. However, TEM based identification is expensive and time consuming. After the phage isolation, before analysis TEM, a DNA based rapid method could be introduced. The DNA based method could dramatically reduce the number of samples analyzed by TEM and thereby increase the speed and reduce the cost of identification. In the present work, four environmental phage isolates were identified based on TEM studies and genome size. The identification of these four phages was validated using DNA based method. The taxon-specific DNA markers were identified through multiple sequence alignments. The primers were designed at conserved genes (DNA polymerase or integrase) of 4 different phage taxa viz. family Ackermannviridae, genus Jerseyvirus, genus T4virus, and genus P22virus. These primers were evaluated using both in vitro and in silico approach for the amplification of the target taxons. Majority of the primer sets were found to amplify member species of the targeted taxa in vitro. In In silico analysis, six primer sets intended for identification of family Ackermannviridae showed positive amplification of ≥86.7% classified species. Further, the primers targeting the genus Jerseyvirus and T4virus showed the amplification of 53.8% and ≥84.6% species, respectively. The present work is a case study performed to explore the possibility of use of taxon-specific primers for identification and taxonomic studies of newly isolated phages to supplement the TEM.

Prophage Integrase Typing Is a Useful Indicator of Genomic Diversity in Salmonella enterica

Salmonella enterica is a bacterial species that is a major cause of illness in humans and food-producing animals. S. enterica exhibits considerable inter-serovar diversity, as evidenced by the large number of host adapted serovars that have been identified. The development of methods to assess genome diversity in S. enterica will help to further define the limits of diversity in this foodborne pathogen. Thus, we evaluated a PCR assay, which targets prophage integrase genes, as a rapid method to investigate S. enterica genome diversity. To evaluate the PCR prophage integrase assay, 49 isolates of S. enterica were selected, including 19 clinical isolates from clonal serovars (Enteritidis and Heidelberg) that commonly cause human illness, and 30 isolates from food-associated Salmonella serovars that rarely cause human illness. The number of integrase genes identified by the PCR assay was compared to the number of integrase genes within intact prophages identified by whole genome sequencing and phage finding program PHASTER. The PCR assay identified a total of 147 prophage integrase genes within the 49 S. enterica genomes (79 integrase genes in the food-associated Salmonella isolates, 50 integrase genes in S. Enteritidis, and 18 integrase genes in S. Heidelberg). In comparison, whole genome sequencing and PHASTER identified a total of 75 prophage integrase genes within 102 intact prophages in the 49 S. enterica genomes (44 integrase genes in the food-associated Salmonella isolates, 21 integrase genes in S. Enteritidis, and 9 integrase genes in S. Heidelberg). Collectively, both the PCR assay and PHASTER identified the presence of a large diversity of prophage integrase genes in the food-associated isolates compared to the clinical isolates, thus indicating a high degree of diversity in the food-associated isolates, and confirming the clonal nature of S. Enteritidis and S. Heidelberg. Moreover, PHASTER revealed a diversity of 29 different types of prophages and 23 different integrase genes within the food-associated isolates, but only identified four different phages and integrase genes within clonal isolates of S. Enteritidis and S. Heidelberg. These results demonstrate the potential usefulness of PCR based detection of prophage integrase genes as a rapid indicator of genome diversity in S. enterica.

Finding a needle in the virus metagenome haystack--micro-metagenome analysis captures a snapshot of the diversity of a bacteriophage armoire

Viruses are ubiquitous in the oceans and critical components of marine microbial communities, regulating nutrient transfer to higher trophic levels or to the dissolved organic pool through lysis of host cells. Hydrothermal vent systems are oases of biological activity in the deep oceans, for which knowledge of biodiversity and its impact on global ocean biogeochemical cycling is still in its infancy. In order to gain biological insight into viral communities present in hydrothermal vent systems, we developed a method based on deep-sequencing of pulsed field gel electrophoretic bands representing key viral fractions present in seawater within and surrounding a hydrothermal plume derived from Loki's Castle vent field at the Arctic Mid-Ocean Ridge. The reduction in virus community complexity afforded by this novel approach enabled the near-complete reconstruction of a lambda-like phage genome from the virus fraction of the plume. Phylogenetic examination of distinct gene regions in this lambdoid phage genome unveiled diversity at loci encoding superinfection exclusion- and integrase-like proteins. This suggests the importance of fine-tuning lyosgenic conversion as a viral survival strategy, and provides insights into the nature of host-virus and virus-virus interactions, within hydrothermal plumes. By reducing the complexity of the viral community through targeted sequencing of prominent dsDNA viral fractions, this method has selectively mimicked virus dominance approaching that hitherto achieved only through culturing, thus enabling bioinformatic analysis to locate a lambdoid viral "needle" within the greater viral community "haystack". Such targeted analyses have great potential for accelerating the extraction of biological knowledge from diverse and poorly understood environmental viral communities.

Sensitive and specific detection of phaseolotoxigenic and nontoxigenic strains of Pseudomonas syringae pv. phaseolicola by TaqMan real-time PCR using site-specific recombinase gene sequences

Pseudomonas syringae pv. phaseolicola, the causative agent of halo blight, is the most important bacterial pathogen of bean. Both nontoxigenic (Tox(-)) and toxigenic (Tox+) strains of this pathogen cause halo blight in beans. However, nontoxigenic strains cannot be detected by currently available molecular and serological tools. In this study, a TaqMan probe and primer set were designed based on the phage integrase family site-specific recombinase of P. s. pv. phaseolicola 1448A because it is known that most site-specific recombinases are structurally and functionally diverse. The specificity of the probe and primers was evaluated using purified DNA from 29 isolates of 3 different pathovars of P. syringae. The probe and primer set were able to detect Tox(-) and Tox+ isolates of P. s. pv. Phaseolicola, but no other phytopathogenic bacteria. The assay was also able to detect at least 5 genome equivalents of cloned amplified target DNA, using purified DNA, or 7 colony forming unit (CFU) per reaction when using calibrated cell suspensions. Thus, the TaqMan real-time PCR-based method can be used for the rapid detection of both types of P. s. pv. Phaseolicola, and will potentially simplify and facilitate the diagnosis and monitoring of this pathogen, and guide plant disease management.

Bacteriophage lambda: a paradigm revisited

Bacteriophage lambda has an archetypal immunity system, which prevents the superinfection of its Escherichia coli lysogens. It is now known that superinfection can occur with toxigenic lambda-like phages at a high frequency, and here we demonstrate that the superinfection of a lambda lysogen can lead to the acquisition of additional lambda genomes, which was confirmed by Southern hybridization and quantitative PCR. As many as eight integration events were observed but at a very low frequency (6.4 x 10(-4)) and always as multiple insertions at the established primary integration site in E. coli. Sequence analysis of the complete immunity region demonstrated that these multiply infected lysogens were not immunity mutants. In conclusion, although lambda superinfection immunity can be confounded, it is a rare event.

Structures of the arm-type binding domains of HPI and HAI7 integrases

The structures of the N-terminal domains of two integrases of closely related but not identical asn tDNA-associated genomic islands, Yersinia HPI (high pathogenicity island; encoding siderophore yersiniabactin biosynthesis and transport) and an Erwinia carotovora genomic island with yet unknown function, HAI7, have been resolved. Both integrases utilize a novel four-stranded beta-sheet DNA-binding motif, in contrast to the known proteins that bind their DNA targets by means of three-stranded beta-sheets. Moreover, the beta-sheets in Int(HPI) and Int(HAI7) are longer than those in other integrases, and the structured helical N terminus is positioned perpendicularly to the large C-terminal helix. These differences strongly support the proposal that the integrases of the genomic islands make up a distinct evolutionary branch of the site-specific recombinases that utilize a unique DNA-binding mechanism.

PCR and partial sequencing of bacteriophage genomes

PCR is a quick and effective way of identifying the presence and 'affiliation' of bacteriophages, or phage-encoded genes from environmental samples, bacterial cells or purified viruses. The limitations are that you have to know what you are looking for in order to find it. Although the bacteriophage world does not have the advantage of a conserved gene, present in all members, there are many phage genes that do show nucleotide conservation even between phages which infect fairly divergent taxa. As more sequence data become available through both metagenomic approaches and the sequencing of complete bacteriophage genomes, PCR primers can be further refined and thus it should be an increasingly useful tool for bacteriophage biology.

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.

Sequence characterization and comparative analysis of three plasmids isolated from environmental Vibrio spp

The horizontal transfer of genes by mobile genetic elements such as plasmids and phages can accelerate genome diversification of Vibrio spp., affecting their physiology, pathogenicity, and ecological character. In this study, sequence analysis of three plasmids from Vibrio spp. previously isolated from salt marsh sediment revealed the remarkable diversity of these elements. Plasmids p0908 (81.4 kb), p23023 (52.5 kb), and p09022 (31.0 kb) had a predicted 99, 64, and 32 protein-coding sequences and G+C contents of 49.2%, 44.7%, and 42.4%, respectively. A phylogenetic tree based on concatenation of the host 16S rRNA and rpoA nucleotide sequences indicated p23023 and p09022 were isolated from strains most closely related to V. mediterranei and V. campbellii, respectively, while the host of p0908 forms a clade with V. fluvialis and V. furnissii. Many predicted proteins had amino acid identities to proteins of previously characterized phages and plasmids (24 to 94%). Predicted proteins with similarity to chromosomally encoded proteins included RecA, a nucleoid-associated protein (NdpA), a type IV helicase (UvrD), and multiple hypothetical proteins. Plasmid p0908 had striking similarity to enterobacteria phage P1, sharing genetic organization and amino acid identity for 23 predicted proteins. This study provides evidence of genetic exchange between Vibrio plasmids, phages, and chromosomes among diverse Vibrio spp.

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.