Phenotypic and genotypic variations within a single bacteriophage species

  • Pieter-Jan Ceyssens1Email author,

    Affiliated with

    • Thea Glonti2Email author,

      Affiliated with

      • ndrew M Kropinski3, 4,

        Affiliated with

        • Rob Lavigne1Email author,

          Affiliated with

          • Nina Chanishvili2,

            Affiliated with

            • Leonid Kulakov5,

              Affiliated with

              • Nino Lashkhi2,

                Affiliated with

                • Marina Tediashvili2 and

                  Affiliated with

                  • Maya Merabishvili2, 6

                    Affiliated with

                    Virology Journal20118:134

                    DOI: 10.1186/1743-422X-8-134

                    Received: 19 January 2011

                    Accepted: 23 March 2011

                    Published: 23 March 2011

                    Abstract

                    Background

                    Although horizontal gene transfer plays a pivotal role in bacteriophage evolution, many lytic phage genomes are clearly shaped by vertical evolution. We investigated the influence of minor genomic deletions and insertions on various phage-related phenotypic and serological properties.

                    Findings

                    We collected ten different isolates of Pseudomonas aeruginosa bacteriophage ϕKMV. All sequenced genomes (42-43 kb, long direct terminal repeats) are nearly identical, which intuitively implied strongly similar infections cycles. However, their latent periods vary between 21 and 28 minutes and they are able to lyse between 5 and 58% of a collection of 107 clinical P. aeruginosa strains. We also noted that phages with identical tail structures displayed profound differences in host spectra. Moreover, point mutations in tail and spike proteins were sufficient to evade neutralization by two phage-specific antisera, isolated from rabbits.

                    Conclusion

                    Although all analyzed phages are 83-97% identical at the genome level, they display a surprisingly large variation in various phenotypic properties. The small overlap in host spectrum and their ability to readily escape immune defences against a nearly identical phage are promising elements for the application of these phages in phage therapy.

                    The potential use of lytic bacteriophages as alternative therapeutic agents against antibiotic-resistant bacteria has been widely documented [13]. This renewed interest in phage therapy resulted in human clinical trials [4, 5] and a considerable rise of commercial interest [6]. Potential downsides of phage applications, i.e. their relatively narrow host range and emerging bacterial resistance, are often countered by the proposed use of mixtures (cocktails) of phages targeting different hosts and/or host receptors. A popular argument in favour of these cocktails is the sheer abundance of phages in nature (1031 particles worldwide) and the ease of isolating new phages infecting phage-resistant bacteria [7].

                    In spite of their ubiquitous character, numerous studies suggest the existence of only a limited number of virulent phage 'types' targeting a specific bacterial species. For example, despite decades of intensive research, only 17 and 13 distinct phage species infecting the model organisms Escherichia coli and Pseudomonas aeruginosa are known [8, 9]. Nearly all newly isolated lytic phages infecting these organisms cluster within existing species and carry genomes which are often over 90% identical at the nucleotide level [1012]. As this limited global diversity might hamper the composition of truly diverse phage cocktails, we were interested in the impact of these subtle "intraspecies" genomic insertions and deletions on infection-related and serological properties of lytic bacteriophages.

                    As a model species, we choose the virulent P. aeruginosa phage ϕKMV which resembles the classic coliphage T7 in morphology and overall genome architecture [13]. It was also the first known T7-related phage encoding a single-subunit RNA polymerase gene downstream its DNA metabolism genes instead of in the early genomic region [13]. Up to now, three phages infecting P. aeruginosa (LKD16, ϕKF77 and LUZ19) have been reported to be 83-90% identical to ϕKMV at the nucleotide level [1416]. In addition, phage LUZ2 was reported to be closely related to ϕKMV based on de novo analysis of structural phage proteins [12].

                    Phages PT2, PT5, PT6, PNC101 and PNM were isolated from Mtkvari River in Georgia [17] and added to our existing collection of ϕKMV-related isolates (Table 1). We also included the more distantly related Pseudomonas phage LKA1 in this analysis; LKA1 shares significant sequence similarity for only 48% of its gene products but has a conserved functional genome organization and gene order [14]. With the exception of LKA1 which forms small plaques (1 mm), all of these phages produce large (4-10 mm diameter) clear plaques. As members of the family Podoviridae, these viruses possess icosahedral heads (diameter of 50-60 nm) and 10-12 nm noncontractile tails. With the exception of LKA1, all the analyzed phages lyse P. aeruginosa cells 21 to 28 minutes after infection at 37°C upon the release of between 116 and 166 newly produced virions (Table 1).
                    Table 1

                    Overview of ϕKMV-related bacteriophages included in this study, all infecting P. aeruginos a.

                    Phage

                    Isolated in

                    Genome (bp)

                    # ORFs

                    DTR length

                    GC (%)

                    ϕKMV identitya

                    Accession number

                    Kadsorption(ml/min)b

                    Latent period (min)

                    Burst size (PFU/cell)

                    KAPS/PT5 c

                    KAPS/PNC101 c

                    ϕKMV

                    Moscow, Russia

                    42,519

                    52

                    414

                    62.3

                    100

                    NC_005045

                    1.6 × 10-8

                    28

                    116

                    259

                    849

                    LKD16

                    Leuven, Belgium

                    43,200

                    54

                    428

                    62.6

                    83

                    NC_009935

                    -

                    27

                    120

                    87.5

                    -

                    LUZ19

                    Leuven, Belgium

                    43,548

                    54

                    472

                    62.2

                    89

                    NC_010326

                    1.2 × 10-7

                    22

                    122

                    -

                    -

                    ϕKF77

                    Moscow, Russia

                    43,152

                    53

                    454

                    62.8

                    90

                    NC_012418

                    1.5 × 10-8

                    23

                    129

                    56.2

                    -

                    PNM

                    Tbilisi, Georgia

                    42,721

                    51

                    422

                    62.3

                    94

                    -

                    ND

                    ND

                    ND

                    -

                    23

                    PT2

                    Tbilisi, Georgia

                    42,961

                    54

                    488

                    62.0

                    97

                    NC_011107

                    9.2 × 10-9

                    21

                    116

                    737

                    51

                    PT5

                    Tbilisi, Georgia

                    42,954

                    52

                    413

                    62.0

                    95

                    NC_011105

                    7.3 × 10-9

                    21

                    134

                    598

                    679

                    PT6

                    Tbilisi, Georgia

                    ND

                    ND

                    ND

                    ND

                    ND

                    -

                    2.1 × 10-8

                    25

                    166

                    257

                    857

                    PNC101

                    Tbilisi, Georgia

                    ND

                    ND

                    ND

                    ND

                    ND

                    -

                    ND

                    ND

                    ND

                    60

                    753

                    LUZ2

                    Leuven, Belgium

                    ND

                    ND

                    ND

                    ND

                    ND

                    -

                    ND

                    ND

                    ND

                    270

                    -

                    LKA1

                    Leuven, Belgium

                    41,593

                    56

                    298

                    60.9

                    ND

                    NC_009936

                    3.9 × 10-9

                    41

                    225

                    -

                    -

                    a Identity with ϕKMV, expressed as % identical nucleotides according to the Stretcher algorithm (EMBOSS suite).

                    b The adsorption constant k was calculated from the unabsorbed phage at a multiplicity of infection of 1:1 by applying the formula ka = 2.3/Bt log P0/P, where P0 is initial number of phage, P is number of unabsorbed phage particles after time t and B is number of bacteria (4). ND, not determined.

                    c Serological neutralization assay using PT5 and PNC101 APS. Dilutions (1/50- 1/1000) of APS were mixed with 106 bacteriophage, samples were titered after different incubation times (5, 15 and 30 min) at 37°C. K, rate constant of neutralization, defined as k = 2.3 D/t logP0/P, where P0 - count of bacteriophage plaques at zero time, P - count of bacteriophage plaques after t time, D - dilution of APS.

                    Infectivity screens on various outer membrane mutants of P. aeruginosa unambiguously identified the necessity for Type IV pili for successful infection of all ϕKMV-like phages except for LKA1, which depends on algC expression. The presence/absence of bacterial flagella (encoded by fliC) has no influence on phage infection (data not shown). Despite this common type IV pili dependence, subsequent screens on a large collection of 114 clinical P. aeruginosa strains showed profound differences in host range. Between 18 (ϕKMV) and 62 (PT6) P. aeruginosa strains were sensitive to phage infection, while the more distantly related phage LKA1 can only lyse six strains from our collection (Figure 1A, Additional file 1). This observation pointed out variations in secondary adsorption mechanisms and intracellular development, and provoked interest in the genome sequences of these new ϕKMV-related isolates.
                    http://static-content.springer.com/image/art%3A10.1186%2F1743-422X-8-134/MediaObjects/12985_2011_1247_Fig1_HTML.jpg
                    Figure 1

                    Host range analysis of the ϕKMV-like viruses. A. Host range analysis of 'phiKMV-like viruses' using a collection of 114 clinical P. aeruginosa strains. In this assay, 106 pfu were spotted on a bacterial lawn and checked for lysis. See additional files for more information on the P. aeruginosa strains and for detailed spectra. B. Host range comparison of four phages displaying >94% nucleotide identity throughout their genomes, stating the number of strains which are uniquely or commonly lysed by PNM, PT2, PT5 and phiKMV.

                    Using traditional Sanger dideoxy methods for PNM and 454 pyrosequencing for PT2 and PT5, the complete genome sequences of these phages were determined (Table 1, Figure 1). All three phages display high nucleotide similarity (94-97% identity) to ϕKMV, except for a 382 bp region located between ORF3 and 4 (Figure 1), which was earlier associated with a localized single-stranded nick on the non-coding strand in the related phages LUZ19 and ϕKF77 [16]. As ϕKMV is the only phage of this species which has lost this nick, it is most probably a more recent representative of this phage species.

                    When comparing these newly sequenced genomes to all known ϕKMV-related isolates at the protein level, the extraordinary conservation throughout their genomes is evident. With the exception of occasional insertions (e.g., the unrelated ORFs 17.1 in LKD16 and ϕKF77) and deletions (e.g., ORF2 in ϕKF77, ORF20 in PNM), the genomes have been stably maintained over time and location. The only phage encoding an additional structural protein is LUZ19 (gp49), which might cause subtle differences in virion. The tail fiber regions of ϕKMV, PT2 and PNM (and of the partially sequenced phage PT6) are 99.5-100% identical, strongly suggesting identical reception/adhesion mechanisms. A similar conservation of tail genes is also present between LUZ19 and ϕKF77 (Figure 2). In contrast, these genes are far less conserved in PT5 (25-52% amino acid identity to ϕKMV), whose tail fibers are more related to phage LKD16 (61-94% amino acid identity). In contrast to clear differences in the tail fibers genes of PT5 and PNM, their early coding region is completely identical (Figure 2). Intriguingly, the fact that ϕKMV, PT2, PNM and PT5 have nearly identical virions and/or early regions does not correlate to their respective host spectra. As shown in Figure 1B, only six P. aeruginosa strains are lysed by all four phages and 30 strains can only be targeted by a single phage.
                    http://static-content.springer.com/image/art%3A10.1186%2F1743-422X-8-134/MediaObjects/12985_2011_1247_Fig2_HTML.jpg
                    Figure 2

                    Comparative genomic analysis of the ϕKMV-like viruses. The predicted open reading frames of all sequenced ϕKMV-like phages are shown. Their amino acid identity to the corresponding ORF of ϕKMV is indicated in different shades of grey. ORFs unique to each phage are not colored, while predicted functions of the annotated genes are indicated on top. Functional equivalent genes are connected with broken lines. Gene annotation was performed as described elsewhere (Ceyssens et al., 2006).

                    In a final step, the serological relatedness of these phages was studied by using a classical cross neutralization method [18]. Two anti-phage sera (APS) were obtained by a two-step immunization of three Swiss white rabbits with phages PT5 and PNC101. Each rabbit was injected intramuscularly with a 1:1 (vol/vol) mixture of phage (1010 pfu/ml) and adjuvant (Freund's complete adjuvant, Difco). One month later, a secondary immunization was performed by injecting phage (1 ml, 1011 pfu/ml) without adjuvant. After two weeks, blood was taken from the left ventricle of the heart. This blood was allowed to clot, the APS was filtered (0.45 μm pore size) and stored at 4°C. Subsequent neutralization assays with these APS showed significant variation in serological response among the tested phages, as only ϕKMV and PT5 were efficiently neutralized by both antisera (Table 1). Phages LKA1 and LUZ19 were not neutralized to any significant degree by either APS, while the inactivation of LKD16, ϕKF77 and PNM was at least tenfold less efficient than of the phages used for the APS production.

                    When looking into these results in more detail, some surprising observations were made. One would expect that neutralizing antibodies primarily function through interaction with capsid and/or tail (fibers) proteins of these phages. It is therefore curious that PT2 but not PNM is neutralized by APSPT5, since the structural proteins of PT2 and PNM are all >99.5% identical (Figure 2). Comparing the genomes of these three phages, the only non-silent structural mutations present in PT2 and PT5 in comparison to PNM are found in the tail tubular proteins A (P20-L20, K55-R55, T102-P102, R115-H115) and B (D26-G26, T111-N111, V142-I142, A298-G298, G329-D329, R363-S363, A644-V644, P660-S660, A784-V784). These amino acids substitutions are conserved in ϕKMV, explaining the sensitivity of this phage for APSPT5. A similar finding was observed using the anti-phage serum produced against the non-sequenced phage PNC101. This serum was over tenfold more effective in neutralizing ϕKMV and PT5 than PT2 and PNM, despite 100% identity of their capsid and connector proteins (Table 1, Figure 2). This time, the difference might be correlated to four amino acid substitutions in the spike protein (gp48) of PT2 and PNM.

                    Discussion

                    Although horizontal gene transfer plays a pivotal role in bacteriophage evolution, many lytic phage genomes are clearly shaped by vertical evolution. The high level of genomic conservation (83-97% identity) observed within the 'phiKMV-like viruses' is a common theme among virulent phages infecting the same bacterial host; for example, the genome of Roseophage SIO1 was found to be completely conserved in isolates sampled over many years in various places [11]. Typically, strongly related phages attach to the same receptor and carry out a highly comparable infection cycle. At first glace, this is indeed the case for all analyzed ϕKMV-like viruses. They are all pili-dependent and quickly lyse their host upon the release of up to 166 newly produced particles.

                    During this study we noticed that small "intraspecies" genomic variations can have essential phenotypic consequences towards the applications of these phages in therapeutic settings. First of all, only limited overlaps in host spectrum exist between isolates with identical tail fibers and/or early genome regions, implying that minor genomic changes can cause a significant shift in infectivity range. Although these differences could be attributed to the evasion of CRISPR repeats [19], a recent survey did not find a single spacer matching a lytic bacteriophage sequence in 122 clinical P. aeruginosa strains [20]. Alternatively, small point mutations could help in the evasion of host restriction-modification systems [21]. In any case, this illustrates the versatility of a phage genome to evade host defenses.

                    Moreover, phage therapy seems not be hindered by the adaptive immune system, since minor variations in tail/capsid structures seem sufficient to evade antibody binding. This confirms observations made by Vitiello and colleagues [22], who reported a long-circulating mutant of phage λ which carried only a single mutation in the major capsid protein. As phages have a virtually unlimited potential of subtly modifying their virion, the potential is present to quickly overcome host immune responses and bacterial resistance during phage therapy.

                    Declarations

                    Acknowledgements

                    A.M.K. is supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. P.J.C. is a postdoctoral fellow supported by the 'Fonds Wetenschappelijk Onderzoek (F.W.O)' of the Flemish Government. This research was performed by members of the 'PhageBiotics' research community, funded by the FWO Vlaanderen (WO.022.09).

                    Authors’ Affiliations

                    (1)
                    Laboratory of Gene Technology (LoGT), Katholieke Universiteit Leuven
                    (2)
                    Eliava Institute of Bacteriophage, Microbiology and Virology (EIBMV)
                    (3)
                    Laboratory for Foodborne Zoonoses, Public Health Agency of Canada
                    (4)
                    Department of Molecular & Cellular Biology, University of Guelph
                    (5)
                    School of Biological Sciences, The Queen’s University of Belfast
                    (6)
                    Laboratory for Molecular and Cellular Technology (LabMCT), Burn Wound Center, Queen Astrid Military Hospital

                    References

                    1. Housby JN, Mann NH: Phage therapy. Drug Discov Today 2009, 14:536–540.PubMedView Article
                    2. Fischetti VA: Bacteriophage endolysins: a novel anti-infective to control Gram-positive pathogens. Int J Med Microbiol 2010, 300:357–362.PubMedView Article
                    3. Górski A, Miedzybrodzki R, Borysowski J, Weber-Dabrowska B, Lobocka M, Fortuna W, Letkiewicz S, Zimecki M, Filby G: Bacteriophage therapy for the treatment of infections. Curr Opin Investig Drugs 2009, 10:766–774.PubMed
                    4. Wright A, Hawkins CH, Anggård EE, Harper DR: A controlled clinical trial of a therapeutic bacteriophage preparation in chronic otitis due to antibiotic-resistant Pseudomonas aeruginosa ; a preliminary report of efficacy. Clin Otolaryngol 2009, 34:349–357.PubMedView Article
                    5. Rhoads DD, Wolcott RD, Kuskowski MA, Wolcott BM, Ward LS, Sulakvelidze A: Bacteriophage therapy of venous leg ulcers in humans: results of a phase I safety trial. J Wound Care 2009, 18:237–238.PubMed
                    6. Deresinski S: Bacteriophage therapy: exploiting smaller fleas. Clin Infect Dis 2009, 48:1096–1101.PubMedView Article
                    7. Loc-Carrillo C, Abedon ST: Pros and cons of phage therapy. Bacteriophage 1:1–4.
                    8. Casjens SR: Diversity among the tailed-bacteriophages that infect the Enterobacteriaceae. Res Microbiol 2008, 159:340–348.PubMedView Article
                    9. Ceyssens PJ, Lavigne R: Bacteriophages of Pseudomonas. Future Microbiol 2010, 5:1041–1055.PubMedView Article
                    10. Savalia D, Westblade LF, Goel M, Florens L, Kemp P, Akulenko N, Pavlova O, Padovan JC, Chait BT, Washburn MP, Ackermann HW, Mushegian A, Gabisonia T, Molineux I, Severinov K: Genomic and proteomic analysis of phiEco32, a novel Escherichia coli bacteriophage. J Mol Biol 2008, 377:774–789.PubMedView Article
                    11. Angly F, Youle M, Nosrat B, Srinagesh S, Rodriguez-Brito B, McNairnie P, Deyanat-Yazdi G, Breitbart M, Rohwer F: Genomic analysis of multiple Roseophage SIO1 strains. Environ Microbiol 2009, 11:2863–2873.PubMedView Article
                    12. Ceyssens PJ, Noben JP, Ackermann HW, Verhaegen J, De Vos D, Pirnay JP, Merabishvili M, Vaneechoutte M, Chibeu A, Volckaert G, Lavigne R: Survey of Pseudomonas aeruginosa and its phages: de novo peptide sequencing as a novel tool to assess the diversity of worldwide collected viruses. Environ Microbiol 2009, 11:1303–1313.PubMedView Article
                    13. Lavigne R, Burkal'tseva MV, Robben J, Sykilinda NN, Kurochkina LP, Jonckx B, Grymonprez B, Krylov VN, Mesyanzhinov VV, Volckaert G: The genome of bacteriophage ϕKMV, a T7-like virus infecting Pseudomonas aeruginosa . Virology 2003, 312:49–59.PubMedView Article
                    14. Ceyssens PJ, Lavigne R, Mattheus W, Chibeu A, Hertveldt K, Mast J, Robben J, Volckaert G: Genomic analysis of Pseudomonas aeruginosa phages LKD16 and LKA1: establishment of the ϕKMV subgroup within the T7 supergroup. J Bacteriol 2006, 188:6924–6931.PubMedView Article
                    15. Lammens E, Ceyssens PJ, Voet M, Hertveldt K, Lavigne R, Volckaert G: Representational Difference Analysis (RDA) of bacteriophage genomes. J Microbiol Methods 2009, 77:207–213.PubMedView Article
                    16. Kulakov LA, Ksenzenko VN, Shlyapnikov MG, Kochetkov VV, Del Casale A, Allen CC, Larkin MJ, Ceyssens PJ, Lavigne R: Genomes of "phiKMV-like viruses" of Pseudomonas aeruginosa contain localized single-strand interruptions. Virology 2009, 391:1–4.PubMedView Article
                    17. Glonti T, Chanishvili N, Taylor PW: Bacteriophage-derived enzyme that depolymerases the alginic acid capsule associated with cystic fibrosis isolates of Pseudomonas aeruginosa . J Appl Microbiol 2010, 108:695–702.PubMedView Article
                    18. Waddell TE, Franklin K, Mazzocco A, Johnson RP: Preparation and characterization of anti-phage serum. Methods Mol Biol 2009, 501:287–292.PubMedView Article
                    19. Deveau H, Garneau JE, Moineau S: CRISPR/Cas system and its role in phage-bacteria interactions. Annu Rev Microbiol 2010, 64:475–493.PubMedView Article
                    20. Cady KC, White AS, Hammond JH, Abendroth MD, Karthikeyan RS, Prajna L, Zegans ME, O'Toole GA: Prevalence, conservation and functional analysis of Yersinia and Escherichia CRISPR regions in clinical Pseudomonas aeruginosa Isolates. Microbiology 157:430–437.
                    21. Labrie SJ, Samson JE, Moineau S: Bacteriophage resistance mechanisms. Nat Rev Microbiol 2010, 8:317–327.PubMedView Article
                    22. Vitiello CL, Merril CR, Adhya S: An amino acid substitution in a capsid protein enhances phage survival in mouse circulatory system more than a 1000-fold. Virus Res 2005, 114:101–103.PubMedView Article

                    Copyright

                    © Ceyssens et al; licensee BioMed Central Ltd. 2011

                    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.

                    Advertisement