Skip to main content

Genome-wide characterization of vibrio phage ϕpp2 with unique arrangements of the mob-like genes

Abstract

Background

Vibrio parahaemolyticus is associated with gastroenteritis, wound infections, and septicemia in human and animals. Phages can control the population of the pathogen. So far, the only one reported genome among giant vibriophages is KVP40: 244,835 bp with 26% coding regions that have T4 homologs. Putative homing endonucleases (HE) were found in Vibrio phage KVP40 bearing one seg D and Vibrio cholerae phage ICP1 carrying one mob C/E and one seg G.

Results

A newly isolated Vibrio phage ϕpp2, which was specific to the hosts of V. parahaemolyticus and V. alginolyticus, featured a long nonenveloped head of ~90 × 150 nm and tail of ~110 nm. The phage can survive at 50°C for more than one hour. The genome of the phage ϕpp2 was sequenced to be 246,421 bp, which is 1587 bp larger than KVP40. 383 protein-encoding genes (PEGs) and 30 tRNAs were found in the phage ϕpp2. Between the genomes of ϕpp2 and KVP40, 254 genes including 29 PEGs for viral structure were of high similarity, whereas 17 PEGs of KVP40 and 21 PEGs of ϕpp2 were unmatched. In both genomes, the capsid and tail genes have been identified, as well as the extensive representation of the DNA replication, recombination, and repair enzymes. In addition to the three giant indels of 1098, 1143 and 3330 nt, ϕpp2 possessed unique proteins involved in potassium channel, gp2 (DNA end protector), tRNA nucleotidyltransferase, and mob-type HEs, which were not reported in KVP40. The ϕpp2 PEG274, with strong promoters and translational initiation, was identified to be a mob E type, flanked by NrdA and NrdB/C homologs. Coincidently, several pairs of HE-flanking homologs with empty center were found in the phages of Vibrio phages ϕpp2 and KVP40, as well as in Aeromonas phages (Aeh1 and Ae65), and cyanophage P-SSM2.

Conclusions

Vibrio phage ϕpp2 was characterized by morphology, growth, and genomics with three giant indels and different types of HEs. The gene analysis on the required elements for transcription and translation suggested that the ϕpp2 PEG274 was an active mob E gene. The phage was signified to be a new species of T4-related, differing from KVP40.

Background

Vibrio parahaemolyticus is a halophilic gram-negative bacterium that is widely distributed in coastal waters worldwide and is associated with gastroenteritis, wound infections, and septicemia[1]. Since the first report of Fujino et al.[2], numerous investigations of V. parahaemolyticus have been performed using stools of patients and diseased fish. The halophile has been found seasonally in sea water of the continental United States, Germany, the Far East, and Hawaii[36]. V. parahaemolyticus infections are frequently reported to occur due to the consumption of undercooked raw shellfish or direct contact with estuarine waters. In Asia, many recent infections have been caused by serotype O3:K6 of V. parahaemolyticus[7].

The phages can control the population of the pathogen. Among the giant T4-like phages that are specific to V. parahaemolyticus, the vibriophage KVP40 is the only strain for which the genome has been determined[8]. The size of the KVP40 genome is 244,835 bp with an overall G + C content of 42.6%. It contains 381 putative protein-encoding genes (PEG), 30 tRNAs, 33 late promoters, and 57 rho-independent terminators. The genome sequence and organization of KVP40 show a degree of conservation with phage T4. While 65% of the PEGs were unique to KVP40, 99 out of the total 381 putative coding regions have homologs in the T4 genome, which includes DNA replication, recombination, and repair enzymes as well as the viral capsid and tail structural genes. KVP40 lacks enzymes involved in DNA degradation, cytosine modification and group I introns, and it probably utilizes NAD salvage pathway that is unique among bacteriophages[8].

Phages can prompt gene recombination via homing endonucleases (HEs). In genome analyses, putative homing endonucleases (HEs) were found in Vibrio parahaemolyticus phages KVP40 and Vibrio cholerae ICP1[8, 9]. Homing endonucleases might act as possible mediators for the diversity among bacteriophage genomes by the acquisition of a novel DNA to create a new species of phage. Although more than 30 T4-related genomes have been published so far, no other known phage genome comes close to encoding the 15 homing endonucleases in T4 phage[1012]. Intron homing[13] and intronless homing[14, 15] endonucleases both utilize homologous recombination between phages to transfer the genetic elements from the HE-encoding genome to a HE-lacking recipient. The seg and mob subtypes, which are also called freestanding endonucleases, belong to the GIY-YIG and HNH homing endonuclease families, respectively[16], a review]. The seg C, seg F, seg G, mob A, and mob E of T4 endonucleases are polycistronically transcribed with their respective upstream genes, whereas the endonuclease-specific promoters for seg A, seg G, mob C and mob D are immediately upstream of the endonuclease genes[16]. There is as yet no convincing evidence that the HEs can move across the boundary of species or genera. Nevertheless, these transposable genes may leave a trace of their involvement after the transfers. The sequence analysis for the Enterobacteria phage JSE intron revealed that the putative intron contained a truncated derivative of a HE gene[17], very similar to the truncated sequence in the intron of the T4 nrd B gene, suggesting that there is a rarely-detectable trace of the mob/seg elements in contemporary phage genomes[18].

We sequenced the genome of ϕpp2 – a new T4-like Vibrio phage with mob genes – which may be another paradigm in the plausible analysis of evolution of HE families in the bacteriophages and their hosts[8]. In the same host, Vibrio parahaemolyticus, the phage ϕpp2 can complement KVP40 in studying the genome spectra of the giant T4-related Vibrio phages.

Methods

Bacteria strains and growth conditions

Vibrio strains were bought from the Bioresource Collection and Research Center, Taiwan; including V. alginolyticus ATCC 17749, V. carchariae ATCC 35084, V. damsela ATCC 33536, V. harveyi ATCC 14126, V. parahaemolyticus ATCC 17802, V. pelagius ATCC 25916, and V. vulnificus BCRC15431. V. parahaemolyticus ATCC 17802 carries O1 serotype and no tdh/trh genes[19]. The Vibrio strains were maintained in Brain Heart Infusion (BHI) medium, supplemented with 3% NaCl. For long-term preservation, bacteria were frozen in BHI supplemented with 1% NaCl and 25% glycerol. When working, the strains were streaked onto the modified sea water yeast extract (rich MSWYE) agar plates consisting of 23.4 g NaCl, 6.98 g MgSO4.7H2O, and 0.75 g KCl in 1000 ml distilled water[19]. The pH was adjusted to 7.6 with 1 N NaOH, followed by addition of 5.0 g of proteose peptone (Difco), 3.0 g of yeast extract (Difco), and 20.0 g of agar per liter.

Isolation and titer of bacteriophage

The water samples were collected from the aquaculture waterways around southern Taiwan. The enrichment procedure for the target phages has been described elsewhere[20]. In brief, 20% of MSWYE medium and 1% seed culture of Vibrio parahaemolyticus were added the micro-filtrated samples and incubated at 37°C for four hours to enrich the phages. In determining the phage concentrations, the bacterium Vibrio parahaemolyticus was freshly grown to 0.3–0.4 OD600, in about two hours, and 200 μl of cells were added to 10 μl phages in a series of dilutions for infection, followed by the Agar Overlay Technique. The plaques were counted in 3–5 hours; the titers per ml were calculated by 100*(dilution factor)*(plaque counts).

Electron microscopy

Preparation of phage particles for electron microscopy has been described elsewhere[20, 21]. In brief, bacteriophage particles were applied onto parafilm to produce a spherical drop. Carbon-coated nitrocellulose films were fabricated on copper grids and placed face down on the sample drops for 1 min to absorb the particles. The samples were stained with freshly prepared 2% uranyl acetate (UA; Tris–HCl, pH 8.0) for 60 seconds. Images of phage particles were taken at a magnification of 40,000x, defocus of 3 μm, using a 200-kV electron microscope (JEOL JEM-2010, equipped with a Gatan-832 CCD camera).

Analyses of bacteriophage DNA

In phage propagation, ten milliliters of ϕpp2 phage stock were added to 50 ml of V. parahaemolyticus (3x 108 CFU ml−1) cultured in MSWYE, incubated in a shaker at 37°C for 3–5 hours, when the lysate was clear with some cell debris. The remaining cells and debris were removed by two centrifugations at 10000 × g for 30 minutes. With an optimal titer of 4 × 109 PFU ml−1, the supernatant was stored at 4°C as a phage stock. To concentrate phages using a standard protocol with polyethylene glycol precipitation[2224], solid NaCl (0.6 M) and polyethylene glycol 8000 (20%) were added and precipitation was performed overnight at 4°C. After centrifugation, the phage particles were resuspended in 2 ml of SM buffer and treated with DNase I and RNase A to remove contamination of host nucleotides. The polyethylene glycol was extracted by adding an equal volume of chloroform until the interface was clear. The aqueous phase containing phages was treated with Proteinase K and sodium dodecyl sulfate (SDS) at 56°C for 1 h. Phenol extraction was carried out three times at room temperature; the aqueous phase was further extracted with a 1:1 mixture of equilibrated phenol and chloroform. DNA precipitated by 2× volume of cold ethanol was re-dissolved in deionized water.

Thermal stability of phage ϕpp2

Thermal stability tests have been described elsewhere[25, 26]. Briefly, the bacterium Vibrio parahaemolyticus was freshly inoculated at the 1% volume of seed from overnight culture into 20 ml of rich MSWYE broth. When the cell density reached 0.4–0.5 OD600, the treated phages of a series dilution were added to infect the host for 5 minutes, mixed with top agar, and poured onto a solid surface of regular agar plate in order to count the plaques in 3−5 hours. 2 × 109 PFU of phage particles were treated under 37−80°C and samples were taken at 15-min intervals. The supernatants from the centrifugation of 14000 × g for 3 minutes were diluted and titered for phage numbers by Agar-overlay method.

Genome sequencing and annotation

Similar to shotgun sequencing described elsewhere, approximately 5 μg of the bacteriophage genomic DNA was randomly sheared by nebulization, and DNA sequencing was performed at Mission Biotech according to the manufacturer’s protocol for the Genome Sequencer GS Junior System (Roche Diagnostic). Low quality sequences of the reads generated by the GS Junior sequencer were trimmed off. De novo assembly of the shotgun reads was performed with the GS Assembler software. Sequence assembly and analyses were performed essentially as described previously[27]. Protein-coding genes (PEG) were predicted using The RAST Server (Rapid Annotations using Subsystems Technology;http://rast.nmpdr.org/)[28] and analyzed with the SEED- Viewer (http://www.theseed.org/wiki/Main_Page)[29]. Protein-coding genes were also checked using the ab initio gene-finding program Glimmer v3.02[30]. rRNA genes of the draft assembly were identified using RNAmmer[31]. tRNA genes for all 20 amino acids that were predicted by the RAST were further verified using tRNAscan-SE[32]. Automatic functional annotation results obtained by the RAST were further compared with the proteins in the GenBank database using PSI-BLAST (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi). The Neural Network Promoter Prediction (NNPP) program was used to find the promoters[33].

Multiple sequence alignments

To determine the taxonomy status of the new phage isolate ϕpp2, the genome sequence data of Enterobacteria phage T4 and Vibrio phages KVP40 were employed to find the high homologous regions with the new phage after PSI-BLAST searches. Complete genome sequences of the Vibrio T4-like phages were acquired from NCBI, including Enterobacteria phage T4 (168903 bp in GenBank access no. NC_000866), Vibrio phage KVP40 (244834 bp in GenBank access no. NC_005083), Aeromonas phage 65 (235229 bp in GenBank access no. NC_015251), Aeromonas phage Aeh1 (233234 bp in GenBank access no. NC_005260), and Prochlorococcus phage P-SSM2 (252401 bp in GenBank access no. NC_006883). T4-like myoviruses also include Enterobacteria phages RB14 (NC_012638), RB16 (NC_014467), RB32 (NC_008515), RB51 (NC_012635), JS10 (NC_012741), and JSE (NC_012740), Aeromonas phages 325(NC_008208), and Vibrio cholerae phage ICP1 (NC_015157). PBCV-1 is the Paramecium bursaria Chlorella virus 1. Sequences of individual target genes retrieved from the genome sets were then aligned using ClustalW with default options[34]. The best alignments of individual genes were analyzed by a neighbor-joining method using the NEIGHBOR program in Phylogeny Inference Package (PHYLIP)[35]. Distances were calculated using the PROTDIST programs of PHYLIP and displayed in TreeView[36]. The ClustalW, PHYLIP, and TreeView were bundled in the BioEdit program version 7.0.5[37].

Results

Phage morphology

The morphology of phage ϕpp2 was observed by transmission electron microscopy, which is traditionally one of the most frequently used methods to classify phages. As Figure1 shows, ϕpp2 was a large phage with nonenveloped head, neck, collar, and tail: the head was approximately 90–95 nm wide by 150–160 nm long and the tail was about 110–120 nm long with 20–25 nm in diameter. A baseplate and tail pins were observed under different focus, while long tail-fibers were threading randomly.

Figure 1
figure 1

Transmission electron micrograph of phage ϕpp2 particles with several structural proteins. The phage particles were purified with three times of centrifugations by PEG-NaCl precipitation method mentioned in the Material section. Virion particles were negatively stained with uranyl acetate for EM. The bars represent a length of 100 nm.

Host range

The susceptibility of seven Vibrio strains to the phage ϕpp2 was also investigated with the Agar-overlay method. Among them, V. parahaemolyticus, V. damsela and V. alginolyticus were found susceptible to phage ϕpp2 while the other four species (V. carchariae, V. damsela, V. harveyi, V. pelagius, and V. vulnificus) could not be infected even at high MOI.

Viability of phage ϕpp2 in the thermal environment

Thermal stability test was carried out to analyze the heat-resistant capability of phage ϕpp2 at pH7.5–8.0. The phage was incubated at 37, 50, 61, 70, and 80°C for one hour, respectively. As Figure2 shows, the phage titers at different time intervals demonstrated that phage ϕpp2 stock solution retained almost 100% infection activity after incubation at temperatures lower than 37°C for one hour. When the temperatures rose above 50°C, viability of phage ϕpp2 declined; about 60% phages remained alive after being heated for 60 minutes. At temperatures over 60°C, nearly all phages were inactivated after 15 minutes of incubation.

Figure 2
figure 2

Thermal stability tests of the phage ϕpp2. Samples were taken at different time intervals to titer the phage particles of infectivity.

Genome organization and annotation

The genome sequence of Vibrio phage ϕpp2 was determined using the Roche Genome Sequencer system (454 Life Sciences, Branford, CT). A total of 21,452 reads and 7,985,781 bases, with an average length of 372.3 bases, were obtained. After de novo assembly among at least 40 nucleotide overlap with minimum overlap identity of 90%, the whole genome was aligned to one single contig, with coverage of 32-fold and the Q40 Plus Bases of 98.89% (where Q40 represents an error rate of 99.99%). Currently, the draft genome has a total of 246,421 bp, which includes 270 nt of Q39 Minus Bases (0.11%). The GenBank access number for this new genome is assigned to be JN849462.

The genome size of the Vibrio phage ϕpp2 is 1587 bp larger than 244,834 bp of KVP40 bp and far bigger than the 168,903 bp of T4, while its average G + C content was 42.55%, which is the same as the 42.60% of KVP40 but not as the 35.3% of T4. No rRNA genes of the draft assembly were identified using RNAmmer. Sixty tRNA overlap genes that were preliminarily predicted by the RAST were further verified to be 30 using tRNAscan-SE. In annotation for protein coding regions, 30 subsystem features were predicted by the SEED-RAST server, including 15 features which were relevant to phage structure proteins, 2 for phage DNA synthesis, 7 for nucleotide reactions, and one each for fluoroquinolone resistance, protein degradation and RNA metabolism. One possible gene for resistance of beta-lactamase was not included by the auto-annotation.

Large indels (insertion/deletions)

Overall of Vibrio phage ϕpp2 was similar to the genome organization of vibriophage KVP40 and Enterobacteria phage T4 (1). In comparison with KVP40, 15 deletions and 19 insertions were found in ϕpp2, of which 25 indels only affected one single ORF. It is noteworthy that a single deletion occurred in the seg D-type HE (PEG145 of KVP40), at the junction of KVP40.0145 (at 84923.85078) and KVP40.0146 (complement 85073.85768), implying that ϕpp2 had lost this HE. Most of the indel sizes were in the range around 100–400 nt; nevertheless, some large replacements existed, i.e., 621 nt at KVP40.0102 (61372.61992), 702 nt at KVP40.0121 (70639.71307), 687 nt at KVP40.0147 (85926.86240), 664 nt at KVP40.0172 (98546.98713), 672 nt at KVP40.0277 (nrd A, 146553.148778), and 693 nt at KVP40.0315 (178766.178930). Additionally, three KVP40 genes were replaced by giant inserts in ϕpp2: 1098 nt of ϕpp2 replaced the gene near KVP40.0363 (gp23, 224506.226050, 1545 nt), 1143 nt of ϕpp2 replaced the gene at KVP40.0263 (137878.138114, 237 nt), and 3330 nt of ϕpp2 replaced the gene at KVP40.0297 (complement 160413.160988,576 nt). The three giant indels signified that the Vibrio phage ϕpp2 was a new species from KVP40.

Table 1 Gene functions of the Vibrio phage ϕpp2

Gene functions

With the extracting plausible protein sequences encoded by the genomic DNAs, 383 PEGs were found in Vibrio phage ϕpp2, in contrast to 381 PEGs for KVP40 found with the same RAST method (1). Functions were identified by sequence similarity (1). 104 (27.2%) out of 383 PEGs were matched to known functions of T4-like phage genes and assorted bacteria genomes, while functions of 279 (72.8%) PEGs were still unknown. Among these, as Figure3 and1 show, 67 PEGs were matched to both T4 and KVP40 (green arrows), 29 PEGs to KVP40 alone (yellow ticks), 7 PEGs to other T4-like (purple and red ticks), and one to assorted bacteria (cyan). Between the genomes of ϕpp2 and KVP40, the similarity of 254 genes was greater than 94%, whereas 17 PEGs of KVP40 and 21 PEGs of ϕpp2 were unmatched to any known, in addition to 15 genes with lower similarity ( Additional file1). At least 29 PEGs (7.6%) were directly related to phage particle structures, such as head, tail, and baseplate. ϕpp2 uniquely possessed the proteins involved in potassium channel, gp2 (DNA end protector), tRNA nucleotidyltransferase, and mob-type HEs, which were not reported in the case of KVP40. Several genes were split: in ϕpp2, PEG297 shared paralogs with PEG296, as the same pattern for PEG119 sharing with PEG274, while KVP40.0089 (54956.55189) and KVP40.0090 (55200.56117) paralogs were matched to one single ϕpp2 PEG88.

Figure 3
figure 3

Genome map of the Vibrio phage ϕpp2. Green arrows indicate the genes matched to both Enterobacteria phage T4 and Vibrio phage KVP40. Yellow ticks on the circle indicate that the genes fitted to KVP40 only while the yellow triangle indicates the absent site for KVP40.0146 HE gene. Purple represents that the genes only aligned well with T4. The cyan is for one gene matched to GTP cyclohydrolase I from Bdellovibrio bacteriovoru HD100, Vibrio angustu S14, and Cytophaga hutchinsonii ATCC 33406. Red bars with the number indicate the PEG numbers of potential HE.

Transfer RNAs

The RAST predicted 60 pieces of potential tRNAs, spanning in the range of 9175 bp in Vibrio phage ϕpp2, while in KVP40 29 tRNAs were found in the range of 8702 bp. Using tRNAscanSE to recalculate the structures with overlapping sequences, 30 tRNAs in the cluster were double verified for Vibrio phage ϕpp2 while 29 tRNAs remained for KVP40; both contained three pseudo-forms of low score for GCA (two) and TGC (one) anticodons. The Vibrio phage ϕpp2 tRNA cluster encoded for 17 amino acid codons, but there were no anticodons for alanine, glutamine, and tyrosine. The KVP40 tRNA region was 475 bp shorter than the ϕpp2 but shared 97% similarity over the cluster. A big insert of 465 nt in the middle of the cluster created no putative tRNA structure in the range of insert. In the Vibrio phage ϕpp2, one extra met-tRNA, which formed from the 28 nt mutation out of 72 nt, was created at the upstream of junction that was 6 nt upstream from the 465-nt insert.

Searching mob-like genes and neighbors

In sequence similarity analysis by PSI-BLAST, three paralog genes of homing endonucleases were found in the Vibrio phage ϕpp2: PEG79, PEG119, and PEG274, in which the number of amino acid residues was 209, 234, and 224 aa, respectively. The PEG119 and PEG274 were aligned to neighborhood of T4 MobE and close to MobD (Figure4A). The PEG79 were situated next to the group of MobA (Figure4A). The PEG119 shared 37% similarity with PEG274, while PEG79 shared 27% and 35% with the other two in pair-wide alignment of amino acid sequences. In Bootstrap analysis with 1000 replicates, the branch percentage showed that the three PEGs in ϕpp2 were all Mob-like homing endonucleases, least likely to be a GIY-YIG type (Figure4B). Although low overall similarity was found between them, all three PEGs aligned the H-N-H motif very well in their N-termini (Figure4C). First two His-32 and His-33 in PEG79 were highly conserved within the motif of ExHH ILPK for PEG119 and PEG274. The second Asn-50 of PEG79 was situated in the motif of SDExN LV, and the third His was paired as HxxxH found in the motif of LTAREH---H xLLxK.

Figure 4
figure 4

Phylogenetic analyses and similarity of the HE genes from different T4-related phages. ϕpp2 is a Vibrio phage isolated in this study. PEG numbers without dash are Enterobacteria phage T4. The homing endonucleases are named with gene product numbers followed by the dash lines for the hosts of the phages: Enterobacteria phages include RB14, RB16, RB32, RB51, JS10, and JSE; Aeromonas phages, Aeh1, 25 and 65; PBCV-1 is Paramecium bursaria Chlorella virus 1; and Vibrio phage ICP1. (A) Rooted phylogenetic tree for the homing endonucleases of Vibrio and T4-like phages by PROTDIST-neighbor joining method; the amino acid sequences were aligned with BLOSUM62 matrix, gap penalty = 8 and extension penalty = 2. (B) Bootstrap analysis for the Mob-type HEs of the Vibrio phages against T4 phages. The bootstrap values of percentages in 1000 replicates are placed on the branch for the nodes defining each monophyletic clade. The scale bars represent distance length. (C) H-N-H alignment of three HE genes from ϕpp2 with T4 mobE and ICP1 ORF28 (a phage in Vibrio cholerae).

We identified the Mob types of ϕpp2 around the genome according to the orientation similarity to the neighbor ORFs of 15 homing endonucleases in Enterobacteria phage T4: gt (glucosyl transferase) and nrd (ribonucleotide reductase) orthologs. The details of the search methods are described in Additional file2. No match to T4 α.gt or β.gt was found in entire genomes of ϕpp2 and KVP40 (NC_005083); therefore, the mob B-like gene could not exist in ϕpp2. Four nrd-like genes were found in ϕpp2: one was found explicitly by the RAST and three others were implicit but manually confirmed with PSI-BLAST searches. The PEG12 protein was similar to the large subunit of anaerobic ribonucleotide reductase of class III (EC 1.17.4.2), with 52.05% similarity to T4 nrd G, while PEG15 was assumed to be the activating protein (EC 1.97.1.4) for the ribonucleotide reductase with 52.74% similarity to T4 nrd D. PEG132, which was matched to T4p232 in the boundary of MobE and downstream close-by seg D, was denoted as nrd B.1. The fourth nrd-like gene, 1041 nt of PEG148 (347 aa, 89176.90216) in ϕpp2, was mapped to T4 nrdC.11.

Using the neighbor-indirect method (details in Additional file2), the neighbors of T4 mob genes were mapped back to ϕpp2 genome. The neighbor gene T4p074 (nrd G) of mob C (T4p075) was back-projected to ϕpp2 PEG15 with a similarity of 52.05%; another neighbor gene T4p076 (nrd D) was matched to ϕpp2 PEG12 with a similarity of 52.74%. The distance of the PEG12/15 pair was at least 37874 nt apart from the PEG79 – it was even farther to PEG 119 and PEG274. Similarly, the PEG132 and PEG148 were still too far to be adjunction neighbors for all three potential ϕpp2 mob genes, i.e., PEG132 was 7895 nt apart from PEG119.

Alternatively, using the so-called neighbor-direct method ( Additional file2), the mob-neighbor genes of ϕpp2 PEG79, PEG119, and PEG274 were manually de novo searched with PSI-BLAST. Neither neighbors of PEG79 (upstream PEGs 70 ~ 78 and downstream PEGs 80 ~ 90) nor of PEG119 (upstream PEGs 110 ~ 117 and downstream PEGs 120 ~ 125) were in any way close to nrd-like genes (Figure5A and B).

Figure 5
figure 5

The best aligned T4-related phage genes for the neighbor genes of ϕpp2 and KVP40 HEs using the neighbor-direct method. The approach is described in the text and Additional file2. The same color arrows represents the homologous genes. The cyan arrows indicate the HE genes for ϕpp2 and KVP40. (A) Neighbors of ϕpp2 PEG79: 78 for RegA translational repressor of early genes, 80 for phage hypothetical protein, and 81 for DNA polymerase. (B) Neighbors of ϕpp2 PEG119: 118 for GTP cyclohydrolase I from Bdellovibrio bacteriovorus HD100, Vibrio angustum S14, and Cytophaga hutchinsonii ATCC 33406; 120 for phage hypothetical protein. (C) Neighbors of ϕpp2 PEG274: 273 & 275 for Nrd, ribonucleotide reductase Ia; 276 for NrdC thioredoxin. (D) Neighbors of KVP40.146: 143 for rIIA protector; 144 for rIIB protector.

Vibrio phage ϕpp2 PEG274 with mob E-type neighbors

In de novo identification of a mob-type for ϕpp2 PEG274 (672 nt) using the neighbor-direct method, ϕpp2 PEG273 (2226 nt) of the upstream neighbor gene was blasted to NrdA of Aeromonas phages (PX29 and phiAS5), Enterobacteria phages (JSE, RB49, phi1, and T4), and Shigella phage SP18 (Figure5C). The downstream neighbor PEG275 (1125 nt) was blasted to NrdB of Aeromonas phages phiAS5 and Aeh1, Klebsiella phage KP15, and Enterobacteria phage RB16. Another neighbor, PEG276 (300 nt), was also blasted to the NrdC thioredoxin; it aligned well as 86% homologous to NrdC thioredoxin in Aeromonas phages phiAS5, Aeh1, and 65, as well as to Klebsiella phage KP15, Shigella phage SP18, and Enterobacteria phages RB16, RB43, and ime09. With the matches of upstream and downstream of nrd-like genes which complemented the full organization of MobE neighbors, the ϕpp2 PEG274 can be annotated as a MobE-type HE, without the existence of I-Tev III intron.

Expression of ϕpp2 PEG 274 gene

All homing endonucleases of ϕpp2 and KVP40 started with an AUG initiation codon. For ϕpp2 PEG 274, AGGA as a ribosome binding site (RBS) was optimally situated 6 nt upstream of the PEG start codon while translation initiation regions are not positioned at the optimal distance of 6–9 nucleotides from the AUG codon for PEG79 and PEG 119. The AAGAGAG for PEG79 was not a good match to antisense of small rRNA while the predicted PEG119 RBS AGGA is immediately adjunct to the AUG codon, which is rarely considered to be a good initiation site for translation.

As shown in Figure3, the direction of transcription for PEG79 and PEG119 was counter-clockwise, whereas the PEG274 promoter was clockwise, which was the same as most T4 homolog genes. The NNPP predicted several promoters of high scores (>0.95) in the upstream of these three genes. It is worth noting that three promoters were identified around PEG274, the aforementioned MobE-type homing endonuclease. In contrast to the translational initiation AUG position of PEG274 at 149293–149964 in the genome of ϕpp2, the nearby promoters were also positioned at 148783 (510 nt upstream; pR148783), 149272 (21 nt immediately upstream; pR149272), and 149974 (10 nt downstream; pR149974). pR149272 was the best fit to the promoter consensus, which consisted of TTGTGA for −35 box and ATGTAAAAT for −10 box. Accompanying this promoter, some weak binding sites for transcription factors were also observed: TGTAAAAT for rpo D17 at position 149258, ATATAAAT for arg R2at 149264, and GTTCATAT for tor R at 149273.

Discussion

Electron microscopy revealed that the phage ϕpp2 particles were morphologically similar to T4 phage and vibriophage KVP40, which is a long head (~140 nm long and ~70 nm wide) with a prolate icosahedral capsid and a contractile tail with associated baseplate and extended tail fibers. ϕpp2 is most likely type A phage in Bradley’s classification of Myoviridae[38], based on the morphological characteristics (Figure1). The protein profiles in ϕpp2 contain a heavy band of ~50 kD, which is similar to known T4 structure proteins of major capsid protein (data not shown). With hourly heat-tolerance at 50°C (Figure2), this phage could infect aquaculture pathogens, V. parahaemolyticus, V. damsela and V. alginolyticus. The complete genome of the new Vibrio phage ϕpp2 was sequenced (GenBank access_no JN849462), which was a sibling phage of KVP40 but with different HE genes (Figure3).

In the phylogenetic tree (Figure4A), the PEG79 was distantly situated next to the group of MobA. Although their overall similarity was low, the N-termini of all three PEGs aligned well with the H-N-H motif (Figure4C). The first His-32/33 in PEG79 was highly conserved within the motif of ExH HILPK for PEG119 and PEG274. The second Asn was situated in the motif of SDExN LV and the third His-pair was in the paired form of HxxxH found in the motif of LTAREH---H xLLxK. This reveals that the ϕpp2 HE genes belong to Mob-type because the H-N-H is the critical motif for the enzyme activity[10]. The vibriophage KVP40 carries segD/C (KVP40.0146)[8]. V. cholerae ORF80 in ICP1 belongs to segG (data not shown) while another ICP1-ORF28 is closely related to MobC (Figure5B and C)[9].

By PSI-BLAST directly from the neighbor genes of ϕpp2 PEG274 (the neighbor-direct method), PEG273, PEG275, and PEG276 were highly homologous to NrdA, NrdB and NrdC thioredoxin in Aeromonas phages, Enterobacteria phages, Klebsiella phage KP15 and Shigella phage SP18, respectively. With match of both up- and downstream, together with the conserved motif of HE in Figure4C, the PEG274 can be annotated as MobE-type HE. For PEG274 protein expression, we found a good promoter (pR149272) immediately upstream of the PEG274 gene; thus, the promoter was considered as endonuclease-specific. The transcript of PEG274 mRNA was also equipped with a good consensus of ribosome binding site AGGA at 6 nt upstream of the start codon AUG.

Sequence of ϕpp2 PEG79 was comparatively similar to MobA gene, but PEG79 was flanked by DNApol and reg A (phage endoribonuclease translational repressor of early genes; Figure5A), where they do not neighbor any mob genes in T4. The PEG119 and PEG79 genes were similar to T4p232 and T4p233 (mob E), respectively. The landmark of T4p131 (e.8, complement 70360.70623) is also very similar to PEG275. In other words, three ϕpp2 mob-like genes (PEG79, PEG119, and PEG274) would be mapped onto the cluster of I-Tev III-nrd B1-mob E located at T4p130 to T4p133 in T4 genome[16], a review]. This implies the characteristics of HE mobility.

KVP40.0146 (696 nt) encoding 231 aa was PSI-BLAST to GIY-YIG endonuclease genes, including Aeromonas phages (phage 25 and phiAS5), Acinetobacter phages (Acj61 and Ac42), Chlorella virus FR483, Enterobacteria phages (RB51, RB16, T4), Klebsiella phage KP15, and Staphylococcus phage PH15. As shown in Figure4A, the phylogenetic analysis plotted KVP40.0146 to be a seg C/D type. Using the neighbor-direct method ( Additional file2), KVP40.0144 and KVP40.0145 could not be matched to any protein of known function (Figure5D) while KVP40.142 and KVP40.143 could be similar to rIIA/B lysis protectors. Both were too distant to bracket the KVP40.0146 of GIY-YIG endonuclease gene for mimicking the T4 segD neighbor. In T4 HEs, types of mob C, mob D, and mob E can be classified by neighbor elements as well as different arrangements of their promoters: nrdD-mobC-nrdG, mobD-nrdC.11, and nrdA-(I-Tev III)-mobE-nrdB, respectively[16]. In KVP40, there are seven nrd-like genes that have been identified: nrd A, B, C, C.11, D, G, and H. The closest one for KVP40.0146 HE was nrdC.11 (KVP40.0153; 88930.89970), but it was still too distant to be a neighbor of KVP40.0146 to form a good setting as the T4 mob C/D/E. Similarly, four nrd genes were found in Vibrio cholerae ICP1 but without any HE insertion. Therefore, KVP40 and ICP1 did not have the same organization of T4 HEs.

KVP40, sharing the same host as ϕpp2, has only one putative seg C/D-type KVP40.0146 (complement 85073.85768), which was also similar in part to T4 seg B/E and I-Tev III, even nrd B.1[9]. Therefore, the two giant Vibrio phages could partially cross the boundary line at nrd B.1 (Figure4A), in the same host of V. parahaemolyticus, to catch the genes and evolve for the future form as the Enterobacteria phage T4 did. The mechanism for the gene exchange and/or evolution may also be similar to the PEG79, PEG119, and PEG 274 in the ϕpp2 as mentioned above.

Conclusions

In summary, the phage ϕpp2 was characterized by the morphology, growth, and genomics. In the complete genome sequence analysis in this study, three giant indels and the mob E-type HE signified the Vibrio phage ϕpp2 to be a new species of T4-related phages, different from KVP40. Our analysis suggested that ϕpp2 PEG274 was an active mob E gene with transcriptional and translational elements. In the same host, Vibrio parahaemolyticus, the new phage ϕpp2 can complement its mob-type HE functions with KVP40 that only carries a seg-type HE gene. This spectrum of genome datasets of T4-related Vibrio phages that can co-infect the same host will be useful to investigate the hypothesis that a lateral transfer of freestanding HEs with self-mobility may result in genomic mosaicism by recombining a variety of genetic sequences in phage genomes[18].

References

  1. Daniels NA, MacKinnon L, Bishop R, Altekruse S, Ray B, Hammond RM, Thompson S, Wilson S, Bean NH, Griffin PM, Slutsker L: Vibrio parahaemolyticus infection in the United States, 1973–1998. J Infect Dis. 2000, 181: 1661-1666. 10.1086/315459.

    Article  CAS  PubMed  Google Scholar 

  2. Fujino T, Okuno Y, Nakada D, Aoyama A, Fukai K, Mukai T, Ueho T: On the bacteriological examination of Shirasu-food poisoning. Med J Osaka Univ. 1953, 4: 299-304.

    Google Scholar 

  3. Levine WC, Griffin PM: Vibrio infections on the Gulf Coast: results of first year of regional surveillance. J Infect Dis. 1993, 167: 479-483.

    Article  CAS  PubMed  Google Scholar 

  4. Morris JG, Black RE: Cholera and other vibrioses in the United States. N Engl J Med. 1985, 312: 343-350. 10.1056/NEJM198502073120604.

    Article  PubMed  Google Scholar 

  5. Johnson DE, Weinberg L, Ciarkowski J, West P, Colwell RR: Wound infection caused by Kanagawa-negative Vibrio parahaemolyticus. J Clin Microbiol. 1984 October, 20 (4): 811-812.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. Alam M, Chowdhury WB, Bhuiyan NA, Islam A, Hasan NA, Nair GB, Watanabe H, Siddique AK, Huq A, Sack RB, Akhter MZ, Grim CJ, Kam K-M, Luey CKY, Endtz HP, Cravioto A, Colwell RR: Serogroup, Virulence, and Genetic Traits of Vibrio parahaemolyticus in the Estuarine Ecosystem of Bangladesh. Appl Environ Microbiol. 2009, 75: 6268-6274. 10.1128/AEM.00266-09.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Marano NN, Daniels NA, Easton AN, McShan A, Ray B, Wells JG, Griffin PM, Angulo FJ: A survey of stool culturing practices for Vibrio species at clinical laboratories in Gulf Coast states. J Clin Microbiol. 2000, 38: 2267-2270.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Miller ES, Heidelberg JF, Eisen JA, Nelson WC, Durkin AS, Ciecko A, Feldblyum TV, White O, Paulsen IT, Nierman WC, Lee J, Szczypinski B, Fraser CM: Complete genome sequence of the broad-host-range vibriophage KVP40: comparative genomics of a T4-related bacteriophage. J Bacteriol. 2003, 185 (17): 5220-5233. 10.1128/JB.185.17.5220-5233.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Seed KD, Bodi KL, Kropinski AM, Ackermann H-W, Calderwood SB, Qadri F, Camilli A: Evidence of a Dominant Lineage of Vibrio cholerae-Specific Lytic Bacteriophages Shed by Cholera Patients over a 10-Year Period in Dhaka. Bangladesh. mBio. 2011, 2 (1): e00334-10.

    PubMed  Google Scholar 

  10. Corina LE, Qiu W, Desai A, Herrin DL: Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs. Nucleic Acids Res. 2009, 37 (17): 5810-5821. 10.1093/nar/gkp624.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Wilson GW, Edgell DR: Phage T4 mobE promotes trans homing of the defunct homing endonuclease I-TevIII. Nucleic Acids Res. 2009, 37 (21): 7110-7123. 10.1093/nar/gkp769.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Mueser TC, Hinerman JM, Devos JM, Boyer RA, Williams KJ: Structural analysis of bacteriophage T4 DNA replication: a review in the Virology Journal series on bacteriophage T4 and its relatives. Virol J. 2010, 7: 359-375. 10.1186/1743-422X-7-359.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Landthaler M, Shub DA: The nicking homing endonuclease I-Bas I is encoded by a group I intron in the DNA polymerase gene of the Bacillus thuringiensis phage Bastille. Nucleic Acids Res. 2003, 31: 3071-3077. 10.1093/nar/gkg433.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Sharma M, Hinton DM: Purification and characterization of the SegA protein of bacteriophage T4, an endonuclease related to proteins encoded by group I introns. J Bacteriol. 1994, 176: 6439-6448.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Tseng MJ, He P, Hilfinger JM, Greenberg GR: Bacteriophage T4 nrd A and nrd B genes, encoding ribonucleotide reductase, are expressed both separately and coordinately: characterization of the nrdB promoter. J Bacteriol. 1990, 172: 6323-6332.

    PubMed Central  CAS  PubMed  Google Scholar 

  16. Edgell DR, Gibb EA, Belfort M: Mobile DNA elements in T4 and related phages. Virol J. 2010, 7: 290-304. 10.1186/1743-422X-7-290.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Eddy SR, Gold L: The phage T4 nrd B intron: a deletion mutant of a version found in the wild. Genes Dev. 1991, 5 (6): 1032-1041. 10.1101/gad.5.6.1032.

    Article  CAS  PubMed  Google Scholar 

  18. Petrov VM, Ratnayaka S, Nolan JM, Miller ES, Karam JD: Genomes of the T4-related bacteriophages as windows on microbial genome evolution. Virol J. 2010, 7: 292-320. 10.1186/1743-422X-7-292.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Cabrera-García ME, Vázquez-Salinas C, Quiñones-Ramírez EI: Serologic and molecular characterization of Vibrio parahaemolyticus strains isolated from seawater and fish products of the Gulf of Mexico. Appl Environ Microbiol. 2004, 70 (11): 6401-6406. 10.1128/AEM.70.11.6401-6406.2004.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Lin Y-R, Chiu C-W, Chang F-Y, Lin C-S: Characterization of a new phage, termed ϕA318, which is specific for Vibrio alginolyticus. Arch Viol. 2012, 10.1007/s00705-012-1244-8. xx (in press; Feb 11, 2012)

    Google Scholar 

  21. Lu M-W, Liu W, Lin C-S: Infection competition against grouper nervous necrosis virus by virus-like particles produced in Escherichia coli. J Gen Virol. 2003, 84: 1577-1582. 10.1099/vir.0.18649-0.

    Article  CAS  PubMed  Google Scholar 

  22. Sambrook J, Russell DW: Molecular cloning: a laboratory manual. 2001, Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press, 3

    Google Scholar 

  23. Shafia F, Thompson TL: Calcium ion requirement for proliferation of bacteriophage Phi Mu-4. J Bacteriol. 1964, 88: 293-296.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Adams MH: Bacteriophages. 1959, New York: Interscience

    Google Scholar 

  25. Abedon ST, Herschler TD, Stopar D: Bacteriophage latent-period evolution as a response to resource availability. Appl Env Microbiol. 2001, 67: 4233-4241. 10.1128/AEM.67.9.4233-4241.2001.

    Article  CAS  Google Scholar 

  26. Mitra S, Basu S: Some biophysical properties of a vibriophage and its DNA. Biochim Biophys Acta. 1968, 155 (1): 143-149. 10.1016/0005-2787(68)90344-4.

    Article  CAS  PubMed  Google Scholar 

  27. Lind PA, Andersson DI: Whole-genome mutational biases in bacteria. Proc. Natl. Acad. Sci. USA. 2008, 105 (46): 7878-17883.

    Article  Google Scholar 

  28. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics. 2008, 9: 75-89. 10.1186/1471-2164-9-75.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, de Crécy-Lagard V, Diaz N, Disz T, Edwards R, Fonstein M, Frank ED, Gerdes S, Glass EM, Goesmann A, Hanson A, Iwata-Reuyl D, Jensen R, Jamshidi N, Krause L, Kubal M, Larsen N, Linke B, McHardy AC, Meyer F, Neuweger H, Olsen G, Olson R, Osterman A, Portnoy V, Pusch GD, Rodionov DA, Rückert C, Steiner J, Stevens R, Thiele I, Vassieva O, Ye Y, Zagnitko O, Vonstein V: The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005, 33 (17): 5691-702. 10.1093/nar/gki866.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Delcher AL, Bratke KA, Powers EC, Salzberg SL: Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics. 2007, 23 (6): 673-679. 10.1093/bioinformatics/btm009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Lagesen K, Hallin PF, Roland EA, Stafeldt HH, Rognes T, Ussery DW: RNammer: consistent annotation of rRNA genes in genomic sequences. Nucleic Acids Res. 2007, 35 (9): 3100-3108. 10.1093/nar/gkm160.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Reese MG: PhD Thesis (PDF). Computational prediction of gene structure and regulation in the genome of Drosophila melanogaster. 2000, UC Berkeley/University of Hohenheim

    Google Scholar 

  34. Thompson JD, Higgins DG, Gibson TJ, CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Felsenstein J: PHYLIP--Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.

    Google Scholar 

  36. Page RDM: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.

    CAS  PubMed  Google Scholar 

  37. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.

    CAS  Google Scholar 

  38. Bradley DE: Ultrastructure of bacteriophages and bacteriocins. Bacteriol Rev. 1967, 31 (4): 230-314.

    PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This research fund is partially supported by the grants from the National Science Council, Taiwan (NSC96-2313-B-110-002-MY3 and NSC99-2313-B-110-002-MY3), and the Ministry of Education, Taiwan (NSYSU95 ~ 99C031701; the second term of Top University Program: NSYSU 00C030205 and NCHU 100-S05-09) under the ATU plan. We thank Professor Long-Huw Lee (National Chung-Hsing University) as the grant organizer of intercampus ATU plan, as well as Chi-Wen Chiu and Feng-Yi Chang in helping initial phage screening, Yu-Tin Liu in helping gel preparations, Professor Y. W. Chiang and Mr. S-C Lin for EM operation, and Kenneth B. Lin and Dr. Simon White for comments and editing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chan-Shing Lin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CSL conceived and designed the study. YRL and CSL did the experiments, analyzed the sequence and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lin, YR., Lin, CS. Genome-wide characterization of vibrio phage ϕpp2 with unique arrangements of the mob-like genes. BMC Genomics 13, 224 (2012). https://doi.org/10.1186/1471-2164-13-224

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-224

Keywords