Comparative genomics of the T4-Like Escherichia coli phage JS98: implications for the evolution of T4 phages

About 130 kb of sequence information was obtained from the coliphage JS98 isolated from the stool of a pediatric diarrhea patient in Bangladesh. The DNA shared up to 81% base pair identity with phage T4. The most conserved regions between JS98 and T4 were the structural genes, but their degree of conservation was not uniform. The head genes showed the highest sequence conservation, followed by the tail, baseplate, and tail fiber genes. Many tail fiber genes shared only protein sequence identity. Except for the insertion of endonuclease genes in T4 and gene 24 duplication in JS98, the structural gene maps of the two phages were colinear. The receptor-recognizing tail fiber proteins gp37 and gp38 were only distantly related to T4, but shared up to 83% amino acid identity to other T6-like phages, suggesting lateral gene transfer. A greater degree of variability was seen between JS98 and T4 over DNA replication and DNA transaction genes. While most of these genes came in the same order and shared up to 76% protein sequence identity, a few rearrangements, insertions, and replacements of genes were observed. Many putative gene insertions in the DNA replication module of T4 were flanked by intron-related endonuclease genes, suggesting mobile DNA elements. A hotspot of genome diversification was located downstream of the DNA polymerase gene 43 and the DNA binding gene 32. Comparative genomics of 100-kb genome sequence revealed that T4-like phages diversify more by the accumulation of point mutations and occasional gene duplication events than by modular exchanges.     

Campylobacter group II phage CP21 is the prototype of a new subgroup revealing a distinct modular genome organization and host specificity

Background

The application of phages is a promising tool to reduce the number of Campylobacter along the food chain. Besides the efficacy against a broad range of strains, phages have to be safe in terms of their genomes. Thus far, no genes with pathogenic potential (e.g., genes encoding virulence factors) have been detected in Campylobacter phages. However, preliminary studies suggested that the genomes of group II phages may be diverse and prone to genomic rearrangements.

Results

We determined and analysed the genomic sequence (182,761 bp) of group II phage CP21 that is closely related to the already characterized group II phages CP220 and CPt10. The genomes of these phages are comprised of four modules separated by very similar repeat regions, some of which harbouring open reading frames (ORFs). Though, the arrangement of the modules and the location of some ORFs on the genomes are different in CP21 and in CP220/CPt10. In this work, a PCR system was established to study the modular genome organization of other group II phages demonstrating that they belong to different subgroups of the CP220-like virus genus, the prototypes of which are CP21 and CP220. The subgroups revealed different restriction patterns and, interestingly enough, also distinct host specificities, tail fiber proteins and tRNA genes. We additionally analysed the genome of group II phage vB_CcoM-IBB_35 (IBB_35) for which to date only five individual contigs could be determined. We show that the contigs represent modules linked by long repeat regions enclosing some yet not identified ORFs (e.g., for a head completion protein). The data suggest that IBB_35 is a member of the CP220 subgroup.

Conclusion

Campylobacter group II phages are diverse regarding their genome organization. Since all hitherto characterized group II phages contain numerous genes for transposases and homing endonucleases as well as similar repeat regions, it cannot be excluded that these phages are genetically unstable. To answer this question, further experiments and sequencing of more group II phages should be performed.

Electronic supplementary material

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

T4-Like genome organization of the Escherichia coli O157:H7 lytic phage AR1

We report the genome organization and analysis of the first completely sequenced T4-like phage, AR1, of Escherichia coli O157:H7. Unlike most of the other sequenced phages of O157:H7, which belong to the temperate Podoviridae and Siphoviridae families, AR1 is a T4-like phage known to efficiently infect this pathogenic bacterial strain. The 167,435-bp AR1 genome is currently the largest among all the sequenced E. coli O157:H7 phages. It carries a total of 281 potential open reading frames (ORFs) and 10 putative tRNA genes. Of these, 126 predicted proteins could be classified into six viral orthologous group categories, with at least 18 proteins of the structural protein category having been detected by tandem mass spectrometry. Comparative genomic analysis of AR1 and four other completely sequenced T4-like genomes (RB32, RB69, T4, and JS98) indicated that they share a well-organized and highly conserved core genome, particularly in the regions encoding DNA replication and virion structural proteins. The major diverse features between these phages include the modules of distal tail fibers and the types and numbers of internal proteins, tRNA genes, and mobile elements. Codon usage analysis suggested that the presence of AR1-encoded tRNAs may be relevant to the codon usage of structural proteins. Furthermore, protein sequence analysis of AR1 gp37, a potential receptor binding protein, indicated that eight residues in the C terminus are unique to O157:H7 T4-like phages AR1 and PP01. These residues are known to be located in the T4 receptor recognition domain, and they may contribute to specificity for adsorption to the O157:H7 strain.     

Phylogeny of the major head and tail genes of the wide-ranging T4-type bacteriophages

We examined a number of bacteriophages with T4-type morphology that propagate in different genera of enterobacteria, Aeromonas, Burkholderia, and Vibrio. Most of these phages had a prolate icosahedral head, a contractile tail, and a genome size that was similar to that of T4. A few of them had more elongated heads and larger genomes. All these phages are phylogenetically related, since they each had sequences homologous to the capsid gene (gene 23), tail sheath gene (gene 18), and tail tube gene (gene 19) of T4. On the basis of the sequence comparison of their virion genes, the T4-type phages can be classified into three subgroups with increasing divergence from T4: the T-evens, pseudoT-evens, and schizoT-evens. In general, the phages that infect closely related host species have virion genes that are phylogenetically closer to each other than those of phages that infect distantly related hosts. However, some of the phages appear to be chimeras, indicating that, at least occasionally, some genetic shuffling has occurred between the different T4-type subgroups. The compilation of a number of gene 23 sequences reveals a pattern of conserved motifs separated by sequences that differ in the T4-type subgroups. Such variable patches in the gene 23 sequences may determine the size of the virion head and consequently the viral genome length. This sequence analysis provides molecular evidence that phages related to T4 are widespread in the biosphere and diverged from a common ancestor in acquiring the ability to infect different host bacteria and to occupy new ecological niches.     

Genome analysis of phage JS98 defines a fourth major subgroup of T4-like phages in Escherichia coli

Numerous T4-like Escherichia coli phages were isolated from human stool and environmental wastewater samples in Bangladesh and Switzerland. The sequences of the major head gene (g23) revealed that these coliphages could be placed into four subgroups, represented by the phages T4, RB69, RB49, and JS98. Thus, JS98 defines a new major subgroup of E. coli T4-like phages. We conducted an analysis of the 169-kb JS98 genome sequence. Overall, 198 of the 266 JS98 open reading frames (ORFs) shared amino acid sequence identity with the reference T4 phage, 41 shared identity with other T4-like phages, and 27 ORFs lacked any database matches. Genes on the plus strand encoded virion proteins, which showed moderate to high sequence identity with T4 proteins. The right genome half of JS98 showed a higher degree of sequence conservation with T4 and RB69, even for the nonstructural genes, than did the left genome half, containing exclusively nonstructural genes. Most of the JS98-specific genes were found in the left genome half. Two came as a hypervariability cluster, but most represented isolated genes, suggesting that they were acquired separately in multiple acquisition events. No evidence for DNA exchange between JS98 phage and the E. coli host genome or coliphages other than T4 was observed. No undesired genes which could compromise its medical use were detected in the JS98 genome sequence.     

Primary Isolation Strain Determines Both Phage Type and Receptors Recognised by Campylobacter jejuni Bacteriophages

In this study we isolated novel bacteriophages, infecting the zoonotic bacterium Campylobacter jejuni. These phages may be used in phage therapy of C. jejuni colonized poultry to prevent spreading of the bacteria to meat products causing disease in humans. Many C. jejuni phages have been isolated using NCTC12662 as the indicator strain, which may have biased the selection of phages. A large group of C. jejuni phages rely on the highly diverse capsular polysaccharide (CPS) for infection and recent work identified the O-methyl phosphoramidate modification (MeOPN) of CPS as a phage receptor. We therefore chose seven C. jejuni strains each expressing different CPS structures as indicator strains in a large screening for phages in samples collected from free-range poultry farms. Forty-three phages were isolated using C. jejuni NCTC12658, NCTC12662 and RM1221 as host strains and 20 distinct phages were identified based on host range analysis and genome restriction profiles. Most phages were isolated using C. jejuni strains NCTC12662 and RM1221 and interestingly phage genome size (140 kb vs. 190 kb), host range and morphological appearance correlated with the isolation strain. Thus, according to C. jejuni phage grouping, NCTC12662 and NCTC12658 selected for CP81-type phages, while RM1221 selected for CP220-type phages. Furthermore, using acapsular ∆kpsM mutants we demonstrated that phages isolated on NCTC12658 and NCTC12662 were dependent on the capsule for infection. In contrast, CP220-type phages isolated on RM1221 were unable to infect non-motile ∆motA mutants, hence requiring motility for successful infection. Hence, the primary phage isolation strain determines both phage type (CP81 or CP220) as well as receptors (CPS or flagella) recognised by the isolated phages.     

Plasticity of the gene functions for DNA replication in the T4-like phages

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.     

Functional interactions of the genomes of Shigella sonnei phages and Escherichia coli phage T4 in mixed infection]

A comparative study of Shigella sonnei phages U and G and Escherichia coli phage T4 has shown that enzymes coded for by the Sh. sonnei phages can functionally substitute for some T4-coded products. This finding in indicative of an evolutionary relationship between T-even phages and disenteric phages U and G. The U phage is uncapable to compensate amber mutants for the genes that control the conversion of cytosine into 5-hydroxymethyl cytosine (5-HMC) and the glucosylation of the latter, which agrees with our earlier finding that the U phage DNA contains no 5-HMC. U and G phages are also found to exclude the T4 phage in the course of mixed infection.     

Genome and proteome of Campylobacter jejuni bacteriophage NCTC 12673

Campylobacter jejuni continues to be the leading cause of bacterial food-borne illness worldwide, so improvements to current methods used for bacterial detection and disease prevention are needed. We describe here the genome and proteome of C. jejuni bacteriophage NCTC 12673 and the exploitation of its receptor-binding protein for specific bacterial detection. Remarkably, the 135-kb Myoviridae genome of NCTC 12673 differs greatly from any other proteobacterial phage genome described (including C. jejuni phages CP220 and CPt10) and instead shows closest homology to the cyanobacterial T4-related myophages. The phage genome contains 172 putative open reading frames, including 12 homing endonucleases, no visible means of packaging, and a putative trans-splicing intein. The phage DNA appears to be strongly associated with a protein that interfered with PCR amplification and estimation of the phage genome mass by pulsed-field gel electrophoresis. Identification and analyses of the receptor-binding protein (Gp48) revealed features common to the Salmonella enterica P22 phage tailspike protein, including the ability to specifically recognize a host organism. Bacteriophage receptor-binding proteins may offer promising alternatives for use in pathogen detection platforms.     

Genome of Xanthomonas oryzae bacteriophage Xp10: an odd T-odd phage

Xp10 is a lytic bacteriophage of the phytopathogenic bacterium Xanthomonas oryzae. Though morphologically Xp10 belongs to the Syphoviridae family, it encodes its own single-subunit RNA polymerase characteristic of T7-like phages of the Podoviridae family. Here, we report the determination and analysis of the 44,373 bp sequence of the Xp10 genome. The genome is a linear, double-stranded DNA molecule with 3' cohesive overhangs and no terminal repeats or redundancies. Half of the Xp10 genome contains genes coding for structural proteins and host lysis functions in an arrangement typical for temperate dairy phages that are related to the Escherichia coli lambda phage. The other half of the Xp10 genome contains genes coding for factors of host gene expression shut-off, enzymes of viral genome replication and expression. The two groups of genes are transcribed divergently and separated by a regulatory region, which contains divergent promoters recognized by the host RNA polymerase. Xp10 has apparently arisen through a recombination between genomes of widely different phages. Further evidence of extensive gene flux in the evolution of Xp10 includes a high fraction (10%) of genes derived from an HNH-family endonuclease, and a DNA-dependent DNA polymerase that is closer to a homolog from Leishmania than to DNA polymerases from other phages or bacteria.     

The genome of bacteriophage K1F, a T7-like phage that has acquired the ability to replicate on K1 strains of Escherichia coli

Bacteriophage K1F specifically infects Escherichia coli strains that produce the K1 polysaccharide capsule. Like several other K1 capsule-specific phages, K1F encodes an endo-neuraminidase (endosialidase) that is part of the tail structure which allows the phage to recognize and degrade the polysaccharide capsule. The complete nucleotide sequence of the K1F genome reveals that it is closely related to bacteriophage T7 in both genome organization and sequence similarity. The most striking difference between the two phages is that K1F encodes the endosialidase in the analogous position to the T7 tail fiber gene. This is in contrast with bacteriophage K1-5, another K1-specific phage, which encodes a very similar endosialidase which is part of a tail gene "module" at the end of the phage genome. It appears that diverse phages have acquired endosialidase genes by horizontal gene transfer and that these genes or gene products have adapted to different genome and virion architectures.     

Characterization of a novel T4-type Stenotrophomonas maltophilia virulent phage Smp14

Stenotrophomonas maltophilia (Sm), with most of the isolates being resistant to multidrugs, is an opportunistic bacterium causing nosocomial infections. In this study, a novel virulent Sm phage, Smp14, was characterized. Electron microscopy showed that Smp14 resembled members of Myoviridae and adsorbed to poles of the host cells during infection. It lysed 37 of 87 clinical Sm isolates in spot test, displayed a latent period of ca. 20 min, and had a burst size of ca. 150. Its genome (estimated to be 160 kb by PFGE), containing m4C and two unknown modified bases other than m5C and m6A as identified by HPLC, resisted to digestion with many restriction endonucleases except MseI. These properties indicate that it is a novel Sm phage distinct from the previously reported phiSMA5 which has a genome of 250 kb digestible with various restriction enzymes. Sequencing of a 16 kb region revealed 12 ORFs encoding structural proteins sharing 15-45% identities with the homologues from T4-type phages. SDS-PAGE displayed 20 virion proteins, with the most abundant one being the 39 kDa major capsid protein (gp23), which had the N-terminal 52 amino acids removed. Phylogenetic analysis based on gp23 classified Smp14 into a novel single-membered T4-type subgroup.     

Methyltransferases acquired by lactococcal 936-type phage provide protection against restriction endonuclease activity

Background

So-called 936-type phages are among the most frequently isolated phages in dairy facilities utilising Lactococcus lactis starter cultures. Despite extensive efforts to control phage proliferation and decades of research, these phages continue to negatively impact cheese production in terms of the final product quality and consequently, monetary return.

Results

Whole genome sequencing and in silico analysis of three 936-type phage genomes identified several putative (orphan) methyltransferase (MTase)-encoding genes located within the packaging and replication regions of the genome. Utilising SMRT sequencing, methylome analysis was performed on all three phages, allowing the identification of adenine modifications consistent with N-6 methyladenine sequence methylation, which in some cases could be attributed to these phage-encoded MTases. Heterologous gene expression revealed that M.Phi145I/M.Phi93I and M.Phi93DAM, encoded by genes located within the packaging module, provide protection against the restriction enzymes HphI and DpnII, respectively, representing the first functional MTases identified in members of 936-type phages.

Conclusions

SMRT sequencing technology enabled the identification of the target motifs of MTases encoded by the genomes of three lytic 936-type phages and these MTases represent the first functional MTases identified in this species of phage. The presence of these MTase-encoding genes on 936-type phage genomes is assumed to represent an adaptive response to circumvent host encoded restriction-modification systems thereby increasing the fitness of the phages in a dynamic dairy environment.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-831) contains supplementary material, which is available to authorized users.     

Genomic and structural analysis of Syn9, a cyanophage infecting marine Prochlorococcus and Synechococcus

Cyanobacteriophage Syn9 is a large, contractile-tailed bacteriophage infecting the widespread, numerically dominant marine cyanobacteria of the genera Prochlorococcus and Synechococcus . Its 177 300 bp genome sequence encodes 226 putative proteins and six tRNAs. Experimental and computational analyses identified genes likely involved in virion formation, nucleotide synthesis, and DNA replication and repair. Syn9 shows significant mosaicism when compared with related cyanophages S-PM2, P-SSM2 and P-SSM4, although shared genes show strong purifying selection and evidence for large population sizes relative to other phages. Related to coliphage T4 – which shares 19% of Syn9's genes – Syn9 shows evidence for different patterns of DNA replication and uses homologous proteins to assemble capsids with a different overall structure that shares topology with phage SPO1 and herpes virus. Noteworthy bacteria-related sequences in the Syn9 genome potentially encode subunits of the photosynthetic reaction centre, electron transport proteins, three pentose pathway enzymes and two tryptophan halogenases. These genes suggest that Syn9 is well adapted to the physiology of its photosynthetic hosts and may affect the evolution of these sequences within marine cyanobacteria.     

Phage T4 SegB protein is a homing endonuclease required for the preferred inheritance of T4 tRNA gene region occurring in co-infection with a related phage

Homing endonucleases initiate nonreciprocal transfer of DNA segments containing their own genes and the flanking sequences by cleaving the recipient DNA. Bacteriophage T4 segB gene, which is located in a cluster of tRNA genes, encodes a protein of unknown function, homologous to homing endonucleases of the GIY-YIG family. We demonstrate that SegB protein is a site-specific endonuclease, which produces mostly 3′ 2-nt protruding ends at its DNA cleavage site. Analysis of SegB cleavage sites suggests that SegB recognizes a 27-bp sequence. It contains 11-bp conserved sequence, which corresponds to a conserved motif of tRNA TψC stem-loop, whereas the remainder of the recognition site is rather degenerate. T4-related phages T2L, RB1 and RB3 contain tRNA gene regions that are homologous to that of phage T4 but lack segB gene and several tRNA genes. In co-infections of phages T4 and T2L, segB gene is inherited with nearly 100% of efficiency. The preferred inheritance depends absolutely on the segB gene integrity and is accompanied by the loss of the T2L tRNA gene region markers. We suggest that SegB is a homing endonuclease that functions to ensure spreading of its own gene and the surrounding tRNA genes among T4-related phages.