pmc logo imageJournal ListSearchpmc logo image
Logo of viroljBioMed Central web siteReference to the article.Search.Manuscript submission.Registration.Journal front page.
Virol J. 2008; 5: 122.
Published online 2008 October 20. doi: 10.1186/1743-422X-5-122.
PMCID: PMC2579914
The genome and proteome of the Kluyvera bacteriophage Kvp1 – another member of the T7-like Autographivirinae
Erika J Lingohr,1 Andre Villegas,1 Yi-Min She,2 Pieter-Jan Ceyssens,3 and Andrew M Kropinskicorresponding author1,4
1Public Health Agency of Canada, Laboratory for foodborne Zoonoses, Guelph, ON N1G 3W4, Canada
2Department of Chemistry, Queen's University, Kingston, Ontario K7L 3N6, Canada
3Laboratory of Gene Technology, Katholieke Universiteit Leuven, Kasteelpark Arenberg 21, Leuven, B-3001, Belgium
4Department of Molecular and Cellular Biology, University of Guelph, Guelph, Ontario, N1G 2W1, Canada
corresponding authorCorresponding author.
Erika J Lingohr: Erika_Lingohr/at/phac-aspc.gc.ca; Andre Villegas: Andre_Villegas/at/phac-aspc.gc.ca; Yi-Min She: YiMin.She/at/chem.queensu.ca; Pieter-Jan Ceyssens: PieterJan.Ceyssens/at/biw.kuleuven.be; Andrew M Kropinski: Andrew_Kropinski/at/phac-aspc.gc.ca
Received September 18, 2008; Accepted October 20, 2008.
Abstract

Background
Kluyvera, a genus within the family Enterobacteriaceae, is an infrequent cause of human infections. Bacteriophage Kvp1, the only bacteriophage isolated for one of its species, Kluyvera cryocrescens, is a member of the viral family Podoviridae.

Results
The genome of Kvp1, the first Kluyvera cryocrescens-specific bacteriophage, was sequenced using pyrosequencing (454 technology) at the McGill University and Genome Québec Innovation Centre. The two contigs were closed using PCR and the sequence of the terminal repeats completed by primer walking off the phage DNA. The phage structural proteome was investigated by SDS-PAGE and mass spectrometry.

Conclusion
At 39,472 bp, the annotated genome revealed a closer relationship to coliphage T3 than T7 with Kvp1 containing homologs to T3 early proteins S-adenosyl-L-methionine hydrolase (0.3) and protein kinase (0.7). The quantitative nature of the relationships between Kvp1 and the other members of the T7-like virus genus (T7, T3, [var phi]A1122, [var phi]YeO3-12, Berlin, K1F, VP4 and gh-1) was confirmed using CoreGenes.

Background

The T7-like bacterial viruses are members of the Podoviridae – phages with short tails – and are characterized by a simple but elegant temporal transcriptional control system [1]. The early genes are transcribed by the host RNA polymerase while the middle and late regions are transcribed by a single subunit phage-encoded RNA polymerase which recognizes unique 23 bp promoters sequences [2]. These viruses are one of the most common types of bacteriophages with 26–29 defined or tentative species according to the VIII report of the International Committee on the Taxonomy of Viruses [3,4]. Most of the host species are members of the γ-Proteobacteria (Erwinia, Escherichia, Klebsiella, Morganella, Pseudomonas, Salmonella, Vibrio, and Yersinia) but viral isolates also infecting α-Proteobacteria (Caulobacter, and Rhizobium) have been isolated. Fifteen T7-like phages have been sequenced and deposited with GenBank. As a result of a reanalysis, at the protein level, of relationships within the "T7-like viruses" this group of bacteriophages have been classified into the subfamily Autographivirinae which currently possesses three genera: T7-like, Sp6-like and [var phi]KMV-like viruses [5]. Kvp1, the first Kluyvera cryocrescens-specific bacteriophage, was isolated from the Matanza River in Buenos Aires (Argentina) by Gadaleta and Zorzopulos [6]. Morphologically this phage is a member of the Podoviridae. Eleven clones derived from AluI or HaeIII digestion of the viral DNA were sequenced, by these authors, revealing strong sequence similarity to coliphage T7. To further analyze the correct taxonomic position of this virus we have completed the sequence of its genome noting its very close similarity to Yersinia phage Berlin and coliphage T3.

Results and discussion

Pyrosequencing (454 technology) has been used to determine the sequence of the genomes of Bacillus thuringiensis phage 0305[var phi]8-36 [7] and coliphage JK98 [8], and, in this incidence, the genome of Kvp1. Sequencing resulted in 2 contigs with 53-fold coverage. While this type of sequencing can result in potential errors at oligonucleotide runs, none were observed in the data on Kvp1. The gap, representing 0 bp, was closed by PCR amplification and ABI sequencing; while the nature of the termini were verified by primer walking off phage DNA template. The total genome is 39,472 bp with 194 bp terminal direct repeats, and a base composition of 48.6 mol%G+C – characteristics remarkably consistent with other T7-like phages. By comparison, the genomes of T7-like phages range from 37.4 kb (Pseudomonas phage gh-1) to 45.4 kb (Erwinia phage Era103) while the reported terminal repeats range from Yersinia phage [var phi]A1122 at 148 bp to Pseudomonas aeruginosa phage LKD16 at 428 bp.

No tRNA genes were discovered, which was not an unexpected observation since no T7-like phages have been found to harbour them; 46 ORFs were delineated encoding protein products which show the strongest sequence similarity to gene products (Gps) from Yersinia phage Berlin (NC_008694). To investigate the relationships further we employed two homology tools, one of which function at the DNA sequence level (Mauve) and one, CoreGenes, which compares proteins.

Several regions of dissimilarity (indicates by areas of white in Figure 1) centred at genes 1.05, 4.7–2.8, 5.3–5.5, 17–17.2, 18.2 and at the left end of the genome are noted. Several of these genes are not found in phage Berlin. The most interesting difference is in gene 17 which encodes the tail fibre. As with other Gp17 homologs sequence similarity is only found at the N-terminus, the part of the protein which is associated with the tail structure. The C-terminus is involved in ligand interactions and exhibits considerable differences.

Figure 1Figure 1
Comparison of the genomes of Yersinia phage Berlin and Kluyvera phage Kvp1 using Mauve. Underneath the name (Kluyvera phage Kvp1) is ruler in kb, the degree of sequence similarity, indicated by the intensity of the red region, and, the gene map with the (more ...)

Using CoreGenes Kvp1 shares 37 (61.7%), 12 (23.1%) and 9 (18.4%) homologs with the type phages of the three Autographivirinae genera – T7, Sp6 and [var phi]KMV-like viruses. The results indicate that Kvp1 is a member of the T7-like virus genus. A comparison with the proteome of phage Berlin indicates 37 homologs – 82.2% common proteins. While the percentage of common proteins is less when compared with coliphage T3 (70.9% similarity) the early regions of T3 and Kvp1 are very similar in that Kvp1 encodes T3-like Gp0.3, 0.6 and 0.7 homologs. The product of early gene 0.3 (Ocr) is a small protein which mimics B-form DNA and binds to, and inhibits, type I restriction endonucleases [9,10]; and, possesses S-adenosyl-L-methionine hydrolase activity [11]. Gp0.7 produces functions in host gene shutoff [12] and as a protein kinase which phosphorylates host elongation factors G and P and ribosomal protein S6 [13]. A major dissimilarity between Kvp1 and T3/T7 is that while the latter phages possess multiple strong promoters recognized by the host RNA polymerase, only a single promoter showing homology to the consensus was found in the Kvp1 genome. In keeping with the protein similarity to Yersinia phage Berlin, the phage specific promoters are also most closely related in sequence to those of this bacterial virus (Fig. 2).

Figure 2Figure 2
Weblogos[27]of some T7-like phage-specific promoters created online at http://weblogo.berkeley.edu/logo.cgi showing that the Kvp1 promoters are most closely related to those of phage Berlin.

Typical of this type of bacteriophage, Kvp1 displays a simple protein profile (Fig. 3) in which most of the protein bands can be assigned based upon the extensive knowledge of these phages, and the mass of the protein bands compared with the in silico analysis of the proteins based upon genomic analysis.

Figure 3Figure 3
Denaturing SDS-PAGE of bacteriophage Kvp1 structural proteins (LaneB) alongside the protein marker (Lane A). The masses of the proteins are indicated in the adjacent lanes. The tentative identification based on in silico analysis of the properties of (more ...)

One band of interest, with a mass of 43.6 kDa, was noted migrating just above that of the major capsid protein (Gp10) which was not immediately linked, on the basis of its mass, to a product of one of the morphogenesis genes. One of the characteristics of some T7-like phages is that they display, on SDS-PAGE, two "versions" of the major capsid protein which are designated as 10A and 10B [14]. The sequences of the amino termini of these proteins are identical, but during translation a rare ribosomal slippage occurs permitting the elongation of the protein product. The features of this system are a protein slightly larger than the capsid, a slippery site in the DNA/RNA and a downstream stem-loop structure capable of forming a pseudoknot [15,16]. We obtained evidence for a potential pseudoknot using pknotsRG at http://bibiserv.techfak.uni-bielefeld.de/pknotsrg/submission.html[17] located 144 bp downstream from the end of the capsid gene, but no typical slippage site was observed. The nature of this protein was investigated by in-gel enzymatic digestion and high-resolution mass spectrometry. MALDI QqTOF MS analysis on a tryptic digest has yielded 70% sequence coverage of the protein Gp10 (Figure 4), and three unique peptides were present at m/z 2244.123, m/z 2372.219 and m/z 2692.310 which revealed the distinct C-terminal amino acid residues 327–372 from the protein sequences of Gp10. This indicates that Kvp1 produces a major capsid protein (10A) and a minor protein (10B) through programmed -1 frameshifting at TTTTCA. The Gp10B protein is predicted to have a calculated mass of 42.1 kDa, consistent with the estimated value of 43.6 kDa by SDS-PAGE gel.

Figure 4Figure 4
MALDI QqTOF MS and MS/MS analyses on the in-gel tryptic digest of the 43.6 kDa protein band. (A) MS spectrum of the trypsin digested peptides. The corresponding Gp10B peptides are shown in the parenthesis. (B) MS/MS sequencing of a typical peptide at (more ...)
Conclusion

Our data conclusively demonstrate that Kluyvera virus Kvp1 is a member of T7-like virus genus of the Podoviridae subfamily Autographivirinae. It differs from phages such as T3 and [var phi]YeO3-12 which exhibit capsid frameshifting at lysyl residues, by ribosomal slippage at polyU residues (phenylalanine) – a property it shares with Yersinia phage [var phi]A1122.

Materials and methods

Purification of phage wV8
Bacteriophage Kvp1 (HER400) and its host K. cryocrescens strain HER1400 were received from the Felix d'Hérelle Reference Center for Bacterial Viruses at Université Laval (Québec, QC, Canada). The phage was propagated at 30°C using standard protocols, precipitated using polyethylene glycol 8000 and purified through two rounds of CsCl equilibrium gradient centrifugation [18].

DNA sequencing
The DNA was isolated using the SDS-proteinase K protocol of Sambrook and Russell (2001) and was submitted to the McGill University and Génome Québec Innovation Centre (Montréal, QC, Canada) for DNA sequencing. This resulted in two contigs which were closed using PCR with custom primers and, standard dideoxy sequencing of the amplicon (University of Guelph, Laboratory Services, Guelph, ON, Canada). The termini were determined by primer walking.

Genome annotation
The genome was screened for tRNA-encoding genes using Aragorn [19] and tRNAScan [20]; and, for protein encoding genes using Kodon (Applied Maths, Austin, TX) and PSI-BLAST [21]. Rho-independent terminators identified using TransTerm [22] at http://nbc11.biologie.uni-kl.de/framed/left/menu/auto/right/clusterinfo2. Phage-specific promoters were discovered using PHIRE [23] and displayed using WebLogo [27]. The sequence of this bacteriophage has been deposited with GenBank (accession no. FJ194439).

Whole genome comparisons
These were carried out using Mauve [24], and CoreGenes [25].

Proteomics
SDS-PAGE [26] was carried out on CsCl-purified phage particles along with the PageRuler Unstained Protein Ladder (Fermentas, Burlington, ON, Canada) stained with Coomassie brilliant blue R250 and characterized using Bionumerics software (Applied Maths). Bands were further characterized by in situ trypsin digestion and mass spectrometry. Briefly, the excised gel bands were destained until colorless, and dried using a SpeedVac. Following reduction with DTT and alkylation with iodoacetamide, the protein was digested with 10 ng of sequencing grade trypsin (Calbiochem) in 25 mM NH4HCO3 (pH 7.6) at 37°C overnight. The proteolytic peptides were extracted, and cleaned up by a C18 Ziptip (Millipore). MALDI data were acquired using an Applied Biosystems/MDS Sciex QStar XL quadrupole time-of-flight (QqTOF) mass spectrometer under a nitrogen laser (337 nm), and 2,5-dihydroxybenzoic acid was used as the matrix. All peptide fingerprinting masses (m/z) on the MS spectrum were compared with the theoretical values generated in-silico by MS-Digest, a ProteinProspector program developed in the UCSF Mass Spectrometry Facility http://prospector.ucsf.edu/. The individual peptide sequence was identified by MALDI MS/MS measurements on the same instrument using argon as the collision gas.

Abbreviations

MALDI: matrix-assisted laser desorption ionization; QqTOF MS: quadrupole time-of-flight mass spectrometry; MS/MS: tandem mass spectrometry.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AMK planned the experiments and prepared the manuscript, EJL propagated and purified the phage; and together with YS contributed to the proteomics. P-JC sequenced the ends of the genome; and AV contributed to the genome annotation.

Acknowledgements

A.K. is supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. We thank Rob Lavigne for his critical review of the MS. P-J.C. holds a predoctoral fellowship from the Instituut voor de Aanmoediging van Innovatie door Wetenschap en Technologie in Vlaanderen (I.W.T., Belgium).

References
  • Molineux, IJ. The T7 Group. In: Calendar R. , editor. The Bacteriophages. New York, NY: Oxford University Press; 2006. pp. 277–301.
  • Chen, Z; Schneider, TD. Information theory based T7-like promoter models: classification of bacteriophages and differential evolution of promoters and their polymerases. Nucleic Acids Research. 2005;33:6172–6187. doi: 10.1093/nar/gki915. [PubMed]
  • Fauquet, CM; Mayo, MA; Maniloff, J; Desselberger, U; Ball, A. Virus Taxonomy. In: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball A. , editor. VIIIth Report of the International Committee on Taxonomy of Viruses. New York, NY: Elsevier Academic Press; 2005. pp. 35–85.
  • Hendrix, RW; Casjens, SR. Myoviridae, Siphoviridae, Podoviridae. In: Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA. , editor. Virus Taxonomy VIIIth Report of the International Committee on Taxonomy of Viruses. New York: Elsevier Academic Press; 2005. pp. 43–47.
  • Lavigne, R; Seto, D; Mahadevan, P; Ackermann, H-W; Kropinski, AM. Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools. Research in Microbiology. 2008;159:406–414. doi: 10.1016/j.resmic.2008.03.005. [PubMed]
  • Gadaleta, P; Zorzopulos, J. Kluyvera bacteriophage Kvp1: a new member of the Podoviridae family phylogenetically related to the coliphage T7. Virus Research. 1997;51:43–52. doi: 10.1016/S0168-1702(97)00078-6. [PubMed]
  • Thomas, JA; Hardies, SC; Rolando, M; Hayes, SJ; Lieman, K; Carroll, CA, et al. Complete genomic sequence and mass spectrometric analysis of highly diverse, atypical Bacillus thuringiensis phage 0305f8-36. Virology. 2007;368:405–421. doi: 10.1016/j.virol.2007.06.043. [PubMed]
  • Zuber, S; Ngom-Bru, C; Barretto, C; Bruttin, A; Brussow, H; Denou, E. Genome analysis of phage JS98 defines a fourth major subgroup of T4-like phages in Escherichia coli. Journal of Bacteriology. 2007;189:8206–8214. doi: 10.1128/JB.00838-07. [PubMed]
  • Walkinshaw, MD; Taylor, P; Sturrock, SS; Atanasiu, C; Berge, T; Henderson, RM, et al. Structure of Ocr from bacteriophage T7, a protein that mimics B-form DNA. Molecular Cell. 2002;9:187–194. doi: 10.1016/S1097-2765(02)00435-5. [PubMed]
  • Sturrock, SS; Dryden, DT; Atanasiu, C; Dornan, J; Bruce, S; Cronshaw, A, et al. Crystallization and preliminary X-ray analysis of ocr, the product of gene 0.3 of bacteriophage T7. Acta Crystallographica Section D-Biological Crystallography. 2001;57:1652–1654. doi: 10.1107/S0907444901011623.
  • Studier, FW; Movva, NR. SAMase gene of bacteriophage T3 is responsible for overcoming host restriction. J Virol. 1976;19:136–145. [PubMed]
  • Marchand, I; Nicholson, AW; Dreyfus, M. High-level autoenhanced expression of a single-copy gene in Escherichia coli: overproduction of bacteriophage T7 protein kinase directed by T7 late genetic elements. Gene. 2001;262:231–238. doi: 10.1016/S0378-1119(00)00526-6. [PubMed]
  • Robertson, ES; Aggison, LA; Nicholson, AW. Phosphorylation of elongation factor G and ribosomal protein S6 in bacteriophage T7-infected Escherichia coli. Molecular Microbiology. 1994;11:1045–1057. doi: 10.1111/j.1365-2958.1994.tb00382.x. [PubMed]
  • Condron, BG; Atkins, JF; Gesteland, RF. Framshifting in gene 10 of bacteriophage T7. Journal of Bacteriology. 1991;173:6998–7003. [PubMed]
  • Alam, SL; Atkins, JF; Gesteland, RF. Programmed ribosomal frameshifting: much ado about knotting! Proceedings of the National Academy of Sciences of the United States of America. 1999;96:14177–14179. doi: 10.1073/pnas.96.25.14177. [PubMed]
  • Chandler, M; Fayet, O. Translational frameshifting in the control of transposition in bacteria. Molecular Microbiology. 1993;7:497–503. doi: 10.1111/j.1365-2958.1993.tb01140.x. [PubMed]
  • Reeder, J; Giegerich, R. Design implementation and evaluation of a practical pseudoknots folding algorithm based upon thermodynamics. BMC Bioinformatics. 2004;5:104–115. doi: 10.1186/1471-2105-5-104. [PubMed]
  • Sambrook, J; Russell, DW. Molecular Cloning: A Laboratory Manual. Third. Cold Spring Harbor, New York: Cold Spring Harbor Press; 2001.
  • Laslett, D; Canback, B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Research. 2004;32:11–16. doi: 10.1093/nar/gkh152. [PubMed]
  • Lowe, TM; Eddy, SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Research. 1997;25:955–964. doi: 10.1093/nar/25.5.955. [PubMed]
  • Altschul, SF; Madden, TL; Schaffer, AA; Zhang, J; Zhang, Z; Miller, W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997;25:3389–4022. doi: 10.1093/nar/25.17.3389. [PubMed]
  • Ermolaeva, MD; Khalak, HG; White, O; Smith, HO; Salzberg, SL. Prediction of transcription terminators in bacterial genomes. Journal of Molecular Biology. 2000;301:27–33. doi: 10.1006/jmbi.2000.3836. [PubMed]
  • Lavigne, R; Sun, WD; Volckaert, G. PHIRE, a deterministic approach to reveal regulatory elements in bacteriophage genomes. Bioinformatics. 2004;20:629–6135. doi: 10.1093/bioinformatics/btg456. [PubMed]
  • Darling, AC; Mau, B; Blattner, FR; Perna, NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome. 2004;14:1394–1403. doi: 10.1101/gr.2289704.
  • Zafar, N; Mazumder, R; Seto, D. CoreGenes: a computational tool for identifying and cataloging "core" genes in a set of small genomes. BMC Bioinformatics. 2002;3:12. doi: 10.1186/1471-2105-3-12. [PubMed]
  • Laemmli, UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970;227:680–685. doi: 10.1038/227680a0. [PubMed]
  • Crooks, GE; Hon, G; Chandonia, JM; Brenner, SE. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–1190. doi: 10.1101/gr.849004. [PubMed]