- Short genome report
- Open Access
Complete genome sequence of bacteriophage P8625, the first lytic phage that infects Verrucomicrobia
Standards in Genomic Sciencesvolume 10, Article number: 96 (2015)
Bacteriophage P8625 is a lytic bacteriophage that infects the verrucomicrobial strain IMCC8625, a marine bacterium affiliated with Verrucomicrobia subdivision 4. Both the bacteriophage and the host bacterial strain were isolated from surface seawater samples collected off the east coast of Korea. The phage particle has an icosahedral capsid with a diameter of ~47 nm and a long tail of ~75 nm in length, showing the distinctive morphology of the Siphoviridae family. The complete genome sequence of phage P8625 is 32,894 bp long with 51.0 % G + C content. This is the first report of the complete genome sequence of a lytic phage that infects the Verrucomicrobia, for which the name “verrucophage” is proposed.
Marine viruses are the most abundant biological components in the ocean, with ~107 virus-like particles per milliliter in surface seawater [1–3]. Marine bacteriophages lyse specific bacterial hosts, controlling bacterial abundance and diversity and influencing biogeochemical cycles, which makes the study of marine viruses ecologically important [2, 4, 5]. Recent metagenomic studies have demonstrated immense genetic diversity among marine viruses [6–8]. However, these studies had difficulty in the phylogenetic interpretation of metagenomes due to the shortage of genetic data from representative marine viral isolates [7, 9]. Therefore, isolation of bacteriophages is also important for a better understanding of marine virome data [10, 11].
Bacterial members of the phylum Verrucomicrobia are distributed widely in the ocean environment, albeit in low densities . Based on 16S rRNA gene sequence analyses, the phylum Verrucomicrobia was found to comprise an average of 2 % of the water column and 1.4 % of the sediment bacterial population . Compared to marine environments, in terrestrial ecosystems such as soil environments, the phylum Verrucomicrobia is highly abundant and often exceed 20 % of the total rRNA gene sequences . To our knowledge, however, no lytic phage infecting verrucomicrobial strains has been reported , perhaps largely due to the scarcity of bacterial hosts. Interestingly, a few recent studies based on fosmid sequencing or single cell genomics have hinted at the presence of verrucomicrobial phages in marine environments [15–17].
In this study, we report the isolation and genomic characterization of bacteriophage P8625, a novel marine siphovirus that infects the marine bacterial strain IMCC8625. Because the bacterial strain belongs to the class Opitutae  of the phylum Verrucomicrobia , phage P8625, isolated off the east coast of Korea, is regarded as the first isolated lytic phage of the Verrucomicrobia , for which the name “verrucophage” is proposed.
Classification and features
Bacterial strain IMCC8625 was isolated using high-throughput cultivation based on dilution to extinction , during a survey of microbial assemblages inhabiting coastal seawater of the East Sea. A comparison of 16S rRNA gene sequences indicated that strain IMCC8625 was closely related to Coraliomargarita akajimensis DSM 45221T (95.4 % similarity), a marine verrucomicrobial strain [20, 21]. Strain IMCC8625 was used as the bacterial host for the screening of marine lytic bacteriophages, which resulted in the isolation of verrucophage P8625.
Verrucophage P8625 was isolated from a surface seawater sample collected off the east coast of South Korea, in the East Sea (Sea of Japan). It is a lytic phage, forming plaques with a diameter of 1 to 2 mm after 5 days of infection of strain IMCC8625. Transmission electron microscopy of purified phage particles showed that P8625 had an icosahedral capsid approximately 47 nm in diameter and a long tail 71 ~ 75 nm in length (Fig. 1a). Bacteriophage P8625 attaches itself to and infects the host cell (Fig. 1b). The tail appears to be long and non-contractile, which is typically seen in members of the Siphoviridae family . The genome of P8625 is comprised of dsDNA with a length of 32,894 bp and 51.0 % G + C content. The morphological evidence (Fig. 1), together with the dsDNA genome, has led to the tentative classification of P8625 as a member of the Siphoviridae family, with an unassigned genus. A summary of the isolation and general phylogenetic features of phage P8625 are shown in Table 1.
Genome sequencing information
Genome project history
Bacteriophage P8625 was selected for genome sequencing because of its importance to understand bacteriophages infecting members of the Verrucomicrobia that are known to play crucial ecological roles in the ocean . This bacteriophage that was shown to infect bacterium IMCC8625 in this study is the first lytic phage infecting Verrucomicrobia . This genome sequence has been submitted to GenBank, and the project information is also available in GOLD. A summary of the project information is shown in Table 2.
Growth conditions and genomic DNA preparation
Verrucophage P8625 was isolated from a surface seawater sample collected at the station where the host strain IMCC8625 had been isolated, using the standard plaque assay after being enriched with the host. For the enrichment, the seawater sample was filtered using a 0.22-μm polyethersulfone membrane filter (Durapore, Millipore) to remove bacterial particles. To 400 ml filtered seawater, 100 ml of 5× R2A broth (BD Difco) and 30 ml of exponentially grown IMCC8625 culture were added. During incubation of this enrichment culture at 20 °C for 2 weeks, 10 ml of the culture was removed and treated with 2 ml of chloroform at 3 days interval. After centrifugation of chloroform-treated cultures to remove bacterial debris, bacteria-free supernatants were used for double agar overlay plaque assay. From an isolated plaque on the final assay plate, a single strain of phage was established and designated P8625.
The purification of phage DNA followed the method outlined in Molecular Cloning: A Laboratory Manual  with minor modifications. Approximately 2 l of phage lysates were prepared for DNA purification. DNase I and RNase A were added at a final concentration of 1 μg per ml. Then, 116.9 g of NaCl was dissolved in the lysates and cooled at 4 °C. After about 1 h, the mixed lysates were centrifuged at 10,000 × g for 30 min at 4 °C to remove the debris. Phage particles in the supernatant were precipitated with 10 % (w/v) PEG 8000 (Sigma-Aldrich). After an overnight incubation at 4 °C, the mixture was pelleted at 10,000 × g for 25 min at 4 °C and the pellet was gently resuspended in 2 ml SM buffer (50 mM Tris–HCl, pH 7.5; 100 mM NaCl; 10 mM MgSO4 · 7H2O; 0.01 % gelatin). PEG was extracted by treatment with an equal volume of chloroform. The phages were then ultracentrifuged at 246,000 × g for 2 h at 4 °C in an L-90 K ultracentrifuge (Beckman) with an SW 55 Ti rotor. Pelleted phage particles were resuspended with 100 μl of SM buffer overnight at 4 °C, and the purified phages were stored in the dark at 4 °C. The genomic DNA of P8625 was extracted using a silica based spin column (Qiagen DNeasy Blood and Tissue Kit) according to the manufacturer’s instructions.
Genome sequencing and assembly
The genomic DNA of P8625 was sequenced at ChunLab, Inc. using an Illumina Miseq system with 2 × 300 bp paired-end reads. Assembly of the resulting reads was performed using SPAdes version 3.1.1 . The Illumina platform provided 3085 × fold coverage of the genome. The genome was assembled into one contig through PCR-based gap closing.
The prediction of genes in the genome was performed using a combination of three gene calling methods: the RAST server , Genemark.hmm 3.25 , and GLIMMER version 3.02 . Assignment of protein function to ORFs was performed manually using BLASTp against the NCBI nonredundant database and RPS-BLAST or HMMER search against the COG database , Pfam database , and TIGRFam database . InterProScan was also used . Search results were combined to assign a putative function for each predicted protein. TMHMM  and SignalP  were used to predict transmembrane helices and signal peptides, respectively.
The properties and statistics of the genome are summarized in Table 3. The genome of P8625 was assembled as a 32,894 bp dsDNA, with a G + C content of 51.0 %. In total, 52 putative CDSs were predicted in the phage genome (Fig. 2). Of the 52 predicted protein-coding genes, 13 genes were assigned putative functions, whereas the remaining genes were annotated as coding for hypothetical proteins. Four proteins with transmembrane helices were identified, but signal peptides were not detected in any protein. The distribution of genes into COG functional categories is presented in Table 4.
Insights from the genome sequence
When all 52 CDSs predicted from the P8625 genome were subjected to functional annotation based mainly on a conserved domain/motif search, only 13 CDSs were found to have specific functions. These functions were related to head and tail assembly, DNA packaging, or DNA replication and metabolism (Fig. 2 and Additional file 1: Table S1). Head and tail structural proteins were similar to those of the Siphoviridae family, while some proteins involved in DNA packaging and DNA replication/metabolism were similar to those of the Myoviridae family, suggesting mosaicism of the genome. The genome contained a putative RecT protein (CDS34) and a protein containing a PAPS reductase domain (CDS45) (Additional file 1: Table S1 and Fig. 2). RecT protein is involved in recombination repair of DNA , and homologs of RecT have been found in many phages . The function of the protein containing the PAPS reductase domain (CDS45) is unclear. Although we cannot exclude the possibility of this protein being involved in sulfate metabolism, the location of this protein in the genome (Fig. 2) and the affiliation of PAPS reductase domain with the adenine nucleotide alpha hydrolase superfamily (cl00292) suggest that this protein could play a role in DNA metabolism.
A TerL, a commonly-used viral phylogenetic marker [37, 38], was predicted in the phage P8625 genome and used for phylogenetic analyses together with other representative tailed phages and Verrucomicrobia prophages (Fig. 3). In the resulting phylogenetic tree, the terminases from P8625 and Verrucomicrobia prophages formed a robust cluster supported by high bootstrap values (Fig. 3), which suggested that P8625 has a strong relationship with Verrucomicrobia prophages, without a clear affiliation with any of the three families of dsDNA tailed phages.
Several proteins predicted from the P8625 genome, including DNA methylase (CDS3), TerL (CDS6), and head structure proteins (CDS9–11), were similar to proteins annotated from a short (~14 kb) contig of a verrucomicrobial single cell genome SCGC AAA027-I19. This contig might have originated from a prophage in the host genome or might have been derived from a lytic phage during an infection cycle when the host cell was selected for multiple displacement amplification.
The virulent verrucophage P8625, a siphovirus isolated from seawater in the East Sea, represents the first identification of a Verrucomicrobia phage. The dsDNA genome of P8625 was 32,894 bp long, with a 51.0 % G + C content. Among the 52 protein-encoding genes predicted in the genome, only 13 genes were functionally annotated, while the remaining 75 % of the genes were annotated as encoding hypothetical proteins. Verrucophage P8625 and its genome sequence can improve our understanding of the interaction of marine verrucomicrobial strains and their phages.
Terminase large subunit
Suttle CA. Viruses in the sea. Nature. 2005;437:356–61.
Fuhrman JA. Marine viruses and their biogeochemical and ecological effects. Nature. 1999;399:541–8.
Breitbart M. Marine viruses: truth or dare. Ann Rev Mar Sci. 2012;4:425–48.
Suttle CA. Marine viruses-major players in the global ecosystem. Nat Rev Microbiol. 2007;5:801–12.
Weinbauer MG. Ecology of prokaryotic viruses. FEMS Microbiol Rev. 2004;28:127–81.
Ma YF, Allen LZ, Palenik B. Diversity and genome dynamics of marine cyanophages using metagenomic analyses. Env Microbiol Rep. 2014;6:583–94.
Angly FE, Felts B, Breitbart M, Salamon P, Edwards RA, Carlson C, et al. The marine viromes of four oceanic regions. PLoS Biol. 2006;4:e368.
Hurwitz BL, Sullivan MB. The Pacific Ocean Virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS ONE. 2013;8:e57355.
Williamson SJ, Allen LZ, Lorenzi HA, Fadrosh DW, Brami D, Thiagarajan M, et al. Metagenomic exploration of viruses throughout the Indian Ocean. PLoS ONE. 2012;7:e42047.
Kang I, Oh HM, Kang D, Cho JC. Genome of a SAR116 bacteriophage shows the prevalence of this phage type in the oceans. Proc Natl Acad Sci U S A. 2013;110:12343–8.
Zhao Y, Temperton B, Thrash JC, Schwalbach MS, Vergin KL, Landry ZC, et al. Abundant SAR11 viruses in the ocean. Nature. 2013;494:357–60.
Freitas S, Hatosy S, Fuhrman JA, Huse SM, Welch DB, Sogin ML, et al. Global distribution and diversity of marine Verrucomicrobia. ISME J. 2012;6:1499–505.
Bergmann GT, Bates ST, Eilers KG, Lauber CL, Caporaso JG, Walters WA, et al. The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol Biochem. 2011;43:1450–5.
Wommack KE, Nasko DJ, Chopyk J, Sakowski EG. Counts and sequences, observations that continue to change our understanding of viruses in nature. J Microbiol. 2015;53:181–92.
Mizuno CM, Rodriguez-Valera F, Kimes NE, Ghai R. Expanding the marine virosphere using metagenomics. PLoS Genet. 2013;9, e1003987.
Rodriguez-Valera F, Mizuno CM, Ghai R. Tales from a thousand and one phages. Bacteriophage. 2014;4:e28265.
Swan BK, Tupper B, Sczyrba A, Lauro FM, Martinez-Garcia M, Gonzalez JM, et al. Prevalent genome streamlining and latitudinal divergence of planktonic bacteria in the surface ocean. Proc Natl Acad Sci U S A. 2013;110:11463–8.
Choo YJ, Lee K, Song J, Cho JC. Puniceicoccus vermicola gen. nov., sp. nov., a novel marine bacterium, and description of Puniceicoccaceae fam. nov., Puniceicoccales ord. nov., Opitutaceae fam. nov., Opitutales ord. nov. and Opitutae classis nov. in the phylum ‘Verrucomicrobia’. Int J Syst Evol Microbiol. 2007;57:532–7.
Connon SA, Giovannoni SJ. High-throughput methods for culturing microorganisms in very-low-nutrient media yield diverse new marine isolates. Appl Environ Microbiol. 2002;68:3878–85.
Yoon J, Yasumoto-Hirose M, Katsuta A, Sekiguchi H, Matsuda S, Kasai H, et al. Coraliomargarita akajimensis gen. nov., sp. nov., a novel member of the phylum ‘Verrucomicrobia’ isolated from seawater in Japan. Int J Syst Evol Microbiol. 2007;57:959–63.
Mavromatis K, Abt B, Brambilla E, Lapidus A, Copeland A, Deshpande S, et al. Complete genome sequence of Coraliomargarita akajimensis type strain (04OKA010-24). Stand Genomic Sci. 2010;2:290–9.
Plisson C, White HE, Auzat I, Zafarani A, Sao-Jose C, Lhuillier S, et al. Structure of bacteriophage SPP1 tail reveals trigger for DNA ejection. EMBO J. 2007;26:3720–8.
Martinez-Garcia M, Brazel DM, Swan BK, Arnosti C, Chain PS, Reitenga KG, et al. Capturing single cell genomes of active polysaccharide degraders: an unexpected contribution of Verrucomicrobia. PLoS ONE. 2012;7, e35314.
Green MR, Sambrook J. Molecular cloning: a laboratory manual. New York: Cold Spring Harbor Laboratory Press; 2012.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.
Besemer J, Borodovsky M. Heuristic approach to deriving models for gene finding. Nucleic Acids Res. 1999;27:3911–20.
Delcher AL, Bratke KA, Powers EC, Salzberg SL. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics. 2007;23:673–9.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6.
Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40:D290–301.
Haft DH, Selengut JD, Richter RA, Harkins D, Basu MK, Beck E. TIGRFAMs and Genome Properties in 2013. Nucleic Acids Res. 2013;41:D387–95.
Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Hall SD, Kolodner RD. Homologous pairing and strand exchange promoted by the Escherichia coli RecT protein. Proc Natl Acad Sci U S A. 1994;91:3205–9.
Iyer LM, Koonin EV, Aravind L. Classification and evolutionary history of the single-strand annealing proteins, RecT, Redbeta, ERF and RAD52. BMC Genomics. 2002;3:8.
Casjens SR, Gilcrease EB, Winn-Stapley DA, Schicklmaier P, Schmieger H, Pedulla ML, et al. The generalized transducing Salmonella bacteriophage ES18: complete genome sequence and DNA packaging strategy. J Bacteriol. 2005;187:1091–104.
Sullivan MB, Krastins B, Hughes JL, Kelly L, Chase M, Sarracino D, et al. The genome and structural proteome of an ocean siphovirus: a new window into the cyanobacterial ‘mobilome’. Environ Microbiol. 2009;11:2935–51.
King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. Virus taxonomy: classification and nomenclature of viruses: ninth report of the international committee on taxonomy of viruses. San Diego: Elsevier; 2012. p. 855–80.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
This study was supported by a grant from the Marine Biotechnology Program (to J-CC, PJT200620, Genome Analysis of Marine Organisms and Development of Functional Applications) funded by the Ministry of Oceans and Fisheries, Korea, and also by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (to IK, NRF-2013R1A1A2063660).
The authors declare that they have no competing interests.
S-JY isolated the host, Verrucomicrobia strain IMCC8625, and collected seawater samples. AC performed laboratory experiments, drafted the manuscript, and analyzed the data. IK and J-CC together organized this study and wrote the manuscript. All authors read and approved the final manuscript
Verrucophage P8625 gene annotation. (PDF 77 kb)