Characterization of the ELPhiS prophage from Salmonella enterica serovar Enteritidis strain LK5

Phages are a primary driving force behind the evolution of bacterial pathogens by transferring a variety of virulence genes into their hosts. Similar to other bacterial genomes, the Salmonella enterica serovar Enteritidis LK5 genome contains several regions that are homologous to phages. Although genomic analysis demonstrated the presence of prophages, it was unable to confirm which phage elements within the genome were viable. Genetic markers were used to tag one of the prophages in the genome to allow monitoring of phage induction. Commonly used laboratory strains of Salmonella were resistant to phage infection, and therefore a rapid screen was developed to identify susceptible hosts. This approach showed that a genetically tagged prophage, ELPhiS (Enteritidis lysogenic phage S), was capable of infecting Salmonella serovars that are diverse in host range and virulence and has the potential to laterally transfer genes between these serovars via lysogenic conversion. The rapid screen approach is adaptable to any system with a large collection of isolates and may be used to test the viability of prophages found by sequencing the genomes of various bacterial pathogens.     

Comparative genomics groups phages of Negativicutes and classical Firmicutes despite different Gram-staining properties

Negativicutes are gram-negative bacteria characterized by two cell membranes, but they are phylogenetically a side-branch of gram-positive Firmicutes that contain only a single membrane. We asked whether viruses (phages) infecting Negativicutes were horizontally acquired from gram-negative Proteobacteria, given the shared outer cell structure of their bacterial hosts, or if Negativicute phages co-evolved vertically with their hosts and thus resemble gram-positive Firmicute prophages. We predicted and characterized 485 prophages (mostly Caudovirales) from gram-negative Firmicute genomes plus 2977 prophages from other bacterial clades, and we used virome sequence data from 183 human stool samples to support our predictions. The majority of identified Negativicute prophages were lambdoids closer related to prophages from other Firmicutes than Proteobacteria by sequence relationship and genome organization (position of the lysis module). Only a single Mu-like candidate prophage and no clear P2-like prophages were identified in Negativicutes, both common in Proteobacteria. Given this collective evidence, it is unlikely that Negativicute phages were acquired from Proteobacteria. Sequence-related prophages, which occasionally harboured antibiotic resistance genes, were identified in two distinct Negativicute orders (Veillonellales and Acidaminococcales), possibly suggesting horizontal cross-order phage infection between human gut commensals. Our results reveal ancient genomic signatures of phage and bacteria co-evolution despite horizontal phage mobilization.     

Bacteriophage crosstalk: coordination of prophage induction by trans-acting antirepressors

Many species of bacteria harbor multiple prophages in their genomes. Prophages often carry genes that confer a selective advantage to the bacterium, typically during host colonization. Prophages can convert to infectious viruses through a process known as induction, which is relevant to the spread of bacterial virulence genes. The paradigm of prophage induction, as set by the phage Lambda model, sees the process initiated by the RecA-stimulated self-proteolysis of the phage repressor. Here we show that a large family of lambdoid prophages found in Salmonella genomes employs an alternative induction strategy. The repressors of these phages are not cleaved upon induction; rather, they are inactivated by the binding of small antirepressor proteins. Formation of the complex causes the repressor to dissociate from DNA. The antirepressor genes lie outside the immunity region and are under direct control of the LexA repressor, thus plugging prophage induction directly into the SOS response. GfoA and GfhA, the antirepressors of Salmonella prophages Gifsy-1 and Gifsy-3, each target both of these phages' repressors, GfoR and GfhR, even though the latter proteins recognize different operator sites and the two phages are heteroimmune. In contrast, the Gifsy-2 phage repressor, GtgR, is insensitive to GfoA and GfhA, but is inactivated by an antirepressor from the unrelated Fels-1 prophage (FsoA). This response is all the more surprising as FsoA is under the control of the Fels-1 repressor, not LexA, and plays no apparent role in Fels-1 induction, which occurs via a Lambda CI-like repressor cleavage mechanism. The ability of antirepressors to recognize non-cognate repressors allows coordination of induction of multiple prophages in polylysogenic strains. Identification of non-cleavable gfoR/gtgR homologues in a large variety of bacterial genomes (including most Escherichia coli genomes in the DNA database) suggests that antirepression-mediated induction is far more common than previously recognized.     

DNA Sequence Heterogeneity of Campylobacter jejuni CJIE4 Prophages and Expression of Prophage Genes

Campylobacter jejuni carry temperate bacteriophages that can affect the biology or virulence of the host bacterium. Known effects include genomic rearrangements and resistance to DNA transformation. C. jejuni prophage CJIE1 shows sequence variability and variability in the content of morons. Homologs of the CJIE1 prophage enhance both adherence and invasion to cells in culture and increase the expression of a specific subset of bacterial genes. Other C. jejuni temperate phages have so far not been well characterized. In this study we describe investigations into the DNA sequence variability and protein expression in a second prophage, CJIE4. CJIE4 sequences were obtained de novo from DNA sequencing of five C. jejuni isolates, as well as from whole genome sequences submitted to GenBank by other research groups. These CJIE4 DNA sequences were heterogenous, with several different insertions/deletions (indels) in different parts of the prophage genome. Two variants of a 3–4 kb region inserted within CJIE4 had different gene content that distinguished two major conserved CJIE4 prophage families. Additional indels were detected throughout the prophage. Detection of proteins in the five isolates characterized in our laboratory in isobaric Tags for Relative and Absolute Quantitation (iTRAQ) experiments indicated that prophage proteins within each of the two large indel variants were expressed during growth of the bacteria on Mueller Hinton agar plates. These proteins included the extracellular DNase associated with resistance to DNA transformation and prophage repressor proteins. Other proteins associated with known or suspected roles in prophage biology were also expressed from CJIE4, including capsid protein, the phage integrase, and MazF, a type II toxin-antitoxin system protein. Together with the results previously obtained for the CJIE1 prophage these results demonstrate that sequence variability and expression of moron genes are both general properties of temperate bacteriophages in C. jejuni.     

Phage regulatory circuits and virulence gene expression

In many pathogenic bacteria, genes that encode virulence factors are located in the genomes of prophages. Clearly bacteriophages are important vectors for disseminating virulence genes, but, in addition, do phage regulatory circuits contribute to expression of these genes? Phages of the lambda family that have genes encoding Shiga toxin are found in certain pathogenic Escherichia coli (known as Shiga toxin producing E. coli) and the filamentous phage CTXphi, that carries genes encoding cholera toxin (CTX), is found in Vibrio cholerae. Both the lambda and CTXphi phages have repressor systems that maintain their respective prophages in a quiescent state, and in both types of prophages this repressed state is abolished when the host cell SOS response is activated. In the lambda type of prophages, only binding of the phage-encoded repressor is involved in repression and this repressor ultimately controls Shiga toxin production and/or release. In the CTXphi prophage, binding of LexA, the bacterial regulator of SOS, in addition to binding of the repressor is involved in repression; the repressor has only limited control over CTX production and has no influence on its release.     

Evolution along the parasitism-mutualism continuum determines the genetic repertoire of prophages

Integrated into their bacterial hosts’ genomes, prophage sequences exhibit a wide diversity of length and gene content, from highly degraded cryptic sequences to intact, functional prophages that retain a full complement of lytic-function genes. We apply three approaches—bioinformatics, analytical modelling and computational simulation—to understand the diverse gene content of prophages. In the bioinformatics work, we examine the distributions of over 50,000 annotated prophage genes identified in 1384 prophage sequences, comparing the gene repertoires of intact and incomplete prophages. These data indicate that genes involved in the replication, packaging, and release of phage particles have been preferentially lost in incomplete prophages, while tail fiber, transposase and integrase genes are significantly enriched. Consistent with these results, our mathematical and computational approaches predict that genes involved in phage lytic function are preferentially lost, resulting in shorter prophages that often retain genes that benefit the host. Informed by these models, we offer novel hypotheses for the enrichment of integrase and transposase genes in cryptic prophages. Overall, we demonstrate that functional and cryptic prophages represent a diversity of genetic sequences that evolve along a parasitism-mutualism continuum.     

Marine phage genomics

Marine phages are the most abundant biological entities in the oceans. They play important roles in carbon cycling through marine food webs, gene transfer by transduction and conversion of hosts by lysogeny. The handful of marine phage genomes that have been sequenced to date, along with prophages in marine bacterial genomes, and partial sequencing of uncultivated phages are yielding glimpses of the tremendous diversity and physiological potential of the marine phage community. Common gene modules in diverse phages are providing the information necessary to make evolutionary comparisons. Finally, deciphering phage genomes is providing clues about the adaptive response of phages and their hosts to environmental cues.     

Bacteriophage Mu genome sequence: analysis and comparison with Mu-like prophages in Haemophilus, Neisseria and Deinococcus

We report the complete 36,717 bp genome sequence of bacteriophage Mu and provide an analysis of the sequence, both with regard to the new genes and other genetic features revealed by the sequence itself and by a comparison to eight complete or nearly complete Mu-like prophage genomes found in the genomes of a diverse group of bacteria. The comparative studies confirm that members of the Mu-related family of phage genomes are genetically mosaic with respect to each other, as seen in other groups of phages such as the phage lambda-related group of phages of enteric hosts and the phage L5-related group of mycobacteriophages. Mu also possesses segments of similarity, typically gene-sized, to genomes of otherwise non-Mu-like phages. The comparisons show that some well-known features of the Mu genome, including the invertible segment encoding tail fiber sequences, are not present in most members of the Mu genome sequence family examined here, suggesting that their presence may be relatively volatile over evolutionary time.The head and tail-encoding structural genes of Mu have only very weak similarity to the corresponding genes of other well-studied phage types. However, these weak similarities, and in some cases biochemical data, can be used to establish tentative functional assignments for 12 of the head and tail genes. These assignments are strongly supported by the fact that the order of gene functions assigned in this way conforms to the strongly conserved order of head and tail genes established in a wide variety of other phages. We show that the Mu head assembly scaffolding protein is encoded by a gene nested in-frame within the C-terminal half of another gene that encodes the putative head maturation protease. This is reminiscent of the arrangement established for phage lambda.     

Morphological and genetic diversity of temperate phages in Clostridium difficile

Eight temperate phages were characterized after mitomycin C induction of six Clostridium difficile isolates corresponding to six distinct PCR ribotypes. The hypervirulent C. difficile strain responsible for a multi-institutional outbreak (NAP1/027 or QCD-32g58) was among these prophage-containing strains. Observation of the crude lysates by transmission electron microscopy (TEM) revealed the presence of three phages with isometric capsids and long contractile tails (Myoviridae family), as well as five phages with long noncontractile tails (Siphoviridae family). TEM analyses also revealed the presence of a significant number of phage tail-like particles in all the lysates. Southern hybridization experiments with restricted prophage DNA showed that C. difficile phages belonging to the family Myoviridae are highly similar and most likely related to previously described prophages phiC2, phiC5, and phiCD119. On the other hand, members of the Siphoviridae phage family are more genetically divergent, suggesting that they originated from distantly related ancestors. Our data thus suggest that there are at least three genetically distinct groups of temperate phages in C. difficile; one group is composed of highly related myophages, and the other two groups are composed of more genetically heterogeneous siphophages. Finally, no gene homologous to genes encoding C. difficile toxins or toxin regulators could be identified in the genomes of these phages using DNA hybridization. Interestingly, each unique phage restriction profile correlated with a specific C. difficile PCR ribotype.     

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.     

Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences

Phage_Finder, a heuristic computer program, was created to identify prophage regions in completed bacterial genomes. Using a test dataset of 42 bacterial genomes whose prophages have been manually identified, Phage_Finder found 91% of the regions, resulting in 7% false positive and 9% false negative prophages. A search of 302 complete bacterial genomes predicted 403 putative prophage regions, accounting for 2.7% of the total bacterial DNA. Analysis of the 285 putative attachment sites revealed tRNAs are targets for integration slightly more frequently (33%) than intergenic (31%) or intragenic (28%) regions, while tmRNAs were targeted in 8% of the regions. The most popular tRNA targets were Arg, Leu, Ser and Thr. Mapping of the insertion point on a consensus tRNA molecule revealed novel insertion points on the 5' side of the D loop, the 3' side of the anticodon loop and the anticodon. A novel method of constructing phylogenetic trees of phages and prophages was developed based on the mean of the BLAST score ratio (BSR) of the phage/prophage proteomes. This method verified many known bacteriophage groups, making this a useful tool for predicting the relationships of prophages from bacterial genomes.     

Isolation of Phages for Phage Therapy: A Comparison of Spot Tests and Efficiency of Plating Analyses for Determination of Host Range and Efficacy

Phage therapy, treating bacterial infections with bacteriophages, could be a future alternative to antibiotic treatment of bacterial infections. There are, however, several problems to be solved, mainly associated to the biology of phages, the interaction between phages and their bacterial hosts, but also to the vast variation of pathogenic bacteria which implies that large numbers of different phages are going to be needed. All of these phages must under present regulation of medical products undergo extensive clinical testing before they can be applied. It will consequently be of great economic importance that effective and versatile phages are selected and collected into phage libraries, i.e., the selection must be carried out in a way that it results in highly virulent phages with broad host ranges. We have isolated phages using the Escherichia coli reference (ECOR) collection and compared two methods, spot testing and efficiency of plating (EOP), which are frequently used to identify phages suitable for phage therapy. The analyses of the differences between the two methods show that spot tests often overestimate both the overall virulence and the host range and that the results are not correlated to the results of EOP assays. The conclusion is that single dilution spot tests cannot be used for identification and selection of phages to a phage library and should be replaced by EOP assays. The difference between the two methods can be caused by many factors. We have analysed if the differences and lack of correlation could be caused by lysis from without, bacteriocins in the phage lysate, or by the presence of prophages harbouring genes coding for phage resistance systems in the genomes of the bacteria in the ECOR collection.     

Comparative genomics and evolution of the tailed-bacteriophages

The number of completely sequenced tailed-bacteriophage genomes that have been published increased to more than 125 last year. The comparison of these genomes has brought their highly mosaic nature into much sharper focus. Furthermore, reports of the complete sequences of about 150 bacterial genomes have shown that the many prophage and parts thereof that reside in these bacterial genomes must comprise a significant fraction of Earth's phage gene pool. These phage and prophage genomes are fertile ground for attempts to deduce the nature of viral evolutionary processes, and such analyses have made it clear that these phage have enjoyed a significant level of horizontal exchange of genetic information throughout their long histories. The strength of these evolutionary deductions rests largely on the extensive knowledge that has accumulated during intensive study into the molecular nature of the life cycles of a few 'model system' phages over the past half century. Recent molecular studies of phages other than these model system phages have made it clear that much remains to be learnt about the variety of lifestyle strategies utilized by the tailed-phage.     

Prophages integrating into prophages: A mechanism to accumulate type III secretion effector genes and duplicate Shiga toxin-encoding prophages in Escherichia coli

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.     

The pKO2 linear plasmid prophage of Klebsiella oxytoca

Temperate bacteriophages with plasmid prophages are uncommon in nature, and of these only phages N15 and PY54 are known to have a linear plasmid prophage with closed hairpin telomeres. We report here the complete nucleotide sequence of the 51,601-bp Klebsiella oxytoca linear plasmid pKO2, and we demonstrate experimentally that it is also a prophage. We call this bacteriophage phiKO2. An analysis of the 64 predicted phiKO2 genes indicate that it is a fairly close relative of phage N15; they share a mosaic relationship that is typical of different members of double-stranded DNA tailed-phage groups. Although the head, tail shaft, and lysis genes are not recognizably homologous between these phages, other genes such as the plasmid partitioning, replicase, prophage repressor, and protelomerase genes (and their putative targets) are so similar that we predict that they must have nearly identical DNA binding specificities. The phiKO2 virion is unusual in that its phage lambda-like tails have an exceptionally long (3,433 amino acids) central tip tail fiber protein. The phiKO2 genome also carries putative homologues of bacterial dinI and umuD genes, both of which are involved in the host SOS response. We show that these divergently transcribed genes are regulated by LexA protein binding to a single target site that overlaps both promoters.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.     

Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.