Skip to content


Standards in Genomic Sciences

Open Access

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

Standards in Genomic Sciences201510:111

Received: 8 April 2015

Accepted: 29 October 2015

Published: 24 November 2015


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.


Bacteriophage Rhodoferax Freshwater Podoviridae Genome


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

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

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




Evidence codea



Genome group: dsDNA viruses, no RNA stage



Phylum: unassigned


Class: unassigned


Order: Caudovirales

TAS [16]


Family: Podoviridae

TAS [16]


Genus: unassigned


Species: unassigned


Strain: P26218


Particle shape

Icosahedral head with a short tail




Freshwater lake, surface



Biotic relationship

Intracellular parasite of Rhodoferax sp. IMCC26218




Virulent phage of Rhodoferax sp. IMCC26218



Geographic location

Soyang Lake, Gangwon-do, South Korea



Sample collection

Oct. 17, 2014




37°57’11” N




127°49’02” E




1 m





aIDA: Inferred from Direct Assay and TAS: Traceable Author Statement. The evidence codes are from the Gene Ontology project [39]

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





Finishing quality



Number of contigs



Libraries used

One paired-end Illumina library


Sequencing platforms

Illumina Miseq

MIGS 31.2

Fold coverage




SPAdes version 3.1.1


Gene calling method

RAST version 2.0, GeneMark.hmm version 3.25 and GLIMMER version 3.02


GenBank ID



GenBank Date of Release

April, 2015








Source Material Identifier



Project relevance

Diversity of freshwater bacteriophage

aNot available

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



% of Totala

Genome size (bp)



DNA coding (bp)



DNA G + C (bp)



DNA scaffolds



Total genes



Protein coding genes



RNA genes



Pseudo genes



Genes in internal clusters



Genes with function prediction



Genes assigned to COGs



Genes with Pfam domains



Genes with signal peptides



Genes with transmembrane helices



CRISPR repeats



aThe total percentage is based on the total number of protein coding genes in the genome

Table 4

Number of genes associated with general COG functional categories



% of totala





Translation, ribosomal structure and biogenesis




RNA processing and modification








Replication, recombination and repair




Chromatin structure and dynamics




Cell cycle control, Cell division, chromosome partitioning




Defense mechanisms




Signal transduction mechanisms




Cell wall/membrane biogenesis




Cell motility




Intracellular trafficking and secretion




Posttranslational modification, protein turnover, chaperones




Energy production and conversion




Carbohydrate transport and metabolism




Amino acid transport and metabolism




Nucleotide transport and metabolism




Coenzyme transport and metabolism




Lipid transport and metabolism




Inorganic ion transport and metabolism




Secondary metabolite biosynthesis, transport, and catabolism




General function prediction only




Function unknown




Mobilome: prophages, transposons




Not in COGs

aThe percentage is based on the total number of protein-coding genes in the genome

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

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


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.



Polyethylene glycol



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.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

School of Biological Sciences, Seoul National University, Seoul, Republic of Korea
Department of Biological Sciences, Inha University, Incheon, Republic of Korea


  1. Casas V, Rohwer F. Phage metagenomics. Methods Enzymol. 2007;421:259–68.View ArticlePubMedGoogle Scholar
  2. Edwards RA, Rohwer F. Viral metagenomics. Nat Rev Microbiol. 2005;3:504–10.View ArticlePubMedGoogle Scholar
  3. Rohwer F. Global phage diversity. Cell. 2003;113:141.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
  6. Suttle CA. Viruses in the sea. Nature. 2005;437:356–61.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticleGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
  18. Breitbart M, Miyake JH, Rohwer F. Global distribution of nearly identical phage-encoded DNA sequences. FEMS Microbiol Lett. 2004;236:249–56.View ArticlePubMedGoogle Scholar
  19. Grabow W. Bacteriophages: update on application as models for viruses in water. Water Sa. 2004;27:251–68.View ArticleGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
  23. Lukashin AV, Borodovsky M. GeneMark.hmm: new solutions for gene finding. Nucleic Acids Res. 1998;26:1107–15.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
  27. Söding J. Protein homology detection by HMM–HMM comparison. Bioinformatics. 2005;21:951–60.View ArticlePubMedGoogle Scholar
  28. Eddy SR. Profile hidden Markov models. Bioinformatics. 1998;14:755–63.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
  32. Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31:371–3.PubMed CentralView ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar
  37. Thompson JD, Gibson T, Higgins DG. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002;Chapter 2:Unit 2.3.PubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle 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.PubMed CentralView ArticlePubMedGoogle Scholar


© Moon et al. 2015