Major spontaneous genomic rearrangements in Haemophilus influenzae S2 and HP1c1 bacteriophages

Physical maps constructed by the localization of the cleavage site of several restriction endonucleases have shown that the genomes of the Haemophilus bacteriophages S2 and HP1c1 exist in variant forms which differ in the molecular organization of the genomes. At least three regions of different organization of the bacteriophage chromosomes have been identified. The different types of molecular organization can be detected both in the DNA isolated from the mature phage particles and after integration of the phage DNA into the bacterial chromosome.     

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.     

Recombination between bacteriophage T4 and plasmid pBR322 molecules containing cloned T4 DNA

Reciprocal recombination between T4 DNA cloned in plasmid pBR322 and homologous sequences in bacteriophage T4 genomes leads to integration of complete plasmid molecules into phage genomes. Indirect evidence of this integration comes from two kinds of experiments. Packaging of pBR322 DNA into mature phage particles can be detected by a DNA--DNA hybridization assay only when a T4 restriction fragment is cloned in the plasmid. The density of the pBR322 DNA synthesized after phage infection is also consistent with integration of plasmid vector DNA into vegetative phage genomes. Direct evidence of plasmid integration into phage genomes in the region of DNA homology comes from genetic and biochemical analysis of cytosine-containing DNA isolated from mature phage particles. Agarose gel electrophoresis of restriction endonuclease-digested DNA, followed by Southern blot analysis with nick-translated probes, shows that entire plasmid molecules become integrated into phage genomes in the region of T4 DNA homology. In addition, this analysis shows that genomes containing multiple copies of complete plasmid molecules are also formed. Among phage particles containing at least one integrated copy, the average number of integrated plasmid molecules is almost ten. A cloning experiment done with restricted DNA confirms these conclusions and illustrates a method for walking along the T4 genome.     

Genetic Variation of Lactobacillus delbrueckii subsp. lactis Bacteriophages Isolated from Cheese Processing Plants in Finland

The genomes of four Lactobacillus delbrueckii subsp. lactis bacteriophages were characterized by restriction endonuclease mapping, Southern hybridization, and heteroduplex analysis. The phages were isolated from different cheese processing plants in Finland between 1950 and 1972. All four phages had a small isometric head and a long noncontractile tail. Two different types of genome (double-stranded DNA) organization existed among the different phages, the pac type and the cos type, corresponding to alternative types of phage DNA packaging. Three phages belonged to the pac type, and a fourth was a cos-type phage. The pac-type phages were genetically closely related. In the genomes of the pac-type phages, three putative insertion/deletions (0.7 to 0.8 kb, 1.0 kb, and 1.5 kb) and one other region (0.9 kb) containing clustered base substitutions were discovered and localized. At the phenotype level, three main differences were observed among the pac-type phages. These concerned two minor structural proteins and the efficiency of phage DNA packaging. The genomes of the pac-type phages showed only weak homology with that of the cos-type phage. Phage-related DNA, probably a defective prophage, was located in the chromosome of the host strain sensitive to the cos-type phage. This DNA exhibited homology under stringent conditions to the pac-type phages.     

Analysis of High-Throughput Sequencing and Annotation Strategies for Phage Genomes

Background

Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage.

Methodology/Principal Findings

To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling.

Conclusions/Significance

These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.     

Filamentous Bacteriophages of Vibrio parahaemolyticus as a Possible Clue to Genetic Transmission

We have previously reported the isolation and characterization of two filamentous bacteriophages of Vibrio parahaemolyticus, designated Vf12 and Vf33. In this study, to understand the potential of these phages as tools for genetic transmission, we investigated the gene structures of replicative-form (RF) DNAs of their genomes and the distribution of these DNAs on chromosomal and extrachromosomal DNAs. The 7,965-bp nucleotide sequences of Vf12 and Vf33 were determined. An analysis of the overall gene structures revealed that Vf12 and Vf33 had conserved regions and distinctive regions. The gene organization of their conserved regions was similar to that of CTX phage of Vibrio cholerae and coliphage Ff of Escherichia coli, while their distinctive regions were characteristic of Vf12 and Vf33 phage genomes. Southern blot hybridization testing revealed that the filamentous phage genomes integrated into chromosomal DNA of V. parahaemolyticus at the distinctive region of the phage genome and were also distributed on some plasmids of V. parahaemolyticus and total cellular DNAs of one Vibrio damsela and one nonagglutinable Vibrio strain tested. These results strongly suggest the possibilities of genetic interaction among the bacteriophage Vf12 and Vf33 genomes and chromosomal and plasmid-borne DNAs of V. parahaemolyticus strains and of genetic transmission among strains through these filamentous phages.     

Integration of phage Mu into the RP4::D3112A-plasmid in Escherichia coli cells

Hybrid plasmids obtained as a result of Mu phage insertions into the RP4::D3112 plasmid in Escherichia coli cells were studied. Stable maintenance of RP4::D3112 plasmid in E. coli cells was provided by using the D3112 phage genome with a point polar mutation in the A gene which prevented early genes' expression. The presence of D3112A- in the RP4 plasmid has been shown to have no effect on efficiency of phage Mu transposition into this plasmid. Moreover, RP4 and D3112 genomes were equivalent targets for Mu integration. The integration of transposable phage into genome of nonrelated phage can be used as one of the approaches to construct recombinant phage genomes in vivo in the absence of DNA homology.     

Comparison of the physical maps and redundant ends of the chromosomes of phages 2C, SP01, SP82 and phi e

The physical map of 2C DNA (cf. following paper in this journal) was compared to the maps of SP01, SP82 and phi e (three other Bacillus subtilis phages containing hydroxymethyluracil in place of thymine in their DNA). The overall organization of the four genomes was remarkably similar, as indicated by the topology of HaeIII and SalI cleavage segments. The proof was gathered for the presence in the four phage DNAs of large redundant ends carrying a single HaeIII recognition site. The location of the latter proved identical for 2C and SP01, but was shifted in the DNAs of SP82 and phi e. Since the redundant end components of these hydroxymethyluracil genomes are colinear, as shown by cross-hybridization studies, the shifting of the HaeIII cleavage site is presumably due to two base substitutions, suppressing an endonuclease recognition site and establishing a new site elsewhere. Relatedness between the genomes of this family of viruses was evaluated from the fraction of conserved restriction fragments. According to these calculations, 6% base substitutions have occurred within the four viral DNAs, in the course of evolution. However, specific segments of 2C DNA were not present in SP01 and phi e DNA, as shown by cross-hybridization with restriction fragments. These data indicate the occurrence of deletions, in addition to base substitutions, as evolutionary mechanisms prevailing in the genomes of this family of phages.     

Protein-primed DNA replication: role of inverted terminal repeats in the Escherichia coli bacteriophage PRD1 life cycle.

Escherichia coli bacteriophage PRD1 and its relatives contain linear double-stranded DNA genomes, the replication of which proceeds via a protein-primed mechanism. Characteristically, these molecules contain 5'-covalently bound terminal proteins and inverted terminal nucleotide sequences (inverted terminal repeats [ITRs]). The ITRs of each PRD1 phage species have evolved in parallel, suggesting communication between the molecule ends during the life cycle of these viruses. This process was studied by constructing chimeric PRD1 phage DNA molecules with dissimilar end sequences. These molecules were created by combining two closely related phage genomes (i) in vivo by homologous recombination and (ii) in vitro by ligation of appropriate DNA restriction fragments. The fate of the ITRs after propagation of single genomes was monitored by DNA sequence analysis. Recombinants created in vivo showed that phages with nonidentical genome termini are viable and relatively stable, and hybrid phages made in vitro verified this observation. However, genomes in which the dissimilar DNA termini had regained identical sequences were also detected. These observations are explained by a DNA replication model involving two not mutually exclusive pathways. The generality of this model in protein-primed DNA replication is discussed.     

Genomic organization of the related Bacillus subtilis bacteriophages SPP1, 41c, rho 15, and SF6.

The genomes of the related virulent Bacillus subtilis bacteriophages SPP1, 41c, rho 15, and SF6 are partially circularly permuted and terminally redundant. Heteroduplex molecules were produced with various combinations of these DNAs. Their electron-microscopic analyses showed a consistent pattern of homologous and heterologous regions of DNA. Restriction maps of the phage DNAs were established. A comparison of these maps showed a pattern of conserved and variable DNAs compatible with the electron-microscopic analyses. In all phage genomes, regions specifying early and late functions were conserved. In each phage genome, such regions were separated by short segments of heterologous DNA characteristic for each phage.     

Evaluation of a transposase protocol for rapid generation of shotgun high-throughput sequencing libraries from nanogram quantities of DNA

Construction of DNA fragment libraries for next-generation sequencing can prove challenging, especially for samples with low DNA yield. Protocols devised to circumvent the problems associated with low starting quantities of DNA can result in amplification biases that skew the distribution of genomes in metagenomic data. Moreover, sample throughput can be slow, as current library construction techniques are time-consuming. This study evaluated Nextera, a new transposon-based method that is designed for quick production of DNA fragment libraries from a small quantity of DNA. The sequence read distribution across nine phage genomes in a mock viral assemblage met predictions for six of the least-abundant phages; however, the rank order of the most abundant phages differed slightly from predictions. De novo genome assemblies from Nextera libraries provided long contigs spanning over half of the phage genome; in four cases where full-length genome sequences were available for comparison, consensus sequences were found to match over 99% of the genome with near-perfect identity. Analysis of areas of low and high sequence coverage within phage genomes indicated that GC content may influence coverage of sequences from Nextera libraries. Comparisons of phage genomes prepared using both Nextera and a standard 454 FLX Titanium library preparation protocol suggested that the coverage biases according to GC content observed within the Nextera libraries were largely attributable to bias in the Nextera protocol rather than to the 454 sequencing technology. Nevertheless, given suitable sequence coverage, the Nextera protocol produced high-quality data for genomic studies. For metagenomics analyses, effects of GC amplification bias would need to be considered; however, the library preparation standardization that Nextera provides should benefit comparative metagenomic analyses.     

Weigle reactivation of diploid M13 phage

Summary The limited ability of ultraviolet (UV)-irradiated E. coli cells to W-reactivate UV-irradiated, single-stranded DNA phages fd and M13 was investigated. The kinetics of induction for W-reactivation of UV-irradiated fd phage are different from that for other SOS functions. W-reactivation of UV-irradiated M13 phage was studied using phage particles that contain at least two single-stranded DNA genomes. No effect on the extent of W-reactivation of diploid phage was observed, compared to that of normal haploid phage, indicating that the mechanism of W-reactivation of single-stranded DNA phages does not involve recombination between partially replicated genomes.     

Two types of recombination hotspots in bacteriophage T4: one requires DNA damage and a replication origin and the other does not

Recombination hotspots have previously been discovered in bacteriophage T4 by two different approaches, marker rescue recombination from heavily damaged phage genomes and recombination during co-infection by two undamaged phage genomes. The phage replication origin ori(34) is located in a region that has a hotspot in both assays. To determine the relationship between the origin and the two kinds of hotspots, we generated phage carrying point mutations that should inactivate ori(34) but not affect the gene 34 reading frame (within which ori(34) is located). The mutations eliminated the function of the origin, as judged by both autonomous replication of plasmids during T4 infection and two-dimensional gel analysis of phage genomic replication intermediates. As expected from past studies, the ori(34) mutations also eliminated the hotspot for marker rescue recombination from UV-irradiated genomes. However, the origin mutations had no effect on the recombination hotspot that is observed with co-infecting undamaged phage genomes, demonstrating that some DNA sequence other than the origin is responsible for inflated recombination between undamaged genomes. The hotspots for marker rescue recombination may result from a replication fork restart process that acts upon origin-initiated replication forks that become blocked at nearby DNA damage. The two-dimensional gel analysis also revealed phage T4 replication intermediates not previously detected by this method, including origin theta forms.     

Genome characteristics of a novel phage from Bacillus thuringiensis showing high similarity with phage from Bacillus cereus

Bacillus thuringiensis is an important entomopathogenic bacterium belongs to the Bacillus cereus group, which also includes B. anthracis and B. cereus. Several genomes of phages originating from this group had been sequenced, but no genome of Siphoviridae phage from B. thuringiensis has been reported. We recently sequenced and analyzed the genome of a novel phage, BtCS33, from a B. thuringiensis strain, subsp. kurstaki CS33, and compared the gneome of this phage to other phages of the B. cereus group. BtCS33 was the first Siphoviridae phage among the sequenced B. thuringiensis phages. It produced small, turbid plaques on bacterial plates and had a narrow host range. BtCS33 possessed a linear, double-stranded DNA genome of 41,992 bp with 57 putative open reading frames (ORFs). It had a typical genome structure consisting of three modules: the "late" region, the "lysogeny-lysis" region and the "early" region. BtCS33 exhibited high similarity with several phages, B. cereus phage Wβ and some variants of Wβ, in genome organization and the amino acid sequences of structural proteins. There were two ORFs, ORF22 and ORF35, in the genome of BtCS33 that were also found in the genomes of B. cereus phage Wβ and may be involved in regulating sporulation of the host cell. Based on these observations and analysis of phylogenetic trees, we deduced that B. thuringiensis phage BtCS33 and B. cereus phage Wβ may have a common distant ancestor.     

Molecular characterization of new group A streptococcal bacteriophages containing the gene for streptococcal erythrogenic toxin A (speA)

Summary Bacteriophage T12 is the prototype phage carrying the streptococcal erythrogenic toxin A (speA) gene. To examine more closely the phages involved in lysogenic conversion, we examined 300 group A streptococcal strains, and identified and isolated two new phages that carry the speA gene. The molecular sizes of these phage genomes were between 32 and 40 kb, similar to that of phage T12 (35 kb). However, as ascertained by restriction analysis, the physical maps of the new phage genomes were different from phage T12 and from each other. Hybridization analysis also showed that all of these phages were only partially related to one another and the speA gene was always located close to the phage attachment site. Additionally, colony hybridization showed that whereas phage T12 or one of its close relatives is the most common phage associated with the group A streptococci, phage 49 has a much stronger association with the speA gene. A defective phage was also found following pulsed field gel electrophoresis of total phage DNA. This phage appears to be a resident of strain T25₃c and is found only following induction of a T25₃c lysogen. Restriction enzyme analysis of the isolated defective phage DNA suggests that it is the source of the submolar amounts of DNA previously found in association with phage T12 digestion patterns. Additionally, the defective phage may serve as the site of integration of the speA gene-carrying phages described above.     

Viral DNA-synthesizing intermediate complex isolated during assembly of bacteriophage phi X174.

A DNA protein complex that is a precursor of mature phi X174 phage was isolated. The complex sedimented with an S value of 50 in a sucrose gradient and contained phage DNA consisting of a replicative form molecule with an extended tail of single-stranded viral DNA. The viral-strand DNA ranged from one to two genomes in length. Proteins coded on the phi X174 genome as well as the host genome were associated with the viral DNA in the 50S precursor complex. Our results indicated that both viral DNA synthesis and cleavage of the growing viral-strand DNA occurred in the 50S complex.     

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.     

A plasmid cloning system utilizing replication and packaging functions of the filamentous bacteriophage fd

DNA cloning vectors were developed which utilize the replication origin (ori) of bacteriophage fd for their propagation. These vectors depend on the expression of viral gene 2 that was inserted into phage lambda, which in turn was integrated into the host genome. The constitutive expression of gene 2 in the host cells is sufficient for the propagation of at least 100 pfd plasmids per cell. In addition to the fd ori, the pfd vectors carry various antibiotic-resistance genes and unique restriction sites. Some of these vectors have no homologies to commonly used pBR plasmids or to lambda DNA. The nucleotide sequence of the vectors can be deduced from published sequences. Large DNA inserts can be stably propagated in pfd vectors; these are more stable than similar DNA fragments cloned in intact genomes of filamentous bacteriophage. Inclusion of phage sequences required for efficient phage packaging and infection with a helper phage resulted in formation of phage particles containing single-stranded plasmid genomes. Growth at 42 degrees C without selective pressure results in loss of pfd plasmids.     

Ultraviolet reactivation of lambda phage: assay of infectivity of DNA molecules by spheroplast transfection

Prior irradiation of non-lysogenic bacteria by ultraviolet light leads to an increase in the viability of infecting irradiated λ phage (ultraviolet reactivation). Similarly, u.v. irradiation of wild type or uvrD bacteria lysogenic for λcIind^? increased the fraction of closed circular duplex phage DNA molecules formed after infection with u.v.-irradiated λ phage. The closed circular molecules isolated from the irradiated lysogens were shown to be free from u.v. damage by a spheroplast transfection assay. The increase of closed circular molecules is sufficient to explain the ultraviolet reactivation observed by the increase of viability of irradiated phage.In ultraviolet reactivation, damage must be erased on irradiated DNA molecules and the repair is independent of total replication of phage genomes, exchange of sister chromatids or recombination between phage genomes. Protein synthesis is necessary to increase the level of closed circular molecules of irradiated λ phage after irradiation of bacteria.