Comparative whole genome analysis of six diagnostic brucellaphages

Whole genome sequencing of six diagnostic brucellaphages, Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C, was followed with genomic comparisons including recently described genomes of the Tb phage from Mexico (Tb_M) and Pr phage to elucidate genomic diversity and candidate host range determinants. Comparative whole genome analysis revealed high sequence homogeneity among these brucellaphage genomes and resolved three genetic groups consistent with defined host range phenotypes. Group I was composed of Tb and Fz phages that are predominantly lytic for Brucella abortus and Brucella neotomae; Group II included Bk, R/C, and Pr phages that are lytic mainly for B. abortus, Brucella melitensis and Brucella suis; Group III was composed of Wb and S708 phages that are lytic for B. suis, B. abortus and B. neotomae. We found that the putative phage collar protein is a variable locus with features that may be contributing to the host specificities exhibited by different brucellaphage groups. The presence of several candidate host range determinants is illustrated herein for future dissection of the differential host specificity observed among these phages.     

Genomic organization and molecular analysis of virulent bacteriophage 2972 infecting an exopolysaccharide-producing Streptococcus thermophilus strain

The Streptococcus thermophilus virulent pac-type phage 2972 was isolated from a yogurt made in France in 1999. It is a representative of several phages that have emerged with the industrial use of the exopolysaccharide-producing S. thermophilus strain RD534. The genome of phage 2972 has 34,704 bp with an overall G+C content of 40.15%, making it the shortest S. thermophilus phage genome analyzed so far. Forty-four open reading frames (ORFs) encoding putative proteins of 40 or more amino acids were identified, and bioinformatic analyses led to the assignment of putative functions to 23 ORFs. Comparative genomic analysis of phage 2972 with the six other sequenced S. thermophilus phage genomes confirmed that the replication module is conserved and that cos- and pac-type phages have distinct structural and packaging genes. Two group I introns were identified in the genome of 2972. They interrupted the genes coding for the putative endolysin and the terminase large subunit. Phage mRNA splicing was demonstrated for both introns, and the secondary structures were predicted. Eight structural proteins were also identified by N-terminal sequencing and/or matrix-assisted laser desorption ionization-time-of-flight mass spectrometry. Detailed analysis of the putative minor tail proteins ORF19 and ORF21 as well as the putative receptor-binding protein ORF20 showed the following interesting features: (i) ORF19 is a hybrid protein, because it displays significant identity with both pac- and cos-type phages; (ii) ORF20 is unique; and (iii) a protein similar to ORF21 of 2972 was also found in the structure of the cos-type phage DT1, indicating that this structural protein is present in both S. thermophilus phage groups. The implications of these findings for phage classification are discussed.     

High frequency of a novel filamentous phage, VCY φ, within an environmental Vibrio cholerae population

Environmental Vibrio cholerae strains isolated from a coastal brackish pond (Oyster Pond, Woods Hole, MA) carried a novel filamentous phage, VCY, which can exist as a host genome integrative form (IF) and a plasmid-like replicative form (RF). Outside the cell, the phage displays a morphology typical of Inovirus, with filamentous particles ～1.8 μm in length and 7 nm in width. Four independent RF isolates had identical genomes, except for 8 single nucleotide polymorphisms clustered in two regions. The overall genome size is 7,103 bp with 11 putative open reading frames organized into three functional modules (replication, structure and assembly, and regulation). VCY shares sequence similarity with other filamentous phages (including cholera disease-associated CTX) in a highly mosaic manner, indicating evolution by horizontal gene transfer and recombination. VCY integrates in the vicinity of the putative translation initiation factor Sui1 in chromosome II of V. cholerae. A screen of 531 closely related host isolates showed that ～40% harbored phages, with 27% and 13% carrying the IF and RF, respectively. The relative frequencies of the RF and IF differed among strains isolated from the pond or lagoon of Oyster Pond, suggesting that the host habitat influences intracellular phage biology. The overall high prevalence within the host population shows that filamentous phages can be an important component of the environmental biology of V. cholerae.     

Lysis to Kill: Evaluation of the Lytic Abilities,and Genomics of Nine Bacteriophages Infective for Gordonia spp. and Their Potential Use in Activated Sludge Foam Biocontrol

Nine bacteriophages (phages) infective for members of the genus Gordonia were isolated from wastewater and other natural water environments using standard enrichment techniques. The majority were broad host range phages targeting more than one Gordonia species. When their genomes were sequenced, they all emerged as double stranded DNA Siphoviridae phages, ranging from 17,562 to 103,424 bp in size, and containing between 27 and 127 genes, many of which were detailed for the first time. Many of these phage genomes diverged from the expected modular genome architecture of other characterized Siphoviridae phages and contained unusual lysis gene arrangements. Whole genome sequencing also revealed that infection with lytic phages does not appear to prevent spontaneous prophage induction in Gordonia malaquae lysogen strain BEN700. TEM sample preparation techniques were developed to view both attachment and replication stages of phage infection.     

Genomic analysis of bacteriophages SP6 and K1-5, an estranged subgroup of the T7 supergroup

We have determined the genome sequences of two closely related lytic bacteriophages, SP6 and K1-5, which infect Salmonella typhimurium LT2 and Escherichia coli serotypes K1 and K5, respectively. The genome organization of these phages is almost identical with the notable exception of the tail fiber genes that confer the different host specificities. The two phages have diverged extensively at the nucleotide level but they are still more closely related to each other than either is to any other phage currently characterized. The SP6 and K1-5 genomes contain, respectively, 43,769 bp and 44,385 bp, with 174 bp and 234 bp direct terminal repeats. About half of the 105 putative open reading frames in the two genomes combined show no significant similarity to database proteins with a known or predicted function that is obviously beneficial for growth of a bacteriophage. The overall genome organization of SP6 and K1-5 is comparable to that of the T7 group of phages, although the specific order of genes coding for DNA metabolism functions has not been conserved. Low levels of nucleotide similarity between genomes in the T7 and SP6 groups suggest that they diverged a long time ago but, on the basis of this conservation of genome organization, they are expected to have retained similar developmental strategies.     

Spontaneous mutation and recombination among brucellaphages

Merz, George S. (University of Wisconsin, Madison), and Joe B. Wilson. Spontaneous mutation and recombination among brucellaphages. J. Bacteriol. 91:2356-2361. 1966.-Two plaque morphology variants, as seen on Brucella abortus 544A, termed c (for clear plaque) and lc (for late-clearing plaque) were isolated from stocks of wild-type brucellaphage and from colonies of B. abortus 544A which had undergone an alteration in colonial morphology associated with the establishment of the phage carrier state. Single-burst experiments showed that the phage variants arise by spontaneous mutation of the wild-type phage during its replication on B. abortus strain R19. Two-factor crosses of independently occurring c mutant phages showed the presence of wild-type recombinants among the progeny. Control experiments showed that there are no strong selective forces against either wild-type or c mutant phage inherent in the cross-procedure. Other control experiments ruled out the possibility that wild-type phage in the cross-progeny resulted from either back mutation of the c mutants or the presence of wild-type phage among the input c mutants.     

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

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

Antisense RNA targeting of primase interferes with bacteriophage replication in Streptococcus thermophilus

The putative primase gene and other genes associated with the Sfi21-prototype genome replication module are highly conserved in Streptococcus thermophilus bacteriophages. Expression of antisense RNAs complementary to the putative primase gene (pri3.1) from S. thermophilus phage kappa 3 provided significant protection from kappa 3 and two other Sfi21-type phages. Expression of pri3.10-AS, an antisense RNA that covered the entire primase gene, reduced the efficiency of plaquing (EOP) of kappa 3 to 3 x 10(-3) and reduced its burst size by 20%. Mutant phages capable of overcoming antisense inhibition were not recovered. Thirteen primase-specific antisense cassettes of different lengths (478 to 1,512 bp) were systematically designed to target various regions of the gene. Each cassette conferred some effect, reducing the EOP to between 0.8 and 3 x 10(-3). The largest antisense RNAs (1.5 kb) were generally found to confer the greatest reductions in EOP, but shorter (0.5 kb) antisense RNAs were also effective, especially when directed to the 5' region of the gene. The impacts of primase-targeted antisense RNAs on phage development were examined. The expression of pri3.10-AS resulted in reductions in target RNA abundance and the number of phage genomes synthesized. Targeting a key genome replication function with antisense RNA provided effective phage protection in S. thermophilus.     

Bioinformatic analysis of the Acinetobacter baumannii phage AB1 genome

As one of the pathogens of hospital-acquired infections, Acinetobacter baumannii poses great challenges to the public health. A. baumannii phage could be an effective way to fight multi-resistant A. baumannii. Here, we completed the whole genome sequencing of the complete genome of A. baumannii phage AB1, which consists of 45,159bp and is a double-stranded DNA molecule with an average GC content of 37.7%. The genome encodes one tRNA gene and 85 open reading frames (ORFs) and the average size of the ORF is 531bp in length. Among 85 ORFs, only 14 have been identified to share significant sequence similarities to the genes with known functions, while 28 are similar in sequence to the genes with function-unknown genes in the database and 43 ORFs are uniquely present in the phage AB1 genome. Fourteen function-assigned genes with putative functions include five phage structure proteins, an RNA polymerase, a big sub-unit and a small sub-unit of a terminase, a methylase and a recombinase and the proteins involved in DNA replication and so on. Multiple sequence alignment was conducted among those homologous proteins and the phylogenetic trees were reconstructed to analyze the evolutionary courses of these essential genes. From comparative genomics analysis, it turned out clearly that the frame of the phage genome mainly consisted of genes from Xanthomonas phages, Burkholderia ambifaria phages and Enterobacteria phages and while it comprises genes of its host A. baumannii only sporadically. The mosaic feature of the phage genome suggested that the horizontal gene transfer occurred among the phage genomes and between the phages and the host bacterium genomes. Analyzing the genome sequences of the phages should lay sound foundation to investigate how phages adapt to the environment and infect their hosts, and even help to facilitate the development of biological agents to deal with pathogenic bacteria.     

Complete genomic sequence of the virulent Salmonella bacteriophage SP6

We report the complete genome sequence of enterobacteriophage SP6, which infects Salmonella enterica serovar Typhimurium. The genome contains 43,769 bp, including a 174-bp direct terminal repeat. The gene content and organization clearly place SP6 in the coliphage T7 group of phages, but there is approximately 5 kb at the right end of the genome that is not present in other members of the group, and the homologues of T7 genes 1.3 through 3 appear to have undergone an unusual reorganization. Sequence analysis identified 10 putative promoters for the SP6-encoded RNA polymerase and seven putative rho-independent terminators. The terminator following the gene encoding the major capsid subunit has a termination efficiency of about 50% with the SP6-encoded RNA polymerase. Phylogenetic analysis of phages related to SP6 provided clear evidence for horizontal exchange of sequences in the ancestry of these phages and clearly demarcated exchange boundaries; one of the recombination joints lies within the coding region for a phage exonuclease. Bioinformatic analysis of the SP6 sequence strongly suggested that DNA replication occurs in large part through a bidirectional mechanism, possibly with circular intermediates.     

The P-SSP7 cyanophage has a linear genome with direct terminal repeats

P-SSP7 is a T7-like phage that infects the cyanobacterium Prochlorococcus MED4. MED4 is a member of the high-light-adapted Prochlorococcus ecotypes that are abundant in the surface oceans and contribute significantly to primary production. P-SSP7 has become a model system for the investigation of T7-like phages that infect Prochlorococcus. It was classified as T7-like based on genome content and organization. However, because its genome assembled as a circular molecule, it was thought to be circularly permuted and to lack the direct terminal repeats found in other T7-like phages. Here we sequenced the ends of the P-SSP7 genome and found that the genome map is linear and contains a 206 bp repeat at both genome ends. Furthermore, we found that a 728 bp region of the genome originally placed downstream of the last ORF is actually located upstream of the first ORF on the genome map. These findings suggest that P-SSP7 is likely to use the direct terminal repeats for genome replication and packaging in a similar manner to other T7-like phages. Moreover, these results highlight the importance of experimentally verifying the ends of phage genomes, and will facilitate the use of P-SSP7 as a model for the correct assembly and end determination of the many T7-like phages isolated from the marine environment that are currently being sequenced.     

Genome organization of temperate phage 11143 from emetic Bacillus cereus NCTC11143

A temperate phage was isolated from emetic Bacillus cereus NCTC 11143 by mitomycin C and characterized by transmission electron microscopy and DNA and protein analyses. Whole genome sequencing of Bacillus phage 11143 was performed by GS-FLX. The phage has a dsDNA genome of 39,077 bp and a 35% G+C content. Bioinformatic analysis of the phage genome revealed 49 putative ORFs involved in replication, morphogenesis, DNA packaging, lysogeny, and host lysis. Bacillus phage 11143 could be classified as a member of the Siphoviridae family by morphology and genome structure. Genomic comparisons at the DNA and protein levels revealed homologous genetic modules with patterns and morphogenesis proteins similar to those of other Bacillus phages. Thus, Bacillus phages might have a mosaic genetic relationship.     

Manipulating or Superseding Host Recombination Functions: A Dilemma That Shapes Phage Evolvability

Phages, like many parasites, tend to have small genomes and may encode autonomous functions or manipulate those of their hosts''. Recombination functions are essential for phage replication and diversification. They are also nearly ubiquitous in bacteria. The E. coli genome encodes many copies of an octamer (Chi) motif that upon recognition by RecBCD favors repair of double strand breaks by homologous recombination. This might allow self from non-self discrimination because RecBCD degrades DNA lacking Chi. Bacteriophage Lambda, an E. coli parasite, lacks Chi motifs, but escapes degradation by inhibiting RecBCD and encoding its own autonomous recombination machinery. We found that only half of 275 lambdoid genomes encode recombinases, the remaining relying on the host''s machinery. Unexpectedly, we found that some lambdoid phages contain extremely high numbers of Chi motifs concentrated between the phage origin of replication and the packaging site. This suggests a tight association between replication, packaging and RecBCD-mediated recombination in these phages. Indeed, phages lacking recombinases strongly over-represent Chi motifs. Conversely, phages encoding recombinases and inhibiting host recombination machinery select for the absence of Chi motifs. Host and phage recombinases use different mechanisms and the latter are more tolerant to sequence divergence. Accordingly, we show that phages encoding their own recombination machinery have more mosaic genomes resulting from recent recombination events and have more diverse gene repertoires, i.e. larger pan genomes. We discuss the costs and benefits of superseding or manipulating host recombination functions and how this decision shapes phage genome structure and evolvability.     

Genomic characterization of four novel bacteriophages infecting the clinical pathogen Klebsiella pneumoniae

Bacteriophages are an invaluable source of novel genetic diversity. Sequencing of phage genomes can reveal new proteins with potential uses as biotechnological and medical tools, and help unravel the diversity of biological mechanisms employed by phages to take over the host during viral infection. Aiming to expand the available collection of phage genomes, we have isolated, sequenced, and assembled the genome sequences of four phages that infect the clinical pathogen Klebsiella pneumoniae: vB_KpnP_FBKp16, vB_KpnP_FBKp27, vB_KpnM_FBKp34, and Jumbo phage vB_KpnM_FBKp24. The four phages show very low (0–13%) identity to genomic phage sequences deposited in the GenBank database. Three of the four phages encode tRNAs and have a GC content very dissimilar to that of the host. Importantly, the genome sequences of the phages reveal potentially novel DNA packaging mechanisms as well as distinct clades of tubulin spindle and nucleus shell proteins that some phages use to compartmentalize viral replication. Overall, this study contributes to uncovering previously unknown virus diversity, and provides novel candidates for phage therapy applications against antibiotic-resistant K. pneumoniae infections.     

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: C_L1, C_L2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (C_L1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages C_L1, C_L2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages C_L1, C_L2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.