Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2009, p. 2659-2667, Vol. 75, No. 9
0099-2240/09/$08.00+0 doi:10.1128/AEM.02483-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Department of Microbiology, Soochow University, Taipei, Taiwan 111, Republic of China,1 Institute of BioMedical Informatics, National Yang-Ming University, Taipei, Taiwan 112, Republic of China,2 Department of Biotechnology and Laboratory Science in Medicine, National Yang-Ming University, Taipei, Taiwan 112, Republic of China,3 Institute of Botany, Academia Sinica, Taipei, Taiwan 115, Republic of China4
Received 29 October 2008/ Accepted 6 March 2009
|
|
|---|
HAP-1 were highly similar, while N15,
KO2, and PY54 made up another closely related group. The presence of DNA adenine methylase and quorum-sensing transcriptional regulators in VP882 may play a specific role in this phage or regulate physiological or virulence-associated traits of the hosts. These genes may also be remnants from the bacterial chromosome following transduction. |
|
|---|
Since 1996, gastroenteritis that is caused by particular O3:K6 strains (designated as pandemic O3:K6) of V. parahaemolyticus has been widespread in Southeast and East Asia (55), North America (27), South America (20), Europe (34), and Africa (3). It has been regarded as the first pandemic of this pathogen (43). These widespread pandemic O3:K6 strains have been found to be genetically closely related based on various molecular methods and appear to constitute a clone that differs markedly from the nonpandemic O3:K6 strains that were isolated before 1996 (4, 43, 55). The acquisition of the traits associated with these strains that account for the pandemic is of great concern.
In pandemic Vibrio cholerae, virulence factors are encoded in the genome of an integrated filamentous phage, CTX
(53). The integration of the filamentous phage CTX
into the bacterial genome conveys pathogenicity to V. cholerae and to Vibrio mimicus (8). Filamentous phage f237, which is similar to the previously described phage VF33, has been identified in a pandemic O3:K6 V. parahaemolyticus strain (38, 39). ORF8 in this phage is unique; no homologous sequence is present in the databases, and no biological function has hitherto been identified (39). This unique sequence has been used as a genetic marker of the presence of this filamentous phage (25), and the homologous sequence has been demonstrated in other pandemic O3:K6 and closely related strains (12, 37).
The role of this filamentous phage in the virulence of the pandemic O3:K6 strain has not been established. Along with the filamentous phage, the presence of temperate phages (VP882 and VP1092) in pandemic O3:K6 strains 882 and 1092, respectively, was demonstrated in our laboratory. In this study, VP882 was characterized using electron microscopy and genomic sequencing to provide more information on the possible roles of temperate phages on lysogenic conversion and on the diversity of this pathogen. Specific traits of the VP882 phage were identified by comparative genomic analysis with other vibriophages and plasmid-like prophages.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. V. parahaemolyticus, V. cholerae, and V. vulnificus strains used in this study
|
Isolation of phage particles and viral DNA.
The lysed bacterial culture was centrifuged at 5,500 x g for 15 min, and the supernatant was filtered through a 0.22-µm membrane (Millipore Co., Billerica, MA). Sodium chloride was added to the filtrate to a final concentration of 1 M, and the filtrate was stored statically for 30 min on ice. Next, the phage particles were precipitated by adding 10% polyethylene glycol 6000 (Sigma) to the filtrate and incubating it at 4°C for at least 24 h (47). The phage particles were recovered by centrifugation at 5,500 x g for 20 min at 4°C, washed with phosphate-buffered saline containing 3% NaCl, and resuspended in a buffer that contained 10 mM of Tris-HCl and 1 mM of EDTA at pH 8.0, and the viral DNA was isolated and purified by the procedure of Sambrook et al. (47).
Electron microscopy.
The phage particles were fixed with 5% formaldehyde, washed in distilled water, and then negatively stained with 2% phosphotungstic acid, before being examined using a Phillips CM100 (FEI Company, Hillsboro, OR) transmission electron microscope (19).
Molecular techniques.
Quantification, electrophoresis, enzyme restriction, nucleotide sequencing, and other general manipulations of DNA were performed by the standard procedures as described by Sambrook et al. (47). The size of the phage genome was determined by pulsed-field gel electrophoresis (PFGE), following published methods (55). Briefly, the purified genome DNA was resolved using a CHEF II apparatus (Bio-Rad) at 200 V with a 0.5- to 1.5-min pulse for 18 h. Lambda DNA monocut mix (New England BioLabs, Beverly, MA) was used as a size marker in PFGE. 2-Log DNA ladder (New England BioLabs) and HindIII-digested lambda DNA were used as size markers in agarose gel electrophoresis.
The primers used for PCR were designed using Primer Design Assistant software (14). PCR was conducted using a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA) as follows: an initial step of 2 min at 94°C; followed by 30 cycles, with 1 cycle consisting of 1 min at 94°C, 1 min at 60°C, and 2 min at 72°C; and a final step of 10 min at 72°C.
Whole-genome shotgun sequencing.
Genome sequencing was conducted at the VGH National Yang-Ming University Genome Research Center by the shotgun sequencing method (13). A 10-fold genome coverage was achieved using procedures briefly described herein. The genomic DNA of phage VP882 was fragmented by sonication, the single strands were digested by BalIII, and blunt-end fragments were made using T4 DNA polymerase and separated by agarose gel electrophoresis. Fragments (2.5 to 3.0 kilobase pairs [kbp]) were eluted, digested by SmaI, ligated to the pUC18 vector, and transformed into Max Efficiency DH5
competent cells (Invitrogen, Carlsbad, CA). Ampicillin-resistant transformants were selected and sequenced using a BigDye terminator cycle sequencing ready reaction kit (Applied Biosystems, Foster City, CA) with an automated sequencer (Applied Biosystems 3730). Sequences were jointly assembled using Phred/Phrap/Consed software (obtained from the University of Washington, Seattle). The accuracy, order, and orientation of contigs were examined based on linking information from forward and reverse sequence ends of each clone. Sequence gaps were closed by editing the end sequences of each contig, by primer walking on linking clones, and by sequencing PCR products from shotgun clones (13).
Gene prediction and annotation.
A combination of three gene prediction programs, Glimmer (version 3.02), GeneMark HMM (version 2.0), and Zcurve_V (version 1.0), were used to predict the protein-coding regions (17, 22, 32). The final designation of protein-coding sequence (CDS) is a unified set of the three prediction results. A protein sequence homology search of the predicted coding regions was performed using NCBI BLAST (version 2.2.15) against NCBI nonredundant protein, KEGG genes, and Uniprot Trembl protein databases to assign a hypothetical role or function to every CDS and thus generate a putative gene map of VP882 genome. CDSs with no blast hit were designated as unknown, and those with blast hits but unknown function were designated as putative phage-related proteins.
Bioinformatics analysis.
Inverted repeat sequence (IRS) analysis of the six plasmid-like prophage genomes was carried out with the Inverted Repeat Finder program (version 3.05) (54). Nucleotide sequence comparison of the six prophage genomes was analyzed with the DOTTER program (50). These prophage genomes included the Vibrio phage VP882 (GenBank accession no. EF057797), Halomonas phage
HAP-1 (NC_010342), Vibrio phage VHML (NC_004456), Yersinia phage PY54 (NC_005069), enterobacterial phage N15 (NC_001901), and Klebsiella phage
KO2 (NC_005857) (NC numbers in parentheses are GenBank nucleotide database accession numbers). IRSs were aligned and visualized by using the Weblogo program (version 2.8.2) (16). The multiple-sequence alignment program Clustalw2 (29) was used to align protelomerase protein sequences. The aligned sequences were taken from prophages VP882 (YP_001039865),
HAP-1 (YP_001686770), VHML (NP_758894), PY54 (NP_892077), N15 (NP_046924), and
KO2 (YP_006606) (YP and NP numbers are GenBank protein database accession numbers). Comparison of the important residues of the protelomerase catalytic domain was also performed.
Plaque formation.
The susceptibility of various strains of V. parahaemolyticus, V. cholerae, and V. vulnificus to phage VP882 was determined by spot assay. Fifty microliters of bacterial culture (about 104 CFU/ml) was spread on different agar plates (TSB with 3% NaCl for V. parahaemolyticus, LB with 2% NaCl for V. cholerae, and BHI with 0.85% NaCl for V. vulnificus) on which 20 µl of viral suspension (about 104 PFU/ml) was spotted. The plates were incubated at 37°C for 16 h, and the characteristics of the plaques were observed.
Nucleotide sequence accession number.
The genome sequence of the VP882 phage was deposited in the GenBank database under accession number EF057797.
|
|
|---|
Electron microscopy revealed that the VP882 and VP1092 particles were morphologically similar; only VP882 was further characterized in this work. VP882 phage exhibited a polyhedral head (48 ± 2 nm) with a rigid long tail (156 ± 8 nm), which was loosely associated with a special sheath-like structure (length, 67 ± 5 nm; width, 20 ± 1 nm) (Fig. 1). VP882 was identified as a Myoviridae phage based on the morphological characteristics (a polyhedral head, a neck/collar region, and a sheathed rigid tail) (Fig. 1). The tail sheath of VP882 was probably loosely associated with the tail and attached to it at various positions or even dislodged in the viral particles (Fig. 1). The heads of VP882 (48 nm) were smaller than those previously described for Myoviridae, and they had diameters of 53 to 160 nm (21, 41, 46). The shape and size of VP882 were similar to those of the VHML phage of Vibrio harveyi (41), while VP882 had a shorter tail than
HAP-1, a similar marine phage, did (36).
![]() View larger version (163K): [in a new window] |
FIG. 1. Electron micrograph of the VP882 particles. Bar, 50 nm.
|
Phages are commonly found in marine environments and have been isolated from 28% of marine molluscan shellfish, crustaceans, seawater, and sediments. The predominant phage types are specific for some strains of V. parahaemolyticus and some other agar-digesting and psychrophilic Vibrio species (6). These marine phages do not infect Vibrio alginolyticus or V. cholerae (5, 6). Recently, a survey of the water column in northwestern America revealed that V. parahaemolyticus phages, mostly siphoviruses, are frequently isolated and infect high percentages of the host strains, especially those isolated from oysters and sediments, and half of these phages are not species specific and can infect V. alginolyticus strains (15). This investigation demonstrated that VP882 was probably a marine phage that exhibits interspecies infection in marine vibrios and V. cholerae, which proliferates well in freshwater and grows well in lake water samples amended with different concentrations of NaCl (52).
Genome properties.
The size of the genomic DNA molecule was determined by PFGE to be about 38 kbp. The genome sequencing of phage VP882 revealed that the genome size was 38,197 base pairs (bp) (GenBank accession no. EF057797), which was close to the size determined by PFGE. The G+C percentage was 56.95%.
Analysis of the genome sequence identified 71 putative open reading frames (ORFs). Based upon translated sequence homology, 27 are putative phage or bacterial genes with hypothetical functions; 15 are putative phage-related genes or other genes with unknown functions (ORF5, -14, -17, -22, -25, -26, -30, -38, -45, -51, -52, -65, -69, -70, and -71), and 29 are unknown (exhibited no BLASTp hit). The table in the supplemental material presents details of position, size, and nearest match with the results of a BLASTp search.
Structural module.
The ORFs of the VP882 genome are arranged in a typical modular structure. The major structural genes are arranged in the order of portal protein, capsid protein, tail plate assembly genes, tail fiber, tail sheath, tail tube, and tail length protein. ORF6 to ORF34 encode the structural proteins (Fig. 2), including portal protein (ORF6), capsid (ORF9), tail proteins (ORF12, -20, -21, -29, -31, -32, and -34), baseplate and assembly proteins (ORF16, -18, and -19), and tail sheath protein (ORF28). These structural genes are organized in a manner similar to those of lambdoid phages and other plasmid-like phages, with head genes near the left end and tail genes clustered immediately to their right. All of these structural genes are transcribed in a rightward direction (9).
|
View larger version (9K): [in a new window] |
FIG. 2. Schematic genome map of phage VP882. ORFs were numbered arbitrarily starting from the left end. The orientation of each gene is denoted by the direction of the arrow.
|
ORF35 putatively encodes a late-control gene D protein which is similar to late-control proteins of many phages. ORF35 locates at the end of the ORFs for tail proteins, and it also shares amino acid identity with a few tail proteins, such as
HAP-1 (54% identity) (36) and a phage of Xylella fastidiosa (GenBank accession no. NP_779294; 29% identity). No other lytic genes have been identified in VP882.
Genome maintenance module.
ORFs associated with maintenance of the VP882 genome in V. parahaemolyticus are similar to those of other plasmid-like temperate phages. The putative ParA protein of VP882 (ORF53) exhibits high identity with those of marine phage
HAP-1 (60%) (36) and V. harveyi phage VHML (65%) (40). ParA, working with ParB and ParS, is a partition protein associated with the movement of the correct copy number of bacterial chromosome and plasmids into daughter cells (49).
Protelomerase (ORF54) was identified, and it is responsible for the maintenance of the linear plasmid-like phages. The IRS, reaction site of protelomerase, in phage VP882 was searched by the Inverted Repeat Finder program (54), and an inverted palindromic repeat sequence was found between the regions encoding ParA (ORF53) and the protelomerase (ORF54) (Table 2). The IRS of VP882 (region of the genome from positions 26135 to 26247) consists of 113 bp and is located 121 bp upstream from the protelomerase (ORF54).
|
View this table: [in a new window] |
TABLE 2. Inverted repeat sequences in phage VP882 and five other phages with a plasmid-like genome
|
HAP-1 (51% identity) (36), PY54 (30% identity) (24), and N15 (30% identity) (45). These replication proteins exhibit multifunctional activities of primase and helicase, similar to Rep protein of plasmids replicating by the
-mechanism (44).
Lysogeny module.
Putative lysogeny-related genes were identified in phage VP882. Repressor protein (ORF59) will bind to the phage DNA; it physically inhibits the binding of any transcription machinery and may be responsible for the lysogeny of VP882. The antirepressor (ORF60) probably binds to and acts against the main phage repressor (33). Antitermination protein Q (ORF62) is involved in late expression and allows RNA polymerase to override the termination signal and transcribe the lysis genes, which are the capsid protein and processing genes.
The putative conjugative transfer protein of VP882 (ORF61) is similar to TraR of Salmonella sp. (47% identity) (GenBank accession no. NP_490571), Escherichia coli (45% identity) (GenBank accession no. YP_538713) and other species, and it regulates genes that are required for conjugation (10) and may also associate with the transfer of the plasmid-like genome of VP882.
Other ORFs.
Some of the putative functional genes of VP882, namely, the DNA adenine methylase (DAM) (ORF37), exonuclease (ORF46), and the transcriptional regulator (LuxR family; ORF56), are often present in bacterial genomes and also in some phages. These genes may be the remains of a bacterial genome after viral transduction or play specific roles in the phages. DAM has been identified in many bacteria, and the putative DAM of phage VP882 (ORF37) exhibits high identity with similar functional proteins of Halomonas aquarmarina phage
HAP-1 (72%) (36), Methylobacillus flagellatus (67%) (GenBank accession no. YP_546786), Chromobacterium violaceum ATCC 12472 (67%) (GenBank accession no. NP_901781), Burkholderia cepacia AMMD (66%) (GenBank accession no. YP_773740) and other Burkholderia species, and low identity with the DAMS of Vibrio species, such as V. parahaemolyticus (24%) (GenBank accession no. NP_799121), V. vulnificus YJ016 (24%) (GenBank accession no. NP_935779), V. harveyi (24%) (GenBank accession no. ABO14793), V. alginolyticus (23%) (GenBank accession no. ZP_01261704) and others. DAM has multiple roles in cell physiology. DAM modifies the affinity of regulatory proteins toward DNA and is a transcriptional regulator. DAM regulates the expression of pilus-adhesins and other virulence-associated factors in E. coli, Salmonella spp., Yersinia pseudotuberculosis and V. cholerae (31).
The putative transcriptional regulator of VP882 (LuxR family; ORF56) exhibits high identity with those of V. cholerae (43%) (GenBank accession no. NP_233459) and Burkholderia oklahomensis (42%) (GenBank accession no. ZP_02358706) and with the DNA-binding helix-turn-helix domain-containing protein of V. vulnificus (41%) (GenBank accession no. NP_763032). Similar helix-turn-helix or LuxR proteins are present in the VHML phage of V. harveyi (40), VP16C and VP16T and in the VpV262 phages of V. parahaemolyticus (23, 48). LuxR is a quorum sensor and regulates the expression of virulence genes (1, 30).
The DAM and quorum-sensing transcriptional regulators in phage VP882 may regulate physiological or virulence-associated traits of the host, V. parahaemolyticus. VP882 infected and lysed several strains of V. parahaemolyticus, V. vulnificus, and V. cholerae (Table 1); may probably be involved in horizontal genetic shift among these three Vibrio species.
The putative exonuclease of VP882 (ORF46) exhibits a high identity with exonuclease RNase T and DNA polymerase III of Yersinia pseudotuberculosis (58%) (GenBank accession no. ACC88872), Thauera sp. (53%) (GenBank accession no. EDS56474), and Vibrio harveyi (48%) (GenBank accession no. ZP_01987444) and also in some phages, including Klebsiella
KO2 (41%) (GenBank accession no. YP_006619) and enterobacterial phage N15 (39%) (GenBank accession no. NP_046936). These putative exonucleases are associated with the processing of single-stranded DNA or RNA and in the replication of bacterial DNA.
Demonstration of plasmid-like genome.
Based on the presence of typical ORFs of linear plasmid-like phages like those of N15 of E. coli (44) and
HAP-1 of H. aquarmarina (36), VP882 most likely belongs to the same group of viruses. After infection of its host, the phage molecule will become circularized via its cohesive ends, and it is then digested by protelomerase to form a linear molecule with covalently closed ends (Fig. 3) (44).
![]() View larger version (28K): [in a new window] |
FIG. 3. Demonstration of the plasmid-like genome in phage VP882. The locations of the primers used in demonstrating the plasmid-like genomes are indicated, and details of the experimental procedures are addressed in Materials and Methods. (Top) After infection of the host, the virion genome of a linear plasmid-like phage will become circularized via its cohesive ends. (Middle) The protelomerase will then cut at the hairpin region and form covalently closed ends (36, 44). Pt, protelomerase; Pa, ParA. (Bottom) In addition, the host chromosomal DNA preparation (lanes 1 and 2) and the VP882 viral DNA (lanes 3 and 4) were separately amplified with the 37698R and 236L primers, and the amplicons were separated with (lanes 2 and 4) or without (lanes 1 and 3) prior digestion with BfucI or SalI. Lane M contains molecular size markers (Fermentas 1-kb DNA ladder; BioLabTech, Kyiv, Ukraine).
|
For further confirmation, the bacterial chromosomal DNA preparation and the purified VP882 phage DNA were separately amplified by PCR using primers 37698R and 236L (5'-CAGGGTTTTCTTTGTGACAGTGATACCCAT, located at the region of the genome from positions 236 to 208). The amplified fragments were separated directly or subjected to digestion with BfuCI or SalI before separation by agarose gel electrophoresis. Both the host chromosomal DNA and phage DNA preparations yielded a single band of 735 bp (Fig. 3, bottom). This result showed the presence of a circular DNA (Fig. 3, top) or a linearized molecule (Fig. 3, middle) in these preparations.
Comparative analysis of protelomerases and their reaction sites.
Protelomerase is in a dimer form, and the core of its monomer consists of an N-terminal domain, a catalytic domain, and a C-terminal DNA binding domain named the stirrup which is required for the resolution of hairpin telomeres (50). In this study, the amino acid sequences of the protelomerases in VP882 and five other plasmid-like prophages, namely,
HAP-1, VHML,
KO2, N15, and PY54, were analyzed.
The putative protelomerases of phages VP882,
HAP-1, VHML,
KO2, N15, and PY54 consist of 538, 520, 509, 640, 631, and 628 amino acid residues, respectively. The entire amino acid sequence of the putative protelomerase of VP882 was aligned with the protelomerases of five other prophages. Alignment analysis revealed that the protelomerase of VP882 has 54, 26, 33, 33, and 32% identity with that of
HAP-1, VHML, N15,
KO2, and PY54, respectively. The sequence conservation, quality, and consensus of these six aligned protelomerases revealed that they could be divided into two groups, except for VHML prophage. The VP882 and
HAP-1 phages were highly similar, while the other three were in another group (data not shown).
The active site residues of the protelomerases have been analyzed and showed mixed features of the tyrosine recombinases and type IB topoisomerases (50). Our results showed that the Arg223, Lys248, Lys325, Arg328, His361, and Tyr370 residues of the catalytic active site of VP882 protelomerase are characteristic residues also identified in five other plasmid-like prophages (Fig. 4). These results suggested that these plasmid-like prophages probably share a common mechanism in using protelomerase to resolve the telomere in the lysogenic cycle.
|
View larger version (31K): [in a new window] |
FIG. 4. Conserved active site of the protelomerase catalytic domain from six plasmid-like prophages. The protelomerase sequences were aligned and presented by Clustalw2 (29). The asterisk below the sequences means that the residues in that column are identical in all sequences in the alignment. Residues with 70% or more sequence identity are shown on a gray background. The numbers above the sequences indicate the corresponding amino acid residue number in the VP882 protelomerase. Gaps introduced to maximize alignment are indicated by dashes.
|
KO2, and PY54 (2, 18). In this study, we examined the IRSs in VP882 and five other plasmid-like prophages. Our results revealed that VP882 phage possesses a diverse IRS pattern, which is more similar to that of
HAP-1 than the traditional IRS pattern of N15 phage (2, 36). The hairpin structure and catalytic binding domain of the IRS region of VP882 are identical to those of
HAP-1 (Fig. 5B and Table 2), whereas the IRS pattern of the catalytic active site of VP882 is only similar to that of N15 phage (Fig. 5B). Phage VHML also exhibits a diverse pattern compared with five other prophages (Fig. 5B).
![]() View larger version (45K): [in a new window] |
FIG. 5. Diverse IRS patterns found in the VP882 phage genome compared with those of five other plasmid-like phages. (A) Illustration of IRS patterns by using the Weblogo tool (16). The predicted cutting site (center hairpin region) of the protelomerase is indicated by a short black vertical line below the graphs. The nucleotide sequences shown inside the boxes are the consensus binding sites of protelomerase genes. (B) Using the binding site of the protelomerase from N15 phage as a reference template (2, 18) to compare the consensus binding patterns among the six plasmid-like phages. The solid arrow lines and the broken arrow lines indicate that the regions have the same sequence or exhibit a different sequence from that of the reference template, respectively. The hairpin region (H) and the binding regions for the stirrup domain (S) and the catalytic domain (C) are indicated.
|
KO2, and PY54, which have three conserved patterns of hairpin structure and catalytic [(c/g)CAT(t/a)(a/c)TACGC] and stirrup binding domains [CAC(c/a)(t/c)A(t/a)(t/c)]. The other group consisted of phage VP882 and
HAP-1, which have only two conserved sequences in hairpin structure and catalytic binding domains, CCCATACTATAC (Fig. 5A and Table 2).
Comparative genomics of six plasmid-like prophages.
Following the addition of VP882, a total of 18 completed genome sequences of vibriophages are available in the NCBI genome database; these are 10 phages from V. cholerae (fs1, VSKK, VSK, KSF-1
, VGJ
, fs2, K139, VP4, VP5, and VP2), seven phages from V. parahaemolyticus (VfO4K68, Vf12, Vf33, VfO3K6, VP882, VpV262, and KVP40), and one phage (VHML) from V. harveyi. While 10 of them are filamentous phages (Inoviridae), 4 are Myoviridae, and another 4 are Podoviridae.
The genome of phage VP882 has certain specific properties; for example, its G+C percentage (56.95) is among the highest of these vibriophages with complete genome sequences. The average G+C percentage of these vibriophages is 46.61, which is close to the values of their hosts. The G+C percentages of V. cholerae and V. parahaemolyticus are 47.48 (GenBank accession no. AE003852.1 and AE003853.1, respectively) and 46% (GenBank accession no. BA000032.2 and BA000031.2, respectively), respectively. Phages VP16T (GenBank accession no. AY328852) and VP16C (GenBank accession no. AY328853) with incomplete genome sequence both have G+C content of 59% (48). A difference of about 11% exists between the VP882 genome and the genome DNA of V. parahaemolyticus. The sequence variations in these vibriophages were further analyzed using the DOTTER program, which is a pair-wise sequence comparison tool (50). Results revealed that the vibriophage VP882 does not exhibited relatedness with other vibriophages (data not shown).
Since there is little sequence similarity of VP882 genome with other complete sequenced vibriophages, we investigated whether VP882 has sequence similarity with five selected plasmid-like prophage genomes by the DOTTER program. The results indicated that the genomic sequence VP882 phage is more similar to the
HAP-1 genome, while the N15,
KO2, and PY54 genomes are more similar to each other than to VP882 (Fig. 6). Moreover,
HAP-1 was the only BLASTN match when the whole VP882 genome was used to search against the NCBI genomic refseq database. From the genomic sequence and IRS and protelomerase sequence comparative analysis, both VP882 and
HAP-1 are rare plasmid-like marine prophages, unlike N15 and other typical plasmid-like prophages, and may be classified into a new subtype of the plasmid-like phages.
![]() View larger version (23K): [in a new window] |
FIG. 6. Genomic alignment of the VP882 and five other plasmid-like phage genomes. The sequence similarity was visualized with the DOTTER program (50) by using a sliding window of 25 bp. The arrangement of the phages was ordered according to their genome sizes, from smallest to largest.
|
Published ahead of print on 13 March 2009. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
and the vibrio pathogenicity island prophage in Vibrio mimicus: evidence for recent horizontal transfer between V. mimicus and V. cholerae. Infect. Immun. 68:1507-1513.
HAP-1 of Halomonas aquamarina possesses a linear plasmid-like prophage genome. J. Virol. 82:6618-6630.This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»