Skip to main content
  • Short genome report
  • Open access
  • Published:

Complete genome sequence of bacteriophage P26218 infecting Rhodoferax sp. strain IMCC26218

Abstract

Bacteriophage P26218 is a virus that thrives in freshwater and infects Rhodoferax sp. strain IMCC26218, both of which were isolated from Soyang Lake, Korea. The bacterial host, IMCC26218, belongs to the genus Rhodoferax and is closely related to R. saidenbachensis, with 98.7 % 16S rRNA gene sequence similarity. Bacteriophage P26218 has an icosahedral head structure with a diameter of ~52 nm and short tail of ~9 nm, which is a typical morphology of the Podoviridae family. Its complete dsDNA genome was 36,315 bp with 56.7 % G + C content. This is the first genome sequence reported for a lytic phage of the genus Rhodoferax.

Introduction

Bacteriophages, which are obligate parasites of bacterial cells, are the most abundant biological entities that can be found in all biospheres [13]. Considering the fact that phages heavily influence the bacterial community structure [4] and various biochemical cycles such as the carbon cycle [5], understanding the genetic potential and diversity of phages would be important in the study of microbial community dynamics. Due to the lack of a universal phylogenetic marker gene to help understand phage diversity, several studies have been reported that include a survey of entire phage populations via metagenomics, from various environments including seawater, hot springs, soil, and freshwater [68]. These viral metagenomic studies demonstrate the extremely diverse nature and novel genetic repertoire of viruses, but the limited number of phage genomes poses a challenge for interpretation of virome data. Difficulty in phage isolation and genome sequencing is simply due to the lack of available bacterial hosts, since many of bacteria in natural environments are yet to be cultured [2]. Therefore, isolation of phages infecting major groups of bacteria and unveiling their genomic information are required to provide detailed information about each phage and enable meaningful interpretation of virome data.

The class Betaproteobacteria is often the most abundant group in freshwater environments, though less abundant in marine environments [9, 10]. Metagenomic studies on several freshwater bacteria revealed that the family Comamonadaceae , arbitrarily named betI [9], is the most frequently found family [11] within this class. The genus Rhodoferax [12], belonging to the family Comamonadaceae , is found in diverse habitats including ditch water, activated sludge, Antarctic microbial mats, and water reservoirs [10, 1214]. Additionally, this is one of the most abundant genera within the 16S rRNA gene database [15]. Therefore, understanding the ecology of the genus Rhodoferax and its lytic phage will contribute to the understanding of freshwater microbial dynamics and help in further freshwater phage genomic studies. To isolate bacteriophages infecting Rhodoferax spp., we successfully isolated phage P26218, which infects Rhodoferax sp. IMCC26218 and further details of its genome features and annotations are described below.

Virus information

Classification and features

A bacteriophage, designated P26218 that infects the bacterial strain IMCC26218 was isolated from Soyang Lake, located inland of Gangwon-do, Korea, in October 2014. A bacterial strain, IMCC26218, was also isolated from the same site, using standard dilution plating technique on R2A agar (Becton, Dickenson and Company, Franklin Lakes, NJ, USA) in April 2014. Based on a comparative 16S rRNA gene sequence analyses, strain IMCC26218 was found to belong to the genus Rhodoferax with 98.7 % sequence similarity to R. saidenbachensis ED16T. To screen a representative lytic phage infecting representatives of the class Betaproteobacteria , Rhodoferax sp. IMCC26218 was used as the bacterial host.

Phage P26218 is a lytic phage that forms plaques of 1 to 2 mm in diameter, on Rhodoferax sp. IMCC26218 culture plates. Transmission electron microscopy of purified phage particles revealed its icosahedral-shaped head (52.1 nm in diameter) with a 9.4-nm long short tail (Fig. 1). The capsid encapsulates a linear dsDNA with length of 36,315 bp with 56.7 % G + C content. The morphology of the viral particle, including a uniform, icosahedral-shaped head with a short tail indicated that this phage belonged to the family Podoviridae of the order Caudovirales [16]. However, when its genomic characteristics were considered, no similar genomic architecture was found among the known viral genera, leaving P26218 without an assigned genus. The amino acid sequence of DNA polymerase I (encoded by polA) of P26218, one of the widely used viral phylogenetic markers [17, 18], was aligned with that of representative strains of the families Podoviridae and Siphoviridae and the aligned sequences were used for phylogenetic analysis. The phylogenetic tree based on DNA polymerase I revealed that P26218 formed a clade with a marine metagenome sequence, parted from previously known type species, confirming limitations in its assignment to a known genus (Fig. 2). A summary of the general phylogenetic features and isolation information are shown in Table 1.

Fig. 1
figure 1

Transmission electron micrographs of phage P26218 particles infecting Rhodoferax sp. IMCC26218. The TEM images were obtained using Philips CM200 electron microscope. Scale bars represent 100 nm in (a) and 20 nm in (b)

Fig. 2
figure 2

Phylogenetic tree highlighting the relationship of phage P26218 infecting Rhodoferax sp. IMCC26218 with representatives of the families Podoviridae and Siphoviridae. Sequences of DNA polymerase I (polA) collected from NCBI were aligned using CLUSTALW software [37], with Bacillus phage SPO1 (NC011421.1) as an outgroup. The phylogenetic tree was generated using the neighbor-joining method implemented in MEGA6 [38]. Bootstrap values representing over 60 % in 1,000 replicates are shown in the tree

Table 1 Classification and general features of phage P26218 infecting Rhodoferax sp. IMCC26218

Genome sequencing information

Genome project history

Compared to phage genomics and viromics in marine environments, relatively fewer studies with phage isolation and viral metagenome have been conducted in freshwater environments. Bacteriophage P26218 is the first lytic phage identified that infects the genus Rhodoferax , one of the representatives of the class Betaproteobacteria in freshwater environments. In this study, both virus and host were isolated from Soyang Lake in Korea. This phage was selected for genome sequencing as an initial approach to understand phages infecting members of the Betaproteobacteria isolated from surface freshwater in Korea. Genomic DNA was sequenced by the ChunLab Inc. Genome assembly, annotation, and submission to GenBank were performed at the Department of Biological Sciences, Inha University. This genome project was registered in Genomes Online Database, with accession ID, Gp0111341 as well as GenBank, with an accession ID of KP792623. A summary of the project information is described in Table 2.

Table 2 Project information

Growth conditions and genomic DNA preparation

The bacterial host, Rhodoferax sp. IMCC26218, was routinely cultured and maintained on R2A agar or in R2A broth (MB Cell, Los Angeles, CA, USA) at 20 °C. To screen lytic phages infecting this bacterial host, 10 l of water sample was collected from Soyang Lake at a depth of 1 m. The water sample was initially filtered using a 0.2-μm polyvinylidene difluoride membrane filter (Merck Millipore, Darmstadt, Germany) to remove bacterial-sized particles. To 400 ml filtrate, 100 ml of 5× R2A broth and 20 ml of IMCC26218 culture in the exponential phase were added, followed by incubation at 20 °C for 2 weeks for enrichment of bacteriophages that specifically infect Rhodoferax sp. IMCC26218. During the incubation period, 10 ml of the enrichment culture was sub-sampled 5 times at a 3-day interval. Each sub-sample was treated with approximately 3 ml chloroform to inactivate the bacterial cells. The treated samples were used for spot-double agar layer plaque assay on a Rhodoferax sp. IMCC26218 lawn plate for phage screening via appearance of plaques [19], resulting in the isolation of phageP26218.

The purification of phage P26218 genomic DNA was performed as per the ‘Molecular Cloning: A Laboratory Manual’ [20] with minor modifications. To 200 ml of phage lysate prepared for DNA purification, 1 μg ml−1 of DNase I and RNase A were added, followed by 11.7 g of NaCl. The obtained mixture was transferred to centrifuge bottles, to which PEG 8000 (Sigma-Aldrich, St. Louis, MO, USA) was added to attain a concentration of 10 % (w/v). After overnight incubation at 4 °C, the mixture was centrifuged at 11,000 × g for 40 min, supernatant was discarded by gentle inversion of the bottle, and the pellet was resuspended in 3–5 ml of SM buffer (50 mM Tris–HCl, pH 7.5; 100 mM NaCl; 10 mM MgSO4 · 7H2O; 0.01 % gelatin). PEG was removed from the liquid by treating it with equal volume of chloroform. The aqueous phase was then collected and further concentrated by ultracentrifugation at 246,000 × g for 2 h using a Beckman Coulter L-90 K ultracentrifuge with a SW 50 Ti swinging-bucket rotor. The phage pellet was resuspended in 100 μl SM buffer and used for genomic DNA extraction, using Qiagen DNeasy Blood and Tissue Kit, according to the manufacturer instructions.

Genome sequencing and assembly

The genome of phage P26218 was sequenced at ChunLab Inc. using Illumina MiSeq system with 2 × 300-bp paired-end reads. The Illumina platform produced a total of 2 × 798,245 reads. The initial total reads were split by 2 × 50,000 reads into 16 sets [21] to facilitate the assembly process. Each set of sequence reads was independently assembled using SPAdes-3.1.1 [22], yielding a single contig but with different start points. Gap-closing PCR was performed with primers designed within the end region of a contig, which resulted in the circularization of the genome sequence. Circular assembly of the genome sequence suggested that the phage genome is terminally redundant or circularly permuted. This procedure for genome sequencing and assembly finally produced 36,315 bp with approximately 2,500× fold-coverage of the genome.

Genome annotation

The ORFs were predicted using 3 gene prediction programs: GeneMark.hmm version 3.25 [23], Rapid Annotation using Subsystem Technology server version 2.0 [24], and NCBI Gene Locator and Interpolated Markov ModelER version 3.02 [25]. Only the ORFs that were identified by 2 of the 3 gene-prediction programs were included in the annotation. Each predicted ORF was translated and used to search for its homologous proteins and predict its domains using the NCBI BLASTP [17, 26], HHpred server [27] and HMMER [28] upon NCBI non-redundant database [26], the Conserved Domain Database [29], Pfam database [30], COG [31], PRK [29], and TIGRFam [32]. Then, TMHMM [33] and SignalP [34] were used to predict transmembrane helices and signal peptides.

Genome properties

The properties and statistics of P26218 genome are summarized in Tables 3 and 4. The total length of the P26218 dsDNA genome was found to be 36,315 bp with 56.7 % G + C content. All 44 predicted ORFs were protein-coding sequences. However, only 15 of them were assigned to putative protein functions, while 29 were assigned to hypothetical proteins. One gene with a signal peptide was identified but none were found to have transmembrane helices.

Table 3 Nucleotide content and gene-count levels of the genome
Table 4 Number of genes associated with general COG functional categories

Insights from the genome sequence

According to the genome annotation, bacteriophage P26218 is a unique phage, with no closely related phages. Therefore, this phage could only be classified based on its morphological characteristics, which attributed it to the family Podoviridae. Out of 44 predicted ORFs, only 15 (34 %) were assigned with a known function. As shown in Fig. 3, four ORFs were predicted to be related to DNA replication, 2 to DNA metabolism, 5 to packaging and structural functions, and 4 to other known functions (Additional file 1). BLASTP analyses showed that each ORF with an identified function was homologous to ORFs from different phages belonging to different viral families. All ORFs encoding viral packaging function were closely related to those of other viruses in the family Podoviridae. The ORFs encoding DNA polymerase I, ATPase component, thymidylate synthase, and hydrolase-like protein were similar to those of the family Siphoviridae, while the genes for DnaB-like ATP-dependent helicase and ParB-like nuclease domain showed a high degree of homology to those of the family Myoviridae This genomic architecture of P26218 confirmed the mosaic genome structure, known to be a result of lateral gene transfer usually predicted in viral genomes in attempts to enhance their genetic diversity [35, 36] and often observed in species of the order Caudovirales such as phages P22 and lambda.

Fig. 3
figure 3

Genome map of Rhodoferax phage P26218. Total length of the genome is 36,315 bp and contig functions are color-coded as follows: light grey represents hypothetical proteins, yellow represents DNA metabolism, red represents DNA replication, blue represents structural and packaging genes and green represents other known functional genes

Conclusions

Lytic bacteriophage P26218, isolated from a freshwater lake is the first virus identified that infects the genus Rhodoferax . Based on its morphology, this phage was identified to be a member of the family Podoviridae, with an icosahedral-shaped head and short tail. All predicted ORFs from this phage genome were protein-coding, with 3 specifically coding for DNA replication, 7 for DNA metabolism, and 5 for packaging and structural proteins. The group of ORFs with similar function was postulated to originate from different groups of viral families (Podoviridae, Siphoviridae, and Myoviridae), which was indicative of the mosaic property of the P26218 genome. It is expected that phage P26218 isolated in this research and its genome sequence would be further used to study bacteria-phage interactions in freshwater environments, to reveal the evolutionary role of phage lateral gene transfer and to interpret freshwater virome data.

Abbreviations

PEG:

Polyethylene glycol

References

  1. Casas V, Rohwer F. Phage metagenomics. Methods Enzymol. 2007;421:259–68.

    Article  CAS  PubMed  Google Scholar 

  2. Edwards RA, Rohwer F. Viral metagenomics. Nat Rev Microbiol. 2005;3:504–10.

    Article  CAS  PubMed  Google Scholar 

  3. Rohwer F. Global phage diversity. Cell. 2003;113:141.

    Article  CAS  PubMed  Google Scholar 

  4. Rodriguez-Brito B, Li L, Wegley L, Furlan M, Angly F, Breitbart M, et al. Viral and microbial community dynamics in four aquatic environments. ISME J. 2010;4:739–51.

    Article  PubMed  Google Scholar 

  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.

    Article  CAS  PubMed  Google Scholar 

  6. Suttle CA. Viruses in the sea. Nature. 2005;437:356–61.

    Article  CAS  PubMed  Google Scholar 

  7. Breitbart M, Wegley L, Leeds S, Schoenfeld T, Rohwer F. Phage community dynamics in hot springs. Appl Environ Microbiol. 2004;70:1633–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Srinivasiah S, Bhavsar J, Thapar K, Liles M, Schoenfeld T, Wommack KE. Phages across the biosphere: contrasts of viruses in soil and aquatic environments. Res Microbiol. 2008;159:349–57.

    Article  CAS  PubMed  Google Scholar 

  9. Newton RJ, Jones SE, Eiler A, McMahon KD, Bertilsson S. A guide to the natural history of freshwater lake bacteria. Microbiol Mol Biol R. 2011;75:14–49.

    Article  CAS  Google Scholar 

  10. Zwart G, Crump BC, Kamst-van Agterveld MP, Hagen F, Han S-K. Typical freshwater bacteria: an analysis of available 16S rRNA gene sequences from plankton of lakes and rivers. Aquat Microbiol Ecol. 2002;28:141–55.

    Article  Google Scholar 

  11. Cottrell MT, Waidner LA, Yu L, Kirchman DL. Bacterial diversity of metagenomic and PCR libraries from the Delaware River. Environ Microbiol. 2005;7:1883–95.

    Article  CAS  PubMed  Google Scholar 

  12. Kaden R, Spröer C, Beyer D, Krolla-Sidenstein P. Rhodoferax saidenbachensis sp. nov., a psychrotolerant, very slowly growing bacterium within the family Comamonadaceae, proposal of appropriate taxonomic position of Albidiferax ferrireducens strain T118T in the genus Rhodoferax and emended description of the genus Rhodoferax. Int J Syst Evol Microbiol. 2014;64:1186–93.

    Article  CAS  PubMed  Google Scholar 

  13. Madigan MT, Jung DO, Woese CR, Achenbach LA. Rhodoferax antarcticus sp. nov., a moderately psychrophilic purple nonsulfur bacterium isolated from an Antarctic microbial mat. Arch Microbiol. 2000;173:269–77.

    Article  CAS  PubMed  Google Scholar 

  14. Hiraishi A, Hoshino Y, Satoh T. Rhodoferax fermentans gen. nov., sp. nov., a phototrophic purple nonsulfur bacterium previously referred to as the “Rhodocyclus gelatinosus-like” group. Arch Microbiol. 1991;155:330–6.

    Article  Google Scholar 

  15. Salka I, Čuperová Z, Mašín M, Koblížek M, Grossart HP. Rhodoferax‐related pufM gene cluster dominates the aerobic anoxygenic phototrophic communities in German freshwater lakes. Environ Microbiol. 2011;13:2865–75.

    Article  CAS  PubMed  Google Scholar 

  16. King AM, Adams MJ, Carstens EB, Lefkowitz EJ. The double stranded DNA viruses. In: Virus Taxonomy: Classification and Nomenclature of Viruses; Ninth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier Academic Press; 2012. p. 63–85.

    Google Scholar 

  17. Adriaenssens EM, Cowan DA. Using signature genes as tools to assess environmental viral ecology and diversity. Appl Environ Microbiol. 2014;80:4470–80.

    Article  PubMed Central  PubMed  Google Scholar 

  18. Breitbart M, Miyake JH, Rohwer F. Global distribution of nearly identical phage-encoded DNA sequences. FEMS Microbiol Lett. 2004;236:249–56.

    Article  CAS  PubMed  Google Scholar 

  19. Grabow W. Bacteriophages: update on application as models for viruses in water. Water Sa. 2004;27:251–68.

    Article  Google Scholar 

  20. Green MR, Sambrook J. Molecular cloning: a laboratory manual. New York: Cold Spring Harbor Laboratory Press; 2012.

    Google Scholar 

  21. Lonardi S, Mirebrahim H, Wanamaker S, Alpert M, Ciardo G, Duma D, et al. When less is more: “slicing” sequencing data improves read decoding accuracy and de novo assembly quality. bioRxiv. 2015; doi:10.1101/013425

  22. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Lukashin AV, Borodovsky M. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Delcher AL, Harmon D, Kasif S, White O, Salzberg SL. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999;27:4636–41.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. 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. Res Microbiol. 2008;159:406–14.

    Article  CAS  PubMed  Google Scholar 

  27. Söding J. Protein homology detection by HMM–HMM comparison. Bioinformatics. 2005;21:951–60.

    Article  PubMed  Google Scholar 

  28. Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14:755–63.

    Article  CAS  PubMed  Google Scholar 

  29. Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011;39:D225–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:D222–30.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31:371–3.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. 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.

    Article  CAS  PubMed  Google Scholar 

  34. Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.

    Article  CAS  PubMed  Google Scholar 

  35. Yoshida M, Yoshida-Takashima Y, Nunoura T, Takai K. Genomic characterization of a temperate phage of the psychrotolerant deep-sea bacterium Aurantimonas sp. Extremophiles. 2015;19:49–58.

    Article  CAS  PubMed  Google Scholar 

  36. Swanson MM, Reavy B, Makarova KS, Cock PJ, Hopkins DW, Torrance L, et al. Novel bacteriophages containing a genome of another bacteriophage within their genomes. PLoS ONE. 2012;7:e40683.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Thompson JD, Gibson T, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002;Chapter 2:Unit 2.3.

    PubMed  Google Scholar 

  38. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by Mid-Career Research Program through National Research Foundation (NRF) funded by the Ministry of Science, ICT and Future Planning (to J-CC, NRF-2013R1A2A2A01068004) and partially by the General Research Program through NRF (to S-JK; 2012R1A1B3003609), Korea.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jang-Cheon Cho or Sang-Jong Kim.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

IK and J-CC designed the study. SK isolated the bacterial host and collected freshwater samples. KM performed laboratory experiments, analyzed the data, and drafted the manuscript. KM, IK, J-CC, and S-JK wrote the manuscript together and finalized the study. We all authors have read and approved the final manuscript.

Additional file

Additional file 1: Table S1.

Description of data: Gene annotation table of bacteriophage P26218. (PDF 183 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Moon, K., Kang, I., Kim, S. et al. Complete genome sequence of bacteriophage P26218 infecting Rhodoferax sp. strain IMCC26218. Stand in Genomic Sci 10, 111 (2015). https://doi.org/10.1186/s40793-015-0090-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40793-015-0090-1

Keywords