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

Lytic Capsule-Specific Acinetobacter Bacteriophages Encoding Polysaccharide-Degrading Enzymes

The genus Acinetobacter comprises both environmental and clinically relevant species associated with hospital-acquired infections. Among them, Acinetobacter baumannii is a critical priority bacterial pathogen, for which the research and development of new strategies for antimicrobial treatment are urgently needed. Acinetobacter spp. produce a variety of structurally diverse capsular polysaccharides (CPSs), which surround the bacterial cells with a thick protective layer. These surface structures are primary receptors for capsule-specific bacteriophages, that is, phages carrying tailspikes with CPS-depolymerizing/modifying activities. Phage tailspike proteins (TSPs) exhibit hydrolase, lyase, or esterase activities toward the corresponding CPSs of a certain structure. In this study, the data on all lytic capsule-specific phages infecting Acinetobacter spp. with genomes deposited in the NCBI GenBank database by January 2024 were summarized. Among the 149 identified TSPs encoded in the genomes of 143 phages, the capsular specificity (K specificity) of 46 proteins has been experimentally determined or predicted previously. The specificity of 63 TSPs toward CPSs, produced by various Acinetobacter K types, was predicted in this study using a bioinformatic analysis. A comprehensive phylogenetic analysis confirmed the prediction and revealed the possibility of the genetic exchange of gene regions corresponding to the CPS-recognizing/degrading parts of different TSPs between morphologically and taxonomically distant groups of capsule-specific Acinetobacter phages.

Adaptive Capacity of a DNA Polymerase Clamp-loader ATPase Complex

The ability of mutations to facilitate adaptation is central to evolution. To understand how mutations can lead to functional adaptation in a complex molecular machine, we created a defective version of the T4 clamp-loader complex, which is essential for DNA replication. This variant, which is ∼5,000-fold less active than the wild type, was made by replacing the catalytic domains with those from another phage. A directed-evolution experiment revealed that multiple substitutions to a single negatively charged residue in the chimeric clamp loader-Asp 86-restore fitness to within ∼20-fold of wild type. These mutations remove an adventitious electrostatic repulsive interaction between Asp 86 and the sliding clamp. Thus, the fitness decrease of the chimeric clamp loader is caused by a reduction in affinity between the clamp loader and the clamp. Deep mutagenesis shows that the reduced fitness of the chimeric clamp loader is also compensated for by lysine and arginine substitutions of several DNA-proximal residues in the clamp loader or the sliding clamp. Our results demonstrate that there is a latent capacity for increasing the affinity of the clamp loader for DNA and the sliding clamp, such that even single-point mutations can readily compensate for the loss of function due to suboptimal interactions elsewhere.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Isolation, Characterization and Genomic Analysis of a Novel Jumbo Phage, AerS_266, That Infects Aeromonas salmonicida

Aeromonas salmonicida is the causative agent of septicemia in fish, and it is associated with significant economic losses in the aquaculture industry. While piscine Aeromonas infections are mainly treated with antibiotics, the emergence of resistance in bacterial populations requires the development of alternative methods of treatment. The use of phages can be one of them. A novel A. salmonicida jumbo phage, AerS_266, was isolated and characterized. This phage infects only mesophilic A. salmonicida strains and demonstrates a slow lytic life cycle. Its genome contains 243,674 bp and 253 putative genes: 84 encode proteins with predicted functions, and 3 correspond to tRNAs. Genes encoding two multisubunit RNA polymerases, chimallin and PhuZ, were identified, and AerS_266 was thus defined as a phiKZ-like phage. While similar phages with genomes >200 kb specific to Aeromonas hydrophila and Aeromonas veronii have been previously described, AerS_266 is the first phiKZ-like phage found to infect A. salmonicida.

An OmpW-dependent T4-like phage infects Serratia sp. ATCC 39006

Serratia sp. ATCC 39006 is a Gram-negative bacterium that has been used to study the function of phage defences, such as CRISPR-Cas, and phage counter-defence mechanisms. To expand our phage collection to study the phage-host interaction with Serratia sp. ATCC 39006, we isolated the T4-like myovirus LC53 in Ōtepoti Dunedin, Aotearoa New Zealand. Morphological, phenotypic and genomic characterization revealed that LC53 is virulent and similar to other Serratia, Erwinia and Kosakonia phages belonging to the genus Winklervirus. Using a transposon mutant library, we identified the host ompW gene as essential for phage infection, suggesting that it encodes the phage receptor. The genome of LC53 encodes all the characteristic T4-like core proteins involved in phage DNA replication and generation of viral particles. Furthermore, our bioinformatic analysis suggests that the transcriptional organization of LC53 is similar to that of Escherichia coli phage T4. Importantly, LC53 encodes 18 tRNAs, which likely compensate for differences in GC content between phage and host genomes. Overall, this study describes a newly isolated phage infecting Serratia sp. ATCC 39006 that expands the diversity of phages available to study phage-host interactions.

Comparative genomic analysis of dwarf Vibrio myoviruses defines a conserved gene cluster for successful phage infection

Tailed bacteriophages have been at the center of attention, not only for their ability to infect and kill pathogenic bacteria but also due to their peculiar and intriguing complex contractile tail structure. Tailed bacteriophages with contractile tails are known to have a Myoviridae morphotype and are members of the order Caudovirales. Large bacteriophages with a genome larger than 150 kbp have been studied for their ability to use multiple infection and lysis strategies to replicate more efficiently. On the other hand, smaller bacteriophages with fewer genes are represented in the GenBank database in greater numbers, and have several genes with unknown function. Isolation and molecular characterization of a newly reported bacteriophage named Athena1 revealed that it is a strongly lytic bacteriophage with a genome size of 39,826 bp. This prompted us to perform a comparative genomic analysis of Vibrio myoviruses with a genome size of no more than 50 kbp. The results revealed a pattern of genomic organization that includes sets of genes responsible for virion morphogenesis, replication/recombination of DNA, and lysis/lysogeny switching. By studying phylogenetic gene markers, we were able to draw conclusions about evolutionary events that shaped the genomic mosaicism of these phages, pinpointing the importance of a conserved organization of the genomic region encoding the baseplate protein for successful infection of Gram-negative bacteria. In addition, we propose the creation of new genera for dwarf Vibrio myoviruses. Comparative genomics of phages infecting aquatic bacteria could provide information that is useful for combating fish pathogens in aquaculture, using novel strategies.

Characterization and genome sequence of the genetically unique Escherichia bacteriophage vB_EcoM_IME392

In this study, a novel Escherichia coli-specific bacteriophage, vB_EcoM_IME392, was isolated from chicken farm sewage in Qingdao, China. The genome of IME392 was found by next-generation sequencing to be 116,460 base pairs in length with a G+C content of 45.4% (GenBank accession number MH719082). BLASTn results revealed that only 2% of the genome sequence of IME392 shows sequence similarity to known phage sequences in the GenBank database, which indicates that IME392 is a novel bacteriophage. Transmission electron microscopy showed that IME392 belongs to the family Myoviridae. The host range, the multiplicity of infection, and a one-step growth curve were also determined.

Influence of Non-canonical DNA Bases on the Genomic Diversity of Tevenvirinae

The Tevenvirinae viruses are some of the most common viruses on Earth. Representatives of this subfamily have long been used in the molecular biology studies as model organisms – since the emergence of the discipline. Tevenvirinae are promising agents for phage therapy in animals and humans, since their representatives have only lytic life cycle and many of their host bacteria are pathogens. As confirmed experimentally, some Tevenvirinae have non-canonical DNA bases. Non-canonical bases can play an essential role in the diversification of closely related viruses. The article performs a comparative and evolutionary analysis of Tevenvirinae genomes and components of Tevenvirinae genomes. A comparative analysis of these genomes and the genes associated with the synthesis of non-canonical bases allows us to conclude that non-canonical bases have a major influence on the divergence of Tevenvirinae viruses within the same habitats. Supposedly, Tevenvirinae developed a strategy for changing HGT frequency in individual populations, which was based on the accumulation of proteins for the synthesis of non-canonical bases and proteins that used those bases as substrates. Owing to this strategy, ancestors of Tevenvirinae with the highest frequency of HGT acquired genes that allowed them to exist in a certain niche, and ancestors with the lowest HGT frequency preserved the most adaptive of those genes. Given the origin and characteristics of genes associated with the synthesis of non-canonical bases in Tevenvirinae, one can assume that other phages may have similar strategies. The article demonstrates the dependence of genomic diversity of closely related Tevenvirinae on non-canonical bases.

Galaxy and Apollo as a biologist-friendly interface for high-quality cooperative phage genome annotation

In the modern genomic era, scientists without extensive bioinformatic training need to apply high-power computational analyses to critical tasks like phage genome annotation. At the Center for Phage Technology (CPT), we developed a suite of phage-oriented tools housed in open, user-friendly web-based interfaces. A Galaxy platform conducts computationally intensive analyses and Apollo, a collaborative genome annotation editor, visualizes the results of these analyses. The collection includes open source applications such as the BLAST+ suite, InterProScan, and several gene callers, as well as unique tools developed at the CPT that allow maximum user flexibility. We describe in detail programs for finding Shine-Dalgarno sequences, resources used for confident identification of lysis genes such as spanins, and methods used for identifying interrupted genes that contain frameshifts or introns. At the CPT, genome annotation is separated into two robust segments that are facilitated through the automated execution of many tools chained together in an operation called a workflow. First, the structural annotation workflow results in gene and other feature calls. This is followed by a functional annotation workflow that combines sequence comparisons and conserved domain searching, which is contextualized to allow integrated evidence assessment in functional prediction. Finally, we describe a workflow used for comparative genomics. Using this multi-purpose platform enables researchers to easily and accurately annotate an entire phage genome. The portal can be accessed at https://cpt.tamu.edu/galaxy-pub with accompanying user training material.

Recent Insights into Aeromonas salmonicida and Its Bacteriophages in Aquaculture: A Comprehensive Review

The emergence and spread of antimicrobial resistance in pathogenic bacteria of fish and shellfish have caused serious concerns in the aquaculture industry, owing to the potential health risks to humans and animals. Among these bacteria, Aeromonas salmonicida, which is one of the most important primary pathogens in salmonids, is responsible for significant economic losses in the global aquaculture industry, especially in salmonid farming because of its severe infectivity and acquisition of antimicrobial resistance. Therefore, interest in the use of alternative approaches to prevent and control A. salmonicida infections has increased in recent years, and several applications of bacteriophages (phages) have provided promising results. For several decades, A. salmonicida and phages infecting this fish pathogen have been thoroughly investigated in various research areas including aquaculture. The general overview of phage usage to control bacterial diseases in aquaculture, including the general advantages of this strategy, has been clearly described in previous reviews. Therefore, this review specifically focuses on providing insights into the phages infecting A. salmonicida, from basic research to biotechnological application in aquaculture, as well as recent advances in the study of A. salmonicida.

Quantification of T4-Like and T7-Like Cyanophages Using the Polony Method Show They Are Significant Members of the Virioplankton in the North Pacific Subtropical Gyre

The North Pacific Subtropical Gyre (NPSG) is one of the largest biomes on Earth, with the cyanobacterium Prochlorococcus being the most abundant primary producer year-round. Viruses that infect cyanobacteria (cyanophages) influence cyanobacterial mortality, diversity and evolution. Two major cyanophage families are the T4-like cyanomyoviruses and T7-like cyanopodoviruses, yet their abundances and distribution patterns remain unknown due to difficulty in quantifying their populations. To address this limitation, we previously adapted the polony method (for PCR colony) to quantify T7-like cyanophages and applied it to spring populations in the Red Sea. Here, we further adapted the method for the quantification of T4-like cyanophages and analyzed the abundances of T4-like and T7-like cyanophage populations in the photic zone of the NPSG in summer 2015 and spring 2016. Combined, the peak abundances of these two cyanophage families reached 2.8 × 10⁶ and 1.1 × 10⁶ cyanophages ⋅ ml^–1 in the summer and spring, respectively. They constituted between 3 and 16% of total virus-like particles (VLPs), comprising a substantial component of the virioplankton in the NPSG. While both cyanophage families were highly abundant, the T4-like cyanophages were generally 1.3–4.4 fold more so. In summer, cyanophages had similar and reproducible distribution patterns with depth. Abundances were relatively low in the upper mixed layer and increased to form a pronounced subsurface peak at 100 m (1.9 × 10⁶ and 9.1 × 10⁵ phages ⋅ ml^–1 for the T4-like and T7-like cyanophages, respectively), coincident with the maximum in Prochlorococcus populations. Less vertical structure in cyanophage abundances was apparent in the spring profile, despite a subsurface peak in Prochlorococcus numbers. In the summer upper mixed layer, cyanophages constituted a smaller proportion of VLPs than below it and cyanophage to cyanobacteria ratios were considerably lower (1.3–2.8) than those of VLPs to bacteria (8.1–21.2). Differences in abundances between the two families and their contribution to VLPs with depth suggest differences in cyanophage production and/or decay processes relative to other members of the virioplankton in the upper mixed layer. These findings highlight the importance of quantifying distinct populations within the virioplankton to gain accurate understanding of their distribution patterns.

Characterizing the Biology of Lytic Bacteriophage vB_EaeM_φEap-3 Infecting Multidrug-Resistant Enterobacter aerogenes

Carbapenem-resistant Enterobacter aerogenes strains are a major clinical problem because of the lack of effective alternative antibiotics. However, viruses that lyze bacteria, called bacteriophages, have potential therapeutic applications in the control of antibiotic-resistant bacteria. In the present study, a lytic bacteriophage specific for E. aerogenes isolates, designated vB_EaeM_φEap-3, was characterized. Based on transmission electron microscopy analysis, phage vB_EaeM_φEap-3 was classified as a member of the family Myoviridae (order, Caudovirales). Host range determination revealed that vB_EaeM_φEap-3 lyzed 18 of the 28 E. aerogenes strains tested, while a one-step growth curve showed a short latent period and a moderate burst size. The stability of vB_EaeM_φEap-3 at various temperatures and pH levels was also examined. Genomic sequencing and bioinformatics analysis revealed that vB_EaeM_φEap-3 has a 175,814-bp double-stranded DNA genome that does not contain any genes considered undesirable for the development of therapeutics (e.g., antibiotic resistance genes, toxin-encoding genes, integrase). The phage genome contained 278 putative protein-coding genes and one tRNA gene, tRNA-Met (AUG). Phylogenetic analysis based on large terminase subunit and major capsid protein sequences suggested that vB_EaeM_φEap-3 belongs to novel genus “Kp15 virus” within the T4-like virus subfamily. Based on host range, genomic, and physiological parameters, we propose that phage vB_EaeM_φEap-3 is a suitable candidate for phage therapy applications.

Enhanced antibacterial effect of the novel T4-like bacteriophage KARL-1 in combination with antibiotics against multi-drug resistant Acinetobacter baumannii

The continuing rise of infections caused by multi-drug resistant bacteria has led to a renewed interest in bacteriophage therapy. Here we characterize phage vB_AbaM-KARL-1 with lytic activity against multi-drug resistant clinical isolates of Acinetobacter baumannii (AB). Besides genomic and phenotypic phage analysis, the objective of our study was to investigate the antibacterial outcome when the phage acts in concert with distinct antibiotics. KARL-1 belongs to the family of Myoviridae and is able to lyse 8 of 20 (40%) tested clinical isolates. Its double-stranded DNA genome consists of 166,560 bp encoding for 253 open reading frames. Genome wide comparison suggests that KARL-1 is a novel species within the subfamily Tevenvirinae, sharing 77% nucleotide identity (coverage 58%) with phage ZZ1. The antibacterial efficacy at various multiplicities of infection (MOI) was monitored either alone or in combination with meropenem, ciprofloxacin, and colistin. A complete clearance of liquid cultures was achieved with KARL-1 at an MOI of 10^-1 and meropenem (>128 mg/l). KARL-1 was still effective at an MOI of 10^-7, but antibacterial activity was significantly augmented with meropenem. While ciprofloxacin did generally not support phage activity, the application of KARL-1 at an MOI of 10^-7 and therapeutic doses of colistin significantly elevated bacterial suppression. Hence, KARL-1 represents a novel candidate for use against multi-drug resistant AB and the therapeutic outcome may be positively influenced by the addition of traditional antibiotics.

Environmental T4-Family Bacteriophages Evolve to Escape Abortive Infection via Multiple Routes in a Bacterial Host Employing "Altruistic Suicide" through Type III Toxin-Antitoxin Systems

Abortive infection is an anti-phage mechanism employed by a bacterium to initiate its own death upon phage infection. This reduces, or eliminates, production of viral progeny and protects clonal siblings in the bacterial population by an act akin to an "altruistic suicide." Abortive infection can be mediated by a Type III toxin-antitoxin system called ToxINPa consisting of an endoribonuclease toxin and RNA antitoxin. ToxINPa is a heterohexameric quaternary complex in which pseudoknotted RNA inhibits the toxicity of the toxin until infection by certain phages causes destabilization of ToxINPa, leading to bacteriostasis and, eventually, lethality. However, it is still unknown why only certain phages are able to activate ToxINPa. To try to address this issue we first introduced ToxINPa into the Gram-negative enterobacterium, Serratia sp. ATCC 39006 (S 39006) and then isolated new environmental S 39006 phages that were scored for activation of ToxINPa and abortive infection capacity. We isolated three T4-like phages from a sewage treatment outflow point into the River Cam, each phage being isolated at least a year apart. These phages were susceptible to ToxINPa-mediated abortive infection but produced spontaneous "escape" mutants that were insensitive to ToxINPa. Analysis of these resistant mutants revealed three different routes of escaping ToxINPa, namely by mutating asiA (the product of which is a phage transcriptional co-activator); by mutating a conserved, yet functionally unknown, orf84; or by deleting a 6.5-10 kb region of the phage genome. Analysis of these evolved escape mutants may help uncover the nature of the corresponding phage product(s) involved in activation of ToxINPa.

Jumbo Bacteriophages: An Overview

Tailed bacteriophages with genomes larger than 200 kbp are classified as Jumbo phages, and are rarely isolated by conventional methods. These phages are designated “jumbo” owing to their most notable features of a large phage virion and large genome size. However, in addition to these, jumbo phages also exhibit several novel characteristics that have not been observed for phages with smaller genomes, which differentiate jumbo phages in terms of genome organization, virion structure, progeny propagation, and evolution. In this review, we summarize available reports on jumbo phages and discuss the differences between jumbo phages and small-genome phages. We also discuss data suggesting that jumbo phages might have evolved from phages with smaller genomes by acquiring additional functional genes, and that these additional genes reduce the dependence of the jumbo phages on the host bacteria.

Comparative Functional Genomic Analysis of Two Vibrio Phages Reveals Complex Metabolic Interactions with the Host Cell

Sequencing and annotation was performed for two large double stranded DNA bacteriophages, φGrn1 and φSt2 of the Myoviridae family, considered to be of great interest for phage therapy against Vibrios in aquaculture live feeds. In addition, phage–host metabolic interactions and exploitation was studied by transcript profiling of selected viral and host genes. Comparative genomic analysis with other large Vibrio phages was also performed to establish the presence and location of homing endonucleases highlighting distinct features for both phages. Phylogenetic analysis revealed that they belong to the “schizoT4like” clade. Although many reports of newly sequenced viruses have provided a large set of information, basic research related to the shift of the bacterial metabolism during infection remains stagnant. The function of many viral protein products in the process of infection is still unknown. Genome annotation identified the presence of several viral open reading frames (ORFs) participating in metabolism, including a Sir2/cobB (sirtuin) protein and a number of genes involved in auxiliary NAD⁺ and nucleotide biosynthesis, necessary for phage DNA replication. Key genes were subsequently selected for detail study of their expression levels during infection. This work suggests a complex metabolic interaction and exploitation of the host metabolic pathways and biochemical processes, including a possible post-translational protein modification, by the virus during infection.

Bacteriophage application to control the contaminated water with Shigella

Shigella is one of the most important waterborne and foodborne pathogens around the world. Emergence of antibiotic-resistant Shigella has made the development of alternatives to conventional antibiotics necessary. In this study, a virulent Myoviridae bacteriophage, pSs-1 was isolated from environmental water in South Korea and showed infectivity to S. flexneri as well as S. sonnei strains. One-step growth analysis showed that pSs-1 has a short latent period (25 min) and a large burst size (97 PFU/cell). According to the genomic analysis, pSs-1 contains 164,999 bp of genome with a G + C content of 35.54% and it is considered as a member of the T4-like bacteriophage group. These results showed that pSs-1 may have potential as a biocontrol agent instead of conventional antibiotics for shigellosis.

Isolation and Characterization of an Aeromonas punctata Bacteriophage

An Aeromonas punctata bacteriophage, named as DH1, was isolated from East Lake, Wuhan city, China. Morphologically, phage DH1 showed a typical Myoviridae structure consisting of an isometric head (50 nm in diameter) and a visible tail. The bacteriophage had a latent period of about 90 minutes and an average burst size of about 125 PFU•Cell^-1. Restriction enzyme pattern of the bacteriophage’s genome showed that the genome is a double-stranded DNA and about 34kb in size. The sequenced genomic fragments showed highly similarities to gp04 and gp16 sequence of other Myoviridae bacteriophages at protein level.

The Caulobacter crescentus transducing phage Cr30 is a unique member of the T4-like family of myophages

Bacteriophage Cr30 has proven useful for the transduction of Caulobacter crescentus. Nucleotide sequencing of Cr30 DNA revealed that the Cr30 genome consists of 155,997 bp of DNA that codes for 287 proteins and five tRNAs. In contrast to the 67 % GC content of the host genome, the GC content of the Cr30 genome is only 38 %. This lower GC content causes both the codon usage pattern and the amino acid composition of the Cr30 proteins to be quite different from those of the host bacteria. As a consequence, the Cr30 mRNAs probably are translated at a rate that is slower than the normal rate for host mRNAs. A phylogenetic comparison of the genome indicates that Cr30 is a member of the T4-like family that is most closely related to a new group of T-like phages exemplified by фM12.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

Genome organisation of the Acinetobacter lytic phage ZZ1 and comparison with other T4-like Acinetobacter phages

Background

Phage ZZ1, which efficiently infects pathogenic Acinetobacter baumannii strains, is the fifth completely sequenced T4-like Acinetobacter phage to date. To gain a better understanding of the genetic characteristics of ZZ1, bioinformatics and comparative genomic analyses of the T4 phages were performed.

Results

The 166,687-bp double-stranded DNA genome of ZZ1 has the lowest GC content (34.4%) of the sequenced T4-like Acinetobacter phages. A total of 256 protein-coding genes and 8 tRNA genes were predicted. Forty-three percent of the predicted ZZ1 proteins share up to 73% amino acid identity with T4 proteins, and the homologous genes generally retained the same order and transcriptional direction. Beyond the conserved structural and DNA replication modules, T4 and ZZ1 have diverged substantially by the acquisition and deletion of large blocks of unrelated genes, especially in the first halves of their genomes. In addition, ZZ1 and the four other T4-like Acinetobacter phage genomes (Acj9, Acj61, 133, and Ac42) share a well-organised and highly conserved core genome, particularly in the regions encoding DNA replication and virion structural proteins. Of the ZZ1 proteins, 70, 64, 61, and 56% share up to 86, 85, 81, and 83% amino acid identity with Acj9, Acj61, 133, and Ac42 proteins, respectively. ZZ1 has a different number and types of tRNAs than the other 4 Acinetobacter phages, although some of the ZZ1-encoded tRNAs share high sequence similarity with the tRNAs from these phages. Over half of ZZ1-encoded tRNAs (5 out of 8) are related to optimal codon usage for ZZ1 proteins. However, this correlation was not present in any of the other 4 Acinetobacter phages.

Conclusions

The comparative genomic analysis of these phages provided some new insights into the evolution and diversity of Acinetobacter phages, which might elucidate the evolutionary origin and host-specific adaptation of these phages.

Electronic supplementary material

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

Characterisation of host growth after infection with a broad-range freshwater cyanopodophage

Freshwater cyanophages are poorly characterised in comparison to their marine counterparts, however, the level of genetic diversity that exists in freshwater cyanophage communities is likely to exceed that found in marine environments, due to the habitat heterogeneity within freshwater systems. Many cyanophages are specialists, infecting a single host species or strain; however, some are less fastidious and infect a number of different host genotypes within the same species or even hosts from different genera. Few instances of host growth characterisation after infection by broad host-range phages have been described. Here we provide an initial characterisation of interactions between a cyanophage isolated from a freshwater fishing lake in the south of England and its hosts. Designated ΦMHI42, the phage is able to infect isolates from two genera of freshwater cyanobacteria, Planktothrix and Microcystis. Transmission Electron Microscopy and Atomic Force Microscopy indicate that ΦMHI42 is a member of the Podoviridae, albeit with a larger than expected capsid. The kinetics of host growth after infection with ΦMHI42 differed across host genera, species and strains in a way that was not related to the growth rate of the uninfected host. To our knowledge, this is the first characterisation of the growth of cyanobacteria in the presence of a broad host-range freshwater cyanophage.

Complete Genome Sequence of Escherichia Phage Lw1, a New Member of the RB43 Group of Pseudo T-Even Bacteriophages

RB43-related bacteriophages have a specific genome type that clearly distinguishes them from other T4-like viruses. Here, we present the complete genome sequence of a new virulent phage, Lw1, isolated as an Escherichia coli BL21(DE3) contaminant. Lw1 shares an RB43-like genome organization, but it does not contain putative AP2-domain endonuclease genes.

Klebsiella phage vB_KleM-RaK2 - a giant singleton virus of the family Myoviridae

At 346 kbp in size, the genome of a jumbo bacteriophage vB_KleM-RaK2 (RaK2) is the largest Klebsiella infecting myovirus genome sequenced to date. In total, 272 out of 534 RaK2 ORFs lack detectable database homologues. Based on the similarity to biologically defined proteins and/or MS/MS analysis, 117 of RaK2 ORFs were given a functional annotation, including 28 RaK2 ORFs coding for structural proteins that have no reliable homologues to annotated structural proteins in other organisms. The electron micrographs revealed elaborate spike-like structures on the tail fibers of Rak2, suggesting that this phage is an atypical myovirus. While head and tail proteins of RaK2 are mostly myoviridae-related, the bioinformatics analysis indicate that tail fibers/spikes of this phage are formed from podovirus-like peptides predominantly. Overall, these results provide evidence that bacteriophage RaK2 differs profoundly from previously studied viruses of the Myoviridae family.

Long tail fibres of the novel broad-host-range T-even bacteriophage S16 specifically recognize Salmonella OmpC

We report isolation and characterization of the novel T4-like Salmonella bacteriophage vB_SenM-S16. S16 features a T-even morphology and a highly modified 160 kbp dsDNA genome with 36.9 mol % G+C, containing 269 putative coding sequences and three tRNA genes. S16 is a virulent phage, and exhibits a maximally broad host range within the genus Salmonella, but does not infect other bacteria. Synthesis of functional S16 full-length long tail fibre (LTF) in Escherichia coli was possible by coexpression of gp37 and gp38. Surface plasmon resonance analysis revealed nanomolar equilibrium affinity of the LTF to its receptor on Salmonella cells. We show that OmpC serves as primary binding ligand, and that S16 adsorption can be transferred to E. coli by substitution of ompC with the Salmonella homologue. S16 also infects 'rough' Salmonella strains which are defective in lipopolysaccharide synthesis and/or its carbohydrate substitution, indicating that this interaction does not require an intact LPS structure. Altogether, its virulent nature, broad host range and apparent lack of host DNA transduction render S16 highly suitable for biocontrol of Salmonella in foods and animal production. The S16 LTF represents a highly specific affinity reagent useful for cell decoration and labelling, as well as bacterial immobilization and separation.

Complete genome sequence of Cronobacter sakazakii bacteriophage vB_CsaM_GAP161

Cronobacter sakazakii is an opportunistic pathogen that causes infant meningitis and is often associated with milk-based infant formula. We have fully sequenced the genome of a newly isolated lytic C. sakazakii myovirus, vB_CsaM_GAP161, briefly named GAP161. It consists of 178,193 bp and has a G+C content of 44.5%. A total of 277 genes, including 275 open reading frames and two tRNA-encoding genes, were identified. This phage is closely related to coliphages RB16 and RB43 and Klebsiella pneumoniae phage KP15.

Isolation and characterization of a lytic Myoviridae bacteriophage PAS-1 with broad infectivity in Aeromonas salmonicida

To search for candidate control agents against Aeromonas salmonicida subsp. salmonicida infections in aquaculture, one bacteriophage (phage), designated as PAS-1, was isolated from the sediment samples of the rainbow trout (Oncorhynchus mykiss) culture farm in Korea. The PAS-1 was morphologically classified as Myoviridae and possessed approximately 48 kb of double-strand genomic DNA. The phage showed broad host ranges to other subspecies of A. salmonicida as well as A. salmonicida subsp. salmonicida including antibiotic-resistant strains. Its latent period and burst size were estimated to be approximately 40 min and 116.7 PFU/cell, respectively. Furthermore, genomic and structural proteomic analysis of PAS-1 revealed that the phage was closely related to other Myoviridae phages infecting enterobacteria or Aeromonas species. The bacteriolytic activity of phage PAS-1 was evaluated using three subspecies of A. salmonicida strain at different doses of multiplicity of infection, and the results proved to be efficient for the reduction of bacterial growth. Based on these results, PAS-1 could be considered as a novel Aeromonas phage and might have potentiality to reduce the impacts of A. salmonicida infections in aquaculture.

Prophage carriage and diversity within clinically relevant strains of Clostridium difficile

Prophages are encoded in most genomes of sequenced Clostridium difficile strains. They are key components of the mobile genetic elements and, as such, are likely to influence the biology of their host strains. The majority of these phages are not amenable to propagation, and therefore the development of a molecular marker is a useful tool with which to establish the extent and diversity of C. difficile prophage carriage within clinical strains. To design markers, several candidate genes were analyzed including structural and holin genes. The holin gene is the only gene present in all sequenced phage genomes, conserved at both terminals, with a variable mid-section. This allowed us to design two sets of degenerate PCR primers specific to C. difficile myoviruses and siphoviruses. Subsequent PCR analysis of 16 clinical C. difficile ribotypes showed that 15 of them are myovirus positive, and 2 of them are also siphovirus positive. Antibiotic induction and transmission electron microscope analysis confirmed the molecular prediction of myoviruses and/or siphovirus presence. Phylogenetic analysis of the holin sequences identified three groups of C. difficile phages, two within the myoviruses and a divergent siphovirus group. The marker also produced tight groups within temperate phages that infect other taxa, including Clostridium perfringens, Clostridium botulinum, and Bacillus spp., which suggests the potential application of the holin gene to study prophage carriage in other bacteria. This study reveals the high incidence of prophage carriage in clinically relevant strains of C. difficile and correlates the molecular data to the morphological observation.

Complete genome sequence of bacteriophage phiAS7, a T7-like virus that infects Aeromonas salmonicida subsp. salmonicida

To date, a number of Myoviridae bacteriophages that infect Aeromonadaceae have been identified and characterized. However, the genome sequences of Aeromonas phages that not belong to the Myoviridae have not been investigated yet. Herein, we report the complete genome sequence of Aeromonas phage phiAS7, which belongs to the Podoviridae and infects Aeromonas salmonicida subsp. salmonicida.

Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes

The last two decades have revealed that phages (viruses that infect bacteria) are abundant and play fundamental roles in the Earth System, with the T4-like myoviruses (herein T4-like phages) emerging as a dominant 'signal' in wild populations. Here we examine 27 T4-like phage genomes, with a focus on 17 that infect ocean picocyanobacteria (cyanophages), to evaluate lateral gene transfer (LGT) in this group. First, we establish a reference tree by evaluating concatenated core gene supertrees and whole genome gene content trees. Next, we evaluate what fraction of these 'core genes' shared by all 17 cyanophages appear prone to LGT. Most (47 out of 57 core genes) were vertically transferred as inferred from tree tests and genomic synteny. Of those 10 core genes that failed the tree tests, the bulk (8 of 10) remain syntenic in the genomes with only a few (3 of the 10) having identifiable signatures of mobile elements. Notably, only one of these 10 is shared not only by the 17 cyanophages, but also by all 27 T4-like phages (thymidylate synthase); its evolutionary history suggests cyanophages may be the origin of these genes to Prochlorococcus. Next, we examined intragenic recombination among the core genes and found that it did occur, even among these core genes, but that the rate was significantly higher between closely related phages, perhaps reducing any detectable LGT signal and leading to taxon cohesion. Finally, among 18 auxiliary metabolic genes (AMGs, a.k.a. 'host' genes), we found that half originated from their immediate hosts, in some cases multiple times (e.g. psbA, psbD, pstS), while the remaining have less clear evolutionary origins ranging from cyanobacteria (4 genes) or microbes (5 genes), with particular diversity among viral TalC and Hsp20 sequences. Together, these findings highlight the patterns and limits of vertical evolution, as well as the ecological and evolutionary roles of LGT in shaping T4-like phage genomes.

Complete genome sequence and characterization of a broad-host range T4-like bacteriophage phiAS5 infecting Aeromonas salmonicida subsp. salmonicida

In this study, we report one lytic Myoviridae bacteriophage (phage) infecting Aeromonas salmonicida subsp. salmonicida. The phage (named as phiAS5) was isolated from environmental river waters in Korea, and showed broad infectivity to other bacterial species in the family Aeromonadaceae as well as antibiotic-resistant A. salmonicida subsp. salmonicida strains. The biological properties and complete genome of phiAS5 were simultaneously investigated. The complete genome of phiAS5 composed of linear double-stranded DNA of 225,268 bp with G+C content of 43.0%, and encoded 343 putative ORFs, 69 putative promoters, 33 transcriptional terminator regions and 24 tRNA-encoding genes. A high degree of similarity to other T4-like Aeromonas phage was found in most ORFs of phiAS5. Therefore, the genome of phiAS5 was further compared with T4 phage and the closest relative, Aeromonas phage Aeh1, and the result demonstrated that it could be classified as a new member of the T4-like group. The bacteriolytic activity of phiAS5 against A. salmonicida subsp. salmonicida was evaluated at different doses of multiplicity of infection using one each of virulent strain that possesses the ascV gene and multi-drug resistant strain, and the results proved to be efficient for the reduction of bacterial growth. Based on these results, phiAS5 may have the potential for reducing the impacts of virulent or antibiotic-resistant A. salmonicida subsp. salmonicida in aquaculture and may also advance our understanding of the biodiversity of T4-like Aeromonas phages.

A homing endonuclease and the 50-nt ribosomal bypass sequence of phage T4 constitute a mobile DNA cassette

Since its initial description more than two decades ago, the ribosome bypass (or "hop") sequence of phage T4 stands out as a uniquely extreme example of programmed translational frameshifting. The gene for a DNA topoisomerase subunit of T4 has been split by a 1-kb insertion into two genes that retain topoisomerase function. A second 50-nt insertion, beginning with an in-phase stop codon, is inserted near the start of the newly created downstream gene 60. Instead of terminating at this stop codon, approximately half of the ribosomes skip 50 nucleotides and continue translation in a new reading frame. However, no functions, regulatory or otherwise, have been imputed for the truncated peptide that results from termination at codon 46 or for the bypass sequence itself. Moreover, how this unusual mRNA organization arose and why it is maintained have never been explained. We show here that a homing endonuclease (MobA) is encoded in the insertion that created gene 60, and the mobA gene together with the bypass sequence constitute a mobile DNA cassette. The bypass sequence provides protection against self-cleavage by the nuclease, whereas the nuclease promotes horizontal spread of the entire cassette to related bacteriophages. Group I introns frequently provide protection against self-cleavage by associated homing endonucleases. We present a scenario by which the bypass sequence, which is otherwise a unique genetic element, might have been derived from a degenerate group I intron.

Characterization of a ViI-like phage specific to Escherichia coli O157:H7

Phage vB_EcoM_CBA120 (CBA120), isolated against Escherichia coli O157:H7 from a cattle feedlot, is morphologically very similar to the classic phage ViI of Salmonella enterica serovar Typhi. Until recently, little was known genetically or physiologically about the ViI-like phages, and none targeting E. coli have been described in the literature. The genome of CBA120 has been fully sequenced and is highly similar to those of both ViI and the Shigella phage AG3. The core set of structural and replication-related proteins of CBA120 are homologous to those from T-even phages, but generally are more closely related to those from T4-like phages of Vibrio, Aeromonas and cyanobacteria than those of the Enterobacteriaceae. The baseplate and method of adhesion to the host are, however, very different from those of either T4 or the cyanophages. None of the outer baseplate proteins are conserved. Instead of T4's long and short tail fibers, CBA120, like ViI, encodes tail spikes related to those normally seen on podoviruses. The 158 kb genome, like that of T4, is circularly permuted and terminally redundant, but unlike T4 CBA120 does not substitute hmdCyt for cytosine in its DNA. However, in contrast to other coliphages, CBA120 and related coliphages we have isolated cannot incorporate 3H-thymidine (3H-dThd) into their DNA. Protein sequence comparisons cluster the putative "thymidylate synthase" of CBA120, ViI and AG3 much more closely with those of Delftia phage φW-14, Bacillus subtilis phage SPO1, and Pseudomonas phage YuA, all known to produce and incorporate hydroxymethyluracil (hmdUra).

Genome of low-temperature T4-related bacteriophage vB_EcoM-VR7

The complete genome sequence of the T4-related low-temperature Escherichia coli bacteriophage vB_EcoM-VR7 was determined. The genome sequence is 169,285 bp long, with a G+C content of 40.3%. Overall, 95% of the phage genome is coding. It encodes 293 putative protein-encoding open reading frames (ORFs) and tRNA(Met). More than half (59%) of the genomic DNA lacks significant homology with the DNA of T4, but once translated, 72% of the vB_EcoM-VR7 genome (211 ORFs) encodes protein homologues of T4 genes. Overall, 46 vB_EcoM-VR7 ORFs have no homologues in T4 but are derived from other T4-related phages, nine ORFs show similarities to bacterial or non-T4-related phage genes, and 27 ORFs are unique to vB_EcoM-VR7. This phage lacks several T4 enzymes involved in host DNA degradation; however, there is extensive representation of the DNA replication, recombination and repair enzymes as well as the viral capsid and tail structural genes.

The gp38 adhesins of the T4 superfamily: a complex modular determinant of the phage's host specificity

The tail fiber adhesins are the primary determinants of host range in the T4-type bacteriophages. Among the indispensable virion components, the sequences of the long tail fiber genes and their associated adhesins are among the most variable. The predominant form of the adhesin in the T4-type phages is not even the version of the gene encoded by T4, the archetype of the superfamily, but rather a small unrelated protein (gp38) encoded by closely related phages such as T2 and T6. This gp38 adhesin has a modular design: its N-terminal attachment domain binds at the tip of the tail fiber, whereas the C-terminal specificity domain determines its host receptor affinity. This specificity domain has a series of four hypervariable segments (HVSs) that are separated by a set of highly conserved glycine-rich motifs (GRMs) that apparently form the domain’s conserved structural core. The role of gp38’s various components was examined by a comparative analysis of a large series of gp38 adhesins from T-even superfamily phages with differing host specificities. A deletion analysis revealed that the individual HVSs and GRMs are essential to the T6 adhesin’s function and suggests that these different components all act in synergy to mediate adsorption. The evolutionary advantages of the modular design of the adhesin involving both conserved structural elements and multiple independent and easily interchanged specificity determinants are discussed.

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.

Campylobacter jejuni group III phage CP81 contains many T4-like genes without belonging to the T4-type phage group: implications for the evolution of T4 phages

CP81 is a virulent Campylobacter group III phage whose linear genome comprises 132,454 bp. At the nucleotide level, CP81 differs from other phages. However, a number of its structural and replication/recombination proteins revealed a relationship to the group II Campylobacter phages CP220/CPt10 and to T4-type phages. Unlike the T4-related phages, the CP81 genome does not contain conserved replication and virion modules. Instead, the respective genes are scattered throughout the phage genome. Moreover, most genes for metabolic enzymes of CP220/CPt10 are lacking in CP81. On the other hand, the CP81 genome contains nine similar genes for homing endonucleases which may be involved in the attrition of the conserved gene order for the virion core genes of T4-type phages. The phage apparently possesses an unusual modification of C or G bases. Efficient cleavage of its DNA was only achieved with restriction enzymes recognizing pure A/T sites. Uncommonly, phenol extraction leads to a significant loss of CP81 DNA from the aqueous layer, a property not yet described for other phages belonging to the T4 superfamily.

Genomic analysis of oceanic cyanobacterial myoviruses compared with T4-like myoviruses from diverse hosts and environments

T4-like myoviruses are ubiquitous, and their genes are among the most abundant documented in ocean systems. Here we compare 26 T4-like genomes, including 10 from non-cyanobacterial myoviruses, and 16 from marine cyanobacterial myoviruses (cyanophages) isolated on diverse Prochlorococcus or Synechococcus hosts. A core genome of 38 virion construction and DNA replication genes was observed in all 26 genomes, with 32 and 25 additional genes shared among the non-cyanophage and cyanophage subsets, respectively. These hierarchical cores are highly syntenic across the genomes, and sampled to saturation. The 25 cyanophage core genes include six previously described genes with putative functions (psbA, mazG, phoH, hsp20, hli03, cobS), a hypothetical protein with a potential phytanoyl-CoA dioxygenase domain, two virion structural genes, and 16 hypothetical genes. Beyond previously described cyanophage-encoded photosynthesis and phosphate stress genes, we observed core genes that may play a role in nitrogen metabolism during infection through modulation of 2-oxoglutarate. Patterns among non-core genes that may drive niche diversification revealed that phosphorus-related gene content reflects source waters rather than host strain used for isolation, and that carbon metabolism genes appear associated with putative mobile elements. As well, phages isolated on Synechococcus had higher genome-wide %G+C and often contained different gene subsets (e.g. petE, zwf, gnd, prnA, cpeT) than those isolated on Prochlorococcus. However, no clear diagnostic genes emerged to distinguish these phage groups, suggesting blurred boundaries possibly due to cross-infection. Finally, genome-wide comparisons of both diverse and closely related, co-isolated genomes provide a locus-to-locus variability metric that will prove valuable for interpreting metagenomic data sets.

Initiation of bacteriophage T4 DNA replication and replication fork dynamics: a review in the Virology Journal series on bacteriophage T4 and its relatives

Bacteriophage T4 initiates DNA replication from specialized structures that form in its genome. Immediately after infection, RNA-DNA hybrids (R-loops) occur on (at least some) replication origins, with the annealed RNA serving as a primer for leading-strand synthesis in one direction. As the infection progresses, replication initiation becomes dependent on recombination proteins in a process called recombination-dependent replication (RDR). RDR occurs when the replication machinery is assembled onto D-loop recombination intermediates, and in this case, the invading 3' DNA end is used as a primer for leading strand synthesis. Over the last 15 years, these two modes of T4 DNA replication initiation have been studied in vivo using a variety of approaches, including replication of plasmids with segments of the T4 genome, analysis of replication intermediates by two-dimensional gel electrophoresis, and genomic approaches that measure DNA copy number as the infection progresses. In addition, biochemical approaches have reconstituted replication from origin R-loop structures and have clarified some detailed roles of both replication and recombination proteins in the process of RDR and related pathways. We will also discuss the parallels between T4 DNA replication modes and similar events in cellular and eukaryotic organelle DNA replication, and close with some current questions of interest concerning the mechanisms of replication, recombination and repair in phage T4.

Engineered split in Pfu DNA polymerase fingers domain improves incorporation of nucleotide gamma-phosphate derivative

Using compartmentalized self-replication (CSR), we evolved a version of Pyrococcus furiosus (Pfu) DNA polymerase that tolerates modification of the γ-phosphate of an incoming nucleotide. A Q484R mutation in α-helix P of the fingers domain, coupled with an unintended translational termination-reinitiation (split) near the finger tip, dramatically improve incorporation of a bulky γ-phosphate-O-linker-dabcyl substituent. Whether synthesized by coupled translation from a bicistronic (−1 frameshift) clone, or reconstituted from separately expressed and purified fragments, split Pfu mutant behaves identically to wild-type DNA polymerase with respect to chromatographic behavior, steady-state kinetic parameters (for dCTP), and PCR performance. Although naturally-occurring splits have been identified previously in the finger tip region of T4 gp43 variants, this is the first time a split (in combination with a point mutation) has been shown to broaden substrate utilization. Moreover, this latest example of a split hyperthermophilic archaeal DNA polymerase further illustrates the modular nature of the Family B DNA polymerase structure.

Mobile DNA elements in T4 and related phages

Mobile genetic elements are common inhabitants of virtually every genome where they can exert profound influences on genome structure and function in addition to promoting their own spread within and between genomes. Phage T4 and related phage have long served as a model system for understanding the molecular mechanisms by which a certain class of mobile DNA, homing endonucleases, promote their spread. Homing endonucleases are site-specific DNA endonucleases that initiate mobility by introducing double-strand breaks at defined positions in genomes lacking the endonuclease gene, stimulating repair and recombination pathways that mobilize the endonuclease coding region. In phage T4, homing endonucleases were first discovered as encoded within the self-splicing td, nrdB and nrdD introns of T4. Genomic data has revealed that homing endonucleases are extremely widespread in T-even-like phage, as evidenced by the astounding fact that ~11% of the T4 genome encodes homing endonuclease genes, with most of them located outside of self-splicing introns. Detailed studies of the mobile td intron and its encoded endonuclease, I-TevI, have laid the foundation for genetic, biochemical and structural aspects that regulate the mobility process, and more recently have provided insights into regulation of homing endonuclease function. Here, we summarize the current state of knowledge regarding T4-encoded homing endonucleases, with particular emphasis on the td/I-TevI model system. We also discuss recent progress in the biology of free-standing endonucleases, and present areas of future research for this fascinating class of mobile genetic elements.

Genomes of the T4-related bacteriophages as windows on microbial genome evolution

The T4-related bacteriophages are a group of bacterial viruses that share morphological similarities and genetic homologies with the well-studied Escherichia coli phage T4, but that diverge from T4 and each other by a number of genetically determined characteristics including the bacterial hosts they infect, the sizes of their linear double-stranded (ds) DNA genomes and the predicted compositions of their proteomes. The genomes of about 40 of these phages have been sequenced and annotated over the last several years and are compared here in the context of the factors that have determined their diversity and the diversity of other microbial genomes in evolution. The genomes of the T4 relatives analyzed so far range in size between ~160,000 and ~250,000 base pairs (bp) and are mosaics of one another, consisting of clusters of homology between them that are interspersed with segments that vary considerably in genetic composition between the different phage lineages. Based on the known biological and biochemical properties of phage T4 and the proteins encoded by the T4 genome, the T4 relatives reviewed here are predicted to share a genetic core, or "Core Genome" that determines the structural design of their dsDNA chromosomes, their distinctive morphology and the process of their assembly into infectious agents (phage morphogenesis). The Core Genome appears to be the most ancient genetic component of this phage group and constitutes a mere 12-15% of the total protein encoding potential of the typical T4-related phage genome. The high degree of genetic heterogeneity that exists outside of this shared core suggests that horizontal DNA transfer involving many genetic sources has played a major role in diversification of the T4-related phages and their spread to a wide spectrum of bacterial species domains in evolution. We discuss some of the factors and pathways that might have shaped the evolution of these phages and point out several parallels between their diversity and the diversity generally observed within all groups of interrelated dsDNA microbial genomes in nature.

Mobile regulatory cassettes mediate modular shuffling in T4-type phage genomes

Coliphage phi1, which was isolated for phage therapy in the Republic of Georgia, is closely related to the T-like myovirus RB49. The ∼275 open reading frames encoded by each phage have an average level of amino acid identity of 95.8%. RB49 lacks 7 phi1 genes while 10 phi1 genes are missing from RB49. Most of these unique genes encode functions without known homologs. Many of the insertion, deletion, and replacement events that distinguish the two phages are in the hyperplastic regions (HPRs) of their genomes. The HPRs are rich in both nonessential genes and small regulatory cassettes (promoter_early stem-loops [PeSLs]) composed of strong σ⁷⁰-like promoters and stem-loop structures, which are effective transcription terminators. Modular shuffling mediated by recombination between PeSLs has caused much of the sequence divergence between RB49 and phi1. We show that exchanges between nearby PeSLs can also create small circular DNAs that are apparently encapsidated by the virus. Such PeSL “mini-circles” may be important vectors for horizontal gene transfer.

Genetic insertions and diversification of the PolB-type DNA polymerase (gp43) of T4-related phages

In Escherichia coli phage T4 and many of its phylogenetic relatives, gene 43 consists of a single cistron that encodes a PolB family (PolB-type) DNA polymerase. We describe the divergence of this phage gene and its protein product (gp43) (gene product 43) among 26 phylogenetic relatives of T4 and discuss our observations in the context of diversity among the widely distributed PolB enzymes in nature. In two T4 relatives that grow in Aeromonas salmonicida phages 44RR and 25, gene 43 is fragmented by different combinations of three distinct types of DNA insertion elements: (a) a short intercistronic untranslated sequence (IC-UTS) that splits the polymerase gene into two cistrons, 43A and 43B, corresponding to N-terminal (gp43A) and C-terminal (gp43B) protein products; (b) a freestanding homing endonuclease gene (HEG) inserted between the IC-UTS and the 43B cistron; and (c) a group I intron in the 43B cistron. Phage 25 has all three elements, whereas phage 44RR has only the IC-UTS. We present evidence that (a) the split gene of phage 44RR encodes a split DNA polymerase consisting of a complex between gp43A and gp43B subunits; (b) the putative HEG encodes a double-stranded DNA endonuclease that specifically cleaves intron-free homologues of the intron-bearing 43B site; and (c) the group I intron is a self-splicing RNA. Our results suggest that some freestanding HEGs can mediate the homing of introns that do not encode their own homing enzymes. The results also suggest that different insertion elements can converge on a polB gene and evolve into a single integrated system for lateral transfer of polB genetic material. We discuss the possible pathways for the importation of such insertion elements into the genomes of T4-related phages.

Classification of Myoviridae bacteriophages using protein sequence similarity

Background

We advocate unifying classical and genomic classification of bacteriophages by integration of proteomic data and physicochemical parameters. Our previous application of this approach to the entirely sequenced members of the Podoviridae fully supported the current phage classification of the International Committee on Taxonomy of Viruses (ICTV). It appears that horizontal gene transfer generally does not totally obliterate evolutionary relationships between phages.

Results

CoreGenes/CoreExtractor proteome comparison techniques applied to 102 Myoviridae suggest the establishment of three subfamilies (Peduovirinae, Teequatrovirinae, the Spounavirinae) and eight new independent genera (Bcep781, BcepMu, FelixO1, HAP1, Bzx1, PB1, phiCD119, and phiKZ-like viruses). The Peduovirinae subfamily, derived from the P2-related phages, is composed of two distinct genera: the "P2-like viruses", and the "HP1-like viruses". At present, the more complex Teequatrovirinae subfamily has two genera, the "T4-like" and "KVP40-like viruses". In the genus "T4-like viruses" proper, four groups sharing >70% proteins are distinguished: T4-type, 44RR-type, RB43-type, and RB49-type viruses. The Spounavirinae contain the "SPO1-"and "Twort-like viruses."

Conclusion

The hierarchical clustering of these groupings provide biologically significant subdivisions, which are consistent with our previous analysis of the Podoviridae.

Non-canonical RNA arrangement in T4-even phages: accommodated ribosome binding site at the gene 26-25 intercistronic junction

Translational initiation region of bacteriophage T4 gene 25 contains three potential Shine and Dalgarno sequences: SD1, SD2 and SD3. Mutational analysis has predicted that an mRNA stem-loop structure may include SD1 and SD2, bringing the most typical sequence SD3, GAGG, to the initiation codon. Here, we report physical evidence demonstrating that previously predicted mRNA stem-loop structure indeed exists in vivo during gene 25 expression in T4-infected Escherichia coli cells. The second mRNA stem-loop structure is identified 14 nucleotides upstream of the stem-loop I, while the SD3 sequence, as well as the start codon of the gene, are proved to be within an unfolded stretch of mRNA. Phylogenetic comparison of 38 T4-like phages reveals that the T-even and some pseudoT-even phages evolve a similar structural strategy for the translation initiation of 25, while pseudoT-even, schizoT-even and exoT-even phages use an alternative mRNA arrangement. Taken together, the results indicate that a specific mRNA fold forms the split ribosome binding site at the gene 26-25 intercistronic junction, which is highly competent in the translational initiation. We conclude that this ribosome binding site has evolved after T-even diverged from other T4-like phages. Additionally, we determine that the SD sequence GAGG is most widespread in T4.