Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genomic and Proteomic Analyses of the Terminally Redundant Genome of the Pseudomonas aeruginosa Phage PaP1: Establishment of Genus PaP1-Like Phages

  • Shuguang Lu,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Shuai Le,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Yinling Tan,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Junmin Zhu,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Ming Li,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Xiancai Rao,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Lingyun Zou,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Shu Li,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Jing Wang,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Xiaolin Jin,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Guangtao Huang,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Lin Zhang,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Xia Zhao,

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

  • Fuquan Hu

    hufuquan2009@yahoo.cn

    Affiliation Department of Microbiology, College of Basic Medical Science, Third Military Medical University, Chongqing, China

Abstract

We isolated and characterized a new Pseudomonas aeruginosa myovirus named PaP1. The morphology of this phage was visualized by electron microscopy and its genome sequence and ends were determined. Finally, genomic and proteomic analyses were performed. PaP1 has an icosahedral head with an apex diameter of 68–70 nm and a contractile tail with a length of 138–140 nm. The PaP1 genome is a linear dsDNA molecule containing 91,715 base pairs (bp) with a G+C content of 49.36% and 12 tRNA genes. A strategy to identify the genome ends of PaP1 was designed. The genome has a 1190 bp terminal redundancy. PaP1 has 157 open reading frames (ORFs). Of these, 143 proteins are homologs of known proteins, but only 38 could be functionally identified. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and high-performance liquid chromatography-mass spectrometry allowed identification of 12 ORFs as structural protein coding genes within the PaP1 genome. Comparative genomic analysis indicated that the Pseudomonas aeruginosa phage PaP1, JG004, PAK_P1 and vB_PaeM_C2-10_Ab1 share great similarity. Besides their similar biological characteristics, the phages contain 123 core genes and have very close phylogenetic relationships, which distinguish them from other known phage genera. We therefore propose that these four phages be classified as PaP1-like phages, a new phage genus of Myoviridae that infects Pseudomonas aeruginosa.

Introduction

Bacteriophages (phages) are ubiquitous in the biosphere [1]. Estimations of phage numbers, ranging from 1030 to 1032 in total, are approximately tenfold higher than those of bacteria [2]. Numerous phage investigations have been performed worldwide since Frederick William Twort and Felix dHerelle first reported the discovery of phages in 1915 and 1917, respectively [3], [4]. Phages are potential antimicrobial agents in various clinical or agricultural settings [5], [6] and have become important molecular and biological tools in facilitating the development of bioscience. Approximately 6300 different phages have been examined by electron microscopy [7]; however, only 759 of these (721 infecting bacteria and 38 infecting archaea) have been completely sequenced based on data from the National Center for Biotechnology Information (NCBI; http://www.ncbi.nlm.nih.gov/; Bethesda, MA, USA, Oct. 2012). This number is far lower than the number of sequenced bacteria (3433 complete genomes of bacteria and 199 complete genomes of archaea as of 28 Oct. 2012). A detailed dissection of phage genomes would add valuable data to our knowledge of phages and help us understand the evolutionary relationships between phages and bacteria.

Pseudomonas aeruginosa (P. aeruginosa), an opportunistic pathogen, is ubiquitous in the environment and often resistant to a large number of antibiotics. As such, the treatment of P. aeruginosa infections is very difficult [8]. Investigating the biology of phages is important for humans to fight multiple-drug resistant pathogens [9]. Sixty-three complete genome sequences of Pseudomonas phages, most of which infect P. aeruginosa, have become available in GenBank as of 28 Oct. 2012. Among these genome sequences, 18 P. aeruginosa phages belong to the Myoviridae family. Members of this family are efficient killers of bacteria and can affect many aspects of bacterial ecology and evolution. P. aeruginosa phages have been studied for decades for use as therapeutics and typing agents [10]. This group of phages seems to be taxonomically diverse and genetically dissimilar [11]. Currently, most characterized myoviruses of P. aeruginosa are classified into four genera, namely, phiKZ-like phages, P2-like phages, PB1-like phages, and KPP10-like phages [12][15]. Some P. aeruginosa phages have been characterized but remain unclassified. Detailed characterizations of novel P. aeruginosa phages will be significant for understanding the interactions between P. aeruginosa and its phages and the exploration of useful therapeutic reagents against P. aeruginosa.

The interests of our group are focused on the dissection of phage genomes and their biological issues. We previously isolated and identified three P. aeruginosa phages from the sewages of our affiliated hospitals and designated them as PaP1, PaP2, and PaP3. PaP1 is a virulent phage, whereas both PaP2 and PaP3 are temperate phages. The genome of PaP3 can integrate into the chromosome of the host bacteria through a tRNA gene locus [16].

The present work focuses on the genomic and proteomic analyses of the phage PaP1. The results indicate that PaP1, in taxonomy, belongs to Myoviridae, and its genome has terminally redundant ends of 1190 bp. High-performance liquid chromatography-mass spectrometry (HPLC-MS) identified 12 PaP1 structural protein coding genes. Based on comparative genomic analysis, we propose that P. aeruginosa phages PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 be classified as a new genus named “PaP1-like phages” within the Myoviridae family.

Results

Morphology of PaP1

PaP1 is in structure and dimensions identical to P. aeruginosa phage PB1 [14]. The head is an icosahedron, as evidenced by the simultaneous presence of hexagonal and pentagonal capsids, and measures 68–70 nm between opposite apices (Figure 1). Shallow depressions in uranyl acetate indicate the presence of capsomers. The head is separated from the tail sheath by an 8 nm-long neck. Uncontracted tails measure 138–140 × 17–20 nm. Contracted tails measure 55 × 22 nm. Base plates are poorly visible on extended tails. Upon contraction, they separate from the sheath and appear as disks of 23 × 3 nm. There are at least 4 straight tail fibers (Figure 1B). In the quiescent tail, the fibers are folded along the sheath.

thumbnail
Figure 1. Electron micrographs of purified PaP1 phage particles.

A shows two PaP1 particles with uncontracted tails and an empty head (uranyl acetate). B shows the contracted tail with straight tail fibers (phosphotungstate). C shows a pentagonal head (uranyl acetate). The scale bar represents 100 nm. These micrographs were taken by Hans-Wolfgang Ackermann, School of Medicine, Laval University, Quebec, Canada.

https://doi.org/10.1371/journal.pone.0062933.g001

Biological characteristics of PaP1

PaP1 forms clear plaques (∼3 mm in diameter) surrounded by a small semitransparent halo on the lawns of the host bacteria. In rich liquid medium, PaP1 is amplified to high titers (∼1011 PFU/ml). PaP1 particles are stable for over two months of storage at 4°C and resistant to chloroform. Based on the one-step growth curve of PaP1 (Figure 2), its latent period is about 20 min, its burst period is about 40 min, and the average number of PaP1 progeny produced from one host bacterium is about 65. PaP1 can lyse two other strains of P. aeruginosa (PA4 and PA6), aside from PA1. This result may provide a potential basis for the application of PaP1 in phage therapy.

thumbnail
Figure 2. One-step growth curve of phage PaP1.

Experiments were repeated three times with duplicate samples. The Y-axis shows the log of plaque forming units per milliliter (PFU/mL).

https://doi.org/10.1371/journal.pone.0062933.g002

General genomic characteristics of PaP1

The genome of PaP1 consists of 91,715 base pairs (bp) with a G+C content of 49.36%, which is significantly less than that of its host (66.34%). The GC skew of the PaP1 genome is shown in Figure 3. In viral genomes, the lowest point on the GC skew curve is typically the origin of replication [17]. Therefore, the putative replication origin of the PaP1 genome is at the end of the genome sequence (Figure 3).

thumbnail
Figure 3. Circular conformation of GC skew and GC plot of the PaP1 genome.

The GC skew is calculated as (G-C)/(G+C) and the GC plot shows GC% content plot. The outer circle represents the GC skew (red for positive and blue for negative); the inner circle represents the GC plot (pink for above-average and powder blue for below-average). The red pane indicates the putative replication origin of the PaP1 genome.

https://doi.org/10.1371/journal.pone.0062933.g003

No direct repeats of more than 50 bp, inverted repeats of more than 26 bp, or mirror repeats of more than 17 bp were found in the PaP1 genome, indicating that it does not contain complicated secondary structures. Twelve tRNA genes were found in the PaP1 genome (Table 1). Among the 12 encoded tRNAs, tRNAArg, tRNAGln, and tRNAGly are used preferentially by PaP1, but not by the host. Codon usage analysis indicated that the three tRNAs are important for the protein expression of PaP1, since phage tRNA genes can overcome differences in codon usage between the phage and the host [18].

Determination of PaP1 genome termini

The PaP1 genome was assembled as a circular molecule when sequencing was completed. The restriction endonucleases NarI and NotI, both of which have only one cut site in the genomic DNA, were selected to digest the DNA and released two short bands (about 2.5 and 6.5 kb) in the gel (Figure 4A) This result indicates that the PaP1 genome is a linear molecule. Figure 4B shows the recognition sites of NarI and NotI within the linear PaP1 genome. The DNA band indicated by the red arrow in Figure 4A is the 3′ end fragment of the PaP1 genome. The restriction endonuclease FspI was used to digest the PaP1 genome and release the 5′ end fragment, as indicated by the red arrow in Figure 4C. Both the 3′ and 5′ end fragments were purified and used as templates for terminal run-off sequencing with primers P1 and P2 (Figure 4D), respectively.

thumbnail
Figure 4. Identification of PaP1 genome ends.

(A) Digestion of the PaP1 genome DNA by NarI and NotI. The red arrow indicates the 3′ end fragment of the PaP1 genome. (B) NarI and NotI restriction map of the PaP1 genome. (C) PaP1 DNA was digested by FspI and the recovered 5′ end fragment (5′ Frag., indicated by the red arrow) was digested by S1 nuclease. (D) Strategy designed to identify PaP1 genome ends. Primers P1 and P2 are annealed to 3′ and 5′ end fragments of the PaP1 genome DNA, respectively. Terminal run-off sequencing of the two ends (also shown in Figure 5C) is then performed. Case 1 (blunt end): The two sequences obtained by P1 and P2 do not have repeated regions and they can be assembled to the PaP1 genome sequence with no gap between them. Case 2 (3′-protruded end): The two sequences obtained by P1 and P2 also do not have repeated regions; however, a gap is observed between the sequences once assembled to the PaP1 genome sequence. The sequence within the gap is the 3′-protruded cohesive sequence. Case 3 (5′-protruded end or terminal redundancy): A repeat between the two obtained sequences is observed. If the repeated sequence is less than 100 bp, it is regarded as the 5′-protruded cohesive sequence [16], [19], [20]; however, if the repeated sequence is over 100 bp, it is regarded as a terminally redundant sequence [12], [62], [63].

https://doi.org/10.1371/journal.pone.0062933.g004

The results of terminal run-off sequencing coincided with case 3 in Figure 4D. The two sequences obtained by P1 and P2 have a long repeat of 1190 bp, suggesting that the repeated sequence is terminally redundant rather than cohesive, since the latter is usually less than 100 bp [16], [19], [20]. We used S1 nuclease, a single-strand digesting enzyme, to digest the 5′ end fragment. The cut and uncut 5′ end fragments had the same size (Figure 4C), indicating that the 1190 bp terminal sequence is double-stranded. Thus, the natural structure of the PaP1 genome DNA molecule can be described as shown in Figure 5A. Two functionally unknown genes (g156 and g157) were found in the terminally redundant region (Figure 5B).

thumbnail
Figure 5. Structure of PaP1 terminal redundancy.

(A) PaP1 genome DNA with terminal redundancy at both ends. The sequence data of the PaP1 genome are calculated from 1 bp to 91,715 bp. (B) Putative genes within the PaP1 terminal redundancy. (C) Terminal run-off sequencing chromatograms of both ends of the PaP1 genome. The sequencing direction of primers P1 and P2 is indicated by the black arrows.

https://doi.org/10.1371/journal.pone.0062933.g005

Identification and organization of PaP1 genes

The open reading frames (ORFs) of the PaP1 genome were identified using ORF Finder [21] with ATG, GTG, and TTG as start codons. A total of 541 ORFs (>100 bp) were predicted from the PaP1 genome, among which 157 ORFs were identified as protein coding genes. The average length of a gene is about 504 bp. Only about 13.8% of the PaP1 genome belongs to non-coding regions (Table 1), which is in concordance with the fact that the ORFs of tailed dsDNA phages are tightly and efficiently organized, with little space between genes. The space between genes is usually occupied by putative regulatory sequences, such as promoters and terminators. Putative promoters and terminators of the PaP1 genome are listed in Table 2 and 3, respectively.

Sequence homology analyses revealed putative functions for some of the 157 coded proteins. A total of 143 PaP1 proteins have homologs in other species, only 38 of which are functionally identified. Both the genome sequence of PaP1 and its annotation have been deposited in the GenBank database, under accession number HQ832595. The detailed annotation and organization of the PaP1 genome is listed in Table S1 and illustrated in Figure 6. The data shows that the PaP1 genome may be divided into several functional modules, revealing an apparent mosaic structure, which is one of the striking characteristics of phage genomes. This finding also suggests that tailed phage genomes evolve from combinations of modules from different species [22].

thumbnail
Figure 6. Diagram of the PaP1 genome with functional annotations.

Blue arrowheads show genes on the plus strand; brownish red arrowheads show genes on the minus strand; green words refer to the amino acid transferred by the corresponding tRNA. NPR: nicotinamide phosphoribosyl, PRP: phosphoribosylpyrophosphate, RDR: ribonucleotide diphosphate reductase.

https://doi.org/10.1371/journal.pone.0062933.g006

Two functionally unknown modules reside near the 5′ and 3′ ends of the PaP1 genome, respectively. A large number of small genes with unknown functions cluster in the two modules. The 3′ end module consisting of g114 to g157 may play critical roles in early transcriptional events. PaP1 codes for its own RNA and DNA polymerases, which are probably involved in the synthesis of its RNA and DNA molecules. At least 10 genes cluster on the minus strand (Figure 6), which probably controls the nucleotide metabolism system of PaP1. These genes can convert the metabolism of the host cell to produce progeny phages [23]. The products of g032 and g110 are dCMP deaminase and thymidylate synthase, respectively, both of which are involved in dTTP synthesis [24]. Phage-encoded thymidylate synthase appears to have diverged from the precursor of the host bacteria. No sequence homologs to phage integrases, repressors, transposases, or excisionases were found, supporting our conclusion that PaP1 is a lytic phage.

Genes closely related to phage morphogenesis also cluster together, among which 12 genes encoding structural proteins were identified (as below). The products of g029 and g072 are cell wall hydrolase and endolysin, respectively, both of which belong to enzybiotic factors [25], [26]. Although endolysin and holin usually constitute a two-component lysis system for the liberation of phage progenies from the host cell, no holin homolog was found in the PaP1 genome. The endolysin gene (g072) of PaP1 has been expressed in our laboratory, and we have shown that the purified product can hydrolyze the cell wall peptidoglycan of P. aeruginosa [27].

Identification of phage PaP1 structural proteins

A dsDNA phage particle is made up of a series of structural proteins and a single DNA molecule containing the entire genome. The PaP1 structural protein coding genes cluster in a module of its genome and are preceded by a terminase gene (Figure 6). The terminase plays an important role in DNA packaging. To identify the structural proteins of PaP1, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) was used to visualize each structural protein in the gel (Figure 7). At least 17 proteins with molecular weights ranging from 6 kDa to 80 kDa were resolved. Each protein band was then excised for HPLC-MS, permitting the allocation of 15 protein bands to 12 corresponding PaP1 genes (Figure 7). The detailed parameters and results of the mass spectrometry are shown in Table 4. The sequence coverage reaches up to 58%. The sequence coverages of gp067 and gp071 are 3% and 4%, respectively, which are relatively low compared with other proteins; hence, the identification of these proteins as structural components of the phage must be confirmed further.

thumbnail
Figure 7. SDS-PAGE analysis of the structural proteins of phage PaP1.

Proteins were visualized in a 15% (w/v) gel and identified by HPLC-MS analysis (Table 4).

https://doi.org/10.1371/journal.pone.0062933.g007

thumbnail
Table 4. Mass spectrometry identification of PaP1 structural proteins.

https://doi.org/10.1371/journal.pone.0062933.t004

The predominant band is, as predicted, the major capsid protein (gp051, ∼40 kDa); the band (∼43 kDa) just above it was also identified as a major capsid protein by mass spectrometry (Figure 7). Peptides corresponding to gp050 (∼15 kDa) were found in three bands (Figure 7), and similar to gp051, a small band (∼16 kDa) just above it was also identified as gp050. This observation may be explained as the result of the known carry-over effect [28], [29] of the massively overrepresented major capsid protein and gp051 bands. The gp050 band at the bottom of the gel (∼6 kDa) suggests posttranslational processing [30] of gp050. An unusual finding was that gp104 (Figure 7), a p09 homolog of the phage PaP3 (Table S1), is not located in the late gene cluster for phage morphogenesis but among the genes probably involved in DNA replication and control (Figure 6). The four identified PaP1 structural proteins, tape measure, tail fiber, baseplate, and major capsid, with molecular weight ranging from 40 kDa to 80 kDa, are important for phage PaP1 particle formation. Interestingly, the major capsid protein of PaP1 shares molecular weights and amino acid sequences identical to those of the P. aeruginosa phages JG004 [31], PAK_P1 [32] and vB_PaeM_C2-10_Ab1 (Table 5), indicating a close relationship among these four phages.

thumbnail
Table 5. Comparison of phages with BlastN scores of more than 200 against the PaP1 genome.

https://doi.org/10.1371/journal.pone.0062933.t005

Comparative genomic analysis of PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1

Comparative analysis of genomes.

Many PaP1 proteins share great similarity with the homologs of the P. aeruginosa phages PAK_P1, JG004, and vB_PaeM_C2-10_Ab1, among which 15 PaP1 proteins share 100% sequence identity with PAK_P1 homologs (Table S1). When the PaP1 genome sequence was compared with the nucleotide database, 11 P. aeruginosa phages were found with BlastN scores of more than 200 (Table 5). In particular, PAK_P1, JG004, and vB_PaeM_C2-10_Ab1 showed query coverages of more than 90% (Table 5).

A graphical comparison (Figure 8) was performed to illustrate the genomic similarities of PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1. This figure is a visualized description of the corresponding data listed in Table 5. The whole genome sequences of phage PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 show great similarities and most of their DNA sequences appear to have descended from a single common ancestral phage. We also performed a dot plot comparison of the genome sequences of PaP1, JG004, PAK_P1, vB_PaeM_C2-10_Ab1, PAK_P3, and KPP10 (Figure 9). The results are in concordance with Table 5 and Figure 8.

thumbnail
Figure 8. Pairwise nucleotide sequence comparison of phages closely related to PaP1.

Comparisons were conducted using BLAST 2.25 and displayed using ACT [56]. Highly related sequences are shown by the blue shadings. The intensity of the blue coloration indicates the level of sequence similarity. The minimum score cutoff is 100 and the minimum identity cutoff is 50%. #The full name is vB_PaeM_C2-10_Ab1.

https://doi.org/10.1371/journal.pone.0062933.g008

thumbnail
Figure 9. Dot plot of genome sequences of six phages using the program Gepard.

The word length used is 9 bp; other parameters are set by default. The black dots indicate that the corresponding genome regions of the abscissa and the ordinate show similarity to each other. #The full name is vB_PaeM_C2-10_Ab1.

https://doi.org/10.1371/journal.pone.0062933.g009

At the protein level, PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 show striking similarities. Intriguingly, the major capsid proteins of these four phages share 100% sequence identity with each other, which explains their similar morphologies [31,32, and the present work]. As shown in Figure 8, the regions with no blue shading represent minor insertions or deletions among the genome sequences of these four phages. These DNA regions are probably the accumulated mutations for phage adaptation to the host bacteria. The main differences between these four phage genomes are located in their tail fiber encoding genes, indicating that these phages may have evolved different host cell adsorption mechanisms.

Protein homology analysis.

The genomic comparison indicated that phages PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 are closely related. To investigate this further, protein homology analysis was performed. The result reveals that the PaP1 genome shares 123 (78.34%) homologs with JG004, PAK_P1, and vB_PaeM_C2-10_Ab1, and shares 55 (35.03%) homologs with KPP10 and PAK_P3, and shares even fewer homologs with other phage groups. These results strongly indicate that PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 are closely associated and distinguishable from other phage groups. Therefore, these four phages do not belong to any known phage genus and may have descended from a common ancestor.

Phylogenetic analysis.

Phylogenetic analysis was performed based on the major capsid proteins. Since related phages are considered to have similar head structural components, they may be clustered based on their major capsid proteins [33]. We chose phages that are listed in Tables 5 and 6 to analyze the phylogenetic relationships between them, and a phylogenetic tree was constructed (Figure 10). Phage PaP2 was included in the tree because it was identified from the same sample (hospital sewage) from which PaP1 and PaP3 were obtained.

thumbnail
Figure 10. Phylogenetic analysis of major capsid protein amino acid sequences.

The diagram was constructed using the MEGA5 program [58]. The relative distances of each main branch are shown in the figure. Both PaP2 and LUZ24-like phages belong to the family of Podoviridae; all other phage groups belong to the Myoviridae phage family. The phages in the same group are marked with the same color. The group of “PaP1-like phages” is first presented in this work. #The full name is vB_PaeM_C2-10_Ab1.

https://doi.org/10.1371/journal.pone.0062933.g010

thumbnail
Table 6. Overview of genomic characteristics of Myoviridae that infect P. aeruginosa.

https://doi.org/10.1371/journal.pone.0062933.t006

Figure 10 shows that different P. aeruginosa phage genera cluster in the phylogenetic tree based on the major capsid proteins. This observation is in accordance with the data listed in Table 6, in which myoviruses of P. aeruginosa have been assigned to several phage genera, except for PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1. As expected, these four phages are closely clustered, distinguishing them from other P. aeruginosa phages. This finding reinforces the idea that these four phages descend from a common ancestor. Thus, PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 should be grouped as a new phage genus: “PaP1-like phages” (Figure 10).

Discussion

In the present work, the newly isolated phage PaP1 was assigned as a member of the Myoviridae family. Over 96% of the investigated phages belong to the tailed phages. A total of 6054 tailed phages are known, among which 1558 are from the family Myoviridae [7]. Although members of Myoviridae have many common features, they actually represent a diverse collection of phages.

Phage PaP1 has a linear genome consisting of 91,715 bp with a terminal redundancy of 1190 bp. We designed a new strategy to determine the genome ends of PaP1 (Figure 4D) and this strategy is useful for the identification of genome ends of many other phages. Many known phages have terminally redundant genomes, such as phages T3, T7, P22, SPP1, and T4. Similar to PaP1, several Myoviridae phages of P. aeruginosa (e.g., phiKZ, EL, JG004, and KPP10) also have terminal redundancies (Table 6). These phages employ a variety of mechanisms to generate long DNA concatemers with terminal redundancy, which ensures phage replication without any loss of genetic information [34]. Intriguingly, the terminally redundant region of phage SPO1 contains a “host take-over module” composed of a cluster of 24 genes. This region is responsible for shutting off transcription and translation of the host genes [35]. The terminally redundant region of the PaP1 genome contains two genes with unknown functions (Figure 5B). Identification of the two genes, the PaP1 terminase, and the terminase recognition (pac) site [36], may provide a basis for a understanding of phage PaP1 morphogenesis.

The 91,715 bp PaP1 genome encodes 157 putative proteins. Based on the predicted functions of these proteins, the genome can be divided into several functional modules, showing an apparent mosaic structure that is characteristic of the phage genomes [22], [37]. The PaP1 genome also contains 12,660 bp worth of non-coding regions. Non-coding cannot be interpreted as an indicator of no biological function because some ncRNAs (non-coding RNAs) or DNA binding motifs may be observed within the non-coding regions in phage genomes, which may be useful for the phages but toxic to their hosts [38]. Less than a quarter of the 157 putative proteins have homologs with known functions. Twelve ORFs of the PaP1 genome were identified as structural protein coding genes (Figure 7). The majority of the phage ORFs have unknown functions, which hinders phage studies. The sequence data of these 12 structural proteins have been added to the phage proteomic pool. As similar data emerge, the collective information will be valuable for future phage studies.

Comparative genome analysis revealed that the PaP1 genome shows great similarity with JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 both at the DNA and protein levels, distinguishing them from other Myoviridae. In essence, the phages have very few relationships with other phage genera. Some similarities consistently exist among the tailed phages, suggesting that phages may undergo genetic material exchange from a large shared pool [22]. Substantial evidence suggests that tailed phages may be of very ancient origin and it has been proposed that all of the dsDNA tailed phages share common ancestry [39]. The comparison of coat protein structure and virion architecture can provide a sound basis for grouping viruses together [33]. The major capsid protein of PaP1 shares 100% identity with those of JG004, PAK_P1, and vB_PaeM_C2-10_Ab1; their particles also share identical morphology. This observation supports the idea that PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 can be grouped together, as shown in Figure 10.

According to similarities in biological characteristics and both DNA and protein sequences, prominent core genes, and the close phylogenetic relationships among PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1, we propose that these four phages, with PaP1 as the type virus, can be grouped as a new genus (PaP1-like phages) of myovirus bacteriophages, as shown in Figure 10. We predict that other newly characterized phages will be assigned to “PaP1-like genus” in the near future, thereby contributing to our understanding of phage biology.

Materials and Methods

Pseudomonas aeruginosa strains

Six P. aeruginosa strains (PA1, PA2, PA3, PA4, PA5, and PA6) were isolated at the second affiliated hospital of the Third Military Medical University, Chongqing, China, and cultivated in our laboratory. These strains belong to serogroups 9, 20, 6, 6, 20, and 11 of the P. aeruginosa international antigenic typing system, respectively. All six strains were cultivated at 37°C in LB medium with shaking for ∼5 h to reach the log phase.

Phage propagation and purification

Phage PaP1 was isolated from the hospital sewage using P. aeruginosa PA1 as host bacterium, based on a standard lambda phage isolation protocol [40]. A liquid culture of the PA1 strain of the log phase growth was infected with PaP1 (MOI of 1/100) and incubated at 37°C with shaking. After ∼5 h, the culture showed signs of lysis, and a few drops of chloroform were added to it. The culture was then centrifuged at 10,000 g for 5 min, and the supernatant was stored at 4°C for subsequent experiments. After storage at 4°C for over two months, the supernatant was diluted, plated onto a Petri dish overlaid with the PA1 stain, and then cultured at 37°C until individual plaques could be picked to test the titers of PaP1 in the supernatant. Another 5 P. aeruginosa strains (PA2, PA3, PA4, PA5, and PA6) were used as host bacteria to test whether or not PaP1 could lyse them. One-step growth experiments of PaP1 were performed, as previously described [41], to determine phage growth characteristics. Crude phage suspensions of PaP1 were concentrated and purified by PEG8000 precipitation according to the method of Govind et al. [42]. The purified PaP1 particles were further purified using CsCl gradient ultracentrifugation [43].

Transmission electron microscopy (TEM)

Filtered phage lysates (about 1011 PFU/mL) were sedimented for 60 min at 25,000 g in a Beckman J2-21 centrifuge (Palo Alto, CA, USA) equipped with a JA1.1 fixed-angle rotor, followed by washing in neutral ammonium acetate buffer (0.1 M) under the same conditions. Phage particles were deposited on carbon-coated copper grids, stained with uranyl acetate (2%, pH 4.5) or potassium phosphotungstate (2%, pH 7.0), and examined under a Philips EM 300 electron microscope. Magnification was monitored with T4 phage tails. Dimensions of PaP1 particles are calculated from 20 particles.

DNA extraction and sequencing

EDTA to a final concentration of 20 mM, proteinase K at 50 µg mL−1, and sodium dodecyl sulfate at 0.5% (w/v) were added to the purified phage PaP1 stock solution. The mixture was incubated at 56°C for 1 h, after which an equal volume of phenol-chloroform-isoamyl alcohol (25∶24∶1) was added to it, followed by centrifugation at 5000 g for 10 min. The aqueous layer was extracted with chloroform at 5000 g for 10 min. The aqueous layer was collected, mixed with 0.6 volumes of isopropanol, and then stored at −20°C for one night. The mixture was centrifuged at 4°C and 12,000 g for 10 min, and the precipitated DNA was collected and washed with 70% and 100% ethanol, respectively. The obtained PaP1 DNA was suspended in TE buffer (pH 8.0) and stored at −20°C for use. DNA sequencing was carried out at the Chinese National Human Genome Center (Shanghai, China) using the Roche/454 GS FLX Titanium system [44]. Roche/454 sequence reads were assembled using the Phred/Phrap/Consed software package [45].

Analysis of PaP1 genome ends

Simulation of the restriction enzyme mapping of the PaP1 genome sequence was performed using the software package DNAStar [46]. The PaP1 DNA was digested by selected restriction endonucleases (NarI, NotI, and FspI, purchased from New England Biolabs, Ipswich, MA, USA). For a reaction system of 20 µL, 10 units of the restriction endonuclease (NarI or NotI) and 200 ng of PaP1 DNA were used. The mixture was incubated at 37°C for 120 min and then used to perform agarose gel electrophoresis. For a reaction system of 100 µL, 1 µg of PaP1 genome DNA and 50 units of restriction endonuclease (NotI or FspI) were used. The mixture was incubated at 37°C for 100 min. Agarose gel electrophoresis was subsequently performed to separate the restriction fragments containing the 5′ and 3′ ends of the PaP1 genome. The restriction fragments containing the 5′ and 3′ ends of the PaP1 genome were purified using Wizard SV Gel and PCR Clean-up System (Promega, Fitchburg, WI, USA), respectively. Terminal run-off sequencing was carried out by BGI-Shenzhen (Shenzhen, China). The 3′ end fragment was sequenced using primer P1 (5′-CGTTCGACGATCCGATGC-3′), and terminal run-off sequencing of the 3′ end fragment was performed by P1. The 5′ end fragment was sequenced using primer P2, which represents three primers (P2a, P2b, and P2c). P2a (5′-CGCCGATGGTCTAGCTGTTG-3′) was the first primer used to sequence the 5′ end fragment. P2b (5′-TCGCCTTCTGCCAGTTATG-3′) was designed based on the DNA sequence acquired by P2a. P2c (5′-ATGCCTTGTCGCAGTTGG-3′) was designed based on the DNA sequence acquired by P2b, and terminal run-off sequencing of the 5′ end fragment was performed by P2c. These primers were prepared by BGI-Shenzhen (Shenzhen, China). We used a strategy to explore the terminal sequence of the PaP1 DNA (Figure 4D). Digestion of the 5′ end fragment with S1 nuclease (Takara Bio, Shiga, Japan) at 23°C for 20 min was carried out to further identify the terminally redundant genome of phage PaP1.

Sequence analysis and genome annotation

The software packages DNAStar [46] and DNAMAN (http://www.lynnon.com/) were used to analyze the basic features of the PaP1 genome sequence. The GC skew of the PaP1 genome was analyzed using DNAPlotter [47]. The internet tool tRNAscan-SE 1.21 [48] was used to predict tRNA genes in the DNA sequence with a cove score cutoff of 20. ORFs were analyzed using NCBI ORF Finder [21], and phage genes were predicted using the software GeneMark.HMM [49] with a length threshold of 100 bp. DNA sequences and protein sequences were scanned for homologs using BLAST [50]. Predicted promoter regions were identified using neural network promoter prediction [51], and putative terminator structures were identified using the web tool FindTerm (http://linux1.softberry.com/berry.phtml).

SDS-PAGE and HPLC-MS of the PaP1 structural proteins

The purified phage particles were resuspended in SDS-PAGE loading buffer [52] and boiled for 5 min before loading onto a 15% (w/v) polyacrylamide gel to identify structural proteins of phage PaP1. We also performed 12% (w/v) and 10% (w/v) SDS-PAGE to better separate proteins of different molecular weight ranges. Protein bands were visualized by staining with Coomassie Brilliant Blue R250 dye for 1 h with shaking and washing with methanol-acetic acid-H2O (5∶1∶4). Slices were excised from the gel and digested as described previously [53]. HPLC-MS was performed using an HPLC-CHIP-MS/MS ION TRAP 6330 system (Agilent, Santa Clara, CA, USA). The acquired data (mass signals) were compared with all of the putative protein sequences of PaP1 using Mill proteomics software (Rev A.03.02.060; Agilent, Santa Clara, CA, USA) to determine genes with products corresponding to the selected protein bands.

Comparative genome analysis

Sequences of the P. aeruginosa phages JG004 [31], vB_PaeM_C2-10_Ab1 (GenBank accession No.: HE983845), PAK_P1 [32], KPP10 [15], and PAK_P3 [54] were compared with PaP1. Dot plot comparison of the phage genomes was performed and displayed using Gepard [55]. Paired genome comparison was conducted using BLAST [50] and the results were displayed using ACT [56]. Multiple sequence alignments of major head proteins (also called major capsid proteins) were conducted using ClustalW [57] with default parameters. A phylogenetic tree was also constructed and displayed using MEGA5 [58] via the neighbor-joining method [59]. Protein homology analysis of PaP1, JG004, PAK_P1, and vB_PaeM_C2-10_Ab1 was performed by CoreGenes [60], [61] with a BlastP threshold score of 100.

Supporting Information

Table S1.

Predicted genes and proteins of phage PaP1.

https://doi.org/10.1371/journal.pone.0062933.s001

(DOCX)

Acknowledgments

We wish to thank Prof. Hans-Wolfgang Ackermann for his professional help in all aspects of our electron microscopy studies, Wenqi Huang for her assistance in our electron microscopy studies, and Ying Wang for her expert help in HPLC-MS. We are also grateful to Qingshan Ni for his help in the bioinformatic analysis and Dr. Lixin Zheng for reading and checking our manuscript.

Author Contributions

Conceived and designed the experiments: FH S. Lu YT. Performed the experiments: S. Lu S. Le. Analyzed the data: S. Lu ML XR L. Zou S. Li. Contributed reagents/materials/analysis tools: JZ JW XJ GH L. Zhang XZ. Wrote the paper: S. Lu FH.

References

  1. 1. Lima-Mendez G, Toussaint A, Leplae R (2007) Analysis of the phage sequence space: the benefit of structured information. Virology 365: 241–249.
  2. 2. Hendrix RW (2003) Bacteriophage genomics. Curr Opin Microbiol 6: 506–511.
  3. 3. Twort A (1993) In focus, out of step: a biography of Frederick William Twort F.R.S., 1877–1950. Phoenix Mill; Dover, NH: A. Sutton. xi, 40 p.
  4. 4. dHerelle F (1917) Sur un microbe invisible antagoniste des bacilles dysentériques. C R Acad Sci Paris 165: 373–375.
  5. 5. Debarbieux L (2008) Experimental phage therapy in the beginning of the 21st century. Med Mal Infect 38: 421–425.
  6. 6. Abedon S (2011) Phage therapy pharmacology: calculating phage dosing. Adv Appl Microbiol 77: 1–40.
  7. 7. Ackermann HW, Prangishvili D (2012) Prokaryote viruses studied by electron microscopy. Arch Virol 157: 1843–1849.
  8. 8. Bonomo RA, Szabo D (2006) Mechanisms of multidrug resistance in Acinetobacter species and Pseudomonas aeruginosa. Clin Infect Dis 43 Suppl 2S49–56.
  9. 9. Skurnik M, Strauch E (2006) Phage therapy: facts and fiction. Int J Med Microbiol 296: 5–14.
  10. 10. Ceyssens PJ, Lavigne R (2010) Bacteriophages of Pseudomonas. Future Microbiol 5: 1041–1055.
  11. 11. Kwan T, Liu J, Dubow M, Gros P, Pelletier J (2006) Comparative genomic analysis of 18 Pseudomonas aeruginosa bacteriophages. J Bacteriol 188: 1184–1187.
  12. 12. Mesyanzhinov VV, Robben J, Grymonprez B, Kostyuchenko VA, Bourkaltseva MV, et al. (2002) The genome of bacteriophage phiKZ of Pseudomonas aeruginosa. J Mol Biol 317: 1–19.
  13. 13. Nakayama K, Kanaya S, Ohnishi M, Terawaki Y, Hayashi T (1999) The complete nucleotide sequence of phi CTX, a cytotoxin-converting phage of Pseudomonas aeruginosa: implications for phage evolution and horizontal gene transfer via bacteriophages. Mol Microbiol 31: 399–419.
  14. 14. Ceyssens PJ, Miroshnikov K, Mattheus W, Krylov V, Robben J, et al. (2009) Comparative analysis of the widespread and conserved PB1-like viruses infecting Pseudomonas aeruginosa. Environ Microbiol 11: 2874–2883.
  15. 15. Uchiyama J, Rashel M, Takemura I, Kato S, Ujihara T, et al. (2012) Genetic characterization of Pseudomonas aeruginosa bacteriophage KPP10. Arch Virol 157: 733–738.
  16. 16. Tan Y, Zhang K, Rao X, Jin X, Huang J, et al. (2007) Whole genome sequencing of a novel temperate bacteriophage of P. aeruginosa: evidence of tRNA gene mediating integration of the phage genome into the host bacterial chromosome. Cell Microbiol 9: 479–491.
  17. 17. Grigoriev A (1999) Strand-specific compositional asymmetries in double-stranded DNA viruses. Virus Res 60: 1–19.
  18. 18. Diard M, Garry L, Selva M, Mosser T, Denamur E, et al. (2010) Pathogenicity-associated islands in extraintestinal pathogenic Escherichia coli are fitness elements involved in intestinal colonization. J Bacteriol 192: 4885–4893.
  19. 19. Wu R, Taylor E (1971) Nucleotide sequence analysis of DNA. II. Complete nucleotide sequence of the cohesive ends of bacteriophage lambda DNA. J Mol Biol 57: 491–511.
  20. 20. Bowden DW, Modrich P (1985) In vitro maturation of circular bacteriophage P2 DNA. Purification of ter components and characterization of the reaction. J Biol Chem 260: 6999–7007.
  21. 21. Wheeler DL, Church DM, Federhen S, Lash AE, Madden TL, et al. (2003) Database resources of the National Center for Biotechnology. Nucleic Acids Res 31: 28–33.
  22. 22. Hendrix RW, Smith MC, Burns RN, Ford ME, Hatfull GF (1999) Evolutionary relationships among diverse bacteriophages and prophages: all the world's a phage. Proc Natl Acad Sci U S A 96: 2192–2197.
  23. 23. Pajunen MI, Kiljunen SJ, Soderholm ME, Skurnik M (2001) Complete genomic sequence of the lytic bacteriophage phiYeO3-12 of Yersinia enterocolitica serotype O:3. J Bacteriol 183: 1928–1937.
  24. 24. Zhang Y, Maley F, Maley GF, Duncan G, Dunigan DD, et al. (2007) Chloroviruses encode a bifunctional dCMP-dCTP deaminase that produces two key intermediates in dTTP formation. J Virol 81: 7662–7671.
  25. 25. Manoharadas S, Witte A, Blasi U (2009) Antimicrobial activity of a chimeric enzybiotic towards Staphylococcus aureus. J Biotechnol 139: 118–123.
  26. 26. Wu H, Lu H, Huang J, Li G, Huang Q (2012) EnzyBase: a novel database for enzybiotic studies. BMC Microbiol 12: 54.
  27. 27. Sun WZ, Tan YL, Jia M, Hu XM, Rao XC, et al. (2010) Functional characterization of the endolysin gene encoded by Pseudomonas aeruginosa bacteriophage PaP1. African Journal of Microbiology Research 4: 933–939.
  28. 28. Williams EA, Degnan SM (2009) Carry-over effect of larval settlement cue on postlarval gene expression in the marine gastropod Haliotis asinina. Mol Ecol 18: 4434–4449.
  29. 29. Ciprandi G, Sormani MP, Filaci G, Fenoglio D (2008) Carry-over effect on IFN-gamma production induced by allergen-specific immunotherapy. Int Immunopharmacol 8: 1622–1625.
  30. 30. Eyer L, Pantucek R, Zdrahal Z, Konecna H, Kasparek P, et al. (2007) Structural protein analysis of the polyvalent Staphylococcal bacteriophage 812. Proteomics 7: 64–72.
  31. 31. Garbe J, Bunk B, Rohde M, Schobert M (2011) Sequencing and characterization of Pseudomonas aeruginosa phage JG004. BMC Microbiol 11: 102.
  32. 32. Debarbieux L, Leduc D, Maura D, Morello E, Criscuolo A, et al. (2010) Bacteriophages Can Treat and Prevent Pseudomonas aeruginosa Lung Infections. Journal of Infectious Diseases 201: 1096–1104.
  33. 33. Bamford DH, Grimes JM, Stuart DI (2005) What does structure tell us about virus evolution? Curr Opin Struct Biol 15: 655–663.
  34. 34. Keppel F, Fayet O, Georgopoulos C (1988) Strategies of bacteriophage DNA replication. New York: Plenum Press. 145–262 p.
  35. 35. Stewart CR, Gaslightwala I, Hinata K, Krolikowski KA, Needleman DS, et al. (1998) Genes and regulatory sites of the “host-takeover module” in the terminal redundancy of Bacillus subtilis bacteriophage SPO1. Virology 246: 329–340.
  36. 36. Bravo A, Alonso JC, Trautner TA (1990) Functional analysis of the Bacillus subtilis bacteriophage SPP1 pac site. Nucleic Acids Res 18: 2881–2886.
  37. 37. Casjens S, Hatfull G, Hendrix R (1992) Evolution of dsDNA tailed-bacteriophage genomes. Semin Virol 3: 383–397.
  38. 38. Kimelman A, Levy A, Sberro H, Kidron S, Leavitt A, et al. (2012) A vast collection of microbial genes that are toxic to bacteria. Genome Res 22: 802–809.
  39. 39. Monod C, Repoila F, Kutateladze M, Tetart F, Krisch HM (1997) The genome of the pseudo T-even bacteriophages, a diverse group that resembles T4. J Mol Biol 267: 237–249.
  40. 40. Sambrook J, Russel DW (2001) Molecular Cloning: A Laboratory Manual, 3rd ed. New York: Cold Spring Harbor Laboratory Press. 170 p.
  41. 41. Garbe J, Wesche A, Bunk B, Kazmierczak M, Selezska K, et al. (2010) Characterization of JG024, a Pseudomonas aeruginosa PB1-like broad host range phage under simulated infection conditions. BMC Microbiol 10: 301.
  42. 42. Govind R, Fralick JA, Rolfe RD (2006) Genomic organization and molecular characterization of Clostridium difficile bacteriophage PhiCD119. J Bacteriol 188: 2568–2577.
  43. 43. Casas V, Rohwer F (2007) Phage metagenomics. Methods Enzymol 421: 259–268.
  44. 44. Zheng ZL, Advani A, Melefors O, Glavas S, Nordstrom H, et al. (2011) Titration-free 454 sequencing using Y adapters. Nature Protocols 6: 1367–1376.
  45. 45. de la Bastide M, McCombie WR (2007) Assembling genomic DNA sequences with PHRAP. Curr Protoc Bioinformatics Chapter 11: Unit11 14.
  46. 46. Rosseel T, Scheuch M, Hoper D, De Regge N, Caij AB, et al. (2012) DNase SISPA-next generation sequencing confirms Schmallenberg virus in Belgian field samples and identifies genetic variation in Europe. PLoS One 7: e41967.
  47. 47. Carver T, Thomson N, Bleasby A, Berriman M, Parkhill J (2009) DNAPlotter: circular and linear interactive genome visualization. Bioinformatics 25: 119–120.
  48. 48. Schattner P, Brooks AN, Lowe TM (2005) The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res 33: W686–689.
  49. 49. Besemer J, Borodovsky M (1999) Heuristic approach to deriving models for gene finding. Nucleic Acids Res 27: 3911–3920.
  50. 50. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.
  51. 51. Reese MG (2001) Application of a time-delay neural network to promoter annotation in the Drosophila melanogaster genome. Comput Chem 26: 51–56.
  52. 52. Moak M, Molineux IJ (2004) Peptidoglycan hydrolytic activities associated with bacteriophage virions. Mol Microbiol 51: 1169–1183.
  53. 53. Welker M, Erhard M (2007) Consistency between chemotyping of single filaments of Planktothrix rubescens (cyanobacteria) by MALDI-TOF and the peptide patterns of strains determined by HPLC-MS. Journal of Mass Spectrometry 42: 1062–1068.
  54. 54. Morello E, Saussereau E, Maura D, Huerre M, Touqui L, et al. (2011) Pulmonary bacteriophage therapy on Pseudomonas aeruginosa cystic fibrosis strains: first steps towards treatment and prevention. PLoS One 6: e16963.
  55. 55. Krumsiek J, Arnold R, Rattei T (2007) Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics 23: 1026–1028.
  56. 56. Carver T, Berriman M, Tivey A, Patel C, Bohme U, et al. (2008) Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics 24: 2672–2676.
  57. 57. Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, et al. (2003) Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res 31: 3497–3500.
  58. 58. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731–2739.
  59. 59. Som A, Fuellen G (2009) The effect of heterotachy in multigene analysis using the neighbor joining method. Mol Phylogenet Evol 52: 846–851.
  60. 60. Zafar N, Mazumder R, Seto D (2002) CoreGenes: a computational tool for identifying and cataloging “core” genes in a set of small genomes. BMC Bioinformatics 3: 12.
  61. 61. Mahadevan P, King JF, Seto D (2009) Data mining pathogen genomes using GeneOrder and CoreGenes and CGUG: gene order, synteny and in silico proteomes. Int J Comput Biol Drug Des 2: 100–114.
  62. 62. Klumpp J, Dorscht J, Lurz R, Bielmann R, Wieland M, et al. (2008) The terminally redundant, nonpermuted genome of Listeria bacteriophage A511: a model for the SPO1-like myoviruses of gram-positive bacteria. J Bacteriol 190: 5753–5765.
  63. 63. Just W, Klotz G (1990) Terminal redundancy and circular permutation of mycoplasma virus L3 DNA. J Gen Virol 71 (Pt 9): 2157–2162.
  64. 64. Alemayehu D, Casey PG, McAuliffe O, Guinane CM, Martin JG, et al. (2012) Bacteriophages phiMR299-2 and phiNH-4 can eliminate Pseudomonas aeruginosa in the murine lung and on cystic fibrosis lung airway cells. MBio 3: e00029–00012.
  65. 65. Ceyssens PJ, Hertveldt K, Ackermann HW, Noben JP, Demeke M, et al. (2008) The intron-containing genome of the lytic Pseudomonas phage LUZ24 resembles the temperate phage PaP3. Virology 377: 233–238.
  66. 66. Hertveldt K, Lavigne R, Pleteneva E, Sernova N, Kurochkina L, et al. (2005) Genome comparison of Pseudomonas aeruginosa large phages. J Mol Biol 354: 536–545.